含孔洞和裂隙混合缺陷的干燥水泥砂浆导热系数相互作用直推预测模型

李云波; 李宗利; 姚希望; 肖帅鹏; 刘士达; 童涛涛

doi:10.13801/j.cnki.fhclxb.20210328.002

含孔洞和裂隙混合缺陷的干燥水泥砂浆导热系数相互作用直推预测模型

李云波^1,,
李宗利^{1, 2, ,},
姚希望^1,,
肖帅鹏^1,,
刘士达^1,,
童涛涛^1,

1.
西北农林科技大学　水利与建筑工程学院，杨凌 712100
2.
西北农林科技大学　旱区农业水土工程教育部重点实验室，杨凌 712100

基金项目: 国家重点研发计划(2017YFC0405101)

详细信息

通讯作者:
李宗利，博士，教授，博士生导师，研究方向为混凝土力学性能及耐久性　E-mail：bene@nwsuaf.edu.cn

中图分类号: TU52；TV43
计量
- 文章访问数: 690
- HTML全文浏览量: 376
- PDF下载量: 39
出版历程
- 收稿日期: 2021-01-26
- 修回日期: 2021-03-09
- 录用日期: 2021-03-22
- 网络出版日期: 2021-03-28
- 刊出日期: 2022-01-14

Interaction direct deduction prediction model of thermal conductivity of dry cement mortar with mixed defects of cavities and cracks

1.
College of Water Resources and Architectural Engineering, Northwest A & F University, Yangling 712100, China
2.
Key Laboratory of Agricultural Soil and Water Engineering in Arid and Semiarid Areas of Ministry of Education, Northwest A&F University, Yangling 712100, China

摘要

摘要: 水泥砂浆在浇筑、养护和受荷载后易出现缺陷，会显著改变其热传导性能。将水泥砂浆中缺陷分为孔洞和裂隙两种基本类型，基于细观力学理论中的相互作用直推(IDD)估计法，建立了孔洞和裂隙单独存在或同时存在时干燥水泥砂浆等效导热系数预测模型；基于数值模拟结果反推法，提出裂隙随机分布影响函数，改进了IDD模型预测精度；与数值模拟结果进行对比验证其合理性。结果表明，推导出的预测模型充分体现了孔洞、裂隙不同含量及特征对水泥砂浆等效导热系数的影响，预测精度高；提出的裂隙随机分布影响函数具有较高精度和良好适用性；在缺陷率一定的情况下，孔洞对水泥砂浆等效导热系数的影响比裂隙更显著。该预测模型物理基础明确，形式简单，便于工程应用。
- 细观力学 /
- 缺陷 /
- 导热系数 /
- 预测模型 /
- 相互作用直推(IDD)理论 /
- 水泥砂浆
Abstract: Cement mortar is prone to product defects during of pouring, curing and loading, which will significantly change its heat conducting properties. The defects in cement mortar were divided into two basic types of cavities and cracks according to its external characteristics and effects. Based on the interaction direct deduction (IDD) estimation method of mesomechanics theory, the prediction models of equivalent thermal conductivity of dry cement mortar with both cavities and cracks alone or in combination were established. Based on the inverse method of numerical simulation results, the random distribution influence function of the cracks was proposed, which improved the prediction accuracy of the proposed model. The IDD calculation results were compared with the numerical simulation results to verify its rationality and accuracy. The results show that the prediction model fully reflects the influence of different contents and characteristics of cavities and cracks on the equivalent thermal conductivity of cement mortar. The model has high prediction accuracy. The proposed random distribution influence function of cracks has high accuracy and good applicability. In the case of a certain defect rate, the effect of cavities on the equivalent thermal conductivity of cement mortar is more significant than that of cracks. The prediction model has a clear physical basis and a simple form, which is convenient for engineering applications.
- mesomechanics /
- defect /
- thermal conductivity /
- prediction model /
- IDD theory /
- cement mortar

HTML全文

水泥基材料的等效导热系数是影响其温度场的重要参数之一，对研究材料的热传导性能和温度效应问题具有重要意义和价值。在拌和、浇筑和振捣的过程中由于工艺和方法的限制，水泥砂浆内部会产生一定气泡，凝固后则形成孔洞；加上水泥砂浆在硬化干燥的过程中由于干缩效应及水化热产生的温度应力及服役期间受到各种荷载，内部会产生一定量微小裂隙，这些孔洞和微小裂隙均会对水泥砂浆的热传导性能产生一定的影响。在细观层次，孔洞和微小裂隙主要存在于水泥基材料的砂浆相。因此，研究不同形式、不同特征的缺陷对砂浆等效导热系数的影响规律是预测水泥砂浆等效导热系数的关键。随着现代材料微观检测手段日趋完善，工业CT及图像处理技术已广泛应用于多孔材料的微观结构检测，基于复合材料细观力学理论建立混凝土或砂浆的导热系数预测模型显得极为必要。

含孔洞和随机裂隙的砂浆可以认为是均质复合材料，其热传导性能受各组分的导热系数影响。在复合材料方面，陈则韶等^[1]基于最小热阻理论和比等效导热系数相等法则，得出考虑复合材料组分的形状及取向对导热系数影响的具体式。部分学者以Maxwell模型为基础^[2-6]，在复合材料等效导热系数预测方面进行了研究。还有一些学者从细观力学角度，应用细观力学模型对复合材料有效导热系数进行了预测研究^[7-10]。孙海浩等^[11]采用随机生成结构法和链表数据结构得到了纤维体等效导热系数理论预测方法。张枫等^[12]通过实验考察了朱伯芳模型等几个模型在混凝土导热系数上应用的合理性。唐世斌等^[13]和孙颖颖等^[14]研究了混凝土细观结构对砂浆和混凝土有效导热系数的影响。李守巨等^[15]采用有限元方法研究了复合材料的等效导热系数。目前对水泥基材料导热性能研究中，也有很多学者考虑了孔隙率的影响^[3,4,6,15]，但并未将孔洞与裂隙分开，分别考虑含量及特征影响。

