水化热参数化分析

1、midas FEA Training Series水化热参数化分析midas FEA Training Series一概要1. 水化热分析浇筑混凝土时，水泥在水化过程中产生大量热量会使混凝土的温度升高。虽然随时间的推移混凝土的温度会慢慢冷却，但结构各个位置的温度下降速度不均匀，结构不同位置将发生相对温差，此温差会使混凝土发生温度应力。温度裂缝发生类型混凝土浇筑初期，因内部温度升高将发生膨胀，但混凝土表面的温度下降较快，相对应变较小，从而使混凝土表面产生拉应力。混凝土内部不同的温度分布引起的不同的体积变化而导致的应力称为内部约束应力。此类拉应力裂缝主要发生在构件尺寸比较大的结构。| 内部约束产生

2、的裂缝（放热时）| 外部约束产生的裂缝（冷却时）| 混凝土在高温状态下温度下降会发生收缩，但受到与其接触的已浇筑混凝土或者地基等的约束而产生的拉力，像这样变形受外部边界约束的状态称为外部约束。此类应力主要发生在像墙这样约束度比较大的结构中。 利用温度裂缝指数预测温度裂缝韩国混凝土规范中使用温度裂缝指数（抗拉强度与发生的温度应力之比）i 值预测是否发生裂缝。一般采用下面的值。 |温度裂缝指数与裂缝发生几率 |裂缝指数(i) =混凝土抗拉强度发生的温度应力 防止裂缝发生时：1.5 以上 限制裂缝发生时：1.2 1.5 限制有害裂缝发生时：0.7 1.2FEA程序的水化热分析水化热分析主要分为热传导

3、分析和热应力分析。. 热传导分析主要计算水泥的水化过程中发热、传导、对流等引起的随时间变化的节点温度。将得到的节点温度作为荷载加载后，计算随时间变化的应力称为热应力分析。 因此通过查看温度分布可以看出输入数据是否有误，如果温度分布没有问题可说明输出的应力结果也是正确的。 2. 水化热参数化分析 输入混凝土的散热特性及浇筑条件等混凝土的温度应力裂缝指数 OK?ENDYesNo水化热分析必须进行反复计算大体积混凝土的温度裂缝可以利用温度裂缝指数(Crack Ratio, Icr) 来验算。温度裂缝指数要满足结构的重要性、功能、环境条件等因素的要求。 温度裂缝指数受水泥的类型、浇筑温度、养生方法等多

4、因素的影响，所以需要对多种条件进行反复分析以找出最佳的浇筑方法。 参数化分析功能为比较多种条件的分析结果需要建立多个模型进行分析，分析结束后需要整理大量的分析结果、还要进行结果保存、对比等工作。 通过FEA的水化热参数化分析功能，可以实现一个模型多种条件分析。可以大大减少单纯繁琐的反复分析过程，从而提高工作效率。参数化分析的使用方法首先建立一个基本模型，在基本模型里使用替换变量的方式定义分析工况。下图是把材料作为变量条件的示例，“Case I”为将混凝土C24变更为C30的工况，“Case II”为将混凝土C35变更为C40的工况。 | 参数化分析的构成 |参数化分析里可以考虑的变量在水化热参

5、数化分析的功能里可以调整的变量有五个，较常用的调整方法具体如下。 施工阶段： 降低浇筑高度缩小各阶段的温度差。浇筑间距过小的话很难达到分段浇筑的效果，但如果太大分界面会产生较大的温差。. 对流边界：对流系数较低时，热量不容易对外流失，可以减少内外温差。 材料：使用弹性模量大的材料时，抗拉强度也较大，可增大裂缝指数。 发热特性：是变量中最为敏感的因素, 定义水化过程中发生的热量。 是否考虑自重：使混凝土产生压应力的荷载，在一定程度上可以减少拉应力，但效果不明显。74二建立基本模型1. 结构分析所需的数据水化热参数化分析如前面图形所示。首先建立基本模型，通过在基本模型里定义“Used”和“New”

