Měření vlhkosti ve stavebních hmotách pomocí mikrovln (III)

V minulých číslech [8], [9] byla popsána experimentální aparatura na měření vlhkosti stavebních hmot na základě propustnosti mikrovlnného záření 2,4 GHz a ukázán experimentální výstup z měření pro nelehčený cihelný střep a porobeton. V tomto pokračování jsou výsledky měření použity pro součinitele kapilární vodivosti a jsou diskutovány výhody stanovení vlhkosti pomocí mikrovln.

Úvod

Měřící aparatura ke sledování transportu kapalné vlhkosti v inertních pórovitých hmotách pomocí změny intenzity elektromagnetického mikrovlnného záření byla sestavena experimentálně na základě fyzikálních teoretických předpokladů, které byly popsány v předchozích kapitolách. Byly tedy aplikovány známé teoretické předpoklady o proudění vody, popisující dynamiku transportu vlhkosti v kombinaci s využitím interakce mikrovln s vlhkým prostředím. V daném případě se jedná o využití poznatků o rychlosti proudění popisované Darcyho rovnicí a plnění pórů vyjádřené rovnicí kontinuity, dosazením do Lykovova vztahu jako popis hustoty toku pro výpočet součinitele kapilární vodivosti pomocí Matanovy metody na základě zjištění rozložení vlhkosti v v nestacionárním stavu průběhu navlhání. Pro použití elektromagnetického mikrovlnného záření v daném případě byly uplatněny poznatky Lambertova - Beerova zákona, jako matematické vyjádření závislosti absorpce elektromagnetického záření od vlastností materiálu, souvztažnost mezi absorbancí, optickou délkou a koncentrací z pohlcující substance.

V dalším pokračování uvádíme výsledky prvotních výstupů měření na vzorcích materiálu z keramického páleného střepu a z pórobetonu, získané jako prvotní naměřené hodnoty, jako vstupní údaje pro stanovení navlhacích křivek a následně výpočet součinitele kapilární vodivosti.

6. Metodika měření na sestavené aparatuře

Na základě testovacích měření na vzorcích porézních materiálů s odlišnou objemovou hmotností byla sestavena metodika postupu měření na sestavené měřící aparatuře.


Obr. 14: Schema znázorňující postup při měření na experimentálně sestavené měřící aparatuře [2], kde je součinitel kapilární vodivosti; um hmotnostní vlhkost a x je délková souřadnice od zdroje navlhání, viz obr. 1

Na obr. 14 je uvedeno schema postupu při získávání podkladů k výpočtu součinitele kapilární vodivosti jak Matanovou metodou, tak i metodou integrální. Uvedeným postupem je možné získat průběžně bez přerušení měření několik navlhacích křivek v nestacionárním stavu ve zvolených časových intervalech pro inertní porézní hmoty bez jejich destrukce.

Jako vstupní údaje o materiálu je třeba zjistit jejich hmotnost, relativní vlhkost v ustáleném stavu a funkční závislost změny intenzity elektromagnetického mikrovlnného záření na jejich hmotnostní vlhkosti.

Navlhavost vyjadřuje časový průběh navlhání a slouží pro orientační stanovení časových intervalů měření před dosažením celkového provlhnutí materiálu, viz kap. II.

Mikrovlnná propustnost materiálu v oblasti 2,4 GHz je vyjádřena jako pokles detekovaných hodnot napětí, které vrací mikrovlnné čidla před a po projití EMVZ vzorkem - viz část II. Všechny tyto vstupní informace jsou zjišťovány pomocí sestavené aparatury.

Vlastní měření sledování transportu vlhkosti je prováděno pomocí napojení měřící aparatury na PC, který umožňuje grafické vyjádření prvotních výstupů měření potřebné pro následné matematické zpracování a výpočet součinitele kapilární vodivosti.

Stanovené souřadnice polohy profilu vlhkostního čela je možné v závěru měření ověřit gravimetrickou metodou.

Výpočet součinitele kapilární vodivosti vychází ze stanovených navlhacích křivek pomocí vzorce (10) – viz část I.

7. Podklady pro určování součinitele vlhkostní vodivosti

