热应力分析在工程中十分常见,之前跟大家一起学习了用ansys进行热应力分析的方法,今天和大家一起学习如何采用Abaqus进行热应力分析。
壹
问题描述
如下图所示的圆筒体带球形封头结构,其几何参数如下:
圆筒及封头外径:1000mm
圆筒及封头厚度:100mm
圆筒长度:1000mm
材料参数:
名称 | 弹性模量 | 泊松比 | 导热系数 | 线膨胀系数 |
单位 | MPa | / | W/(m·℃) | 1/℃ |
值 | 200000 | 0.3 | 35 | 1e-5 |
结构内部温度边界对流换热边界:换热系数 1000 W/(m^2·℃),环境温度为400℃。结构外表面温度边界为对流换热边界:换热系数 50 W/(m^2·℃),环境温度为25℃。
采用Abaqus进行热应力分析。由于该问题中温度只有边界条件,没有温度场分布,因此,在进行热应力分析之前,需先进行稳态温度场分析。
贰
稳态温度场分析
01
定义材料
分别定义材料的导热系数、弹性模量、泊松比、热膨胀系数,如下图所示(由于使用的是mm制单位,相应的数值需要注意单位协调)。
02
创建截面属性、赋予截面属性、模型装配
过程略(可参考文章《Abaqus入门(5)-弹塑性分析实例》)
03
创建稳态传热分析步(Step-1)
04
定义对流换热边界
在Interaction中分别定义两个对流换热边界(过程参考《Abaqus入门(6)-稳态传热分析》),注意对流换热系数的单位制换算。
05
划分网格
需要在Asign Element Type中指定单元类型为Heat Transfer。
06
创建分析并提交计算
创建一个名为HeatTransfer的工作,并提交分析。
07
查看温度场分布
温度场分布
叁
热应力分析
01
创建结构分析步
在Step Manager中将Step-1(Heat Transfer)抑制掉。
然后创建一个Step-2(Static,General)
02
创建边界条件
分别在三个方向的对称面上创建对称约束。
03
施加温度场
通过Create Predefined Field,选择Temperature。
Distribution选择From results or output database file。
点击File name后的打开文件目录按钮,选择之前分析生成的温度场结果文件HeatTransfer.odb。
04
修改单元类型
在Asign Element Type中指定单元类型为3D Stress。
05
创建热应力分析工作
创建名为Structure的工作,提交计算。
06
查看热应力、位移
应力云图
位移云图