ANSYS热分析指南与经典案例.docx

ANSYS热分析指南与经典案例.docx

《ANSYS热分析指南与经典案例.docx》由会员分享，可在线阅读，更多相关《ANSYS热分析指南与经典案例.docx（44页珍藏版）》请在冰豆网上搜索。

ANSYS热分析指南与经典案例

第一章简介

一、热分析的目的

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

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

二、ANSYS的热分析

∙在ANSYS/Multiphysics、ANSYS/Mechanical、ANSYS/Thermal、ANSYS/FLOTRAN、ANSYS/ED五种产品中包含热分析功能,其中ANSYS/FLOTRAN不含相变热分析。

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

∙ANSYS热分析包括热传导、热对流及热辐射三种热传递方式。

此外，还可以分析相变、有内热源、接触热阻等问题。

三、ANSYS热分析分类

∙稳态传热：

系统的温度场不随时间变化

∙瞬态传热：

系统的温度场随时间明显变化

四、耦合分析

∙热－结构耦合

∙热－流体耦合

∙热－电耦合

∙热－磁耦合

∙热－电－磁－结构耦合等

第二章基础知识

一、符号与单位

项目

国际单位

英制单位

ANSYS代号

长度

m

ft

时间

s

s

质量

Kg

lbm

温度

℃

oF

力

N

lbf

能量（热量）

J

BTU

功率（热流率）

W

BTU/sec

热流密度

W/m2

BTU/sec-ft2

生热速率

W/m3

BTU/sec-ft3

导热系数

W/m-℃

BTU/sec-ft-oF

KXX

对流系数

W/m2-℃

BTU/sec-ft2-oF

HF

密度

Kg/m3

lbm/ft3

DENS

比热

J/Kg-℃

BTU/lbm-oF

C

焓

J/m3

BTU/ft3

ENTH

二、传热学经典理论回顾

热分析遵循热力学第一定律，即能量守恒定律:

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

式中:

Q——热量;

W——作功;

——系统内能;

——系统动能；

——系统势能；

●对于大多数工程传热问题：

；

●通常考虑没有做功：

则：

；

●对于稳态热分析：

，即流入系统的热量等于流出的热量；

●对于瞬态热分析：

，即流入或流出的热传递速率q等于系统内能的变化。

三、热传递的方式

1、热传导

　　热传导可以定义为完全接触的两个物体之间或一个物体的不同部分之间由于温度梯度而引起的内能的交换。

热传导遵循付里叶定律：

，式中

为热流密度（W/m2），

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

2、热对流

　　热对流是指固体的表面与它周围接触的流体之间，由于温差的存在引起的热量的交换。

热对流可以分为两类：

自然对流和强制对流。

热对流用牛顿冷却方程来描述：

，式中h为对流换热系数（或称膜传热系数、给热系数、膜系数等），

为固体表面的温度，

为周围流体的温度。

3、热辐射

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

物体温度越高，单位时间辐射的热量越多。

热传导和热对流都需要有传热介质，而热辐射无须任何介质。

实质上，在真空中的热辐射效率最高。

在工程中通常考虑两个或两个以上物体之间的辐射，系统中每个物体同时辐射并吸收热量。

它们之间的净热量传递可以用斯蒂芬—波尔兹曼方程来计算：

，式中

为热流率，

为辐射率（黑度），

为斯蒂芬－波尔兹曼常数，约为×10-8W/m2.K4，A1为辐射面1的面积，

为由辐射面1到辐射面2的形状系数，

为辐射面1的绝对温度，

为辐射面2的绝对温度。

由上式可以看出，包含热辐射的热分析是高度非线性的。

四、稳态传热

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

q流入+q生成-q流出=0，则系统处于热稳态。

在稳态热分析中任一节点的温度不随时间变化。

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

式中：

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

为节点温度向量；

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

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

、

以及

。

五、瞬态传热

　　瞬态传热过程是指一个系统的加热或冷却过程。

在这个过程中系统的温度、热流率、热边界条件以及系统内能随时间都有明显变化。

根据能量守恒原理，瞬态热平衡可以表达为（以矩阵形式表示）：

式中:

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

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

为节点温度向量；

为温度对时间的导数;

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

六、线性与非线性

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

①、材料热性能随温度变化，如K（T）,C（T）等；

②、边界条件随温度变化，如h（T）等；

③、含有非线性单元；

4、考虑辐射传热

非线性热分析的热平衡矩阵方程为：

七、边界条件、初始条件

ANSYS热分析的边界条件或初始条件可分为七种：

温度、热流率、热流密度、对流、辐射、绝热、生热。

八、热分析误差估计

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

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

∙基于单元边界的热流密度的不连续；

∙仅对一种材料、线性、稳态热分析有效；

