竹炭濕物性測試及其對外墻熱濕性能影響的模擬

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

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

竹炭(BC)的研究與應(yīng)用興起于20世紀(jì)90年代，主要發(fā)生在東亞及非洲國家.竹炭由纖維素、半纖維素和少量木質(zhì)素經(jīng)不同程度地分解炭化而成[1]，應(yīng)用于空氣凈化[2]、土壤改良、水凈化[3]，以及生活用品調(diào)濕性填充物[4]和復(fù)合材料中[5-6].未經(jīng)活化處理的竹炭孔隙結(jié)構(gòu)中以大孔為主，微孔形成和發(fā)展在500℃以后開始；竹炭內(nèi)部孔徑分布隨溫度升高從0.25～50.00μm收縮到0.55～5.50μm；比表面積在600～700℃之間急劇發(fā)展，到700℃時可達(dá)385m2/g，是普通木炭的2～5倍，此時竹炭具有最好的調(diào)濕機(jī)能[7-8].蔣正武等[9]以平衡吸放濕量、吸放濕效率和調(diào)濕穩(wěn)定性作為指標(biāo)評價竹炭調(diào)濕能力，并研究工藝參數(shù)對竹炭調(diào)濕性能的影響.盧克陽等[10]對竹炭吸濕性能的研究表明，顆粒大小對竹炭吸濕量影響很小，比表面積對吸濕量有不成比例關(guān)系的影響，環(huán)境溫度升高會導(dǎo)致竹炭的吸濕性能降低.

日本對調(diào)濕材料基礎(chǔ)理論有較為充實(shí)的研究積累，1949年，西藤等[11]提出“調(diào)濕材料”的概念，指不需要借助人工能源和機(jī)械設(shè)備，依靠材料自身的吸放濕性能，感應(yīng)所調(diào)空間空氣溫濕度的變化，自動調(diào)節(jié)空氣相對濕度的材料；Maki 等[12]提出用材料水蒸氣擴(kuò)散系數(shù)和最大平衡含濕量來衡量材料的調(diào)濕性能，兩者都大的材料調(diào)濕性能好.Ikeda[13]基于材料平衡吸放濕曲線u(x,T)的分析，提出以k=(?u/?x)T和v=(?u/?T)x分別表示空氣溫度T一定和空氣絕對濕度x一定時，材料平衡含水率u隨空氣相對濕度變化的能力.Qin等[14]對香港和巴黎的案例進(jìn)行模擬計(jì)算，表明調(diào)濕材料對室內(nèi)空氣濕度和建筑能耗具有改善作用.在建筑調(diào)濕材料體系中，多孔類材料包括竹炭、沸石、高嶺土等，它們均具有透過性良好、吸放濕速度快、放濕滯后小的特點(diǎn)，且生產(chǎn)工藝簡單、成本低，因此在民間應(yīng)用最為廣泛[15].在中國和日本，調(diào)濕性竹炭被用于調(diào)節(jié)室內(nèi)空氣濕度和建筑構(gòu)造含濕量，并起抑菌殺菌作用[16].

依托豐富的竹林資源，中國成為主要竹炭生產(chǎn)國和出口國之一.將竹炭作為調(diào)濕性填充材應(yīng)用于建筑工業(yè)，具有利用竹資源和提高建筑性能的雙重意義.本文基于建筑圍護(hù)結(jié)構(gòu)熱濕過程模型，系統(tǒng)地對竹炭進(jìn)行濕物性測試，研究竹炭對建筑構(gòu)件和圍護(hù)空間熱濕性能的影響，為其合理應(yīng)用提供基礎(chǔ).

1 試驗(yàn)

1.1 原材料

BC樣品為600℃炭化工藝所得的調(diào)濕性竹炭，填充密度為270kg/m3，單元直徑1.0～1.5cm.其表面形貌如圖1所示(采用蔡司EVO18掃描式電子顯微鏡獲得圖像).

圖1 BC樣品表面形貌Fig.1 BC sample surface morphology

1.2 平衡含水率測試

