ABAQUS热应力分析实例详解
ABAQUS热应力分析实例详解
热应力分析实例详解学习要点通过实例分析，学习如何进行热应力分析，并掌握ABAQUS/CAE 的以下功能：1）在Material 功能模块中，定义线胀系数；2）在Load 功能模块中，使用预定义场（predefined field）来定义温度场；实例1：带孔平板的热应力分析定义材料属性——Property Property——Material——Edit——steelMechanical——Elastic, 输入弹性模量和泊松比定义材料属性——Property Property——Material——Edit——steelMechanical——Expansion, 输入线胀系数定义边界条件——Load定义边界条件——Load定义边界条件——Load固支边界条件使用预定义场定义初始温度Load——PredefinedField Manager使用预定义场使模型温度升高至120℃网格划分——Mesh结果分析——Visualization小结在ABAQUS中进行热应力分析的基本步骤：⏹定义线胀系数⏹定义初始温度场⏹定义分析步中的温度场实例2：法兰盘感应淬火的残余应力场模拟问题描述：◆表面感应淬火是一种工程中常用的热处理工艺，其原理是使用感应器来对工件的局部进行加热，然后迅速冷却，从而使工件表面产生残余压应力，抵消工作载荷所产生的一部分拉应力。
◆表面感应淬火可显著提高工件弯曲疲劳抗力和扭转疲劳抗力，工件表面产生的马氏体具有良好的耐磨性。
实例2：法兰盘感应淬火的残余应力场模拟 本例中的法兰盘经淬火后，由试验测得法拉盘的内圆角表面残余压应力约为-420MPa。
法拉盘的一端固定，另一端的整个端面受向下的面载荷p=100MPa，法拉盘内孔直径为24mm，材料的弹性模量为210000MPa，泊松比为0.3，线胀系数为1.35e-5/ ℃。
要求：模拟分析感应淬火所产生的残余应力场，并分析此残余应力场在缓和应力集中方面所起的作用。
ABAQUS热应力分析实例详解
热应力分析实例详解学习要点通过实例分析，学习如何进行热应力分析，并掌握ABAQUS/CAE 的以下功能：1）在Material 功能模块中，定义线胀系数；2）在Load 功能模块中，使用预定义场（predefined field）来定义温度场；实例1：带孔平板的热应力分析定义材料属性——Property Property——Material——Edit——steelMechanical——Elastic, 输入弹性模量和泊松比定义材料属性——Property Property——Material——Edit——steelMechanical——Expansion, 输入线胀系数定义边界条件——Load定义边界条件——Load定义边界条件——Load固支边界条件使用预定义场定义初始温度Load——PredefinedField Manager使用预定义场使模型温度升高至120℃网格划分——Mesh结果分析——Visualization小结在ABAQUS中进行热应力分析的基本步骤：⏹定义线胀系数⏹定义初始温度场⏹定义分析步中的温度场实例2：法兰盘感应淬火的残余应力场模拟问题描述：◆表面感应淬火是一种工程中常用的热处理工艺，其原理是使用感应器来对工件的局部进行加热，然后迅速冷却，从而使工件表面产生残余压应力，抵消工作载荷所产生的一部分拉应力。
◆表面感应淬火可显著提高工件弯曲疲劳抗力和扭转疲劳抗力，工件表面产生的马氏体具有良好的耐磨性。
实例2：法兰盘感应淬火的残余应力场模拟 本例中的法兰盘经淬火后，由试验测得法拉盘的内圆角表面残余压应力约为-420MPa。
法拉盘的一端固定，另一端的整个端面受向下的面载荷p=100MPa，法拉盘内孔直径为24mm，材料的弹性模量为210000MPa，泊松比为0.3，线胀系数为1.35e-5/ ℃。
要求：模拟分析感应淬火所产生的残余应力场，并分析此残余应力场在缓和应力集中方面所起的作用。
基于ABAQUS的自由辊温度场及热应力场分析
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!（2）多种零件排样选择part1、part2、part3，其中，part1的长宽不变为100×70，part2的长宽为200×100，part3的长宽为150×90，板材大小仍然为设800×500，X 间距、Y 间距和零件与板材的边距都设为8，排样对话框设置和在AutoCAD 中的排样结果如图4、图5。
5结论通过排样程序可以看出，ObjectARX 编程环境提供的与MFC 相关的用户界面类AdUi 和AcUi 使开发的应用程序能与AutoCAD 风格界面很好地融合在一起，能在同一个界面对AutoCAD 和MFC 对话框进行操作，同时应用Visual.C ++设计较为复杂的对话框。
虽然用ObjectARX 开发应用程序功能强大，但掌握ObjectARX编程方法并不容易，尤其是运行到AutoCAD 界面出现的错误，因为没有具体的错误提示，很难从程序中找到错误，需要经过不停的调试才能找到原因，这给程序设计带来很大困难。
［参考文献］［1］秦洪现，崔惠岚，孙剑，等．Autodesk 系列产品开发培训教程［M ］．北京：化学工业出版社，2008．［2］江思敏，曹默，胡春江.AutoCAD2000开发工具———ObjectARX开发工具与应用实例［M ］．北京：人民邮电出版社，1999．［3］刘蓉梅，姜秀萍，华徐勇，等．ObjectARX 二次开发及应用实例［J ］．机械设计与制造，2002（3）：27-29.（编辑昊天）作者简介：谢友宝（1968-），男，教授，硕士研究生导师，主要研究方向为机电一体化设备研制、数控技术、CAD/CAM 技术、计算机软硬件系统开发等。
收稿日期：2009-06-18图3排样结果图2设置对话框图5排样结果图4对话框设置基于ABAQUS 的自由辊温度场及热应力场分析杨桂芳1，罗会信1，林刚2，代宗岭2（1.武汉科技大学机械自动化学院，武汉430081；2.中冶京诚工程技术有限公司，北京100081）自由辊是连铸机中重要的零件，在结晶器、支撑导向段、扇形段中都有使用。
ABAQUS应用培训-07 热分析实例
7. Mesh：seed part，mesh part, *** Assign element type: heat tr.. 根据傅里叶定律，解析解：
14
2015/7/5
7
ABAQUS传热分析-瞬态热分析
在稳态分析模型的接触上，做出如下修改
1. 2. 3. Material：conductivity 80W/(Cmm),Density,2.7e-9,Special heat 0.8 Step：Heat transfer Transient, time 0.00001
求解传热问题的提法为，在满足边界条件及初始条件的许可温度场中，真 实的温度场使以下泛函 I 取极小值，即
BC S1 , S 2 , S 3 T IC
min
I
1 T 2 T T T x ( ) y ( ) 2 z ( ) 2 2( Q cT )T d 2 x y z t
Fourier传热定律 能量守恒定律
控制方程
6
3
传热过程的基本方程
• Fourier定律：单位时间内通过单位截面积所传递的热量（热流密度， Heat Flux）正比于当地垂直于截面方向上的温度变化率（梯度），热 量传递的方向与温度升高的方向相反。
q k
t x
k为导热系数或热导率[W/(oC▪m)]
1
T 2 T 2 T 2 1 2 x ( x ) y ( y ) z ( z ) d S2 q f TdA 2 S3 hc T T dA
11
稳态传热过程的有限元分析列式
将物体离散为单元体，将单元的温度场表示为节点温度的插值关系，有
16
ABAQUS热应力分析解析实例详解
ABAQUS热应力分析解析实例详解ABAQUS是一种常用的有限元分析软件，可以进行各种不同类型的分析，包括热应力分析。
热应力分析是通过模拟材料受热后发生的变形来评估材料的热稳定性和耐久性。
在这篇文章中，我们将详细介绍ABAQUS热应力分析的步骤和实例。
首先，我们需要创建一个ABAQUS模型。
模型包括几何形状、材料属性和边界条件。
在热应力分析中，我们通常需要定义一个热源，以及材料的热传导、热膨胀和热辐射等属性。
在这个实例中，我们将模拟一个烤箱的加热过程。
模型是一个简单的长方体，材料是钢铁，边界条件是恒定的热流。
下一步是定义材料属性。
我们需要定义钢铁的热传导系数，热膨胀系数和热辐射系数。
这些属性通常可以从材料手册或实验中获得。
我们将使用以下参数：-热传导系数：40W/mK-热膨胀系数：12e-61/°C-热辐射系数：0.8接下来，我们需要定义边界条件。
在这个实例中，我们将模拟一个恒定的热流输入。
我们可以通过选择“控制模拟”菜单中的“载荷”选项来定义边界条件。
在强制边界条件下选择“热流”载荷，然后指定热流的大小和方向。
我们将选择1000W的热流输入。
然后，我们需要定义分析步骤。
在这个实例中，我们将使用一个稳态热分析步骤。
在强制模式下选择“热”分析步骤，然后指定步骤的参数，包括时间步长和总时间。
我们将选择0.1s的时间步长和10s的总时间。
在模拟之前，我们需要定义网格划分。
网格划分是将模型分解为多个小元素的过程，以便于进行数值计算。
ABAQUS中有多种网格划分方法可供选择。
我们可以通过选择“网格”菜单中的“划分”选项来进行网格划分，然后选择适当的网格划分方法和参数。
当所有定义都完成后，我们可以点击“开始模拟”按钮开始进行热应力分析。
ABAQUS将使用已定义的模型、材料属性、边界条件和分析步骤来进行数值计算。
计算结果将显示在ABAQUS的图形界面中。
在热应力分析完成后，我们可以查看结果并进行后处理。
基于ABAQUS的热应力分析
1.1基于ABAQUS的热应力分析1.1.1 温度场数据处理(1)打开INP_Generator.exe，出现如下软件界面：图1.数据处理软件(2)点击“浏览”按钮，选择由FLUENT导出的inp文件所在路径，如下图所示：图2.路径选择(3)点击“生成”按钮，则在inp文件所在路径下自动生成包含多个温度场的ABAQUS输入文件ABAQUSinputfile.inp。
图3.生成包含连续温度场INP文件1.1.2 复材工装模板热应力分析(1)打开ABAQUS，导入inp文件后，打开Tools菜单下“Set - Manager”，如下图所示。
检查是否有名为“PID6”的set，若没有则创建一个名为“PID*”的set，set为模板整体。
（“*”为任意数字或字母）图4.创建SET(2)打开Plug-ins菜单下“CAC Project - Composite Analyse”，弹出如下界面。
在Step1标签中输入用到的材料名称并选择工作路径；在Step2中定义铺层信息，可通过右键删除或添加行；按照Step3和Step4的提示，使用ABAQUS/CAE自身功能完成剩余分析工作。
(a)(b)(c)图5.定义材料及铺层(3)进入Load模块，定义垂直于模板表面平面部分的局部坐标系。
选择“Tools”菜单下“Datum”，Type选择“CSYS”Method选择“3Points”，然后默认点击“Continue”按钮。
依次在模板表面选择坐标原点、X轴上点和XY面上的点，生成局部坐标。
图6.定义模板局部坐标系(4)点击“Create Boundary Condition”按钮，弹出边界条件定义对话框。
Step设为“Initial”，Category选择为“Mechanical”，Types for Selected Step 选择为“Displacement/Rotation”，点击“Continue”，如下图所示：图7.选择约束类型(5)将“Select regions for the boundary condition”选为“by angle”，选中模板下表现所有结点（按住Shift键可多选），点击鼠标中键，弹出如下边界条件编辑对话框，给模板施加U3和UR3约束，CSYS选择为模板局部坐标系。
abaqus热应力分析实例_200105
Abaqus热应力分析实例1 说明：本例通过简单的杆状零件，介绍abaqus热分析的基本步骤。
利用abaqus/CAE分析图1所示的杆状零件，四面加热条件下（随时间升温T=20+5t）的温度场，并以该温度为初始条件，分析零部件受力状况。
图1为杆状零件截面的图2传热分析2.1创建part进入part模块，点击创建部件，name输入bar，模型所在空间选择3维，类型选择可变性，shape选择Solid，Type选择Extrusion，Approximate size 输入200，设置如下图，点击Continue，进入二维截面创建，分别输入（25,25）、（-25,-25）两两点，完成草图绘制，Depth（长度）输入500，完成部件的创建，如下图所示。
2.2 创建材料和截面切换到property模块，Density输入7.74e-09，Conductivity（传热率）、Specific Heat （比热）与温度有关，输入如下：2.3点击，弹出Create Section对话框，name输入Section-1，Categeory选择Solid，type选择Homogeneous，点击continue，弹出Edit Section，选择刚创建的材料Steel。
2.4赋予属性点击，选择部件，中键确定，完成材料赋予。
2.5创建分析步创建一个Heat Transfer（热传递）分析步，点击Continue，basic工具栏设置，选择Transient(瞬态分析)，time period设置为100，切换到incrementation,设置如下图。
2.6 热传递与热辐射设置在杆四周面加载一个随时间变化的的温度T=20+5t，切换到interation模块，创建温度曲线，Tools》Amplitude》create，name输入Amp-1，Type选择Tabular，列表设置如下左图。
点击，分析步选择step-1，选择surface file condition，点击continue，film coefficient 设置为0.4，Sink temperature 为1，Sink amplitude 选择上述创建的温度曲线。
ABAQUS热分析
h *FILM 二维情况下施加在单边上,三维情况下施加在单元面上 *FILM ELSET, F3., 450, 2.3E-3 温度 *SFILM 二维情况下施加在面上
h
页码，5/14
/view/d16dff05cc17552707220897?pn=2&vw=all&ssid=&from=&bd_page_type=1&uid=806880DB5C59C1... 2013-10-9
/view/d16dff05cc17552707220897?pn=2&vw=all&ssid=&from=&bd_page_type=1&uid=806880DB5C59C1... 2013-10-9-Boltzmann常数和绝对零度 *PHYSICAL CONSTANTS, ABSOLUTE ZERO = -273.16 STEFAN BOLTZMANN = 5.6697E-8 边界条件与载荷 4. 向环境的辐射 辐射率 emissivity 是衡量一个表面有多接近理想黑体的指标
页码，6/14 第3/6页
NT13 NT12 NT11
n
-- 壳单元表面下方的温度自由度为11(输出变量为NT11) -- 在正表面的温度自由度为 10+n, n 为壳截面上使用截面点的数量 -- 在单层(均质)壳 中,截面点在厚度上均匀分布,默认为5个点 -- 每层壳必须是奇数个截面点,这是由 ABAQUS/standard 在厚度上使用分段抛物线 型插值方法决 定的.6页
的. 热传导分析中,这对共轭变量是 温度 --- 热率(单位时间的能量流) 默认情况下温度是未知的,热率是已知的 -- 已知的热率 = 0, 相当于绝 热边界条件; -- 没有外部的能量流进或流出节点. ABAQUS 中的几种热边界条件和热载荷 1. 在某些节点上预设温度, *BOUNDARY, 自由度1 1 2. 在某些点上或者某些表面上或者体积内预设热率 q *CFLUX, *DFLUX, *DSFLUX 3. 在某些点上或者某些表面上的边界层(薄膜)条件 *C FILM, *FILM 和 *SFILM 4. 在某些点上或者某些表面上的辐射条件 *CRADIATE, *RADIATE, 和 *SRADIATE 5. 自然边界条件(默认)
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
设置分析步——step
1）第一个分析步HighTemper-Noload: 令淬硬层区域温度升 高至120℃，其余区域温度仍保持20℃，不施加外载荷； 2）第二个分析步HighTemper-WithLoad:保持上述温度场不 变（相应的残余应力也不会变），施加外载荷； 3）第三个分析步LowTemper-WithLoad：令整个工件的温 度都变为20℃（即去掉残余应力），保持外载荷不变，从 而得到没有残余应力时的应力场，用来与第二个分析步的 结果相比较。
18
建模要点说明
上述方法的优点是比较简便，不必进行复杂的传热分析和 热弹塑性分析，并且通用性强，可用于模拟各种不同工艺 所产生的残余应力场，但其缺点是模拟精度不高，通过选 择Thigh只能保证工件局部区域的压应力值较准确。 一种改进的方法就是为淬硬层的不同区域设定不同的温度 值Thigh，从而得到与试验结果更加接近的残余应力场。 本实例中，为简单起见，只为整个淬硬层设定单一的温度 值Thigh =120℃。
热应力分析实例详解
热应力分析的主要问题
设定材料的线胀系数 设定模型的初始温度场 可以直接给出温度值 也可以读入传热分析的结果文件 修改在分析步中的温度场 可以直接给出温度值 也可以读入传热分析的结果文件
2
学习要点
通过实例分析，学习如何进行热应力分析，并掌 握 ABAQUS/CAE 的以下功能： 1）在 Material 功能模块中，定义线胀系数； 2）在 Load 功能模块中，使用预定义场 （predefined field）来定义温度场；
要求：模拟分析感应淬火所产生的残余
应力场，并分析此残余应力场在缓和应 力集中方面所起的作用。
17
建模要点说明
使用ABAQUS可以模拟感应淬火的完整过程，即通过分析工件 与感应器之间以及工件和冷却液之间的传热过程来确定工件的 温度场，从而得到相应的塑性应变场和冷却后的残余应力场。 这一模拟过程较复杂，下面介绍一种模拟残余应力的简化方 法。 设置整个模型的初始温度为20℃，在分析步令淬硬层的温度升 高至某一温度值Thigh(例如120℃)，其余区域温度仍保持20℃。 这种温度差异会使高温区产生压应力，相当于所要模拟的残余 压应力。经几次试算即可找到合适的Thigh，使法兰盘内圆角表 面的压应力与试验结果大致吻合。施加工作载荷时，仍保持上 述温度场不变，就可以模拟在残余应力作用下的应力场了。
13
结果分析——Visualization
14
小结
在ABAQUS中进行热应力分析的基本步骤：
定义线胀系数 定义初始温度场 定义分析步中的温度场
15
实例2：法兰盘感应淬火的残余应力场模拟 问题描述：
表面感应淬火是一种工程中常用的热处理工艺， 其原理是使用感应器产生残余压应 力，抵消工作载荷所产生的一部分拉应力。 表面感应淬火可显著提高工件弯曲疲劳抗力和扭 转疲劳抗力，工件表面产生的马氏体具有良好的 耐磨性。
3
实例1：带孔平板的热应力分析
问题描述：
去掉原有的拉伸载荷，在平板顶部添加固支边界条件。 整个平板的初始温度为20℃，当温度升高为120℃时，平 板会发生热膨胀，而平板顶部的固支约束会限制模型的变 形，模型的应力场发生相应的改变。材料的线胀系数为 100 1.35×10-5/℃。
要求：分析模型在120℃下的应力场。
Mechanical——Expansion, 输入线胀系数
7
定义边界条件——Load
8
定义边界条件——Load
9
定义边界条件——Load
固支边界条件
10
使用预定义场定义初始温度
Load——Predefined Field Manager
11
使用预定义场使模型温度升高至120℃
12
网格划分——Mesh
565号节点
与残余压应力的试验结果420MPa基本吻合！
33
2）分析残余压应力在缓解应力集中方面的作用
第一个分析步
σminP= 276MPa
560号节点
34
2）分析残余压应力在缓解应力集中方面的作用
第二个分析步：存在残余应力！
σmaxP= 412MPa
560号节点
35
2）分析残余压应力在缓解应力集中方面的作用
19
几何建模
CAD平面图
20
几何建模
Part
12mm
旋转轴
21
几何建模
22
几何建模
23
定义材料属性——Property
E=210 000 =0.3 Expansion=1.03e-5
24
网格划分——Mesh
淬硬层以外区域: Element type: Hex-dominated Technique: Sweep Algorithm: Advancing front 淬硬层: Element type: Hex Technique: Sweep Algorithm: Medial axis (选中Minimize the mesh transition)
y 0 x
50 50
R5 100
4
几何建模——Part
100
50
y 0 x
R5 100 50
5
定义材料属性——Property
Property——Material——Edit——steel
Mechanical——Elastic, 输入弹性模量和泊松比
6
定义材料属性——Property
Property——Material——Edit——steel
26
边界条件——Load
一端固定
27
边界条件——Load
28
施加载荷——Load
p=100MPa
29
定义温度场——Load
30
定义温度场——Load
20 120 120 20
淬硬层
31
后处理
1）查看残余应力的模拟结果
32
1）查看残余应力的模拟结果
第一个分析步
σminiP= 416.17MPa
16
实例2：法兰盘感应淬火的残余应力场模拟
本例中的法兰盘经淬火后，由试验测得法拉盘的内圆角表 面残余压应力约为 -420MPa。 法拉盘的一端固定，另一端的整个端面受向下的面载荷 p=100MPa，法拉盘内孔直径为24mm，材料 的弹性模量为210000MPa，泊松比为0.3， 线胀系数为1.35e-5/ ℃。