∙使用自适应网格划分可以对误差进行控制。

第三章　稳态传热分析

一、稳态传热的定义

　　稳态传热用于分析稳定的热载荷对系统或部件的影响。

通常在进行瞬态热分析以前，进行稳态热分析用于确定初始温度分布。

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

二、热分析的单元

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

线性：

LINK32两维二节点热传导单元

LINK33三维二节点热传导单元

LINK34二节点热对流单元

LINK31二节点热辐射单元

二维实体：

PLANE55四节点四边形单元

PLANE77八节点四边形单元

PLANE35三节点三角形单元

PLANE75四节点轴对称单元

PLANE78八节点轴对称单元

三维实体SOLID87六节点四面体单元

SOLID70八节点六面体单元

SOLID90二十节点六面体单元

壳SHELL57四节点

点MASS71

　　有关单元的详细解释，请参阅《ANSYSElementReferenceGuide》

三、ANSYS稳态热分析的基本过程

　　ANSYS热分析可分为三个步骤：

∙前处理：

建模

∙求解：

施加载荷计算

∙后处理：

查看结果

1、建模

①、确定jobname、title、unit；

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

③、定义单元实常数；

④、定义材料热性能参数，对于稳态传热，一般只需定义导热系数，它可以是恒定的，也可以随温度变化；

⑤、创建几何模型并划分网格，请参阅《ANSYSModelingandMeshingGuide》。

2、施加载荷计算

1、定义分析类型

●如果进行新的热分析：

Command:

ANTYPE,STATIC,NEW

GUI:

Mainmenu>Solution>-AnalysisType->NewAnalysis>Steady-state

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

Command:

ANTYPE,STATIC,REST

GUI:

Mainmenu>Solution>AnalysisType->Restart

2、施加载荷

　　可以直接在实体模型或单元模型上施加五种载荷（边界条件）：

a、恒定的温度

　　通常作为自由度约束施加于温度已知的边界上。

CommandFamily:

　D

GUI：

MainMenu>Solution>-Loads-Apply>-Thermal-Temperature

b、热流率

　　热流率作为节点集中载荷，主要用于线单元模型中（通常线单元模型不能施加对流或热流密度载荷），如果输入的值为正，代表热流流入节点，即单元获取热量。

如果温度与热流率同时施加在一节点上则ANSYS读取温度值进行计算。

注意：

如果在实体单元的某一节点上施加热流率，则此节点周围的单元要密一些，在两种导热系数差别很大的两个单元的公共节点上施加热流率时，尤其要注意。

此外，尽可能使用热生成或热流密度边界条件，这样结果会更精确些。

CommandFamily:

F

GUI：

MainMenu>Solution>-Loads-Apply>-Thermal-HeatFlow

c、对流

　　对流边界条件作为面载施加于实体的外表面，计算与流体的热交换，它仅可施加于实体和壳模型上，对于线模型，可以通过对流线单元LINK34考虑对流。

CommandFamily:

SF

GUI：

MainMenu>Solution>-Loads-Apply>-Thermal-Convection

d、热流密度

　　热流密度也是一种面载。

当通过单位面积的热流率已知或通过FLOTRANCFD计算得到时，可以在模型相应的外表面施加热流密度。

如果输入的值为正，代表热流流入单元。

热流密度也仅适用于实体和壳单元。

热流密度与对流可以施加在同一外表面，但ANSYS仅读取最后施加的面载进行计算。

CommandFamily:

F

GUI：

MainMenu>Solution>-Loads-Apply>-Thermal-HeatFlux

e、生热率

　　生热率作为体载施加于单元上，可以模拟化学反应生热或电流生热。

它的单位是单位体积的热流率。

CommandFamily:

BF

GUI：

MainMenu>Solution>-Loads-Apply>-Thermal-HeatGenerat

③、确定载荷步选项

　　对于一个热分析，可以确定普通选项、非线性选项以及输出控制。

a.普通选项

∙时间选项：

虽然对于稳态热分析，时间选项并没有实际的物理意义，但它提供了一个方便的设置载荷步和载荷子步的方法。

Command:

TIME

GUI:

MainMenu>Solution>-LoadStepOpts-Time/Frequenc>Time-TimeStep/TimeandSubstps

∙每载荷步中子步的数量或时间步大小：

对于非线性分析，每一载荷步需要多个子步。

Command:

NSUBST

GUI:

MainMenu>Solution>-LoadStepOpts->Time/Frequenc>TimeandSubstps

Command:

DELTIM

GUI:

MainMenu>Solution>-LoadStepOpts->Time/Frequenc>Time-TimeStep

∙递进或阶越选项：

如果定义阶越（stepped）选项，载荷值在这个载荷步内保持不变；如果为递进（ramped）选项，则载荷值由上一载荷步值到本载荷步值随每一子步线性变化。

