Trwa remanent po dużej awarii portalu. Przepraszam za nieczytelną część treści do zakończenia remanentu !
W wielu zagadnieniach projektowania budowli występuje potrzeba przeprowadzenia obliczeń cieplno-wilgotnościowych przegród. Dotyczy to nie tylko budynków, co jest uregulowane w Rozporządzeniu [1], ale również budowli przemysłowych takich jak: zbiorniki do przechowywania wody i cieczy zamarzających oraz silosy do przechowywania cementu, wapna, gipsu itd. Ochrona może dotyczyć zarówno niskich lub wysokich temperatur jak i wykraplania się pary wodnej z powietrza. Szczególnie złożone są problemy pojawiąjące się wówczas, gdy przegroda pozostaje w styku z gruntem a także problem mostków cieplnych.
Przewodnictwo cieplne przez przegrody
Rozwiązanie klasycznego zagadnienia przewodnictwa ciepła
Rozwiążemy zagadnienie brzegowe jednowymiarowego, jednorodnego przepływu ciepła przez przegrodę wielowarstwową. Przegroda składa się z n warstw o grubościach di i współczynnikach przewodnictwa ciepła λi , ( i =1, …n). Zakładamy, że na zewnątrz przegrody stała temperatura jest zadana: po lewej stronie jest tw , a po prawej tz (rys.1). Zakładamy również, że można pominąć efekty odbicia na powierzchniach zewnętrznych, a powierzchnie kontaktu między warstwami przegrody są idealne (tzn. temperatura jest na nich ciągła).
Aby określić wymianę ciepła i rozkład temperatur pomiędzy powierzchniami przegrody rozwiązujemy układ równań:
div→q=0;→q=−→λ⋅gradt
Równanie (1 jest spełnione, gdy spełnione są warunki (2) i (3):
dqdx=0→q=const
q=−λdtdx→t−i−1=t+i−1;t−i=t+i;(i=1,2,..n)t−0=tw;t+n=tz
gdzie ” + ” oznacza wartość po prawej stronie powierzchni kontaktu, a ” − ” – po lewej. Tym samym rozkład temperatury w każdej warstwie jest liniowy.
Dla temperatur na powierzchniach kontaktowych otrzymujemy:
t1=tw–q⋅d1λ1t2=t1–q⋅d2λ1t1=tw–q⋅d1λ1….tn=tz=tn−1−q⋅dnλn
skąd wynika
tz=tw−q⋅n∑1diλi=n∑1Ri
czyli
q=U ⋅Δt
gdzie wprowadzono oznaczenia:
U=1R
R=n∑1Ri
Ri=diλi
Δt=tw–tz
Temperatura wokół przegrody
Temperatura wokół przegrody wg normy konstrukcyjnej [2]
Temperatura wokół budynku jest traktowana jak obciążenie klimatyczne, a jej wartość charakterystyczna powinna być przyjmowana zgodnie z normą PN-EN 1991-1-5 [2]. Rozróżnia się temperaturę maksymalną Tmax (latem) i minimalną Tmin (zimą). Temperatury te zostały ustalone z ekstremalnego rozkładu prawdopodobieństwa Gumbela z okresem powrotu 50 lat.
Na Rys.2. pokazano szacunkowe strefy temperatur charakterystycznych w Polsce, które zostały uzyskane z programu SPECBUD [3], który można pobrać również z tego panelu. Precyzyjne wartości dla miejscowości w Polsce na poziomie 0 m npm. można odczytać z programu.

Rys.2. Szacunkowe strefy temperatur charakterystycznych w Polsce: a) t_{min} , b) t_{max} wg normy [2]
W przypadku rozpatrywania miejsca położonego na innym niż zerowy poziomie – H metrów nad poziom morza – należy zastosować współczynniki korekcyjne:
t(H)={tmin−0,0035⋅H dla temperatury minimalnej tmax−0,0053⋅H dla temperatury maksymalnej
Temperatura powietrza na zewnątrz przegrody t__{out} zależy od tego , czy przegroda jest położona powyżej terenu, czy też pod nim.
W przypadku powierzchni zewnętrznych przegród wystawionych na działanie promieni słonecznych ( położonych powyżej powierzchni terenu) temperatura powietrza zewnętrznego latem zależy od współczynnika absorpcji ciepła przez powierzchnię, czyli koloru powierzchni wg zależności [2]– tab 5.2.
Dla powierzchni wystawionych na działanie słońca od strony północno-wschodniej (Pn-Wch) – południowo-zachodniej (Pd-Zach) dla powierzchni:
J – jasnej, lśniącej,
K – kolorowej lśniącej,
C- ciemnej,
mamy odpowiednio:
tout,nadziemne,lato=tmax+{0oC dla J : Pn-Wsch2oC dla K : Pn-Wsch4oC dla C : Pn-Wsch18oC dla J : Pd-Zach 30oC dla K : Pd-Zach42oC dla C : Pd-Zach
Od strony północnej osiągane są zwykle najmniejsze temperatury tout,nadziemne,lato,Pn= ok 50% tmax
Temperatura powietrza zewnętrznego zimą, przyjmuje się niezależną od obsorpcji promieni przez powierzchnię oraz od strony świata
tout,nadziemne,zima=tmin
Temperaturę powietrza wokół podziemnych części przegród przyjmuje się dla przegród zagłębionych o H poniżej poziomu terenu [2]– tab 5.3.:
tout,podziemne={−3oCdla H>1m−5oCdla H<1m+5oCdla H>1m+8oCdla H<1m
Temperaturę środowiska wewnętrznego (powietrza wewnątrz pomieszczeń) przy braku danych dokładniejszych przyjmuje się [2]– tab 5.1:
tin={20oC latem25oC zimą
Temperatura wokół przegrody wg normy ogrzewnictwa [4]
Do projektowania instalacji grzewczych zgodnie z normą PN-EN 12831 [4]
w budynkach przyjmuje się odmienną od pokazanej w poprzednim punkcie metodologię oraz wartości temperatur wokół przegród budowlanych. Na rys. 3 pokazano mapę klimatyczną Polski zgodnie z tą normą i temperatury projektowe (obliczeniowe) dla pięciu stref termicznych.

Rys.3. Szacunkowe strefy temperatur charakterystycznych w Polsce: a) t_{min} , b) t_{max} wg normy [2]
Od razu zauważamy ogromną różnice temperatur klimatycznych zgodnie z normą konstrukcyjną [2] i instalacji grzewczych [4] . Na przykład dla kilku miast w Polsce mamy:
Szczecin:
−16oC wg [4],
−29oC wg [4] :,
z rys.2a −28÷−30oC ; średnio −29oC ,
kt = 16/29 = 0,55;
Warszawa:
−20oC wg [4],
−29oC wg [4];
kt = 20/29 = 0,69;
Suwałki:
−24oC wg [4],
−34oC wg [2] ,
kt = 24/34 = 0,71;
Zauważmy, że mapa zawarta w normie ogrzewniczej [4] pochodzi z roku 2006 i odpowiada dokładnie podziałowi Polski przyjętemu w starej z roku 1982 „postradzieckiej” normie PN-82/B-02403 [5]. Zgodnie z normą [4] projektowa temperatura zewnętrzna jest najniższą wartością temperatury średniej z dwóch dni, zarejestrowanej dziesięć razy w okresie min 20 lat. Temperatura projektowa jest więc wyznaczona jako wartość szczytowa z prawdopodobieństwem przekroczenia p=2⋅10/20=0,1, co odpowiada okresowi powrotu T=1/0,1=10lat. W starej normie [5] nie ma niestety podanego okresu powrotu, więc nie jest jasne, czy podane tam projektowe temperatury były właściwe – przypuszczamy, że okres powrotu wynosił T=5lat , bo tak zakładały ówczesne normy klimatyczne.
Przeliczenie temperatur z okresu powrotu 10 lat (czyli prawdopodobieństwa przewyższenia p= 1/10= 0,1) na 50 lat (czyli p=1/50=0,02) można dokonać na podstawie współczynników stosunku Tmin,p/Tmin zgodnie z rys. NB.1 zamieszczonym w normie [2]. Stosunek Tmin,0,1/Tmin,0,02=0,78. Oznacza to, że w celu uzyskania temperatury obliczeniowej (projektowej) dla potrzeb ogrzewania z temperatury charakterystycznej powinien polegać na zmniejszeniu współczynnikiem kt=0,78
Sugerujemy, by do czasu korekty map klimatycznych zawartych w normie grzewczej [4] korzystać z normy konstrukcyjnej [2] po przemnożeniu odczytanej wartości przez współczynnik kt=0,78.
Na przykład dla Warszawy mamy tprojektowe,zima=−29⋅0,78=−23oC ( a nie −20oC).
W ten sam sposób można wyznaczyć temperatury projektowe latem, istotne przy projektowaniu klimatyzacji.
W tab. NB.2 normy [4] podano też projektowe temperatury w pomieszczeniach w budynkach tin . Zalecenia te przedstawiono w tab. 1 na podstawie rozporządzenia „Warunki techniczne” [6]), w którym powtórzono wymogi normowe. Wymogi te są wobec tego prawnie obowiązujące.
Tab.1. Projektowa temperatura wewnętrzna tin wg § 134 [6])
Współczynnika przenikania ciepła przez przegrody U
Ograniczenie współczynnika U dla budynków
Zgodnie z [1]– zał 2. ogranicza się od dołu współczynniki przenikalności ciepła przez przegrody wszelkich budynków z uwzględnieniem poprawek na pustki powietrzne, mostki cieplne w tym poprawki na łączniki mechaniczne zgodnie z tab.2. :
Tab.2. Maksymalne współczynniki przewodzenia ciepła UCmax [1]– zał 2
Obliczanie współczynnika przenikania ciepła U
Obliczenia współczynnika U należy przeprowadzić zgodnie z normą [7] w sposób wynikający z rozwiązania zagadnienia podstawowego przedstawionego wyżej.
Współczynnik przenikania ciepła Uc [W/(m2⋅K)] (7) jest odwrotnością oporu cieplnego całej przegrody (9) uzupełnionego o opory przejmowania i oddawania ciepła przez powierzchnie przegrody:
RT=Rsi+n∑1+Rse
gdzie:
Rsi opór przejmowania ciepła przez przegrodę, to znaczy opór powietrza otulającego przegrodę od strony wewnętrznej „i” (na rys 1 „w”)
Rse opór oddawania ciepła przez przegrodę, to znaczy opór powietrza otulającego przegrodę od strony zewnętrznej „e” (na rys.1 „z”)
Opory cieplne przejmowania ciepła przez przegrodę
Wartości oporów napływu i odpływu ciepła z przegród należy przyjmować zgodnie z tab.3:
Tab.3 Opory przejmowania ciepła Rs [(m2K)/W]
Oznaczenia (si) oraz (se) należy traktować jako oznaczenie typu powierzchni przegrody:
(si) – powierzchnia, przy której warstwa powietrza jest stabilna (nie jest omywana wiatrem lub przeciągiem)
(se) -powierzchnia, przy której warstwa powietrza pozostaje w ruchu (jest omywana wiatrem lub przeciągiem)
Oznacza to, że dla danej przegrody obie powierzchnie mogą być typu (si) – dotyczy w szczególności przegród położonych wewnątrz budowli, lub obie powierzchnie mogą być (se) dla przegród wystawionych na zewnątrz, tak, że obie mogą być omywane wiatrem.
Opór cieplny warstwy „i” przegrody o grubości di i obliczeniowym współczynniku przewodzenia ciepła λi wyznacza się z zależności (9).
Współczynnik przewodzenia ciepła materiałów budowlanych
Współczynniki λi występujące w (9) należy przyjmować na podstawie tablic w których zestawiono wyniki pomiarów dokonywane przez producentów oraz niezależne instytuty ( w Polsce ITB). Dotychczas zalecano tablice z normy PN-EN ISO 5946:1999 [8] – zał krajowy NC. Norma ta została zastąpiona normą pod tym samym numerem ale z roku 2008 [7], z której usunięto załącznik krajowy NC, W tab. 4 podajemy wyciąg z tego załącznika obejmujący najczęściej stosowne materiały Bardziej szczegółowe tablice współczynników przewodności podano w: normie [9], w podręcznikach fizyki budowli, np. [10] oraz badaniach producentów. Wartościowe są tablice opublikowane na stronie Tablice, zebrane przez Karolinę Kurtz. Tablice te są udostępnione również w tym artykule Tablice_materialowe.
Tab. 4. Średnie współczynniki przewodności wybranych materiałów budowlanych
λ w [(m⋅K)/W]
Opór cieplny pustek powietrznych
Często przegrody budowlane mają wbudowane pustki powietrzne. W zależności od pola powierzchni otworów łączących szczelinę powietrzną ze środowiskiem zewnętrznym mogą to być warstwy powietrza:
- niewentylowane ,
- słabo wentylowane
- wentylowane
W przypadku warstw powietrza niewentylowanych opór cieplny w zależności od grubości warstwy można przyjmować z tab.5. Dla warstw o innej grubości d można stosować interpolację
Tab. 5 Opór w [(m2K)/W] niewentylowanych warstw powietrza o wysokiej emisyjności powietrza
Dla warstw w pełni wentylowanych opór cieplny można pominąć, a opór przejmowania ciepła na powierzchni zewnętrznej przyjmuje się równy Rsi.
W sytuacji gdy pustka jest słabo wentylowana to całkowity opór warstwy wyznaczany jest z interpolacji liniowej między wartościami dla szczeliny słabo i dobrze wentylowanej z zależności:
RT=(1,5−Av)⋅RT,u+(Av−0,5)⋅RT,v
gdzie:
Av – powierzchnia otworów wentylujących w [cm2]
RT,u – całkowity opór cieplny z niewentylowaną warstwą powietrza
RT,v -całkowity opór cieplny bez uwzględnienia warstwy powietrza
Opór cieplny warstw niejednorodnych
Gdy przegroda składa się z warstw jednorodnych i niejednorodnych cieplnie, to jej całkowity opór wyznacza się jako średnią arytmetyczną:
RT=(RT′+R′′T)/2
z kresu górnego RT′ oporu cieplnego obliczonego wg [7]– pkt.6.2.3,
oraz kresu dolnego RT′′ oporu cieplnego obliczonego wg [7]– pkt.6.2.4.
Opór cieplny przestrzeni dachowych
Dach stromy z płaskim izolowanym stropem przestrzeni poddasza można uznać za warstwę jednorodną cieplnie o oporze podanym w tab.6 :
Tab.6. Opór cieplny w [(m2K)/W] przestrzeni dachowych
Mostki cieplne
Mostek termiczny (cieplny) powstaje w wyniku naruszenia ciągłości struktury wewnętrznej przegrody w związku z występowaniem materiałów budowlanych różniących się, najczęściej znacznie, wielkością współczynników przewodności cieplnej i może występować w kontekście zarówno geometrii miejsca występowania mostka, jak i charakter jego oddziaływania na przegrodę lub ewentualne złącza przegród.
Występowanie mostków termicznych jest często niedostrzegane przez projektantów, architektów i konstruktorów. Tymczasem jest to zjawisko, które w istotny sposób wpływa na parametry cieplne budynku, a tym samym na jego charakterystykę energetyczną.
Zagadnienie mostków cieplnych omawiają normy [11] oraz [12]. Oszacowanie wpływu mostka cieplnego na przenikalność cieplną przegrody powszechnie wykonuje się metodami numerycznymi , np z wykorzystaniem programu TRISCO [13].
W niniejszym artykule problem ten tylko sygnalizujemy. Jeśli nie są przeprowadzane dokładne obliczenia, to zalecamy, aby współczynnik przewodności każdej przegrody wyznaczony z równania (7) zwiększyć o 3% do 5%,
Podłogi na gruncie
Zagadnienie przenikalności ciepła z pomieszczeń przez posadzki do gruntu jest złożonym zagadnieniem Do obliczeń uproszczonych zwykle używa się metody [8]– wg załącznika krajowego NB. Niestety nie uzyskuje się ta metodą zadawalających wyników [14]. Obecnie powszechnie stosuje się nowoczesną metody opublikowaną w normie PN-EN ISO 13370.. Na użytek szacunkowych obliczeń w w ogrzewnictwie opracowano uproszczoną metodę PN-EN 12831.
Poniżej omawiamy wszystkie trzy metody.
Metoda PN-EN-ISO 6946:1999 (NB) [8]– zał krajowy NB
Uwaga: Metoda PN-EN ISO 6946 – zał NB nie jest aktualnie zalecana i wspierana.
Zgodnie z normą [8]– zał krajowy NB – podłogę stykającą się z gruntem dzieli się dwie strefy:
strefę I obwodową, o szerokości 1 m przy ścianach zewnętrznych (rys.4)
strefę II – wewnętrzną – pozostała powierzchnia podłogi.

