第6章ANSYS热分析.ppt

《第6章ANSYS热分析.ppt》由会员分享，可在线阅读，更多相关《第6章ANSYS热分析.ppt（49页珍藏版）》请在三一文库上搜索。

1、第6章 ANSYS热分析,李达 西南交通大学材料学院,1 简介 一、热分析的目的 二、ANSYS的热分析 三、ANSYS热分析分类 四、耦合分析 2 基础知识 一、符号与单位 二、传热学经典理论回顾 三、热传递的方式 四、稳态传热 五、瞬态传热 六、线性与非线性 七、边界条件、初始条件 八、热分析误差估计,主要内容,3 稳态传热分析 (1)、稳态传热的定义 (2)、热分析的单元 (3)、ANSYS稳态热分析的基本过程 练习 4 瞬态传热分析 (1)、瞬态传热分析的定义 (2)、瞬态热分析的单元及命令 (3)、ANSYS瞬态热分析的主要步骤 1、建模 2、加载求解 3、后处理 (4)、相变问题

2、练习,6-1. 热分析的目的 热分析用于计算一个系统或部件的温度分布及其它热物理参数，如热量的获取或损失、热梯度、热流密度（热通量等。 热分析在许多工程应用中扮演重要角色，如内燃机、换热器、管路系统、电子元件等。,6-2. ANSYS的热分析 在ANSYS/Multiphysics、ANSYS/Mechanical、ANSYS/Thermal、ANSYS/FLOTRAN、ANSYS/ED五种产品中包含热分析功能。 ANSYS热分析基于能量守恒原理的热平衡方程，用有限元法计算各节点的温度，并导出其它热物理参数。 ANSYS热分析包括热传导、热对流及热辐射三种热传递方式。此外，还可以分析相变、有内

3、热源、接触热阻等问题。,6-3. ANSYS的热分析分类 稳态传热：系统的温度场不随时间变化 瞬态传热：系统的温度场随时间明显变化 9-4.与热有关的耦合分析 热结构耦合 热流体耦合 热电耦合 热磁耦合 热电磁结构耦合等,热分析 基 础 知 识,1、符号与单位 2、传热学经典理论回顾 3、热传递的方式 4、稳态传热 5、瞬态传热 6、线性与非线性 7、边界条件、初始条件 8、热分析误差估计,1、符号与单位,2、传热学经典理论回顾,1、热传导 热传导可以定义为完全接触的两个物体之间或一个物体的不同部分之间由于温度梯度而引起的内能的交换。热传导遵循付里叶定律：qn=-k*(dT/dx)，式中qn为

4、热流密度（W/m2），k为导热系数（W/m-），“-”表示热量流向温度降低的方向。 2、热对流 热对流是指固体的表面与它周围接触的流体之间，由于温差的存在引起的热量的交换。热对流可以分为两类：自然对流和强制对流。热对流用牛顿冷却方程来描述： qn= h*(TS-TB)，式中h为对流换热系数（或称膜传热系数、给热系数、膜系数等）， TS为固体表面的温度， TB为周围流体的温度。,3、热传递的方式,3、热辐射 热辐射指物体发射电磁能，并被其它物体吸收转变为热的热量交换过程。物体温度越高，单位时间辐射的热量越多。热传导和热对流都需要有传热介质，而热辐射无须任何介质。实质上，在真空中的热辐射效率最高。

5、 在工程中通常考虑两个或两个以上物体之间的辐射，系统中每个物体同时辐射并吸收热量。它们之间的净热量传递可以用斯蒂芬波尔兹曼方程来计算：q=A1F12(T14-T24)，式中q为热流率， 为辐射率（黑度）， 为斯蒂芬波尔兹曼常数，约为5.6710-8W/m2.K4，A1为辐射面1的面积，F12为由辐射面1到辐射面2的形状系数，T1为辐射面1的绝对温度，T2为辐射面2的绝对温度。由上式可以看出，包含热辐射的热分析是高度非线性的。,如果系统的净热流率为，即流入系统的热量加上系统自身产生的热量等于流出系统的热量：q流入+q生成-q流出=0，则系统处于热稳态。在稳态热分析中任一节点的温度不随时间变化。稳

