黃祖堅(jiān)， 孫一民， MUSSO Florian

(1.華南理工大學(xué) 建筑學(xué)院， 廣東 廣州 510640； 2.華南理工大學(xué) 亞熱帶建筑科學(xué)國家重點(diǎn)實(shí)驗(yàn)室， 廣東 廣州 510640； 3.慕尼黑工業(yè)大學(xué) 建筑構(gòu)造與材料教研所， 慕尼黑 80333)

對(duì)于竹原纖維(BF)的研究始于20世紀(jì)90年代，主要由某些東亞國家主持，且面向紡織業(yè)和復(fù)合材料應(yīng)用[1].竹原纖維的化學(xué)組分包含纖維素、半纖維素、木質(zhì)素果膠和脂蠟質(zhì)等[2-3]，提取方法主要有機(jī)械制備和閃爆工藝[4-5].已有研究關(guān)注的是竹原纖維的結(jié)構(gòu)特性和制備、作為服飾用品原料的熱濕舒適性以及作為復(fù)合材料增強(qiáng)元素的力學(xué)性能等.田慧敏等[6]對(duì)竹原纖維表面形貌的研究表明，竹原纖維的表面及孔隙結(jié)構(gòu)有利于其快速吸放濕分.萬玉芹等[7]研究了竹原纖維吸放濕特征曲線并推算了吸放濕速率方程，表明竹原纖維吸放濕速度較快.趙春梅[8]研究表明竹原纖維是高吸濕纖維，在環(huán)境溫度θ為20℃，環(huán)境相對(duì)濕度φ為65%時(shí)的回潮率為12%.López[9]研究表明竹材的纖維飽和點(diǎn)與其品種有關(guān)，大約在13%(Phyllostachys pubescens,毛竹)至20%(Dendrocalamus strictus,印度實(shí)竹)之間，小于木材.Liese[10]研究表明，竹材的纖維長(zhǎng)度大約在1.04mm(Phyllostachys nigra，紫竹)至2.64mm(Bambusa vulgaris cv. Wamin，大佛肚竹)之間，大于硬木而小于軟木.因?yàn)榫哂辛己玫牧W(xué)特性，竹原纖維作為復(fù)合材料的力學(xué)增強(qiáng)元素在國內(nèi)外得到了較多研究[11-12].為降低制備成本，有學(xué)者和企業(yè)研發(fā)了初級(jí)加工的竹原纖維，使其可適用于量大面廣的建筑工業(yè).姚文斌等[13-14]發(fā)明了一種麻形竹原纖維制法.魏書元[15]采用竹原纖維制成建筑保溫材料并申請(qǐng)了專利.Grosser等[16-17]研究表明竹纖維基本上在縱向平行排布，僅在竹節(jié)處有少量徑向纖維，與木材有很大區(qū)別.這一各向異性強(qiáng)烈、缺乏橫向纖維拉結(jié)的特點(diǎn)，使得初級(jí)竹原纖維的提取耗能較小[18]，有利于大規(guī)模生產(chǎn).

在建筑圍護(hù)結(jié)構(gòu)分層構(gòu)造體系中，填充材作為必要的組成元素，承擔(dān)了提高構(gòu)造整體聲學(xué)和熱工性能的作用.植物纖維是填充材體系中的重要一支，具有可再生、可降解等優(yōu)勢(shì)[19-20].Boukhattem等[21-26]分別研究了稻稈、棗棕櫚纖維、椰子纖維、香蕉纖維等作為保溫材料的熱工性質(zhì)，發(fā)現(xiàn)這些植物纖維的導(dǎo)熱系數(shù)為0.038～0.065W/(m·K)，其中木纖維的導(dǎo)熱系數(shù)為0.038～0.050W/(m·K).

