ANSYS非稳态热分析及实例详解

ANSYS流体及热场分析

第 7 章 非稳态热分析及实例详解

本章向读者介绍非稳态热分析的基本知识，主要包括非稳态热分析的应用、非稳态热分析单元、非稳态热分析的基本步骤。

本章要点

? 非稳态导热的基本概念 ? 非稳态热分析的应用 ? 非稳态热分析单元 ?

分析的基本步骤

本章案例

? 钢球非稳态传热过程分析

? 不同材料金属块水中冷却的非稳态传热过程分析 ? 高温铜导线冷却过程分析

1

第7章 非稳态热分析

7.1 非稳态热分析概述

物体的温度随时间而变化的导热过程称为非稳态导热。根据物体温度随着时间的推移而变化的特性可以区分为两类非稳态导热：物体的温度随时间的推移逐渐趋于恒定的值以及物体的温度随时间而作周期性的变化。无论在自然界还是工程实际问题中，绝大多数传热过程都是非稳态的。许多工程实际问题需要确定物体内部的温度场随时间的变化，或确定其内部温度达到某一限定值所需要的时间。例如：在机器启动、停机及变动工况时，急剧的温度变化会使部件因热应力而破坏，因此需要确定物体内部的瞬时温度场；钢制工件的热处理是一个典型的非稳态导热过程，掌握工件中温度变化的速率是控制工件热处理质量的重要因素。再例如，金属在加热炉内加热时，需要确定它在加热炉内停留的时间，以保证达到规定的中心温度。可见，非稳态热分析是有相当大的应用价值的。ANSYS 11.0及其相关的下属产品均支持非稳态的热分析。非稳态热分析确定了温度以及其它随时间变化的热参数。

7.1.1 非稳态热分析特性

瞬态热分析用于计算一个系统的随时间变化的温度场及其它热参数。在工程上一般用瞬态热分析计算温度场，并将之作为热载荷进行应力分析。

瞬态热分析的基本步骤与稳态热分析类似。主要的区别是瞬态热分析中的载荷是随时间变化的。为了表达随时间变化的载荷，首先必须将载荷-时间曲线分为载荷步。对于每一个载荷步，必须定义载荷值及时间值，同时必须选择载荷步为渐变或阶越。

7.1.2 非稳态热分析的控制方程

热储存项的计入将稳态系统变为非稳态系统，计入热储存项的控制方程的矩阵形式如下：

????K??T???Q? ?C??T?为热储存项。 其中，?C?T在非稳态分析时，载荷是和时间有关的函数，因此控制方程可表示如下：

??????K??T???Q?t?? ?C??T若分析为分线性，则各参数除了和时间有关外，还和温度有关。非线性的控制方程可表示如

下：

???C?T????T????K?T????T???Q?T,t??

7.1.3 时间积分与时间步长

1、时间积分

从求解方法上来看，稳态分析和非稳态分析之间的差别就是时间积分。利用ANSYS 11.0分析

问题时，只要在后续载荷步中将时间积分效果打开，稳态分析即转变为非稳态分析；同样，只要在后续载荷步中将时间积分关闭，非稳态分析也可转变为稳态分析。 2、时间步长

两次求解之间的时间称为时间步，一般来说，时间步越小，计算结果越精确。确定时间步长的方法有两种：

（1）指定裕度较大的初始时间步长，然后使用自动时间步长增加时间步。 2

ANSYS流体及热场分析

（2）大致估计初始时间步长。

在非稳态热分析中估计初始时间步长，可以使用Biot数和Fourier数。 Biot数是不考虑尺寸的热阻对流和传导比例因子，其定义为：

Bi?h?x K式中：?x——名义单元宽度； h——平均表面换热系数； K——平均导热系数。

Fourier数是不考虑尺寸的时间（?t/t），其定义为：

Fo?式中：?——平均密度；

K?t

?c(?x)2 c——比热容；

如果Bi?1，可将Fourier数设为常数并求解?t来预测时间步长：

?c(?x)2(?x)2??t???? ??

???c式中：?——热耗散。

如果Bi?1，时间步长可应用Fourier数和Biot数的乘积预测：

???t??h?x??h?t?Fo?Bi???????? 2????c(?x)??K???c?x?求解?t得到：

?t??其中，0.1???0.5

?c?xh

时间步长的预测精度随单元宽度的取值、平均的方法、比例因子?的变化而变化。

7.1.4 数值求解过程

当前温度矢量?Tn?假设为已知，可以是初始温度或由前面的求解得到的。定义下一个时间点的温度矢量为：