平衡吸放濕試驗(yàn)采用國際標(biāo)準(zhǔn)ISO 12571:2012和美國標(biāo)準(zhǔn)ASTM-04a(2016)推薦的方法[17-18]，用于測試BC的平衡含水率u值，如圖2所示.將烘干并稱量好的BC置于玻璃皿中(其質(zhì)量為md)，1式6份，吸濕、放濕過程各3份.將3份烘干試件靜置養(yǎng)護(hù)至20℃，然后置入真空干燥器中，真空干燥器底部盛有不同飽和鹽溶液以控制預(yù)定的相對濕度.將真空干燥器連同溫濕度記錄儀整體移入恒溫恒濕箱中，在20℃，相對濕度φ=50%條件下進(jìn)行養(yǎng)護(hù)，每天監(jiān)測真空干燥器內(nèi)溫濕度的變化，直至穩(wěn)定到預(yù)定值.4周后，每隔3d 打開真空干燥器對試件進(jìn)行稱重，完成稱重后立刻放回并關(guān)閉真空干燥器，重復(fù)以上操作，直至連續(xù)3次稱重所得結(jié)果差異不超過0.1%時，認(rèn)為材料已經(jīng)達(dá)到了吸濕平衡.記錄潮濕試件質(zhì)量mw.根據(jù)式(1)計(jì)算吸濕試件的平衡含水率u，得出其等溫吸濕曲線圖.

(1)

將3份在恒溫恒濕箱(20℃，φ=90%)環(huán)境中養(yǎng)護(hù)2周直至吸濕平衡的試件置入底部盛有除K2SO4飽和鹽溶液外的其余8個真空干燥器中，重復(fù)以上操作，得出放濕試件的等溫放濕曲線圖.

圖2 平衡吸放濕試驗(yàn)Fig.2 Absorption and desorption test

1.3 蒸汽滲透試驗(yàn)

蒸汽滲透試驗(yàn)用于測試BC材料的蒸汽滲透系數(shù)δ值，采用國際標(biāo)準(zhǔn)ISO 12572:2001(E)和美國標(biāo)準(zhǔn)ASTM E 96/E 96M-2005等推薦的方法[19-21]，如圖3所示，采用亞克力管和玻璃碗制作符合標(biāo)準(zhǔn)的干、濕杯試驗(yàn)裝置，其中上方的亞克力管用于盛放BC試件，下方的玻璃碗用于盛放干燥劑/飽和鹽溶液，分別形成干、濕杯.將烘干并稱量好的BC填充于亞克力管容器中(亞克力管圓形開口底部預(yù)先粘好不銹鋼絲網(wǎng)，保證竹炭不掉落，并且不影響竹炭吸濕)，亞克力管固定于塑料碗蓋上表面，碗蓋預(yù)留直徑100mm的圓形開口，記錄試件及塑料碗蓋整體干燥質(zhì)量md,lid，再整體移入恒溫恒濕箱中，在23℃，φ=50%條件下進(jìn)行養(yǎng)護(hù).2周后，試件達(dá)到吸濕平衡，將試件及塑料碗蓋密封于玻璃碗上，玻璃碗內(nèi)盛有粒徑Φ=3mm的CaCl2干燥劑(干杯)或K2SO4飽和鹽溶液(濕杯)以控制試件下方相對濕度，得到干、濕杯試件各3份.將干、濕杯移入恒溫恒濕箱中，在23℃，相對濕度50%條件下進(jìn)行養(yǎng)護(hù).每隔3d對干、濕杯進(jìn)行稱重，當(dāng)總質(zhì)量連續(xù)5次以上呈線性變化后停止養(yǎng)護(hù)，立刻取下碗蓋進(jìn)行稱重得潮濕質(zhì)量mw,lid，對記錄的質(zhì)量和時間數(shù)據(jù)進(jìn)行線性擬合，所得質(zhì)量變化速率即為蒸汽傳遞速率G(kg/s).根據(jù)式(2)計(jì)算試件蒸汽滲透系數(shù)δ[kg/(m·s·Pa)]：

(2)

式中：d為試件厚度，m；A為碗蓋開口面積，m2；Tk為絕對溫度，K；Δφ為試件兩側(cè)空氣相對濕度差，%；R0為理想氣體常數(shù)，R0=461.5N·m/(kg·K)；p為氣壓，Pa；p0為標(biāo)準(zhǔn)大氣壓，p0=101325Pa；hair為試件下表面空氣厚度，m.

