CN113468638B - 一种基于热湿耦合的装配式建筑热桥分析方法 - Google Patents

一种基于热湿耦合的装配式建筑热桥分析方法 Download PDF

Info

Publication number
CN113468638B
CN113468638B CN202110685229.3A CN202110685229A CN113468638B CN 113468638 B CN113468638 B CN 113468638B CN 202110685229 A CN202110685229 A CN 202110685229A CN 113468638 B CN113468638 B CN 113468638B
Authority
CN
China
Prior art keywords
heat
thermal
building
method based
analysis method
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110685229.3A
Other languages
English (en)
Other versions
CN113468638A (zh
Inventor
郭兴国
季姜君
刘向伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanchang University
Original Assignee
Nanchang University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanchang University filed Critical Nanchang University
Priority to CN202110685229.3A priority Critical patent/CN113468638B/zh
Publication of CN113468638A publication Critical patent/CN113468638A/zh
Application granted granted Critical
Publication of CN113468638B publication Critical patent/CN113468638B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Architecture (AREA)
  • Civil Engineering (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

本发明的一种基于热湿耦合的装配式建筑热桥分析方法,主要考虑了装配式建筑外墙连接件部位热传递、湿迁移及其耦合作用的问题,通过COMSOL用有限元方法对模型进行数值求解,对装配式建筑热桥进行合理的分析,提高准确度,为装配式建筑节能效率不达标的问题提供了解答。本发明可为装配式建筑热桥部位的保温节能措施提供理论指导,有利于提高资源利用效率。适用于各种装配式建筑的热桥分析。

Description

一种基于热湿耦合的装配式建筑热桥分析方法
技术领域
本发明涉及建筑工程技术领域,更具体的,涉及一种基于热湿耦合的装配式建筑热桥分析方法。
背景技术
目前,装配式建筑在节能设计方面存在的问题主要可归结于以下几个方面:一是我国关于装配式建筑施工技术标准以及技术指导方面存在缺口,使得设计人员在进行节能设计时过多的参考普通建筑的保温技术设计规定与标准,造成了装配式建筑在施工完成后出现节能效率不达标和建筑材料大量浪费的后果;二是装配式建筑在现场施工和设计当中处理欠妥,如构造节点保温处理不当等,导致建筑产生大量热桥现象从而造成能耗增大。
除此之外,围护结构暴露在室内外热湿环境中,在温度和水蒸气分压力梯度的共同作用下,室内外的水分会通过围护结构发生迁移。围护结构各层之间的水分迁移不仅会通过降解材料和增加霉菌生长的风险而缩短系统的使用寿命,还会对室内空气湿度,热负荷和湿负荷产生重大影响。多孔材料(如外墙面板的隔热层和结构层)中的蒸汽凝结也会提高整体传热。
发明内容
本发明的目的在于提出一种基于热湿耦合的装配式建筑热桥分析方法,在各种复杂环境条件下,对不同装配式建筑连接节点热桥进行热湿分析,为装配式建筑保温设计提供理论指导。
本发明是通过以下技术方案实现的:
本发明所述的一种基于热湿耦合的装配式建筑热桥分析方法,其特征在于包括以下步骤:
步骤(1)获取墙体材料的性能参数,包括密度(kg/m3)、恒压热容(J/(kg·K))、导热系数(W/(m·K))、含水量(kg/m3)、扩散系数(m2/s)、蒸汽渗透性(kg/(m·s·Pa));
(以上材料参数可由查找文献或实验得出)
步骤(2)查询当地典型年份气象参数,包括月平均温度Text(℃)、月平均相对湿度φext、月平均风速u(m/s);
步骤(3)在COMSOL Multiphysics选择三维建筑材料热湿传递模型,其湿控制方程和热控制方程分别为
Figure GDA0003668210060000021
其中(ρc)eff=ρmcp,m+ωcp,l(方程3)
式中ξ为湿容(等温吸放曲线的斜率)(kg/m3),φ为空气相对湿度(%),Dw为液态水扩散系数(m2/s),δP为水蒸气渗透率(kg/(m·s·Pa)),Psat为饱和蒸汽压力(Pa),keff为当量导热系数,Lv为水蒸气的汽化潜热(J/kg),ρm为干材料密度(kg/m3),cp,m为干材料的比热(J/(kg·K)),ω为含湿量(kg/m3),cp,l为液态水的比热(J/(kg·K)),T为墙体温度(K),t为时间(s),
Figure GDA0003668210060000022
为梯度算子;
步骤(4)导入三维物理模型;
步骤(5)输入参数,包括气象参数以及材料参数;
步骤(6)根据已知条件为各个区域选择合适的初始条件与边界条件,边界条件包括热边界条件与湿边界条件,用户可根据已知条件选择温度、热通量、热绝缘、湿度、水汽通量、防潮层边界条件;
步骤(7)划分网格,在网格划分方面COMSOL自带预定义的自由四面体网格单元划分方式,用户只需选择合理的网格尺寸(极端粗化、特别粗化、较粗化、粗化、正常、细化、较细化、特别细化、极端细化)对研究对象进行网格划分,网格尺寸越密,计算越精确,但所需计算时间越长,设定求解方式和时间步长;
步骤(8)计算结果;
步骤(9)选择几种不同的装配式建筑连接方式以及材料,根据各个月份的气象参数分别计算每个月的墙体温度与相对湿度参数。
进一步地,所述气象参数包括风速、温度、湿度、太阳辐射照度、天空有效温度、降雨以及盐雾浓度中的至少一种。
进一步地,所述网格尺寸包括至少对以下极端粗化、特别粗化、较粗化、粗化、正常、细化、较细化、特别细化、极端细化中的一种对研究对象进行网格划分。
进一步地,所述墙体材料可由MatWeb材料库导入。
进一步地,所述三维物理模型可由SolidWorks或CAD导入。
本发明的有益效果为:在各种复杂环境条件下,对不同装配式建筑连接节点热桥进行热湿分析,为装配式建筑热桥部位的保温节能措施提供了理论指导,有利于提高资源利用效率。
附图说明
图1装配式外围护结构组件和钢连接的热物理性质,
图2装配式外围护结构3D模型;
图3装配式外围护结构边界条件;
图4装配式外围护结构网格划分;
图5求解器设置;
图6连接件热桥截面的温度分布[℃];
图7连接件热桥截面的热通量分布[W/m2];
图8连接件热桥截面的相对湿度分布[%];
图9连接件热桥截面的水汽通量分布[kg/(m2·s)];
图10装配式建筑连接件热桥三视图及3D模型;
图11COMSOL求解流程图。
具体实施方式
以下结构附图对本发明的实施例进行详细说明。但本发明可以根据权利要求限定和覆盖的多种不同方式实施。
参考图11基于热湿耦合的装配式建筑热桥分析方法,以下为具体实施例:
本例以南昌地区的螺栓连接装配式建筑围护结构为例,某墙体为加气混凝土,保温层为XPS(挤塑式聚苯乙烯隔热保温板),墙面抹灰材料为的水泥砂浆,楼板为钢筋混凝土。
获取墙体材料的性能参数。
表1装配式外围护结构组件和钢连接的热物理性质
Figure GDA0003668210060000041
Figure GDA0003668210060000051
注意:f(x)表示该参数为参数x的函数,如图1所示。
以南昌地区二月份为例:
月平均温度Text为6.8℃,月平均相对湿度φext为81.5%,月平均风速2.4m/s,取室内设计温度Tint为21℃,相对湿度φint为60%。
在COMSOL Multiphysics选择三维建筑材料热湿传递模型,其湿控制方程和热控制方程分别为
Figure GDA0003668210060000052
Figure GDA0003668210060000053
其中(ρc)eff=ρmcp,m+ωcp,l(方程3)
式中ξ为湿容(等温吸放曲线的斜率)(kg/m3),φ为空气相对湿度(%),Dw为液态水扩散系数(m2/s),δP为水蒸气渗透率(kg/(m·s·Pa)),Psat为饱和蒸汽压力(Pa)。keff为当量导热系数,Lv为水蒸气的汽化潜热(J/kg),ρm为干材料密度(kg/m3),cp,m为干材料的比热(J/(kg·K)),ω为含湿量(kg/m3),cp,l为液态水的比热(J/(kg·K))。
使用SketchUp或其他3D建模软件对连接件热桥部位进行建模并导入COMSOL,如图2所示。
根据设计条件为各个区域选择材料并定义初始条件,如图3所示。
其中外表面边界条件为
Figure GDA0003668210060000061
内表面边界条件为
Figure GDA0003668210060000062
所有面板材料的初始条件设置为相对湿度60%,温度21℃(294.15K)。
使用COMSOL自动网格生成器生成了非结构化非均匀的四面体网格,选择网格尺寸为细化,如图4所示。为实现精确且独立于网格的模拟,将细孔(约10毫米)分布在钢连接附近的区域中,在这些区域中热量和水分传递的潜在变化可能性很大。大网格(约80毫米)分布在非关键区域,以使总网格数保持在合理的范围内。
本模拟案例使用的是月平均温湿度,因此选择基于有限元方法(FEM)的稳态求解器来数值求解耦合的控制方程(方程式1和方程式2)。采用全耦合方法来迭代求解稳态偏微分方程(PDE),设置如图5所示。设置完成后进行计算。
计算完成后从COMSOL结果方案中提取模型中的温度,热通量,相对湿度和水汽通量等,分别如图6所示、如图7所示、如图8所示、如图9所示。
同理可计算其他月份或其他连接方式的装配式建筑热桥的热湿性能,对不同装配式建筑连接节点热桥进行热湿分析。

Claims (5)

1.一种基于热湿耦合的装配式建筑热桥分析方法,其特征在于包括以下步骤:
步骤(1)获取墙体材料的性能参数,包括密度(kg/m3)、定压热容(J/(kg·K))、导热系数(W/(m·K))、含水量(kg/m3)、扩散系数(m2/s)、蒸汽渗透性(kg/(m·s·Pa));
步骤(2)查询当地典型年份气象参数,包括月平均温度Text(℃)、月平均相对湿度
Figure FDA0003726373600000011
月平均风速u(m/s);
步骤(3)在COMSOL Multiphysics选择三维建筑材料热湿传递模型,其湿控制方程和热控制方程分别为
Figure FDA0003726373600000012
Figure FDA0003726373600000013
其中(ρc)eff=ρmcp,m+ωcp,l (方程3)
式中ξ为湿容(kg/m3),φ为空气相对湿度(%),Dw为液态水扩散系数(m2/s),δP为水蒸气渗透率(kg/(m·s·Pa)),Psat为饱和蒸汽压力(Pa),keff为当量导热系数,Lv为水蒸气的汽化潜热(J/kg),ρm为干材料密度(kg/m3),cp,m为干材料的比热(J/(kg·K)),ω为含湿量(kg/m3),cp,l为液态水的比热(J/(kg·K)),T为墙体温度(K),t为时间(s),
Figure FDA0003726373600000014
为梯度算子;
步骤(4)导入三维物理模型;
步骤(5)输入参数,包括气象参数以及材料参数;
步骤(6)根据已知条件为各个区域选择合适的初始条件与边界条件,边界条件包括热边界条件与湿边界条件,用户可根据已知条件选择温度、热通量、热绝缘、湿度、水汽通量、防潮层边界条件;
步骤(7)划分网格,在网格划分方面COMSOL自带预定义的自由四面体网格单元划分方式,用户只需选择合理的网格尺寸对研究对象进行网格划分,网格尺寸越密,计算越精确,但所需计算时间越长,根据需求设定求解方式和时间步长;
步骤(8)计算结果;
步骤(9)选择几种不同的装配式建筑连接方式以及材料,根据各个月份的气象参数分别计算每个月的墙体温度与相对湿度参数。
2.如权利要求1所述的一种基于热湿耦合的装配式建筑热桥分析方法,其特征在于:
所述气象参数包括风速、温度、湿度、太阳辐射照度、天空有效温度、降雨以及盐雾浓度中的至少一种。
3.如权利要求1所述的一种基于热湿耦合的装配式建筑热桥分析方法,其特征在于:
所述网格尺寸包括至少对以下极端粗化、特别粗化、较粗化、粗化、正常、细化、较细化、特别细化、极端细化中的一种对研究对象进行网格划分。
4.如权利要求1所述的一种基于热湿耦合的装配式建筑热桥分析方法,其特征在于:
所述墙体材料可由MatWeb材料库导入。
5.如权利要求1所述的一种基于热湿耦合的装配式建筑热桥分析方法,其特征在于:
所述三维物理模型可由SolidWorks或CAD导入。
CN202110685229.3A 2021-06-21 2021-06-21 一种基于热湿耦合的装配式建筑热桥分析方法 Active CN113468638B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110685229.3A CN113468638B (zh) 2021-06-21 2021-06-21 一种基于热湿耦合的装配式建筑热桥分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110685229.3A CN113468638B (zh) 2021-06-21 2021-06-21 一种基于热湿耦合的装配式建筑热桥分析方法

Publications (2)

Publication Number Publication Date
CN113468638A CN113468638A (zh) 2021-10-01
CN113468638B true CN113468638B (zh) 2022-09-16

Family

ID=77868828

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110685229.3A Active CN113468638B (zh) 2021-06-21 2021-06-21 一种基于热湿耦合的装配式建筑热桥分析方法

Country Status (1)

Country Link
CN (1) CN113468638B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108009309A (zh) * 2017-10-17 2018-05-08 西安建筑科技大学 一种建筑材料导热系数计算方法
CN111812147A (zh) * 2020-06-24 2020-10-23 浙江大学 一种含热源土壤热湿耦合传递模化实验的设计方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10612834B2 (en) * 2016-07-26 2020-04-07 Whirlpool Corporation Method for manufacturing an insulated structure for a refrigerator
CN111305377A (zh) * 2020-03-06 2020-06-19 河北省建筑科学研究院有限公司 一种装配式钢结构超低能耗建筑处理方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108009309A (zh) * 2017-10-17 2018-05-08 西安建筑科技大学 一种建筑材料导热系数计算方法
CN111812147A (zh) * 2020-06-24 2020-10-23 浙江大学 一种含热源土壤热湿耦合传递模化实验的设计方法

Also Published As

Publication number Publication date
CN113468638A (zh) 2021-10-01

Similar Documents

Publication Publication Date Title
Wu et al. Phase change humidity control material and its impact on building energy consumption
Thiers et al. Thermal and environmental assessment of a passive building equipped with an earth-to-air heat exchanger in France
Kumar et al. Dynamic earth-contact building: A sustainable low-energy technology
CN103114661B (zh) 防冷凝外保温墙体结构
CN113468638B (zh) 一种基于热湿耦合的装配式建筑热桥分析方法
dos Santos et al. Numerical analysis of passive cooling using a porous sandy roof
Biplab et al. Comparative assessment of thermal comfort with insulation and phase change materials utilizations in building roofs and walls
Ashraf Ali et al. A Review on Innovative Nanomaterials for Enhancing Energy Performance of the Building Envelope
Mao et al. Experiment and simulation for the effect of window opening behavior on indoor air status points in wuhan area in summer
Ponikiewski et al. Moisture transport in cellular concrete walls with the connector for thermal insulation
Rahadian et al. The Evaluation of Thermal Comfort using a BIM-based Thermal Bridge Simulation (case Study: Itenas Mosque Building Bandung)
Gagliano et al. Computational fluid dynamic simulations of natural convection in ventilated facades
Lajimi et al. Numerical study of the thermal behavior of bi-zone buildings
AP BD et al. Experimental and simulation analysis of the performance of a low cost and space preserving ventilated façade for hot and humid climates
Vytchikov et al. Analyzing energy consumption while heating one-layer building envelopes in conditions of intermittent heating
Rode et al. Evaluation of moisture buffer effects by performing whole-building simulations
Siddique et al. Optimum insulation thickness for walls and roofs for reducing peak cooling loads in residential buildings in Lahore
Siegele et al. Modelling and simulation of façade integrated active components with Matlab/Simulink
Charai et al. Morocco’s housing challenge: Effect of current marketed building materials on the energy performance of residential buildings
Vacek et al. Case study of application of capillary active thermal insulation systems used as an interior insulation for historical buildings
Urzędowski et al. Numerical and Experimental Investigation of the Influence of External Wall Moisture on Building Temperature Distributions
Selvaraj et al. Numerical investigation of optimal wall materials for effective thermal performance in warm and hot climatic regions
Juras et al. Parametric study: Impact of selected factors on the energy demand of a family house
Lu et al. Analytical model for predicting whole building heat transfer
Ma et al. Finite element modelling of coupled heat and moisture transfer in typical earth-sheltered building envelope

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant