Skip to content

Note

Click here to download the full example code

Bauphysik

Formeln für Bauphysikalische Berechnungen

1 Randbedingungen (Klima)

Link enbau

A1.1 gleitenden Mittelwerts der Aussentemperatur

Die Komfort-Anforderungen an die empfundene Temperatur werden in der Norm SIA 180:2014 [1.4] neu in Abhängigkeit des gleitenden Mittelwerts der Aussentemperatur θrm über 48 Stunden gemäss folgender Formel festgelegt.

\(\theta_{rm}(t) = \frac{ 1 }{ N } * \sum_{j=0}^{N-1} \theta_e(t-j)\)

  • gleitender Mittelwert der Aussentemperatur fiir die Stunde t
  • Aussenlufttemperatur der Stunde t-j
  • N Anzahl in den Mittelwert einbezogener Stunden N=48

A1.2 maximale Luftgeschwindigkeit

\(V_{max} = 0.54 * (\frac{g* \Delta T}{T_i})^{0.5} * H^{0.5}\)

\(V_{max} = k * \sqrt{\Delta \theta * H}\)

  • Temperaturdifferenz zwischen Raumluft und kalter Oberflache
  • Höhe der vertikalen Flache
  • Regressionskoeffizient I

A1.1

# A1.2
def max_air_flow_speed(delta_temp, air_temperature, height):
    """Vmax = m*s^-1"""
    GRAV = 9.81
    return 0.54 * (GRAV * delta_temp / air_temperature) ** 0.5 * height**0.5

print(max_air_flow_speed.__doc__)
print(max_air_flow_speed(10, 20, 1))

Out:

Vmax = m*s^-1
1.1959506678789056

2 Wärme

Link Enbau Wärme

wärme_strahlung

wärme_leitung

wärme_konvektion

docs/generated/gallery/images/wärme_konvektion.excalidraw.png

A2.1 Wärmestromdichte q [\(W/m^{2}\)]

\(q = h_c * (\theta_s - \theta_a) + h_r * (\theta_s - \theta_r)\)

  • \(q\) Warmestromdichte \(W/m^{2}\)
  • \(h_c\) Warmeiibergangskoeffizient Konvektion \(W/(m^{2}*K)\)
  • \(h_r\) Warmeiibergangskoeffizient Strahlung \(W/(m^{2}*K)\)
  • \(\theta_s\) Oberflaichentemperatur \(°C\)
  • \(\theta_a\) Lufttemperatur \(°C\)
  • \(\theta_r\) Strahlungstemperatur \(°C\)

Die beiden Wärmeübergangkoeffizienten Konvektion und Strahlung können zu einem kombinierten Wärmeübergangskoeffizienten \(h = h_c + h_r\) zusammengefasst werden. Die Formel zur Bestimmung der Wärmestromdichte einer Oberfläche zu seiner Umgebung vereinfacht sich durch die Verwendung der Umgebungstemperatur θn zur Formel A2.2:

A2.2 Wärmestromdichte mittels Umgebungstemperatur

\(q= h * (\theta_s - \theta_n)\)

  • \(h\) kombinierter Warmetibergangskoeffizient \(W/(m^{2}*K)\)
  • \(\theta_n\) Umgebungstemperatur \(°C\)

Die Umgebungstemperatur θn wird gemäss Formel A2.3 ermittelt und kann als eine Gewichtung der Luft- und Strahlungstemperatur betrachtet werden.

A2.3 Umgebungstemperatur θn

\(\theta_n = \frac{h_c * \theta_a + h_r * \theta_r}{h_c + h_r} = F_c * \theta_a + (1- F_c) * \theta_r\)

  • \(F_c\) = konvektiver Anteil
  • \(F_c\) = \(h_c/(h_c+h_r)\)

A2.4 Umgebungstemperatur aussen, der Aussentemperatur θe:

\(\theta_e = \frac{h_c * \theta_{a,e} + h_{r,e} * \theta_{r,e}}{h_{c,e} + h_{r,e}}\)

A2.5 Umgebungstemperatur innen, der Raumtemperatur θi:

\(\theta_i = \frac{h_{c,i} * \theta_{a,i} + h_{r,i} * \theta_{r,i}}{h_{c,i} + h_{r,i}}\)

A2.6 Wärmeübergangswiderstände \(R_{se}\)

\(R_{se} = \frac{1}{h_c} = \frac{1}{h_{r,e} + h_{c,e}}\)

\(R_{se} = \frac{1}{\varepsilon * 4 * \sigma * T^{3}_{me} + (a + b * v)}\)

A2.7 Wärmeübergangswiderstände \(R_{si}\)

\(R_{si} = \frac{1}{h_i} = \frac{1}{h_{r,i} + h_{c,i}}\)

\(R_{si} = \frac{1}{\varepsilon * 4 * \sigma * T^{3}_{mi} + c * \sqrt[3]{\Delta \theta }}\)

  • \(R_{se}\) Warmeiibergangswiderstand aussen \((m^{2}*K)/W\)
  • \(R_{si}\) Warmeiibergangswiderstand innen \((m^{2}*K)/W\)
  • \(\Delta \theta\) Temperaturdifferenz \(\Delta \theta = \theta_{si} - \theta_{a,i}\) \(K\)
  • \(T_{me}\) Mitteltemperatur aussen \((T_{r,e}+T_{se})/2\) \(K\)
  • \(T_{mi}\) Mitteltemperatur innen \((T_{r,e} + T_{si}) /2\) \(K\)
  • \(a\) Regressionskoeffizient a = 4.0 \(W/(m^{2}*K)\)
  • \(b\) Regressionskoeffizient b = 4.0 \((W*s)/(m^{3}*K)\)
  • \(v\) Windgeschwindigkeit \(m/s\)
  • \(c\) Regressionskoeffizient c = 1.31 (bei Warmestr6mung horizontal) \(W/(m^{2}*K^{4/3})\)

2.2 Wärmespeicherung

A2.8 statische, flächenbezogene Wärmespeicherfähigkeit

\(k_{stat} = \sum d_j*\rho_j*c_j\)

  • \(d_j\) Dicke der Schicht j
  • \(\rho_j\) Rohdichte der Schicht j
  • \(c_j\) spezifische Warmekapazitat der Schicht j

A2.9 speicherwirksame Dicke \(d_{T,max}\)

\(d_{T,max} = \sqrt{\frac{a*T}{2*\pi}} = \sqrt{\frac{\lambda * T}{\rho*c*2*\pi}}\)

  • \(a\) Temperaturleitfahigkeit \(m^{2}*s^{-1}\)
  • \(T\) Periodenlange \(s\)
  • \(\lambda\) Warmeleitfahigkeit \(W*(m*K)^{-1}\)
  • \(\rho\) Rohdichte \(kg*m^{-3}\)
  • \(c\) spezifische Warmekapazitat \(kJ*(kg*K)^{-1}\)

2.3 Wärmebrücken

A2.10 thermische Gesamtleitwert \(L_{i,j}\)

\(L_{3D,i,j} = \sum^{N_k}_{k=1}U_{k,(i,j)} * A_k + \sum^{N_m}_{m=1} \psi_{m(i,j)}*l_m + \sum^{N_n}_{n=1} \chi_{n(i,j)}\)

  • \(L_{3D,i,j}\) thermischer Leitwert aus einer 3D-Berechnung \(W*K^{-1}\)
  • \(U_{k,(i,j)}\) Warmedurchgangskoeffizient des 1D-Bauteils \(W*(m^{2}*K)^{-1}\)
  • \(A_k\) Flache, über die der \(U_k\)-Wert gilt \(m^{2}\)
  • \(\psi_{m(i,j)}\) längenbezogener Warmedurchgangskoeffizient \(W\)
  • \(l_m\) Länge, über die der \(U_k\)-Wert gilt \(m\)
  • \(X_{n,(i,j)}\) punktbezogener Warmedurchgangskoeffizient \(W*K^{-1}\)
  • \(N_k\) Anzahl der Warmedurchgangskoeffizienten \(-\)
  • \(N_m\) Anzahl der langenbezogenen Warmedurchgangskoeffizienten \(-\)
  • \(N_n\) Anzahl der punktbezogenen Warmedurchgangskoeffizienten \(-\)