多孔介质的热传导特性是线性物质的一种传导性质，Shafiro等^[16]给出了线性物理问题中含椭球形夹杂基体的热传导问题稀疏解。Zheng等^[17]在建立有效自洽法的基础上进一步提出了相互作用直推(Interaction direct deduction，IDD)估计。Zhou等^[18-19]采用IDD法计算裂缝作用下的固体渗透率，并研究了裂纹特征对多孔材料渗透特性的影响。IDD估计与其他现有的细观力学估计相比，结构更简单明确，普遍适用于各种夹杂类型的多相复合材料，可考虑夹杂物的分布、形状、取向和夹杂密度等影响^[17]。在物理现象的本质层面，热传导过程与渗透过程、电导问题等均属于线性物质传导问题，数学方程的表达形式相同^[20]。

本文基于IDD估计框架分别推导出含有孔洞与裂隙缺陷夹杂的干燥水泥砂浆等效导热系数预测模型。提出裂隙随机分布影响函数，来考虑裂隙交叉、连通和近场相互作用等因素造成的有效密度变化。建立便于工程应用的考虑不同形式缺陷夹杂的干燥水泥砂浆导热系数预测模型，可为损伤对水泥基材料等效导热系数影响提供理论支撑。

1. 基于IDD理论的水泥砂浆等效导热系数预测模型的建立

图1为含有孔洞和裂隙的水泥砂浆混合缺陷模型。将孔洞和裂隙认为是椭圆形夹杂的特例，以最基本的椭圆夹杂推导其等效导热系数预测模型。

图 1 含孔洞和裂隙的水泥砂浆混合缺陷模型

Figure 1. Mixed defect model with cavities and cracks of cement mortar

下载: 全尺寸图片幻灯片

1.1 稀疏解

稀疏解作为细观力学方法中的一种基础解法，是应用IDD理论框架研究细观夹杂对砂浆等效导热系数影响的理论基础。将基体相的导热系数记为λ₀，第i个夹杂相的导热系数记为λ_i，复合材料整体导热系数记为λ，整体复合材料有效介质的导热系数的增量记为K，夹杂相对基体相的导热系数的增量记为K_i，则它们具有以下关系：

${{\boldsymbol{K}}} = {{\boldsymbol{\lambda }}} - {{\boldsymbol{\lambda }}_{\rm{0}}}$

(1)

${{\boldsymbol{K}}_{{i}}} = {{\boldsymbol{\lambda }}_{{i}}} - {{\boldsymbol{\lambda }}_{\rm{0}}}$

(2)

线性物理问题的椭球型夹杂稀疏解H^d表示为^[9]

${{\boldsymbol{H}}^{\rm{d}}} = \sum\limits_{{i}} {{\boldsymbol{H}}_{{i}}^{\rm{d}}}$

(3)

${\boldsymbol{H}}_{{i}}^{\rm{d}}{\rm{ = }}{c_{{i}}}{{\boldsymbol{K}}_{{i}}}{\boldsymbol{A}}$

(4)

式中： ${{{ {\boldsymbol{H}}}}}_{{i}}^{\rm{d}}$ 为第i个夹杂所对应的导热系数稀疏解；c_i是第i个夹杂的体积分数；对于各向同性传热的夹杂问题，A等价于弹性问题中的局部化张量^[18]，可按下式计算：

${\boldsymbol{A}} = {\left[ {{\boldsymbol{I}} + {{\boldsymbol{K}}_{{i}}}{\boldsymbol{\lambda }}_0^{ - 1}{\boldsymbol{J}}} \right]^{ - 1}}$

(5)

${{ \boldsymbol I}}{\rm{ = }}{{\boldsymbol e}_1}{{\boldsymbol e}_1}{\rm{ + }}{{\boldsymbol e}_2}{{\boldsymbol e}_2}{\rm{ + }}{{\boldsymbol e}_3}{{\boldsymbol e}_3}$

(6)

式中：I为单位张量； ${{\boldsymbol e}_1}{{\boldsymbol e}_1}、{{\boldsymbol e}_2}{{\boldsymbol e}_2}和{{\boldsymbol e}_3}{{\boldsymbol e}_3}$ 分别代表空间中三个方向的单位矢量；J是对角形式的几何特征张量，对应于基体氛围的取向及大小。其分量J₁、J₂和J₃可用标准化椭圆积分进行计算：

${{ J}_{ 1}} = \frac{{{a_1}{a_2}{a_3}}}{2}\int_0^\infty {\frac{{{\rm{d}}u}}{{\left( {a_1^2 + u} \right)\sqrt {\left( {a_1^2 + u} \right)\left( {a_2^2 + u} \right)a_3^2 + u} }}}$

(7)

式中：a₁、a₂和a₃分别为椭圆三个轴的半轴长(图2)。J₂和J₃可通过指标轮换进行计算。

图 2 椭球夹杂示意图

Figure 2. Schematic diagram of ellipsoid inclusion

下载: 全尺寸图片幻灯片

将式(5)代入式(4)，复合材料导热系数增量的稀疏解可表示为

${{\boldsymbol{H}}^{\rm{d}}} = \sum\limits_{{i}} {{c_{{i}}}{{\left( {{{\boldsymbol{K}}_{{i}}}^{ - 1} + {{\boldsymbol{\lambda }}_0}^{ - 1}{\boldsymbol{J}}} \right)}^{ - 1}}}$

(8)

1.2 IDD估计框架

IDD理论框架是针对复合材料的弹性性质提出来的。它的结构简洁、统一，可以将各个量的物理意义清晰地表示出来，适用于含有各种几何形状夹杂、且任意分布的多相复合材料有效性质的计算^[21]。考虑复合材料夹杂间相互作用的有效导热系数IDD理论解H^idd表示如下：

${{\boldsymbol{H}}^{{\rm{idd}}}} = {\left( {{\boldsymbol{I}} - \sum\limits_{{i}} {{\boldsymbol{H}}_{{i}}^{\rm{d}}{\boldsymbol{\varOmega }}_{{\rm{D}}i}^0} } \right)^{ - 1}}{{\boldsymbol{H}}^{\rm{d}}}$

(9)

