­
Analiza cieplno-wilgotnościowa przegród ⋆ Chodor-Projekt ⋆ Architekci i Inżynierowie. Encyklopedia PiWiki
Loading [MathJax]/jax/output/HTML-CSS/jax.js
A B D E F G H I K Ł M N O P R S T U W Z
Al An Ar

Analiza cieplno-wilgotnościowa przegród

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.

Spis treści ukryj

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ń:

divq=0;q=λgradt

Przepływ ciepła przez przegrodę

Rys.1. Schemat przewodzenia ciepła przez przegrodę

Równanie (1 jest spełnione, gdy spełnione są warunki (2) i  (3):

dqdx=0q=const

q=λdtdxti1=t+i1;ti=t+i;(i=1,2,..n)t0=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=twqd1λ1t2=t1qd2λ1t1=twqd1λ1.tn=tz=tn1qdnλn

skąd wynika

tz=twqn1diλi=n1Ri

czyli

q=U Δt

gdzie wprowadzono oznaczenia:

U=1R

R=n1Ri

Ri=diλi

Δt=twtz

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)={tmin0,0035H dla temperatury minimalnej tmax0,0053H 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>1m5oCdla 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.

Strefy temperatur w Polsce

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=210/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=290,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ł 2Przewodności graniczne ciepłą przez przegrody

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/(m2K)] (7) jest odwrotnością oporu cieplnego całej przegrody  (9) uzupełnionego o opory przejmowania i oddawania ciepła przez powierzchnie przegrody:

RT=Rsi+n1+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]Opory napływu

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 [(mK)/W]  Współczynniki przewodności

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,5Av)RT,u+(Av0,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+RT)/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

Opory dachów

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,57Z+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 gruntu II strefaOpó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.

Rys. 6. Schemat interakcji cieplnej budynku i gruntu

Tab. 9 Właściwości cieplne gruntów [15]– tab.1