圖3 蒸汽滲透試驗(yàn)Fig.3 Water vapor transmission test

1.4 外墻熱濕性能模擬

建筑圍護(hù)結(jié)構(gòu)熱濕空氣流動耦合模型(HAM模型)用于描述圍護(hù)結(jié)構(gòu)熱量、濕分及空氣的儲存和傳遞過程，通常采用高度耦合的非線性偏微分方程來同時描述圍護(hù)結(jié)構(gòu)中的熱濕傳遞過程，再通過數(shù)值方法求解.這一計(jì)算方法由于對材料參數(shù)、計(jì)算時間等要求更高，因此在實(shí)際應(yīng)用中受到局限.Künzel[22]通過理論推算和實(shí)測檢驗(yàn)將熱濕傳遞計(jì)算所需材料參數(shù)進(jìn)行了簡化，選擇蒸汽壓和相對濕度作為氣態(tài)水和液態(tài)水傳遞的驅(qū)動勢，并形成計(jì)算機(jī)程序WUFI(w?rme-und feuchtetransport instation?r，非穩(wěn)態(tài)熱濕過程模擬軟件，由德國Fraunhofer建筑物理研究所開發(fā))的基礎(chǔ)模型.

將本文測試所得BC材料參數(shù)輸入WUFI Plus軟件中，設(shè)置建筑外墻模型組進(jìn)行全年熱濕性能模擬，如圖4所示.圖中BC代表含BC構(gòu)造層的模型組，NC代表不含BC構(gòu)造層的模型組；BFB為竹重組材，BSB為竹集成材，BPB為竹刨花板，這3種竹板的性質(zhì)參數(shù)是筆者通過相同的試驗(yàn)條件和方法測試所得；A concrete代表加氣混凝土，L S brick為石灰磚，這2種材料的性質(zhì)參數(shù)是采用WUFI Plus軟件中的默認(rèn)值.模型外部條件設(shè)置時考慮到中國濕氣候參數(shù)中缺乏降雨數(shù)據(jù)，無法形成墻體液態(tài)水的驅(qū)動勢，會給模擬結(jié)果造成誤差[23]，而根據(jù)K?ppen氣候分類法，類型豐富的北美氣候區(qū)對緯度相近的中國典型氣候區(qū)具有較為完整的參考價值，因此使用北美的3個熱濕及溫濕性氣候區(qū)中6個代表城市的氣候數(shù)據(jù).內(nèi)部條件設(shè)置中建立尺寸為開間×進(jìn)深×層高=3.0m×3.0m×3.0m的空間單元.設(shè)置8:00～17:00為房間工作時間，并賦予標(biāo)準(zhǔn)辦公間使用強(qiáng)度的室內(nèi)熱濕負(fù)荷；邊界條件選擇四面外墻作為研究對象，將頂面天花及底面地板設(shè)為絕緣模式.構(gòu)造上設(shè)置常用的3種類型，包括L型(外面板+填充層+內(nèi)面板)、M型(砌體層+填充層+內(nèi)面板)和H型(砌體層+內(nèi)面板).各組構(gòu)造傳熱系數(shù)U值控制在1.0W/(m2·K)左右，采用40mm厚的BC材料作為填充層靠外側(cè)布置.模型運(yùn)行過程中首先關(guān)閉暖通空調(diào)系統(tǒng)(HVAC)，模擬無設(shè)備條件下室內(nèi)熱濕環(huán)境和外墻熱濕性能；然后打開HVAC，以理想的HVAC維持室內(nèi)溫濕度，計(jì)算全年HVAC需求；運(yùn)行時間為2a，對第2年全年數(shù)據(jù)進(jìn)行采集.

圖4 WUFI Plus建筑外墻熱濕性能模擬Fig.4 WUFI Plus hygrothermal performance simulation for exterior walls

2 結(jié)果與分析

2.1 材料測試結(jié)果

BC試件平衡含水率及蒸汽滲透系數(shù)測試結(jié)果分別見表1,2.由表1,2可知：BC吸濕量高，在φ=96.3%時達(dá)到平衡含水率u=16.67%，體積含濕量w=45.01kg/m3；在測試的相對濕度內(nèi)(11.2%～85.4%)，放濕過程所得u值比吸濕過程所得u值高0.25%～1.75%，放濕滯后小，有利于發(fā)揮BC的調(diào)濕作用；BC蒸汽滲透系數(shù)干、濕杯測試結(jié)果比較接近，取均值為9.32×10-11kg/(m·s·Pa).