${\boldsymbol{\varOmega }}_{{\rm{D}}i}^0{\rm{ = }}{\boldsymbol{\lambda }}_0^{ - 1}{\boldsymbol{J}}$

(10)

式中： ${\boldsymbol{\varOmega }}_{{\rm{D}}i}^0$ 为基体氛围所对应的特征刚度张量。

1.3 不同类型缺陷夹杂的IDD解

1.3.1 裂隙夹杂

本文认为缺陷内部充满干燥空气，无孔砂浆基体导热系数为1.65 W/(m·K)^[6]，干燥空气的导热系数为0.0259 W/(m·K)^[22]。相对于砂浆基体的导热系数，空气的导热系数仅为其1.6%。为了简化处理，充填干燥空气的裂隙夹杂可认为其绝热(λ_i = 0)。整理式(8)可得裂隙夹杂时砂浆导热系数的稀疏解 ${\boldsymbol{ H}}_{\rm{c}}^{\rm{d}}$ 为

${\boldsymbol{H}}_{\rm{c}}^{\rm{d}} = \sum\limits_{{i}} {{c_{{i}}}{{\boldsymbol{\lambda }}_0}{{\left[ {{\boldsymbol{I}} - {\boldsymbol{J}}} \right]}^{ - 1}}}$

(11)

根据文献[18]，裂隙夹杂可认为由椭球退化成椭圆，再进一步退化成为裂隙。如图2所示，为一个椭球夹杂在立方体基体中，对于椭圆夹杂可认为椭球的a₁趋于无穷大。椭圆夹杂的几何张量J的各个分量分别为^[18]

${J_1} = 0$

(12)

${J_2} = \frac{{{a_3}}}{{{a_2} + {a_3}}}$

(13)

${J_3} = \frac{{{a_2}}}{{{a_2} + {a_3}}}$

(14)

将式(12)~(14)代入式(11)整理可得

${\boldsymbol{H}}_{\rm{c}}^{\rm{d}} = \sum\limits_{{i}} {\frac{{{\rm{{\text{π}} }}{a_2}{a_3}}}{A}} {{\boldsymbol{\lambda }}_0}\left[ \begin{gathered} {{\boldsymbol{e}}_1}{{\boldsymbol{e}}_1} + \frac{{{a_2} + {a_3}}}{{{a_2}}}{{\boldsymbol{e}}_2}{{\boldsymbol{e}}_2} + \frac{{{a_2} + {a_3}}}{{{a_3}}}{{\boldsymbol{e}}_3}{{\boldsymbol{e}}_3} \\ \end{gathered} \right]$

(15)

对于裂隙夹杂，可认为a₃→0，将其代入式(15)，则裂隙夹杂的稀疏解为

${\boldsymbol{H}}_{\rm{c}}^{\rm{d}} = \sum\limits_{{i}} {\frac{{{\rm{{\text{π}} }}a_{{i}}^2}}{A}} {{\boldsymbol{\lambda }}_0}{{\boldsymbol{e}}_3}{{\boldsymbol{e}}_3}$

(16)

式中：a_i为裂隙的半长。裂隙密度^[18] $\rho = \displaystyle\sum a_{{i}}^2{\rm{/}}A$ ，则上式可进一步整理为

${\boldsymbol{H}}_{\rm{c}}^{\rm{d}} = {\rm{{\text{π}} }}\rho {{\boldsymbol{\lambda }}_0}{{\boldsymbol{e}}_3}{{\boldsymbol{e}}_3}$

(17)

对于二维的裂隙夹杂，其特征刚度矩阵 ${\boldsymbol{\varOmega }}_{{\rm{D}}i}^0$ 可按下式计算^[18]：

${\boldsymbol{\varOmega }}_{{\rm{D}}i}^0{\rm{ = }}{\boldsymbol{\lambda }}_0^{ - 1}\left[ {\gamma {{\boldsymbol{e}}_2}{{\boldsymbol{e}}_2} - \left( {1 - \gamma } \right){{\boldsymbol{e}}_3}{{\boldsymbol{e}}_3}} \right]$

(18)

$\gamma {\rm{ = }}\dfrac{{\sqrt {1{\rm{ + }}\dfrac{4}{{{\text{π}} \rho }}} {\rm{ - }}1}}{{\sqrt {1{\rm{ + }}\dfrac{4}{{{\text{π}} \rho }}} {\rm{ + }}1}}$

(19)

式中，γ为长宽比。

将式(17)~(18)代入式(9)便可得到平行裂隙垂直于温度梯度方向夹杂时的IDD解增量 ${\boldsymbol{H}}_{\rm{c}}^{{\rm{idd}}}$ 为

${\boldsymbol{H}}_{\rm{c}}^{{\rm{idd}}} = - \frac{{{\rm{{\text{π}} }}\rho {\lambda _0}}}{{1 - {\rm{{\text{π}} }}\rho \left( {1 - \gamma } \right)}}{{\boldsymbol{e}}_3}{{\boldsymbol{e}}_3}$

(20)

则平行裂隙垂直于温度梯度方向夹杂时，等效导热系数IDD解 ${{\boldsymbol{\lambda}} _{{\rm{cp}}}}$ 为

${{\boldsymbol{\lambda}} _{{\rm{cp}}}} = \left[ {1 - \frac{{{\rm{{\text{π}} }}\rho }}{{1 - {\rm{{\text{π}} }}\rho \left( {1 - \gamma } \right)}}} \right]{{\boldsymbol{\lambda}} _0}$

(21)

对于随机裂隙夹杂时，等效导热系数IDD解 ${{\boldsymbol{\lambda}} _{{\rm{cr}}}}$ 为

${{\boldsymbol{\lambda}} _{{\rm{cr}}}} = \left[ {1 - \frac{{{\rm{{\text{π}} }}\rho }}{{2 - {\rm{{\text{π}} }}\rho \left( {1 - \gamma } \right)}}} \right]{{\boldsymbol{\lambda}} _0}$

(22)