Właściwości cieplne gruntów

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,457B+dtdla podłóg dobrze izolowanych:dtB2λπB+dtln(πBdt+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(ngeDdt+1)ln(ngeDdt+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=Rndnλ

gdzie:
opór cieplny dodatkowej izolacji krawędziowej (lub fundamentu) w  m2K/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=2hUwB+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=twUc(twtz)Rsw

gdzie:
Rsw=0,167m2K/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=φwpn

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,130,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:

ϑwts

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.

Rys.8  Przykład: warstwy podłogi

Grubość ścian wynosi w=0,25\,m $.

Na rys.9 pokazano plan podłogi (w rzucie z góry).
Pole powierzchni podłogi A=10936=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$

Rys.9 Kształt podłogi z przykładu

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/(mK)

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,018m2K/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,948m2K/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/m2K.(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=204=16oC
Strumień ciepła  (6) wynosi q=0,51316=8,2W/m2

Temperatura we wnętrzu pomieszczenia   20oC,
Temperatura na powierzchni wewnętrznej przegrody  200,178,2=18,6oC,
Temperatura na granicy betonu i styropianu 18,60,0188,2=18,5oC,
Temperatura na granicy styropianu i betonu 18,51,58,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 tw16oC
$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/(m2K) 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/(mK)
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=21,5π3,789+0,565ln(π3,7890,565+1)=0,744W/(m2K)
(23): U=0,744 \, W/(m^2 \cdot K)$

Ponieważ U=0,744>Umax=0,3W/(m2K) 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/(mK)

Rys.10. Posadowienie na ławie z lekkiego betonu

Opór cieplny dodatkowej izolacji krawędziowej (fundamentu):
(30) Rn=0,40,25=1,600(m2K)/W,
(29) R=1,6000,41,5=1,333(m2K)/W,
Dodatkowa grubość ekwiwalentna, wynikająca z izolacji krawędziowej (fundamentu)
(28d=1,3331,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(20,60,565+1)ln(20,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/(mK)
    Opór cieplny izolacji :
    (30) Rn=0,080,05=1,600m2K)/W
    (29 R=1,6000,081,5=1,547m2K)/,W
    Grubość ekwiwalentna  pionowej izolacji krawędziowej:
    (28d=1,5471,5=2,320m,
    Liniowy współczynnik przenikania ciepła:
    (26Ψge=1,5π[ln(20,50,565+1)ln(20,50,565+2,320+1)]=0,344W/(mK)
  • Izolację poziomą :grubość  dn=0,08m ;  długość D=0,7m  ; przewodność λizol=0,05W/(mK)
    Opór cieplny izolacji :
    (30) Rn=0,080,05=1,600m2K)/W
    (29 R=1,6000,081,5=1,547m2K)/,W
    Grubość ekwiwalentna  pionowej izolacji krawędziowej:
    (28d=1,5471,5=2,320m,
    Liniowy współczynnik przenikania ciepła:
    (26Ψge=1,5π[ln(20,70,565+1)ln(10,70,565+2,320+1)]=0,281W/(mK)

Rys. 11 Przykład izolacji krawędziowej

(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/(m2K) 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/(mK)
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=21,5π3,789+0,565ln(π3,7890,565+1)=0,2894W/(m2K)
(23): U=0,744 \, W/(m^2 \cdot K)$

Ponieważ U=0,2894Umax=0,3W/(m2K) podłogę  ocieplono prawidłowo

 Literatura

  1. 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)
  2. PN-EN 1991-1-5, Eurokod 1: Oddziaływania na konstrukcje – Część 1-5: Oddziaływania ogólne – Oddziaływania termiczne
  3. SPECBUD, Mapa stref oddziaływań klimatycznych v.1.2, [ https://www.computerworld.pl/ftp/pobierz/pc/109084.html ]
  4. PN-EN 12831:2006, Instalacje grzewcze w budynkach, Metoda obliczania projektowego obciążenia cieplnego
  5. PN-82/B-02403, Ogrzewnictwo – Temperatury obliczeniowe zewnętrzne
  6. 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
  7. PN-EN ISO 6946: 2008, Komponenty budowlane i elementy budynku. Opór cieplny i współczynnik przenikania ciepła. Metoda obliczania
  8. PN-EN ISO 6946:1999, Komponenty budowlane i elementy budynku. Opór cieplny i współczynnik przenikania ciepła. Metoda obliczania
  9. PN-EN ISO 12524:2003, Materiały i wyroby budowlane – Właściwości cieplno-wilgotnościowe. Stabelaryzowane wartości obliczeniowe
  10. Pogorzelski J., A., (2003), Przewodnik po PN-EN ochrony cieplnej budynków, Wydawnictwa ITB, Warszawa
  11. PN-EN ISO 14683: 2007,  Mostki cieplne w budynkach – Liniowy współczynnik przenikania ciepła –  Metody uproszczone i wartości orientacyjne
  12. PN-EN ISO 10211: 2008 „Mostki cieplne w budynkach. Strumień cieplny i temperatura powierzchni. Obliczenia szczegółowe
  13. Physibel, TRISCO Software, Belgia, [https://www.physibel.be/en/products/trisco]
  14. 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
  15. PN-EN ISO 13370: 2008, Cieplne właściwości użytkowe budynków. Przenoszenia ciepła przez grunt. Metody obliczania
  16. Physibel, TRISCO Software, Belgia, https://www.physibel.be/en/products/trisco]
  17. Robakiewicz M., (2007), Obliczanie strat ciepła przez grunt na podstawie nowych norm, Energia i budynek 4 (4), s.29-31
  18. 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

________________________________

Comments : 0
O autorze
* dr inż. Leszek Chodor. Architekt i Inżynier Konstruktor; Rzeczoznawca budowlany. Autor wielu projektów budowli, w tym nagrodzonych w konkursach krajowych i zagranicznych, a między innymi: projektu wykonawczego konstrukcji budynku głównego Centrum "Manufaktura" w Łodzi, projektu budowlanego konstrukcji budynku PSE w Konstancinie Bielawa, projektów konstrukcji "Cersanit" ( Starachowice, Wałbrzych, Nowograd Wołyński-Ukraina), projektu konstrukcji hali widowiskowo-sportowej Arena Szczecin Autor kilkudziesięciu prac naukowych z zakresu teorii konstrukcji budowlanych, architektury oraz platformy BIM w projektowaniu.

Wyślij

Translate »