Pro výpočet součinitele vlhkostní vodivosti potřebujeme nejdříve získat křivky navlhání – tj. funkce vyjadřující závislost mezi vlhkostí a vzdáleností od zdroje vlhkosti. Postup pro určení křivek navlhání je následující.

7.1. Stanovení závislosti mezi množstvím záření, které prochází vzorkem a obsahem vlhkosti ve vzorku

Tuto závislost jsme určili použitím gravimetrické metody – gravimetrickou metodou určíme hmotnostní vlhkost vzorků a zároveň vždy určíme množství záření, které prochází vzorkem. Naměřené grafy funkčních závislostí změny intenzity elektromagnetického mikrovlnného záření na hmotnostní vlhkosti pro keramický pálený střep s objemovou hmotností cca 1 800 kg.m-3 (cihelna Štíty na Moravě) a pórobeton s objemovou hmotností cca 520 kg.m-3 (výrobce Hrušovany), byly uvedeny v předchozích částech I. a II.

Z hodnot naměřených pro 6 vzorků potřebujeme stanovit závislost vlhkosti u na množství záření, které projde vzorkem. K vyjádření této závislosti lze využít program Maple, který metodou nejmenších čtverců vypočítá rovnici funkce:

Tab. 4: Naměřené hodnoty (závislost hmotnostní vlhkost - změna intenzity EMWZ) [2]

vzorek 1
EMWZ (mV)386094197265336383424435440
Hmotn. vlhkost (%)18.8516.0014.649.266.033.682.601.731.130.4
vzorek 2
EMWZ (mV)3345103143195250273301343425
Hmotn. vlhkost (%)17.6015.0513.159.546.874.773.813.062.351.4
vzorek 3
EMWZ (mV)325980124187230279301360420
Hmotn. vlhkost (%)18.3515.9311.5810.387.184.803.652.952.181.4
vzorek 4
EMWZ (mV)367088144197251285331435
Hmotn. vlhkost (%)17.4815.0613.199.607.084.903.862.420.7
vzorek 5
EMWZ (mV)323866125186224260290341445
Hmotn. vlhkost (%)18.4815.7313.659.936.884.663.522.802.051.4
vzorek 6
EMWZ (mV)364281134191255296318345439
Hmotn. vlhkost (%)17.6214.9413.069.526.954.933.993.342.541.4

Z hodnot naměřených jako funkční závislost změny intenzity EMVZ na hmotnostní vlhkosti pro 6 vzorků byla metodou nejmenších čtverců v programu Maple stanovena rovnice funkce vyjadřující závislost vlhkosti um na intenzitě EMW záření z, které projde vzorkem (podle rovnice regrese vložené křivky funkční závislosti):

um = -1,342033167·10-7 · z3 + 0,0001936510773 · z2 - 0,1038753765 · z + 20,78641097

kde z je veličina udaná v mV, která vyjadřuje intenzitu elmg. mikrovlnného záření po průchodu vzorkem, přičemž

z0 = 500 mV odpovídá intenzitě záření, které vstupuje do vzorku.
Obr. 15: Vyjádření funkční závislosti změny intenzity elektromagnetického mikrovlnného záření na hmotnostní vlhkosti pro keramický pálený střep v programu Maple [2].

7.2. Určení závislosti mezi množstvím mikrovlnného záření, které prochází vzorkem a vzdáleností od zdroje vlhkosti

Hodnoty souřadnice polohy profilu vlhkostního čela získáme měřením na mikrovlnné měřící aparatuře. Ve zvolených časových intervalech od počátku navlhání vzorku materiálu připraveného ke sledování transportu vlhkosti (viz kap.II. odst. 4.3) je prováděna detekce množství vlhkosti, pronikající do materiálu. Na obr. 15 je vyobrazen vzorek materiálu keramického páleného střepu v průběhu navlhání.

Obr. 16: Vzorek materiálu keramický pálený střep v průběhu navlhání a po ukončení experimentu po 30ti minutách [2].

Na obr. 16 jsou grafické výstupy měření detekce koncentrace vlhkosti, vyjádřené jako závislost změny intenzity elektromagnetického záření v čase pojezdu vlnovodů v délce sledovaného vzorku materiálu keramického páleného střepu.