與無機(jī)材料相比，植物纖維填充材會(huì)與環(huán)境，尤其是濕物理環(huán)境發(fā)生相互作用.Boukhattem等[21-22]研究表明含水率對(duì)棕櫚纖維、大麥秸稈、稻草、亞麻、麻纖維和木纖維的導(dǎo)熱系數(shù)具有顯著影響.植物纖維這一特性一方面有助于發(fā)揮其“呼吸性”[23]，對(duì)建筑構(gòu)件及其圍護(hù)的建筑空間起到熱濕調(diào)節(jié)作用；另一方面由于自身含濕量的變化，可能對(duì)其熱工性能產(chǎn)生不利影響.

植物纖維在歐美國家被廣泛應(yīng)用于低層民用建筑，相比之下中國在這一領(lǐng)域具有廣闊的發(fā)展空間.充實(shí)的竹林資源和研究基礎(chǔ)使得竹原纖維可以成為中國研發(fā)植物纖維填充材的一個(gè)優(yōu)選對(duì)象.對(duì)竹原纖維的熱濕物理性質(zhì)進(jìn)行系統(tǒng)、關(guān)聯(lián)性的研究，是將其合理應(yīng)用的基礎(chǔ).

1 材料與方法

1.1 測(cè)試對(duì)象

測(cè)試對(duì)象為初級(jí)竹原纖維(BF)，是通過碾壓出纖工藝從毛竹(Phyllostachys pubescens，剛竹屬，是中國主要的工業(yè)用竹種)莖部初步取得的韌皮纖維，可進(jìn)一步加工為可紡竹原纖維和再生竹纖維.BF的單元長(zhǎng)度為2.0～10.0cm.

1.2 測(cè)試方法

1.2.1平衡吸放濕試驗(yàn)

平衡吸放濕試驗(yàn)用于測(cè)試BF的等溫吸放濕曲線，以獲得表征其濕儲(chǔ)存性質(zhì)的平衡含水率u值.試驗(yàn)時(shí)，稱取預(yù)定質(zhì)量的干燥BF填充于玻璃皿中，獲得填充密度ρ為110kg/m3的BF試件，1式6份，吸濕、放濕過程各用3份.主要試驗(yàn)裝置有：真空干燥器9個(gè)，內(nèi)徑為30cm；溫濕度記錄儀TH10R-EX 1個(gè)及外置傳感器9個(gè)，溫濕度精度分別為±0.2℃和±2%；玻璃皿51個(gè)，內(nèi)徑為100mm，高度為20mm.

1.2.2干燥試驗(yàn)

干燥試驗(yàn)用于測(cè)試BF的干燥曲線，據(jù)此計(jì)算材料的干燥速率，并以干燥速率作為表征材料濕傳遞性質(zhì)的吸放濕速率U值.試驗(yàn)時(shí)稱取預(yù)定質(zhì)量的干燥BF，填充于亞克力透明管容器中，獲得填充密度ρ為110kg/m3的BF試件，1式3份.主要試驗(yàn)裝置有：SHIMADZU UX6200H電子天平,精度0.01g；HWS-250B恒溫恒濕箱，溫濕度精度分別為±0.5℃ 和±5%；亞克力透明管容器12個(gè)，底面密封，預(yù)留5個(gè)直徑5mm的氣孔，容器內(nèi)徑100mm，高度h分別為40，60，80，100mm.

1.2.3導(dǎo)熱系數(shù)測(cè)試

導(dǎo)熱系數(shù)λ值用于表征材料的熱傳遞性質(zhì).試驗(yàn)時(shí)稱取預(yù)定質(zhì)量的干燥BF，參考木纖維保溫板的常用密度范圍，將填充密度70～170kg/m3作為BF的測(cè)試范圍，故將干燥BF填充于30cm×30cm的EPS模板中(模板厚度為28mm，在中部預(yù)留邊長(zhǎng)為20cm的正方形孔洞用于填充BF)，獲得填充密度ρ分別為70，90，110，130，150，170kg/m3的BF試件，1式2份.主要試驗(yàn)裝置有：CD-DR3030平板導(dǎo)熱儀，熱板θ=35℃，冷板θ=15℃，中心傳熱區(qū)域尺寸為15cm×15cm，精度±2%；HWS-250B恒溫恒濕箱，溫濕度精度分別為±0.5℃，±5%(養(yǎng)護(hù)過程中采用精度分別為±0.2℃，±2%的TH10R-EX溫濕度記錄儀對(duì)箱內(nèi)溫濕度環(huán)境進(jìn)行實(shí)時(shí)監(jiān)測(cè)).