6、的对应关系来定义分析工况。 | 进行参数化分析时输入变量的示意图|但是利用这种方法很难输入多个变量，虽然输入热源函数或对流系数函数很简单，但反复定义施工阶段和对流边界面的过程较繁琐。 为避免这种繁琐的定义过程，在定义水化热分析变量时，先不定义此两项。而是先定义阶段工况，即按不同的施工方案定义不同的施工阶段工况，然后再对不同的阶段工况定义各自的材料、对流系数等。| 进行参数化分析所输入的数据 |1) 建立模型此操作例题主要介绍“水化热参数化分析”的方法，仅对于相关变量的输入、各种分析条件的定义、查看分析结果部分进行详细说明。导入附件里的“HYD_Pier_Mesh.feb” 模型文件。n

7、88; 对称模型施工阶段水化热分析模型一般单元数量较多，所以分析所需的时间也较长，而且还要进行多条件分析工况的分析，所以需要更多的分析时间。如果模型属于对称模型，可以只建立¼模型以减少分析时间。 这样不但可以减少分析时间也有利于查看模型中心部位的分析结果。2) 特性1. 单位体系 :确认指定为 (N, m, J)分析 > 分析控制.操作步骤Procedure1分析 > 时间依存性材料> 徐变/收缩操作步骤Procedure1. 名称 : (C30)2. 规范 : (中国(JTG D62-2004)3. 混凝土28天材龄抗压强度 : (3.0e7)4. 点击 适用5.

8、 名称 : (C45)6. 混凝土28天材龄抗压强度 : (4.5e7)7. 点击 确认时间依存材料特性在FEA里可以考虑徐变/收缩特性和混凝土的抗压强度变化。 徐变计算方法可以选择用户定义或者使用规范的计算方法。 弹性模量折减方法是假设发生徐变，然后折减混凝土弹性模量的简易计算方法，一般的水化热分析里均使用这种方法。因弹性模量折减方法只适用于水化热分析，为了避免在一般的施工阶段分析中误用，在水化热分析控制里单独定义。 采用弹性模量折减方法时，为了要指定计算徐变的单元（通过材料对话框中的相应选项将徐变函数与材料连接起来，然后通过给单元赋予材料将徐变函数与单元连接起来），需要任意定义一个徐变计算

9、方法，但在这里定义的徐变计算方法并不参与水化热分析的徐变的计算。操作步骤分析 > 时间依存性材料 > 抗压强度Procedure1. 名称 : (C30_Normal)2. 类型 : 设计规范3. 规范 : ACI4. 混凝土28天抗压强度 (f28) : (3.0e7)5. a : (4.5)6. b : (0.95)7. 点击适用8. 重复上面的步骤定义C45_Normal 的强度发展函数。计算裂缝指数用的抗拉强度类型里选择“用户定义”后，用户可以通过表格自定义随时间变化的弹性模量、抗压强度、抗拉强度。 选择“设计规范”时，利用规范提供的公式计算弹性模量和抗拉强度。 利用抗拉强

10、度计算温度裂缝指数。 混凝土抗压强度系数与水泥的种类有关，一般硅酸盐水泥a为4.5, b为0.95。 序列号名称抗压强度ab1C30_Normal3.0e74.50.952C45_Normal4.5e74.50.951. 名称 : (Soil)2. 弹性模量 : (1e7)3. 泊松比 : (0.2)4. 膨胀系数. : (1e-5)5. 重量密度 : (26000)6. 点击 适用7. 参考下表输入C30_Normal和C45_Normal特性。分析 > 材料过程Procedure序列号名称弹性模量泊松比膨胀系数重量密度徐变/收缩抗压强度1Soil1e70.21e-526000-2C3

11、0_Normal3.0e100.21e-524517C30C30_Normal3C45_Normal3.35e100.21e-524517C45C45_Normal分析 > 特性操作步骤1. 点击 下拉菜单选择3D。2. 参考下表定义三个特性。Procedure2序列号名称材料1Soil1: Soil2Foundation2: C30_Normal3Pier3: C45_Normal1. 点击 “前视图”2. 特性 : (3D)3. 选择 “1:Soil”4. 选择 “3D Element(R)5. 选择地基1710个单元6. 点击适用7. 重复上述过程定义基础和桥墩的特性。分析 >

12、; 材料Procedure操作步骤145 地基-1: Soil 基础-2: Foundation 桥墩-3: Pier3) 边界条件分析 > 边界条件 > 约束操作步骤Procedure1. 点击 “前视图”2. 边界组 : Sym3. 选择Y-Z 对称面上的360个节点4. DOF : (T1)5. 点击 适用6. 点击 “左视图”7. 选择X-Z对称面上的705个节点8. 自由度 : (T2)9. 点击确认对称面边界条件对于Y-Z平面上的所有节点约束DX。37需要查看输入的对称边界条件时，可将网格显示为特征边线后查看。如下图所示。分析 > 边界 > 约束操作步骤Pr

