Github Repository
4. Flächenberechnung und Newtonverfahren (1D)
Contents
4. Flächenberechnung und Newtonverfahren (1D)¶
Aufgabe dieser Übung ist es, die optimale Flächennutzung von Solarpanelen für ein kreisförmiges Dach zu berechnen.
Es stehen zwei Module zur Auswahl. Die Berechnung erfolgt mittels des eindimensionalen Newtonverfahren.
Situation¶
Auf das kreisförmige Dach eines Lagergebäudes sollen Photovoltaik-Module installiert werden. Zur Auswahl stehen zwei verschiedene Module:
Modul 1: 0,85m x 1,74m, 133€
Modul 2: 0,98m x 1,62m, 143€
Ziel ist eine maximale Ausnutzung der Dachfläche.
Hinweis
Im Rahmen der Aufgabe wird angenommen, dass die Module alle entlang eines Rechteck-Rasters montiert werden. Eine verschobene Anordnung (ähnlich einem Mauerwerk) wird nicht betrachtet.
Tipp
Berechnen Sie zunächst alle nötigen Werte für einen Viertelkreis und rechnen Sie anschließend alles auf einen vollen Kreis hoch.
1. Kreis-Plot¶
Erzeugen Sie ein neues Skript und legen Sie dort die Functions main und circle an. Ergänzen Sie nun die Funktion circle so, dass Sie die y-Werte eines Kreises abhängig von x ausgibt. Der Radius des Dachs beträgt 10 m.
Plotten Sie nun einen Viertelkreis.
%% file main.m
function main
end
function circle
axis equal % Legt gleiche Seitenverhältnisse für x- und y-Achse im Plot fest
end
Hinweis
Die Formel eines Kreises lautet \(y=\sqrt{r^2-x^2}\).
2. Berechnung der Modulverteilung¶
Zur Berechnung der Modulverteilung sollten Sie wie folgt vorgehen (Modulbreite: \(xw\), Modulhöhe: \(yw\)):
Start einer Schleife bei \(x(i)\), z.B. \(x = 0\).
\(x\) um \(xw\) erhöhen.
Prüfung, ob Kreis bei \(x(i) = x(i-1) + xw\) noch oberhalb von \(yw\) liegt.
Wenn ja: Anzahl an Module in diesem Intervall ist gleich y-Wert des Kreises bei \(x(i)\) geteilt durch Modulhöhe \(yw\) (auf ganze Zahl abrunden).
Wenn nein: Anzahl Module ist gleich Null.
Wieder bei a beginnen. Modulanzahlen bei jedem Schleifendurchlauf aufaddieren. Schleife solange durchlaufen lassen, bis \(x(i)\) größer als \(r\) ist.
% your code here
3. Plotten der Modulverteilung¶
Plotten Sie die Modulverteilung für einen Viertelkreis mit der MATLAB-Funktion rectangle.
% Plottet ein Rechteck mit unterer, linker Ecke bei x und y sowie den Kantenlängen xw und yw:
rectangle(‘Position‘,[x,y,xw,yw])
Tipp
Ergänzen Sie die Berechnung aus 2. so, dass Sie pro Durchlauf die Anzahl der Module auf einem x-Intervall speichern. Somit wissen Sie für den Plot, wie viele Module pro x-Schritt erzeugt werden müssen.
4. Auswertung¶
Geben Sie an, wie viele Module jeweils maximal auf der gesamten Dachfläche platziert werden können und welche Fläche diese Module einnehmen würden. Ist das Modul für mehr Gesamtfläche die günstigere oder teurere Wahl?
5. Berechnung Einspeisedauer¶
Ein beispielhafter Tagesgang der PV-Leistung kann vereinfacht über folgende Funktion dargestellt werden (\(t\) in Stunden von \(0…24\), \(P\) normalisiert als \(0…1\)):
Plotten Sie die Funktion über 24h.
% your code here
Die Leistung der PV-Anlage ist so dimensioniert, dass in der Woche (Mo-Fr) nur etwa 85%, am Wochenende nur ca. 30% der maximal verfügbaren Leistung benötigt werden. Der Rest wird gespeichert oder ins Netz eingespeist.
Bestimmen Sie numerisch, mit dem Newtonverfahren, von wann bis wann jeweils überschüssige Leistung produziert wird. Ermitteln Sie daraus die Anzahl an Stunden pro Woche, an welchen Leistung eingespeist wird.
% your code here
Hinweis
Die Vorschrift des eindimensionalen Newtonverfahrens lautet:
Die Ableitung \(f'(x)\) können Sie zuvor analytisch bestimmen.
\(a\) ist ein Dämpfungsparameter (im Normalfall \(= 1\)). Falls Sie Probleme bei der Lösung bekommen, können Sie \(0<a<1\) setzen.