Ansys热分析教程

第一章介绍先决条件•热分析讲座的先决条件包括：–完成ANSYS第一部分的培训课程或具备相应的ANSYS经验。–理工科学士或具有相等程度的热传递知识。综述：本培训教程相关章节的主要内容•第1章－介绍•第2章－基本概念–术语–符号和单位–热传递的类型–能量守恒定律(热力学第一定律)–瞬态热传导的控制微分方程–有限元方法–有限元热分析的基本特征–如何使热传递分析包括非线性？–什么时候需要定义比热和密度–与结构分析的比较–单元描述：功能和限制–例题1－基本的热传递分析章节内容概述•第3章－稳态热传导–稳态热传递的控制方程–热边界条件类型–热分析模板–GUI和ANSYS命令–借助例题详细讲解分析过程的每一步(带散热片的钢管热分析)–前处理—建立模型–求解分析–后处理－“后处理101”–例题2－晶体管的稳态热分析章节内容概述•第4章－非线性分析的一些特殊考虑–时间，载荷步，子步和平衡迭代–收敛准则–初始温度–阶跃或渐变载荷–其他非线性选项–输出控制选项–监测/查看非线性分析–例题3－晶体管的非线性热分析•第5章－瞬态分析–控制方程–与稳态分析的区别–查看瞬态结果–例题4－晶体管的瞬态热分析章节内容概述•第6章－随时间和空间变化的复杂边界条件–表格化的热边界条件(载荷)–基本变量–用户定义的独立变量–典型应用–优点和缺点–定义关键时间和输出时间–ANSYS表格和数组的回顾–函数工具–例题5－有表格化载荷的瞬态练习章节内容概述•第七章－对流选项和简单的热/流单元–对流作为面载荷施加–对流连接单元–表面效应单元–模拟接触热阻–1D热/流单元–用户对流子程序–例题6－低压燃气轮机匣的热分析章节内容概述•第8章－热辐射分析–辐射概念的回顾–基本定义–辐射模拟的可选方法–辐射矩阵–辐射通量密度(Radiation)法–辐射分析例题－使用辐射矩阵(隐藏和非隐藏方法)进行散热器分析–例题7-使用辐射通量法计算辐射传热•第9章－相变分析–基本概念/术语–在ANSYS中求解相变–相变例题－飞轮铸造分析章节内容概述•第10章－耦合场分析(重点是热应力)–直接和序列耦合场–物理环境–热应力分析–例题8－带散热片的钢管的热应力分析•第11章－热接触–基本概念–热传导、生热、热对流、热辐射、摩擦生热–变化热接触系数的传导问题–例题9－导热系数是接触压力的函数的接触传热•例题2的附录分步讲述分析过程相关的讲座和培训•ANSYS介绍，第2部分•CFD(计算流体动力学)分析•CFD高级分析•电子封装其它信息来源•ANSYS在线文档–例题手册，校验手册，热分析指南，耦合场分析指南•ANSYS自学例题•ANSYS参考论文和白皮书•ANSYS新闻和“AnalysisSolutions”杂志•其它ANSYS课程的培训手册•ANSYS会议论文集•ANSYS网站：(中文)第2章基本概念热传导基础符号下列符号在全文中的意义如下：tTchKQqqEtimetemperaturedensityspecificheatfilmcoefficientemissivityStefan-Boltzmannconstantthermalconductivityheatflow(rate)heatfluxinternalheatgeneration/volumeenergy*f热传导基础ANSYS中标准单位(英制)•温度•热流量•热传导率•密度•比热•对流换热系数•热流•温度梯度•内部热生成•ºF•BTU/hr•BTU/(hr-in-ºF)•lbm/(in3)•BTU/(lbm-ºF)•BTU/(hr-in2-ºF)•BTU/(hr-in2)•ºF/in•BTU/(hr-in3)注：对于结构热容量,密度/Gc和比热*Gc经常使用该单位。其中Gc为重力加速度=386.4(lbm-in)/(lbf-sec2)热传导基础ANSYS中标准单位(国际单位制)•温度•热流量•热传导率•密度•比热•对流换热系数•热流•温度梯度•内部热生成•ºC(orK)•Watts•Watts/(m-ºC)•Kg/(m3)•(Watt-sec)/(Kg-ºC)•Watt/(m2-ºC)•Watt/(m2)•ºC/m•Watt/(m3)注：对于结构热容量,密度/Gc和比热*Gc经常使用该单位。其中Gc为重力加速度=9.8(m/sec2)热传导基础ANSYS中的标准单位–汇总热传导基础热传递的类型•热传递有三种基本类型:–传导-两个良好接触的物体之间的能量交换或一个物体内由于温度梯度引起的内部能量交换。–对流-在物体和周围介质之间发生的热交换。–辐射-一个物体或两个物体之间通过电磁波进行的能量交换。•在绝大多数情况下，我们分析的热传导问题都带有对流和/或辐射边界条件。热传导基础传导•传导的热流由传导的傅立叶定律决定：•负号表示热沿梯度的反向流动(例如,热从高温部分流向低温部份).Tnq*dTdn热传导基础对流•对流的热流由冷却的牛顿准则得出:•对流一般作为面边界条件施加TBTs热传导基础辐射•从平面i到平面j的辐射热流由施蒂芬-玻斯曼定律得出：•在ANSYS中将辐射按平面现象处理(即：体都假设为不透明的)。ij热传导基础热力学第一定律•能量守恒要求系统的能量改变与系统边界处传递的热和功数值相等。•能量守恒在一个短的时间增量下可以表示为方程形式：•将其应用到一个微元体上，就可以得到热传导的控制微分方程。热传导基础控制微分方程•热传导的控制微分方程热传导基础有限元方法•将控制微分方程转化为等效的积分形式(参阅ANSYS理论手册第6.1节)。cTTtvLTLTDLTdvolTqdSThTTdSTqdvolvolLxyzqhTqTTvolSfBSvolTfBFHGIKJFHGIKJLNMOQPzzzzlqlqlqlqchlq()()()()()()**2323where=volumeoftheelement=heatflux,filmcoefficient,bulkfluidtemp.heatgenerationperunitvolumeT=anallowablevirtualtemperaturesurfaceswithappliedfluxsurfaceswithappliedconvectionSS23热传导基础有限元方法(续)•将区域分解(也称“划分”)为简单的形状;2-D模型中的四边形和/或三角形,3-D模型中的四面体，金字塔形或六面体。1热传导基础有限元方法(续)•假设单元内温度变化可以用多项式表示。一般情况下，根据单元类型的不同，应当包含不同的一次项,平方和混合的立方项。多项式假设保证了温度在单元内部和单元边界上都是连续的。•写出以单元节点温度为未知数的多项式：热传导基础有限元热分析的基本特点•求解连续性–温度在一个单元内和单元边界上是连续的(即，单值的)；–温度剃度和热通量在一个单元内是连续的，在单元边界上是不连续的。•能量平衡在每个节点上都能够满足，因为基本方程就表示了节点能量平衡。•由于热传导的傅立叶定律用于推导基本方程并用于从单元温度梯度中求解单元热通量，因而自然得到满足。•一般来说，稳态分析中网格上节点温度比实际温度要低。也就是说，如果加密网格，温度将增加，但加密到一定程度，结果将不显著增加(也就是说,结果达到收敛)。T网格密度热传导基础有限元热分析的基本特点(续)•引起温度奇异性的原因–整体求解的奇异性•在稳态分析中当有热量输入(比如,施加节点热流，热流，内部热源)而无热流流出(指定的节点温度，对流载荷等)，稳态的温度将是无限大的。•等同于结构分析中的刚体位移。–温度梯度/热流奇异性•如果对点热源处的网格细分下去的话，梯度/热流将无限增加。•凹角和网格中的“裂缝”。•形状不好的单元。热传导基础有限元热分析的基本特点(续)•实际上任何产生不连续热通量区域的有限元模型都是有误差的。在单元边界上的热通量不连续的大小将作为ANSYS误差估计的基础。•网格划分误差估计一般用于实体和壳单元，而且单元所在区域的单元类型是相同的(具有共同的特性)，热通量在该区域中也就是连续的。•在ANSYS理论手册中对误差的计算有详细的叙述。热传导基础误差估计热传导基础ANSYS误差估计•ANSYS计算了几个数值，可以用来评估网格划分误差。误差计算可以用于线性和非线性的稳态分析，在通用后处理器-POST1中进行(FullGraphics设置为ON)。•ANSYS中的网格划分误差度量功能:–TEPC–能量范数百分误差，表示由于特定的网格划分而引起的相对误差。要想知道应该在什么地方细化网格，可绘制TERR(详见下面描述)–TERR-估计选定单元中的热耗散能。单位是能量单位，比如,BTU,焦耳等。在POST1中可以使用ETABLE命令存储，排序和列表。TERR的云图可以使用：ContourPlotElementSolution来显示。–TDSG－单元中最大的热流偏差。计算单元中每个节点在任意方向上平均热流和非平均热流之间的最大差值。单位是热流单位，比如，BTU/(hour-in2；J/(s-m2))。存储，排序，列表和绘图方法与TERR类似。热传导基础ANSYS误差估计(续)•网格划分误差度量(续)–误差限SMNB和SMXB-当用云图绘制不连续数值(温度梯度和热流)时(误差估计功能处于打开状态),SMNB和SMXB将出现在图例区域，表示出该数值不连续的范围。热传导基础ANSYS误差估计(续)•网格划分误差度量(续)–例如:假如云图显示的是X方向的平均节点热流(PLNS,TF,X),SMNB和SMXB将显示在图例中，其计算方法如下:SMXBqTDSGSMNBqTDSGqTDSGixixixmax()min()maxmaxmaxforallselectednodesforallselectednodeswhereaveragenodalfluxinthex-directionatnodeiandmaximumTDSGofanyselectedelementthatconnectstonodei.热传导基础如何使热传递分析包括非线性?•当比热矩阵，热传导率矩阵和/或等效节点热流向量至少一个是温度的函数时，分析就是非线性的，需要迭代求解平衡方程。如果所有三项都是与温度有关的，那么控制方程可以写为如下形式:•下面几项都可以使得分析包括非线性:–与温度有关的材料特性–与温度有关的对流换热系数–使用辐射单元–与温度有关的热源(热流或热流矢量)–使用耦合场单元(假设载荷向量耦合)CTTKTTQT()()()mrlqlq热传导基础何时需要定义比热和密度?•瞬态问题,这些数值用于形成比热矩阵(该矩阵表示瞬态分析中需要的热能存储效果)；•稳态分析中包括有热质量传递效果(例如,模型中有流动导体介质)。热传导基础与结构分析的比较结构•位移•力•压力•应变•应力•温度分布•惯性载荷•弹性基础•无•接触热•温度•热流率•热通量(施加的)•温度梯度•热通量(计算的)•内部热生成(热/体积)•无•对流•辐射•恒温器对于熟悉结构分析的人来说，下面的表格将是非常有帮助的:CTKTQmrlqlqMUCUKUFmrmrlqlq热传导基础单元库•1-D热单元(可以与1-D,2-D和3-D单元联合使用)。LINK31辐射连接单元LINK322-D传导杆LINK333-D传导杆•1-D热单元(续)LINK34节点－节点对流连接单元MASS71集中热质量单元可以用于定义与温度有关的热源热传导基础单元库•控制单元-允许用户在有限元模型中加入反馈。最简单的方法-恒温器!–根据控制节点K或L的温度或温度差，一阶或二阶导数，温度积分，或时间，程序可以打开或关闭节点I和J之间的热流。COMBIN37节点－节点控制单元热传导基础单元库•2-D实体-传导–平面或轴对称(几何,载