Rys.4. Podłoga na gruncie. Podział na strefy. Podział podłogi na gruncie na strefy, Oznaczenia: te – obliczeniowa (projektowa) temperatura zewnętrzna, tg – obliczeniowa temperatura gruntu dla II strefy podłogi, tg=8 °C
Przy zagłębieniu górnej powierzchni podłogi więcej niż 1 m poniżej poziomu terenu cała podłogę traktuje się jako strefę II.
Wartość oporu cieplnego gruntu Rgr należy przyjmować:
w strefie I Rgr=0,5(m2K)/W
w strefie II wg tab. 7, przy czym wartość ta nie może przekraczać wartości granicznej
Rg,max=0,57⋅Z+0,09
gdzie Z- odległość w metrach poziomu górnego poziomu podłogi od zwierciadła wody gruntowej.:
Tab.7. Opór cieplny gruntu przylegającego do podłogi w strefie II [8]– Tab. NB.1
Opór cieplny przegrody posadzki na gruncie wyznacza się z zależności:
RT,gr=RT+Rgr
gdzie RT opór warstw posadzki wg zależności (2) ,
gdzie należy przyjąć: Rsi=0,17 , Rse=0.
Opór cieplny gruntu stykającego się ze pionową ścianą przyjmuje się wg tab 8.:
Tab.8. Opór cieplny gruntu przylegającego do ściany [8]– Tab. NB.1
Opór ściany wraz z gruntem oblicza się wg formuły (7), przy czym Rsi=0,13, Rse=0.
Metoda PN-EN ISO 13370: 2008 [15]
Pole temperatury w gruncie poza budynkiem zależy praktycznie nie zależy od gradientu geotermicznego w skorupie ziemskiej, a zależy od czynników atmosferycznych (temperatury powietrza zewnętrznego, nasłonecznienia, prędkości wiatru i pokrywy śnieżnej). W analizie przyjmuje się, że temperatura powietrza zewnętrznego, zmienia się według funkcji harmonicznej. z okresem niezależnym od głębokości. Amplituda wahań temperatury w gruncie zanika ekspotencjalnie (według funkcji wykładniczej) ze wzrostem głębokości. W rezultacie można wykazać (np. [14]), że na głębokości gruntu równej kilku metrom wahania dobowe i wszystkie dzienne krótkotrwałe zmiany temperatury można zaniedbać w porównaniu z wahaniami rocznymi. W praktyce na głębokości od 7 m do 10 m możemy zakładać stałą temperaturę gruntu równą wieloletniej średniej rocznej temperaturze powietrza (w Polsce środkowej przeciętnie 8 °C), podwyższonej o około 1,5 K przez promieniowanie słoneczne. Temperatura powierzchni będzie superpozycją średniej rocznej temperatury i funkcji harmonicznej o okresie 1 rok i amplitudzie rocznej.
Izotermy w gruncie przebiegają po liniach sferycznych, w sposób pokazany na rys. 5 na przykładzie podłogi o temperaturze 20oC ułożonej na gruncie według przykładu wzorcowego programu TRISCO [16].

