ansys热分析.docx

1、ansys热分析ANSYS热分析指南(第一章)第一章简介1.1热分析的目的热分析用于计算一个系统或部件的温度分布及其它热物理参数，我们一般关心的参数有：温度的分布热量的增加或损失热梯度热流密度热分析在许多工程应用中扮演着重要角色，如内燃机、涡轮机、换热器、管路系统、电子元件等等。通常在完成热分析后将进行结构应力分析，计算由于热膨胀或收缩而引起的热应力。1.2ANSYS中的热分析ANSYS/Multiphysics、ANSYS/Mechanical、ANSYS/Professional、ANSYS/FLOTRAN四种产品中支持热分析功能。ANSYS热分析基于由能量守恒原理导出的热平衡方程，有关细

2、节，请参阅ANSYS Theory Reference。ANSYS使用有限元法计算各节点的温度，并由其导出其它热物理参数。ANSYS可以处理所有的三种主要热传递方式：热传导、热对流及热辐射。1.2.1对流热对流在ANSYS中作为一种面载荷，施加于实体或壳单元的表面。首先需要输入对流换热系数和环境流体温度，ANSYS将计算出通过表面的热流量。如果对流换热系数依赖于温度，可以定义温度表，以及在每一个温度点处的对流换热系数。1.2.2辐射ANSYS提供了四种方法来解决非线性的辐射问题：辐射杆单元(LINK31)使用含热辐射选项的表面效应单元(SURF151-2D,或SURF152-3D)在AUX12

3、中，生成辐射矩阵，作为超单元参与热分析使用Radiosity求解器方法有关辐射的详细描述请阅读本指南第四章。1.2.3特殊的问题除了前面提到的三种热传递方式外，ANSYS热分析还可以解决一些诸如：相变(熔融与凝固)、内部热生成(如焦耳热)等的特殊问题。例如，可使用热质点单元MASS71模拟随温度变化的内部热生成。1.3热分析的类型ANSYS支持两种类型的热分析：1.稳态热分析确定在稳态的条件下的温度分布及其他热特性，稳态条件指热量随时间的变化可以忽略。2.瞬态热分析则计算在随时间变化的条件下，温度的分布和热特性。1.4耦合场分析ANSYS中可与热分析进行耦合的方式有热结构、热-电磁等。耦合场分

4、析可以使用ANSYS中的矩阵耦合单元，或者在独立的物理环境中使用序惯荷载耦合。有关耦合场分析的详细描述，请参阅ANSYS Coupled-Field Analysis Guide。1.5关于菜单路径和命令语法在本指南中，您将会看到相关的ANSYS命令及其等效的菜单路径。这些参考的命令仅仅包括命令名，因为并不总是需要指定所有的参数，而且不同的参数组合会有不同的作用。有关ANSYS命令的更多的叙述，请参考ANSYS Commands Reference。菜单路径将近可能完整得列出。对于多数情况，选择菜单就能够完成所需要的功能;但还有一些情况，选择文中所示菜单后会弹出一个菜单或是对话框，由此定义其他

5、的选项来执行一些特定的任务。第二章 基础知识2.1符号与单位项目国际单位英制单位ANSYS代号长度mft时间ss质量Kglbm温度oFTEMP力Nlbf能量（热量）JBTU功率（热流率）WBTU/secHEAT热流密度W/m2BTU/sec-ft2HFLUX生热速率W/m3BTU/sec-ft3HGEN导热系数W/m-BTU/sec-ft-oFKXX对流系数 W/m2-BTU/sec-ft2-oFHF密度Kg/m3lbm/ft3DENS比热J/Kg-BTU/lbm-oFC焓J/m3BTU/ft3ENTH2.2传热学经典理论回顾热分析遵循热力学第一定律，即能量守恒定律。对于一个封闭的系统（没有质