当两个或两个以上的微裂隙随机分布时，由于裂隙之间的交叉，连通和近场相互作用，它们会放大或缩小传导效应，这种现象称为聚类的放大或缩小效应^[23]。对于本文研究的传热问题，式(21)和式(22)预测随机分布的裂隙对热传导产生的影响时，由于未能考虑聚类的放大或缩小效应，会使得导热系数整体计算结果偏小。文献[18]在研究裂隙的随机分布特征对渗透性的影响时引入裂隙连通系数 $\phi$ 和缩放系数α，通过计算有效裂隙密度ρ′，综合反映裂隙的随机分布对其结果的影响：

$\rho '{\rm{ = }}\left[ {1{\rm{ + }}\left( {{\alpha ^2}{\rm{ + }}2\alpha } \right)\phi } \right]\rho$

(23)

文献[18]根据精心设计的算法，对随机裂隙模型中的裂隙的连通数量和相交数量进行精确控制，来对其进行定量分析。但裂隙连通系数 $\phi$ 和缩放系数α计算较为繁琐，且在实际工程中很难实现对裂隙连通数量和交叉数量进行精确统计分析，因此通过后续数值模拟结果反推，将其统一归于裂隙随机分布影响函数ζ(ρ)，并综合验证了反推结果合理性和普遍适用性。对于裂隙随机分布的情况，材料近似被视为各向同性；裂隙随机平行分布时，由于裂隙的方向性使材料变为正交各向异性，其随机分布影响函数分别为

$\zeta \left( \rho \right){\rm{ = }}\left\{ {\begin{array}{*{20}{l}} { - 3.429{\rho ^2} - 2.069\rho + 0.704}&{{\rm{RPC}}} \\ {0.606}&{{\rm{RC}}} \end{array}} \right.$

(24)

式中：RPC为随机平行分布裂隙；RC为随机分布裂隙。

根据裂隙的随机分布影响函数计算有效密度：

$\rho '{\rm{ = }}\xi \left( \rho \right)\rho$

(25)

随机平行裂隙垂直于温度梯度方向时，等效导热系数修正IDD解 ${{\boldsymbol{\lambda}} '_{{\rm{cp}}}}$ 为

${{\boldsymbol{\lambda}} '_{{\rm{cp}}}} = \left[ {1 - \frac{{{\rm{{\text{π}} }}\rho '}}{{1 - {\rm{{\text{π}} }}\rho '\left( {1 - \gamma '} \right)}}} \right]{{\boldsymbol{\lambda}} _0}$

(26)

式中，γ′为有效密度对应的长宽比(式(19))。

根据最小热阻原理和裂隙密度定义，随机平行裂隙方向与温度梯度方向夹角为θ时的等效导热系数修正IDD解 ${{\boldsymbol{\lambda}} '_{{\rm{cp\theta }}}}$ 为

${{\boldsymbol{\lambda}} '_{{\rm{cp\theta }}}} = \left[ {1 - \frac{{{\rm{{\text{π}} }}\rho '{{\sin }^2}\theta }}{{1 - {\rm{{\text{π}} }}\rho '\left( {1 - \gamma '} \right){{\sin }^2}\theta }}} \right]{{\boldsymbol{\lambda}} _0}$

(27)

随机分布裂隙夹杂时，等效导热系数修正IDD解 ${{\boldsymbol{\lambda}} '_{{\rm{cr}}}}$ 为

${{\boldsymbol{\lambda}} '_{{\rm{cr}}}} = \left[ {1 - \frac{{{\rm{{\text{π}} }}\rho '}}{{2 - {\rm{{\text{π}} }}\rho '\left( {1 - \gamma '} \right)}}} \right]{{\boldsymbol{\lambda}} _0}$

(28)

1.3.2 孔洞夹杂

与裂隙夹杂同理，孔洞夹杂项为热绝缘体(λ_i= 0)。对于圆形孔洞夹杂可认为是椭圆夹杂中的a₂ = a₃。根据式(12)~(14)，圆孔夹杂的几何张量J的各个分量J₁ = 0，J₂ = 0.5和J₃ = 0.5。将其代入式(11)，可得到圆孔夹杂的稀疏解 ${\boldsymbol{H}}_{\rm{h}}^{\rm{d}}$ ：

${\boldsymbol{H}}_{\rm{h}}^{\rm{d}} = 2{{\boldsymbol{\lambda}} _0}c$

(29)

式中，c为孔洞率，按下式计算：

$c{\rm{ = }}\sum\limits_{{i}} {\frac{{{\rm{{\text{π}} }}r_{{i}}^2}}{A}}$

(30)

式中，r_i为第i个圆孔的半径。

将式(29)代入式(9)可得到圆孔夹杂时的IDD解增量 ${\boldsymbol{H}}_{\rm{h}}^{{\rm{idd}}}$ 为

${\boldsymbol{H}}_{\rm{h}}^{{\rm{idd}}} = - \frac{{2c{{\boldsymbol{\lambda}} _0}}}{{1 + c}}$

(31)

则圆孔夹杂的复合材料等效导热系数IDD解 ${{\boldsymbol{\lambda}} _{\rm{h}}}$ 为

${{\boldsymbol{{\boldsymbol{\lambda}}}} _{\rm{h}}} = \left( {1 - \frac{{2c}}{{1 + c}}} \right){{\boldsymbol{\lambda}} _0}$

(32)

1.3.3 多种形式缺陷夹杂

在实际工程中，砂浆内部缺陷往往不是单一存在，而是裂隙和孔洞并存。可通过叠加法，得到等效导热系数预测公式。

随机平行裂隙方向与温度梯度方向夹角为θ时修正IDD解 ${{\boldsymbol{\lambda}} '_{{\rm{p\theta }}}}$ 为

${{\boldsymbol{\lambda}} '_{{\rm{p\theta }}}} = \left( {1 - \frac{{2c}}{{1 + c}} - \frac{{{\rm{{\text{π}} }}\rho '{{\sin }^2}\theta }}{{1 - {\rm{{\text{π}} }}\rho '\left( {1 - \gamma '} \right){{\sin }^2}\theta }}} \right){{\boldsymbol{\lambda}} _0}$

(33)