A2.11 Wärmebrückenzuschläge 2D

\(\psi = L_{2D} - \sum^{N_k}_{k=1}U_k*l_k\)

  • \(L_{2D}\) thermischer Leitwert aus einer 2D-Berechnung \(W*(m*K)^{-1}\)
  • \(U_k\) Warmedurchgangskoeffizient des 1D-Bauteils \(W*(m^{2}*K)^{-1}\)
  • \(l_k\) Lange, über die der \(U_k\)-Wert gilt \(m\)

A2.12 Wärmebrückenzuschläge 3D

\(\chi = L_{3D} - \sum^{N_k}_{k=1}U_k*A_k - \sum^{N_m}_{m=1} \psi_m*l_m\)

  • \(L_{3D}\) thermischer Leitwert aus einer 3D-Berechnung \(W*K^{-1}\)
  • \(U_k\) Warmedurchgangskoeffizient des 1D-Bauteils \(W*(m^{2}*K)^{-1}\)
  • \(A_k\) Flache, über die der \(U_k\)-Wert gilt \(m^{2}\)
  • \(\psi_m\) langenbezogener Warme- durchgangskoeffizient \(W*(m*K)^{-1}\)
  • \(l_m\) L$nge, über die der \(U_k\)-Wert gilt \(m\)

A2.13 Mittlerer U-Wert

\(U_m = \frac{L_{}i,j}{A_{tot}}\)

  • \(L_{i,j}\) Gesamtleitwert aller betrachteten Gebaudehiillenbereiche \(W*K^{-1}\)
  • \(A_{tot}\) Gesamtfliche der betrachteten Gebiudehiillenbereiche (Projektionsflaiche) \(m^{2}\)

3 Feuchte

Enbau Formelbuch

3.1 Oberflächentemperatur

Die innere Oberflächentemperatur \(\theta_{si}\) einer Wand wird nach Formel A3.1 berechnet:

\(\theta_{si} = \theta_i - (R_{si} * U) * (\theta_i - \theta_e)\)

  • \(\theta_e\) Aussentemperatur \(°C\)
  • \(\theta_i\) Innentemperatur \(°C\)
  • \(R_{si}\) Warmeiibergangswiderstand innen \(m^{2}*K*W^{-1}\)
  • \(U\) Warmedurchgangskoeffizient der Wand \(W*(m^{2}*K)^{-1}\)
"""
Rsi:
windows, doors  = 0.13 [m^2*K/W]
other           = 0.25 [m^2*K/W]
"""


def surface_temperature(temp_out, temp_in, thermal_res_in, heat_transfer_coeff):
    return temp_in - (thermal_res_in * heat_transfer_coeff) * (temp_in - temp_out)


def surface_temperature_from_factor(
    temp_out,
    temp_in,
    surface_temp_factor,
):
    return temp_out + surface_temp_factor * (temp_in - temp_out)


def surface_temperature_unheated(g1, temp_1, g2, temp_2, g3, temp_3):
    """g1 + g2 + g3 = 1"""
    return (g1 * temp_1) + (g2 * temp_2) + (g3 * temp_3)


(
    surface_temperature(-10, 21, 0.13, 0.2),
    surface_temperature_from_factor(-10, 21, 0.974),
)

Out:

(20.194, 20.194)

Oberflächentemperaturfaktor

Oberflächentemperaturfaktor fRsi. Er ist als Verhältnis zwischen der Differenz der inneren Oberflächentemperatur eines Aussenbauteils und der Lufttemperatur aussen sowie der Differenz zwischen den Lufttemperaturen innen und aussen, bei vorgegebenem innerem Wärmeübergangswiderstand Rsi, definiert (siehe Formel A3.2).

Für flächige Bauteile (eindimensionale Wärmestromsituation) gilt: fRsi = 1 – (Rsi . U)

def surface_temperature_factor(surface_temperature, temp_out, temp_in):
    return (surface_temperature - temp_out) / (temp_in - temp_out)