6、量的流入或流出）：式中：热量作功系统内能 系统动能 系统势能对大多数工程传热问题：；通常不考虑做功：，则；对于稳态热分析：，即流入的热量等于流出的热量；对于瞬态热分析：，即流入流出的热传递速率等于系统内能的变化。2.3热传递的方式2.3.1热传导热传导可以定义为完全接触的两个物体之间或一个物体的不同部分之间由于温度梯度而引起的内能的交换。热传导遵循傅立叶定律：，式中为热流密度（W/m2），为导热系数(W/m-)，负号表示热量流向温度降低的方向。2.3.2热对流热对流是指固体的表面与它周围接触的流体之间，由于温差的存在引起的热量的交换。热对流可以分为两类：自然对流和强制对流。热对流用牛顿冷却方程

7、来描述：，式中为对流换热系数（或称膜传热系数、给热系数、膜系数等）；为固体表面的温度，为周围流体的温度。2.3.3热辐射热辐射指物体发射电磁能，并被其它物体吸收转变为热的热量交换过程。物体温度越高，单位时间辐射的热量越多。热传导和热对流都需要有传热介质，而热辐射无须任何介质。实质上，在真空中的热辐射效率最高。在工程中通常考虑两个或两个以上物体之间的辐射，系统中每个物体同时辐射并吸收热量。它们之间的净热量传递可以用斯蒂芬波尔兹曼方程来计算：，式中为热流率，为辐射率（黑度），为斯蒂芬波尔兹曼常数，约为约为5.6710-8W/m2.K4，为辐射面1的面积，为由辐射面1到辐射面2的形状系数，为辐射面1

8、的绝对温度，为辐射面2的绝对温度，由上式可以看出，包含热辐射的热分析是高度非线性的。2.4稳态传热如果系统的净流滤为0，即流入系统的热量加上系统自身产生的热量等于流出系统的热量：，则系统热稳态。在稳态热分析中，任一节点的温度不随时间变化。稳态热分析的能量平衡方程为（以矩阵形式表示）：式中：为传导矩阵，包含热系数、对流系数及辐射和形状系数；为节点温度向量；为节点热流率向量，包括热生成；ANSYS利用模型几何差数、材料热性能参数以及所施加的边界条件，生成、及。2.5瞬态传热瞬态传热过程是指一个系统的加热或冷却过程。在这个过程中系统的温度、热流率、热边界条件以及系统内能随时间都有明显变化。根据能量守

9、恒原理，瞬态热平衡可以表达为（以矩阵形式表示）：式中：为传导矩阵，包含热系数、对流系数及辐射和形状系数；为比热矩阵，考虑系统内能的增加；为节点温度向量；为温度对时间的导数；为节点热流率向量，包括热生成；2.6线性与非线性如果有下列情况产生，则为非线性热分析：材料热性能随温度变化，如K(T)，C(T)等；边界条件随温度变化，如h(T)等；含有非线性单元；考虑辐射传热；非线性热分析的热平衡方程为：2.7边界条件和初始条件ANSYS热分析的边界条件或初始条件可分为七种：温度，热流率、热流密度、对流、辐射、绝热、生热。在本指南中，您将会看到相关的ANSYS命令及其等效的菜单路径。这些参考的命令仅仅包括

10、命令名，因为并不总是需要指定所有的参数，而且，不同的参数的组合会有不同的作用。有关ANSYS命令的更多的叙述，请参考ANSYS Commands Reference。菜单路径将近可能完整得列出，2.8热分析误差估计仅用于评估由于网格密度不够带来的误差；仅适用于SOLID或SHELL的热单元（只有一个温度自由度）；基于单元边界的热流密度的不连续；仅对一种材料、线性、稳态热分析有效；使用自适应网格划分可对误差进行控制。第三章 稳态热分析3.1稳态传热的定义ANSYS/Multiphysics，ANSYS/Mechanical，ANSYS/FLOTRAN和ANSYS/Professional这些产品