13、ocedure1. 边界组 : Support2. 点击 “前视图”3. 选择地基的外部轮廓4. 点击 “左视图”5. 选择地基的外部轮廓6. 自由度 : (T1, T2, T3)实体单元没有旋转自由度所以只选择上述三个自由度。 7. 点击确认固结边界条件水化热分析模型里建立地基时，一般将地基下部的边界条件设为完全固结。混凝土产生的热量将充分地传递给地基，后续不再传递温度，也不存在温差，所以也不会发生相对位移。为了让混凝土产生的热量充分地传递给地基，需要建立足够大的地基模型。 552. 热传导分析所需的数据1) 材料的热特性数据定义水化过程产生的热传递的特性。模型单元内的热传导由比热和热传导率

14、确定，单元外的散热由下一节定义的对流系数来确定。分析 > 材料操作步骤1. 定义的材料列表里选择1:Soil 2. 点击 修改3. 点击 热工参数4. 传导率 : (3.45)5. 比热 : (784)6. 点击 确认7. 参考下表输入混凝土的热特性值。Procedure| 一般岩体和混凝土的热工系数 |热传导率 (W / m2 ·°C)比热 (kJ / kg · °C)岩体1.7 5.20.71 0.88混凝土2.6 2.81.05 1.26| 热工特性 |序列号材料名称热传导率比热1Soil3.457842C30_Normal2.711763C

15、45_Normal2.711762) 固定温度为了定义实体单元表面与外界的热传递关系，需要定义相应的对流边界条件。 固定温度边界、对流边界、绝热边界如下图所示。对流边界（热传导边界）通过对流系数定义各时刻与外部交换的热量。对于不同的模板输入不同的对流系数。固定温度隔热边界对流边界绝热边界输入DX，DY对称边界的部分。不输入与外界的热传递关联数据时，自动认为是绝热边界。仅单元内部进行热传递，不对外散热。 固定温度边界固定温度边界条件意味着该边界处的温度为恒定值，不随时间发生变化，该边界处输入的温度会一直不变。 1. 边界组 : (Prescribed Temp)2. 点击 “前视图”3. 选择地

16、基的外部轮廓4. 点击 “左视图”5. 选择地基的外部轮廓6. 温度 : (20)7. 点击 确认分析 > 水化热分析 > 强制温度操作步骤Procedure3) 发热特性输入混凝土的绝热温度上升曲线。 混凝土的发热特性应根据材料配合比条件不同会有所不同，所以对于实际工程项目要进行绝热温度上升试验，然后在自定义类型里输入实验结果数据。 没有实验数据时可以根据混凝土规范里的说明输入数据。分析 >水化热 > 热源操作步骤1. H/S 荷载组 : (Heat)2. 点击“前视图”3. 选择混凝土的3549个单元4. 点击 “ 热源函数”5. 函数名称 : (Normal)6.