Rys. 5. Przykład rozkładu temperatur w układzie posadzka -grunt-ściana (z programu [13]
Współczynnik przenikania cieplnego przez grunt
W normie [15] została przyjęta koncepcja, aby rzeczywisty kształt podłogi na gruncie zastąpić figurą o obwodzie i polu powierzchni A zgodnymi z danymi rzeczywistego rzutu podłogi. obwodowi budynku P przypisujemy pas podłogi o szerokości zwanej „wymiarem charakterystycznym”, zdefiniowanej wzorem
B′=AP/2
Wielkością charakteryzującą przenikanie ciepła przez grunt w stanie ustalonym jest stacjonarny współczynnik sprzężenia cieplnego definiowany jako iloraz strumienia cieplnego przez różnicę średniej w czasie temperatury powietrza wewnętrznego i zewnętrznego. W celu uproszczenia jego wyrażania wprowadza się pojęcie „grubości równoważnej”, to jest grubości gruntu o tym samym oporze cieplnym, jak zastępowany nią opór cieplny. Wyrażenie na „grubość równoważną” dt dla podłóg jest następujące:
dt=w+λ⋅(Rsi+Rf+Rse)
gdzie:
w – grubość całkowita ściany budynku (Rys.6)
λ – współczynnik przewodności cieplnej gruntu wg tab.9,
Rsi , Rse – współczynniki napływu wg tab.3,
Rf – opór cieplny płyty podłogi, łącznie z każdą warstwą izolacyjną, poniżej lub powyżej podłogi i każdym pokryciem podłogi.
Tab. 9 Właściwości cieplne gruntów [15]– tab.1
W kol.3 tab.9 podano szacunkowe wartości współczynnika przewodzenia ciepła przez grunty, a w kol. 4 iloczyn ρc, gdzie: : ρ – masa właściwa gruntu ( średnio 2000 kg/m3), c- ciepło właściwe niezamarzniętego gruntu ( średnio 1000 J/(kg K ) ).
Podłoga typu „płyta na gruncie” może być:
• nieizolowana,
• równomiernie izolowana na całej powierzchni (powyżej, poniżej lub wewnątrz płyty).
W przypadku płyty izolowanej izolacja może być:
a) pozioma – rys 7a,
b) pionowa – rys. 7b.
Ponadto rozróżnia się:
c) podłogi podniesione – rys. 7 c,
d) podłogi opuszczone poniżej poziomu gruntu – rys. 7d

