北美典型气候区建筑围护结构HM模拟及分析

黄祖坚 孙一民

摘   要：为明确气候参数、材料参数对建筑围护结构模拟结果的影响，以及热湿耦合传递模型（HM）和被普遍采用的Glaser模型之间的差异，在热湿过程模拟软件WUFI Plus中选择北美典型气候区的20个代表城市，设置15组外墙构件，进行空间单元HVAC需求、室内热湿环境和外墙热湿性能的全年模拟.模拟结果与室外气候参数进行相关性分析，表明室外气温、太阳辐射、相对湿度和风驱雨不同程度地影响相关结果.通过材料参数控制所得2组HM和3组Glaser模型的对比，表明液态水相关参数、含湿量引起的材料参数的变物性取值对模拟结果有显著影响.北美典型气候区HM模拟对我国建筑围护结构热湿过程相关的气候参数、材料参数和评估方法的研究具有借鉴意义.

关键词：典型气候区;建筑围护结构;热湿耦合传递模型;Glaser模型

中图分类号：TU111.19                          文献标志码：A

Abstract：To determine the impact of climate data and material parameter on building envelope simulation results，and to investigate the difference between coupled heat and moisture transfer model （HM model） and the widely used Glaser model，15 groups of exterior walls were constructed in WUFI Plus with 20 representative cities of the North America typical climate zones as external conditions，and simulated for the annual HVAC demand，indoor hygrothermal environment and the exterior walls hygrothermal performance. Correlation analyses between the simulation results and the climate parameters showed that the exterior air temperature，solar radiation，relative humidity and drying rain affected the related simulation results to varying degrees. Comparison among the 2 HM and 3 Glaser model groups distinguished with specific material parameters，showed that liquid water dependent parameters，variation of moisture content dependent heat and moisture transport properties had significant influence on the simulation results. The HM model simulation in North America typical climate zones could contribute implications to the studies on climate data，material parameters and evaluation methods for the heat and moisture process in building envelope.

Key words：typical climate zone;building envelope;coupled heat and moisture transfer model;Glaser model

準确理解建筑围护结构的热湿过程，并进行合理控制，有助于延长建筑构件寿命[1]、降低建筑运行能耗[2]、提高室内环境舒适性与卫生质量[3-4].在建筑热工学领域通常采用热湿过程模型进行计算，在1990年以前由Glaser提出的稳态蒸汽渗透模型，采用蒸汽压作为驱动势，基于Fick定律确定纯蒸汽的湿分迁移速率[5]，这一方法操作性强，可以对建筑构件进行初步的湿性能评价，预测冷凝部位和冷凝程度，因而被广泛使用，如EN ISO 13788、DIN 4108、GB 50176-93/2016[6].然而实际的湿过程常常是非稳态、多相同时发生的，对建筑围护结构性能的精确认识要求从Glaser的简单评估转变为对湿过程的仿真模拟.

1975年Luikov基于Fourier 定律、Fick 定律和Darcy 定律，根据能量守恒定律、质量定律和动量守恒定律建立起控制方程，从机理上较真实地描述湿分在多孔介质中的迁移过程[7].近20年，在Luikov热湿耦合传递理论的基础上，Molenda、Pedersen、Künzel、H？覿upl、Janssen、Mendes、Steeman和Tariku等人相继提出分析建筑围护结构中热量和湿分耦合传递模型（Coupled heat and moisture transfer model，简称HM模型）[8-15].这些模型已经得到不同角度的验证，2008年完成的IEA Annex 41项目对17种HM模型及数值模拟工具进行详细比较，并完成7项联合对比任务，证实主要HM模型的计算结果、及其与实测结果之间非常接近[16-17].HM模型基础理论在我国也得到许多学者的关注，陈友明、郭兴国、刘向伟、王莹莹、孔凡红、郑茂余、张华玲等学者从不同角度开展HM模型研究[18-27]，其中湖南大学的李念平团队推导竹胶合板外墙热湿过程模型并进行实测检验[28].