Command:

KBC

GUI:

MainMenu>Solution>-LoadStepOpts-Time/Frequenc>Time-TimeStep/TimeandSubstps

b.非线性选项

∙迭代次数：

本选项设置每一子步允许的最多的迭代次数。

默认值为25，对大数热分析问题足够。

Command:

NEQIT

GUI:

MainMenu>Solution>-LoadStepOpts-Nolinear>EquilibriumIter

∙自动时间步长:

对于非线性问题，可以自动设定子步间载荷的增长，保证求解的稳定性和准确性。

Command:

AUTOTS

GUI:

MainMenu>Solution>-LoadStepOpts-Time/Frequenc>Time-TimeStep/TimeandSubstps

∙收敛误差：

可根据温度、热流率等检验热分析的收敛性。

Command:

CNVTOL

GUI:

MainMenu>Solution>-LoadStepOpts-Nolinear>ConvergenceCrit

∙求解结束选项：

如果在规定的迭代次数内，达不到收敛，ANSYS可以停止求解或到下一载荷步继续求解。

Command:

NCNV

GUI:

MainMenu>Solution>-LoadStepOpts-Nolinear>CriteriatoStop

∙线性搜索：

设置本选项可使ANSYS用Newton-Raphson方法进行线性搜索。

Command:

LNSRCH

GUI:

MainMenu>Solution>-LoadStepOpts-Nolinear>LineSearch

∙预测矫正：

本选项可激活每一子步第一次迭代对自由度求解的预测矫正。

Command:

PRED

GUI:

MainMenu>Solution>-LoadStepOpts-Nolinear>Predictor

c.输出控制

∙控制打印输出：

本选项可将任何结果数据输出到*.out文件中。

Command:

OUTPR

GUI:

MainMenu>Solution>-LoadStepOpts-OutputCtrls>SoluPrintout

∙控制结果文件：

控制*.rth的内容。

Command:

OUTRES

GUI:

MainMenu>Solution>-LoadStepOpts-OutputCtrls>DB/ResultsFile

④、确定分析选项

a.Newton-Raphson选项（仅对非线性分析有用）

Command:

NROPT

GUI:

MainMenu>Solution>AnalysisOptions

b.选择求解器：

可选择如下求解器中一个进行求解：

∙Frontalsolver（默认）

∙JacobiConjugateGradient（JCG）solver

∙JCGout-of-memorysolver

∙IncompleteCholeskyConjugateGradient（ICCG）solver

∙Pre-ConditionedConjugateGradientSolver（PCG）

∙Iterative（automaticsolverselectionoption）

Command:

EQSLV

GUI:

MainMenu>Solution>AnalysisOptions

注意：

热分析可选用Iterative选项进行快速求解，但如下情况除外：

∙热分析包含SURF19或SURF22或超单元；

∙热辐射分析；

∙相变分析

∙需要restartananalysis

c.确定绝对零度：

在进行热辐射分析时，要将目前的温度值换算为绝对温度。

如果使用的温度单位是摄氏度，此值应设定为273；如果使用的是华氏度，则为460。

Command:

TOFFST

GUI:

MainMenu>Solution>AnalysisOptions

5、保存模型:

点击ANSYS工具条SAVE_DB。

⑥、求解

Command:

SOLVE

GUI:

MainMenu>Solution>CurrentLS

3、后处理

　　ANSYS将热分析的结果写入*.rth文件中，它包含如下数据：

　　基本数据：

∙节点温度

　　导出数据：

∙节点及单元的热流密度

∙节点及单元的热梯度

∙单元热流率

∙节点的反作用热流率

∙其它

　　对于稳态热分析，可以使用POST1进行后处理，关于后处理的完整描述，可参阅《ANSYSBasicAnalysisProceduresGuide》。

　　进入POST1后，读入载荷步和子步：

Command:

SET

GUI:

MainMenu>GeneralPostproc>-ReadResults-ByLoadStep

　　可以通过如下三种方式查看结果：

∙彩色云图显示

Command:

PLNSOL,PLESOL,PLETAB等

GUI:

MainMenu>GeneralPostproc>PlotResults>NodalSolu,ElementSolu,ElemTable

∙矢量图显示

Command:

PLVECT

GUI:

MainMenu>GeneralPostproc>PlotResults>Pre-definedorUserdefined

∙列表显示

Command:

PRNSOL,PRESOL,PRRSOL等

GUI:

MainMenu>GeneralPostproc>ListResults>NodalSolu,ElementSolu,ReactionSolu

　　详细过程请参阅《ANSYSBasicAnalysisProceduresGuide》。

实例1：

　　某一潜水艇可以简化为一圆筒，它由三层组成，最外面一层为不锈钢，中间为玻纤隔热层，最里面为铝层，筒内为空气，筒外为海水，求内外壁面温度及温度分布。