Rychlost pojezdu vlnovodů je konstantní a jejím přepočtem na délkové údaje v programu Linregrese Excel jsou stanoveny souřadnice polohy pohlceného elektromagnetického mikrovlnného záření v závislosti na obsažené hmotnostní vlhkosti viz obr. 17.

Obr. 17 (3 grafy): Prvotní naměřené hodnoty pro stanovení polohy profilu vlhkostního čela při sledování navlhání vzorku materiálu keramického páleného střepu [2].

Obr. 18 (3 grafy): Převedení časových údajů podle rychlosti pojezdu vlnovodů na délkové údaje souřadnice X polohy profilu vlhkostního čela pro keramický pálený střep v časových intervalech po 10 minutách v programu Linregrese Excel [2].

Obr. 19: Vzorek materiálu z pórobetonu v průběhu navlhání [2].

7.3. Sledování polohy vlhkosti v pórobetonu

Pro srovnání průběhu navlhání jsou uvedeny také výsledky měření postupu profilu vlhkostního čela u vzorku materiálu z pórobetonu.

Obr. 20: Během 10 minut příjem vody navlháním +7,60 g; detekované hodnoty prošlého EMVZ - porobeton; souřadnice polohy profilu vlhkostního čela A = 5 mm; B = 15 mm

Obr. 21: Po 20 minutách navlhání detekováno vlhkostní čelo s příjmem vody +8,25 g a s vizuálním postupem vzlínající vody na úrovni cca 15 mm souřadnice polohy vlhkostního profilu A = 10 mm; B = 18 mm, viz obr. 1

Obr. 22: Po 30 minutách navlhání detekováno vlhkostní čelo s příjmem vody +9,35 g a s vizuálním postupem vzlínající vody na úrovni cca 18 mm; souřadnice polohy vlhkostního profilu A = 15 mm B = 18 mm - viz obr. 1.

Obr. 23: Po 40 minutách navlhání detekováno vlhkostní čelo s příjmem vody +12,90 g a s vizuálním postupem vzlínající vody na úrovni cca 22 mm; souřadnice polohy vlhkostního profilu A = 25 mm; B = 20 mm, viz obr. 1.

Na obr. 20 až 23 [2] je převedení časových údajů z prvotních podkladů měření podle rychlosti pojezdu vlnovodů na délkové údaje souřadnice X polohy profilu vlhkostního čela pro pórobeton v časových intervalech po 10ti minutách v programu Linregrese Excel. Z těchto grafů je zřejmé, že je možné průběžně sledovat transport vlhkosti v průběhu navlhání sledovaného vzorku kontinuálně a hodnoty udávající informace o rozložení vlhkosti zapisovat ve zvolených časových intervalech.

Tab. 5: Závislost hodnoty prošlého EMVZ v mV na vzdálenosti x v metrech od zdroje vlhkosti pro vzorek keramický pálený střep [2]. Vstupní hodnota EMVZ je vyjádřena napětím 500 mV.

vzdálenost x od zdroje vlhkosti (m) množ. záření (mV) (po 10 min. navlhání) množství záření (po 20 min. navlhání) množství záření (po 30 min. navlhání)
0432
0.00333462
0.006668126
0.01211411
0.01333-1614
0.016661372913
0.02216-14
0.0233327814830
0.02666315224156
0.03333282238
0.03333343323-
0.03666350341303
0.04357351341
0.04333361355356
0.04666363354354
0.05362352351
0.05333361354352
0.05666364359359
0.06368366368

8. Závislost EMWZ na poloze vlhkosti ve vzorku materiálu

Z hodnot naměřených pro tři různé časy (10, 20, 30 min.) od začátku navlhání, byly stanoveny metodou nejmenších čtverců v programu Maple rovnice závislosti hustoty toku prošlého záření z (udaného v mV) na vzdálenosti od zdroje vlhkosti x (m), předpoklad navlhacích křivek a „funkční závislost detekované veličiny změny intenzity EMVZ na hmotnostní vlhkosti um na délce vzorku ve vybraném časovém intervalu jeho navlhání.