HM模型通常采用高度耦合的非线性偏微分方程来同时描述建筑围护结构中的热湿传递过程[29].

由于运算过程中需要不断地通过有限差分法对方程组进行求解，并通过积分更新构造内热湿场分布，其计算时间可达普通能耗模型的102～104倍.计算机工具是HM模型投入实际应用的前提条件，已有部分计算机软件与HM模型结合进行特定分析，如BES-HAM的建筑能耗模拟分析，WUFI系列软件的建筑热湿环境、建筑构造热湿传递和霉菌增长分析等[29-30].常用建筑热湿过程模拟软件还有Delphin（德国）、CHAMPS（美国），部分综合能耗模拟软件如EnergyPlus（美国）、TRNSYS（美国&欧洲）也增加HM功能.计算机工具的支持使未来以HM模型的模拟仿真取代Glaser模型的简单评估具有可行性.

在实验室和实践工程中，HM模型及其计算机工具体现出的可信度为其赢来越来越多从业者的接受，DIN4108—3标准已经承认这些方法[10].相比之下，我国目前GB 50176—2016《民用建筑热工设计规范》对建筑围护结构的计算偏重于热过程，而湿过程仍基于Glaser模型的简单评估，与此相应的气候和材料研究也相对滞后.在建筑气候研究上，采用气温和太阳辐射为分区依据，未给予环境湿度和降雨参数足够重视;在材料数据库上，目前对常用建筑材料的参数积累只能用于Glaser模型的简单计算，无法支撑完整的HM模型.

为明确气候参数、材料参数对建筑围护结构模拟结果的影响，以及HM模型和Glaser模型之间的差异，本文采用德国Fraunhofer IBP开发的热湿过程模拟软件WUFI Plus，选择气候参数完善的北美典型气候区进行模拟和分析，为我国建筑围护结构热湿过程相关的气候参数、材料参数和评估方法的研究提供参考.

1   基于Künzel方程的HM模型

HM模型与Glaser模型存在许多不同.HM模拟采用热湿耦合的方式进行，而非让热过程和湿过程分离开来计算;除气温、相对湿度、太阳辐射外，HM气候参数还包括降雨因子;除表观密度、比热容、导热系数、蒸汽渗透阻力因子外，HM模型还需要等温吸放湿曲线、孔隙率、液态水传递系数等材料参数;考虑多相湿过程，以及材料湿储存性质（含湿量）对热传递性质、液态水和气态水传递性质的影响.

不同学者在研究HM模型中，由于基础描述方程的差异，需要不同的材料参数作为输入.而对于同一材料参数，由于对驱动势选择和取值方法等方面的差异，以及对环境因素、材料含湿量等影响的不同判断，也会有不同测试要求.采用温度作为热量传递的驱动势得到比较一致认可，但对于湿分传递的驱动势则没有形成统一的观点.常用的湿分传递驱动势有温度、含湿量、蒸汽压和抽吸应力.

1995年，德国建筑物理学家H M Künzel通过理论推算和实测检验将热湿传递计算所需材料参数进行简化，提出以蒸汽压和相对湿度作为湿传递驱动势的改进的Luikov模型，将热量和湿分过程描述为两个耦合微分方程[10].两等式左侧由储存部分构成：热储存包括干燥材料的热容和材料中湿分的热容，湿储存根据材料等温吸放湿曲线来描述.等式右侧由传递部分构成：热传递根据含湿量相关的导热系数和蒸汽焓流来描述，后者是由于湿分从一处蒸发吸热，然后扩散到另一处冷凝放热的潜热效果.液态水传递包含受相对湿度差驱动的表面扩散和毛细传导，受温度影响相对较小.气态水扩散则受温度场强烈影响.以Künzel方程为基础，德国Fraunhofer IBP开发出热湿过程模拟软件WUFI Plus（W？覿rme- Und Feuchte- transport Instation？覿r非稳态瞬时热量和水分传输）.如图1所示.

2   模拟方案

WUFI Plus模型的运行需要设置外部条件、内部条件、边界条件、HVAC条件，和数据网格、时间步长等控制参数，如图2所示.