11、支持稳态热分析。稳态传热用于分析稳定的热载荷对系统或部件的影响。通常在进行瞬态热分析以前，进行稳态热分析用于确定初始温度分布。也可以在所有瞬态效应消失后，将稳态热分析作为瞬态热分析的最后一步进行分析。稳态热分析可以计算确定由于不随时间变化的热载荷引起的温度、热梯度、热流率、热流密度等参数。这些热载荷包括：对流辐射热流率热流密度（单位面积热流）热生成率（单位体积热流）固定温度的边界条件稳态热分析可用于材料属性固定不变的线性问题和材料性质随温度变化的非线性问题。事实上，大多数材料的热性能都随温度变化，因此在通常情况下，热分析都是非线性的。当然，如果在分析中考虑辐射，则分析也是非线性的。3.2热分析

12、的单元ANSYS和ANSYS/Professional中大约有40种单元有助于进行稳态分析。有关单元的详细描述请参考ANSYS Element Reference，该手册以单元编号来讲述单元，第一个单元是LINK1。单元名采用大写，所有的单元都可用于稳态和瞬态热分析。其中SOLID70单元还具有补偿在恒定速度场下由于传质导致的热流的功能。这些热分析单元如下：表3-1二维实体单元单元维数形状及特点自由度PLANE35二维六节点三角形单元温度（每个节点）PLANE55二维四节点四边形单元温度（每个节点）PLANE75二维四节点谐单元温度（每个节点）PLANE77二维八节点四边形单元温度（每个节点）

13、PLANE38二维八节点谐单元温度（每个节点）表3-2三维实体单元单元维数形状及特点自由度SOLID70三维八节点六面体单元温度（每个节点）SOLID87三维十节点四面体单元温度（每个节点）SOLID90三维二十节点六面体单元温度（每个节点）表3-3辐射连接单元单元维数形状及特点自由度LINK31二维或三维二节点线单元温度（每个节点）表3-4传导杆单元单元维数形状及特点自由度LINK32二维二节点线单元温度（每个节点）LINK33三维二节点线单元温度（每个节点）表3-5对流连接单元单元维数形状及特点自由度LINK34三维二节点线单元温度（每个节点）表3-6壳单元单元维数形状及特点自由度SHEL

14、L57三维四节点四边形单元温度（每个节点）表3-7耦合场单元单元维数形状及特点自由度PLANE13二维四节点热应力耦合单元温度、结构位移、电位、磁矢量位CONTACT48二维三节点热应力接触单元温度、结构位移CONTACT49三维热应力接触单元温度、结构位移FLUID116三维二或四节点热流单元温度、压力SOLID5三维八节点热应力和热电单元温度、结构位移、电位、磁标量位SOLID98三维十节点热应力和热电单元温度、结构位移、电位、磁矢量位PLANE67二维四节点热电单元温度、电位LINK68三维两节点热电单元温度、电位SOLID69三维八节点热电单元温度、电位SHELL157三维四节点热电单

15、元温度、电位表3-8特殊单元单元维数形状及特点自由度MASS71一维到三维一个节点的质量单元温度COMBINE37一维四节点控制单元温度、结构位移、转动、压力SURF151二维二到四节点面效应单元温度SURF152三维四到九节点面效应单元温度MATRIX50由包括在超单元中的单元类型决定没有固定形状的矩阵或辐射矩阵超单元由包括在超单元中的单元类型决定INFIN9二维二节点无限边界单元温度、磁矢量位INFIN47三维四节点无限边界单元温度、磁矢量位COMBINE14一维到三维两节点弹簧阻尼单元温度、结构位移、转动、压力COMBINE39一维两节点非线性弹簧单元温度、结构位移、转动、压力COMBI

16、NE40一维两节点组合单元温度、结构位移、转动、压力3.3热分析的基本过程ANSYS热分析包含如下三个主要步骤：前处理：建模求解：施加荷载并求解后处理：查看结果以下的内容将讲述如何执行上面的步骤。首先，对每一步的任务进行总体的介绍，然后通过一个管接处的稳态热分析的实例来引导读者如何按照GUI路径逐步完成一个稳态热分析。最后，本章提供了该实例等效的命令流文件。3.4建模建立一个模型的内容包括：首先为分析指定jobname和title；然后在前处理器（PREP7）中定义单元类型，单元实常数，材料属性以及建立几何实体。ANSYS Modeling and Meshing Guide中对本部分有详细说

