CN111678945A - 一种考虑骨料形状的混凝土导热系数的计算方法 - Google Patents
一种考虑骨料形状的混凝土导热系数的计算方法 Download PDFInfo
- Publication number
- CN111678945A CN111678945A CN202010584348.5A CN202010584348A CN111678945A CN 111678945 A CN111678945 A CN 111678945A CN 202010584348 A CN202010584348 A CN 202010584348A CN 111678945 A CN111678945 A CN 111678945A
- Authority
- CN
- China
- Prior art keywords
- concrete
- formula
- phase
- temperature gradient
- aggregate
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
- G01N25/20—Investigating or analyzing materials by the use of thermal means by investigating the development of heat, i.e. calorimetry, e.g. by measuring specific heat, by measuring thermal conductivity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
- G01N25/18—Investigating or analyzing materials by the use of thermal means by investigating thermal conductivity
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
Description
技术领域
本发明属于材料和土木、水利工程技术领域,具体涉及一种考虑骨料形状的混凝土导热系数的计算方法。
背景技术
掌握混凝土的导热性能对于地热工程、核电站工程等的安全运行具有至关重要的意义。一般情况下,根据毛细管中是否充满了水,可将混凝土分为干燥、非饱和及饱和状态。为了定量研究混凝土中的热传递过程,减少复杂、费力、耗时的试验和数值模拟,建立准确预测混凝土有效导热系数,为工程的安全运行以及高性能混凝土材料的设计具有重要意义和必要性。
现有混凝土有效导热系数主要有三方面的计算方案:(1)基于实验的半经验模型,该方法包含经验参数,具有较强的经验性;(2)基于数值模拟的数学模型,该方法的精确性在一定程度上依赖于数值模拟结果的准确性,模型缺乏物理意义,不能很好的揭示机理;(3)基于不同理论的分析模型。有效导热系数的理论模型主要可以分为两类:其一,毛细管束模型。该方法是将孔视为连通上下边界的管道,由于分形几何相对于欧几里得传统几何能更为真实的反应毛细管的自相似性,因此,基于分形几何的毛细管束模型更常被采用。其二,细观力学模型,该方法是将混凝土视为夹杂-基体的系统。总的来说,现有理论模型没有将骨料与饱和度的统一考虑,例如,仅针对于饱和问题或非饱和不考虑骨料问题等。由于现有理论模型的局限性,不能很好的适应实际工程的真实性及高性能混凝土设计的复杂性。
已有的最接近的技术方案是:Li等通过Mori-Tanaka方法,推导了非饱和水泥基(不考虑骨料的存在)的有效导热系数模型,另外,该课题组基于多尺度细观力学和Mori-Tanaka方法建立了干燥与饱和情况下考虑骨料的有效导热系数模型。本发明模型的主要优势:本发明模型兼顾了最接近方案的两种情况,能综合分析骨料存在且混凝土处在不同饱和度情况下(包括Li等考虑的干燥与饱和情况)的热传导问题。
发明内容
发明目的:本发明提供的一种考虑骨料形状的混凝土导热系数的计算方法,可节省大量时间、人力和物力,且预测结果与实验吻合度高。
技术方案:本发明所述的一种考虑骨料形状的混凝土导热系数的计算方法,包括如下步骤:
(1)计算各相的温度梯度和体积分数;
(2)引入了二阶浓度因子Ar,建立各相的温度梯度与体积平均温度梯度的联系;
(4)从单夹杂推广到多相夹杂,考虑饱和度和各相(骨料相、液相和气相) 及其形状等因素的耦合作用;
(5)得到考虑骨料存在时不同饱和度的混凝土有效导热系数。
进一步地,所述步骤(1)的实现过程如下:
对于无源温度场,热通量q和温度梯度T的函数表达式如下:
式中,K为导热系数二阶张量,T为温度,xj为坐标,负号表示热传导是沿温度梯度减小的方向,为了在混凝土中产生统计均匀温度场,整体边界条件为T =-T0xj,其中T0为整体温度梯度,且体积平均温度梯度等于T0,和体积平均热通量可以写为:
式中,下标r、m分别是第r相夹杂与基体;p,l,g分别是骨料、水和气; Vr和Vm是第r相夹杂和基体的体积分数;根据平均场理论和公式(1),各相和混凝土的热通量可表示为:
式中:Keff是混凝土的有效导热系数;将式(4)代入式(3)有:
式中:Kp,Kl,Kg和Km分别是骨料、水、气和基体的导热系数;将式(2) 代入式(5)中有:
式中:I是单位张量。
进一步地,所述步骤(2)通过以下公式实现:
其中,Ar是二阶浓度因子。
4、根据权利要求1所述的一种考虑骨料形状的混凝土导热系数的计算方法,其特征在于,所述步骤(3)实现过程如下:
式中:Pr为去极化因子,可以写为:
式中:χ为椭球的长径比χ=a3/a1,a1和a3是椭球的特征方向的半径,χ>1, χ=1和χ<1表示扁长形、球形和扁圆形,长纤维时(χ→∞)有Qr=1/2,薄片形(χ=0)有Qr=0。
进一步地,所述步骤(4)实现过程如下:
各相的浓度因子为:
非饱和混凝土的有效导热系数为:
由于液相和气相的体积分数之和就是混凝土的孔隙率φ,即:
φ=Vl+Vg (18)
饱和度用来定量分析热传导与体积含水量θ的关系,其定义如下:
式中:θres和θsat为剩余体积含水量和饱和体积含水量;初始情况被视为完全干燥,即θres=0;另外,在这里认为水和气的密度与饱和度无关,则Θ可写为:
Θ=Vl/φ (20a)
1-Θ=Vg/φ (20b) 将式(20)与式(17)中,则有:
进一步地,所述步骤(5)通过以下公式实现:
其中,干燥状态为Θ=0;饱和状态为Θ=1。
有益效果:与现有技术相比,本发明的有益效果:1、有限元软件计算混凝土导热系数的操作繁琐,利用本发明可节省大量时间、人力和物力,且预测结果与实验吻合度高;2、针对现有模型无法兼顾骨料与饱和度问题,提出了一个细观力学理论框架,用于准确预测干燥、非饱和及饱和状态混凝土的有效导热系数,为混凝土的测试和设计提供了一个重要的理论工具,以便有效和全面地了解和分析各种关键控制变量的影响;3、本发明为混凝土有效导热系数的计算和预测提供了理论基础,用于设计高性能混凝土材料和预测混凝土结构的使用周期寿命,具有较高科研和实际应用价值。
附图说明
图1为本发明的流程图;
图2为非饱和混凝土的归一化有效导热系数K+与饱和度Θ的预测效果图;
图3为干燥混凝土的有效导热系数Keff与骨料体积分数Vp的预测效果图;
图4为非饱和含骨料混凝土的有效导热系数Keff与饱和度Θ的预测效果图;
图5为干燥状态混凝土的有效导热系数Keff与骨料体积分数关系预测图;
图6为饱和状态混凝土的有效导热系数Keff与骨料体积分数关系预测图。
具体实施方式
下面结合附图对本发明作进一步详细说明。
针对现有模型无法兼顾骨料与饱和度问题,本发明提出了一个细观力学理论框架下对干燥、非饱和及饱和状态混凝土有效导热系数的预测模型,如图1所示,具体包括以下步骤:
步骤1:计算各相的温度梯度和体积分数。
热传导是由温度梯度所驱动。对于无源温度场,热通量q和温度梯度T的函数表达式如下:
式中:K为导热系数二阶张量,T为温度,xj为坐标,负号表示热传导是沿温度梯度减小的方向。为了在混凝土中产生统计均匀温度场,整体边界条件为T =-T0xj,其中T0为整体温度梯度,且体积平均温度梯度等于T0。和体积平均热通量可以写为:
式中:下标r、m分别是第r相夹杂与基体;p,l,g分别是骨料、水和气; Vr和Vm是第r相夹杂和基体的体积分数。根据平均场理论和公式(1),各相和混凝土的热通量可表示为:
式中:Keff是混凝土的有效导热系数。将式(4)代入式(3)有:
式中:Kp,Kl,Kg和Km分别是骨料、水、气和基体的导热系数。将式(2) 代入式(5)中有:
式中:I是单位张量。
步骤2:引入了二阶浓度因子Ar,建立各相的温度梯度与体积平均温度梯度的联系。
上式表明求解Ar是计算Keff的关键。
考虑一个在r相中热导系数为Kr的椭球夹杂,它被嵌入一个具有温度梯度 T1和热导系数K1的无限均匀基体中。
式中:Pr为去极化因子,可以写为:
式中:χ为椭球的长径比χ=a3/a1,a1和a3是椭球的特征方向的半径。χ>1, χ=1和χ<1表示扁长形、球形和扁圆形,长纤维时(χ→∞)有Qr=1/2,薄片形(χ=0)有Qr=0。
步骤4:从单夹杂推广到多相夹杂,考虑饱和度和各相(骨料相、液相和气相)及其形状等因素的耦合作用。
在本技术方案中,混凝土包含三类夹杂物:骨料、水和气。由于是统计均匀的材料,在同一相中的夹杂物具有相同的温度梯度、导热系数和长径比,因此式 (10)可进一步写为:
依据式(13)和式(15)和Ar的定义,各相的浓度因子可以写为:
将式(16)代入式(8)中,非饱和混凝土的有效导热系数可以写为:
由于液相和气相的体积分数之和就是混凝土的孔隙率φ,即:
φ=Vl+Vg (18)
另一方面,饱和度Θ用来定量分析热传导与体积含水量θ的关系,其定义如下:
式中:θres和θsat为剩余体积含水量和饱和体积含水量。通常情况下,初始情况被视为完全干燥,即θres=0。另外,在这里认为水和气的密度与饱和度无关,则Θ可写为:
Θ=Vl/φ (20a)
1-Θ=Vg/φ (20b) 将式(20)代入式(17)中,则有:
将其写为标量的形式为:
综上:依据式(22)、(23)可以计算考虑骨料存在时不同饱和度的混凝土有效导热系数,其中,干燥状态为Θ=0;饱和状态为Θ=1。
采用Jerman等关于非饱和水泥基材料(无骨料情况)的实验数据进行第一组对比中。首先令K+=Keff/Kg进行归一化处理。依据实验参数输入参数为:Kp= 1.322Wm-1K-1,Kl=0.6065Wm-1K-1,Kg=0.02623Wm-1K-1,χp=0.5,χg=0.1, 0.12和0.18分别对应于φ=0.874,0.819和0.802。预测结果与实验结果对比如图 2所示,可以看到利用本发明模型所预测的有效导热系数与实验结果具有较好的一致性。
采用Zhang等关于干燥混凝土不同骨料体积分数的实验数据,进行第二组对比中。输入参数为:Kp=2.6Wm-1K-1,Km=1.62Wm-1K-1,χp=2.5,χg=0.6,φ=19.74%,Vp分别为0.3,0.398,0.5,0.6。由图3可知,本发明模型针对干燥混凝土的骨料体积分数与气相耦合作用具有较好的预测效果。
采用Zhang等关于非饱和混凝土不同骨料体积分数的实验数据进行第三组对比中。输入参数为:Kp=2.6Wm-1K-1,Km=1.62Wm-1K-1,Kl=0.6065Wm-1K-1, Kg=0.02623Wm-1K-1,φ=19.74%,Vp=0.398,χp=2.5,χl=0.5,χg=0.11。由图4可知,本发明模型针对非饱和混凝土多相(骨料相、气相和液相)耦合作用的总体较好,但在饱和度很高时出现偏差,这是由于骨料的导热系数是随饱和度动态变化的,而饱和度较低时,由于混凝土中存在大量极易吸水的胶凝孔,以及骨料周围存在大孔径的界面过渡区(ITZ),使得骨料接近于干燥状态,此时本发明模型有较高的预测精度,而在饱和度较高时,由于骨料迅速饱和,导致骨料的导热系数Kp大于2.6Wm-1K-1,导致总体有效导热系数的实际值偏高。
采用Kim等关于干燥和饱和混凝土不同骨料体积分数的实验数据,进行第四组对比。值得注意的是,该实验没有测量孔隙率,但是提供了骨料体积分数为0 时的导热系数,即该组混凝土可视为基体和骨料的两相材料。干燥状态的输入参数为:Km=0.768Wm-1K-1,Kp=3.0Wm-1K-1,χp=2.5;饱和状态时的输入参数为:Km=1.163Wm-1K-1,Kp=3.4Wm-1K-1,χp=2.5。由图5、图6可知,本发明模型针对干燥和饱和状态的骨料-基体两相系统具有很有的预测效果,说明了该发明模型的灵活性。另一方面,证明了干燥与饱和状态时的骨料导热系数存在一定差异,为图3中饱和度较高时实验导热系数是由骨料导热系数变化所引起的提供理论了支撑。
综上所述,本发明模型能准确预测考虑骨料存在时的干燥、非饱和及饱和状态混凝土的有效导热系数。另外,还能通过本发明模型还能对骨料的特性(例如:长径比、导热系数和体积分数)和气相、液相的形状以及基体的导热系数等参数进行影响定量分析。
Claims (6)
2.根据权利要求1所述的一种考虑骨料形状的混凝土导热系数的计算方法,其特征在于,所述步骤(1)的实现过程如下:
对于无源温度场,热通量q和温度梯度T的函数表达式如下:
式中,K为导热系数二阶张量,T为温度,xj为坐标,负号表示热传导是沿温度梯度减小的方向,为了在混凝土中产生统计均匀温度场,整体边界条件为T=-T0xj,其中T0为整体温度梯度,且体积平均温度梯度等于T0,和体积平均热通量可以写为:
式中,下标r、m分别是第r相夹杂与基体;p,l,g分别是骨料、水和气;Vr和Vm是第r相夹杂和基体的体积分数;根据平均场理论和公式(1),各相和混凝土的热通量可表示为:
式中:Keff是混凝土的有效导热系数;将式(4)代入式(3)有:
式中:Kp,Kl,Kg和Km分别是骨料、水、气和基体的导热系数;将式(2)代入式(5)中有:
式中:I是单位张量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010584348.5A CN111678945B (zh) | 2020-06-24 | 2020-06-24 | 一种考虑骨料形状的混凝土导热系数的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010584348.5A CN111678945B (zh) | 2020-06-24 | 2020-06-24 | 一种考虑骨料形状的混凝土导热系数的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111678945A true CN111678945A (zh) | 2020-09-18 |
CN111678945B CN111678945B (zh) | 2023-07-21 |
Family
ID=72456515
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010584348.5A Active CN111678945B (zh) | 2020-06-24 | 2020-06-24 | 一种考虑骨料形状的混凝土导热系数的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111678945B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112630254A (zh) * | 2020-12-11 | 2021-04-09 | 江苏华彤建设工程有限公司 | 基于变温下热寿命损耗的母线槽寿命预测方法 |
CN114835508A (zh) * | 2022-04-14 | 2022-08-02 | 湖南工程学院 | 加气混凝土的检测方法及制备方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08240548A (ja) * | 1995-03-03 | 1996-09-17 | Fujita Corp | コンクリート構造物の躯体内部の含水率測定方法 |
CN101609057A (zh) * | 2009-07-15 | 2009-12-23 | 河海大学 | 普通混凝土导热系数随龄期发展预测模型的构建方法 |
CN102162802A (zh) * | 2011-01-11 | 2011-08-24 | 河海大学 | 全级配混凝土导热系数预测模型的构建方法 |
CN104133052A (zh) * | 2014-07-30 | 2014-11-05 | 中国水利水电科学研究院 | 一种大体积混凝土温控效果全过程评价方法 |
CN107402227A (zh) * | 2017-07-14 | 2017-11-28 | 天津城建大学 | 土体矿物导热系数的计算方法 |
RU2647546C1 (ru) * | 2016-12-22 | 2018-03-16 | Федеральное государственное унитарное предприятие "Институт химических реактивов и особо чистых химических веществ Национального исследовательского центра "Курчатовский институт" | Способ определения устойчивости покрытий и устройство для его осуществления |
CN108828002A (zh) * | 2018-04-17 | 2018-11-16 | 辽宁工程技术大学 | 一种固液气三相非饱和混凝土有效导热系数预测方法 |
CN111027236A (zh) * | 2019-10-30 | 2020-04-17 | 武汉大学 | 多相复合材料的等效导热系数的细观尺度研究方法 |
-
2020
- 2020-06-24 CN CN202010584348.5A patent/CN111678945B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08240548A (ja) * | 1995-03-03 | 1996-09-17 | Fujita Corp | コンクリート構造物の躯体内部の含水率測定方法 |
CN101609057A (zh) * | 2009-07-15 | 2009-12-23 | 河海大学 | 普通混凝土导热系数随龄期发展预测模型的构建方法 |
CN102162802A (zh) * | 2011-01-11 | 2011-08-24 | 河海大学 | 全级配混凝土导热系数预测模型的构建方法 |
CN104133052A (zh) * | 2014-07-30 | 2014-11-05 | 中国水利水电科学研究院 | 一种大体积混凝土温控效果全过程评价方法 |
RU2647546C1 (ru) * | 2016-12-22 | 2018-03-16 | Федеральное государственное унитарное предприятие "Институт химических реактивов и особо чистых химических веществ Национального исследовательского центра "Курчатовский институт" | Способ определения устойчивости покрытий и устройство для его осуществления |
CN107402227A (zh) * | 2017-07-14 | 2017-11-28 | 天津城建大学 | 土体矿物导热系数的计算方法 |
CN108828002A (zh) * | 2018-04-17 | 2018-11-16 | 辽宁工程技术大学 | 一种固液气三相非饱和混凝土有效导热系数预测方法 |
CN111027236A (zh) * | 2019-10-30 | 2020-04-17 | 武汉大学 | 多相复合材料的等效导热系数的细观尺度研究方法 |
Non-Patent Citations (5)
Title |
---|
HEDONG LI 等: ""Effect of pore shape on the thermal conductivity of partially saturated cement-based porous composites"", 《CEMENT AND CONCRETE COMPOSITES》 * |
HIROSHI HATTA 等: ""Effective thermal conductivity of a misoriented short fiber composite"", 《JOURNAL OF APPLIED PHYSICS》 * |
刘凯等: "水泥混凝土导热性能影响因素及预估模型研究", 《建筑材料学报》 * |
张仁波等: "混凝土温度相关热传导行为细观分析", 《北京工业大学学报》 * |
张伟平等: "混凝土导热系数的试验研究与预测模型", 《建筑材料学报》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112630254A (zh) * | 2020-12-11 | 2021-04-09 | 江苏华彤建设工程有限公司 | 基于变温下热寿命损耗的母线槽寿命预测方法 |
CN112630254B (zh) * | 2020-12-11 | 2022-03-08 | 江苏华彤建设工程有限公司 | 基于变温下热寿命损耗的母线槽寿命预测方法 |
CN114835508A (zh) * | 2022-04-14 | 2022-08-02 | 湖南工程学院 | 加气混凝土的检测方法及制备方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111678945B (zh) | 2023-07-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ghosh et al. | A multi-level computational model for multi-scale damage analysis in composite and porous materials | |
Karakoti et al. | Nonlinear transient analysis of porous P-FGM and S-FGM sandwich plates and shell panels under blast loading and thermal environment | |
Jiang et al. | Probability-interval hybrid reliability analysis for cracked structures existing epistemic uncertainty | |
Yan et al. | A new 2D continuous-discontinuous heat conduction model for modeling heat transfer and thermal cracking in quasi-brittle materials | |
Sommer et al. | Coupling anisotropic viscosity and fiber orientation in applications to squeeze flow | |
Zhou et al. | Numerical simulation of thermal cracking in rocks based on general particle dynamics | |
CN111678945A (zh) | 一种考虑骨料形状的混凝土导热系数的计算方法 | |
Zhang et al. | Multiscale modeling of the effect of the interfacial transition zone on the modulus of elasticity of fiber-reinforced fine concrete | |
Sun et al. | A new anisotropic thermal conductivity equation for h-BN/polymer composites using finite element analysis | |
Schneider | An FFT‐based method for computing weighted minimal surfaces in microstructures with applications to the computational homogenization of brittle fracture | |
Zhao et al. | Parallelization of the distinct lattice spring model | |
Martinie et al. | Fiber orientation during casting of UHPFRC: electrical resistivity measurements, image analysis and numerical simulations | |
Duan et al. | Numerical simulation of the freezing process of concrete | |
Rukhaiyar et al. | Strength behavior of plain cement concrete subjected to true triaxial compression | |
Bennett et al. | Effective thermoelasticity of polymer-bonded particle composites with imperfect interfaces and thermally expansive interphases | |
Cheng et al. | New technique for frictional contact on crack slip in the extended finite-element method framework | |
Wang et al. | Density functional theory-based database development and CALPHAD automation | |
Sun et al. | A coupled 3D thermo-mechanical peridynamic model for cracking analysis of homogeneous and heterogeneous materials | |
Liu et al. | XFEM for thermal crack of massive concrete | |
Liu et al. | Surface effect on a nano‐elliptical hole or nano‐crack in magnetoelectroelastic materials under antiplane shear | |
Jain et al. | Homogenization based 3d continuum damage mechanics model for composites undergoing microstructural debonding | |
Tran et al. | Extended finite element method for plastic limit load computation of cracked structures | |
Liu et al. | A stochastic multi-scale model for predicting the thermal expansion coefficient of early-age concrete | |
Li et al. | A meshless model for transient heat conduction analyses of 3D axisymmetric functionally graded solids | |
Peralta et al. | Multiscale homogenization analysis of the effective elastic properties of masonry structures |
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 |