采用國際標(biāo)準(zhǔn)ISO 8302—1991推薦的方法[26]，首先采用CD -DR3030平板導(dǎo)熱儀對(duì)干燥試件進(jìn)行導(dǎo)熱系數(shù)測(cè)試，得出干燥材料導(dǎo)熱系數(shù)λd.再將試件置入恒溫恒濕箱中養(yǎng)護(hù)至設(shè)定的平衡含水狀態(tài)，養(yǎng)護(hù)環(huán)境的相對(duì)濕度φ設(shè)置為40%，50%，60%，70%，80%，90%.2周后將試件取出稱重，并測(cè)試對(duì)應(yīng)含濕量的導(dǎo)熱系數(shù)，完成測(cè)試后將試件放回恒溫恒濕箱中，設(shè)置下一梯度濕度環(huán)境進(jìn)行養(yǎng)護(hù).重復(fù)以上操作，直至獲得所有填充密度梯度和濕度梯度試件的λ值.

2 結(jié)果分析

20世紀(jì)中期，Luikov[27-28]提出用于描述多孔材料中熱濕耦合遷移的高度耦合的非線性偏微分方程.1995年，Künzel[29]改進(jìn)了Luikov模型，提出以相對(duì)濕度φ為濕傳遞驅(qū)動(dòng)勢(shì)，并研發(fā)出計(jì)算建筑圍護(hù)結(jié)構(gòu)內(nèi)熱濕耦合傳遞的數(shù)值模擬計(jì)算軟件WUFI.材料的熱濕物理性質(zhì)通常用含濕量、密度、環(huán)境溫度和相對(duì)濕度等因子進(jìn)行關(guān)聯(lián)性描述，有時(shí)1項(xiàng)性質(zhì)對(duì)應(yīng)有數(shù)十種擬合方程.本文的測(cè)試結(jié)果顯示：BF的u，U和λ值與環(huán)境相對(duì)濕度φ值之間存在指數(shù)函數(shù)關(guān)系，可以采用Y=a·eb·φ+c的一般方程進(jìn)行描述.

2.1 u -φ測(cè)試結(jié)果分析

BF的平衡吸放濕試驗(yàn)結(jié)果見表1.BF的等溫吸放濕u-φ曲線近似屬于第2類型，即反“S”型曲線.通常曲線前段采用BET方程，后段采用Kelvin方程加以描述.將表1中的平衡含水率u對(duì)相對(duì)濕度φ進(jìn)行指數(shù)函數(shù)擬合，有u=0.0088e3.363 7φ(擬合度R2=0.9495,顯著性Sig=0.000009).其中吸濕等溫曲線為u=0.0039e4.282 7φ(R2=0.9584,Sig=0.000004)；放濕等溫曲線為u=0.0145e2.836 3φ(R2=0.9274,Sig=0.000031).從所得的R2和Sig值可見，采用指數(shù)函數(shù)可以較好地描述BF的u-φ曲線.

BF吸濕量高，在φ=96.3%時(shí)其u值可以達(dá)到25.84%.在本測(cè)試范圍內(nèi)，φ=85.4%～96.3%對(duì)應(yīng)

表1 真空干燥器內(nèi)溫濕度控制及平衡含水率測(cè)試結(jié)果