從等溫吸放濕曲線(圖5)可見，BC的等溫吸放濕曲線在φ=33.4%～85.4%范圍內(nèi)向上凸起，表明炭化處理使BC在這一范圍內(nèi)的調(diào)濕容量得到擴(kuò)大.此外，由各材料蒸汽滲透系數(shù)對比(圖6)可知，BC蒸汽滲透系數(shù)為竹板蒸汽滲透系數(shù)的10.0～213.2倍，濕傳遞速率大幅提高，表明BC的調(diào)濕速率得到提高.與構(gòu)造中其余材料相比，BC在調(diào)濕容量和調(diào)濕速率方面的優(yōu)勢有利于發(fā)揮其調(diào)濕作用.

2.2 WUFI Plus模擬結(jié)果

以開啟理想HVAC工況下的年HVAC開啟時間HVACperiod，年供暖和制冷量P值以及加濕和除濕量H值來表征圍護(hù)空間的HVAC需求；以關(guān)閉HVAC工況下的室內(nèi)空氣溫度θi的均值θi-mean和振幅θi-amp、外墻內(nèi)表面溫度θis的均值θis-mean和振幅θis-amp、相對濕度φi的均值φi-mean和振幅φi-amp來表征圍護(hù)空間的熱濕環(huán)境質(zhì)量；通過建筑四面外墻的熱流量Hflow和濕流量Mflow，以及L型和M型構(gòu)造的保溫填充層(cellulose fiber)含濕量wCF來表征外墻構(gòu)造熱濕性能.對含BC構(gòu)造層的模型組(BC模型組)和不含BC構(gòu)造層的模型組(NC模型組)的模擬結(jié)果進(jìn)行對比分析.

表1 真空干燥器內(nèi)溫濕度控制及平衡含水率測試結(jié)果Table 1 Temperature(θ) and relative humidity of ambient air control in vacuum dryers and the measured equilibrium moisture content

表2 干濕杯內(nèi)濕度控制及蒸汽滲透系數(shù)測試結(jié)果Table 2 Relative humidity control in cups and the measured water vapor transfer coefficient

圖5 BC、竹板和砌體材料平衡含水率對比Fig.5 Equilibrium moisture content comparison among BC, bamboo panels and masonry

圖6 BC、竹板和砌體材料蒸汽滲透系數(shù)對比Fig.6 Water vapor transfer coefficient comparison among BC, bamboo panels and masonry

2.2.1BC對建筑構(gòu)件熱濕性能的影響

以南側(cè)外墻為例分析L型和M型構(gòu)造填充層的全年含濕量，結(jié)果見圖7.由圖7(a)可見，BC對wCF的均值并不產(chǎn)生明顯影響，但能明顯減小其變化振幅，對于HVAC-on工況，wCF振幅減小了0.71～104.52kg/m3，以BC/NC模型組的模擬結(jié)果比值作為指標(biāo)，評價BC對構(gòu)造熱濕性能的影響.對于wCF振幅，BC模型組為NC模型組的6.44%～85.58%，可見BC對填充層含濕量起到顯著穩(wěn)定作用.由圖7(b),(c)可見，BC由于自身吸放濕特性，對外墻濕流起到緩沖作用，從而可減小外墻濕熱流量.對于濕流量Mflow，BC模型組為NC模型組的84.47%～101.55%；相應(yīng)地，對于熱流量Hflow，BC模型組為NC模型組的67.79%～92.45%.

2.2.2BC對圍護(hù)空間熱濕性能的影響