?????t?T?? ?Tn?1???Tn??(1??)?t?Tnn?1其中?称为欧拉参数，默认为1，下一个时间点的温度为：

????K??T???Q? ?C??Tn?1n?1由上面两式可得：

1?????1??1C?KT?Q?CT?Tn?? ??????n?1???????n???1???t????t? 3

第7章 非稳态热分析

??K???Tn?1???Q?

其中??1 K??C???K??????????t??Q???C?1?????1T?Tn????Q? ???n???t1??欧拉参数?的数值在0.5~1之间。在这个范围内，时间积分算法是不明显而且是不稳定的。因

此，ANSYS 11.0总是忽略时间积分步的幅值来计算。但是，这样的计算结果并不总是准确的。下面是选择积分参数的一些建议：

当?=0.5时，时间积分方法采用“Crank-Nicolson”技术。本设置对于绝大多数热瞬态问题都是精确有效的。

当?=1时，时间积分方法采用“Backward Euler”技术。这是缺省的和最稳定的设置，因为它消除了可能带来严重非线性或高阶单元的非正常振动。 本技术一般需要相对Crank-Nicolson较小的时间积分步得到精确的结果。

7.2 非稳态热分析单元

非稳态热分析和稳态热分析使用的分析单元相同，具体请读者参见本书第6章。

7.3 非稳态热分析基本步骤

非稳态热分析的基本步骤主要包括：建模、加载求解和后处理。下面分别对这三个基本步骤进行具体的阐述。

7.3.1 建立有限元模型

就这一步骤而言，并没有稳态和非稳态之分，可参照稳态分析的建模方法进行。因此，在这里不在赘述。

7.3.2 加载求解

1、定义分析类型

如果第一次进行分析或重新进行分析，操作步骤如下： Command: ANTYPE,TRANSIENT,NEW

GUI：Main Menu>Solution>Analysis Type>New Analysis>Transient 如果接着上次的分析继续进行（例如增加其它载荷），操作步骤如下： Command: ANTYPE,TRANSIENT,REST GUI：Main Menu>Solution>Analysis Type>Restart 2、获得非稳态热分析的初始条件

（1）定义均匀温度场

如果已知模型的起始温度是均匀的，可设定所有节点初始温度，操作步骤如下： Command: TUNIF

4

ANSYS流体及热场分析

GUI: Main Menu>Solution>Loads>Settings>Uniform Temp

如果不在对话框中输入数据，则默认为参考温度，参考温度的值默认为零，但可通过如下方法设定参考温度：

Command: TREF

GUI: Main Menu>Solution>Loads>Settings>Reference Temp

注意：设定均匀的初始温度，与如下的设定节点的温度（自由度）不同，设定节点温度的操作步骤如下：

Command: D

GUI: Main Menu>Solution>Loads>Apply>Thermal>Temperature>On Nodes 初始均匀温度仅对分析的第一个子步有效；而设定节点温度将保持贯穿整个瞬态分析过程，除非通过下列方法删除此约束：

Command: DDELE

GUI: Main Menu>Solution>Loads>Delete>Thermal>Temperature>On Nodes （2）设定非均匀的初始温度

在瞬态热分析中，节点温度可以设定为不同的值，操作步骤如下： Command: C

GUI: Main Menu>Solution>Loads>Apply>Initial Condit'n>Define 如果初始温度场是不均匀的且又是未知的，就必须首先作稳态热分析确定初始条件，步骤如下： ·设定载荷（如已知的温度、热对流等） ·将时间积分设置为OFF： Command: TIMINT, OFF

GUI: Main Menu>Preprocessor>Loads>Load Step Opts>Time/Frequenc>Time Integration ·设定一个只有一个子步的，时间很小的载荷步（例如0.001）： Command: TIME

GUI: Main Menu>Preprocessor>Loads>Load Step Opts>Time/Frequenc>Time and Substps ·写入载荷步文件： Command: LSWRITE

GUI: Main Menu>Preprocessor>Loads>Write LS File 或先求解：

Command: SOLVE

GUI: Main Menu>Solution>Solve>Current LS 3、设定载荷步选项

（1）普通选项

·设置时间和时间步步，操作如下： Command: TIME

GUI: Main Menu>Solution>Load Step Opts>Time/Frequenc>Time-Time Step ·设置每个载荷步的载荷子步数，或时间增量 Command: NSUBST or DELTIM

GUI: Main Menu>Solution>Load Step Opts>Time/Frequenc>Time and Substps （2）非线性选项 ·设置迭代次数：（每个子步默认的次数为25，这对大多数非线性热分析已经足够） Command: NEQIT

