第一章

一、热分析之目的

简 介

热分析用于计算一个系统或部件的温度分布及其它热物理参数，如热量的获取或损失、热梯度、热流密度（热通量〕等。

热分析在许多工程应用中扮演重要角色，如内燃机、涡轮机、换热器、管路系统、电子元件等。 二、ANSYS之热分析

?

在

ANSYS/Multiphysics

、

ANSYS/Mechanical

、

ANSYS/Thermal

、

ANSYS/FLOTRAN、ANSYS/ED五种产品中包含热分析功能,其中ANSYS/FLOTRAN不含相变热分析。 ?

ANSYS热分析基于能量守恒原理的热平衡方程，用有限元法计算各节点的温度，并导出其它热物理参数。 ?

ANSYS热分析包括热传导、热对流及热辐射三种热传递方式。此外，还可以分析相变、有内热源、接触热阻等问题。

三、ANSYS 热分析之分类

? ?

稳态传热：系统的温度场不随时间变化 瞬态传热：系统的温度场随时间明显变化

四、耦合分析之

? ? ? ? ?

热－结构耦合 热－流体耦合 热－电耦合 热－磁耦合

热－电－磁－结构耦合等

第二章

一、符号与单位

项目 长度 时间 质量 温度 力 能量（热量） 功率（热流率） 热流密度 生热速率 导热系数 对流系数 密度 比热 焓 二、传热学经典理论回顾 国际单位 m s Kg ℃ N J W W/m2 W/m3 W/m-℃ W/m2-℃ Kg/m3 J/Kg-℃ J/m3 基础知识 英制单位 ft s lbm oF lbf BTU BTU/sec BTU/sec-ft2 BTU/sec-ft3 BTU/sec-ft-oF BTU/sec-ft2-oF lbm/ft3 BTU/lbm-oF BTU/ft3 ANSYS代号 KXX HF DENS C ENTH 热分析遵循热力学第一定律，即能量守恒定律:

?

对于一个封闭的系统（没有质量的流入或流出〕

W —— 作功;

式中: Q —— 热量;

?U——系统内能; ?KE——系统动能； ?PE——系统势能；

? ? ?

对于大多数工程传热问题：?； KE＝?PE＝0??U通常考虑没有做功：W, 则：Q； ?0??U?0对于稳态热分析：Q，即流入系统的热量等于流出的热量；

对于瞬态热分析：

?

q?dU，即流入或流出的热传递速率q等于系统内能的变化。 dt三、热传递的方式

1、热传导

热传导可以定义为完全接触的两个物体之间或一个物体的不同部分之间由于温度梯度而引起的内能的交换。热传导遵循付里叶定律：（W/m2），

dT????qk，式中q??为热流密度

dxk为导热系数（W/m-℃），“-”表示热量流向温度降低的方向。

2、热对流

热对流是指固体的表面与它周围接触的流体之间，由于温差的存在引起的热量的交换。热对流可以分为两类：自然对流和强制对流。热对流用牛顿冷却方程来描述：

??，式中h为对流换热系数（或称膜传热系数、给热系数、膜系数等），q?h(T?T)SBTS为固体表面的温度，TB为周围流体的温度。

3、热辐射

热辐射指物体发射电磁能，并被其它物体吸收转变为热的热量交换过程。物体温

度越高，单位时间辐射的热量越多。热传导和热对流都需要有传热介质，而热辐射无须任何介质。实质上，在真空中的热辐射效率最高。

在工程中通常考虑两个或两个以上物体之间的辐射，系统中每个物体同时辐射并吸收热量。它们之间的净热量传递可以用斯蒂芬—波尔兹曼方程来计算：

q???A1F12(T14?T24)，式中q为热流率，?为辐射率（黑度），?为斯蒂芬－波尔

兹曼常数，约为5.67×10-8W/m2.K4，A1为辐射面1的面积，

F12为由辐射面1到辐射面