裂缝随机分布时修正IDD解 ${{\boldsymbol{\lambda}} '_{\rm{r}}}$ 为

${{\boldsymbol{\lambda}} '_{\rm{r}}} = \left( {1 - \frac{{2c}}{{1 + c}} - \frac{{{\rm{{\text{π}} }}\rho '}}{{2 - {\rm{{\text{π}} }}\rho '\left( {1 - \gamma '} \right)}}} \right){{\boldsymbol{\lambda}} _0}$

(34)

2. 基于数值模拟对干燥水泥砂浆导热系数IDD预测模型的验证

基于数值模拟验证含有缺陷夹杂的水泥砂浆等效导热系数IDD解的合理性。无孔砂浆基体导热系数为1.65 W/(m·K)^[6]，干燥空气的导热系数为0.0259 W/(m·K)^[22]。

2.1 温度传导方程

稳态热传导满足热力学第一定律^{[22, 24]}，基于能量守恒定律和傅里叶热传导定律，直角坐标系下砂浆的热传导微分方程为

$\nabla \left( {{\boldsymbol{\lambda}} \nabla T} \right) = 0$

(35)

式中： $\nabla$ 为拉普拉斯算子；T为温度。

2.2 初始值和边界条件

采用150×150 mm的含有缺陷的二维几何模型。温度初始值为25℃；上下两侧边界定义为热绝缘；左右两侧为恒温边界，左侧温度为60℃，右侧温度为25℃。网格划分如图3所示，在裂隙和孔洞区域进行加密处理。

图 3 干燥水泥砂浆网格划分图

Figure 3. Meshing diagram of dry cement mortar

下载: 全尺寸图片幻灯片

2.3 干燥水泥砂浆随机缺陷模型建立

为模拟求解缺陷对基体导热系数的影响，应用MATLAB软件编制随机缺陷生成程序。图4为随机缺陷程序流程图，利用随机函数生成孔洞和裂隙等随机缺陷；距离公式判断是否与边界及其他缺陷有交叉；循环语句循环生成孔洞、裂隙等缺陷，以达到所需含量。通过控制缺陷形状、尺寸、含量、角度等特征参数，生成含不同形式、不同尺寸、不同含量、不同角度缺陷的水泥砂浆模型。

图 4 随机缺陷程序流程图

Figure 4. Flow chart of random defect program

下载: 全尺寸图片幻灯片

2.4 数值模拟求解干燥水泥砂浆导热系数

将含有缺陷夹杂的几何模型导入到有限元软件COMSOL中，应用固体传热模块进行细观热传导数值模拟。由于裂隙的长度方向远大于厚度方向，在进行有限元计算时，采用薄层对基体中含有的裂隙进行特殊处理。

3. 不同缺陷特征对水泥砂浆等效导热系数的影响

模型将砂浆中缺陷分为孔洞与裂隙两类，在二维平面中，缺陷分别用圆形孔洞和线状裂隙表示。研究不同缺陷特征情况下，干燥砂浆导热系数的变化规律，并将基于细观力学理论IDD法计算的结果与应用有限元软件求解的导热系数进行对比分析。

3.1 裂隙特征对水泥砂浆等效导热系数的影响

由于裂隙的形状特性，其长度、数量、方向不同，会导致裂隙对砂浆等效导热系数产生影响不同。如图5和图6所示，建立随机平行分布和随机分布两种裂隙缺陷模型，从以下两个方面研究裂隙特征对砂浆导热系数的影响规律。

图 5 平行裂隙模型(ρ=0.05、θ=60°)

Figure 5. Parallel crack model (ρ=0.05, θ=60°)

下载: 全尺寸图片幻灯片

图 6 随机裂隙模型(ρ=0.05)

Figure 6. Random crack model (ρ=0.05)

下载: 全尺寸图片幻灯片

3.1.1 平行裂隙的方向影响

裂隙方向作为裂隙的主要特征之一，对导热系数的影响规律值得探究。宏观导热系数作为材料的属性，研究裂隙水平和垂直两个具有代表性的方向对砂浆导热系数的影响具有重要意义。建立二维定向裂隙模型，固定裂隙数量为200，取裂隙密度ρ分别为0.05、0.10、0.15、0.20和0.25，生成水平裂隙和垂直裂隙夹杂的砂浆模型。并选择其中裂隙密度ρ为0.05、0.15和0.25的模型，使其裂隙方向从水平到垂直之间以10°为步长逐渐增长。进行数值模拟研究，为了减小离散性对数值模拟结果的影响，每个方向的模型生成三组，求其平均值作为最终结果。

基于最小热阻理论，热流在砂浆基体内流动时，会沿着阻力最小的方向流动。如图7和图8所示，裂隙方向平行于温度梯度方向时，对热传导几乎不产生阻碍作用，砂浆导热系数等于基体导热系数1.65 W/(m·K)；当裂隙方向垂直于温度梯度方向时对导热系数影响最大，且随裂隙密度增大等效导热系数减小。基于此通过数值模拟结果反推法对平行裂隙夹杂的砂浆等效导热系数预测模型进行了修正，得到式(24)的裂隙随机分布影响函数。

图 7 裂隙密度对水泥砂浆正交方向导热系数的影响

Figure 7. Influence of crack density on the thermal conductivity of cement mortar in the orthogonal direction

CIDD—Correct IDD solution; IDD—Interaction direct deduction solution; NS—Numerical simulation solution

下载: 全尺寸图片幻灯片

图 8 裂隙方向对水泥砂浆导热系数的影响

Figure 8. Influence of crack direction on thermal conductivity of cement mortar

下载: 全尺寸图片幻灯片

图8为数值模拟计算结果与修正IDD理论解对比。结果显示：不同裂隙密度下，不同角度时数值模拟解与修正IDD解的变化趋势一致且误差均很小。数值模拟解整体略大于修正IDD解，且随着裂隙密度的增加，误差逐渐增大，这主要是由于数值模拟是基于有限单元模型进行求解，而修正IDD解的提出是基于无限大基体。本文建立的考虑夹杂间相互影响的修正IDD公式可以很好地预测不同方向裂隙对砂浆等效导热系数的影响。

3.1.2 随机裂隙长度与数量影响