圍護(hù)空間室內(nèi)熱濕環(huán)境和年HVAC需求模擬結(jié)果見圖8.由圖8可知，BC的添加并不對θi-mean，φi-mean和θis-mean產(chǎn)生明顯影響，但可以減小相應(yīng)的振幅θi-amp，φi-amp和θis-amp.由圖8(a)可見，在HVAC-off工況下，BC模型組φi-amp最大減幅達(dá)4.0%；由圖8(b) 可見，BC模型組θi-amp相比于NC模型組減小了0.7～2.5℃；由圖8(c)可見，以南側(cè)外墻為例，BC模型組θis-amp比NC模型組θis-amp減小了0.9～3.0℃，表明BC更有助于維持室內(nèi)熱濕環(huán)境的穩(wěn)定；由圖8(d)可見，BC的添加有助于縮短HVAC開啟時間HVACperiod，BC模型組HVACperiod為NC模型組的91.12%～96.39%；由圖8(e),(f)可見，在打開理想HVAC設(shè)備維持空間單元一定舒適性的情況下，與NC模型組相比，BC模型組的P值減小為77.37%～86.00%，H值為93.81%～102.07%.H值除受圍護(hù)結(jié)構(gòu)濕流量大小的影響外，還受室內(nèi)濕負(fù)荷、建筑構(gòu)件濕流方向等多方面因素綜合影響，該模擬結(jié)果表明BC并不能對加濕和除濕量起明確作用.

圖7 外墻熱濕性能模擬結(jié)果Fig.7 Exterior wall hygrothermal performance simulation results

圖8 圍護(hù)空間室內(nèi)熱濕環(huán)境和年HVAC需求模擬結(jié)果Fig.8 Enclosed space indoor hygrothermal environment and annual HVAC demand simulation results

3 結(jié)論

(1)與3種竹板相比，在本文測試范圍內(nèi)竹炭平衡含水率在相對濕度為33.4%～85.4%時得到擴(kuò)大；其蒸汽滲透系數(shù)為9.32×10-11kg/(m·s·Pa)，是竹板的10.0～213.2倍，表明炭化處理使BC的調(diào)濕容量和調(diào)濕速率得到擴(kuò)大.

(2)建筑構(gòu)件和圍護(hù)空間熱濕性能模擬結(jié)果顯示，竹炭對構(gòu)造填充層含濕量起顯著穩(wěn)定作用，并可減小外墻熱濕流量和內(nèi)表面溫度振幅.竹炭有助于減小室內(nèi)空氣溫濕度振幅，縮短HVAC開啟時間，降低供暖和制冷量，但對加濕和除濕量影響不明顯.

致謝：浙江農(nóng)林大學(xué)工程學(xué)院張文標(biāo)教授、浙江笙炭控股有限公司陳俊輝先生為本研究工作提供竹炭樣品；Fraunhofer IBP(德國費(fèi)勞恩霍夫建筑物理研究所)提供WUFI Plus軟件使用權(quán)； Matthias Pazold先生、Manuel Lindauer博士和Simon Schmidt博士為模擬方案和軟件使用提供建議與技術(shù)支持，在此鳴謝.

參考文獻(xiàn)：

[1] YANG Li,LIU Hongbo,ZHANG Dongsheng,et al.Electron microscopy study on microstructure of bamboo charcoal[J].Journal of Chinese Electron Microscopy Society,2011，30(2):137-142.

[2] ISA S S M,RAMLI M M,HAMBALI N A M A,et al.Adsorption properties and potential applications of bamboo charcoal:A review[C]//MATEC Web of Conferences,2nd International Conference on Green Design and Manufacture 2016. Phuket,Thailand:EDP Sciences,2016,78:01097.

[3] WU G,ZHANG Q.Adsorptive removal of p-nitroaniline from aqueous solution by bamboo charcoal:Kinetics,isotherms,thermodynamics,and mechanisms[J].Desalination and Water Treatment,2016,57(55):26448-26460.

[4] FUJIWARA S,SHIMA K,CHIBA K.Fundamental characteristics and humidity control capacity of bamboo charcoal[J].Mokuzai Gakkaishi/Journal of the Japan Wood Research Society,2003,49(5):333-341.(in Japanese)

[5] LAOHHASURAYOTIN K,POOKBOONMEE S.Multifunctional properties of Ag/TiO2/bamboo charcoal composites:Preparation and examination through several characterization methods[J].Applied Surface Science,2013,282:236-244.

[6] LI Wenyan,YANG Hui,GUO Xingzhong,et al.Building decoration material containing charcoal prepared by using bamboo charcoal and diatomite as raw materials[J].Journal of Materials Science and Engineering,2011,29(1):7-11.