几何参数：

筒外径30feet

总壁厚2inch

不锈钢层壁厚0.75inch

玻纤层壁厚1inch

铝层壁厚0.25inch

筒长200feet

导热系数不锈钢8.27BTU/hr.ft.oF

玻纤0.028BTU/hr.ft.oF

铝117.4BTU/hr.ft.oF

边界条件空气温度70oF

海水温度44.5oF

空气对流系数2.oF

海水对流系数2.oF

　　沿垂直于圆筒轴线作横截面，得到一圆环，取其中1度进行分析,如图示。

以下分别列出log文件和菜单文件。

/filename,Steady1

/title,Steady-statethermalanalysisofsubmarine

/units,BFT

Ro=15!

外径（ft）

Rss=15-（0.75/12）!

不锈钢层内径ft）

Rins=15-（1.75/12）!

玻璃纤维层内径（ft）

Ral=15-（2/12）!

铝层内径（ft）

Tair=70!

潜水艇内空气温度

Tsea=44.5!

海水温度

Kss=8.27!

不锈钢的导热系数（BTU/hr.ft.oF）

Kins=0.028!

玻璃纤维的导热系数（BTU/hr.ft.oF）

Kal=117.4!

铝的导热系数（BTU/hr.ft.oF）

Hair=2.5!

空气的对流系数（BTU/hr.ft2.oF）

Hsea=80!

海水的对流系数（BTU/hr.ft2.oF）

/prep7

et,1,plane55!

定义二维热单元

mp,kxx,1,Kss!

设定不锈钢的导热系数

mp,kxx,2,Kins!

设定玻璃纤维的导热系数

mp,kxx,3,Kal!

设定铝的导热系数

pcirc,Ro,Rss,-0.5,0.5!

创建几何模型

aglue,all

numcmp,area

lesize,1,,,16!

设定划分网格密度

lesize,4,,,4

lesize,14,,,5

lesize,16,,,2

eshape,2!

设定为映射网格划分

mat,1

amesh,1

mat,2

amesh,2

mat,3

amesh,3

/SOLU

SFL,11,CONV,HAIR,,TAIR!

施加空气对流边界

SFL,1,CONV,HSEA,,TSEA!

施加海水对流边界

SOLVE

/POST1

PLNSOL！

输出温度彩色云图

finish

菜单操作:

1.UtilityMenu>File>changejobename,输入Steady1；

2.UtilityMenu>File>changetitle,输入Steady-statethermalanalysisofsubmarine；

3.在命令行输入:

/units,BFT；

4.MainMenu:

Preprocessor；

5.MainMenu:

Preprocessor>ElementType>Add/Edit/Delete,选择PLANE55；

6.MainMenu:

Preprocessor>MaterialProp>-Constant-Isotropic,默认材料编号为1,在KXX框中输入8.27,选择APPLY,输入材料编号为2,在KXX框中输入0.028,选择APPLY,输入材料编号为3,在KXX框中输入；

7.MainMenu:

Preprocessor>-Modeling->Create>-Areas-Circle>ByDimensions,在RAD1中输入15,在RAD2中输入15-（.75/12）,在THERA1中输入-0.5,在THERA2中输入0.5,选择APPLY,在RAD1中输入15-（.75/12）,在RAD2中输入15-（1.75/12）,选择APPLY,在RAD1中输入15-（1.75/12）,在RAD2中输入15-2/12,选择OK；

8.MainMenu:

Preprocessor>-Modeling->Operate>-Booleane->Glue>Area,选择PICKALL；

9.MainMenu:

Preprocessor>-Meshing-SizeContrls>-Lines-PickedLines,选择不锈钢层短边,在NDIV框中输入4,选择APPLY,选择玻璃纤维层的短边,在NDIV框中输入5,选择APPLY,选择铝层的短边,在NDIV框中输入2,选择APPLY,选择四个长边,在NDIV中输入16；

10.MainMenu:

Preprocessor>-Attributes-Define>PickedArea,选择不锈钢层,在MAT框中输入1,选择APPLY,选择玻璃纤维层,在MAT框中输入2,选择APPLY,选择铝层,在MAT框中输入3,选择OK；

11.MainMenu:

Preprocessor>-Meshing-Mesh>-Areas-Mapped>3or4sided,选择PICKALL；

12.MainMenu:

Solution>-Loads-Apply>-Thermal-Convection>Onlines,选择不锈钢外壁,在VALI框中输入80,在VAL2I框中输入44.5,选择APPLY,选择铝层内壁,在VALI框中输入2.5,在VAL2I框中输入70,选择OK；

13.MainMenu:

Solution>-Solve-CurrentLS；

14.MainMen