Rys. 7: Cieplna izolacja krawędziowa budynku: a) pozioma, b) pionowa, c) podniesioną podłogą, d) opuszczeniem podłogi poniżej poziomu gruntu
Podłoga z izolacją krawędziową (pionową lub poziomą)
W przypadku podłóg bez izolacji krawędziowej (rys.6) współczynnik przenikania ciepła U wynosi
U=U0
a w przypadku podłóg z izolacją krawędziową
U=U0+2⋅ΨgeB′
gdzie:
U0={λ0,457⋅B′+dtdla podłóg dobrze izolowanych:dt≥B′2λπB′+dt⋅ln(πB′dt+1),dla podłóg nieizolowanych i słabo izolowanych:dt<B′
Liniowy współczynnik przenikania ciepła związany z połączeniami z izolacją krawędziową wyznacza się z równania:
Ψge=−λπ⋅[ln(nge⋅Ddt+1)–ln(nge⋅Ddt+d′+1)]
gdzie
nge={1,dla izolacji krawędziowej poziomej2,dla izolacji krawędziowej pionowej
W równaniach (25), (26):
D – szerokość pionowej lub poziomej izolacji krawędziowej zgodnie z rys 3a) i b),
d′ – dodatkowa grubość ekwiwalentna, wynikająca z izolacji krawędziowej wg zależności
d′=R′⋅λ
R′=Rn−dnλ
gdzie:
opór cieplny dodatkowej izolacji krawędziowej (lub fundamentu) w m2⋅K/W
Rn=dnλizolacja
dn – grubość izolacji krawędziowej (lub fundamentu) w m
Jeśli w danym detalu zastosowano kilka rodzajów izolacji krawędziowej (pionową lub poziomą, wewnętrzną lub zewnętrzną), to liniowy współczynnik przenikania ciepła Ψge,i (26) należy obliczyć dla każdej izolaacji (i=1, n) osobno , a do dalszych analiz (24) wybrać tą izolację, która daje największy redukcję przenikania ciepła, tzn:
Ψge=–nmaxi=1|Ψge,i|
Podłoga podniesiona
Podłoga podniesiona (rys. 7c) jest to podłoga zbudowana nad gruntem, a powstała przestrzeń podpodłogowa jest wentylowana grawitacyjnie (naturalnie) lub mechanicznie powietrzem zewnętrznym. W przypadki podłogi wentylowanej naturalnie przenikalność cieplną U z pomieszczenia do gruntu można wyznaczyć zależności
1U=1Uf+1Ug+Ux
gdzie:
Uf w [W/(m2K)] -współczynnik przenikania ciepła podłogi od środowiska wewnętrznego do podpodłogowej przestrzeni powietrznej. Wyznaczamy go zgodnie z normą [7] podług ogólnej zależności (6), biorąc pod uwagę warstwy podłogi.
Ug=U – współczynnik przenikania warstwy podpodłogowej. Wyznaczamy go z zależności (22) i (25) , ale:
w (22) zamieniamy Rf na Rg, a w (25) zamieniamy df na dg ,
gdzie: opór cieplny przestrzeni podpodłogowej zwykle Rg=0, gdy przestrzeń nie jest izolowana i nie traktuje się jej jako pustki powietrznej (jest z byt duża) i wobec tego z (25) otrzymujemy grubość ekwiwalentną przestrzeni podpodłogowej dg=w+λ(Rsi+Rse.
Zwykle głębokość podstawy przestrzeni podpodłogowej z (rys. 7d) nie przekracza 0,5 m i wówczas Ug wyznaczamy z zależności (22) (po zmianie df na dg). Jeśli jednak z>0,5, to stosujemy zależność
dt=w+λ⋅(Rsi+Rf+Rse)
Przenikalności od wentylacji przestrzeni podpodłogowej Uxobliczamy z zaleźnosci
Ux=2⋅h⋅UwB′+1450⋅ε⋅ν⋅fwB′
Metoda PN-EN 12831:2006 [4]
Ze względu na to, że szereg projektantów postrzegało trudności obliczeniowe przy stosowaniu normy [15], więc na ich użytek (projektowania instalacji grzewczych w budynkach) opracowano uproszczoną metodę wyznaczania współczynnika przenikania ciepła przez podłogę na gruncie oraz przez ścianę pomieszczenia zagłębionego w gruncie wg zasad podanych w normie [7].
W celu obliczenia wartości U wyznacza się opór cieplny poszczególnych warstw z uwzględnieniem oporu przejmowania ciepła od strony wewnętrznej budynku Rsi, a pomijając opór przejmowania ciepła od strony gruntu, czyli przyjmując Rse=0. Na podstawie wyliczonych wartości B’ oraz U określa się wartości Uequiv,bf przy użyciu tab.9 lub nomogramów podanych w normie [15]. Korzysta się z wartości, które są funkcją wielkości zagłębienia w terenie oraz wartości B′ (21) i Upodłogi w rozpatrywanym budynku. Dla wartości pośrednich Uequiv,bf należy wyznaczyć stosując metodę interpolacji liniowej.
Tab.10 Wartości Uequiv,bf podłogi ogrzewanego podziemia w funkcji zagłębienia poniżej poziomu terenu. i współczynnika przenikania ciepła podłogi i wartości B′ – na podstawie normy PN-EN 12832:2002 [17]
W dobie powszechnego stosowania arkuszy kalkulacyjnych powyższa, uproszczona metoda jest coraz rzadziej stosowana.
Porównanie metod obliczania strat ciepła przez grunt
W pracy [18] przeprowadzono analizę dokładności metod obliczania strat ciepła do gruntu za pomocą:
A trójwymiarowej symulacji numerycznej zgodnie z normą PN-EN ISO 10211 [12],
B metody przybliżonej podług normy PN-EN ISO 13370 [15] i symulacją 2D mostków wg [12],
C metody orientacyjnej wg normy normy PN-EN ISO 13370 [15] i symulacją mostków zgodnie z normą [11].
Przeprowadzono analizę podłogi na gruncie (rys 8 a) oraz podziemia nieogrzewanego (rys 8b). W obu analizowanych sytuacjach projektowych stwierdzono, że metoda szacunkowa C jest niedokładna (daje znacznie zawyżone wartości), a metoda B( przybliżona) jest mierna.
Rozkład temperatur w przegrodzie
Obliczanie rozkładu temperatury w przegrodzie powinno służyć do sprawdzenia wilgotnościowego przegród jako zadanie pomocnicze.. Wyznacza się go z zależności klasycznej (4). Strumień cieplny q (6) zależy od różnicy ustalonych temperatur w pomieszczeniach po obu stronach przegrody (10) i jest najczęściej znany z wymogów normowych .
W przypadku nieogrzewanych pomieszczeń temperatura tw nie jest znana. Wówczas mamy do czynienia z procesem czasowym, związanym z bezwładnością cieplną obiektu, a nie tylko pomieszczenia.
Sprawdzenie kondensacji powierzchniowej
Temperatura wewnętrzna powierzchni przegrody bez mostków cieplnych liniowych
Temperaturę wewnętrzną ϑw na powierzchni przegrody bez mostków cieplnych liniowych wyznacza się ze wzoru
ϑw=tw–Uc⋅(tw−tz)⋅Rsw
gdzie:
Rsw=0,167m2⋅K/W – opór przejmowania ciepła przyjmowany przy sprawdzaniu minimalnej temperatury wewnętrznej powierzchni przegród nieprzezroczystych
Uc
Ciśnienie cząstkowe pary wodnej w pomieszczeniu
Ciśnienie cząstkowe pary wodnej w pomieszczeniu pw zależy od wilgotności względnej φw oraz ciśnienia cząstkowego pary wodnej nasyconej dla temperatury obliczeniowej w tym pomieszczeniu pn.
Ciśnienie cząstkowe pary wodnej nasyconej w powietrzu pn w temperaturze tw odczytuje się z tab.11. Na przykład dla temperatury w pomieszczeniu tw=4oC odczytujemy pn=8,13hPa
Ciśnienie cząstkowe pary wodnej wylicza się z zależności
pw=φw⋅pn
gdzie φw w ułamku dziesiętnym ( np dla 45% wstawiamy 0,45).
Jeśli z założeń technologicznych lub bilansu wilgoci nie wynikają inne wartości, to wilgotność względna w pomieszczeniach można przyjmować następująco:
φw=55 dla pomieszczeń mieszkalnych (w tym kuchnie, łazienki oraz ustępy), pokoje chorych w szpitalach, sanatoriach, pokoje dziecięce żłobkach, i przedszkolach
φw=45 dla pomieszczeń pozostałych (użyteczności publicznej i produkcyjnych) w których nie wydziela się para wodna z otwartych zbiorników lub urządzeń technologicznych i nie stosuje się nawilżania,
Dla obliczonego ciśnienia cząsteczkowego pary pw ponownie z tab.11 odczytujemy temperaturę powietrza nasyconego ts.
Na przykład dla pw=8,13⋅0,45=3,65hPa mamy temperaturę nasycenia ts=−6,1oC
Warunek ryzyka wykraplania (kondensacji) pary wodnej
Warunek przy którym nie nastąpi wykraplanie pary wodnej na powierzchni przegrody jest następujący:
ϑw≥ts
gdzie: ϑw (35) , ts – wyżej.
Tab.11 Ciśnienie cząstkowe pary wodnej nasyconej w powietrzu pn w funkcji temperatury [8]– Tab. NA.3
Przykłady rachunkowe
Podłoga na gruncie
Dla podłogi ułożonej na gruncie (rys.6) z warstwami jak na Rys.8 wyznaczyć opór cieplny , współczynnik przenikania ciepła oraz rozkład temperatur.
Grubość ścian wynosi w=0,25\,m $.
Na rys.9 pokazano plan podłogi (w rzucie z góry).
Pole powierzchni podłogi A=10⋅9–3⋅6=72m2
Obwód podłogi P=10+9+4+3+6+6=38m
Szerokość charakterystyczna:
(21 \toB^\prime =\cfrac{72}{38/2}=3,789 \,m$
Płyta jest ułożona na gruncie kat.1 (glina lub ił).
tab.9 → współczynnik przewodzenia ciepła przez grunt wynosi
λ=1,5W/(m⋅K)
Opór cieplny posadzki i współczynnik przewodności cieplnej
W tab. 12 zestawiono wyniki obliczeń oporu cieplnego i przewodności cieplnej posadzki dla danych rodzajów, grubości i przewodności cieplnej warstw (kol. 1, 2 i 3).
Tab.12. Zestawienie oporów warstw posadzki z przykładu
Opory napływu ciepła na przegrodę Rsi=0,17m^2\cdot K/ W$ przyjęto z tab.3 jak dla przepływu ciepła w dół przy powierzchni napływu z warstwą powietrza w bezruchu. (stałego). Natomiast opór wychodzenia ciepła z przegrody Rse=0 jak dla przegrody ułożonej na gruncie.
Opór cieplny warstw w kol ostatniej wyliczano ze wzoru (9).
Na przykład dla warstwy drugiej – betonu, obliczono jak następuje R2=0,031,7=0,018m2⋅K/W
Po zsumowaniu oporu cieplnego poszczególnych warstw wg wzoru (8) otrzymano opór przegrody: R=0,170+0,018+1,500+0,029+0,009+0,222+0=1,948m2⋅K/W
Współczynnik przenikalności cieplnej jest odwrotnością oporu cieplnego zgodnie z formułą (7) U=1/1,948=0,513W/(m^2 \cdot K$.
W rzeczywistości współczynnik przewodzenia ciepła przez przegrodę będzie większy ze względu na mostki cieplne w przegrodzie, wywołane przez przestrzenny układ i połączenia (w tym noraża ścian) z innymi przegrodami, łączniki, niejednorodność warstw, spoiny w murach, i inne. Nie przeprowadzamy dokładnych obliczeń, ale powiększyliśmy go o 3% jak dla przegrody dobrej jakości zgodnie z uwagą zamieszoną w punkcie Mostki cieplne.
Rozkład temperatur w przegrodzie
Rozkład temperatur w przegrodzie został nakreślony czerwona linią na rys 8.
Uzyskano go zw wyrażeń (4) dla strumienia cieplnego wyznaczonego ze współczynnika przenikalności cieplnej Unetto=0,513W/m2⋅K.(bez rozrzucania przenikalności wynikającej z mostków na poszczególne warstwy).
Temperaturę powietrza wewnątrz pomieszczenia przyjęto tw=20oC jak dla pomieszczeń przeznaczonych na pobyt ludzi bez okryć zewnętrznych wg tab.1.
Temperaturę powietrza pod przegrodą (od strony gruntu) przyjęto tz=+4oC na podstawie indywidulanych pomiarów minimalnej temperatury pod posadzką zimą.
Różnica temperatur po obu stronach przegrody (10) wynosi Δt=20−4=16oC
Strumień ciepła (6) wynosi q=0,513⋅16=8,2W/m2
Temperatura we wnętrzu pomieszczenia 20oC,
Temperatura na powierzchni wewnętrznej przegrody 20−0,17⋅8,2=18,6oC,
Temperatura na granicy betonu i styropianu 18,6−0,018⋅8,2=18,5oC,
Temperatura na granicy styropianu i betonu 18,5−1,5⋅8,2=6,1oC.
Zestawienie obliczeń podano w tab.13
Tab.13. Obliczenia strumienia cieplnego i rozkładu temperatur w przegrodzie
Przenikanie ciepła przez podłogę z gruntem
Maksymalny współczynnik przewodzenia ciepła UCmax dla temperatury wewnątrz pomieszczenia tw≥16oC→
$tab 2, poz. 6 :U_{max} = 0,3 \, W/(m^2 \cdot K)$
Wariant przy braku izolacji
przy założeniu: opór cieplny podłogi Rf=0W/(m2⋅K) wg tab.
Współczynnik napływu przy przepływie ciepła w dół wg tab.3, Rsi=0,17(m2K)/W
Współczynnik odpływu przy przepływie ciepła w dół wg tab.3, Rse=0,04(m2K)/W
Współczynnik przewodności cieplnej gruntuλ=1,5W/(m⋅K)
Grubość ścian zewnętrznych (Rys.6) w=0,25m
Grubość równoważna posadzki z gruntem
(22) →dt=0,25+1,5⋅(0,17+0+0,04)=0,565m
Ponieważ dt=0,565<B′=3,789m→ podłoga jest nieizolowana (lub słabo izolowana) i:
(25): U0=2⋅1,5π⋅3,789+0,565⋅ln(π⋅3,7890,565+1)=0,744W/(m2⋅K)
(23): U=0,744 \, W/(m^2 \cdot K)$
Ponieważ U=0,744>Umax=0,3W/(m2⋅K)→ podłogę należy ocieplić.
Wariant z zastosowaniem ław z betonu komórkowego γ=500kg/m3 (rys.10)
Rozważmy wariant posadowienia ścian na ławie wykonanej z lekkiego betonu – rys. 10. Sytuację oceniamy procedurą analogicznie do pionowej izolacji krawędziowej dla głębokości
D=0,6m i grubości dn=0,4m.
Tab. 4: dla betonu komórkowego o gęstości γ=500kg/m3 współczynnik przewodzenia ciepła
λizolacja=0,25W/(m⋅K)
Opór cieplny dodatkowej izolacji krawędziowej (fundamentu):
(30) →Rn=0,40,25=1,600(m2K)/W,
(29) →R′=1,600–0,41,5=1,333(m2K)/W,
Dodatkowa grubość ekwiwalentna, wynikająca z izolacji krawędziowej (fundamentu)
(28) →d′=1,333⋅1,5=2,000m,
Liniowy współczynnik przenikania ciepła związany z izolacją krawędziową ( w tym przypadku lekkiego fundamentu):
(26) →
Ψge=–1,5π⋅[ln(2⋅0,60,565+1)–ln(2⋅0,60,565+2,000+1)]=0,477⋅(ln(3,123)–ln(1,468)=–0,361W/(mK),
Współczynnik przenikalności:
(24) →U=0,744+2⋅(−0,361)3,789=0,554W/(m2K).
Ponieważ U=0,554>Umax=0,3W/(m2K),
to podłoga wymaga innego docieplenia.
Wariant z izolacją krawędziową – przeciwmrozową (Rys. 11)
Zastosowano:
- Izolację pionową: grubość dn=0,08m ; długość D=0,6m ; przewodność λizol=0,05W/(m⋅K)
Opór cieplny izolacji :
(30) →Rn=0,080,05=1,600m2K)/W
(29 →R′=1,600–0,081,5=1,547m2K)/,W
Grubość ekwiwalentna pionowej izolacji krawędziowej:
(28) →d′=1,547⋅1,5=2,320m,
Liniowy współczynnik przenikania ciepła:
(26) →Ψge=–1,5π⋅[ln(2⋅0,50,565+1)–ln(2⋅0,50,565+2,320+1)]=–0,344W/(mK) - Izolację poziomą :grubość dn=0,08m ; długość D=0,7m ; przewodność λizol=0,05W/(m⋅K)
Opór cieplny izolacji :
(30) →Rn=0,080,05=1,600m2K)/W
(29 →R′=1,600–0,081,5=1,547m2K)/,W
Grubość ekwiwalentna pionowej izolacji krawędziowej:
(28) →d′=1,547⋅1,5=2,320m,
Liniowy współczynnik przenikania ciepła:
(26) →Ψge=–1,5π⋅[ln(2⋅0,70,565+1)–ln(1⋅0,70,565+2,320+1)]=–0,281W/(mK)
(31) →Ψge=–max|−0,344|;|−0,281|=−0,344
czyli bardziej efektywna jest izolacja krawędziowa pionowa. Izolacji krawędziowej poziomej można nie wykonywać. Należy natomiast wykonać dodatkowe izolacje , wynikające z ochrony przed mrozem w narożach oraz innych mostków cieplnych.
Współczynnik przenikalności:
(24) →U=0,744+2⋅(−0,344)3,789=0,563W/(m2K).
Ponieważ U=0,563>Umax=0,3W/(m2K),
to podłoga wymaga innego docieplenia.
Wariant przy izolacji posadzki 60 mm
Opór cieplny podłogi Rf=1,948W/(m2⋅K) wg tab. 12
Współczynnik napływu przy przepływie ciepła w dół wg tab.3, Rsi=0,17(m2K)/W
Współczynnik odpływu przy przepływie ciepła w dół wg tab.3, Rse=0,04(m2K)/W
Współczynnik przewodności cieplnej gruntuλ=1,5W/(m⋅K)
Grubość ścian zewnętrznych (Rys.6) w=0,25m
Grubość równoważna posadzki z gruntem
(22) →dt=0,25+1,5⋅(0,17+1,948+0,04)=3,487m
Ponieważ dt=3,487<B′=3,789m→ podłoga jest słabo izolowana:
(25): U0=2⋅1,5π⋅3,789+0,565⋅ln(π⋅3,7890,565+1)=0,2894W/(m2⋅K)
(23): U=0,744 \, W/(m^2 \cdot K)$
Ponieważ U=0,2894≤Umax=0,3W/(m2⋅K)→ podłogę ocieplono prawidłowo
Literatura
- Rozporządzenie Ministra Infrastruktury z dnia 12 kwietnia 2002 r. w sprawie warunków technicznych, jakim powinny odpowiadać budynki i ich usytuowanie, ( Dz.U. 2002 nr 75, poz.69, tekst jednolity Dz.U. 2019 poz. 1065 (z późn zmianami)
- PN-EN 1991-1-5, Eurokod 1: Oddziaływania na konstrukcje – Część 1-5: Oddziaływania ogólne – Oddziaływania termiczne
- SPECBUD, Mapa stref oddziaływań klimatycznych v.1.2, [ https://www.computerworld.pl/ftp/pobierz/pc/109084.html ]
- PN-EN 12831:2006, Instalacje grzewcze w budynkach, Metoda obliczania projektowego obciążenia cieplnego
- PN-82/B-02403, Ogrzewnictwo – Temperatury obliczeniowe zewnętrzne
- Rozporządzenie Ministra Infrastruktury z dnia 12 kwietnia 2002 r. w sprawie warunków technicznych, jakim powinny odpowiadać budynki i ich usytuowanie, ( Dz.U. 2002 nr 75, poz.69, tekst jednolity Dz.U. 2019 poz. 1065 (z późn zmianami
- PN-EN ISO 6946: 2008, Komponenty budowlane i elementy budynku. Opór cieplny i współczynnik przenikania ciepła. Metoda obliczania
- PN-EN ISO 6946:1999, Komponenty budowlane i elementy budynku. Opór cieplny i współczynnik przenikania ciepła. Metoda obliczania
- PN-EN ISO 12524:2003, Materiały i wyroby budowlane – Właściwości cieplno-wilgotnościowe. Stabelaryzowane wartości obliczeniowe
- Pogorzelski J., A., (2003), Przewodnik po PN-EN ochrony cieplnej budynków, Wydawnictwa ITB, Warszawa
- PN-EN ISO 14683: 2007, Mostki cieplne w budynkach – Liniowy współczynnik przenikania ciepła – Metody uproszczone i wartości orientacyjne
- PN-EN ISO 10211: 2008 „Mostki cieplne w budynkach. Strumień cieplny i temperatura powierzchni. Obliczenia szczegółowe
- Physibel, TRISCO Software, Belgia, [https://www.physibel.be/en/products/trisco]
- Pogorzelski J., A., (2002), Straty ciepła z budynku przez grunt według PN-EN ISO 13370:2001 Wydawnictwa ITB, kwartalnik nr 3(123), Warszawa
- PN-EN ISO 13370: 2008, Cieplne właściwości użytkowe budynków. Przenoszenia ciepła przez grunt. Metody obliczania
- Physibel, TRISCO Software, Belgia, https://www.physibel.be/en/products/trisco]
- Robakiewicz M., (2007), Obliczanie strat ciepła przez grunt na podstawie nowych norm, Energia i budynek 4 (4), s.29-31
- Dylla A., Pawłowski K., Rożenek P., Analiza metody obliczania strat ciepła do gruntu z wykorzystaniem normy PN-EN ISO 14683, w: Wybrane problemy budownictwa,(Red. Sobczak-Piąstka J.), Wydawnictwa Uczelniane Uniwersytetu Technologiczno-Przyrodniczego w Bydgoszczy, Bydgoszcz 2015
________________________________