def simple_surface_temperature_factor(thermal_res_in, heat_transfer_coeff):
    return 1 - (thermal_res_in * heat_transfer_coeff)


(
    surface_temperature_factor(20.194, -10, 21),
    simple_surface_temperature_factor(0.13, 0.2),
)

Out:

(0.974, 0.974)

Feuchtebelastung

Die Feuchtebelastung in einem Raum ist abhängig von der Feuchteproduktion G, der Aussenluftfeuchte νe und dem Aussenluft-Volumenstrom qv. Der raumseitige Feuchteüberschuss ∆ν kann gemäss Formel A3.5 bestimmt werden.

In Tabelle A3.2 ist die Feuchteproduktion G einer Person in Abhängigkeit zu seiner Tätigkeit angegeben.

(A3.5)

\(\Delta v = v_i - v_e = \frac{G}{q_v}\)

  • \(\Delta v\) raumseitiger Feuchtetiberschuss \(g*m^{-3}\)
  • \(v_i\) volumenbezogene Raumluftfeuchte \(g*m^{-3}\)
  • \(v_e\) volumenbezogene Aussenluftfeuchte \(g*m^{-3}\)
  • \(G\) Feuchteproduktion im Raum \(g*h^{-1}\)
  • \(q_v\) Aussenluft- Volumenstrom \(m^{3}/h\)
Tatigkeit Feuchteproduktion G in \(g/h\)
ruhig liegend, schlafend 45
ruhig sitzend 60
sitzende Tatigkeit (Biro, Schule, Labor), ruhig stehend 70
leichte Tatigkeit, stehend (Laden, Werkbankarbeit, Labor) 95
mittelschwere Tatigkeit, stehend (Haushalt, Werkstatt) 115

Tabelle: SIA 180:2014[3.2]

A3.5

def excess_humidity_room(indoor_hum, outdoor_hum):
    return indoor_hum - outdoor_hum


def excess_humidity_room_from_flow(hum_production_person, outdoor_air_flow):
    return hum_production_person / outdoor_air_flow


# A3.6
def water_vapour_overpressure(excess_humidity_room, temp_in):
    RV = 462  # Pa*m3 * (kg*K)^-1
    temp_in += 273.15
    return excess_humidity_room * RV * temp_in


# A3.7
def water_vapour_pressure_in(water_vapour_overpressure, water_vapour_pressure_out):
    return water_vapour_overpressure + water_vapour_pressure_out

3.6 Wasserdampfüberdruck

Der raumseitige Wasserdampfüberdruck ∆p im Raum lässt sich somit nach Formel A3.6, der raumseitige Wasserdampfdruck pv,i nach Formel A3.7 ermitteln:

(A3.6)

\(\Delta p = \Delta v * R_v * T_i\)

(A3.7)

\(P_{v,i} = P_{v,e} + \Delta p\)

  • \(p_{v,i}\) Wasserdampfdruck innen \(Pa\)
  • \(p_{v,e}\) Wasserdampfdruck aussen \(Pa\)
  • \(\Delta_{P}\) raumseitiger Wasserdampfüberdruck \(Pa\)
  • \(R_v\) spez. Gaskonstante für Wasserdampf = 462 \(Pa * m^{3} * (kg*K)^{-1}\)
  • \(T_i\) absolute Temperatur innen \(K\)

Für maritime Klimate wird im Anhang A.2 die Luftfeuchtelast in fünf Luftfeuchteklassen gemäss Tabelle A3.3 eingeteilt und jeder Klasse in Abbildung A3.4 ein Wert für den raumseitigen Feuchteüberschuss ∆ν bzw. den raumseitigen Wasserdampfüberdruck ∆p in Abhängigkeit des monatlichen Mittelwerts der Aussenlufttemperatur zugeordnet.

Raumseitige Luftfeuchteklassen in Abhängigkeit der Aussenlufttemperatur

Abbildung A3.4: Raumseitige Luftfeuchteklassen in Abhängigkeit der Aussenlufttemperatur (EN ISO 13788 [3.1])

