基于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选择为模板局部坐标系。

图8.定义约束

(6)打开Predefined Field Manager,检查温度场Field-1，若为模板整体温度

场则删除它，否则保留。创建一个初始温度场：Step设为“Initial”，Category 选择为“Other”，Type for Selected Step选择为“Temperature”，点击“Continue”按钮，选中模板整体，输入常温293K。并将其在Initial以后各步聚中设为“Reset to initial”。

图9.检查温度场

(7)进入Job模块，创建作业，提交计算，完成后打开工作目录下对应的odb

文件查看分析结果。

创建并提交作业

(8)打开Plug-ins菜单下“CAC Project – Post Process”,弹出后处理对话框，

各选项作用如下图(b)所示：

(a)

整体变形

变形后云图

变形前后云图

应力云图

应变云图

位移云图

(b)

图10.后处理

(9)打开Plug-ins菜单下“CAC Project –Strength Check”,弹出强度校核对话

框。

图11.强度校核

(10)点击“Select…”选择分析后的结果文件(Odb文件)，输入insurance name，

默认为PART-1-1（应与Odb文件中insurance name保持一致）。最后点击

“OK”，将对模板进行最大应力强度准则校核，并在ABAQUS下面的信

息栏给出校核结果：（Strength Qualified or Strength Unqualified）。

1.1.3 制件精度分析

a)对FLUENT导出的INP文件作处理，去除其中温度点，如下图。

图12.去除温度点

b)使用File-Import-Model功能将去除温度点的INP导入ABAQUS。按5.3.2的方法，设置型面以模板材料从常温升致固化温度。温度场设置如下：点击

Load模块定义升温过程温度场，“Step”选为“Initial”；“Category”选为“Other”；“Type for Selected Step”选为“Temperature”，点击“Continue”，如下图所示：

图13.建立温度场

初始步(Initial)温度设为常温25度，如下图所示：

图14.常温设置

打开“Predefined Field Manager”，双击“Step-1”下面的“Propagated”，如下图：

图15.温度场控制器

将“Step-1”中温度场更改为固化温度197度，设置如下图：

图16.改为固化温度

c)上一步分析完成后，重新打开ABAQUS，导入固化温度下的型面。使用

File-Import-Part功能，如下图所示：

图17.导入高温型面

按5.3.2的方法设置型面以制件材料从固化温度降为常温，分析得到制件的处形。

d)按5.3.2的方法设置型面以制件材料从常温到常温，形状不变的一个过

程，分析得到制件理论外形。

e)制件精度比较：打开Plug-ins菜单下“CAC Project –Dimensional

Accuracy”,弹出制件精度分析对话框，如下图：

图18.制件精度分析

上面选择制件理论外形的结果文件，下面选择制件变形的结果文件，点击OK，将得到包含制件理论外形和变形后所有节点坐标值的文本文件和对应节点最大位移值，并以此判断制件的精度。

ABAQUS顺序热力耦合分析实例

ABAQUS顺序热力耦合分析实例此实例中需要确定一个冷却栅管的温度场分布。温度场的求解采用稳态热分析，在此之后还将进行热应力分析来求出冷却栅管在温度作用下产生的位移和应力分布。由于冷却栅管比较长，并且是轴对称结构，根据上述特点，可以简化有限元分析模型。此实例中使用国际单位制。 1、part中创建轴对称可变形壳体，大致尺寸为1，通过creat line创建一个封闭曲线（0.127,0） （0.304,0）（0.304,0.006）（0.152,0.006）（0.152,0.031）（0.127,0.031）（0.127,0） 使用creat Fillet功能对模型倒角处设置0.005的倒圆角。倒角后，模型并未改变，需要在模型树中，part下的Features右键，Regenerate，最终模型如下图所示。 2、在材料模块中定义密度7800，弹性模量1.93E11，泊松比0.3。所不同的是，热分析还需 要指定热传导系数以及比热。在Thermal里输入参数，热铲刀系数25.96，比热451。 3、创建截面属性以及装备部件，和普通的静力分析设置一样。 4、Step有所不同，分析类型仍为通用分析步，下面要更改为Heat Transfer。在Edit Step窗 口中，使用默认的瞬态分析（Transient），时长设置为3s。切换到Incrementatin进行相应的设置，如下图。