1）外部条件.根据K？觟ppen的气候区划分方法，北美气候区类型丰富，对我国典型气候区具有较为完整的参考价值.选择北美15种典型气候区的20个代表城市的气候数据作为模型的外部条件，如图3所示.

2）内部条件.建立尺寸为 3.0 m×3.0 m×3.0 m（开间×进深×层高）的空间单元.设置9：00-17：00为房间工作时间，并赋予标准办公间使用强度的室内热湿负荷，其中对流热：33.3 W，辐射热：25.2 W，湿分：17.55 g/h，CO2发生量：20.79g /h，人员活动量：1.2 met.

3）边界条件.选择四面外墙作为研究对象，面积均为9 m2，将顶面天花及底面地板设置为相同内部条件的两房间之间的分隔（space with the same interior conditions），以避免頂、底面的热湿传递干扰.外墙构造上类型上设置常用的3种类型，包括L型（外面板+填充层+内面板）、H型（砌体层+内面板）和M型（砌体层 + 填充层 + 内面板）.构造层厚度设置参考北美对相应气候区的U值要求，分为0.15/0.2/0.4/0.6/1.0 W/（m2·K）的5个组，如表1所示.

4）HVAC条件.基于以上设置的外部、内部和边界条件，首先关闭HVAC，模拟在无设备条件下室内热湿环境和外墙热湿性能;然后打开HVAC，以理想的供暖和制冷、加湿和除湿设备维持一定室内温湿度，输出空间单元的供暖和制冷量，以及加湿和除湿量.运行时间2015-01-01至2017-01-01，对2016-01-01 至2017-01-01进行数据收集，时间步长为1 h.

通过材料参数调整，设置2组HM模型（HM和HM0）和3组Glaser模型（Glaser0，Glaser50和Glaser75）.其中HM组采用完整的材料参数，而HM0在HM材料参数的基础上将其毛细吸水速率和液态水传递速率设为0，并将模型外墙降雨负荷关闭.Glaser0、Glaser50和Glaser75分别采用恒定对应RH=0%，RH=50%和RH=75%的热传递和气态水传递性质的参数，且将湿储存性质和液态水相关性质设为0，从而切断热湿过程之间的耦合关系，使模型热过程和湿过程模拟相互独立开.

3   结果分析

3.1   气候参数相关性分析

选择开启理想HVAC情况下年供暖&制冷量P值、供暖峰值Pheating.max和制冷峰值Pcooling.max、加湿&除湿量H值、加湿峰值Hhumid.max和除湿峰值Hdehumid.max表征空间单元HVAC需求;关闭HVAC情况下室内年温、湿度均值Timean和RHimean及相应波幅Tiamp和RHiamp表征室内热湿环境;通过建筑四面外墙的年热、湿流量Hflow和Mflow及相应峰值Hflow.max和Mflow.max表征外墙构造热湿性能.选择室外年空气温度Te、太阳辐射总量SR、年空气相对湿度RHe和年风驱雨量DR作为室外气候参数.将模拟结果与气候参数进行相关性分析.

考虑到模拟结果和气候参数之间不符合双变量正态分布，且变量之间无线性相关，因而选择非参数统计方法的Spearman相关系数进行分析，借助统计分析软件SPSS进行操作.将相关指标分为两组，组1为模拟结果和气候参数之间总量和均值指标的相关性分析，组2则针对波幅和峰值（注：由于所有城市气候参数中年室外相对湿度最大值均为100%，年风驱雨最小值均为0，因此相应地仅采用年室外相对湿度最小值和年风驱雨最大值进行分析）.