的曲線段斜率大幅高于φ=11.2%～85.4%對(duì)應(yīng)的曲線段斜率.這一方面是因?yàn)榍罢吒叩乃魵夥謮旱耐庠谧饔?，另一方面是由于BF在φ=85.4%～96.3%曲線段中達(dá)到了纖維飽和點(diǎn)，進(jìn)入自由水吸附階段，從而使得其平衡含水率大幅增加.

材料的等溫吸放濕曲線是一定的環(huán)境相對(duì)濕度及其所對(duì)應(yīng)的平衡含水率的集合，由于材料達(dá)到其平衡含水率的過程不同，所得結(jié)果會(huì)有差異，通常放濕過程所得u值會(huì)不同程度地高于吸濕過程所得u值，這種情況稱為滯后效應(yīng).滯后效應(yīng)因材而異，Künzel[29]對(duì)4種常用礦物建筑材料進(jìn)行測(cè)試的結(jié)果顯示其滯后效應(yīng)不明顯，認(rèn)為在此類材料中采用等溫吸濕曲線對(duì)于實(shí)際工程而言已經(jīng)足夠準(zhǔn)確.對(duì)于纖維類材料而言，干燥過程中纖維素分子的羥基會(huì)部分借助副價(jià)健直接相連而達(dá)到相互飽和，因而削弱了之后材料對(duì)水分的吸附能力，使其吸濕過程的u值偏小.本文中，在φ=11.2%～85.4%的范圍內(nèi)，等溫吸濕和放濕過程有較為明顯的滯后效應(yīng)，放濕過程所得u值比吸濕過程高0.52%～3.50%.因此建議在BF平衡含水率取值時(shí)，采用吸濕過程和放濕過程所得u值的平均值.

2.2 U-φ測(cè)試結(jié)果分析

BF的吸放濕速率U值的測(cè)試結(jié)果見表2.將表2中的U值與材料平衡含水率對(duì)應(yīng)的相對(duì)濕度φ值進(jìn)行指數(shù)函數(shù)擬合，有U=2.7610×10-10e0.111 0φ(R2=0.9960,Sig=9.0054×10-10).考慮試件高度h的影響，有U=(1.009×10-10-1.396×10-10lnh)·e0.111 0φ(R2=0.9500,Sig=0.026).從所得的R2和Sig值可見，采用指數(shù)函數(shù)可以較好地描述BF的U-φ曲線.

BF的干燥曲線表現(xiàn)出前后兩個(gè)階段的特點(diǎn).在測(cè)試范圍內(nèi)，φ=75%～95%階段的U值大幅高于φ= 50%～75%階段.前一階段U值主要取決于環(huán)境溫度和相對(duì)濕度，后一階段U值主要取決于自身物性.4組不同高度試件的測(cè)試結(jié)果對(duì)比表明，在h= 40～100mm 時(shí)，U值受試件高度影響不明顯.

表2 竹原纖維吸放濕速率U值

2.3 λ-φ測(cè)試結(jié)果分析

圖1為BF導(dǎo)熱系數(shù)測(cè)試結(jié)果.對(duì)圖1中的測(cè)試結(jié)果進(jìn)行Spearman相關(guān)性分析，所得λ-φ和λ-u的相關(guān)系數(shù)R分別為0.8688和0.8406，其雙尾顯著性分別為8.833×10-14和3.2946×10-12.2組數(shù)據(jù)相關(guān)性均在0.01層上顯著(雙尾)，其中λ-φ比λ-u相關(guān)性略強(qiáng).在另一方面，BF的λ值與u值之間存在“S”型曲線關(guān)系，在試件從干燥狀態(tài)直至平衡含水率較低的初始階段，λ值增大速率隨著平衡含水率的增加而增大；在平衡含水率達(dá)到一定值后，在本測(cè)試中為u=11.5%～13.0%時(shí)，λ值增大速率減緩.“S”型曲線較難采用單一函數(shù)進(jìn)行描述，在WUFI等基于HAM建筑圍護(hù)結(jié)構(gòu)熱濕耦合傳遞的計(jì)算模型中，通常采用線性方程結(jié)合含水量補(bǔ)償系數(shù)b值進(jìn)行近似估算[29].在本文中，采用線性函數(shù)對(duì)λ-u測(cè)試結(jié)果進(jìn)行擬合，有λ=0.0421-0.016870u(R2=0.8976,Sig=0.004073)，而采用指數(shù)函數(shù)方程對(duì)λ-φ測(cè)試結(jié)果進(jìn)行擬合，有λ=0.00457e4.397 4φ+0.043 7(R2=0.9955,Sig=0.000008)；從所得的R2和Sig值可見，至少在本文中以φ作為自變量的指數(shù)函數(shù)，可提高BF的導(dǎo)熱系數(shù)與其濕物理性質(zhì)關(guān)系的擬合度.

