稳态热传导问题的有限元法

1、6.稳态热传导问题的有限元法本章的内容如下：6.1热传导方程与换热边界6.2稳态温度场分析的一般有限元列式6.3三角形单元的有限元列式6.4温度场分析举例6.1热传导方程与换热边界在分析工程问题时，经常要了解工件内部的温度分布情况，例如发动机的工作温度、金 属工件在热处理过程中的温度变化、流体温度分布等。物体内部的温度分布取决于物体内部 的热量交换，以及物体与外部介质之间的热量交换，一般认为是与时间相关的。物体内部的 热交换采用以下的热传导方程(Fourier方程)来描述，+旋一6+ 6z(6-1)式中P为密度，kg/m3； c为比热容，J/(kg - K) ； X ,X ,X为导热系数，w.

2、(m- k)； Tx y z为温度，C； t为时间，s； Q为内热源密度，w/m3。对于各向同性材料，不同方向上的导热系数相同，热传导方程可写为以下形式，(6-2)6T 。62T . 62T . 62TPC有云可云除了热传导方程，计算物体内部的温度分布，还需要指定初始条件和边界条件。初始条 件是指物体最初的温度分布情况，T t_0 = T0 (x,y,z)(6-3)边界条件是指物体外表面与周围环境的热交换情况。在传热学中一般把边界条件分为三 类。给定物体边界上的温度，称为第一类边界条件。物体表面上的温度或温度函数为已知，T = T TOC o 1-5 h z 或 T = T (x, y, z,

3、 t)(6-4)给定物体边界上的热量输入或输出，称为第二类边界条件。已知物体表面上热流密度，6T. 6T. 6T . HYPERLINK l bookmark20 o Current Document n +人 n +人 n )= qx 6x x y 6y y z 6z z| = q, (x, y, z, t)(6-5)6T6T6T .(人n + X n + Xn ) x 6x x y 6y y z 6z z给定对流换热条件，称为第三类边界条件。物体与其相接触的流体介质之间的对流换热系数和介质的温度为已知。(6-6)人空n +人咀n +人空n = h(T - T )x dx x y dy y

4、z Qz z / s其中h为换热系数，W/(m2K)；气是物体表面的温度；Tf是介质温度。如果边界上的换热条件不随时间变化，物体内部的热源也不随时间变化，在经过一定时 间的热交换后，物体内各点温度也将不随时间变化，即dt这类问题称为稳态(Steady state)热传导问题。稳态热传导问题并不是温度场不随时间 的变化，而是指温度分布稳定后的状态，我们不关心物体内部的温度场如何从初始状态过渡 到最后的稳定温度场。随时间变化的瞬态(Transient)热传导方程就退化为稳态热传导方程， 三维问题的稳态热传导方程为，Q L QT)Q LcT) Q LQT)不八f 人 | +M 人M +人 l + Q

5、= 0(6-7)Qx k x Qx )Qy k yQy) Qz k zQz )对于各向同性的材料，可以得到以下的方程，称为Poisson方程，(6-8)Q 2T Q 2T Q 2T Q 八 + = 0 Qx2 Qy2 Qz2 人考虑物体不包含内热源的情况，各向同性材料中的温度场满足Laplace方程，(6-9)Q 2T + Q 2T + Q 2T = 0Qx2Qy2Qz2在分析稳态热传导问题时，不需要考虑物体的初始温度分布对最后的稳定温度场的影 响，因此不必考虑温度场的初始条件，而只需考虑换热边界条件。计算稳态温度场实际上是 求解偏微分方程的边值问题。温度场是标量场，将物体离散成有限单元后，每

6、个单元结点上 只有一个温度未知数，比弹性力学问题要简单。进行温度场计算时有限单元的形函数与弹性 力学问题计算时的完全一致，单元内部的温度分布用单元的形函数，由单元结点上的温度来 确定。由于实际工程问题中的换热边界条件比较复杂，在许多场合下也很难进行测量，如何 定义正确的换热边界条件是温度场计算的一个难点。6.2稳态温度场分析的一般有限元列式在前面我们已经介绍了有限元方法可以用来分析场问题，稳态温度场计算是一个典型的 场问题。我们可以采用虚功方程建立弹性力学问题分析的有限元格式，推导出的单元刚度矩 阵有明确的力学含义。在这里，介绍如何用加权余量法 Weighted Residual Method