2的形状系数，T1为辐射面1的绝对温度，T2为辐射面2的绝对温度。由上式可以看出，包含热辐射的热分析是高度非线性的。 四、稳态传热

如果系统的净热流率为０，即流入系统的热量加上系统自身产生的热量等于流出系统的热量：q

流入

+q

生成

-q

流出

=0，则系统处于热稳态。在稳态热分析中任一节点的温度不随时间变化。稳

态热分析的能量平衡方程为（以矩阵形式表示）

式中：?K?为传导矩阵，包含导热系数、对流系数及辐射率和形状系数；

?T?为节点温度向量；

?Q?为节点热流率向量，包含热生成；

ANSYS利用模型几何参数、材料热性能参数以及所施加的边界条件，生成?K?、

?T?以及?Q?。

五、瞬态传热

瞬态传热过程是指一个系统的加热或冷却过程。在这个过程中系统的温度、热流率、热边界条件以及系统内能随时间都有明显变化。根据能量守恒原理，瞬态热平衡可以表达为（以矩阵形式表示）： 式中:

?K?为传导矩阵，包含导热系数、对流系数及辐射率和形状系数；

?C?为比热矩阵,考虑系统内能的增加;

?T?为节点温度向量；

?T??为温度对时间的导数;

?Q?为节点热流率向量，包含热生成。

①、材料热性能随温度变化，如K(T),C(T)等； ②、边界条件随温度变化，如h(T)等； ③、含有非线性单元；

六、线性与非线性

如果有下列情况产生，则为非线性热分析：

④、 考虑辐射传热

非线性热分析的热平衡矩阵方程为： 七、边界条件、初始条件

ANSYS热分析的边界条件或初始条件可分为七种：温度、热流率、热流密度、对流、辐射、绝热、生热。 八、热分析误差估计 ? ?

仅用于评估由于网格密度不够带来的误差；

仅适用于SOLID或SHELL的热单元（只有温度一个自由度）；

? ? ?

基于单元边界的热流密度的不连续； 仅对一种材料、线性、稳态热分析有效； 使用自适应网格划分可以对误差进行控制。

第三章 稳态传热分析

一、稳态传热的定义

稳态传热用于分析稳定的热载荷对系统或部件的影响。通常在进行瞬态热分析以前，进行稳态热分析用于确定初始温度分布。

稳态热分析可以通过有限元计算确定由于稳定的热载荷引起的温度、热梯度、热流率、热流密度等参数 二、热分析的单元

热分析涉及到的单元有大约40种，其中纯粹用于热分析的有14种：

线性：

LINK32

壳 点

两维二节点热传导单元

三维二节点热传导单元 二节点热对流单元 二节点热辐射单元 四节点四边形单元 八节点四边形单元 三节点三角形单元 四节点轴对称单元 八节点轴对称单元 六节点四面体单元 八节点六面体单元 二十节点六面体单元

LINK33 LINK34 LINK31

二维实体： PLANE55

PLANE77 PLANE35 PLANE75 PLANE78 SOLID87 SOLID70 SOLID90

三维实体

SHELL57 MASS71

四节点

有关单元的详细解释，请参阅《ANSYS Element Reference Guide》 三、ANSYS稳态热分析的基本过程 ANSYS热分析可分为三个步骤：

? ? ?

前处理： 后处理：

建模 查看结果

求解： 施加载荷计算

1、建模

①、确定jobname、title、unit；

②、进入PREP7前处理，定义单元类型，设定单元选项； ③、定义单元实常数；

④、定义材料热性能参数，对于稳态传热，一般只需定义导热系数，它可以是恒定

的，也可以随温度变化；

⑤、创建几何模型并划分网格，请参阅《ANSYS Modeling and Meshing Guide》。

2、施加载荷计算 ①、 定义分析类型

?

如果进行新的热分析：

Command: ANTYPE, STATIC, NEW

GUI: Main menu>Solution>-Analysis Type->New Analysis>Steady-state

?

如果继续上一次分析，比如增加边界条件等：

Command: ANTYPE, STATIC, REST