5、Load模块中，设置左边温度为100度，右边及上边温度为20度。Creat BC，类型选择 Other>Temperature。在纯粹的热传导分析方程中，没有位移项，因此不会发生刚体位移，这里也就不需要设置位移边界条件。 6、接下来划分网格，种子尺寸给0.005，单元类型需要在单元族中选择专门用来热分析的 Heat Transfer，查看下面确保使用的单元为DCAX4。使用结构化的全四边形网格划分方法。 7、到此，热分析的设置已经完成，可以提交计算，完成后，查看变量NT11即为节点温度。

abaqus热分析

Abaqus： ABAQUS 是一套功能强大的工程模拟的有限元软件，其解决问题的范围从相对简单的线性分析到许多复杂的非线性问题。ABAQUS 包括一个丰富的、可模拟任意几何形状的单元库。并拥有各种类型的材料模型库，可以模拟典型工程材料的性能，其中包括金属、橡胶、高分子材料、复合材料、钢筋混凝土、可压缩超弹性泡沫材料以及土壤和岩石等地质材料，作为通用的模拟工具。 模块： ABAQUS/?ba:kjus/有两个主求解器模块— ABAQUS/Standard 和ABAQUS/Explicit。ABAQUS 还包含一个全面支持求解器的图形用户界面，即人机交互前后处理模块—ABAQUS/CAE。 ABAQUS 对某些特殊问题还提供了专用模块来加以解决。ABAQUS 被广泛地认为是功能最强的有限元软件，可以分析复杂的固体力学结构力学系统，特别是能够驾驭非常庞大复杂的问题和模拟高度非线性问题。ABAQUS 不但可以做单一零件的力学和多物理场的分析，同时还可以做系统级的分析和研究。ABAQUS 的系统级分析的特点相对于其他的分析软件来说是独一无二的。由于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”按钮，弹出边界条件定义对话框。

abaqus有限元建模例子

问题一:工字梁弯曲 1.1问题描述: 在<<材料力学实验>>中，弯曲实验測定了工字梁弯曲应变大小及其分布，以验证弯曲正应力公式。在这里，採用ABAQUS/CAE建立试验件的有限元模型，ABAQUS/Standard模块进行分析求解，得到应力、应变分布，对比其与理论公式计算值及实验測量值的差別。 弯曲实验的相关数据： 材料：铝合金E=70GPa泊松比0.3 实验装置结构简图如图所示： 结构尺寸测量值：H=50（+/-0.5mm） h=46(+/-0.5mm) B=40（+/-0.5mm） b=2(+/-0.02mm) a=300(+/-1mm) F1=30N Fmax=300N N ? F100 = 1.2ABAQUS有限元建模及分析 一对象: 工字型截面铝合金梁 梁的结构简图如图1所示，結构尺寸、载荷、約束根据1.1设定，L取1600mm，两端各伸出100mm。 二用ABAQUS/CAE建立实验件的有限元模型,效果图如下： 边界条件简化： 左侧固定铰支座简化为下表面左参考点处的约束U1=U2=U3=0

右侧活动铰支座简化为下表面右参考点处的约束U1=U2=UR3=0 几何模型

有限元模型 三ABAQUS有限元分析結果 ①应力云图（Z方向正应力分量）：施加载荷前 F=300N

②应变（Z方向分量）： 中间竖直平面的厚度方向应变分布图： F=100N F=200N

F=300N 由上图可以看出应变沿着厚度方向呈线性比例趋势变化，与实验测得的应变值变化趋势相同。中性轴处应变均接近零值，应变与距离中性轴位移基本为正比关系。 1.3分析结果: 中间竖直截面上下边缘轴向应力数值对比：*10^-6MPa 距中性轴距ABAQUS模拟实验测量值平均理论值 1/2H-96.182*70000-97*70000-6.9165=-70000*98.807 -1/2H95.789*7000092*70000 6.9165

热应力分析

