Chapter-11：ANSYS 8.1热力学分析

圆筒轴线作横截面，得到一圆环，取其中 1 度进行分析，如图 11.3 所示。

3．查看结果 ANSYS 将热分析的结果将写入*.rth 文件中，它包含基本结果数据和导出结果数据。基 本结果数据为节点温度；导出结果数据有：节点及单元的热流密度、节点及单元的热梯度、 单元热流率和节点的反作用热流率等。 对于稳态热分析，可以使用 POST1 进行后处理。和基本的静力学分析一样，它也可以 输出等值线图、矢量图和数据列表。
第十一章：ANSYS 8.1 热力学分析
热流密度载荷） ，如果输入的值为正，代表热流流入节点，即单元获取热量。如果温度与热 流率同时施加在一节点上则 ANSYS 读取温度值进行计算。 单击 Main Menu>Solution>Define Loads>Apply>Thermal>Heat Flow>On Nodes 菜单，可 在节点上施加热流率。
ANSYS 有限元分析实用教程
http://www.lubanren.com
第十一章
ANSYS 8.1 热力学分析
热力学分析是用于计算一个系统或部件的温度分布及其它热物理参数，如热量的获取 或损失、热梯度、热流密度（热通量〕等。热力学分析在许多工程应用中扮演重要角色， 如内燃机、涡轮机、换热器、管路系统、电子元件等。 ANSYS 热力学分析基于能量守恒原理的热平衡方程，用有限元法计算各节点的温度， 并导出其它热物理参数。在 ANSYS/Multiphysics、ANSYS/Mechanical、ANSYS/Thermal、 ANSYS/FLOTRAN、ANSYS/ED 五种产品中包含热分析功能，其中 ANSYS/FLOTRAN 不 含相变热分析。 热力学分析包括热传导、热对流及热辐射三种热传递方式。此外，还可以分析相变、 有内热源、接触热阻等问题。 ANSYS 热力学分析主要分为以下两类： 稳态传热：系统的温度场不随时间变化； 瞬态传热：系统的温度场随时间明显变化。 本章将分别介绍这两种分析。
图11.2
稳态热分析选项
可以直接在实体模型或单元模型上施加五种载荷（边界条件） ： 恒定温度 恒 定 温 度 通 常 作 为 自 由 度 约 束 施 加 于 温 度 已 知 的 边 界 上 。 单 击 Main Menu>Solution>Define Loads>Apply>Thermal>Temperature>On Nodes 菜单， 在节点上施加温 度边界条件。 热流率 热流率作为节点集中载荷，主要用于线单元模型中（通常线单元模型不能施加对流或