9. Vyjádření křivek navlhání

Křivky navlhání jsou určeny jako grafy složených funkcí, které vzniknou složením funkcí z předchozích výpočtů. Sestrojení grafů funkcí vyjadřujících závislost obsahu vlhkosti na vzdálenosti od zdroje vlhkosti, tj. grafy složených funkcí podle předpokladů uvedených v části I [8]:

um,t = f (zt (x) ),

kde t je označení časového intervalu křivky navlhání. Vyjádření rozložení um na délce vzorku jako „navlhací křivky”

u10 = f(z10(x))
u20 = f(z20(x))
u30 = f(z30(x))

10. Součinitel vlhkostní vodivosti v závislosti na hmotnostní vlhkosti a v časových intervalech navlhání

Zjištěné informace o transportu profilu vlhkostního čela měřením na sestavené měřící aparatuře byl Matanovou metodou výpočtem podle vzorce (10) stanoven součinitel kapilární vodivosti v časových intervalech navlhání vzorku materiálu z keramického páleného střepu a pro různé hmotnostní vlhkosti jsou jeho hodnoty uvedeny v tabulce č. 6.

Tab. 6 (6 tabulek): Hodnoty κ jako charakteristické veličiny
pro daný materiál, závislé na hmotnostní vlhkosti pórovitého materiálu
získané výpočtem Matanovou metodou - na základě měření na 6ti vzorcích -
keramický pálený střep [2].

vlhkost (%) κ (10 min.) κ (20 min.) κ (30 min.)
51.627 . 10-78.879 . 10-85.798 . 10-8
101.621 . 10-79.804 . 10-88.116 . 10-8
151.845 . 10-71.129 . 10-71.011 . 10-7
203.098 . 10-71.622 . 10-71.480 . 10-7

vlhkost (%) κ (10 min.) κ (20 min.) κ (30 min.)
51.766 . 10-78.578 . 10-86.569 . 10-8
101.839 . 10-71.059 . 10-79.22 . 10-8
152.076 . 10-71.275 . 10-71.295 . 10-7
202.767 . 10-71.75 . 10-72.433 . 10-7

vlhkost (%) κ (10 min.) κ (20 min.) κ (30 min.)
51.235 . 10-75.68 . 10-87.438 . 10-8
102.409 . 10-71.316 . 10-71.162 . 10-7
153.328 . 10-72.186 . 10-71.747 . 10-7
206.207 . 10-711.85 . 10-74.251 . 10-7

vlhkost (%) κ (10 min.) κ (20 min.) κ (30 min.)
51.332 . 10-78.308 . 10-88.349 . 10-8
101.725 . 10-71.345 . 10-71.441 . 10-7
151.757 . 10-72.5 . 10-72.261 . 10-7
201.741 . 10-739.1 . 10-78.678 . 10-7

vlhkost (%) κ (10 min.) κ (20 min.) κ (30 min.)
51.239 . 10-78.2 . 10-85.598 . 10-8
101.355 . 10-78.31 . 10-85.7 . 10-8
151.253 . 10-79.439 . 10-86.536 . 10-8
201.123 . 10-71.385 . 10-71.076 . 10-7

vlhkost (%) κ (10 min.) κ (20 min.) κ (30 min.)
57.2 . 10-81.633 . 10-77.07 . 10-9
101.683 . 10-8záporné ???5.526 . 10-8
159.333 . 10-81.296 . 10-81.1 . 10-7
202.134 . 10-75.086 . 10-83.563 . 10-7

Na obr. 20 je uvedeno grafické vyjádření závislosti součinitele vlhkostní vodivosti stanoveného na základě experimentálního měření vstupních údajů na sestavené aparatuře pomocí EMVZ a vypočítaného pomocí Matanovy metody na hmotnostní vlhkosti pro keramický pálený střep.

Obr. 24: Vyjádření změny intenzity EMWZ v závislosti na hmotnostní vlhkosti v délce sledovaného vzorku materiálu v časových intervalech v průběhu jeho navlhání v programu Maple [2].