总量和均值相关性分析结果表明，与热过程相关的指标中，Timean、P值和Hflow均与Temean和SR存在0.01层上显著的相关关系，且Timean和Hflow分别和DR在0.05和0.01层上显著相关，可见建筑围护结构热过程不仅取决于热气候参数，也受湿气候参数的影响;与湿过程相关的指标中，RHimean和RHemean存在0.01层上显著的相关关系;Mflow和DR在0.05层上显著相关;而H值则没有表现出与室外气象参数的明显相关关系.波幅和峰值相关性分析结果表明，Tiamp和RHiamp和室外空气温度振幅均存在0.01层上显著的相关关系，且RHiamp还与室外相对湿度和年风驱雨峰值在0.01层上显著相关;在HVAC峰值相关指标中，除年供暖峰值外，制冷、加湿和除湿峰值分别与年均室外气温和年风驱雨峰值存在显著相关关系;外墙Hflow.max和Mflow.max分别和年室外气温和年风驱雨峰值之间在0.01和0.05层上显著相关，如表2所示.

同一热分区，不同湿分区的不同模型组之间模拟结果存在显著差异，室外相对湿度和降雨等湿气候参数的影响不可忽略.因此依据气温、太阳辐射强度进行单纯的热工分区，以及单一的分区传热系数U值规定作为建筑围护结构设计依据存在一定缺陷.美国的建筑气候分区方法采用气温和太阳辐射作为一级区划指标，其下根据相对湿度和降雨量又作二级区划，以此保证热气候参数为主，又兼顾湿气候参数作为设计依据，这对我国建筑气候区划方法有参考意义，如图4所示.

3.2   材料参数对模拟结果的影响

3.2.1   液态水对模拟结果的影响

以HM0模型组（不考虑风驱雨）和HM模型组（考虑风驱雨）模拟结果的比值作为指标，对南面外墙L类构造和M类构造的填充层含湿量，以及H类构造的砌体层含湿量进行比较.HM0 /HM模型组计算结果比值最小分别可达92.84%、31.93%和5.86%.可见液态水传递对模拟结果有不可忽略的影响，通常采用的Glaser评估方法未考虑降雨因素作用，会使构造含湿量被不同程度地低估，如图5所示.

HM0/HM模型组模拟结果对比表明液态水对HVAC需求、室内热湿环境、建筑构件热湿性能计算结果有明显影响.对于L类构造，由于外面板的液态水传递阻力较大，液态水对构造整体及所围护空间的影响受到阻挡，此时HM0/HM模型组模拟结果各指标的比值在98.85%～100.88%之间.对于H类构造，由于构造主体层直接与外部雨水接触，风驱雨对构造影响最大，HM0 /HM模型组模拟结果比值中，有：HVAC相关的P：78.07% ～ 100.00%，Pheating.max：78.77%～ 100.00%，Pcooling.max ： 95.02%～ 125.95%;

Hdehumid.max：51.91% ～ 06.83%;室内热湿环境相关的Timean：100.00%～ 104.92%，Tiamp：96.43% ～ 102.75%，RHimean：94.01% ～ 100.00%，RHiamp：91.61% ～ 106.57%;外墙热湿性能相关的Hflow：94.95%～100.07%，Mflow：92.08%～ 100.00%.M类构造受液态水影响介于L和H之间，相应各指标的比值在59.78%～108.19%之间，如图6所示.

3.2.2   導热系数变物性取值对计算结果的影响

以HM0模型组为基准，将年供暖&制冷量及其峰值、年室内空气平均温度及其振幅和外墙年热流量的Glaser/HM0模拟结果比值r-P、r-Pheating.max、r-Pcooling.max、r-Timean、r-Tiamp和r-Hflow作为指标，相应的值为r-P：78.10%～116.35%、r-Pheating.max：44.02%～106.75%、r- Pcooling.max：52.73%～ 146.66%、r- Timean：100.00%～162.50%、r-Tiamp：78.52%～102.06%和r-Hflow：97.06%～117.38%.该值在3类构造中偏幅不同，其大小与构造自身受含湿量影响程度有关，除r-Pcooling.max外，H类构造偏幅最大，如图7所示.

对Glaser0/50/75 3组计算结果进行差异性分析，表明由含湿量提高引起的导热系数的增大，使计算所得P值和Hflow值增大.值得一提的是对于作为常规建筑热工计算的主要指标P值，除H类构造的个别案例，HM0模型组模拟结果普遍落在Glaser50和Glaser75模型组所得结果之间，这意味着采用室内和室外空气平均相对湿度对应的材料导热系数作简化计算，对P值估算有较好的参考价值.