11.2.3 ANSYS 稳态热分析的基本过程
ANSYS 进行稳态热分析可分为三个步骤： 建模 施加荷载并计算 查看结果 1．建模 这一步主要是在前处理器（PREP7）中进行，包括定义单元类型和单元选项、定义实常 数、定义材料性能参数和创建几何模型并划分网格等。 其中，定义材料热性能参数，对于稳态传热，一般只需定义导热系数，它可以是恒定 的，也可以随温度变化。 2．施加荷载并计算 这一步主要是在求解模块（SOLU）中进行，包括定义分析类型、确定分析选项、施加 荷载及荷载步等。 定义分析类型，可单击 Main Menu>Solution>Analysis Type>New Analysis 菜单，弹出如 图 11.1 所示的对话框。在【Type of analysis】单选框中选中【Steady-State】单选按钮，单击 【OK】按钮即可。
11.2.4 稳态热分析实例
本节将结合一个实例介绍用 ANSYS 进行稳态热分析的基本过程。 1．问题描述 某一潜水艇可以简化为一圆筒，它由三层组成，最外面一层为不锈钢，中间为玻纤隔 热层，最里面为铝层，筒内为空气，筒外为海水，求内外壁面温度及温度分布。沿垂直于
ANSYS 有限元分析实用教程
http://www.lubanren.com
dU ，即流入或流出的热传递速率 q 等于系统内能的变化。 dt
11.1.3 热传递的方式
1．热传导 热传导可以定义为完全接触的两个物体之间或一个物体的不同部分之间由于温度梯度 而引起的内能的交换。 热传导遵循付里叶定律：q k
dT ， 式中 q 为热流密度 （W/m2） ， dx
， “-”表示热量流向温度降低的方向。 k 为导热系数（W/m-℃） 2．热对流 热对流是指固体的表面与它周围接触的流体之间，由于温差的存在引起的热量的交换。 热对流可以分为两类：自然对流和强制对流。热对流用牛顿冷却方程来描述： q h(TS TB ) ，式中 h 为对流换热系数（或称膜传热系数、给热系数、膜系数等） ，TS 为 固体表面的温度， TB 为周围流体的温度。 3．热辐射 热辐射指物体发射电磁能，并被其它物体吸收转变为热的热量交换过程。物体温度越 高，单位时间辐射的热量越多。热传导和热对流都需要有传热介质，而热辐射无须任何介 质。实质上，在真空中的热辐射效率最高。 在工程中通常考虑两个或两个以上物体之间的辐射，系统中每个物体同时辐射并吸收 热量。它们之间的净热量传递可以用斯蒂芬-波尔兹曼方程来计算： ，式中 q 为热流率， 为辐射率（黑度） ， 为斯蒂芬－波尔兹曼常数，约为 5.67×10-8W/m2.K4，A1 为辐射面 1 的面积， F12 为由辐射面 1 到辐射面 2 的形状系数， T1 为辐射面 1 的绝对温度， T2 为辐射 面 2 的绝对温度。由上式可以看出，包含热辐射的热分析是高度非线性的。
11.1 热分析的基础Байду номын сангаас识
要进行热力学分析需要用户具有热力学分析的基础知识，本节简单回顾一下热力学分 析的基础知识。已经掌握了这部分知识的读者跳过本节的学习。
11.1.1 符号与单位
热力学中常用的符号和单元如表 11.1 所示。
表11.1 项目 长度 时间 质量 温度 力 能量（热量） 功率（热流率） 热流密度 生热速率 导热系数 对流系数 密度 比热 焓 热力学分析中的符号与单位 英制单位 ft s lbm o F 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 代号 国际单位 m s Kg ℃ N J W W/m2 W/m3 W/m-℃ W/m2-℃ Kg/m3 J/Kg-℃ J/m3
注意：如果在实体单元的某一节点上施加热流率，则此节点周围的单元要密一些，在两种导热系数差别很 大的两个单元的公共节点上施加热流率时，尤其要注意。此外，尽可能使用热生成或热流密度边界条件， 这样结果会更精确些。
对流 对流边界条件作为面载施加于实体的外表面，计算与流体的热交换，它仅可施加于实 体和壳模型上，对于线模型，可以通过对流线单元 LINK34 考虑对流。 单击 Main Menu>Solution>Define Loads>Apply>Thermal>Convection>On Nodes 菜单， 可 在节点上施加对流。 热流密度 热流密度也是一种面载。当通过单位面积的热流率已知或通过 FLOTRAN CFD 计算得 到时，可以在模型相应的外表面施加热流密度。如果输入的值为正，代表热流流入单元。 热流密度也仅适用于实体和壳单元。热流密度与对流可以施加在同一外表面，但 ANSYS 仅 读取最后施加的面载进行计算。 单击 Main Menu>Solution>Define Loads>Apply>Thermal>Heat Flux>On Nodes 菜单，可 在节点上施加热流密度。 生热率 生热率作为体载施加于单元上，可以模拟化学反应生热或电流生热。它的单位是单位 体积的热流率。 单击 Main Menu>Solution>Define Loads>Apply>Thermal>Heat Generat>On Nodes 菜单， 可在节点上施加生热率。
11.2.2 热分析的单元
在 ANSYS 中，热分析涉及到的单元有大约 40 种，其中纯粹用于热分析的有 14 种： 线性 4 种如下： LINK32：二维二节点热传导单元 LINK33：三维二节点热传导单元 LINK34：二节点热对流单元
第十一章：ANSYS 8.1 热力学分析
LINK31：二节点热辐射单元 二维实体 5 种如下： PLANE55：四节点四边形单元 PLANE77：八节点四边形单元 PLANE35：三节点三角形单元 PLANE75：四节点轴对称单元 PLANE78：八节点轴对称单元 三维实体 3 种如下： SOLID87：六节点四面体单元 SOLID70：八节点六面体单元 SOLID90：二十节点六面体单元 壳单元 1 种：SHELL57 点单元 1 种：MASS71 有关单元的详细解释，请读者参考 ANSYS 的帮助文档。
KXX HF DENS C ENTH
第十一章：ANSYS 8.1 热力学分析
11.1.2 传热学经典理论回顾
热力学分析遵循热力学第一定律，即能量守恒定律，其描述如下： 对于一个封闭的系统（没有质量的流入或流出〕有
Q W U KE PE
其中： Q—— 热量； W—— 作功； U ——系统内能； ——系统动能； PE ——系统势能； 对于大多数工程传热问题： KE ＝ PE ＝ 0 ； 通常考虑没有做功： W 0 ，则： Q U ； 对于稳态热分析： Q U 0 ，即流入系统的热量等于流出的热量； 对于瞬态热分析： q
ANSYS 有限元分析实用教程
http://www.lubanren.com
图11.1
定义分析类型
确认分析选项， 可单击 Main Menu>Solution>Analysis Type>Analysis Options 菜单， 弹出 如图 11.2 所示的 【Static or Steady-State Analysis】 对话框， 如图 11.2 所示。 在 【Newton – Raphson option】下拉列表框中选择牛顿-辛普森算法（当进行非线性分析时才需设置） ；在【Equaton solver】下拉列表框中选择求解器；在【Temperature difference】文本框中输入绝对零度值， 在进行热辐射分析时，要将目前的温度值换算为绝对温度。如果使用的温度单位是摄氏度， 此值应设定为【273】 ，如果使用的是华氏度，则为【460】 。
KT Q
其中： K 为传导矩阵，包含导热系数、对流系数及辐射率和形状系数；
ANSYS 利用模型几何参数、材料热性能参数以及所施加的边界条件，生成 K 、 T 以及 Q 。 5．瞬态传热 瞬态传热过程是指一个系统的加热或冷却过程。在这个过程中系统的温度、热流率、 热边界条件以及系统内能随时间都有明显变化。根据能量守恒原理，瞬态热平衡可以表达 为（以矩阵形式表示） ： 其中： K 为传导矩阵，包含导热系数、对流系数及辐射率和形状系数；
Q 为节点热流率向量，包含热生成；
T 为节点温度向量；
C 为比热矩阵,考虑系统内能的增加; T 为节点温度向量；
为温度对时间的导数;
Q 为节点热流率向量，包含热生成。
11.2 稳态热分析
11.2.1 稳态传热的定义
稳态传热用于分析稳定的热载荷对系统或部件的影响。通常在进行瞬态热分析以前， 进行稳态热分析用于确定初始温度分布。稳态热分析可以通过有限元计算确定由于稳定的 热载荷引起的温度、热梯度、热流率、热流密度等参数。
ANSYS 有限元分析实用教程
http://www.lubanren.com
4．稳态传热 如果系统的净热流率为０，即流入系统的热量加上系统自身产生的热量等于流出系统 的热量：q 流入+q 生成-q 流出=0，则系统处于热稳态。在稳态热分析中任一节点的温度不随时间 变化。稳态热分析的能量平衡方程为（以矩阵形式表示）