GUI: Main Menu>Solution>Load step opts>Nonlinear>Equilibrium Iter ·自动时间步长：（本选项为ON时，在求解过程中将自动调整时间步长）

5

第7章 非稳态热分析

Command: AUTOTS

GUI: Main Menu>Solution>Load Step Opts>Time/Frequenc>Time and Substps ·时间积分效果：（如果将此选项设定为OFF，将进行稳态热分析） Command: TIMINT

GUI: Main Menu>Solution>Load Step Opts>Time/Frequenc>Time Integration （3）输出选项 ·控制打印输出：（本选项可将任何结果数据输出到*.out 文件中） Command: OUTPR

GUI: Main Menu>Solution>Load Step Opts>Output Ctrls>Solu Printout ·控制结果文件：（控制*.rth的内容） Command: OUTRES

GUI: Main Menu>Solution>Load Step Opts>Output Ctrls>DB/Results File （4）存盘求解

7.3.3 后处理

ANSYS提供两种后处理方式：

·POST1，可以对整个模型在某一载荷步（时间点）的结果进行后处理； Command: POST1

GUI: Main Menu>General Postproc.

·POST26，可以对模型中特定点在所有载荷步（整个瞬态过程）的结果进行后处理。 Command: POST26

GUI: Main Menu>TimeHist Postproc 1、用POST1进行后处理

·进入POST1后，可以读出某一时间点的结果： Command: SET

GUI: Main Menu>General Postproc>Read Results>By Time/Freq

如果设定的时间点不在任何一个子步的时间点上，ANSYS会进行线性插值。

·此外还可以读出某一载荷步的结果：

GUI: Main Menu>General Postproc>Read Results>By Load Step

然后就可以采用与稳态热分析类似的方法，对结果进行彩色云图显示、矢量图显示、打印列表等后处理。

2、用POST26进行后处理

·首先要定义变量：

Command: NSOL or ESOL or RFORCE

GUI: Main Menu>TimeHist Postproc>Define Variables ·然后就可以绘制这些变量随时间变化的曲线： Command: PLVAR

GUI: Main Menu>TimeHist Postproc>Graph Variables

6

ANSYS流体及热场分析

或列表输出： Command: PRVAR

GUI: Main Menu>TimeHist Postproc>List Variables

此外，POST26还提供许多其它功能，如对变量进行数学操作等，请参阅《ANSYS Basic Analysis Procedures Guide》

实例7－1：钢球非稳态传热过程分析

一个直径为12cm，温度为1000℃的钢制小球突然被放入了盛满了水的、完全绝热的横截面直径和高度均为60cm的圆柱体水槽内（钢球放在水槽的正中央），水的温度为18℃，材料参数表如表7.1所示。试求10分钟后钢球与水的温度场分布。

表7.1 物性参数表

热性能 导热系数 密度 比热 单位制 W/m℃ Kg/m3 J/kg℃ 铁 70 7800 448 水 0.60 1000 4200 该问题是典型的瞬态传热问题，研究对象为钢球和水。由于对称性，在求解过程中取钢球和水中心纵截面的1/4建立几何模型，如图7-1所示。本例选取PLANE55轴对称单元进行求解。

图7-1 几何模型

——附带光盘“Ch7\\实例7－1_start” ——附带光盘“Ch7\\实例7－1_end” ——附带光盘“AVI\\Ch7\\7－1.avi”

1、定义工作文件名

选择Utility Menu>File>Change Jobname，弹出Change Jobname对话框。在对话框中将工作名改为example7-1，单击OK关闭该对话框。选择Main Menu>Preferences，弹出Preferences for GUI Filtering对话框，选中Thermal复选框，然后单击Ok按钮关闭该对话框。

7

第7章 非稳态热分析

2、定义单元类型

选择Main Menu>Preprocessor>Element Type>Add/Edit/Delete，弹出Element Type对话框，单击Add按钮，弹出Library of Element Types对话框。在Library of Element Types对话框的两个列表框中分别选择Thermal Solid、Quad 4node 55选项，如图7-2所示。单击OK按钮关闭该对话框。

图7-2 单元类型列表对话框

单击Element Type对话框中的Options按钮，弹出PLANE55 element type options对话框，在Element behavior K3下拉列表框中选择Axisymmetric选项，其余选项均采用默认设置，如图7-3所示，单击OK关闭该对话框。单击Element Type对话框中的Close按钮关闭该对话框。

图7-3 PLANE55单元属性设置对话框

3、定义材料性能参数