圖1 竹原纖維導(dǎo)熱系數(shù)Fig.1 Thermal conductivity of BF

不同填充密度BF干燥試件的導(dǎo)熱系數(shù)為0.0423～0.0465W/(m·K)，與木纖維保溫板性質(zhì)相近.在填充密度ρ=70～170kg/m3范圍內(nèi)，λ值隨ρ值增大而增大，且不同填充密度的BF試件之間的導(dǎo)熱系數(shù)差異隨著環(huán)境相對(duì)濕度的增大而增大，如在干燥條件下，ρ=170kg/m3試件和ρ=70kg/m3試件的λ值相差9.46%，而在φ=90%條件下，這2種試件的λ值相差達(dá)到20.88%，表明BF的導(dǎo)熱系數(shù)受環(huán)境相對(duì)濕度的影響與自身填充密度存在正相關(guān)關(guān)系.

2.4 相關(guān)熱工性能參數(shù)推算

選擇填充密度為110kg/m3，高度為60mm的試件作為個(gè)案，設(shè)置3組對(duì)比計(jì)算對(duì)象，其中組1采用本測(cè)試所得BF的u-φ，U-φ和λ-φ函數(shù)，以1h為時(shí)間步長(zhǎng)計(jì)算24h內(nèi)的λ0，S0，D0值，作為計(jì)算參照；組2將組1中的U值減小一半，u值和λ值不變，以1h為時(shí)間步長(zhǎng)計(jì)算24h內(nèi)的λU，SU和DU值，計(jì)算結(jié)果用于評(píng)價(jià)吸放濕速率U值對(duì)試件熱工性能的影響；組3將組1中u值提高2倍，U值和λ值不變，以1h為時(shí)間步長(zhǎng)計(jì)算24h內(nèi)的λemc，Semc和Demc值，計(jì)算結(jié)果用于評(píng)價(jià)平衡含水率u值對(duì)試件熱工性能的影響.所得結(jié)果見圖2.

計(jì)算結(jié)果顯示，在相對(duì)濕度周期性變化的環(huán)境下，BF填充材的熱工性能會(huì)隨之波動(dòng)，采用極值差Dev作為表征其波動(dòng)幅度的指標(biāo)，則本測(cè)試所得BF的導(dǎo)熱系數(shù)λ0，蓄熱系數(shù)S0和熱惰性指標(biāo)D0的Dev值分別為16.60%，12.12%和4.28%.

圖2 24h周期的相對(duì)濕度φ，導(dǎo)熱系數(shù)λ，蓄熱系數(shù)S和熱惰性指標(biāo)D的計(jì)算結(jié)果Fig.2 φ in cycle of 24h and the calculation results of λ, S and D