ABAQUS可以求解以下类型的传热问题： 1.非耦合传热分析：温度场不受应力应变场或电场的影响。应用ABAQUS/Standard可以求 解导热问题、强制对流、边界辐射和空腔辐射问题，其分析类型可以是瞬态或稳态的，也可以是线性或非线性的。 2.顺序耦合热应力分析：应力应变场受温度场的影响，但温度场不受应力应变场的影响。 此类问题用ABAQUS/Standard求解的步骤为：先求解温度场，然后以其作为已知条件，进行热应力分析，得到应力应变场。分析传热问题和热应力分析可以使用不一样的网格，abaqus会自动进行差值处理（此类问题称为热应力分析）。 3.完全耦合热应力分析：温度场和应力应变场之间有着强烈的相互作用。 4.绝热分析：在此类分析中，力学变形会产生热，而且整个过程中时间极短，不发生热扩 散。 5.热电耦合分析：用来求解电流产生的温度场。 7.1热应力分析中的主要问题 设定线胀系数、模型的初始温度场，并可以修改分析步中的温度场。 7.2带孔平板的热应力分析 学习： 在LOAD功能模块中，使用预定义场（predefined field）来定义温度场。 在此模块中可以直接指定温度场或读入分析结果文件中的温度场，可以指定并精确读入某个分析步中某个增量步的温度场 7.3法兰盘感应淬火的残余应力模拟 学习： 使用热应力来模拟残余应力；在LOAD功能模块中，为模型的各个区域定义不同的温度场 表面感应淬火：常用的热处理工艺，使用感应器对工件表面进行局部加热，然后迅速冷却，在工件内部产生残余压应力。它可以提高工件的弯曲疲劳抗力和扭转疲劳抗力，工件表面的

马氏体具有良好的耐磨性。 Abaqus可以完整的模拟淬火的全过程，即通过分析工件和感应器之间以及工件与冷却液之间的热场过程来确定工件的温度场，从而得到相应的塑性应变场和冷却后的残余应变场。 比较简单的模拟方法：先设定整个模型的初始温度场，在分析过程中令淬硬层区域的温度升高至某个温度值，其余区域的温度保持不变。经过几次试算，找到合适的淬硬层温度值，使得法兰盘内圆角处的表面压应力与实验结果吻合。施加工作载荷，保持上述温度场不变，就可以模拟在残余应力作用下的应力场。 优点：通用性强，可以模拟不同工艺所产生的残余应力场 缺点：精确度不高 改进方法：参淬硬层的不同区域设定不同的温度值

abaqus热残余应力分析实例

利用Abaqus的Moldflow接口进行翘曲分析和残余应力分析 Abaqus关键特征和优势 1.力学性质、有限元网格以及残余应力数据都能从Moldflow很简便地传递到Abaqus 2.包含了成型工艺残余应力的Abaqus分析使得注塑模具产品的仿真更加精确 分析方法 对一个注塑模具产品的翘曲和应力分析的过程来说，一开始是利用Moldflow对注塑成型过程进行仿真。Moldflow的分析结果包括材料性质的描述以及固化零件中的残余应力分布。Abaqus的Moldflow接口此时用来将这些数据转换成Abaqus可以应用的格式。特别强调的是，接口产生的文件包含了塑料的网格信息、残余应力结果以及材料的性质。这些数据会在接下来的Abaqus分析中用来进行翘曲和残余应力影响的建模。椅子和手机外壳塑模的离散化模型如图1所示。对于这两个模型，Moldflow分析在模型厚度上分了21层并使用了壳体网格元素。翘曲的仿真运用Abaqus/Standard的静态分析功能分析完成。 图1：椅子和手机外壳模型的网格 结果和讨论 运用Abaqus/Standard进行翘曲分析后，椅子模型和手机外壳模型的变形如图2及图3所示。

图2：椅子模型的翘曲位移[米]分布云图 图3：手机外壳模型的翘曲位移[米]分布云图 由Abaqus/Standard翘曲分析所得到的椅子模型和手机外壳模型的Mises应力分布云图如图4及图5所示。很明显可以看出，由于翘曲引起了变形，原来零件中所储存的Mises 应力大小降低了。

图4：椅子模型的Mises应力[帕]分布分布—翘曲前[左]和翘曲后[后] 图5：手机外壳模型的Mises应力[帕]分布—翘曲前[左]和翘曲后[后] 结论 Abaqus为进行细致的结构分析提供了强大的能力。Moldflow为注塑模具产品提供了运算残余应力和材料性质的能力。Abaqus的Moldflow接口通过提供Moldflow分析结果向Abaqus分析过程传送的方法，使得更加精确、更加高效的设计过程得以实现。

