midasfea水化热参数化分析.docx
《midasfea水化热参数化分析.docx》由会员分享,可在线阅读,更多相关《midasfea水化热参数化分析.docx(15页珍藏版)》请在冰豆网上搜索。
midasfea水化热参数化分析
midasfea-水化热参数化分析
一.概要
1.水化热分析
浇筑混凝土时,水泥在水化过程中产生大量热量会使混凝土的温度升高。
虽然随时间的推移混凝土的温度会慢慢冷却,但结构各个位置的温度下降速度不均匀,结构不同位置将发生相对温差,此温差会使混凝土发生温度应力。
温度裂缝发生类型
混凝土浇筑初期,因内部温度升高将发生膨胀,但混凝土表面的温度下降较快,相对应变较小,从而使混凝土表面产生拉应力。
混凝土内部不同的温度分布引起的不同的体积变化而导致的应力称为内部约束应力。
此类拉应力裂缝主要发生在构件尺寸比较大的结构。
|内部约束产生的裂缝(放热时)||外部约束产生的裂缝(冷却时)|
混凝土在高温状态下温度下降会发生收缩,但受到与其接触的已浇筑混凝土或者地基等的约束而产生的拉力,像这样变形受外部边界约束的状态称为外部约束。
此类应力主要发生在像墙这样约束度比较大的结构中。
利用温度裂缝指数预测温度裂缝
韩国混凝土规范中使用温度裂缝指数(抗拉强度与发生的温度应力之比)i值预测是否发生裂缝。
一般采用下面的值。
因此通过查看温度分布可以看出输入数据是否有误,如果温度分布没有问题可说明输出的应力结果也是正确的。
水化热分析必须进行反复计算
大体积混凝土的温度裂缝可以利用温输入混凝土的散热特性及浇筑条件等度裂缝指数(CrackRatio,Icr)来验
算。
温度裂缝指数要满足结构的重要
混凝土的温度性、功能、环境条件等因素的要求。
温度裂缝指数受水泥的类型、浇筑温度、养生方法等多因素的影响,所以应力需要对多种条件进行反复分析以找出最佳的浇筑方法。
No裂缝指数
Ye参数化分析功能
END为比较多种条件的分析结果需要建立
多个模型进行分析,分析结束后需要整理大量的分析结果、还要进行结果保存、对比等工作。
通过FEA的水化热参数化分析功能,可以实现一个模型多种条件分析。
可以大大减少单纯繁琐的反复分析过程,从而提高工作效率。
参数化分析的使用方法
首先建立一个基本模型,在基本模型里使用替换变量的方式定义分析工况。
下图是把材料作为变量条件的示例,“CaeI”为将混凝土C24变更为C30的工况,“CaeII”为将混凝土C35变更为C40的工况。
|参数化分析的构成|
裂缝指数(i)=
混凝土抗拉强度发生的温度应力
参数化分析里可以考虑的变量
在水化热参数化分析的功能里可以调整的变量有五个,较常用的调整方法具体如下。
施工阶段:
降低浇筑高度缩小各阶段的温度差。
浇筑间距过小的话很难
达到分段浇筑的效果,但如果太大分界面会产生较大的温差。
.对流边界:
对流系数较低时,热量不容易对外流失,可以减少内外温差。
材料:
使用弹性模量大的材料时,抗拉强度也较大,可增大裂缝指数。
发热特性:
是变量中最为敏感的因素,定义水化过程中发生的热量。
是否考虑自重:
使混凝土产生压应力的荷载,在一定程度上可以减少拉应力,但效果不明显。
|温度裂缝指数与裂缝发生几率|防止裂缝发生时:
1.5以上限制裂缝发生时:
1.2~1.5限制有害裂缝发生时:
0.7~1.2
FEA程序的水化热分析
水化热分析主要分为热传导分析和热应力分析。
.
热传导分析主要计算水泥的水化过程中发热、传导、对流等引起的随时间变化的节点温度。
将得到的节点温度作为荷载加载后,计算随时间变化的应力称为热应力分析。
二.建立基本模型
1.结构分析所需的数据
水化热参数化分析如前面图形所示。
首先建立基本模型,通过在基本模型里定义“Ued”和“New”的对应关系来定义分析工况。
|进行参数化分析时输入变量的示意图|
但是利用这种方法很难输入多个变量,虽然输入热源函数或对流系数函数很简单,但反复定义施工阶段和对流边界面的过程较繁琐。
为避免这种繁琐的定义过程,在定义水化热分析变量时,先不定义此两项。
而是先定义阶段工况,即按不同的施工方案定义不同的施工阶段工况,然后再对不同的阶段工况定义各自的材料、对流系数等。
|进行参数化分析所输入的数据|
1)建立模型
此操作例题主要介绍“水化热参数化分析”的方法,仅对于相关变量的输入、各种分析条件的定义、查看分析结果部分进行详细说明。
导入附件里的“HYD_Pier_Meh.feb”模型文件。
对称模型
施工阶段水化热分析模型一般单元数量较多,所以分析所需的时间也较长,而且还要进行多条件分析工况的分析,所以需要更多的分析时间。
如果模型属于对称模型,可以只建立模型以减少分析时间。
这样不但可以减少分析时间也有利于查看模型中心部位的分析结果。
2)特性
操作步骤分析>分析控制...Procedure1.单位体系:
确认指定为(N,m,J)
操作步骤分析>时间依存性材料>徐变/收缩Procedure1.名称:
(C30)2.规范:
(中国(JTGD62-2004)3.混凝土28天材龄抗压强度:
(3.0e7)4.点击[适用]5.名称:
(C45)6.混凝土28天材龄抗压强度:
(4.5e7)7.点击[确认]
时间依存材料特性
在FEA里可以考虑徐变/收缩特性和混凝土的抗压强度变化。
徐变计算方法
可以选择用户定义或者使用规范的计算方法。
弹性模量折减方法是假设发生徐变,然后折减混凝土弹性模量的简易计算方法,一般的水化热分析里均使用这种方法。
因弹性模量折减方法只适用于水化热分析,为了避免在一般的施工阶段分析中误用,在水化热分析控制里单独定义。
采用弹性模量折减方法时,为了要指定计算徐变的单元(通过材料对话框中的相应选项将徐变函数与材料连接起来,然后通过给单元赋予材料将徐变函数与单元连接起来),需要任意定义一个徐变计算方法,但在这里定义的徐变计算方法并不参与水化热分析的徐变的计算。
操作步骤Procedure分析>时间依存性材料>抗压强度1.名称:
(C30_Normal)2.类型:
设计规范3.规范:
ACI4.混凝土28天抗压强度(f28):
(3.0e7)5.a:
(4.5)6.b:
(0.95)7.点击[适用]8.重复上面的步骤定义‘C45_Normal’的强度发展函数。
计算裂缝指数用的抗拉强度
类型里选择“用户定义”后,用户可以通过表格自定义随时间变化的弹性模量、抗压强度、抗拉强度。
选择“设计规范”时,利用规范提供的公式计算弹性模量和抗拉强度。
利用抗拉强度计算温度裂缝指数。
混凝土抗压强度系数
与水泥的种类有关,一般硅酸盐水泥a为4.5,b为0.95。
序列号12名称C30_NormalC45_Normal过程分析>材料序列号123名称SoilC30_NormalC45_Normal弹性模量1e73.0e103.35e10泊松比0.20.20.2膨胀系数重量密度1e-51e-51e-5260002451724517徐变/收缩-C30C45抗压强度-C30_NormalC45_Normal抗压强度3.0e74.5e7a4.54.5b0.950.95Procedure操作步骤Procedure分析>特性1.点击下拉菜单选择3D。
2.参考下表定义三个特性。
序列号123
名称SoilFoundationPier材料1:
Soil2:
C30_Normal3:
C45_Normal21.名称:
(Soil)2.弹性模量:
(1e7)3.泊松比:
(0.2)4.膨胀系数.:
(1e-5)5.重量密度:
(26000)6.点击[适用]7.参考下表输入‘C30_Normal’和‘C45_Normal’特性。
操作步骤分析>材料Procedure1.点击“前视图”2.特性:
(3D)3.选择“1:
Soil”4.选择“3DElement(R)5.选择地基1710个单元6.点击[适用]7.重复上述过程定义基础和桥墩的特性。
对称面边界条件
对于Y-Z平面上的所有节点约束D某。
1347
需要查看输入的对称边界条件时,可将网格显示为特征边线后查看。
如下图所示。
5
地基-1:
Soil基础-2:
Foundation桥墩-3:
Pier
3)边界条件
操作步骤分析>边界条件>约束Procedure1.点击“前视图”2.边界组:
Sym3.选择Y-Z对称面上的360个节点4.DOF:
(T1)5.点击[适用]6.点击“左视图”7.选择某-Z对称面上的705个节点8.自由度:
(T2)9.点击[确认]
操作步骤Procedure分析>边界>约束1.边界组:
Support2.点击“前视图”3.选择地基的外部轮廓4.点击“左视图”5.选择地基的外部轮廓6.自由度:
(T1,T2,T3)实体单元没有旋转自由度所以只选择上述三个自由度。
7.点击[确认]
固结边界条件
水化热分析模型里建立地基时,一般将地基下部的边界条件设为完全固结。
混凝土产生的热量将充分地传递给地基,后续不再传递温度,也不存在温差,所以也不会发生相对位移。
为了让混凝土产生的热量充分地传递给地基,需要建立足够大的地基模型。
5
5
2.热传导分析所需的数据
1)材料的热特性数据
定义水化过程产生的热传递的特性。
模型单元内的热传导由比热和热传导率确定,单元外的散热由下一节定义的对流系数来确定。
操作步骤Procedure分析>材料
|一般岩体和混凝土的热工系数|1.定义的材料列表里选择1:
Soil2.点击[修改]3.点击[热工参数…]4.传导率:
(3.45)5.比热:
(784)6.点击[确认]7.参考下表输入混凝土的热特性值。
岩体混凝土
|热工特性|热传导率(W/m2·°C)1.7~5.22.6~2.8比热(kJ/kg·°C)0.71~0.881.05~1.26序列号123材料名称SoilC30_NormalC45_Normal热传导率3.452.72.7比热78411761176
2)固定温度
为了定义实体单元表面与外界的热传递关系,需要定义相应的对流边界条件。
固定温度边界、对流边界、绝热边界如下图所示。
固定温度边界
固定温度边界条件意味着该边界处的温度为恒定
值,不随时间发生变化,该边界处输入的温度会一
直不变。
绝热边界
输入D某,DY对称边界的部分。
不输入与外界的热
传递关联数据时,自动认为是绝热边界。
仅单元内
部进行热传递,不对外散热。
对流边界(热传导边界)
固定温度
通过对流系数定义各时刻与外部交换的热量。
对于
隔热边界
不同的模板输入不同的对流系数。
对流边界
3)发热特性
输入混凝土的绝热温度上升曲线。
混凝土的发热特性应根据材料配合比条件不同会有所不同,所以对于实际工程项目要进行绝热温度上升试验,然后在自定义类型里输入实验结果数据。
没有实验数据时可以根据混凝土规范里的说明输入数据。
操作步骤分析>水化热>热源Procedure1.H/S荷载组:
(Heat)2.点击“前视图”3.选择混凝土的3549个单元
4.点击“热源函数”5.函数名称:
(Normal)6.勾选“设计标准”7.水泥类型:
(普通硅酸盐水泥)8.温度:
(20)9.单位体积水泥用量.:
(320)10.点击[确认],点击[确认]
操作步骤Procedure分析>水化热分析>强制温度1.边界组:
(PrecribedTemp)2.点击“前视图”3.选择地基的外部轮廓4.点击“左视图”5.选择地基的外部轮廓6.温度:
(20)7.点击[确认]
三.输入参数
1.发热特性
此操作例题里定义了四个工况,工况中调整了施工阶段、热源函数、对流系数。
下边的表格里标记为‘-’的使用与ControlGroup相同的条件。
|工况|序列号1234工况ControlgroupStageHeatConvection施工阶段4段分段浇筑9段分段浇筑--发热函数普通硅酸盐水泥-低热硅酸盐水泥-对流系数钢模板--模板
1)发热函数
‘2:
Stage’工况里添加要查看的低热水泥的发热函数。
操作步骤Procedure分析>水化热分析>热源函数1.名称:
(LowHeat)2.函数类型:
设计标准3.最大绝热温度升高(K):
(34.2)4.导温系数(a):
(0.103)5.点击[重画]6.点击[确认]
|K和a的标准值|水泥低热硅酸盐水泥浇筑温度102030K(C)=aC+ba0.1030.1000.096b1.92.22.6a(C)=gC+hg(某10e-3)0.180.320.62h0.2290.2250.218
操作步骤Procedure分析>时间依存材料>抗压强度1.名称:
(C30_LowHeat)2.规范:
(ACI)3.混凝土28天抗压强度:
(3.0e7)4.a:
(16.2)5.b:
(0.82)6.点击[适用]7.重复上述过程定义“C45_LowHeat”的强度发展函数。
|强度发展系数|水泥的种类普通硅酸盐水泥中热硅酸盐水泥高早强硅酸盐水泥粉煤灰水泥B类型(含有20%)低热硅酸盐水泥a4.56.22.97.916.2b0.950.930.970.900.82d(28)1.111.151.071.401.4
2)材料
材料种类不同时,其抗压强度的发展特性也有所不同,所以需要定义新的抗压强度发展函数
操作步骤Procedure分析>材料1.在列表里选择“2:
C30_Normal”2.点击[复制]3.列表里选择“4:
C30_Normal”4.点击[修改]5.名称:
(C30_LowHeat)6.抗压强度:
(C30_LowHeat)7.点击[确认]8.重复上述过程定义低热硅酸盐水泥的强度发展特性“C45_LowHeat”材料。
2.对流边界
结构温度由高温部位向低温部位传递,传递的热量根据对流系数确定。
混凝土内部温度外部温度对流系数
混凝土的内部温度由程序自动计算,但“外部温度”、“对流系数”以及“对流边界面”需要用户定义。
此例题对于如下图所示的两个工况设置不同的对流边界条件,然后对分析结果进行比较。
1:
Controlgroup4:
Convection外部接触面钢模板模板
2)定义对流边界函数
操作步骤水化热分析>单元对流边界Procedure1.名称:
(Form_A)2.函数类型:
(常量)3.对流系数:
(14)4.点击[适用]5.参照对流系数表格定义模板和外部接触面的对流系数|对流系数|名称钢模板(Form_A)模板(Form_B)外部接触面(E某poedSurface)对流系数146131)定义外部温度函数
操作步骤Procedure分析>水化热>环境温度函数1.名称:
(AmbientTemp.)2.函数类型:
(常量)3.温度:
(20)4.点击:
[确认]
3)对流边界面
一般水化热分析里指定对流边界的方法与前面定义固定温度时的方法比较类似。
一般水化热分析时,同时定义固定温度和相应的固定温度的区域。
但在FEA的水化热参数化分析里分别定义对流系数和对流系数的区域。
即,对于输入相同对流系数的区域指定为一个对流边界面,然后输入多个对流系数分析后进行结果比较。
操作步骤Procedure水化热阶段>参数化分析设定1.点击“前视图”2.确认为1tForm3.选择混凝土部分。
选择过滤里指定为单元-面时只选中一个网格组内的自由面,即只选择网格组的表面。
4.点击“取消选择”5.调整视点取消选择标记的部分。
选择没有对称边界及钢材模板的水平面。
6.取消选择曲线上没有选中的面。
7.点击[确认]
随季节变化的温度变化
像大坝、煤气仓等大规模混凝土结构需要较长时间的浇筑施工,此时用户可以自定义随季节变化的外部温度曲线。
在同一个工况里可以输入三个对流系数
与定义对流边界面的对话框一样,可以对不同的两个区域(第一对流边界、第二对流边界)进行定义。
属于对流边界面区域(未指定固定温度边界或绝热边界的区域),未适用于第一对流边界或第二对流边界时,程序将其默认为与大气接触的边界面。
水化热参数化分析时,再同一个工况可以分别定义三个区域(第一对流边界、第二对流边界、大气接触面)的对流系数。
3.施工阶段
“定义对流边界面”和“定义施工阶段”参数时,采用不同的方法来定义。
通过在各阶段激活单元来定义施工阶段
在进行水化热分析时,一般通过激活钝化三个组(网格/边界/荷载)来定义施工阶段。
对于不同的参数化分析工况需要多次重复定义施工阶段,工作相当繁琐。
且施工阶段数据有变化时,还需要定义不同的网格组、对流边界、边界组等数据。
在进行水化热参数化分析时,为了提高工作效率,程序中添加了仅激活单元就能定义施工阶段的功能。
只要激活单元,相应的边界条件将自动被激活。
在定义工况时,还设置了考虑自重荷载的选项。
一般在进行水化热分析时,很少有钝化结构或钝化边界条件的过程。
所以只要考虑在哪个施工阶段激活哪一部分的结构单元即可。
激活单元的方法1–输入高度(便于建模的方法)输入分段浇筑混凝土的高度来定义施工阶段。
适用于单向浇筑混凝土,且可按一定高度分段施工的结构。
水化热分析时,常用这种方法来定义施工阶段。
例:
桥墩,基础
激活单元的方法2–输入网格组(可灵活应用的方法)
如下图所示,对于无法按照特定方向、特定高度分段施工的结构,可采用此方法。
将需要同时激活的单元定义为同一个网格组,在相应的阶段激活网格组即可。
能够定义多种水化热施工阶段。
例:
大坝,煤气柜
操作步骤水化热分析>水化热参数化分析Procedure1.选择“阶段设定”表单2.阶段名称:
(A_5m)3.选择“应用距离”4.参考如下对话框输入分段浇筑数据5.点击[添加]6.重复上述过程输入“B_2m”施工阶段数据
按特定浇筑高度划分网格
使用“应用距离”选项定义施工阶段时,必须按照施工阶段的特定高度划分网格。
初始温度
定义混凝土的初始温度。
在相应的阶段里输入激活单元时的温度。
|“B_2m”的施工阶段数据|StageWidthDuration1234567896.51.52222222.5120120120120120120120120500Init.emp202020202020202020Step10,30,50,80,12010,30,50,80,12010,30,50,80,12010,30,50,80,12010,30,50,80,12010,30,50,80,12010,30,50,80,12010,30,50,80,12010,30,50,80,120,180,250,350,500四.参数化分析
1.定义工况
基本模型的工况(1:
ControlGroup)为:
“4阶段分段浇筑”、“普通硅酸盐水泥”、“钢模板”条件。
修改上述变量后进行分析,并查看变量对结果的影响。
|定义工况|序列号1234工况ControlgrStageHeatConvection施工阶段4段分段浇筑9段分段浇筑4段分段浇筑4段分段浇筑热源函数普通硅酸盐普通硅酸盐普通->低热普通硅酸盐对流系数钢模板钢模板钢模板模板(整体)操作步骤水化热分析>水化热阶段参数化分析Procedure1.选择“参数”表单2.参考表格、对话框以及帮助文件定义分析工况。
必需输入的参数
定义各个变量工况时,还需要输入基本模型里没有定义的数据。
即,与“阶段工况”、“热源函数组”、“对流边界“相关的“大气温度”,“系数函数1”,“环境系数函数”等参数。
本例题以“4阶段分段浇筑”、“普通硅酸盐水泥”、“钢模板”条件(1:
ControlGroup)的工况为基本模型,分别修改上述变量后进行分析,并查看变量对结果的影响。
|各变量的输入|
2.分析
操作步骤Procedure水化热分析>水化热阶段参数化分析1.选择“分析”表单2.勾选“有效模量”3.运行里确认是否勾选了所有的工况4.点击[运行]
必需输入的参数
用红色标记的区域是必需要输入的参数。
参数化分析时,定义过程较繁琐的施工阶段和对流边界,不在这里输入。
即,这两项不包含在基本模型中,在定义工况时一定要输入。
FEA的发热特性作为荷载来输入,在相应的工况里定义热源荷载组。
1:
ControlGroup
基本模型的工况采用“4阶段分段浇筑”、“普通硅酸盐水泥”、“钢模板”条件,然后再输入必须输入的参数。
2:
Stage
所有的数据都与基本模型工况(1:
ControlGroup)相同,仅修改施工阶段数据为“2:
B_2m”。
3:
Heat
“施工阶段”、“对流边界”与基本模型相同。
基本模型中输入的热源函数“1:
Normal”变更为“2:
LowHeat”,水泥的类型变化后,其刚度发展特性也会发生变化,所以还要注意修改材料数据。
4:
Convection
相对于基本模型仅修改“对流边界”。
基本模型只对第一对流边界面应用了钢模板,在此工况里对所有对流边界应用“Form_B”对流系数。
计算徐变方法
混凝土规范介绍了水化热分析考虑徐变时的有效模量(弹性模量折减法)方法。
选择“一般”选项后可通过迭代计算徐变效应,同时要在定义材料时定义徐变/收缩函数。
在使用“有效模量”的计算方法时,也需要定义徐变/收缩函数,这是为了仅对于与收缩/徐变函数相的构件计算徐变而定义的。
例如模型中同时有混凝土和地基时,告知程序哪一部位是混凝土材料。
运行
对于勾选的多个工况进行分析。
导出
选择最终工况后点击,自动生成相应的独立工况的水化热模型。
五.查看结果
2.应力
允许应力
温度裂缝指数是由发生的抗拉强度和温度应力的比值决定的。
在“节点其他”项里查看计算裂缝系数所用的抗拉强度。
计算抗拉强度时采用随混凝土的水化反应程度而变化的等价材龄,有可能与抗拉强度发展曲线对应的抗拉强度有些差异。
操作步骤Procedure后处理表单>传导系数>Stage1Step4(80)>3D单元应力1.双击“LO-SOLID,P1(V)”2.确认“”是否亮显3.点击“输出滑动结果”4.用鼠标调整滚动条或者用键盘的↓键查看每个步骤的温度结果。
1.温度
分析结束后,在“结果-树形菜单”里同时查看四个工况的分析结果。
首先查看基本模型的温度情况。
操作步骤Procedure后处理表单>Controlgroup>Stage1Step5(120)>节点其他.1.双击“温度”输入基本模型120小时的温度等值线。
2.打开“动态.”动态按钮打开后等值线结果将即时更新。
3.点击“输出滑动结果”会出现可以调整步骤的滚动条。
4.用鼠标调整滚动条或者用键盘的↓键可以查看每个步骤的温度结果。
3.温度裂缝指数
操作步骤网格>单元/节点>显示控制Procedure1.点击“前视图”2.选择3291单元选择除与地基相邻的混凝土部分3.点击[显示]4.点击[关闭]
只查看有意义的裂缝指数
裂缝指数是根据抗拉强度和温度应力的比值计算的。
因为裂缝指数是节点上的结果,所以在地基和混凝土的边界面上的节点受地基的容许拉应力的影响可能会输出不真实的结果。
地基的抗拉强度是根据混凝土里使用的公式计算,所以输出的结果没什么意义。
由于地基的弹性模量比混凝土的弹性模量小很多,计算的抗拉强度也很小。
所以最不利的裂缝指数可能集中在地基和混凝土的边界面上,这样的结果是不真实的。
为了避免输出上述错误结果,程序将自动除去边界面上的单元,仅激活混凝土部分的单元后输出裂缝指数结果。
若想重新查看隐藏的地基单元时,点击即可。
操作步骤特性>等值线Procedure1.反转颜色:
(True)裂缝指数越小越不利所以最小的值标记为红色。
2.最大最小值显示:
(True)虽然裂缝指数输出0~20范围内的值,但是“2”以上就已经是非常安全的值,所以范例变更为0~2后查看。
3.最大值:
(2)4.点击[适用]
操作步骤Procedure后处理表单>Heat>HYDRATIONMIN>节点其他1.双击“裂缝指数(法向)”
HYDRATIONMIN/HYDRATIONMA某
输出相应工况的所有步骤中的最大或最小值。
例如HYDRAIONMIN输出的裂缝指数是相应节点的在所有步骤中的最小值。
水化热分析的目的就是要查看发生最不利裂缝指数的节点位置。
操作步骤Procedure后处