根据裂隙密度定义，裂隙长度的增加会引起裂隙密度呈幂指数增加，对导热系数的影响也会随之增加。同时，裂隙数量的变化，也会导致裂隙缺陷对导热系数的影响发生变化。为了研究长度与数量对砂浆导热系数的影响规律，建立含随机角度裂隙缺陷的砂浆模型，裂隙长度分别为5 mm、7.5 mm、10 mm时，裂隙数量从20到200递增，步长为20，进行数值模拟求解，并与IDD计算结果进行对比分析。

图9为数值模拟计算结果与IDD及修正IDD预测结果对比。结果显示：数值模拟值整体大于随机分布IDD解，主要是由于IDD理论模型在计算的过程中暂未考虑裂隙交叉、连通和近场相互作用等因素对导热系数阻碍作用的缩放效应。为了提高IDD理论模型的预测精度，基于数值模拟对其进行了修正，弥补了这方面的不足。

图 9 裂隙长度与数量对水泥砂浆导热系数的影响

Figure 9. Influence of crack length and quantity on thermal conductivity of cement mortar

下载: 全尺寸图片幻灯片

由图9还可以看出，在裂隙长度一定时，砂浆导热系数随着其数量的增加呈线性递减，且裂隙长度越长，导热系数随数量增加而减小的速率越快。在裂隙数量一定时，导热系数随裂隙长度增加而减小，数量越多，导热系数随长度变化越显著。而长度和数量的变化，本质上会引起裂隙密度的变化，如图10所示，即为不同长度和不同数量所对应的裂隙密度下的导热系数变化。

图 10 随机裂隙密度对水泥砂浆导热系数的影响

Figure 10. Effect of random crack density on thermal conductivity of cement mortar

L—Crack length

下载: 全尺寸图片幻灯片

图9和图10结果表明，应用反推的裂隙随机分布影响系数对于预测含不同长度、不同数量裂隙的砂浆导热系数具有普遍适用性。

3.2 孔洞特征对水泥砂浆等效导热系数的影响

孔洞的特征参数包括孔径与孔洞数量。由于浇筑、振捣等因素影响，不同的砂浆中会产生不同孔径、不同数量的孔洞缺陷。为了简化计算并定量分析其影响，将孔洞缺陷等效化处理，如图11所示，所得结果用以反映影响规律。建立等效孔径分别为2 mm、4 mm、6 mm时，砂浆的孔洞数量在20~200之间变化，步长为20的模型。

图 11 随机孔洞模型(c=0.05)

Figure 11. Random cavity model (c=0.05)

下载: 全尺寸图片幻灯片

图12为IDD理论预测与数值模拟计算的砂浆导热系数结果。由于干燥的水泥砂浆，认为其孔洞内部充满空气，而空气的导热系数远远小于砂浆基体的导热系数，孔洞的存在会对砂浆热传导产生明显阻碍作用，干燥砂浆的导热系数随着孔径与孔洞数量的增加呈现明显的减小趋势。

图 12 不同孔洞特征对水泥砂浆导热系数的影响

Figure 12. Influence of different cavity characteristics on thermal conductivity of cement mortar

d—Cavity diameter

下载: 全尺寸图片幻灯片

由图12还可以看出，数值模拟结果与理论计算结果变化趋势一致，且误差很小，故IDD理论在求解含孔洞缺陷的砂浆导热系数时具有很高的可靠度。计算值存在较小误差，且理论计算值整体略大于数值模拟值，这是由于理论推导过程中，为了简化计算，将充满干燥空气的孔洞视为热绝缘体，因此与数值模拟结果产生了微小误差。随着夹杂含量的增加，将孔洞视为绝缘夹杂带来的微小误差会累积，导致误差逐渐增大，但整体误差依然较小。在孔径一定的情况下，砂浆的导热系数随孔洞数量的增加呈减小趋势，且孔径越大，导热系数随孔洞数量变化越显著。究其原因，孔径和孔洞数量的变化导致孔洞率变化，引起的导热系数变化，如图13所示，即为图12中不同孔径和孔洞数量对应的孔洞率下砂浆的导热系数变化。

图 13 孔洞率对水泥砂浆导热系数的影响

Figure 13. Influence of porosity on thermal conductivity of cement mortar

下载: 全尺寸图片幻灯片

3.3 多种缺陷夹杂的水泥砂浆导热系数模型

在实际砂浆中孔洞和裂隙是同时存在的，应用MATLAB编制程序，建立孔洞与裂隙含量占比例不同、总含量分别为0.05、0.10和0.15的混合缺陷砂浆模型，如图14所示，其等效孔径为2~6 mm，裂隙长度为2~10 mm。

图 14 含孔洞和裂隙缺陷模型(ρ=0.04、c=0.06)

Figure 14. Defect model with cavities and cracks (ρ=0.04, c=0.06)

下载: 全尺寸图片幻灯片

将修正IDD理论预测模型应用到同时含有孔洞与裂隙缺陷的砂浆等效导热系数预测中，并与数值模拟计算结果进行比较。图15为孔隙率一定时，孔洞与裂隙相对含量变化对砂浆导热系数的影响。结果显示，修正IDD解与数值模拟解吻合良好。在孔隙率一定时，随着孔洞率与裂隙密度的比值增加，砂浆的导热系数逐渐减小，说明孔洞对砂浆导热系数的影响较裂隙显著。

图 15 孔洞率 c 与裂隙密度 ρ 之比对水泥砂浆导热系数影响

Figure 15. Effect of the ratio of porosity c to crack density ρ on thermal conductivity of cement mortar

下载: 全尺寸图片幻灯片

4. 结论

(1) 应用细观力学基础理论，基于相互作用直推(IDD)估计框架，分别推导出含有孔洞与裂隙两种形式缺陷夹杂的等效导热系数预测模型。根据数值模拟反推出的裂隙随机分布影响函数具有较高精度和良好适用性。

(2) 该预测模型，充分体现了孔洞、裂隙不同含量及特征对水泥砂浆等效导热系数的影响，预测精度较高；在缺陷率一定的情况下，水泥砂浆等效导热系数随孔洞率增加而减小，孔洞对水泥砂浆等效导热系数的影响比裂隙更显著。