3.2.3   蒸汽渗透阻力因子变物性取值对模拟结果的影响

以HM0模型组为基准，将年加湿&除湿量H值及其峰值、年室内空气平均相对湿度及其峰值和外墙年湿流量的Glaser/HM0比值r-H、r-Hhumid.max、r-Hdehumid.max、r-RHimean、r-RHiamp和r-Mflow作为指标，有r-H：59.79%～121.91%、r-Hhumid.max：24.61%～199.37%、r-Hdehumid.max：0.00%～115.78%、r-RHimean：92.78%～100.45%、r-RHiamp：76.86%～115.70%和r-Mflow：81.72%～143.71%，如图8所示.

對Glaser0/50/75 3组计算结果进行差异性分析，表明由含湿量提高引起的蒸汽渗透阻力的减小，使建筑构件年湿流量Mflow增大，其中Glaser75/ Glaser0的比值达到97.50%～139.58%;RHimean值主要取决于室外空气相对湿度，受构件自身湿传递性质影响不明显，Glaser50和Glaser75对Glaser0的RHimean比值在94.56%～100.60%;年加湿&除湿量H值除受建筑构件湿传递性质外，还受室内外湿负荷、建筑构件湿流方向等多因素综合作用，因此蒸汽渗透系数的变物性取值并未对H值模拟结果产生明确影响.

与液态水引起的结果差异的对比显示，在H类构造中，含湿量引起的导热系数和蒸汽渗透系数变物性取值所带来的模拟结果差异，大幅小于液态水的影响.然而在L和M类构造中，相应的差异则不可忽略.

3.3   长周期动态的湿过程评估方法

Glaser模型对建筑围护结构湿性能采用的是稳态的评估方法，即假设建筑构件内外表面具有一定的蒸汽压差驱动通过建筑构件的稳定蒸汽流.在HM的全年模拟中，受外部气候条件、内部热湿负荷、HVAC条件等因素影响，全年动态模拟显现出和建筑构件稳态评估不同的结果.建筑构件层面稳态评估的各项指标不能直接、准确地反映全年HVAC需求、室内湿环境和建筑构件湿性能.相比于热过程，湿过程模拟结果与外部、内部、边界条件之间相关性更弱，因而更加难以进行预判，因此对建筑围护结构模拟实际使用条件、长周期、动态的评估具有必要性.

4   结 论

目前我国建筑围护结构研究偏重单纯的热过程，湿过程仅基于Glaser模型进行简单评估，本文对北美典型气候区15组外墙构件进行HM模拟，明确气候参数、材料参数对模拟结果的影响，以及HM模型和Glaser模型之间的差异.模拟和分析结果对我国建筑围护结构相关的气候参数、材料参数和评估方法的研究具有参考意义.

1）气候参数的影响.相关性分析表明环境相对湿度和降雨因素对室内热湿环境和建筑热湿性能有显著影响，未考虑风驱雨参数的模拟会造成对构造含湿量、年能耗量和构造热湿流量的低估，对于主体层外露的构造类型，液态水影响大幅高于气态水.

2）材料参数的影响.由构造层含湿量引起的导热系数变物性取值对年能耗量、室内热环境和建筑构件热流量有显著影响;蒸汽渗透阻力因子变物性取值明显影响构件湿传递，但对室内湿环境和年加湿&除湿量的影响不明确.

3）受气候条件、构造类型等因素影响，全年动态模拟显现出和建筑构件稳态评估不同的结果，对建筑围护结构模拟实际使用条件下长周期、动态的评估具有必要性.

致谢   Fraunhofer IBP（德国费劳恩霍夫建筑物理研究所）为本研究提供WUFI Plus计算机程序使用权; Matthias Pazold先生、Manuel Lindauer博士和Simon Schmidt博士为模拟方案和软件使用提供建议与技术支持.在此鸣谢.

参考文献