17、 勾选 “设计标准”7. 水泥类型 : (普通硅酸盐水泥)8. 温度 : (20)9. 单位体积水泥用量. : (320)10. 点击 确认, 点击 确认Procedure三输入参数1. 发热特性此操作例题里定义了四个工况，工况中调整了施工阶段、热源函数、对流系数。下边的表格里标记为-的使用与Control Group相同的条件。 | 工况 |序列号工况施工阶段发热函数对流系数1Control group4段分段浇筑普通硅酸盐水泥钢模板2Stage9段分段浇筑-3Heat-低热硅酸盐水泥-4Convection-模板1) 发热函数2: Stage 工况里添加要查看的低热水泥的发热函数。 1.

18、名称 : (Low Heat)2. 函数类型 : 设计标准3.最大绝热温度升高 (K) : (34.2)4. 导温系数(a) : (0.103)5. 点击重画6. 点击确认分析 > 水化热分析 > 热源函数操作步骤Procedure| K和a的标准值 |水泥浇筑温度K(C) = aC+ba(C) = gC+habg( X 10e-3)h低热硅酸盐水泥100.1031.90.180.229200.1002.20.320.225300.0962.60.620.2182) 材料材料种类不同时，其抗压强度的发展特性也有所不同，所以需要定义新的抗压强度发展函数分析 > 时间依存材料 &

19、gt; 抗压强度操作步骤1. 名称 : (C30_Low Heat)2. 规范 : (ACI)3. 混凝土28天抗压强度 : (3.0e7)4. a : (16.2)5. b : (0.82)6. 点击 适用7. 重复上述过程定义“C45_Low Heat”的强度发展函数。 Procedure| 强度发展系数 |水泥的种类abd(28)普通硅酸盐水泥4.50.951.11中热硅酸盐水泥6.20.931.15高早强硅酸盐水泥2.90.971.07粉煤灰水泥B类型(含有20%)7.90.901.40低热硅酸盐水泥16.20.821.4操作步骤Procedure分析 > 材料1. 在列表里选择

20、“2:C30_Normal” 2. 点击 复制3. 列表里选择“4:C30_Normal” 4. 点击修改5. 名称: (C30_Low Heat)6. 抗压强度 : (C30_Low Heat)7. 点击 确认8. 重复上述过程定义低热硅酸盐水泥的强度发展特性“C45_Low Heat” 材料。2. 对流边界结构温度由高温部位向低温部位传递，传递的热量根据对流系数确定。 混凝土内部温度外部温度对流系数混凝土的内部温度由程序自动计算，但“外部温度”、“对流系数”以及“对流边界面”需要用户定义。 此例题对于如下图所示的两个工况设置不同的对流边界条件，然后对分析结果进行比较。1: Control

21、group4: Convection外部接触面钢模板模板1) 定义外部温度函数1. 名称 : (Ambient Temp.)2. 函数类型 : (常量)3. 温度 : (20)4. 点击 : 确认 分析 > 水化热 > 环境温度函数操作步骤Procedure随季节变化的温度变化像大坝、煤气仓等大规模混凝土结构需要较长时间的浇筑施工，此时用户可以自定义随季节变化的外部温度曲线。 2) 定义对流边界函数1. 名称 : (Form_A)2. 函数类型 : (常量)3. 对流系数 : (14)4. 点击 适用 5. 参照对流系数表格定义模板和外部接触面的对流系数 水化热分析 > 单元

22、对流边界操作步骤Procedure| 对流系数 |名称对流系数钢模板 (Form_A)14模板 (Form_B)6外部接触面 (Exposed Surface)133) 对流边界面一般水化热分析里指定对流边界的方法与前面定义固定温度时的方法比较类似。 一般水化热分析时，同时定义固定温度和相应的固定温度的区域。但在FEA的水化热参数化分析里分别定义对流系数和对流系数的区域。即，对于输入相同对流系数的区域指定为一个对流边界面，然后输入多个对流系数分析后进行结果比较。 1. 点击 “前视图”2. 确认为1st Form 。3. 选择混凝土部分。选择过滤里指定为单元-面时只选中一个网格组内的自由面，即