1、选择Main Menu>Preprocessor>Material Props>Material Model，弹出Define Material Model Behavior对话框，如图7-4所示。

2、在Material Models Available列表中依次双击Thermal、Conductivity、Isotropic选项，弹出Conductivity for Material Number1对话框，在文本框中输入钢的导热系数70，如图7-5所示，点击OK关闭对话框。

3、双击Define Material Model Behavior对话框上的Specific Heat按钮，弹出Specific Heat for Material Number1对话框，在文本框中输入钢的比热448，如图7-6所示，单击OK按钮关闭该对话框。

4、双击Define Material Model Behavior对话框上的Density按钮，弹出Density for Material Number1对话框，在文本框中输入钢的密度7800，如图7-7所示，单击OK按钮关闭该对话框。

8

ANSYS流体及热场分析

图7-4 定义材料属性对话框 图7-5 定义材料导热系数对话框

图7-6 定义材料比热对话框 图7-7 定义材料密度对话框

5、定义水的材料属性：在Define Material Model Behavior对话框中选择Material>New Model，弹出Define Material ID对话框，在文本框中输入材料参数号2，如图7-8所示单击OK按钮关闭该对话框。

图7-8 定义材料编号对话框

6、在Material Models Available列表中依次双击Thermal、Conductivity、Isotropic选项，弹出Conductivity for Material Number2对话框，在文本框中输入水的导热系数0.6，点击OK关闭对话框。

7、双击Define Material Model Behavior对话框上的Specific Heat按钮，弹出Specific Heat for Material Number2对话框，在文本框中输入水的比热4200，单击OK按钮关闭该对话框。

8、双击Define Material Model Behavior对话框上的Density按钮，弹出Density for Material Number2对话框，在文本框中输入水的密度1000，单击OK按钮关闭该对话框。关闭Define Material Model Behavior对话框。 4、创建几何模型

1、选择Main Menu>Preprocessor>Modeling>Create>Areas>Rectangle>By Dimensions，弹出Create Rectangle by Dimensions对话框，如图7-9所示，在X1、X2文本框中分别输入0、0.3，在Y1、Y2文本框中分别输入0、0.3，然后单击OK按钮确认设置。

2、选择Main Menu>Preprocessor>Modeling>Create>Areas>Circle>By Dimensions，弹出Circular Areas by Dimensions对话框。参照图7-10对其进行设置，然后单击OK确认设置。

9

第7章 非稳态热分析

图7-9 创建矩形面对话框 图7-10 创建圆面对话框

3、选择Main Menu>Preprocessor>Modeling>Operate>Booleans>Overlap>Areas，弹出Overlap Areas对话框。单击Pick All按钮选取所有的面。

4、选择Main Menu>Preprocessor>Numbering Ctrls>Compress Numbers，弹出Compress Numbers对话框。在Label Item to be compressed下拉列表框中选择All选项，单击OK按钮确认设置。

5、选择Utility Menu>PlotCtrls>Numbering，弹出Plot Numbering Controls对话框，选择LINE选项，使其状态从Off变为On，其余选项均采用默认设置，单击OK按钮关闭该对话框。完成上述操作后，生成的几何模型如图7-11所示。

图7-11 生成的几何模型

5、划分有限元网格

1、选择Utility Menu>Plot>Lines。

2、选择Main Menu>Preprocessor>Meshing>Size Cntrls>ManualSize>Lines>Picked Lines，弹出Element Size on菜单，在文本框中输入4，5，单击OK按钮，弹出Element Size on Picked Lines对话框，在NDIV文本框中输入划分的网格单元个数30，在SPACE文本框中输入0.1，如图7-12所示，单击OK按钮关闭对话框。

3、选择Main Menu>Preprocessor>Meshing>Size Cntrls>ManualSize>Lines>Picked Lines，弹出Element Size on菜单，在文本框中输入6，7，单击OK按钮，弹出Element Size on Picked Lines对话框，在NDIV文本框中输入划分的网格单元个数32，在SPACE文本框中输入0.1，单击OK按钮关闭对话框。

4、选择Main Menu>Preprocessor>Meshing>Size Cntrls>ManualSize>Lines>Picked Lines，弹出Element Size on菜单，在文本框中输入3，单击OK按钮，弹出Element Size on Picked Lines对话框，10

第7章 非稳态热分析

图7-26 水内部温度场分布等值线云图

8、选择Utility Menu>PlotCtrls>Style>Graphs>Modify Axes，弹出Axes Modification for Graph Plots对话框，参照图7-27进行设置，设置完成后单击OK按钮关闭该对话框。