[1]   GEVING S，HOLME J. The drying potential and risk for mold growth in compact wood frame roofs with built-in moisture[J].Journal of Building Physics，2010，33（3）：249—269.

[2]   MENDES N，WINKELMANN F C，LAMBERTS R，et al. Moisture effects on conduction loads[J].Energy and Buildings，2003，35（7）：631—644.

[3]   ZHANG H，YOSHINO H，HASEGAWA K.Assessing the moisture buffering performance of hygroscopic material by using experimental method[J].Building and Environment，2012，48：27—34.

[4]  ABUKU M，JANSSEN H，ROELS S.Impact of wind-driven rain on historic brick wall buildings in a moderately cold and humid climate：Numerical analyses of mould growth risk，indoor climate and energy consumption[J]. Energy and Buildings，2009，41（1）：101—110.

[5]   GLASER H.Graphical method for investigation of diffusional process[J].Kalteteohnik，1959，39（5）：519—550.

[6]   GB 50176-93/2016民用建筑熱工设计规范[S]. 北京：中国建筑工业出版社，2016：33—37.

GB 50176-93/2016. Thermal design code for civil building[S].Beijing：China Architecture & Building Press，2016：33—37.（In Chinese）

[7]   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（6）：1—14.

[8]    MOLENDA C H A，CRAUSSE P，LEMARCHAND D.The influence of capillary hysteresis effects on the humidity and heat coupled transfer in a non-saturated porous medium[J].International Journal of Heat and Mass Transfer，1992，35（6）：1385—1396.

[9]  PEDERSEN C R. Prediction of moisture transfer in building constructions[J]. Building and Environment，1992，27（3）：387—397.

[10] K？譈NZEL H M. Simultaneous heat and moisture transport in building components [R]. Fraunhofer IRB Verlag Suttgart，1995：38—46.

[11]  H？魧UPL P， GRUNEWALD J， FECHNER H，et al. Coupled heat air and moisture transfer in building structures [J]. International Journal of Heat and Mass Transfer，1997，40（7）：1633—1642.

[12] JANSSEN H， BLOCKEN B， CARMELIET J. Conservative modelling of the moisture and heat transfer in building components under atmospheric excitation[J].International Journal of Heat and Mass Transfer，2007， 50（5/6）： 1128—1140.

[13] MENDES N，PHILIPPI P C. A method for predicting heat and moisture transfer through multilayered walls based on temperature and moisture content gradients [J]. International Journal of Heat and Mass Transfer， 2005， 48（1）： 37—51.

[14] STEEMAN H J， VAN BELLEGHEM M， JANSSENS A， et al. Coupled simulation of heat and moisture transport in air and porous materials for the assessment of moisture related damage[J]. Building and Environment，2009， 44（10）： 2176—2184.

[15] TARIKU F，KUMARAN K，FAZIO P. Transient model for coupled heat， air and moisture transfer through multilayered porous media[J]. International Journal of Heat and Mass Transfer，2010， 53（15/16）： 3035—3044.

[16] KUMARAN M K. IEA Annex 24： heat， air and moisture transfer in insulated envelope parts. Final Report， Volume 3， Task 3： Material properties[R]. Brussels： Laboratory of Building Physics， Civil Engineering Department，KU Leuven，Beligum，1996：14—132.

[17] WOLOSZYN M， RODE C. IEA Annex 41： Whole building heat， air， moisture response. Subtask 1： Modeling principles and common exercises[R].La Rochelle： International Energy Agency，Executive committee on Energy Conservation in Buildings and Community Systems，2008：136—239.

[18]  陳友明，邓永强，郭兴国，等.建筑围护结构热湿耦合传递特性实验研究与分析[J].湖南大学学报（自然科学版）， 2010，37（4）：11—16.

CHEN Y M， DENG Y Q，GUO X G，et al. Experimental study of coupled heat and moisture transfer performance through building envelope[J]. Journal of Hunan University （Natural Sciences），2010，37（4）：11—16.（In Chinese）