Raumseitige Randbedingungen nach SIA 180 [3.2]

Um Feuchteschäden zu vermeiden, darf die Feuchte in Räumen mit Personenbelegung die in Tabelle A3.4 angegebenen Werte im Tagesmittel in Abhängigkeit zur Aussenlufttemperatur nicht übersteigen. Die Angaben zur relativen Luftfeuchte beziehen sich auf eine Raumlufttemperatur von 20 °C und ein fRsi von 0,70. Für abweichende Raum- und Aussenlufttemperaturen ist die maximal zulässige relative Feuchte in Abbildung A3.5 angegeben, berechnet mit Formel A3.8. Die Aussenluft-Volumenströme sind so zu wählen, dass diese Grenzen nicht überschritten werden.

Maximal zulässige relative Feuchte der Raumluft

Abbildung A3.5: Maximal zulässige relative Feuchte der Raumluft (Tagesmittelwerte) (Norm SIA 180:2014 [3.2])

Bei abweichenden Nutzungsbedingungen (Raumlufttemperaturen ≠ 20 °C) und in Räumen mit unvermeidbaren Wärmebrücken mit einem Oberflächentemperaturfaktor unter 0,70 ist eine Berechnung der maximal zulässigen relativen Raumluftfeuchte φi,max mit der Gleichung A3.8 notwendig:

(A3.8)

4 Luftströmungen

Link Enbau

A4.1 Rechtecköffnung

\(q_{v,Rechtecköffnung} = c_d * H * B * \frac{1}{3}*\sqrt{g*H*\frac{(T_i-T_e)}{T_e}}\)

Die empirische Grösse cd ist insbesondere von der Fensterposition und -geometrie abhängig [4.2] und kann – falls keine Resultate aus Experimenten oder CFD-Simulationen vorliegen – in grober Näherung mit etwa 0.6 angenommen werden.

A4.2 Kippfenstern

\(q_{v,Kippflügel}(\alpha) = c_k(\alpha) * q_{v,Rechtecköffnung}\)

5 Lichttechnik

5.1 Lichtstrom Tagsehen

\(\Phi_v = Km * \int_{380nm}^{780nm} \Phi_{e\lambda}(\lambda) * V(\lambda) d\lambda\)

5.2 Lichstrom Nachtsehen