7、)建立稳 态温度场分析的有限元列式。微分方程的边值问题，可以一般地表示为未知函数u满足微分方程组，A1（u ）、A（u）= = 0（在域 Q 内）（6-10）.未知函数u还满足边界条件，B（u ）、B（u） = = 0（在边界r 上）（6-11）.如果未知函数u是上述边值问题的精确解，则在域中的任一点上u都满足微分方程 （6-10）,在边界的任一点上都满足边界条件（6-11）。对于复杂的工程问题，这样的精确解 往往很难找到，需要设法寻找近似解。所选取的近似解是一族带有待定参数的已知函数，一 般表示为u 机 u =咒 Na = Na （6-12）i=1其中a,为待定系数，N,为已知函数，被称为试

8、探函数。试探函数要取自完全的函数序列， 是线性独立的。由于试探函数是完全的函数序列，任一函数都可以用这个序列来表示。采用这种形式的近似解不能精确地满足微分方程和边界条件，所产生的误差就称为余 量。微分方程（6-10）的余量为，R = A （Na）（6-13）边界条件（6-11）的余量为，R = B（Na）（6-14）选择一族已知的函数，使余量的加权积分为零，强迫近似解所产生的余量在某种平均意 义上等于零，j WTRdQ + j WTRr = 0（6-15）W和W称为权函数，通过公式（6-15）可以选择待定的参数ai。这种采用使余量的加权积分为零来求得微分方程近似解的方法称为加权余量法。对权函

9、数的不同选择就得到了不同的加权余量法，常用的方法包括配点法、子域法、最小二乘法、 力矩法和伽辽金法（Galerkin method）。在很多情况下，采用Galerkin法得到的方程组的系 数矩阵是对称的，在这里也采用Galerkin法建立稳态温度场分析的一般有限元列式。在Galerkin法中，直接采用试探函数序列作为权函数，取W. = N ., W =-N,。下面用求解二阶常微分方程为例，说明Galerkin法（参见，王勖成编著“有限元法基本 原理和数值方法”的1.2.3节）。例，求解二阶常微分方程+ u + X = 0（0 X ccc JT Tydy dyx A4A2cjmijmTjmdy第

10、）、 c2ccc c1 T 1iijim1ccc2cc1Tijjjmjc cccc 2T .imjmmmMm 4A单元的刚度矩阵为，b 2bbbbMc 2ccc cMiijijiijimxbbb2bb+ yccc2cc4i jjjm4ijjjmbbbbb 2c cccc 2_ i mmmimjmmK e =显然，单元的导热矩阵是对称的。如果单元的内部热源为常数，由内部热源产生的温度载荷项为，1 N 1111QA1 NdA = 1,j3 IN hmj NyQdA = Qj（6-29）（6-27）（6-28）由Green公式可得j : （NT M g） + （NT M gdA（6-30）A dxx

11、 dxdyy dydT. dT=j （NtM n + NtM n ）dS方便起见，把换热边界统一表示为第三类换热边界，a . dTa . dT(6-30)J 瓦（N以x瓦）+否（NT气哥dA=J hNt（T -T ）dS = J hNtT dS J hNtNTedS sf ssfs如果在单元边上存在热交换，各条边上的边界换热条件在单元刚度矩阵中生成的附加项为，2 10hlKe r61 2 00 0 0(6-31)hlK e =i60000210122 0 1K e =00061 0 2（6-32）(6-33)由边界换热条件生成的温度载荷向量为，1 hT l尸e =j 2 0（6-34）尸e =

12、 hTfljm J1 21（6-35）(6-36)6.4温度场分析举例正方形截面的烟囱如图6-2所示，烟囱由混凝土建造，边长为60cm，通道的边长为20cm, 混凝土的导热系数为k = 1.4W/（m - K）。假定烟囱内表面的温度为100。，烟囱外表面暴露在空气中，空气的温度为30。，换热系数为h = 20W/（m2 K）。计算烟囱截面内的稳态温度场。（参见，Finite Element Method Theory and Application with ANSYS, p279）图6-2烟囱截面图6-3有限元模型图6-4稳态温度分布图6-5热流量分布稳态温度场分布与物体的初始状态无关，那么是否与材料的导热系数相关？我们把烟囱 的模型做些修改，假定烟囱壁由两层材料构成。内层材料为混凝土，外表面的截面尺寸为 30cm x 30cm，烟囱通道的尺寸不变，仍为20cm x 20cm。外层材料的导热系数为k = 0.1W /(m - K)，外部表面的截面尺寸不变，内部表面的截面尺寸为30cm x 30cm。换热边界条件不变，双层烟囱的有限元模型如图6-6所