图7-27 坐标设置对话框

9、选择Utility Menu>PlotCtrls>Style>Graphs>Modify Curves，弹出Curve Modification for Graph Plots对话框，在[/GTHK]下拉列表框中选择Triple选项，单击OK按钮关闭该对话框。选择Utility Menu>Select>Everything。

10、选择Main Menu>TimeHist Postpro>Define Variables，弹出Define Time-History Variables对话框，单击Add按钮，弹出Add Time-History Variable对话框，选中Nodal DOF Result单选按钮，16

ANSYS流体及热场分析

如图7-28所示，单击OK按钮，弹出Define Nodal Data菜单，在文本框中输入1，单击OK按钮，弹出Define Nodal Data对话框，参照图7-29对其进行设置，单击OK按钮关闭该对话框。单击Define Time-History Varables对话框中的Close按钮关闭该对话框。

图7-28 添加时间历程变量对话框 图7-29 定义节点数据对话框

11、选择Main Menu>TimeHist Postpro>Graph Variables，弹出Graph Time-History Varable对话框，在NAVR1文本框中输入2，如图7-30所示，单击OK按钮，ANSYS显示窗口将显示将显示球心温度随时间的变化关系曲线图，如图7-31所示。

图7-30 显示时间历程对话框 图7-31 球心温度随时间变化曲线图

实例7－2：不同材料金属块水中冷却的非稳态传热过程分析

一个30公斤重、温度为1000℃的铜块,以及一个20公斤重、温度为800℃的铁块，突然放入温度为20℃、盛满了300升水的、完全绝热的水箱中，如图7-32所示。试求：一个小时后，铜块与铁块的温度分布(假设忽略水的流动)。材料的物性参数如表7.2所示。

表7.2 材料物性参数表

热性能 导热系数 密度 比热 单位制 W/m℃ Kg/m3 J/kg℃ 铜 383 8889 1390 铁 70 7800 448 水 0.6 1000 4200 17

第7章 非稳态热分析

图7-32 物理模型示意图

该问题是典型的非稳态传热问题，研究对象为3个：铜块、铁块和水。由于分析对象的长度远远大于其横截面积，因此只需对其横截面进行分析即可。本例选取PLANE77单元进行求解。

——附带光盘“Ch7\\实例7－2_start” ——附带光盘“Ch7\\实例7－2_end” ——附带光盘“AVI\\Ch7\\7－2.avi”

1、定义工作文件名

选择Utility Menu>File>Change Jobname，弹出Change Jobname对话框。在对话框中将工作名改为example7-2，单击OK关闭该对话框。选择Main Menu>Preferences，弹出Preferences for GUI Filtering对话框，选中Thermal复选框，然后单击OK按钮关闭该对话框。 2、定义单元类型

选择Main Menu>Preprocessor>Element Type>Add/Edit/Delete，弹出Element Type对话框，单击Add按钮，弹出Library of Element Types对话框。在两个文本框中分别选择Thermal Solid、8node 77,点击OK按钮关闭该对话框。单击Element Type对话框上的Close按钮关闭该对话框。 3、创建几何模型

1、选择Main Menu>Preprocessor>Modeling>Create>Areas>Rectangle>By Dimensions，弹出Create Rectangle by Dimensions对话框，在X1、X2、Y1、Y2文本框中分别输入0、0.6、0、0.5，如图7-33所示。

2、单击Apply按钮，在X1、X2、Y1、Y2文本框中分别输入0.15、0.225、0.225、0.27，如图7-34所示。

3、单击Apply按钮，在X1、X2、Y1、Y2文本框中分别输入0.6-0.2-0.058、0.6-0.2、0.225、0.225+0.044，如图7-35所示。然后单击OK按钮完成设置。

18

ANSYS流体及热场分析

图7-33 建立长方形对话框1 图7-34 建立长方形对话框2

图7-35 建立长方形对话框3

4、选择Main Menu>Preprocessor>Modeling>Operate>Booleans>Overlap>Areas，弹出Overlap Areas菜单，单击Pick All按钮。

5、选择Utility Menu>PlotCtrls>Numbering，弹出Plot Numbering Controls对话框，点选AREA复选框，使其状态由Off变为On，然后单击OK按钮完成设置。

6、选择Main Menu>Preprocessor>Numbering Ctrls>Compress Numbers，弹出Compress Numbers对话框，在Label下拉列表框中选择Areas选项，如图7-36所示，单击OK按钮关闭该对话框。

图7-36 压缩编号对话框

7、选择Utility Menu>Plot>Areas，完成后的平面几何模型如图7-37所示，

