计算混凝土湿热耦合变形的解析-有限元结合解法.doc

计算混凝土湿热耦合变形的解析－有限元结合解法 陈德鹏(,2，钱春香，缪昌文1，刘加平1 (1. 江苏省建筑科学研究院有限公司 博特新材料研发中心, 南京 210008; . 安徽工业大学 建筑工程学院, 马鞍山 243002;. 东南大学 材料科学与工程学院, 南京 211189) 摘 要：。 关键词：中图分类号：TU528.01 1 引言 混凝土变形有其湿热变化本质混凝土的温度变形和干湿变形通常是同时发生的，混凝土内部的温度场和湿度场也是相互影响、相互作用的[1, 2]。国际著名混凝土专家Bazant教授很早就曾明确指出混凝土收缩是一种湿热变形（hygrothermal deformation）[3]。因此，温湿度变化及其耦合作用对混凝土内部温湿度分布及变化规律的影响，决定着混凝土的体积变形的发展，也极大影响着混凝土耐久性及使用寿命，研究混凝土体积变形必须考虑温湿度变化的耦合作用。 混凝土是一种典型的多孔介质材料，基于多孔介质传热传质理论对混凝土的传输、变形及耐久性能进行研究已经是目前混凝土研究的主要热点之一[4-8]。以Bazant[3]、Gawin[9]、Meschke[10]等为代表的众多研究者对混凝土应力及变形的研究也都建立在多孔介质传热传质理论的基础上；刘光廷等对混凝土温度场和湿度场的计算模拟研究也是基于混凝土的多孔介质特性而开展的[11]；李荣涛等[12]基于多孔介质相关原理对高温下混凝土中化学－热－湿－力学耦合过程的数值模拟及破坏进行了研究；陈松等[7]考虑水分对混凝土徐变重要影响，基于多孔介质理论对混凝土的徐变力学行为进行了研究。 实现对混凝土湿热耦合作用进行数值模拟分析的主要有种可能：①基于多孔介质热、质迁移的唯象热力学方法；②以Luikov传热传质耦合方程为基础的数值分析方法；③有限元分析方法。因受试验条件和许多材料参数的影响，唯象学分析方法难以在实际中应用于混凝土[13]；以Luikov传热传质方程为基础的解析法在求解温湿度耦合过程中，要么由于存在很难计算的复杂特征值而得不到结果，要么在控制方程和边界条件同时耦合时无能无力，要么准确性欠佳难以反映出多孔介质内的温湿度分布，用传递函数对耦合传输进行求解的方法则可行[14]，但对于温、湿度同时变化的问题，解析过程将非常复杂；对于第三种，目前用主流的商业有限元软件无法直接计算析混凝土中的湿度应力及变形。而对于号称能够实现多物理场耦合的、基于偏微分方程建立的有限元软件FEMLAB和FEPG（飞箭），能否用来对混凝土中湿热传输耦合作用进行模拟分析、能否简单易行的直接实现湿热耦合变形分析，还有待于相关实践的验证。本文采用数值解析法与有限元分析相结合的方法对混凝土湿热耦合变形进行计算。 解析－有限元结合解法 2.1 解析－有限元结合解法的方法体系 用以计算混凝土湿热耦合变形的解析－有限元结合解法，主要是基于较为成熟的多介质湿热传输理论以及有限元计算混凝土变形时难以处理湿度引起的应力及变形而提出的。在解析－有限元结合解法体系中，主要包括温湿度场的解析求解和湿热耦合变形的有限元分析两部分。此外，由于有限元软件难以处理湿度分布引起的应力与变形，进行有限元分析前采用理论公式将湿度分布进行应力转化。 因此，混凝土变形模拟的整个过程主要分三步：首先以Luikov传热传质方程为基础解析法耦合作用下的温、湿度分布；然后根据混凝土干湿变形机理，通过相关方程计算出湿度分布引起的应力分布；最后根据温度分布及湿度应力分布，通过有限单元法计算混凝土在温、湿度耦合作用下的变形及规律。 图1 混凝土湿热耦合变形数值模拟方法流程图 Fig.1 Schematic of the numerical simulation for hygro-thermal deformation of concrete 2.2 混凝土温湿度分布的解析法求解 根据理论分析[14]，考虑到混凝土内部可能存在热源（热汇）或湿源（湿汇）、描述混凝土中湿扩散时不可忽略的Knudsen扩散作用[15]以及湿热物参数间的相互影响，可建立混凝土湿热耦合传输的微分方程： （） （） 式中， M；为湿源或湿汇，因水泥水化作用降低湿含量或其它非湿迁移过程引起的湿含量变化变化；为考虑Knudsen扩散影响的混凝土湿扩散系数，；为Knudsen扩散影响系数；Dm 为多孔材料的密度；λ为多孔材料的导热系数；T 为温度；r 为相变因子；hlv为蒸发潜热；δ为热梯度系数。 在根据湿热耦合传输方程进行混凝土内部温湿度分布计算时，需根据实际情况确定分析对象的边界条件、初始条件以及Q和W的存在与否。为求解方便，可首先对方称进行简化，然后对基本方程、边界条件初