6、态热分析的能量平衡方程为（以矩阵形式表示） KT=Q 式中： K为传导矩阵，包含导热系数、对流系数及辐射率和形状系数； T为节点温度向量； Q为节点热流率向量，包含热生成； ANSYS利用模型几何参数、材料热性能参数以及所施加的边界条件，生成K 、 T以及Q 。,4、稳态传热,瞬态传热过程是指一个系统的加热或冷却过程。在这个过程中系统的温度、热流率、热边界条件以及系统内能随时间都有明显变化。根据能量守恒原理，瞬态热平衡可以表达为（以矩阵形式表示）： C +KT=Q 式中: K为传导矩阵，包含导热系数、对流系数及辐射率和形状系数； C为比热矩阵,考虑系统内能的增加; T为节点温度向量； 为温度对

7、时间的导数; Q为节点热流率向量，包含热生成。,6、瞬态传热,6、线性与非线性,如果有下列情况产生，则为非线性热分析： 、材料热性能随温度变化，如K(T),C(T)等； 、边界条件随温度变化，如h(T)等； 、含有非线性单元； 、考虑辐射传热 非线性热分析的热平衡矩阵方程为： C(T) +K(T)T=Q (T),7、边界条件、初始条件,ANSYS热分析的边界条件或初始条件可分为七种： 温度; 热流率; 热流密度; 对流; 辐射; 绝热; 生热。, 仅用于评估由于网格密度不够带来的误差； 仅适用于SOLID或SHELL的热单元（只有温度一个自由度）； 基于单元边界的热流密度的不连续； 仅对一种材

8、料、线性、稳态热分析有效； 使用自适应网格划分可以对误差进行控制。,8、热分析误差估计,稳态传热分析,目的：通过学习，介绍稳态传热分析的知 识，以助于大家了解以后的稳态传热分析，并 能对简单的问题进行求解。 1、稳态传热的定义 2、热分析的单元 3、 ANSYS稳态热分析的基本过程 练习,1 稳态传热定义 稳态传热用于分析稳定的热载荷对系统或部件的影响。通常在进行瞬态热分析以前，进行稳态热分析用于确定初始温度分布。 稳态热分析可以通过有限元计算确定由于稳定的热载荷引起的温度、热梯度、热流率、热流密度等参数,热分析涉及到的单元有大约40种，其中纯粹用于热分析的有14种： 线性： LINK32 两

9、维二节点热传导单元 LINK33 三维二节点热传导单元 LINK34 二节点热对流单元 LINK31 二节点热辐射单元 二维实体：PLANE55 四节点四边形单元 PLANE77 八节点四边形单元 PLANE35 三节点三角形单元 PLANE75 四节点轴对称单元 PLANE78 八节点轴对称单元 三维实体SOLID87 六节点四面体单元 SOLID70 八节点六面体单元 SOLID90 二十节点六面体单元 壳 SHELL57 四节点 点 MASS71 有关单元的详细解释，请参阅ANSYS Element Reference Guide,2 热分析单元,ANSYS热分析可分为三个步骤： 前处理

10、： 建模 求解： 施加载荷计算 后处理： 查看结果,3 热分析基本过程,建模,、确定jobname、title、unit； 、进入PREP7前处理，定义单元类型，设定单元选项； 、定义单元实常数； 、定义材料热性能参数，对于稳态传热，一般只需定义导热系数，它可以是恒定的，也可以随温度变化； 、创建几何模型并划分网格，请参阅ANSYS Modeling and Meshing Guide。,加载计算,、定义分析类型 如果进行新的热分析： Command: ANTYPE, STATIC, NEW GUI: Main menuSolution-Analysis Type-New AnalysisSt

11、eady-state 如果继续上一次分析，比如增加边界条件等： Command: ANTYPE, STATIC, REST GUI: Main menuSolutionAnalysis Type-Restart,、施加载荷 可以直接在实体模型或单元模型上施加五种载荷(边界条件) ： a、恒定的温度 通常作为自由度约束施加于温度已知的边界上。 Command Family: D GUI：Main MenuSolution-Loads-Apply-Thermal-Temperature,b、热流率 热流率作为节点集中载荷，主要用于线单元模型中(通常线单元模型不能施加对流或热流密度载荷)，如果输入的

12、值为正，代表热流流入节点，即单元获取热量。如果温度与热流率同时施加在一节点上则ANSYS读取温度值进行计算。 注意：如果在实体单元的某一节点上施加热流率，则此节点周围的单元要密一些，在两种导热系数差别很大的两个单元的公共节点上施加热流率时，尤其要注意。此外，尽可能使用热生成或热流密度边界条件，这样结果会更精确些。 Command Family: F GUI：Main MenuSolution-Loads-Apply-Thermal-Heat Flow,c、对流 对流边界条件作为面载施加于实体的外表面，计算与流体的热交换，它仅可施加于实体和壳模型上，对于线模型，可以通过对流线单元LINK34考虑

13、对流。 Command Family: SF GUI：Main MenuSolution-Loads-Apply-Thermal-Convection,d、热流密度 热流密度也是一种面载。当通过单位面积的热流率已知或通过FLOTRAN CFD计算得到时，可以在模型相应的外表面施加热流密度。如果输入的值为正，代表热流流入单元。热流密度也仅适用于实体和壳单元。热流密度与对流可以施加在同一外表面，但ANSYS仅读取最后施加的面载进行计算。 Command Family: F GUI：Main MenuSolution-Loads-Apply-Thermal-Heat Flux,e、生热率 生热率作为

14、体载施加于单元上，可以模拟化学反应生热或电流生热。它的单位是单位体积的热流率。 Command Family: BF GUI：Main MenuSolution-Loads-Apply-Thermal-Heat Generat,定载荷步选项 对于一个热分析，可以确定普通选项、非线性选项以及输出控制。 a. 普通选项 时间选项：虽然对于稳态热分析，时间选项并没有实际的物理意义，但它提供了一个方便的 设置载荷步和载荷子步的 方法。 Command: TIME GUI: Main MenuSolution-Load Step Opts- Time/Frequenc Time-Time Step/ T

15、ime and Substps,每载荷步中子步的数量或时间步大小：对于非线性分析，每一载荷步需要多个子步。 Command: NSUBST GUI: Main MenuSolution-Load Step Opts- Time/Frequenc Time and Substps Command: DELTIM GUI: Main MenuSolution-Load Step Opts- Time/Frequenc Time-Time Step,递进或阶越选项：如果定义阶越(stepped)选项，载荷值在这个载荷步内保持不变；如果为递进(ramped)选项，则载荷值由上一载荷步值到本载荷步值随每

16、一子步线性变化。 Command: KBC GUI: Main MenuSolution-Load Step Opts-Time/FrequencTime-Time Step/Time and Substps,定载荷步选项 b. 非线性选项 迭代次数：本选项设置每一子步允许的最多的迭代次数。默认值为25，对大数热分析问题足够。 Command: NEQIT GUI: Main MenuSolution-Load Step Opts-NolinearEquilibrium Iter,自动时间步长: 对于非线性问题，可以自动设定子步间载荷的增长，保证求解的稳定性和准确性。 Command: AUT