[19]  郭兴国.热湿气候地区多层墙体热湿耦合迁移特性研究[D].长沙：湖南大学土木工程学院，2010：18—44.

GUO X G.Research on coupled heat and moisture transfer characteristics of multilayer walls in hot and humid climate[D].Changsha： College of Civil Engineering， Hunan University， 2010：18—44.（In Chinese）

[20]  刘向伟.夏热冬冷地区建筑墙体热、空气、湿耦合迁移特性研究[D].长沙：湖南大学土木工程学院，2015：18—82.

LIU X W. Investigation of the coupled heat， air and moisture transport in building walls in hot summer and cold winter zone[D].Changsha： College of Civil Engineering， Hunan University， 2015：18—82.（In Chinese）

[21]  刘向伟，陈国杰，陈友明.墙体热、湿及空气耦合传递非稳态模型及验证 [J].湖南大学学报（自然科学版），2016，43（1）：152—156.

LIU X W，CHEN G J，CHEN Y M. Modeling of the transient heat，air and moisture transfer in building walls[J].Journal of Hunan University （Natural Sciences），2016，43（1）：152—156.（In Chinese）

[22]  LIU X W，CHEN Y M，GE H.Numerical investigation for thermal performance of exterior walls of residential buildings with moisture transfer in hot summer and cold winter zone of China[J].Energy and Buildings， 2015（93）：259—268.

[23] LIU X W，CHEN Y M，GE H.Determination of optimum insulation thickness for building walls with moisture transfer in hot summer and cold winter zone of China[J].Energy and Buildings， 2015（109）：361—368.

[24]  王莹莹.围护结构湿迁移对室内热环境及空调负荷影响关系研究[D].西安：西安建筑科技大学环境与市政工程学院，2013：26—42.

WANG Y Y.Research on the effect of the palisade structure moisture transfer on the indoor thermal environment and air-conditioning load[D]. Xi′an： School of Environmental and Municipal Engineering， Xi′an University of Architecture and Technology，2013：26—42.（In Chinese）

[25]  孔凡红，郑茂余，韩宗伟，等.新建建筑围护结构热质传递对室内温湿度环境的影响[J].建筑科学，2008，24（8）：94-97.

KONG F H， ZHENG M Y， HAN Z W，et al. Influences of envelope heat and moisture transfer of new building on indoor temperature and humidity[J].Building Science，2008，24（8）：94—97.（In Chinese）

[26]  郑茂余，孔凡红，韩宗伟.新建建筑外保温围护结构热质耦合传递[J]. 哈尔滨工业大学学报，2009，41（4）：118—122.

ZHENG M Y，KONG F H，HAN Z W.Simulations on heat and mass coupling transfer in exterior insulated envelope of new building[J]. Journal of Harbin Institute of Technology，2009，41（4）： 118—122.（In Chinese）

[27]  张华玲，刘朝，刘方，等.地下洞室多孔墙体热湿传递的数值模拟[J].暖通空调，2006，36（12）：9—13.

ZHANG H L，LIU C， LIU F，et al.Numerical simulation of heat and moisture transfer in porous wall in underground caves[J].Heating Ventilating & Air Conditioning，2006，36（12）： 9—13.（In Chinese）

[28]  王倩，李念平，曾德军，等.竹材组合墙体构件热湿耦合迁移特性[J].科技导报，2010，28（17）：87—90.

WANG Q，LI N P，ZENG D J，et al. Coupled heat and moisture transfer in bamboo composite wall [J].Science & Technology Review，2010，28（17）：87—90.（In Chinese）

[29] VAN BELLEGHEM M， STEEMAN H J， STEEMAN M，et al. Sensitivity analysis of CFD coupled non-isothermal heat and moisture modelling[J].Building and Environment，2010， 45（11）：2485—2496.

[30] NZEL H M. Verfahren zur ein-und zweidimensionalen berechnung des gekoppelten rme- und feuchtetransports in bauteilen mit einfachen kennwerten[D]. Stuttgart： Lehrstuhl für Konstruktive Bauphysik，Dissertation Universit t Stuttgart， 1994：40—58.