abaqus热分析

ABAQUS 热分析常用概念介绍 热传递通过热传导、对流和热辐射三种方式实现。 热传导是热量重系统的一部分传到另一部分或由一个系统传到另一系统的现象。模型中有两种方式实现，共点网格和接触对。热阻系数=空气热传导率/空气间隙。 对流是液体或气体中较热部分和较冷部分之间通过循环流动使温度趋于均匀的过程，对流是液体和气体中热传递的特有方式，气体的对流现象比液体明显，对流分为自然对流和强迫对流。 辐射是物体因自身的温度而具有向外发射能量的本领。

热传导分析中的基本物理量 ·温度Temperature 单位℃ ·热能Heat energy 单位J ·热率Heat rate power 单位J/t or W ·热流量Heat Flux = Powerper unit area单位J/t/L2 材料参数介绍

1. CONDUCTIVITY 热传导率用于度量热量在材料中流动的难易程度： 单位：W/m/℃ 在热传递分析中，传导率为必需的材料属性。 2. SPECIFIC HEAT 比热用于度量热能在材料中存储的难易程度： 单位：J/Kg/℃ 3. Emissivity辐射率是衡量一个表面有多接近理想黑体的指标 （0~1）。 边界条件与载荷介绍 边界条件与载荷，在热传导分析中，每个自由度的共轭变量为温度-热率（单位时间的能量流）。 1. 预设的温度*BOUNDARY，包括两种，恒定温度和变化的温度，温度的共轭反作用是热率（热能进入一个已经预设温度值的节点的流通率）。

2. 预设热流量（热率），*CFLUX,节点的集中热流量；*DFLUX，施加在面或体上的分布热流量,*DSFLUX施加在面上的分布热流量。 3. 预设边界层条件最常见的一种边界条件为一个自由表明被紧临的流体加热或降温，关键字*CFILM,施加在节点上；*FILM二维中施加在单元边上，三维中施加在单元面上；*SFILM二维中施加在单元边上；边界层系数h是ABAQUS的一个输入参数，量纲：J/L2*T*θ。 4. 辐射条件*CRADIATE,施加在节点上；*RADIATE,施加在单元上；*SRADIATE，施加在面上，定义辐射边界条件，需要定义 Stefan-Boltzmann常数和绝对零度。 5. 自然边界条件（默认）

abaqus热分析

ABAQUS作为最常用的求解器，具有强大的仿真功能和热分析求解能力。ABAQUS 不仅可以用于热传导分析，还可以用于温度场和其他领域的耦合分析 1.传热 2.耦合温度位移 3.耦合热电分析 4.耦合热电结构分析 ①导热分析 对于热分析，准确定义材料和元素尤为重要。ABAQUS为此分析提供了一个单位（dc3d8）。在材料定义方面，ABAQUS提供电导率，比热，密度等。此外，对于某些特殊效果，可以使用以下材料属性：内部发热（仅ABAQUS /标准）和用户定义的本构响应（ABAQUS）/标准）。此外，ABAQUS提供了电导率，比热，密度，弹性模量（Ex），泊松比等的定义。 根据热分析的类型，ABAQUS提供稳态分析，瞬态分析和非线性分析。

ABAQUS提供各种形式的温度指定，热通量指定，对流边界条件设置，对周围环境的辐射的定义以及自然边界条件和初始条件的设置。 对于热分析中的接触问题，ABAQUS提供了一种热“接触”方式，它通过界面传热，热相互作用，间隙传热和间隙辐射来模拟接触位置的传热。 ②热耦合分析 热应力耦合分析是热分析必不可少的部分。ABAQUS提供了两种方法进行热应力耦合分析，顺序耦合分析和完全耦合分析。顺序耦合分析是先进行热传导分析，然后使用热传导分析的结果进行热应力分析。假定温度会导致热应力，但是应力对温度没有响应。完整的耦合分析考虑了两者之间的相互反应。在热耦合分析中，ABAQUS为不同类型的热耦合分析提供了特殊的耦合元素，并且热传导分析前面的材料，载荷和边界的定义适用于耦合分析。 ③ABAUQS胎面制动的热分析 使用ABAQUS耦合温度位移分析步骤执行完全耦合热分析。