图7-37 几何模型

4、定义材料性能参数

1、选择Main Menu>Preprocessor>Material Props>Material Models，弹出Define Material Model Behavior对话框。

2、在Material Models Available列表框依次双击Thermal、Conductivity、Isotropic选项，弹出

19

第7章 非稳态热分析

Conductivity for Material Number 1对话框，在KXX文本框中输入铜的导热系数383，如图7-38所示，单击OK按钮关闭对话框。

3、双击Specific Heat选项，弹出Specific Heat for Material Number 1对话框，在C文本框中输入铜的比热1390，如图7-39所示，单击OK按钮关闭对话框。

4、双击Density选项，弹出Density for Material Number 1对话框，在DENS文本框中输入铜的密度8889，如图7-40所示，单击OK按钮关闭对话框。

图7-38 定义导热系数对话框 图7-39 定义比热对话框

图7-40 定义密度对话框

5、在Define Material Model Behavior对话框中单击Material>New Model，弹出Define Material ID对话框，在Define Material ID文本框中输入2，单击OK按钮关闭该对话框。

6、在Material Models Available列表框依次双击Thermal、Conductivity、Isotropic选项，弹出Conductivity for Material Number 2对话框，在KXX文本框中输入铁的导热系数70，单击OK按钮关闭对话框。

7、双击Specific Heat选项，弹出Specific Heat for Material Number 2对话框，在C文本框中输入铁的比热448，单击OK按钮关闭对话框。

8、双击Density选项，弹出Density for Material Number 2对话框，在DENS文本框中输入铁的密度7800，单击OK按钮关闭对话框。

9、在Define Material Model Behavior对话框中单击Material>New Model，弹出Define Material ID对话框，在Define Material ID文本框中输入3，单击OK按钮关闭该对话框。

10、在Material Models Available列表框依次双击Thermal、Conductivity、Isotropic选项，弹出Conductivity for Material Number 3对话框，在KXX文本框中输入水的导热系数0.61，单击OK按钮关闭对话框。

11、双击Specific Heat选项，弹出Specific Heat for Material Number 3对话框，在C文本框中输入水的比热4200，单击OK按钮关闭对话框。

12、双击Density选项，弹出Density for Material Number 3对话框，在DENS文本框中输入水的密度1000，单击OK按钮关闭对话框。最后关闭Define Material Model Behavior对话框。 5、划分有限元网格

1、选择Main Menu>Preprocessor>Meshing>Size Cntrls>ManualSize>Global>Size，弹出Global

20

第7章 非稳态热分析

图7-51 结果文件写入控制对话框对话框

19、选择Utility Menu>Select>Everything。选择Utility Menu>file>Save as命令,弹出Save Dtabase对话框,在Save Dtabase to文本框中输入7－2_start.db, 保存上述操作过程，单击OK按钮关闭该对话框。

20、选择Main Menu>Solution>Solve>Current LS，出现Solve Current Load Step对话框，同时出现/STAT Command窗口，仔细阅读/STAT Command窗口中的内容，然后点击Close按钮，关闭/STAT Command窗口。点击Solve Current Load Step对话框中的OK按钮，ANSYS开始求解计算。求解结束后，ANSYS显示窗口出现Note提示框，单击Close按钮关闭该对话框。

21、选择Utility Menu>file>Save as命令,弹出Save Dtabase对话框,在Save Dtabase to文本框中输入7－2_end.db, 保存求解结果，单击OK按钮关闭该对话框。 7、后处理

1、选择Main Menu>General Postproc>Read Results>Last Set。选择Utility Menu>Select>Entities，弹出Select Entities对话框。在第1个下拉列表框中选择Elements选项，在第2个下拉列表框中选择By Attributes选项，在第3个选项组中选择Material num单选按钮，在文本框中输入1，单击OK按钮关闭该对话框。

2、选择Utility Menu>Select>Entities，弹出Select Entities对话框。在第1个下拉列表框中选择Nodes选项，在第2个下拉列表框中选择Attached to选项，在第3个选项组中选择Elements单选按钮，单击OK按钮关闭该对话框。 3、选择Main Menu>General Postproc>Plot Results>Contour Plot>Nodal Solu，弹出Contour Nodal Solution Data对话框。选择Nodal Solu>DOF Solution>Nodal Temperature，单击OK按钮，铜块内部的温度场分布如图7-52所示。

4、选择Utility Menu>Select>Entities，弹出Select Entities对话框。在第1个下拉列表框中选择Elements选项，在第2个下拉列表框中选择By Attributes选项，在第3个选项组中选择Material num单选按钮，在文本框中输入2，单击OK按钮关闭该对话框。

