热分析(ansys教程)..

热分析热分析的目的热分析用于计算一个系统或部件的温度分布及其它热物理参数，如热量的获取或损失、热梯度、热流密度(热通量)等热分析在许多工程应用中扮演重要角色，如内燃机、涡轮机、换热器、管路系统、电子元件等ANSYS的热分析在ANSYS/Multiphysics、ANSYS/Mechanical、ANSYS/Thermal、ANSYS/FLOTRAN、ANSYS/ED五种产品中包含热分析功能ANSYS热分析基于能量守恒原理的热平衡方程，用有限元法计算各节点的温度，并导出其它热物理参数ANSYS热分析包括热传导、热对流及热辐射三种热传递方式。此外，还可以分析相变、有内热源、接触热阻等问题ANSYS的热分析分类ANSYS的热分析分类稳态传热：系统的温度场不随时间变化瞬态传热：系统的温度场随时间明显变化与热有关的耦合分析热－结构耦合热－流体耦合热－电耦合热－磁耦合热－电－磁－结构耦合等热分析的符号与单位项目国际单位英制单位ANSYS代号长度mft[英尺]时间ss质量Kglbm[磅质量]温度℃oF力Nlbf能量(热量)JBTU[英制热单位]功率(热流率)WBTU/sec热流密度W/m2BTU/sec-ft2生热速率W/m3BTU/sec-ft3导热系数W/m-℃BTU/sec-ft-oFKXX对流系数W/m2-℃BTU/sec-ft2-oFHF密度Kg/m3lbm/ft3DENS比热J/Kg-℃BTU/lbm-oFC焓J/m3BTU/ft3ENTH表征物体吸收的热量，为一个体系的内能与体系的体积和外界施加于体系的压强的乘积之和传热学经典理论回顾热传递的方式1、热传导热传导可以定义为完全接触的两个物体之间或一个物体的不同部分之间由于温度梯度而引起的内能的交换。热传导遵循付里叶定律：qn=-k*(dT/dx)，式中qn为热流密度（W/m2），k为导热系数(W/m-℃)，“-”表示热量流向温度降低的方向。2、热对流热对流是指固体的表面与它周围接触的流体之间，由于温差的存在引起的热量的交换。热对流可以分为两类：自然对流和强制对流。热对流用牛顿冷却方程来描述：qn=h*(TS-TB)，式中h为对流换热系数(或称膜传热系数、给热系数、膜系数等)，TS为固体表面的温度，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}式中：[K]为传导矩阵，包含导热系数、对流系数及辐射率和形状系数；{T}为节点温度向量；{Q}为节点热流率向量，包含热生成；ANSYS利用模型几何参数、材料热性能参数以及所施加的边界条件，生成[K]、{T}以及{Q}。瞬态传热瞬态传热过程是指一个系统的加热或冷却过程。在这个过程中系统的温度、热流率、热边界条件以及系统内能随时间都有明显变化。根据能量守恒原理，瞬态热平衡可以表达为(以矩阵形式表示)：[C]{}+[K]{T}={Q}式中:[K]为传导矩阵，包含导热系数、对流系数及辐射率和形状系数；[C]为比热矩阵,考虑系统内能的增加;{T}为节点温度向量；{}为温度对时间的导数;{Q}为节点热流率向量，包含热生成。TT线性与非线性如果有下列情况产生，则为非线性热分析：①材料热性能随温度变化，如K(T),C(T)等；②边界条件随温度变化，如h(T)等；③含有非线性单元；④考虑辐射传热非线性热分析的热平衡矩阵方程为：[C(T)]{}+[K(T)]{T}={Q(T)}T边界条件、初始条件ANSYS热分析的边界条件或初始条件可分为七种：温度：模型区温度已知热流率：热流率已知的点对流：表面的热传递给周围的流体通过对流。输入对流换热系数h和环境流体的平均温度Tb热辐射：通过辐射产生热传递的面.输入辐射系数，Stefan-Boltzmann常数，“空间节点”的温度作为可选项输入绝热面：“完全绝热”面，该面上不发生热传递热通量：单位面积上的热流率已知的面热生成率：体的生热率已知的区域热分析误差估计仅用于评估由于网格密度不够带来的误差；仅适用于SOLID或SHELL的热单元(只有温度一个自由度)；基于单元边界的热流密度的不连续；仅对一种材料、线性、稳态热分析有效；使用自适应网格划分可以对误差进行控制。稳态传热分析稳态传热用于分析稳定的热载荷对系统或部件的影响。通常在进行瞬态热分析以前，进行稳态热分析用于确定初始温度分布。稳态热分析可以通过有限元计算确定由于稳定的热载荷引起的温度、热梯度、热流率、热流密度等参数热分析的单元热分析涉及到的单元有大约40种，其中纯粹用于热分析的有14种：线性：LINK32两维二节点热传导单元LINK33三维二节点热传导单元LINK34二节点热对流单元LINK31二节点热辐射单元二维实体：PLANE55四节点四边形单元PLANE77八节点四边形单元PLANE35三节点三角形单元PLANE75四节点轴对称单元PLANE78八节点轴对称单元热分析的单元(续)三维实体：SOLID87六节点四面体单元SOLID70八节点六面体单元SOLID90二十节点六面体单元壳：SHELL57四节点点：MASS71稳态热分析的基本过程ANSYS热分析可分为三个步骤：•前处理：建模•求解：施加载荷计算•后处理：查看结果建模①确定jobname、title、unit；②进入PREP7前处理，定义单元类型，设定单元选项；③定义单元实常数；④定义材料热性能参数，对于稳态传热，一般只需定义导热系数，它可以是恒定的，也可以随温度变化；⑤创建几何模型并划分网格。几何尺寸(模型)既可用ANSYS建立模型，也可用其它方法建好模型后导入模型建好后，以上两种建模方法的具体过程将不再显示划分网格首先定义单元属性:单元类型,实常数,材料属性.单元类型下表给出了常用的热单元类型每个结点只有一个自由度:温度2-DSolid3-DSolid3-DShellLineElementsLinearPLANE55SOLID70SHELL57LINK31,32,33,34QuadraticPLANE77PLANE35SOLID90SOLID87划分网格(续)材料属性必须输入导热系数,KXX如果施加了内部热生成率，则需指定比热(C)ANSYS提供的材料库(/ansys57/matlib)包括几种常用材料的结构属性和热属性,但是建议用户创建、使用自己的材料库把优先设置为“热分析”，使材料模型图形用户界面只显示材料的热属性实常数主要应用于壳单元和线单元划分网格(续)划分网格存储数据文件使用MeshTool划分网格，使用缺省的智能网格划分级别6可以生成很好的初始网格至此完成前处理，下面开始求解施加载荷计算①定义分析类型如果进行新的热分析：Command:ANTYPE,STATIC,NEWGUI:MainmenuSolution-AnalysisType-NewAnalysisSteady-state如果继续上一次分析，比如增加边界条件等：Command:ANTYPE,STATIC,RESTGUI:MainmenuSolutionAnalysisType-Restart施加载荷计算(续)②施加载荷可以直接在实体模型或单元模型上施加五种载荷(边界条件)：a、恒定的温度通常作为自由度约束施加于温度已知的边界上。CommandFamily:DGUI：MainMenuSolution-Loads-Apply-Thermal-Temperature施加载荷计算(续)b、热流率热流率作为节点集中载荷，主要用于线单元模型中(通常线单元模型不能施加对流或热流密度载荷)，如果输入的值为正，代表热流流入节点，即单元获取热量。如果温度与热流率同时施加在一节点上则ANSYS读取温度值进行计算。注意：如果在实体单元的某一节点上施加热流率，则此节点周围的单元要密一些，在两种导热系数差别很大的两个单元的公共节点上施加热流率时，尤其要注意。此外，尽可能使用热生成或热流密度边界条件，这样结果会更精确些。CommandFamily:FGUI：MainMenuSolution-Loads-Apply-Thermal-HeatFlow施加载荷计算(续)c、对流对流边界条件作为面载施加于实体的外表面，计算与流体的热交换，它仅可施加于实体和壳模型上，对于线模型，可以通过对流线单元LINK34考虑对流。CommandFamily:SFGUI：MainMenuSolution-Loads-Apply-Thermal-Convection施加载荷计算(续)d、热流密度热流密度也是一种面载。当通过单位面积的热流率已知或通过FLOTRANCFD计算得到时，可以在模型相应的外表面施加热流密度。如果输入的值为正，代表热流流入单元。热流密度也仅适用于实体和壳单元。热流密度与对流可以施加在同一外表面，但ANSYS仅读取最后施加的面载进行计算。CommandFamily:FGUI：MainMenuSolution-Loads-Apply-Thermal-HeatFlux施加载荷计算(续)e、生热率生热率作为体载施加于单元上，可以模拟化学反应生热或电流生热。它的单位是单位体积的热流率。CommandFamily:BFGUI：MainMenuSolution-Loads-Apply-Thermal-HeatGenerat施加载荷计算(续)③定载荷步选项对于一个热分析，可以确定普通选项、非线性选项以及输出控制。a.普通选项·时间选项：虽然对于稳态热分析，时间选项并没有实际的物理意义，但它提供了一个方便的设置载荷步和载荷子步的方法。Command:TIMEGUI:MainMenuSolution-LoadStepOpts-Time/FrequencTime-TimeStep/TimeandSubstps施加载荷计算(续)·每载荷步中子步的数量或时间步大小：对于非线性分析，每一载荷步需要多个子步。Command:NSUBSTGUI:MainMenuSolution-LoadStepOpts-Time/FrequencTimeandSubstpsCommand:DELTIMGUI:MainMenuSolution-LoadStepOpts-Time/FrequencTime-TimeStep施加载荷计算(续)·递进或阶越选项：如果定义阶越(stepped)选项，载荷值在这个载荷步内保持不变；如果为递进(ramped)选项，则载荷值由上一载荷步值到本载荷步值随每一子步线性变化。Command:KBCGUI:MainMenuSolution-LoadStepOpts-Time/FrequencTime-TimeStep/TimeandSubstps施加载荷计算(续)b.非线性选项·迭代次数：本选项设置每一子步允许的最多的迭代次数。默认值为25，对大数热分析问题足够。Command:NEQITGUI:MainMenu