\(\Phi'_v = K'm * \int_{380nm}^{780nm} \Phi_{e\lambda}(\lambda) * V'(\lambda) d\lambda\)

5.3 Skalierungsfaktoren Km

\(K_m \approx 683 lm/W\)

\(K'_m \approx 1700 lm/W\)

5.4 Lichtstärke

\(\Phi_v(\varphi_1) = lv(\varphi_1)* \Omega\)

5.5 Lichtstärke

\(l_v(\varphi_1) = \frac{d\Phi_v(\varphi_1)}{d\Omega}\)

Iv (φ1) wird als Lichtstärke bezeichnet. Die physiologische Einheit Candela (cd) der Lichtstärke Iv ist zugleich eine Basiseinheit des SI-Einheitensystems. Der Gesamtlichtstrom Φv, tot berechnet sich aus:

5.6 Gesamtlichtstrom Φv, to

\(\Phi_{v,tot} = \int_{4\pi} l_v (\varphi_1) * d\Omega\)

5.7 Leuchtdichte

\(L_{v0} = \frac{l_v(\varphi = 0)}{A}\)

\(L_{v(\varphi1)} = \frac{l_v (\varphi_1)}{A*\cos \varphi_1} * \frac{1}{\cos \varphi_1}\)

\(L_{v(\varphi1)} = \frac{dl_v}{dA} * \frac{1}{\cos \varphi_1}\)

Leuchtdiode Lichstrom, lm
Leuchdiode 10 \(^{-2}\)
LED-Lampe 220 V, 8 W 480
Gluhlampe 220 V, 60 W 730
Gluhlampe 220 V, 100 W 1380
Leuchtstoffrohre 220 V, 40 W 2300
Quecksilberdampflampe 220 V, 125W 5400
Lichtbogen 250 W \(\sim\) 10'000
Quecksilberdampflampe 220 V, 2000 W \(\sim\) 125'000
Xe-Hochstdrucklampe 100 kW 3*10 \(^{6}\)

Tabelle 5.1: Gesamtlichtstrom einiger Lichtquellen

5.8 Leuchtdichte

\(cd*m^{-2}\) \(cd*m^{-2}\)
Nachthimmel 10 \(^{-3}\) Glimmlampe 2+10 \(^2\) - 1 * 10 \(^3\)
graue Wolken 10 \(^{3}\) Leuchtstofflampe 3*10 \(^3\) - 7 * 10 \(^3\)
blauer Himmel 3*10 \(^{3}\) Glühlampe 3*10 \(^4\) - 3 * 10 \(^5\)
weisse Wolken 10 \(^{4}\) Bogenlampe 1.8 * 10 \(^8\)
Vollmond 2.5 * 10 \(^{3}\) Hg-Höchstdrucklampe 6 - 10 \(^8\)
Sonne 1.5 * 10 \(^{9}\) Xe-Höchstdrucklampe 1 * 10 \(^{10}\)

Tabelle 5.2: Leuchdichten ausgewählter Lichtquellen

5.9 Beleuchtungsstärke

\(E_v = \frac{\Phi_v}{A_E}\)

\(E_v = \frac{d\Phi_v}{dA_E}\)

TODO: Bild

\(E_v = l_v(\varphi_1) * \frac{\cos \varphi_2}{r^2} * \Omega_0\)

Beluchtungsstärke zeitpunkt lx
Sonnenlicht Juni
Dezember
bis 100'000
9'000
Tageslicht bei bedecktem Himmel Juni
Dezember
4'000 - 20'000
900 - 2'000
Vollmond 0.25
Sterne ohne Mond, klare Nacht 10 \(^{-3}\)
Grenze der Farbwahrnehmung 3
Arbeitsplatzbeleuchtung, hohe Ansprüche 1'000
Wohnzimmerbeleuchtung 120

Tabelle: Beleuchtungsstärke

Gegenüberstellung von physikalischen und photometrischen Strahlungsgrössen

6 Energie/Leistung

Link Enbau

A6.1 Klimakorrektur mit akkumulierten Temperaturdifferenzen (ATD)

6.2 Heizenergiebedarf

\(E_{H,an} = \frac{\theta_{\sum,an}}{\theta_{\sum,per}}*E_{H,per}\)

6.3 Heizenergiebedarf

\(E_{H,1} = \frac{\theta_{\sum,per1}}{\theta_{\sum,per2}}*E_{H,2}\)

A6.4 Statische Wärmekapazität \(C_{stat}\)

\(C_{stat} = \sum A_i * \sum d_j * \rho_j * c_j\)

  • \(A_i\) Flache des Bautells \(m^{2}\)
  • \(d_j\) Dicke der Schicht j \(m\)
  • \(\psi_j\) Rohdichte der Schicht j \(kg*m^{3}\)
  • \(c_j\) spezifische Warmekapazitat der Schicht \(kJ*(kg*K)^{-1}\)

A6.5 Dynamische Wärmekapazität \(C_{stat}\)

6.5 Kühlungsbedarf

A6.6 Kühlungsbedarf

\(C_r = \sum A_i*k_i\)

  • \(A_i\) Flache des Bauteils i \(m^{2}\)
  • \(K_i\) flachenbezogene Warmespeicherfahigkeit des Bauteils (mit \(R_{si}\) und \(R_{se}\)),
    berechnet nach EN ISO 13786 [6.7] \(kJ * (m^{2}* K)^{-1}\)

Anhang

Link Enbau

9.2 maximal zulässigen U-Wertes

9.3 maximal zulässigen U-Wertes bei max zul. Luftgeschwindigkeit

9.4 Sonnenschutz Gesamtenergiedurchlassgrad des Fensters

Total running time of the script: ( 0 minutes 0.001 seconds)

Download Python source code: plot_1_Bauphysik.py

Download Jupyter notebook: plot_1_Bauphysik.ipynb

Gallery generated by mkdocs-gallery