17、明。对于热分析有：定义单元类型命令：ETGUI：Main MenuPreprocessorElement TypeAdd/Edit/Delete定义固定材料属性命令：MPGUI：Main MenuPreprocessorMaterial PropsMaterial ModelsThermal定义温度相关的材料属性，首先要定义温度表，然后定义对应的材料属性值。通过下面的方法定义温度表命令：MPTEMP或MPTEGN，然后定义对应的材料属性，使用MPDATAGUI：Main MenuPreprocessorMaterial Props Material ModelsThermal对于温度相关的对流

18、换热系数也是通过上述的GUI路径和命令来定义的。注意-如果以多项式的形式定义了与温度相关的膜系数，则在定义其它具有固定属性的材料之前，必须定义一个温度表。创建几何模型及划分划分网格的过程，请参阅ANSYS Modeling and Meshing Guide3.5施加荷载和求解在这一步骤中，必须指定所要进行的分析类型及其选项，对模型施加荷载，定义荷载选项，最后执行求解。3.5.1指定分析类型在这一步中，可以如下指定分析类型：GUI: Main MenuSolutionNew AnalysisSteady-state(static)命令：ANTYPE,STATIC,NEW如果是重新启动以前的分析

19、，比如，附加一个荷载。命令：ANTYPE,STATIC,rest。（条件是先前分析的jobname.ESAV、jobname.DB等文件是可以利用的）3.5.2施加荷载可以直接在实体模型（点、线、面、体）或有限元模型（节点和单元）上施加载荷和边界条件，这些载荷和边界条件可以是单值的，也可以是用表格或函数的方式来定义复杂的边界条件，详见ANSYS基本分析过程指南。可以定义以下五种热载荷：3.5.2.1恒定的温度(TEMP)通常作为自由度约束施加于温度已知的边界上。3.5.2.2热流率（HEAT）热流率作为节点集中载荷，主要用于线单元（如传导杆、辐射连接单元等）模型中，而这些线单元模型通常不能直接

20、施加对流和热流密度载荷。如果输入的值为正，表示热流流入节点，即单元获取热量。如果温度与热流率同时施加在一节点上，则温度约束条件优先。注意-如果在实体单元的某一节点上施加热流率，则此节点周围的单元应该密一些；特别是与该节点相连的单元的导热系数差别很大时，尤其要注意，不然可能会得到异常的温度值。因此，只要有可能，都应该使用热生成或热流密度边界条件，这些热荷载即使是在网格较为粗糙的时候都能得到较好的结果。3.5.2.3对流（CONV）对流边界条件作为面载施加于分析模型的外表面上，用于计算与模型周围流体介质的热交换，它仅可施加于实体和壳模型上。对于线单元模型，可以通过对流杆单元LINK34来定义对流。

21、3.5.2.4热流密度（HEAT）热流密度也是一种面载荷。当通过单位面积的热流率已知或通过FLOTRAN CFD的计算可得到时，可以在模型相应的外表面或表面效应单元上施加热流密度。如果输入的值为正，表示热流流入单元。热流密度也仅适用于实体和壳单元。单元的表面可以施加热流密度也可以施加对流，但ANSYS仅读取最后施加的面载进行计算。3.5.2.5热生成率（HGEN）热生成率作为体载施加于单元上，可以模拟单元内的热生成，比如化学反应生热或电流生热。它的单位是单位体积的热流率。下表总结了在热分析中的载荷类型：表3-9 热荷载类型载荷类型类别命令族GUI 路径温度 (TEMP)约束DMain Menu

22、Solution-Loads-Apply -Thermal-Temperature热流率 (HEAT)力FMain MenuSolution-Loads-Apply -Thermal-Heat Flow对流 (CONV), 热流密度 (HFLUX)面载荷SFMain MenuSolution-Loads-Apply -Thermal-ConvectionMain MenuSolution-Loads-Apply -Thermal-Heat Flux热生成率 (HGEN)体载荷BFMain MenuSolution-Loads-Apply -Thermal-Heat Generat下表详细列出了