[7] 周建斌.竹炭環(huán)境效應(yīng)及作用機(jī)理研究[D].南京：南京林業(yè)大學(xué),2005.

ZHOU Jianbin.Study on mechanism of action and environmental application of bamboo charcoal[D].Nanjing:Nanjing Forestry University,2005.(in Chinese)

[8] 左宋林，高尚愚，徐柏森，等.炭化過程中竹材內(nèi)部形態(tài)結(jié)構(gòu)的變化[J].林產(chǎn)化學(xué)與工業(yè),2004,24(4):56-60.

ZUO Songlin,GAO Shangyu,XU Bosen,et al.Change of inner morphological texture of bamboo during carbonization[J].Chemistry and Industry of Forest Products，2004,24(4):56-60.(in Chinese)

[9] 蔣正武,張靜娟,曾志勇,等.炭化工藝參數(shù)對竹炭調(diào)濕性能的影響研究[J].建筑材料學(xué)報(bào),2017,20(5):717-722.

JIANG Zhengwu,ZHANG Jingjuan,ZENG Zhiyong,et al.Effect of carbonation process parameters on humidity-controlling properties of bamboo charcoal[J].Journal of Building Materials,2017,20(5):717-722.(in Chinese)

[10] 盧克陽,張齊生,蔣身學(xué).竹炭吸濕性能的初步研究[J].木材工業(yè),2006,20(3):20-22.

LU Keyang,ZHANG Qisheng,JIANG Shenxue.A primary study of bamboo charcoal moisture absorption[J].China Wood Industry,2006,20(3):20-22.(in Chinese)

[11] 西藤,宮野,田中.屋內(nèi)濕度變化と壁體材料.日本建築學(xué)會主編.日本建築學(xué)會研究報(bào)告:第 3 集[R].福岡:秀巧社印刷株式會社,1949:21.

[12] MAKI F,NORIMOTO M，YAMADA T.Humidity conditions in the house lined with wood based materials[J].Wood Research & Technical Notes,1985,21:87-95.

[13] IKEDA T.Moisture absorption and desorption characteristics and the effects of porous materials(Inorganic porous material)[J].Gypsum & Lime,1992,240:355-363.

[14] QIN Menghao,RAFIK B,ABDELKARIM A,et al.Simulation of coupled heat and moisture transfer in air-conditioned buildings[J].Automation in Construction,2009,18(5):624-631.

[15] 李國勝.海泡石礦物材料的顯微結(jié)構(gòu)與自調(diào)濕性能研究[D].天津:河北工業(yè)大學(xué),2005.

LI Guosheng.Microstructure and self-humidity-control performance of mineral material sepiolite[D].Tianjin:Hebei University of Technology,2005.(in Chinese).

[16] 侯國艷,冀志江,王靜,等.調(diào)濕材料的國內(nèi)外研究概況[J].材料導(dǎo)報(bào),2008,22(8):78-82.

HOU Guoyan,JI Zhijiang,WANG Jing,et al.Domestic and abroad research status of humidity-control materials[J].Materials Review,2008,22(8):78-82.(in Chinese).

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

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

[19] ISO 12572:2001(E) Hygrothermal performance of building materials and products-Determination of water vapour transmission properties[S].

[20] ASTM E 96/E 96M-2005 Standard Test Methods for Water Vapor Transmission of Materials[S].

[21] MUKHOPADHYAYA P,KUMARAN M K,LACKEY J.Use of the modified cup method to determine temperature dependency of water vapor transmission properties of building materials[J].Journal of Testing and Evaluation,2005,33(5):316-322.

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

[23] 孫立新,周輝,馮馳.風(fēng)驅(qū)雨對巖棉薄抹灰外墻外保溫系統(tǒng)含濕量模擬研究[C]//第十二屆全國建筑物理學(xué)術(shù)會議論文集.北京：中國建筑工業(yè)出版社,2016:49-52.

SUN Lixin,ZHOU Hui,FENG Chi.Research on the simulation of hygrothermal distribution of wind-driven rain on stone wool ETICS[C]//Proceedings of the 12th National Conference on Building Physics.Beijing:China Architecture & Building Press,2016:49-52.(in Chinese).