23、只选择网格组的表面。 4. 点击 “取消选择”5. 调整视点取消选择标记的部分。选择没有对称边界及钢材模板的水平面。 6. 取消选择曲线上没有选中的面。7. 点击确认 水化热阶段> 参数化分析设定操作步骤Procedure在同一个工况里可以输入三个对流系数与定义对流边界面的对话框一样，可以对不同的两个区域（第一对流边界、第二对流边界）进行定义。 属于对流边界面区域（未指定固定温度边界或绝热边界的区域），未适用于第一对流边界或第二对流边界时，程序将其默认为与大气接触的边界面。水化热参数化分析时，再同一个工况可以分别定义三个区域（第一对流边界、第二对流边界、大气接触面）的对流系数。3. 施工

24、阶段“定义对流边界面”和“定义施工阶段”参数时，采用不同的方法来定义。 通过在各阶段激活单元来定义施工阶段在进行水化热分析时，一般通过激活钝化三个组（网格 / 边界 / 荷载）来定义施工阶段。 对于不同的参数化分析工况需要多次重复定义施工阶段，工作相当繁琐。且施工阶段数据有变化时，还需要定义不同的网格组、对流边界、边界组等数据。 在进行水化热参数化分析时，为了提高工作效率，程序中添加了仅激活单元就能定义施工阶段的功能。只要激活单元，相应的边界条件将自动被激活。在定义工况时，还设置了考虑自重荷载的选项。一般在进行水化热分析时，很少有钝化结构或钝化边界条件的过程。所以只要考虑在哪个施工阶段激活哪一

25、部分的结构单元即可。激活单元的方法 1 输入高度（便于建模的方法）输入分段浇筑混凝土的高度来定义施工阶段。适用于单向浇筑混凝土，且可按一定高度分段施工的结构。 水化热分析时，常用这种方法来定义施工阶段。 例：桥墩, 基础 激活单元的方法 2 输入网格组（可灵活应用的方法）如下图所示，对于无法按照特定方向、特定高度分段施工的结构，可采用此方法。将需要同时激活的单元定义为同一个网格组，在相应的阶段激活网格组即可。能够定义多种水化热施工阶段。 例： 大坝, 煤气柜Procedure1. 选择“阶段设定” 表单2. 阶段名称 : (A_5m)3. 选择“应用距离” 4. 参考如下对话框输入分段浇筑数据

26、5. 点击 添加6. 重复上述过程输入“B_2m”施工阶段数据水化热分析> 水化热参数化分析操作步骤按特定浇筑高度划分网格使用“应用距离”选项定义施工阶段时，必须按照施工阶段的特定高度划分网格。 初始温度定义混凝土的初始温度。在相应的阶段里输入激活单元时的温度。 | “B_2m”的施工阶段数据 |StageWidthDurationInit. empStep16.51202010, 30, 50, 80, 12021.51202010, 30, 50, 80, 120321202010, 30, 50, 80, 120421202010, 30, 50, 80, 120521202010

27、, 30, 50, 80, 120621202010, 30, 50, 80, 120721202010, 30, 50, 80, 120821202010, 30, 50, 80, 12092.55002010, 30, 50, 80, 120, 180, 250, 350, 500四参数化分析1. 定义工况基本模型的工况（1：Control Group）为：“4阶段分段浇筑”、“普通硅酸盐水泥”、“钢模板”条件。修改上述变量后进行分析，并查看变量对结果的影响。| 定义工况 |序列号工况施工阶段热源函数对流系数1Control gr4段分段浇筑普通硅酸盐钢模板2Stage9段分段浇筑普通硅酸

28、盐钢模板3Heat4段分段浇筑普通 -> 低热钢模板4Convection4段分段浇筑普通硅酸盐模板 (整体)1. 选择“参数” 表单2. 参考表格、对话框以及帮助文件定义分析工况。 水化热分析 > 水化热阶段参数化分析操作步骤Procedure必需输入的参数定义各个变量工况时，还需要输入基本模型里没有定义的数据。即，与“阶段工况”、“热源函数组”、“对流边界“相关的“大气温度”，“系数函数1”，“环境系数函数”等参数。本例题以“4阶段分段浇筑”、“普通硅酸盐水泥”、“钢模板”条件（1：Control Group）的工况为基本模型，分别修改上述变量后进行分析，并查看变量对结果的影响