5、选择Utility Menu>Select>Entities，弹出Select Entities对话框。在第1个下拉列表框中选择Nodes选项，在第2个下拉列表框中选择Attached to选项，在第3个选项组中选择Elements单选按钮，单击OK按钮关闭该对话框。 6、选择Main Menu>General Postproc>Plot Results>Contour Plot>Nodal Solu，弹出Contour Nodal Solution Data对话框。选择Nodal Solu>DOF Solution>Nodal Temperature，单击OK按钮，铁块内部的温度场分布如图7-53所示。

图7-52 铜块内部温度场分布云图 图7-53 铁块内部温度场分布云图 5、选择Utility Menu>Select>Everything。选择Main Menu>General Postproc>Plot Results>Contour Plot>Nodal Solu，弹出Contour Nodal Solution Data对话框。选择Nodal Solu>DOF Solution>Nodal Temperature，单击OK按钮，整体的温度场分布如图7-54所示。 26

ANSYS流体及热场分析

图7-54 整体温度分布云图

实例7－3：高温铜导线冷却过程分析

如图7-55所示，初始温度为300℉的铜导线置于环境温度为70℉的空气中，其材料的物性参数如表7.3所示。试求5分钟后的温度场分布。

表7.3 材料物性参数表

密度ρ lb/ft3 500 比热C Btu/lb·℉ 0.1

对流换热系数h Btu/ hr·ft2·℉ 2 半径r ft 0.00125 该问题是典型的非稳态传热问题，在计算过程中选用LINK34单轴热对流单元和MASS71质量单元进行求解。MASS71单元需要输入热容实常数Ca，LINK34单元需要输入面积实常数A。Ca和A可分别由下式求解：

Ca???C?V?500?0.1???0.001252?2.45?10?4BTU/oF

A?2?r?2???0.00125?7.85?10?3ft2

——附带光盘“Ch7\\实例7－3_start” ——附带光盘“Ch7\\实例7－3_end” ——附带光盘“AVI\\Ch7\\7－3.avi”

1、定义工作文件名

选择Utility Menu>File>Change Jobname，弹出Change Jobname对话框。在对话框中将工作名改为example7-3，单击OK关闭该对话框。选择Main Menu>Preferences，弹出Preferences for GUI Filtering对话框，选中Thermal复选框，然后单击OK按钮关闭该对话框。

27

第7章 非稳态热分析

2、定义单元类型

1、选择Main Menu>Preprocessor>Element Type>Add/Edit/Delete，弹出Element Type对话框，单击Add按钮，弹出Library of Element Types对话框。在两个列表框中分别选择Thermal Mass、3D mass 71,点击OK按钮关闭该对话框。单击Apply按钮，在Library of Element Types的两个列表框中分别选择Thermal Link、3D convection 34,在Elemental type reference number文本框中输入2，点击OK按钮关闭该对话框。

2、在Element Type对话框上选择Type 1 Mass71选项，单击Option按钮，弹出MASS71 element type options对话框，在Interpret real const CON1 as K3下拉列表框中选择Therm capacitnce选项，如图7-56所示，单击OK按钮关闭该对话框。单击Element Type对话框上的Close按钮关闭该对话框。

图7-56 MASS71单元属性设置对话框

1、选择Main Menu>Preprocessor>Real Constants>Add/Edit/Delete，弹出Real Constants对话框，单击Add按钮，弹出Element Type for R对话框，选择Type 1 MASS71选项，单击OK按钮，弹出Real Constants Set Number 1，for MASS71对话框，参照图7-57进行设置，单击OK关闭该对话框。

2、在Real Constants对话框上单击Add按钮，弹出Element Type for R对话框，选择Type 2 LINK34选项，单击OK按钮，弹出Real Constants Set Number 2，for LINK34对话框，参照图7-58进行设置，单击OK关闭该对话框。单击Real Constants对话框上的Close按钮关闭该对话框。

图7-57 MASS71单元实常数设置对话框 图7-58 LINK34单元实常数设置对话框 3、定义材料性能参数

1、选择Main Menu>Preprocessor>Material Props>Material Models，弹出Define Material Model Behavior对话框。

2、在Material Models Available列表框依次双击Thermal、Convection or Film Coef.选项，弹出Convection or Film coefficient for Material Number 1对话框，在HF文本框中输入热对流系数2，如图7-59所示，单击OK按钮关闭对话框。在Define Material Model Behavior对话框中选择Material>Exit，关闭该对话框。

28

