Energiehaushalt eines Hauses: Energieerzeugung EEZ: Unterschied zwischen den Versionen
Zeile 254: | Zeile 254: | ||
== Komponente Flächenberechnung == | == Komponente Flächenberechnung == | ||
'''Eingänge''' | '''Eingänge und Parameter''' | ||
* PAR_LKT_Grundflaeche = <math>100\,m^2</math> | * PAR_LKT_Grundflaeche = <math>100\,m^2</math> | ||
* PAR_LKT_Deckenhoehe = <math>2{,}5\,m</math> | * PAR_LKT_Deckenhoehe = <math>2{,}5\,m</math> | ||
Zeile 287: | Zeile 287: | ||
== Komponente Wärmeverlust == | == Komponente Wärmeverlust == | ||
'''Eingänge und Parameter''' | |||
* PAR_LKT_Grundflaeche = <math>100\,m^2</math> | |||
* ISO_FlaecheWand = <math>75\,m^2</math> | |||
* ISO_FlaecheFenster = <math>25\,m^2</math> | |||
* PAR_LKT_Temperatur aus Temperaturdaten.mat | |||
* ISO_IstTemp = <math>25\,^{\circ}C</math> | |||
* PAR_ISO_WaendeUWert = <math>0{,}187\,W/(m^2K)</math> | |||
* PAR_ISO_DachUWert = <math>0{,}188\,W/(m^2K)</math> | |||
* PAR_ISO_FensterUWert = <math>0{,}5\,W/(m^2K)</math> | |||
* PAR_ISO_BodenUWert = <math>0{,}073\,W/(m^2K)</math> | |||
'''Ausgänge''' | |||
* ISO_Waermeverlust [W] | |||
* ISO_TempDiff [°C] | |||
{| class="wikitable" | |||
{| class="wikitable mw-collapsible mw-collapsed" | |||
| align="center" style="background:#4472C4; color:white"|'''Testfall-ID''' | |||
| align="center" style="background:#4472C4; color:white"|'''Testfall-Name''' | |||
| align="center" style="background:#4472C4; color:white"|'''Anforderung''' | |||
| align="center" style="background:#4472C4; color:white"|'''Vorbedingungen und Eingänge''' | |||
| align="center" style="background:#4472C4; color:white"|'''Aktionen''' | |||
| align="center" style="background:#4472C4; color:white"|'''Erwartetes Ergebnis''' | |||
| align="center" style="background:#4472C4; color:white"|'''Ergebnis''' | |||
| align="center" style="background:#4472C4; color:white"|'''Bewertung''' | |||
| align="center" style="background:#4472C4; color:white"|'''Kommentar''' | |||
|- | |||
| 1 || Test auf die Aktualisierung der Dateien || 001 || Software SVN vorbereiten. || Update SVN || Keine Fehlermeldung || Keine Fehlermeldung || i. O. || | |||
|- | |||
| 2 || Test auf die Matlab/Simulink Version || 002 || MATLAB R2022a + Simulink Version 10.5 installieren. || MATLAB Skript "start.m" mit MATLAB R2022a öffnen und im Editor "Run" drücken || Haus.slx wird ohne Fehlermeldung geöffnet || Haus.slx wird ohne Fehlermeldung geöffnet || i. O. || | |||
|- | |||
| 3 || Test auf alle Eingänge = 0 || 028 || Parameter geladen. || Simulation eines Zeitschritts || Alle Ausgänge = 0 || Alle Ausgänge = 0 || i. O. || | |||
|- | |||
| 4 || Test auf Laden der Parameter PAR_LKT_Temperatur, ISO_IstTemp, ISO_TempDiff || 047-049 || Parameter geladen. || Run MATLAB Parameterdatei || Parameter im Workspace || Parameter im Workspace || i. O. || | |||
|- | |||
| 5 || Test des Wärmeverlusts der Wandfläche || 023 || Parameter geladen. PAR_FlaecheWand = 75; PAR_Flaeche_Fenster = 0; PAR_Grundflache = 0; PAR_LKT_Temperatur = 10 °C; ISO_IstTemp = 25 °C; || Simulation eines Zeitschritts || Waermeverlust = <math>-210{,}4\,W</math> || Waermeverlust = <math>-210{,}4\,W</math> || i. O. || | |||
|- | |||
| 6 || Test des Wärmeverlusts der Fensterfläche || 016 || Parameter geladen. PAR_FlaecheWand = 0; PAR_Flaeche_Fenster = 25; PAR_Grundflache = 0; PAR_LKT_Temperatur = 10 °C; ISO_IstTemp = 25 °C; || Simulation eines Zeitschritts || Waermeverlust = <math>-187{,}5\,W</math> || Waermeverlust = <math>-187{,}5\,W</math> || i. O. || | |||
|- | |||
| 7 || Test des Wärmeverlusts der Bodenfläche || 022 || Parameter geladen. PAR_FlaecheWand = 0; PAR_Flaeche_Fenster = 0; PAR_Grundflache = 100; PAR_LKT_Temperatur = 10 °C; ISO_IstTemp = 25 °C; || Simulation eines Zeitschritts || Waermeverlust = <math>-391{,}5\,W</math> || Waermeverlust = <math>-391{,}5\,W</math> || i. O. || | |||
|- | |||
| 8 || Test des Wärmeverlusts des Dachs || 014 || Parameter geladen. PAR_FlaecheWand = 0; PAR_Flaeche_Fenster = 0; PAR_Grundflache = 100; PAR_LKT_Temperatur = 10 °C; ISO_IstTemp = 25 °C; || Simulation eines Zeitschritts || Waermeverlust = <math>-391{,}5\,W</math> || Waermeverlust = <math>-391{,}5\,W</math> || i. O. || | |||
|} | |||
== Komponente Wärmekapazität == | == Komponente Wärmekapazität == | ||
'''Eingänge und Parameter''' | |||
* ISO_FlaecheWand = <math>100\,m^2</math> | |||
* ISO_Grundflaeche = <math>100\,m^2</math> | |||
* ISO_FlaecheFenster = <math>25\,m^2</math> | |||
* PAR_ISO_WaermekapazitaetWand = <math>353\,kJ/(m^2K)</math> | |||
* PAR_ISO_WaermekapazitaetDach = <math>448\,kJ/(m^2K)</math> | |||
* PAR_ISO_WaermekapazitaetFenster = <math>21{,}6\,kJ/(m^2K)</math> | |||
* PAR_ISO_WaermekapazitaetBoden = <math>63\,kJ/(m^2K)</math> | |||
'''Ausgänge''' | |||
* ISO_Waermekapazitaet | |||
{| class="wikitable" | |||
{| class="wikitable mw-collapsible mw-collapsed" | |||
| align="center" style="background:#4472C4; color:white"|'''Testfall-ID''' | |||
| align="center" style="background:#4472C4; color:white"|'''Testfall-Name''' | |||
| align="center" style="background:#4472C4; color:white"|'''Anforderung''' | |||
| align="center" style="background:#4472C4; color:white"|'''Vorbedingungen und Eingänge''' | |||
| align="center" style="background:#4472C4; color:white"|'''Aktionen''' | |||
| align="center" style="background:#4472C4; color:white"|'''Erwartetes Ergebnis''' | |||
| align="center" style="background:#4472C4; color:white"|'''Ergebnis''' | |||
| align="center" style="background:#4472C4; color:white"|'''Bewertung''' | |||
| align="center" style="background:#4472C4; color:white"|'''Kommentar''' | |||
|- | |||
| 1 || Test auf alle Eingänge = 0 || 032 || Parameter geladen. || Simulation eines Zeitschritts || alle Ausgänge = 0 || alle Ausgänge = 0 || i. O. || Keine Anforderung im Lastenheft dazu gefunden | |||
|- | |||
| 2 || Test auf korrekte Formel, wenn Außentemperatur größer als Innentemperatur ist || 032 || Parameter geladen. ISO_TempDiff = <math>-1\,K</math> || Simulation eines Zeitschritts || ISO_Waermekapazitaet = <math>24{,}3432\,kWh</math> || ISO_Waermekapazitaet = <math>24{,}3432\,kWh</math> || i. O. || Keine Anforderung im Lastenheft dazu gefunden | |||
|- | |||
| 3 || Test auf korrekte Formel, wenn Innentemperatur größer als Außentemperatur ist || 032 || Parameter geladen. ISO_TempDiff = <math>-1\,K</math> || Simulation eines Zeitschritts || ISO_Waermekapazitaet = <math>-24{,}3432\,kWh</math> || ISO_Waermekapazitaet = <math>-24{,}3432\,kWh</math> || i. O. || Keine Anforderung im Lastenheft dazu gefunden -> macht eine negative Wärmekapazität Sinn? | |||
|- | |||
| 4 || Test auf Einheiten || 032 || || Einheiten der Signale ohne Simulation überprüft || Wärmekapazität laut Wikipedia in J/K || Wärmekapazität in J bzw. kWh, also eine Energie || n. i. O. || Ist der Ausgang dann eher eine Wärmeenergie? | |||
|} | |||
== Komponente IstTemp == | == Komponente IstTemp == | ||
'''Eingänge und Parameter''' | |||
* ISO_Waermekapazitaet = <math>25\,kWh</math> | |||
'''Ausgänge''' | |||
* ISO_IstTemp [<math>^{\circ}C</math>] | |||
{| class="wikitable" | |||
{| class="wikitable mw-collapsible mw-collapsed" | |||
| align="center" style="background:#4472C4; color:white"|'''Testfall-ID''' | |||
| align="center" style="background:#4472C4; color:white"|'''Testfall-Name''' | |||
| align="center" style="background:#4472C4; color:white"|'''Anforderung''' | |||
| align="center" style="background:#4472C4; color:white"|'''Vorbedingungen und Eingänge''' | |||
| align="center" style="background:#4472C4; color:white"|'''Aktionen''' | |||
| align="center" style="background:#4472C4; color:white"|'''Erwartetes Ergebnis''' | |||
| align="center" style="background:#4472C4; color:white"|'''Ergebnis''' | |||
| align="center" style="background:#4472C4; color:white"|'''Bewertung''' | |||
| align="center" style="background:#4472C4; color:white"|'''Kommentar''' | |||
|- | |||
| 1 || Test auf alle Eingänge = 0 || 048 || Parameter geladen. || Simulation eines Zeitschritts || alle Ausgänge = 0 || alle Ausgänge = 0 || i. O. || | |||
|- | |||
| 2 || Test auf verbrauchte Leistung größer als Wärmeverlust || 048, 050 || HZT_Verbrauchteleistung = 200; ISO_Waermeverlust = 100 || Simulation mehrere Zeitschritte, da Integratorblock genutzt wird || Ausgang ISO_IstTemp steigt an || Ausgang ISO_IstTemp steigt an || i. O. || Gain-Block muss angepasst werden. Temperatur steigt aktuell drastisch an auf unrealistische Werte. Evtl. begrenzen mit Saturationsblock. | |||
|- | |||
| 3 || Test auf verbrauchte Leistung größer als verbrauchte Leistung || 048, 050 || HZT_Verbrauchteleistung = 100; ISO_Waermeverlust = 200 || Simulation mehrere Zeitschritte, da Integratorblock genutzt wird || Ausgang ISO_IstTemp sinkt ab || Ausgang ISO_IstTemp sinkt ab || i. O. || Gain-Block muss angepasst werden. Temperatur steigt aktuell drastisch an auf unrealistische Werte. Evtl. begrenzen mit Saturationsblock. | |||
|} | |||
= Integrationstest = | = Integrationstest = |
Version vom 3. Juli 2023, 19:22 Uhr

Autoren: Daniel Gosedopp und Junjie Lyu
Betreuer: Prof. Dr.-Ing. M. Göbel
→ zum Hauptartikel: Systems Design Engineering - Seminaraufgabe SoSe 2023: Energiehaushalt eines Hauses
Einleitung
Im Rahmen der Lehrveranstaltung Systems Design Engineering im Studiengang Business and Systems Engineering soll ein Modell für den Energiehaushalt eines Hauses entwickelt werden. Dazu wird das Gesamtsystem in mehrere Module eingeteilt. Dieser Artikel beschreibt das Modul Energieerzeugung (kurz: EEZ).
Zielsetzung der Seminaraufgabe
Ziel der Seminaraufgabe ist es, das in der Softwareentwicklung weit verbreitete V-Modell (siehe Abb. 2) auf ein Beispielsystem (Haus) anzuwenden. Dabei wird der komplette Entwicklungsprozess von der Anforderungsdefinition über den Systementwurf und die Implementierung bis zum Testen der Software durchlaufen. In jeder Phase wird überprüft, ob die Arbeitsergebnisse den Spezifikationen aus der vorherigen Phase genügen. Dadurch wird die Effizienz und Qualität der Softwareentwicklung erhöht und die Fehleranfälligkeit minimiert. Durch die direkte Anwendung auf ein System sollen die Vorteile bei der Einhaltung des V-Modells "hands-on" von den Studierenden erfahren werden.

Anforderungsdefinition: Lastenheft
ID | Typ (I = Info, A = Anforderung) | Kapitel | Inhalt | Ersteller | Datum | Durchsicht von | am | Status Auftraggeber | Kommentar Auftraggeber | Status Auftragnehmer | Kommentar Auftragnehmer |
Funktionaler Systementwurf
Der funktionale Systementwurf war im Sommersemester 2023 kein Bestandteil der Gruppenaufgabe und wurde von Prof. Göbel, wie in Abb. 3 dargestellt, bereitgestellt.

Technischer Systementwurf
Das Modul EEZ wird in die vier Komponenten "Sonne", "Photovoltaik", "Warmwasserkollektor" und "Wechselrichter" untergliedert. Im technischen Systementwurf werden die Ein- und Ausgänge der Komponenten ebenfalls festgelegt. Die Komponente "Sonne" bekommt als Eingang die Simulationszeit und gibt eine Globalstrahlung aus. Diese dient als Eingang für die Komponenten "Photovoltaik" und "Warmwasserkollektor", welche eine elektrische Leistung (Gleichstrom) bzw. eine Warmwasserleistung ausgeben. Die Komponente "Wechselrichter" macht aus der Gleichstromleistung eine Wechselstromleistung. Abbildung 4 zeigt den Aufbau des Moduls EEZ.

Komponentenspezifikation
In der Komponentenspezifikation werden die Aufgabe, das Verhalten, der innere Aufbau und Schnittstellen für die im technischen Systementwurf erarbeiteten Komponenten definiert. Nachfolgend wird für jede der vier Komponenten eine Tabelle mit den jeweiligen Spezifikationen angelegt.
Sonne
Eingänge
- Simulationszeit [s]
Ausgänge
- EEZ_Globalstrahlung [W/m^2]
Parameter
- PAR_EEZ_StrahlungLippstadt [W/m^2]
ID | Kapitel | Inhalt | Ersteller | Datum |
Photovoltaik
Eingänge
- EEZ_Globalstrahlung [W/m^2]
Ausgänge
- EEZ_PVLeistung_DC [W]
Parameter
- PAR_EEZ_PVNeigungswinkel [°]
- PAR_EEZ_PVWirkungsgrad [-]
- PAR_EEZ_PVAnzahlModule [-]
- PAR_EEZ_PVModullaenge [m]
- PAR_EEZ_PVModulbreite [m]
ID | Kapitel | Inhalt | Ersteller | Datum |
Warmwasserkollektor
Eingänge
- EEZ_Globalstrahlung [W/m^2]
Ausgänge
- EEZ_Warmwasserleistung [W]
Parameter
- PAR_EEZ_WWKNeigungswinkel [°]
- PAR_EEZ_WWKAnzahlModule [-]
- PAR_EEZ_WWKModullaenge [m]
- PAR_EEZ_WWKModulbreite [m]
- PAR_EEZ_WWKTransmissionskoeffizient [-]
- PAR_EEZ_WWKAbsorptionskoeffizient [-]
ID | Kapitel | Inhalt | Ersteller | Datum |
Wechselrichter
Eingänge
- EEZ_PVLeistungDC [W]
Ausgänge
- EEZ_PVLeistungAC [W]
Parameter
- PAR_EEZ_DCACWirkungsgrad [-]
ID | Kapitel | Inhalt | Ersteller | Datum |
Programmierung / Modellierung
Nachdem die Vorbereitungen und Überlegungen in den Phasen des linken Asts des V-Modells abgeschlossen sind, kann die Programmierung bzw. Implementierung der einzelnen Komponenten erfolgen.
Komponente Sonne
In dieser Komponente wird aus dem aktuellen Monat und der aktuellen Stunde per "Lookup-Table"-Block eine Globalstrahlung in W/m^2 ermittelt. Für die aktuelle Stunde gilt:
Fehler beim Parsen (Syntaxfehler): {\displaystyle \text{Aktuelle Stunde} = floor\left(\frac{\text{Simulationszeit}\:mod\:86400}{3600}\right)}
Der aktuelle Monat wird aus einer weiteren "Lookup-Table" ermittelt. Dazu muss der aktuelle Tag im Jahr bekannt sein (also Tag 1 - Tag 365). Für den aktuellen Tag gilt:
Fehler beim Parsen (Syntaxfehler): {\displaystyle \text{Aktueller Monat} = floor\left(\frac{\text{Simulationszeit}}{24 \cdot 60 \cdot 60}\right)\:mod\:365 + 1}
Die Lookup-Table "Tag zu Monat" ermittelt schließlich den aktuellen Monat. Mit diesem und der aktuellen Stunde kann die Globalstrahlung ausgegeben werden.

Komponente Photovoltaik
Hier werden zunächst aus den Abmaßen und der Anzahl der PV-Module für die drei Ausrichtungen Süden, Westen und Osten die Flächen berechnet (also Fläche Süd, Fläche West und Fläche Ost). Diese können zusammen mit dem Wirkungsgrad der PV-Module und der Globalstrahlung zu drei Leistungen verrechnet werden. Außerdem wird ein weiterer Wirkungsgrad für die Ausrichtung und Neigung dazumultipliziert (für Details siehe Photovoltaikanlage: Ausrichtung und Neigung). Dieser wird für die drei Ausrichtungen S, W, O mithilfe des Neigungswinkels aus je einer Lookup-Table ermittelt. Durch Addition der drei Teilleistungen S, W, O ergibt sich die aktuell gelieferte Gesamtleistung der PV-Anlage (Gleichstrom).

Komponente Warmwasserkollektor
Auch hier wird zunächst aus den Abmaßen und der Anzahl der WWK-Module für die Ausrichtung Süden die Fläche berechnet. Mit einem sog. Transmissions- und einem Absorptionskoeffizienten sowie der Globalstrahlung kann durch Multiplikation mit der Fläche eine Warmwasserleistung berechnet werden (für Details siehe die Kollektorgleichung). Diese muss ebenfalls mit einem Wirkungsgrad für den Neigungswinkel, welcher aus einer Lookup-Table ausgelesen wird, multipliziert werden, um die tatsächliche Warmwasserleistung zu erhalten.

Komponente Wechselrichter
In dieser Komponente wird lediglich der Wirkungsgrad eines Wechselrichters berücksichtigt, sodass sich im Ergebnis eine Wechselstromleistung ergibt.

Komponententest
Den Komponententest haben wir für die Gruppe ISO durchgeführt. Die Testergebnisse werden für jede Komponente als Tabelle dargestellt.
Komponente Flächenberechnung
Eingänge und Parameter
- PAR_LKT_Grundflaeche =
- PAR_LKT_Deckenhoehe =
- PAR_LKT_Fensterflaeche =
Ausgänge
- ISO_FlaecheWand []
- ISO_FlaecheFenster []
Testfall-ID | Testfall-Name | Anforderung | Vorbedingungen und Eingänge | Aktionen | Erwartetes Ergebnis | Ergebnis | Bewertung | Kommentar |
Komponente Wärmeverlust
Eingänge und Parameter
- PAR_LKT_Grundflaeche =
- ISO_FlaecheWand =
- ISO_FlaecheFenster =
- PAR_LKT_Temperatur aus Temperaturdaten.mat
- ISO_IstTemp =
- PAR_ISO_WaendeUWert =
- PAR_ISO_DachUWert =
- PAR_ISO_FensterUWert =
- PAR_ISO_BodenUWert =
Ausgänge
- ISO_Waermeverlust [W]
- ISO_TempDiff [°C]
Testfall-ID | Testfall-Name | Anforderung | Vorbedingungen und Eingänge | Aktionen | Erwartetes Ergebnis | Ergebnis | Bewertung | Kommentar |
Komponente Wärmekapazität
Eingänge und Parameter
- ISO_FlaecheWand =
- ISO_Grundflaeche =
- ISO_FlaecheFenster =
- PAR_ISO_WaermekapazitaetWand =
- PAR_ISO_WaermekapazitaetDach =
- PAR_ISO_WaermekapazitaetFenster =
- PAR_ISO_WaermekapazitaetBoden =
Ausgänge
- ISO_Waermekapazitaet
Testfall-ID | Testfall-Name | Anforderung | Vorbedingungen und Eingänge | Aktionen | Erwartetes Ergebnis | Ergebnis | Bewertung | Kommentar |
Komponente IstTemp
Eingänge und Parameter
- ISO_Waermekapazitaet =
Ausgänge
- ISO_IstTemp []
Testfall-ID | Testfall-Name | Anforderung | Vorbedingungen und Eingänge | Aktionen | Erwartetes Ergebnis | Ergebnis | Bewertung | Kommentar |