1、,有限元法基础,1,第11章 传热分析与热应力,11 传热分析与热应力,本章参考书 孔祥谦. 有限单元法在传热学中的应用. 第三版，科学出版社，1998 王勖成. 有限单元法. 清华大学出版社，2003. 第12章 曾攀. 有限元分析及应用. 清华大学出版社，2004 . 第8章,2,11. 传热分析与热应力,传热是广泛存在的自然现象，只要有温度差存在，就会有热量的传递，只要有热量的输入和输出，就会引起温度的变化。 传热分析的对象是固体、液体和气体，其应用包括热量交换、化学反应、材料相变、能量转换等。 温度的变化和不均匀分布，引起结构出现应力变化，称为热应力。 当以应力分析为目的时，为确定温度
2、场，需要对固体进行传热计算，以便确定相关的热应力。,3,11 传热分析与热应力,11.1 传热问题的基本方程 固体热传导的现象,4,11 传热分析与热应力,术语和单位在国际标准单位制中，传热分析的术语和单位,5,11 传热分析与热应力,控制方程对于微元dxdydz，生成的热量为微元体内的净流出热流量为由于热量的储存使内能增加，即由能量守恒定律，在微元内有,6,11 传热分析与热应力,由于方程不封闭，需补充条件。各向同性均匀材料的Fourier定律于是，有 第一类边界条件：边界上已知温度T，即第二类边界条件：边界上已知热流密度，即第三类边界条件：已知与物体相接触的流体介质的温度和换热系数为已知，
3、即,7,11 传热分析与热应力,11.2 变分原理与有限元瞬态热传导问题变分泛函为稳态问题，温度不随时间变化泛函的变分取驻值，可得控制方程和第二类和第三类边界条件第一类边界条件应强制满足，称为本质边界条件；第二、第三类边界条件是自然边界条件。,8,11 传热分析与热应力,有限元法将物体离散为n个单元体，即 ，将单元内的温度场用节点上的温度插值，有代入泛函，泛函称为,9,11 传热分析与热应力,将单元矩阵组合后，得到物体的总体矩阵，对泛函变分取极值，即得有限元方程组KT称为热传导矩阵，C为热容矩阵，RT为等效节点温度载荷列阵。对稳态问题，上述方程是线性的，热传导方程是对称正定的，求解类方法似于结
4、构分析。 对瞬态问题，需采用有限差分法，将 离散，或采用显式时间积分，如中心差分，或采用隐式时间积分，如Newmark差分。,10,11 传热分析与热应力,对线性瞬态问题，还可采用结构分析中的模态迭加法首先求特征问题每个特征向量T i相对于C正则化，即 令 是模态矩阵，它的每一列是正则化的特征向量T i，于是可将节点温度表示为广义温度Z的关系将其代入有限元方程，并左乘 得到n个解耦的方程组积分上述方程组后，得Z(t)，由此可得到节点T(t)。,11,11 传热分析与热应力,11.3热辐射考虑两个无限大的平行平面，由于无限大，不用考虑边界效应。设每个平面都有均匀温度，平面1的温度为T1，平面2的
5、温度为T2，平面都是理想的黑体，因此每个平面都是理想的吸收体和辐射体，平面表面的热流量为是Stefan-Boltzman常数。由于实际的辐射面并非理想黑体，也不是无穷大的平面，也不一定平行，因此把面积为A1和温度T1表面所接收的热流量表示为包含了各种因素引起的辐射折减，包括视图因子和辐射率等。,12,11 传热分析与热应力,由于辐射面是有限的、非平行的，用视图因子表示对于两个无限大的平行面为1，对于两个相互看不见的平面是0,13,11 传热分析与热应力,与面积为A1交换辐射能的表面有多少个，就有多少个式子。如果A1不是很大，可认为Q1在A1上是个常数，因此与对流边界比较，该式与之相同。因此在有
6、限元方程中将h项中换成 就可用有限元分析辐射问题。 需注意， 是与温度有关的，故辐射问题是高度非线性问题。 当材料常数是温度的函数，问题也是非线性的。 非线性方程可采用Newton-Raphson法求解，但有一些特殊的适合传热问题的处理方法。 为了避免辐射的强非线性，实际问题的处理，有时将 也处理为常数，如在热锻时的传热分析。,14,11 传热分析与热应力,11.4伴有相变的导热问题 特点：控制方程是非稳态导热方程，区域内存在一个随时间移动的两相界面，在界面上放出或吸收潜热。1891年J. Stefan关于地极冰层厚度的研究首次讨论这一课题。 当越过相变区间时，热流密度不连续，在数学上是一个强
7、非线性问题，计算发生困难。,15,11 传热分析与热应力,相变界面的两边各自满足非稳态导热控制方程，一般为了简单略去液相区的自然对流或强制对流等作用。 在相变界面S(t)上，满足温度连续条件能量守恒条件设L J/kg为物质的相变潜热，则由于相变界面的移动，给数值方法带来困难。,16,11 传热分析与热应力,焓法模型采用焓（HcT）和温度同时作为待求函数。由于相变界面上温度随时间的变化曲线是间断的，但焓随时间的变化曲线是连续的，因此用数值方法求解焓分布时不需跟踪两相界面，从而使液相区和固相区统一处理称为可能，焓场解出后，温度场可容易得到。,17,11 传热分析与热应力,焓表示的传热方程为其中温度
8、与焓的关系为,18,11 传热分析与热应力,有限元法中的处理 （1）相变潜热作为附加比热用加权余量法推导有限元方程，有使用T 的插值的形函数作为权函数，可得Galerkin法的有限元方程。 式中fs是固相率,为无因次量, , 液相0, 固相为1。,19,11 传热分析与热应力,（2）相变潜热作为源项在前面的公式中，相变潜热项 与内热源项 作相同处理，得有限元方程是与相变项有关的等效载荷。由于相变得复杂性，在实际数值模拟中，尤其在大型结构计算时，不考虑相变的两相问题，相变的影响只考虑潜热，把它作为一个常热源处理。,20,11 传热分析与热应力,11.5 热应力（Thermal stress）大多
9、数情况下，传热问题所确定的温度场将直接影响物体的热应力,而热应力对温度场的耦合影响不大，因而可将物体的热问题看成是单向耦合过程。 （一）热应力问题中的物理方程热膨胀系数（Thermal expansion coefficent）, 为Lam系数,21,11 传热分析与热应力,将热应变看作初始应变，上式写为其中初应变为 ，矩阵形式 （二）虚功原理弹性力学问题的虚功原理为内力虚功等于外力虚功，即将物理方程代入，得,22,11 传热分析与热应力,（三）有限元列式设单元节点列阵为 假设单元内位移由节点位移表示的插值函数为 应变可表示为 虚位移与虚应变为 代入虚功原理，得 其中组装到总体矩阵后，由于 的
10、任意性，得到有限元方程,23,11 传热分析与热应力,（三）求解热应力的方法在有限元分析程序中解热应力问题有两种方法，即直接法和间接法。 直接法直接将传热分析和热应力耦合起来分析的方法。在求解时，直接将传热边界条件、力学边界条件施加在有限元模型上，以节点温度和位移作为未知变量求解。在有限元商业软件中，有多场耦合单元，如ANSYS中Solid5和Solid98有热、电、磁、压电和结构的耦合场单元，每节点有6DOF，即温度、电势、磁场强度势函数，3个方向位移。,24,11 传热分析与热应力,间接法热应力问题是一个单向耦合问题，在多数情况下，温度变化会产生热应变，但热应变部不引起温度变化。这样，可以
11、将热应力问题分成两个过程来计算分析，即传热分析和热应力计算。首先通过传热分析获得结构的温度场分布，然后在已知温度分布的情况下求热应力。优点：1）求解规模比直接法小，有更高的计算效率；2）在瞬态问题时，可先求解温度场，然后在关心的时间点上求热应力，可节省大量的存储空间，即使是非线性问题也可这样处理。,25,11 传热分析与热应力,例 稳态温度分析 一圆柱型容器上，垂直接有一小管道， 容器内装有450F的液体，小管内有100F 的流体。假设容器足够长，远端的温度 为450F。,26,11 传热分析与热应力,例 稳态温度长分析（续）,网格图,温度分布,27,11 传热分析与热应力,例 稳态温度常分析
12、（续）,28,11 传热分析与热应力,例 稳态温度长分析（续）等效应力分布图,29,11 传热分析与热应力,例 液固相变分析有厚度为a得液体，初始温度为T0, 在其上表面突加温度TsT0，分析液体 的温度变化及相变情况，假设液体的 其他边界是绝热的。,30,11 传热分析与热应力,例 液固相变分析（续）相变潜热（latent heat）的影响以焓的快速变化来反映。焓（enthalpy）与温度的变化曲线,31,例 液固相变分析（续）经计算分析，在789s797s间液体开始完全固化，即节点2开始降温至1 。 在501s时温度沿厚度的分布 各节点温度随时间的变化,11 传热分析与热应力,32,11
13、传热分析与热应力,例 热处理分析将工件放在一定介质中加热到适宜的温度、并在此温度中保持一定时间后、又以不同速度冷却的一种工艺方法。热处理不改变工件的形状和整体化学成分，而是通过改变工件内部的显微组织，或改变工件的表面的化学成分，赋予或改变工件的使用性能。特点：1）有热传导2）有相变3）热力耦合,33,11 传热分析与热应力,例 热处理分析(续)淬火分析的网格和 温度分布,34,2 传热分析与热应力,35,11 传热分析与热应力,热处理分析（续）典型的导热系数曲线,36,11 传热分析与热应力,热处理分析（续）典型的质量密度曲线,37,11 传热分析与热应力,热处理分析（续）典型的比热容曲线,3
14、8,11 传热分析与热应力,热处理分析（续）典型的相变焓曲线,39,11传热分析与热应力,热处理分析（续）弹性模量曲线,40,11 传热分析与热应力,热处理分析（续）屈服应力曲线,41,11 传热分析与热应力,热处理分析（续）淬火实验设备及淬火介质的换热系数曲线,42,11 传热分析与热应力,热处理分析（续）几种钢的金相对比,43,2 传热分析与热应力,热处理分析（续）,44,11 传热分析与热应力,热处理分析（续）,45,11 传热分析与热应力,热处理分析（续）,46,11 传热分析与热应力,焊接残余应力通过改变加辅助热源改变残余应力分布,47,11 传热分析与热应力,焊接残余应力通过改变加辅助热源改变残余应力分布,48,11 传热分析与热应力,T形焊接,49,11 传热分析与热应力,横梁的焊接残余应力焊接现场,50,11 传热分析与热应力,横梁的焊接残余应力有限元模型,51,11 传热分析与热应力,横梁的焊接残余应力温度场分布,52,11 传热分析与热应力,横梁的焊接残余应力沿焊缝方向的残余应力分布,53,11 传热分析与热应力,横梁的焊接残余应力垂直焊缝方向的残余应力分布,54,