29、。 | 各变量的输入 |必需输入的参数用红色标记的区域是必需要输入的参数。参数化分析时，定义过程较繁琐的施工阶段和对流边界，不在这里输入。即，这两项不包含在基本模型中，在定义工况时一定要输入。 FEA的发热特性作为荷载来输入，在相应的工况里定义热源荷载组。1: Control Group基本模型的工况采用“4阶段分段浇筑”、“普通硅酸盐水泥”、“钢模板”条件，然后再输入必须输入的参数。2: Stage所有的数据都与基本模型工况（1: Control Group）相同，仅修改施工阶段数据为“2：B_2m”。 3: Heat“施工阶段”、“对流边界”与基本模型相同。 基本模型中输入的热源函数“1：

30、Normal”变更为“2：Low Heat”，水泥的类型变化后，其刚度发展特性也会发生变化，所以还要注意修改材料数据。4: Convection相对于基本模型仅修改“对流边界”。基本模型只对第一对流边界面应用了钢模板，在此工况里对所有对流边界应用“Form_B”对流系数。2. 分析1. 选择“分析” 表单2. 勾选 “有效模量”3. 运行里确认是否勾选了所有的工况4. 点击运行 水化热分析> 水化热阶段参数化分析操作步骤Procedure计算徐变方法混凝土规范介绍了水化热分析考虑徐变时的有效模量（弹性模量折减法）方法。选择“一般”选项后可通过迭代计算徐变效应，同时要在定义材料时定义徐变/

31、收缩函数。在使用“有效模量”的计算方法时，也需要定义徐变/收缩函数，这是为了仅对于与收缩/徐变函数相的构件计算徐变而定义的。例如模型中同时有混凝土和地基时，告知程序哪一部位是混凝土材料。运行对于勾选的多个工况进行分析。 导出选择最终工况后点击，自动生成相应的独立工况的水化热模型。 五查看结果1. 温度分析结束后，在“结果-树形菜单”里同时查看四个工况的分析结果。首先查看基本模型的温度情况。 1. 双击“温度”输入基本模型120小时的温度等值线。 2. 打开“动态.”动态按钮打开后等值线结果将即时更新。3. 点击“输出滑动结果”会出现可以调整步骤的滚动条。 4.用鼠标调整滚动条或者用键盘的键可以

32、查看每个步骤的温度结果。后处理表单 > Control group > Stage 1 Step 5(120) > 节点其他. 操作步骤Procedure2. 应力1. 双击 “LO-SOLID, P1(V)”2. 确认 “ ” 是否亮显 3. 点击“ 输出滑动结果”4. 用鼠标调整滚动条或者用键盘的键查看每个步骤的温度结果。 后处理表单> 传导系数 > Stage 1 Step 4(80) > 3D 单元应力Procedure操作步骤允许应力温度裂缝指数是由发生的抗拉强度和温度应力的比值决定的。在“节点其他”项里查看计算裂缝系数所用的抗拉强度。计算抗拉强度

33、时采用随混凝土的水化反应程度而变化的等价材龄，有可能与抗拉强度发展曲线对应的抗拉强度有些差异。 3. 温度裂缝指数1. 点击“前视图”2. 选择3291单元选择除与地基相邻的混凝土部分 3. 点击 显示 4. 点击 关闭 网格 > 单元/节点 > 显示控制Procedure操作步骤只查看有意义的裂缝指数裂缝指数是根据抗拉强度和温度应力的比值计算的。 因为裂缝指数是节点上的结果，所以在地基和混凝土的边界面上的节点受地基的容许拉应力的影响可能会输出不真实的结果。 地基的抗拉强度是根据混凝土里使用的公式计算，所以输出的结果没什么意义。由于地基的弹性模量比混凝土的弹性模量小很多，计算的抗拉强度也很小。所以最不利的裂缝指数可能集中在地基和