(3) 推导的预测模型物理基础明确，形式简单，便于工程应用，可为研究含有一定损伤的水泥基材料等效导热系数提供理论支撑，同时可为含有线状夹杂和球状夹杂的复合材料导热系数研究提供参考。

图 1 含孔洞和裂隙的水泥砂浆混合缺陷模型

Figure 1. Mixed defect model with cavities and cracks of cement mortar

下载: 全尺寸图片幻灯片

图 2 椭球夹杂示意图

Figure 2. Schematic diagram of ellipsoid inclusion

下载: 全尺寸图片幻灯片

图 3 干燥水泥砂浆网格划分图

Figure 3. Meshing diagram of dry cement mortar

下载: 全尺寸图片幻灯片

图 4 随机缺陷程序流程图

Figure 4. Flow chart of random defect program

下载: 全尺寸图片幻灯片

图 5 平行裂隙模型(ρ=0.05、θ=60°)

Figure 5. Parallel crack model (ρ=0.05, θ=60°)

下载: 全尺寸图片幻灯片

图 6 随机裂隙模型(ρ=0.05)

Figure 6. Random crack model (ρ=0.05)

下载: 全尺寸图片幻灯片

图 7 裂隙密度对水泥砂浆正交方向导热系数的影响

Figure 7. Influence of crack density on the thermal conductivity of cement mortar in the orthogonal direction

CIDD—Correct IDD solution; IDD—Interaction direct deduction solution; NS—Numerical simulation solution

下载: 全尺寸图片幻灯片

图 8 裂隙方向对水泥砂浆导热系数的影响

Figure 8. Influence of crack direction on thermal conductivity of cement mortar

下载: 全尺寸图片幻灯片

图 9 裂隙长度与数量对水泥砂浆导热系数的影响

Figure 9. Influence of crack length and quantity on thermal conductivity of cement mortar

下载: 全尺寸图片幻灯片

图 10 随机裂隙密度对水泥砂浆导热系数的影响

Figure 10. Effect of random crack density on thermal conductivity of cement mortar

L—Crack length

下载: 全尺寸图片幻灯片

图 11 随机孔洞模型(c=0.05)

Figure 11. Random cavity model (c=0.05)

下载: 全尺寸图片幻灯片

图 12 不同孔洞特征对水泥砂浆导热系数的影响

Figure 12. Influence of different cavity characteristics on thermal conductivity of cement mortar

d—Cavity diameter

下载: 全尺寸图片幻灯片

图 13 孔洞率对水泥砂浆导热系数的影响

Figure 13. Influence of porosity on thermal conductivity of cement mortar

下载: 全尺寸图片幻灯片

图 14 含孔洞和裂隙缺陷模型(ρ=0.04、c=0.06)

Figure 14. Defect model with cavities and cracks (ρ=0.04, c=0.06)

下载: 全尺寸图片幻灯片

图 15 孔洞率 c 与裂隙密度 ρ 之比对水泥砂浆导热系数影响

Figure 15. Effect of the ratio of porosity c to crack density ρ on thermal conductivity of cement mortar

下载: 全尺寸图片幻灯片

参考文献(24)