23、热分析中用于施加载荷，删除载荷，对载荷进行操作、列表的所有命令：表3-10 热荷载相关的命令载荷类型实体或有限元模型实体施加删除列表显示运算设置温度实体模型关键点DKDKDELEDKLISTDTRAN-有限元模型节点DDDELEDLISTDSCALEDCUMTUNIF热流率实体模型关键点FKFKDELEFKLISTFTRAN-有限元模型节点FFDELEFLISTFSCALEFCUM对流, 热流密度实体模型线SFLSFLDELESFLLISTSFTRANSFGRAD实体模型面SFASFADELESFALISTSFTRANSFGRAD有限元模型节点SFSFDELESFLISTSFSCALESFGR

24、ADSFCUM有限元模型单元SFESFEDELESFELISTSFSCALESFBEAMSFCUMSFFUNSFGRAD生热率实体模型关键点BFKBFKDELEBFKLISTBFTRAN-实体模型线BFLBFLDELEBFLLISTBFTRAN-实体模型面BFABFADELEBFALISTBFTRAN-实体模型体BFVBFVDELEBFVLISTBFTRAN-有限元模型节点BFBFDELEBFLISTBFSCALEBFCUM单元BFEBFEDELEBFELISTBFSCALEBFCUM3.5.3采用表格和函数边界条件除了一般的使用表格来定义边界条件的方法，本节讨论热分析中特有的一些问题。关于定

25、义表参数的详细叙述，请参考ANSYS APDL Programmers Guide。本节内容对单元类型没有特别的限制。下表列出了热分析中能够用于每一种边界条件的自变量：表3-11荷载边界条件及其自变量热边界条件命令族自变量固定温度DTIME, X, Y, Z热流FTIME, X, Y, Z, TEMP对流换热系数 (对流)SFTIME,X, Y, Z, TEMP, VELOCITY环境温度 (对流)SFTIME, X, Y, Z热流密度SFTIME, X, Y, Z, TEMP热生成BFTIME, X, Y, Z, TEMP流体单元(FLUID116) 边界条件流率SFETIME压力DTIME

26、, X, Y, Z后面有一个例题详细介绍在一个稳态热分析中如何采用表格边界条件。为了使用更加灵活的热传导系数，可以使用函数的方式来定义边界条件。有关这种用法的详细说明，可以参考ANSYS Basic Analysis Procedures Guide。除了上述自变量外，函数边界条件还可用下面的参数作为函数的自变量：表面温度（TS）（SURF151、SURF152单元的表面温度）密度（）（材料属性DENS）比热（材料属性C）导热率（材料属性kxx）导热率（材料属性kyy）导热率（材料属性kzz）粘度（材料属性）辐射率（材料属性）3.5.4定义载荷步选项对于一个热分析，可以确定通用选项、非线性选项

27、以及输出控制。下表列出了热分析中可能用到的载荷步选项：表3-12分析中的载荷步选项选项命令GUI 路径通用选项时间TIMEMain MenuSolution-Load Step Opts-Time/FrequencTime-Time Step时间步数NSUBSTMain MenuSolution -Load Step Opts-Time/FrequencTime and Substps时间步长DELTIMMainMenuSolution -Load Step Opts-Time /FrequencTime-Time Step阶跃或斜坡加载KBCMainMenuSolution -Load St

28、ep Opts-Time /FrequencTime -Time Step非线性选项最大平衡迭代数NEQITMain MenuSolution -Load Step Opts-NonlinearEquilibrium Iter自动时间步长AUTOTSMain MenuSolution -Load Step Opts-Time/FrequencTime-Time Step收敛容差CNVTOLMain MenuSolution -Load Step Opts-NonlinearConvergence Crit求解中断选项NCNVMainMenuSolution -Load Step Opts -NonlinearCriteria to Stop