當(dāng)BF的吸放濕速率U值降為原來的1/2時(shí)，BF填充材的熱工性能波動(dòng)幅度大幅縮小，其中λU，SU和DU的Dev值分別降低至0.80%，0.59%和0.21%，且表征熱傳遞性質(zhì)的導(dǎo)熱系數(shù)λ值得到明顯改善，24h平均值由0.0536W/(m·K)降至0.0518W/(m·K).當(dāng)BF的平衡含水率u值增至原來的2倍時(shí)，BF填充材的熱工性能波動(dòng)幅度也明顯下降，其中λemc，Semc和Demc的Dev值分別降至8.18%，7.72%和1.12%，且同時(shí)表征熱儲(chǔ)存和熱傳遞的蓄熱系數(shù)S值和熱惰性指標(biāo)D值均得到了明顯改善：S的24h平均值從1.8559W/(m2·K)提高到1.9255W/(m2·K)，D的平均值從2.1215提高到了2.2404.

3 結(jié)論

(1)采用指數(shù)函數(shù)Y=a·eb·φ+c，可以將u，U和λ值分別與環(huán)境相對(duì)濕度φ進(jìn)行較好的擬合.

(2)在測(cè)試范圍內(nèi)，受環(huán)境相對(duì)濕度φ及自身特性影響，BF的等溫吸放濕u-φ曲線在φ為85.4%～96.3%、吸放濕速率U-φ曲線在φ為75%～95%時(shí)快速增大.

(3)在填充密度ρ為70～170kg/m3范圍內(nèi)，干燥BF的λ值為0.0423～0.0465W/(m·K)，與木纖維保溫板性質(zhì)相近；潮濕BF的λ值隨φ值增大而增大的速率與其自身填充密度呈現(xiàn)正相關(guān)關(guān)系.

(4)在相對(duì)濕度動(dòng)態(tài)變化的環(huán)境中，λ值，S值和D值會(huì)相應(yīng)波動(dòng)，波幅與U值呈正相關(guān)關(guān)系，且S值和D值隨u值的增大而增大，λ值隨U值的減小及u值的增大而減小，表明增大u值和減小U值均有利于改善BF的熱工性能.

參考文獻(xiàn)：

[1] LI H,ZHOU G Y,ZHANG H Y.Research and utilization status of natural bamboo fiber[J].Advanced Materials Research,2010,159:236-241.

[2] 蔡再生.纖維化學(xué)與物理[M].北京:中國紡織出版社,2004:19-20.

CAI Zaisheng.Fiber chemistry and physics [M].Beijing:China Textile & Apparel Press,2004:19-20.(in Chinese)

[3] WANG Yueping,WANG Ge,CHENG Haitao,et al.Structures of natural bamboo fiber for textiles[J].Textile Research Journal,2010,80(1):334-343.

[4] ZAKIKHANI P,ZAHARI R,SULTAN M T H,et al.Extraction and preparation of bamboo fibre-reinforced composites[J].Materials and Design,2014,63:820-828.

[5] KHALIL H P S A,BHAT I U H,JAWAID M,et al.Bamboo fibre reinforced biocomposites:A review[J].Materials & Design,2012,42:353-368.

[6] 田慧敏,蔡玉蘭.竹原纖維微觀形態(tài)及聚集態(tài)結(jié)構(gòu)的研究[J].棉紡織技術(shù),2008,36(9):544-546.

TIAN Huimin,CAI Yulan.Research of nature bamboo fiber microstructure and aggregation state structure [J].Cotton Textile Technology,2008,36(9):544-546.(in Chinese)

[7] 萬玉芹,吳麗莉,俞建勇.竹纖維吸濕性能研究[J].紡織學(xué)報(bào),2004,25(3):14-16.

WAN Yuqin,WU Lili,YU Jianyong.Study on hygroscopic properties of bamboo fiber[J].Journal of Textile Research,2004,25(3):14-16.(in Chinese)

[8] 趙春梅.竹原纖維性能研究及其產(chǎn)品開發(fā)[D].上海：東華大學(xué),2010.

ZHAO Chunmei.Research on performance of natural bamboo fiber and product development[D].Shanghai:Donghua University,2010.(in Chinese)