17、OTS GUI: Main MenuSolution-Load Step Opts-Time/FrequencTime-Time Step/Time and Substps,收敛误差：可根据温度、热流率等检验热分析的收敛性。 Command: CNVTOL GUI: Main MenuSolution-Load Step Opts-NolinearConvergence Crit,求解结束选项：如果在规定的迭代次数内，达不到收敛，ANSYS可以停止求解或到下一载荷步继续求解。 Command: NCNV GUI: Main MenuSolution-Load Step Opts-Nolinea

18、rCriteria to Stop,线性搜索：设置本选项可使ANSYS用Newton-Raphson方法进行线性搜索。 Command: LNSRCH GUI: Main MenuSolution-Load Step Opts-NolinearLine Search 预测矫正：本选项可激活每一子步第一次迭代对自由度求解的预测矫正。 Command: PRED GUI: Main MenuSolution-Load Step Opts-NolinearPredictor,定载荷步选项 c. 输出控制 控制打印输出：本选项可将任何结果数据输出到*.out 文件中。 Command: OUTPR G

19、UI: Main MenuSolution-Load Step Opts-Output CtrlsSolu Printout,c. 输出控制 控制结果文件：控制*.rth的内容。 Command: OUTRES GUI: Main MenuSolution-Load Step Opts-Output CtrlsDB/Results File,确定分析选项 a. Newton-Raphson选项（仅对非线性分析有用） Command: NROPT GUI: Main MenuSolutionAnalysis Options,b. 选择求解器：可选择如下求解器中一个进行求解： Frontal so

20、lver(默认） Jacobi Conjugate Gradient(JCG) solver JCG out-of-memory solver Incomplete Cholesky Conjugate Gradient(ICCG) solver Pre-Conditioned Conjugate Gradient Solver(PCG) Iterative(automatic solver selection option) Command: EQSLV GUI: Main MenuSolutionAnalysis Options 注意：热分析可选用Iterative选项进行快速求解，但如下

21、情况除外： 热分析包含SURF19或SURF22或超单元； 热辐射分析； 相变分析 需要restart an analysis,c. 确定绝对零度：在进行热辐射分析时，要将目前的温度值换算为绝对温度。如果使用的温度单位是摄氏度，此值应设定为273；如果使用的是华氏度，则为460。 Command: TOFFST GUI: Main MenuSolutionAnalysis Options,保存模型: 点击ANSYS工具条SAVE_DB。 求解 Command: SOLVE GUI: Main MenuSolutionCurrent LS ANSYS将热分析的结果写入*.rth文件中，它包含如下

22、 数据： 基本数据： 节点温度 导出数据： 节点及单元的热流密度 节点及单元的热梯度 单元热流率 节点的反作用热流率 其它,对于稳态热分析，可以使用POST1进行后处理，关于后处理的完整描述，可参阅ANSYS Basic Analysis Procedures Guide。 进入POST1后，读入载荷步和子步： Command: SET GUI: Main MenuGeneral Postproc-Read Results-By Load Step,可以通过如下三种方式查看结果： 彩色云图显示 Command: PLNSOL, PLESOL, PLETAB等 GUI: Main MenuGen

23、eral Postproc Plot ResultsNodal Solu, Element Solu, Elem Table,矢量图显示 Command: PLVECT GUI: Main MenuGeneral PostprocPlot ResultsPre-defined or Userdefined,列表显示 Command: PRNSOL, PRESOL, PRRSOL等 GUI: Main MenuGeneral PostprocList ResultsNodal Solu, Element Solu, Reaction Solu 详细过程请参阅ANSYS Basic Analysis Procedures Guide。,练习,电线生热分析： 电热丝的半径r=1.015mm，电阻率=8010-6*cm，热导率=19.03W/(m.k)，稳态时通过电热丝的电流I=150A，几何模型如图所示，简化为轴对称，取轴线长度1.5mm，求中心线温度比表面温度高多少。,问题分析： 1 可选用平面热分析PLANE55单元和SOLID单元进行分析。 2 将电流产生的热能作为热生成体载荷施加到电热丝上，热能Q的计算,