Obr. 25: Navlhací křivky vyjádřené v programu Maple pro keramický pálený střep [2].
Obr. 26: Závislost součinitele vlhkostní vodivosti na hmotnostní vlhkosti pro keramický pálený střep [2].

Závěry

F.1 Pro výpočet součinitele kapilární vodivosti dává tato metoda, ve srovnání s metodou gravimetrickou, možnost dosažení vyšší četnosti a přesnosti údajů vlhkostních poměrů v detailních řezech na délce měřeného vzorku přepočtem posuvu v měřeném čase ve směru podélné osy vzorku na úseky měřené souřadnice. Hodnoty z kontinuálního měření, jako křivky vyššího stupně, jsou z hlediska exaktnosti měření vhodným podkladem pro matematické zpracování tak, aby bylo dosahováno při modelování vlhkostního pole a výpočtech součinitele kapilární vodivosti co největšího přiblížení pravděpodobnosti skutečnému stavu. Detekce EMWZ záření v mV je detailnější právě v oblastech s nižší vlhkostí.

F.2 Tento způsob zjištění rozložení vlhkosti umožňuje ve srovnání s metodou gravimetrickou získání kontinuálně naměřených hodnot při vyloučení nepřesností vlivem lidského faktoru.

F.3 Měřící aparatura, a vytvořená metodika postupu měření umožňuje návaznost na matematické zpracování získaných výstupů k výpočtu součinitele kapilární vodivosti κ.

F.4 Výhodou je nedestruktivní, bezkontaktní, kontinuální a relativně rychlé získání výsledků měření.

F.5 Nedostatkem této metodiky měření je potřeba stanovení „cejchovní” křivky pro jednotlivé látky, to znamená, stanovení závislosti absorpce EMWZ na hmotnostní vlhkosti pro různé stavební materiály.

Příspěvek vznikl s pomocí výzkumného záměru MSM0021630511 „Progresivní stavební materiály s využitím druhotných surovin a jejich vliv na životnost konstrukcí" na fakultě stavební VUT v Brně.

Použitá literatura:

[1] Moudrý.I.: Aplikace mikrovln. záření pro stanovení změn vlhkosti dřeva a nedestruktivní zjišťování stavu dřevěných konstrukcí, VUT FAST, Brno,1985.
[2] Škramlik. J: Vlhkost v dutinách stavebních konstrukcí, VUT FAST, Disertační práce 2005.
[3] Šťastník S.: Fyzikální předpoklady pro popis komplexního vlhkostního šíření, interní zpráva, VUT FAST, Brno 2003.
[4] Garvan, F.: The Maple Book.
[5] Kutílek M.: Vlhkost pórovitých materiálů, SNTL, 1992.
[6] Mrlík F.: Vlhkostné problémy stavebných materiálov a konštrukcií, Alfa Bratislava, 1985.
[7] Kieβl, K.: Kapillarer und dampfförmiger Feuchtentransport in mehrschichtigen Bauteilen. Rechnerische Erfassung und bauphysikalische Anwendung, Dissertation Universität Essen-Fachbereich Bauwesen, Essen 1983.
[8]Moudrý, I.; Škramlik J.: Měření vlhkosti ve stavebních hmotách pomocí mikrovln (I), Stavebnictví a interiér č. 12/2006, str. 6, Měření vlhkosti ve stavebních hmotách pomocí mikrovln (I).
[9] Moudrý, I.; Škramlik J.: Měření vlhkosti ve stavebních hmotách pomocí mikrovln (II), Stavebnictví a interiér č. 1/2007, str. 38, Měření vlhkosti ve stavebních hmotách pomocí mikrovln (II).

1 Ing. Jan Škramlik, Ph.D., VUT FAST Brno, Ústav pozemního stavitelství, odborný asistent, 602 00 Brno, Veveří 95, tel.: 541 147 430, e-mail: skramlik.j@fce.vutbr.cz.
2 Doc. Ing. Ivan Moudrý, CSc., VUT FAST Brno, Ústav pozemního stavitelství, odborný asistent, 602 00 Brno, Veveří 95, tel.: 541 147 430, e-mail: moudry.i@fce.vutbr.cz.
Autor:
Foto: -