[1]	陈则韶, 钱军, 叶一火. 复合材料等效导热系数的理论推算[J]. 中国科学技术大学学报, 1992, 22(4):416-424. CHEN Z S, QIAN J, YE Y H. Predicting theory of effective thermal conductivity of complex material[J]. Journal of University of Science and Technology of China,1992,22(4):416-424(in Chinese).
[2]	MAXWELL J C. A treatise on electricity and magnetism[M]. New York: Dover Publication, 1954: 16-20.
[3]	张伟平, 童菲, 邢益善, 等. 混凝土导热系数的试验研究与预测模型[J]. 建筑材料学报, 2015, 18(2):183-189. DOI: 10.3969/j.issn.1007-9629.2015.02.001 ZHANG W P, TONG F, XING Y S, et al. Experimental study and prediction model of thermal conductivity of concrete[J]. Journal of Building Materials,2015,18(2):183-189(in Chinese). DOI: 10.3969/j.issn.1007-9629.2015.02.001
[4]	张伟平, 邢益善, 顾祥林. 基于细观复合材料的混凝土导热系数模型[J]. 结构工程师, 2012, 28(2):39-45. DOI: 10.3969/j.issn.1005-0159.2012.02.007 ZHANG W P, XING Y S, GU X L. Theoretical models of effective thermal conductivity of concrete based on composite materials in mesoscale[J]. Structural Engineers,2012,28(2):39-45(in Chinese). DOI: 10.3969/j.issn.1005-0159.2012.02.007
[5]	姚秀鹏, 韩阳, 沈雷, 等. 高温后聚丙烯纤维增强水泥基复合材料导热的多尺度方法[J]. 复合材料学报, 2021, 38(10):3531-3542. YAO X P, HAN Y, SHEN L, et al. Multi-scale method for thermal conductivity of polypropylene fiber reinforced cementitious composites after high temperature[J]. Acta Materiae Compositae Sinica,2021,38(10):3531-3542(in Chinese).
[6]	王立成, 常泽, 鲍玖文. 基于多相复合材料的混凝土导热系数预测模型[J]. 水利学报, 2017, 48(7):765-772. WANG L C, CHANG Z, BAO J W. Prediction model for the thermal conductivity of concrete based on its composite structure[J]. Journal of Hydraulic Engineering,2017,48(7):765-772(in Chinese).
[7]	刘嘉涵, 徐世烺, 曾强. 基于多尺度细观力学方法计算水泥基材料的导热系数[J]. 建筑材料学报, 2018, 21(2):293-298. DOI: 10.3969/j.issn.1007-9629.2018.02.019 LIU J H, XU S L, ZENG Q. An investigation of thermal conductivity of cement-based composites with multi-scale micromechanical method[J]. Journal of Building Materials,2018,21(2):293-298(in Chinese). DOI: 10.3969/j.issn.1007-9629.2018.02.019
[8]	ESHELBY J D. The determination of the elastic field of an ellipsoidal inclusion, and related problems[J]. Proceedings of the Royal Society of London, 1957, 241(1226): 376−396.
[9]	张研, 韩林. 细观力学基础[M]. 北京: 科学出版社, 2014: 78-90. ZHANG Y, HAN L. Foundation of mesomechanics[M]. Beijing: Science Press, 2014: 78-90 (in Chinese).
[10]	钟轶峰, 张亮亮, 周小平, 等. 复合材料热传导性能的变分渐近均匀化细观力学模型[J]. 复合材料学报, 2015, 32(4):1173-1178. ZHONG Y F, ZHANG L L, ZHOU X P, et al. Variational asymptotic homogenization micromechanics model for thermal conductivity of composites[J]. Acta Materiae Compositae Sinica,2015,32(4):1173-1178(in Chinese).
[11]	孙海浩, 石义雷, 刘伟强, 等. 随机结构复合材料等效导热系数的理论预测[J]. 复合材料学报, 2021, 38(9):2918-2926. SUN H H, SHI Y L, LIU W Q, et al. Theoretical prediction for effective thermal conductivity of composite materials with random structure[J]. Acta Materiae Compositae Sinica,2021,38(9):2918-2926(in Chinese).
[12]	张枫, 肖建庄, 宋志文. 混凝土导热系数的理论模型及其应用[J]. 商品混凝土, 2009(2):23-25, 51. ZHANG F, XIAO J Z, SONG Z W. The prediction models of thermal conductivity of concrete and their application[J]. Ready-Mixed Concrete,2009(2):23-25, 51(in Chinese).
[13]	唐世斌, 唐春安, 梁正召, 等. 混凝土热传导与热应力的细观特性及热开裂过程研究[J]. 土木工程学报, 2012, 45(2):732-747. TANG S B, TANG C A, LIANG Z Z, et al. Study of thermal conduction and thermal stress of concrete at mesoscopic level and its thermal cracking processes[J]. China Civil Engineering Journal,2012,45(2):732-747(in Chinese).
[14]	孙颖颖, 周璐瑶, 韩宇, 等. 气泡和气隙影响六方氮化硼/环氧树脂复合材料导热性能的有限元模拟[J]. 复合材料学报, 2020, 37(10):2482-2488. SUN Y Y, ZHOU L Y, HAN Y, et al. Numerical analysis of the effect of air bubbles and gaps on thermal conductivity of hexagonal boron nitride/epoxy composites[J]. Acta Materiae Compositae Sinica,2020,37(10):2482-2488(in Chinese).
[15]	李守巨, 刘迎曦, 于贺. 多孔材料等效导热系数与分形维数关系的数值模拟研究[J]. 岩土力学, 2009, 30(5):1465-1470. DOI: 10.3969/j.issn.1000-7598.2009.05.049 LI S J, LIU Y X, YU H. Numerical simulation of relationship between thermal conductivity of porous material and fractal dimension[J]. Rock and Soil Mechanics,2009,30(5):1465-1470(in Chinese). DOI: 10.3969/j.issn.1000-7598.2009.05.049
[16]	SHAFIRO B, KACHANOV M. Anisotropic effective conductivity of materials with nonrandomly oriented inclusions of diverse ellipsoidal shapes[J]. Journal of Applied Physics,2000,87(12):8561-8569. DOI: 10.1063/1.373579
[17]	ZHENG Q S, DU D X. An explicit and universally applicable estimate for the effective properties of multiphase composites which accounts for inclusion distribution[J]. Journal of the Mechanics and Physics of Solids,2001,49:2765-2788. DOI: 10.1016/S0022-5096(01)00078-3
[18]	ZHOU C S, LI K F, PANG X Y. Effect of crack density and connectivity on the permeability of microcracked solids[J]. Mechanics of Materials,2011,43:969-978. DOI: 10.1016/j.mechmat.2011.08.011
[19]	周春圣, 李克菲. 含裂纹夹杂多孔材料的渗透性理论与数值分析[J]. 工程力学, 2013, 30(4):150-156. DOI: 10.6052/j.issn.1000-4750.2011.11.0753 ZHOU C S, LI K F. Theoretical and numerical analyses of permeability of porous medium with cracking inclusions[J]. Engineering Mechanics,2013,30(4):150-156(in Chinese). DOI: 10.6052/j.issn.1000-4750.2011.11.0753
[20]	HASHIN Z. Assessment of the self-consistent scheme approximation: Conductivity of particulate composites[J]. Journal of Composite Materials,1968,2(3):284-300. DOI: 10.1177/002199836800200302
[21]	DU D X, ZHENG Q S. A further exploration of the interaction direct derivative (IDD) estimate for the effective properties of multiphase composites taking into account inclusion distribution[J]. Acta Mechanics,2002,157(1):61-80.
[22]	杨世铭, 陶文铨. 传热学[M]. 北京: 高等教育出版社, 2006: 15-36. YANG S M, TAO W Q. Heat transfer[M]. Beijing: Higher Education Press, 2006: 15-36 (in Chinese).
[23]	KACHANOV M. Effective elastic properties of cracked solids: Critical review of some basic concepts[J]. Applied Mechanics Reviews,1992,45 (8):304–335.
[24]	朱伯芳. 大体积混凝土温度应力与温度控制[M]. 北京: 中国电力出版社, 2012: 8-12. ZHU B F. Thermal stresses and temperature control of mass concrete[M]. Beijing: China Electric Power Press, 2012: 8-12(in Chinese).

施引文献(6)

期刊类型引用(4)

1.	万佩，夏辉，刘晨，贾吉龙，何学，丁安心. 基于多尺度数值模型的复合材料各向异性热膨胀系数预测. 复合材料学报. 2023(02): 1208-1217 . 本站查看
2.	朱丽华，刘海林，韩伟. 基于细观尺度的再生混凝土多相导热系数理论模型. 材料导报. 2023(12): 117-123 . 百度学术
3.	王兴禹. 不同船型海上水泥运输与过驳的经济可行性分析. 中国航务周刊. 2022(18): 51-53 . 百度学术
4.	张凡，米晶晶，李泽树，陈涛，胡魁. 水泥砂浆外墙装饰施工质量控制探讨. 散装水泥. 2022(06): 18-19 . 百度学术