ANSYS流体及热场分析

图7-59 设置热对流系数对话框

4、加载求解

1、选择Main Menu>Preprocessor>Modeling>Create>Nodes>In Active CS，弹出Create Nodes in Active Coordinate System对话框，在NODE文本框中输入1，单击Apply按钮，在NODE文本框中输入2，单击OK按钮关闭该对话框。

2、选择Main Menu>Preprocessor>Modeling>Create>Elements>Auto Numbered>Thru Nodes，弹出Elements from Nodes菜单，在文本框中输入1，单击OK按钮关闭该菜单。

3、选择Main Menu>Preprocessor>Modeling>Create>Elements>Elem Attributes，弹出Elements Attributes对话框，参照图7-60进行设置，然后单击OK按钮关闭该对话框。

4、选择Main Menu>Preprocessor>Modeling>Create>Elements>Auto Numbered>Thru Nodes，弹出Elements from Nodes菜单，在文本框中输入1，2，单击OK按钮关闭该菜单。

图7-60 设置单元属性对话框

5、定义求解类型

1、选择Main Menu>Solution>Analysis Type>New Analysis，弹出New Analysis对话框，选择分析类型为Transient，单击OK按钮，弹出Transient Analysis对话框，在Solution method选项组中选择Full单选按钮，单击OK关闭该对话框。 7．求解控制

2、选择Main Menu>Solution>Analysis Type>Sol’n Controls，弹出Solution Controls对话框，选

29

第7章 非稳态热分析

择Basic选项卡，参照图7-61进行设置，选择Transient选项卡，参照图7-62进行设置，单击OK按钮关闭该对话框。

图7-61 求解控制基本选项设置选项卡 图7-62 求解控制非稳态选项设置选项卡 3、选择Main Menu>Solution>Define Loads>Apply>Thermal>Temperature>Uniform Temp，弹出Uniform Temperature对话框，在文本框中输入300，如图7-63所示，单击OK按钮关闭该对话框。

图7-63 定义整体温度对话框

4、选择Main Menu>Solution>Define Loads>Apply>Thermal>Temperature>On Nodes，弹出Apply TEMP on Nodes菜单，在文本框中输入2，单击OK按钮，弹出Apply TEMP on Nodes对话框，在Lab2列表框中选择TEMP选项，在文本框中输入70，单击OK按钮关闭该对话框。

5、选择Utility Menu>Select>Everything。选择Utility Menu>file>Save as命令,弹出Save Dtabase对话框,在Save Dtabase to文本框中输入7－3_start.db, 保存上述操作过程，单击OK按钮关闭该对话框。

6、选择Main Menu>Solution>Solve>Current LS，出现Solve Current Load Step对话框，同时出现/STAT Command窗口，仔细阅读/STAT Command窗口中的内容，然后点击Close按钮，关闭/STAT Command窗口。点击Solve Current Load Step对话框中的OK按钮，ANSYS开始求解计算。求解结束后，ANSYS显示窗口出现Note提示框，单击Close按钮关闭该对话框。

7、选择Utility Menu>file>Save as命令,弹出Save Dtabase对话框,在Save Dtabase to文本框中输入7－3_end.db, 保存求解结果，单击OK按钮关闭该对话框。 6、后处理

1、选择Main Menu>General Postproc>Read Result>Last Set，读最后一步的分析结果。选择Utility Menu>PlotCtrls>Style>Graphs>Modify Axes，弹出Axes Modifications for Graph Plots对话框，在[/AXLAB] X-axis label文本框中输入TIME（hour），在[/AXLAB] Y-axis label文本框中输入TEMP，单击OK关闭该对话框。

2、选择Main Menu>TimeHist Postpro>Define Variables，弹出Defined Time-history Variables对话框，单击其上的Add按钮，弹出Add Time-History Variable对话框，选中Nodal DOF result单选按钮，单击OK按钮，弹出Define Nodal Data菜单，在文本框中输入1，单击OK按钮，弹出Define Nodal Data对话框，参照图7-64进行设置，单击OK按钮关闭该对话框。然后单击Defined Time-history Variables对话框上的上的Close按钮关闭该对话框。 30

ANSYS流体及热场分析

图7-64 定义节点变量对话框

3、选择Main Menu>TimeHist Postpro>Graph Variables，弹出Graph Time-History Variables对话框，在文本框中输入2，单击OK按钮，ANSYS显示窗口将显示铜导线温度随时间的变化曲线，如图7-65所示。

图7-65 导线温度随时间变化关系曲线图

31

本文来源：https://www.bwwdw.com/article/e6mg.html