[10] LIESE W.INBAR technical report 18—The anatomy of bamboo culms[C]//International Network for Bamboo and Rattan.Beijing:INBAR,1998:58-62.

[11] FEI P,XIONG H,CAI J,et al.Enhanced the weatherability of bamboo fiber-based outdoor building decoration materials by rutile nano-TiO2[J].Construction and Building Materials,2016,114:307-316.

[12] OUYANG Weijia,HU Weiping,YANG Xiaoping.On the application of bamboo fiber in home furnishing design[J].Advanced Materials Research,2012,460:244-248.

[13] 姚文斌,愈偉鵬,張蔚.一種麻形竹纖維的生產(chǎn)工藝:CN 102206873 B[P].2014-05-14.

YAO Wenbin,YU Weipeng,ZHANG Wei.A hemp production technology shaped bamboo fiber:CN 102206873 B[P].2014-05-14.(in Chinese)

[14] YAO Wenbin,ZHANG Wei.Research on manufacturing technology and application of natural bamboo fibre[C]//International Conference on Intelligent Computation Technology and Automation.[S.l.]:IEEE,2011:143-148.

[15] 魏書元.一種竹纖維建筑保溫材料:CN 103553423 A [P].2014-02-05.

WEI Shuyuan.One kind of bamboo fiber building insulation materials:CN 103553423 A[P].2014-02-05.(in Chinese)

[16] GROSSER D.Beitrag zur Histologie und Klassifikation asiatischer Bambusarten[D].Hamburg:Kommissionsverlag Buchhandlung Max Wiedebusch Hamburg,1971.

[17] ZOU Linhua,JIN Helena,LU Weiyang.Nanoscale structural and mechanical characterization of the cell wall of bamboo fibers[J].Materials Science and Engineering C,2009,29:1375-1379.

[18] DEFOIRDT N,BISWAS S,VRIESE L D,et al.Assessment of the tensile properties of coir,bamboo and jute fibre[J].Composites Part A:Applied Science and Manufacturing,2010,41(5):588-595.

[19] MURPHY R J,NORTON A.Life cycle assessments of natural fibre insulation materials final report[R].London:Imperial College,2008:75-76.

[20] HURTADO P L,ROUILLY A,VANDENBOSSCHE V,et al.A review on the properties of cellulose fibre insulation[J].Building and Environment,2016,96(1):170-177.

[21] BOUKHATTEM L,BOUMHAOUT M,HAMDI H,et al.Moisture content influence on the thermal conductivity of insulating building materials made from date palm fibers mesh [J].Construction and Building Materials,2017,148:811-823.

[22] PALUMBO M,LACASTA A M ,HOLCROFT N,et al.Determination of hygrothermal parameters of experimental and commercial bio-based insulation materials[J].Construction and Building Materials,2016,124:269-275.

[23] HOLCROFT N,SHEA A.Moisture buffering and latent heat effects in natural fibre insulation materials[C]//Portugal SB13-Contribution of Sustainable Building to Meet the EU.[S.l.]:[s.n.],2013:221-228.

[24] ISO 12571:2012 Hygrothermal performance of building materials and products-Determination of hygroscopic sorption properties[S].

[25] ASTM-04a(2016) Standard test method for hygroscopic sorption isotherms of building materials[S].

[26] ISO 8302—1991 Thermal insulation; determination of steady-state thermal resistance and related properties; guarded hot plate apparatus[S].

[27] LUIKOV A V.Heat and mass transfer in capillary-porous bodies[M].Amsterdam:Elsevier,1964:123-184.

[28] LUIKOV A V.Systems of differential equations of heat and mass transfer in capillary-porous bodies(review)[J].International Journal of Heat and Mass Transfer,1975,18(1):1-14.

[29] KüNZEL H M.Simultaneous heat and moisture transport in building components[R].Suttgart:Fraunhofer IRB Verlag,1995:7-8.