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

1、6.稳态热传导问题的有限元法本章的内容如下：6热传导方程与换热边界6.2稳态温度场分析的一般有限元列式6.3三角形单元的冇限元列式6.4温度场分析举例6热传导方程与换热边界在分析工程问题时，经常要了解工件内部的温度分布情况，例如发动机的工作温度、金 属工件在热处理过程屮的温度变化、流体温度分布等。物体内部的温度分布収决于物体内部 的热量交换，以及物体与外部介质z间的热量交换，一般认为是与时间相关的。物体内部的 热交换采用以下的热传导方程（fourier方程）來描述，dt dx i_8_+ 內1dy)d+ 一8z2鲁卜。i qz丿(6-1)式中。为密度，kg/m c为比热容，j/（kgk）；,

2、av, 为导热系数，w/（mk）； t为温度，°c；（为时间，s； 为内热源密度，w/naa2t “ a2t对于各向同性材料，不同方向上的导热系数相同，热传导方程可写为以下形式,(6-2)显:7卑+2卑+2卑+6t 5x2 sy2 az2除了热传导方程，计算物体内部的温度分布，还需要指定初始条件和边界条件。初始条 件是指物体最初的温度分布情况，tt=0 =t0(x,y,z)(6-3)边界条件是指物体外表而与周围环境的热交换情况。在传热学中一般把边界条件分为三 类。1）给定物体边界上的温度，称为第一类边界条件。物体表血上的温度或温度函数为已知，或 t =t,（x,y,z,/）（6-4）

3、s2）给定物体边界上的热量输入或输出，称为第二类边界条件。 已知物体表面上热流密度，3) 给定对流换热条件，称为第三类边界条件。物体与其相接触的流体介质z间的对流换热系数和介质的温度为已知。dt , dt(6-6)dtdz 、其中h为换热系数，w/(m2k)； 7；是物体表面的温度；7是介质温度。如果边界上的换热条件不随时间变化，物体内部的热源也不随时间变化，在经过一定时 间的热交换后，物体内各点温度也将不随时间变化，即竺=0dt这类问题称为稳态(steady state)热传导问题。稳态热传导问题并不是温度场不随时间 的变化，而是指温度分布稳定后的状态，我们不关心物体内部的温度场如何从初始状

4、态过渡 到最后的稳定温度场。随时间变化的瞬态(transient)热传导方程就退化为稳态热传导方程, 三维问题的稳态热传导方程为，a5xdx)乜丿聖、dz丿(6-7)对于各向同性的材料，可以得到以卜-的方程，称为poisson方程,卑*卑*啤+2=odx2 dy2 qz a(6-8)考虑物体不包含内热源的情况，各向同性材料中的温度场满足laplace方程,(6-9)a2t a2t a2t在分析稳态热传导问题吋，不需要考虑物体的初始温度分布对最后的稳定温度场的影 响，因此不必考虑温度场的初始条件，而只需考虑换热边界条件。计算稳态温度场实际上是 求解偏微分方程的边值问题。温度场是标量场，将物体离散

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

6、al method)建立稳 态温度场分析的有限元列式。微分方程的边值问题，可以一燉地表示为未知函数u满足微分方程组，a（w）a（w） = < a2（u） > = 0（在域。内）（6-10） 未知函数u还满足边界条件，bgb（w） = < b2（m） > = 0（在边界上）（6-11） 如果未知函数u是上述边值问题的精确解，则在域屮的任一点上u都满足微分方程 （6-10）,在边界的任一点上都满足边界条件（6-11）。对于复杂的工程问题，这样的精确解 往往很难找到，需耍设法寻找近似解。所选取的近似解是一族带有待定参数的已知函数，一 般表示为« u 帀=工 niai

7、 = nai=(6-12)其小勺为待定系数，n,为己知函数，被称为试探函数。试探函数要取占完全的函数序列,是线性独立的。由于试探函数是完全的函数序列，任一函数都对以用这个序列来表示。 采用这种形式的近似解不能精确地满足微分方程和边界条件，所产生的谋差就称为余fio微分方程（6-10）的余量为，r = a(na)(6-13)边界条件（611）的余量为，r = b(na)(6-14)选择一族已知的函数，使余量的加权积分为零，强迫近似解所产生的余量在某种平均意 义上等于零，£ w：rdq+ wjr/r = 0（6-15）w/和兀称为权函数，通过公式（6-15）可以选择待定的参数勺。这种采用

8、使余量的加权积分为零来求得微分方程近似解的方法称为加权余量法。对权函 数的不同选择就得到了不同的加权余量法，常用的方法包括配点法、了域法、最小二乘法、 力矩法和伽辽金法（galcrkin method）o在很多情况下，采用galcrkin法得到的方程组的系 数矩阵是对称的，在这里也采用galerkin法建立稳态温度场分析的一般有限元列式。在 galerkin法屮，直接采用试探函数序列作为权函数，取= n., wj = -n . 0下而用求解二阶常微分方程为例，说hj galerkin法（参见，王勖成编苦“有限元法基本 原理和数值方法”的123节）。例，求解二阶常微分方程j 2牛 + u + x

9、 = 0（0 < x < 1）dx2边界条件：当x = 0时，w=0；当兀=1时，w=0o 取两项近似解：n、= x(l - x)n2 =x2(1-x)u = n仇 + n2a2 = ax(y 一工)+ a2x (1 一 x)州=m, w2 = n2由公式(6-15)可以得到两个加权积分方程，x( - x)x + ax (-2 + x-x2) + 6z2(2-6x + x2 - x3)dx = 0(1 兀)兀 + % (_2 x (2 6x + 兀2 _ )dx = 0积分后可以得到一个二元一次方程组，解得，% =0.1924,a 2 = 0.1707近似解为,u=x(l-x)(0

10、.1924 + 0707x)该方程的粕确解为，u=-xsinl近似解与楮确解的结果比较见表6-1,表6-1近似解与精确解比较x=025x=0.5x=0.75sinxu =xsin 10.044010.069750.06006m = x(l - x)(0.1924 + 0.1707x)0.044080.069440.06008假定单元的形函数为，n = n,皿nj单元结点的温度为，te = t2 . tnt单元内部的温度分布为，r = 7vte以二维问题为例，说明用galerkin法建立稳态温度场的一般有限元格式的过程。二维问 题的稳态热传导方程为，(6-16a)聖、dx)第一类换热边界为(6-

11、16c)第二类换热边界条件为,oxdy第三类边界条件为,(6-16d)在一个单元内的加权积分公式为,；叩？仇浮)+ £仇普) + odg = o1 ox ex dy dy(6-17)由分部积分得,別讥甞)=譽“乎+時仇乎知认寥)=譽久甞)+時0乎dydy dy dy dy dy应用green定理，一个单元内的加权积分公式写为, r r5yvj of svvj cti叽孰+4轨)dr = 0采jij galerkin方法，选择权函数为，(6-18)将单元内的温度分布函数和换热边界条件代入(6-18)式,£ >(4+讐久警川门掘4 dxoy dy-£ nqd&#

12、163;l-2niqsdv+ l tv""""- nihtfdr = ody dy单元的加权积分公式为,(6-19)换热边界条件代入后，在(6式内相应岀现了第二类换热边界项-j ndt,第三类换热边界项匸n/n7t"-匚ngdr ,但没有出现与笫一类换热边界对应的项。这是因为，釆用tv,作为权函数，第一类换热边界被自动满足。写成矩阵形式有,(叨罗4 oxox oycy(6-20)-£ n! qdsl-n1 cisdr+灿77- nrhtfdv = q公式(6-20)是n个联立的线性方程组，可以确定n个结点的温度7；。按有限元格式将(6

13、-20)表示为,kyte=pe(6-21)其屮炬阵kf为单元的导热矩阵或称为温度刚度矩阵，t)e为单元的结点温度向量，p° 称为单元的温度载荷向量或热载荷向量(thermal load vector)o对于某个特定单元，单元导 热矩阵kf和温度载荷向量pf的元素分别为，(6-22)的 cnjc：n(，nj 十hnjnjdtx dx dx y dy dyz nihtfdt+ £ n©dv(6-23)如果某个单元完全处于物体的内部,k厂(人如殂+人殂殂jx dx dx )dy dyp, = £ n © 在在整个物体上的加权积分方程是单元积分方程的和

14、,ze讐)l7tdq(6-24)-e £【ng工帥15声+ee工 匚灿ntntecir-y £ nthtfcir = 0根据单元结点的局部编号与整休编号的关系，血接求和得到整体刚度矩阵，整体方程纽 为，kt = p6.3三角形单元的有限元列式冋顾第三章的内容可以发现，与计算弹性力学平而问题时所采用的方法一样，二维温度 场问题计算中所采用的三角形单元可以使用相同的形函数，ni =方（也+勺兀+沙nj =云（勺+巧兀+勺刃v” 二丄加 + x + 5）2aai = xjymxmyj bj =yi-ym cj =xtn-xjaj =xmyi-xiymbj = ym - cj =无

15、-兀”am =兀儿一勺几 bm =yi-yj cm =xj-xi1 x： y, 11 xj yj =t j心儿_t =2a单元内的温度分布用结点上的温度值表示为,、t=lnt nj nm t（6-25）t川丿在三角形单元上，采用galerkin法可得，(6-26)可吟 h (血 m) + da = 0 ox ox oy oy灯？(人冬ox ox oxoxdx dx假定单元内的导热系数为常数,4a2(6-27)4a®bjbib,nb®bjb几sbt,町软岛)(心浮dy dy dydydy dyj血皿竺da 山 dy dygm >dajtrri jc?(6-28)a4a单

16、元的刚度矩阵为，b-k b”jb；bbmc.2wjac ；55cjj£显然，单元的导热矩阵是对称的。如果单元的内部热源为常数，由内部热源产生的温度载荷项为,ntqda = q由green公式可得f °“入口7 0t5t dt（"人=）+（"人必ja dxdx dydy(6-29)(6-30)=j （仏轨+可九轨）ds 方便起见，把换热边界统一表示为第三类换热边界，夕(ns浮)+ g(n玄知必oxox dyoy(6-30)灿 nv (t-ts)ds = 灿 nytjds-?灿 nrw7tds如果在单元边上存在热交换，各条边上的边界换热条件在单元刚度矩阵屮生

17、成的附加项为,ke21_0(6-31)2hlk“ =亠 0012(6-32)(6-33)山边界换热条件牛成的温度载荷向量为,py(6-34)pe-/(6-35)py(6-36)6.4温度场分析举例止方形截而的烟囱如图6-2所示，烟囱由混凝土建造，边长为60cm,通道的边长为20cm, 混凝土的导热系数为r=1.4w/(mk)。假定烟囱内表面的温度为100°c,烟囱外表面暴鉤 在空气中，空气的温度为30°c,换热系数为h = 20w /(/h2-)o计算烟囱截面内的稳态温 度场。(参见，finite element method theory and application

18、with ansys, p279 )20cm60cm图62烟囱截面modal solutionstep"丄sub .丄time 丄temp(avg)b3y3-03hm *34.1623hx 丄0034.182图6-3有限元模型46.808ore is zooz23.33.s541.49556.12170.74885.374100图6-4稳态温度分布ansys图6-5热流量分布稳态温度场分布与物体的初始状态无关,那么是否与材料的导热系数相关？我们把烟囱 的模型做些修改，假定烟囱壁由两层材料构成。内层材料为混凝土，外表面的截面尺寸为 30cmx30cm ,烟囱通道的尺寸不变,仍为20cmx20cm o外层材料的导热系数为 k = 0aw /（m-a：）,外部表面的截血尺寸不变，内部表面的截面尺寸为30cm x 30cm o换热边界条件不变，双层烟囱的冇限元模型如图6-6所示。图6-6双层烟囱的有限元模型1hodal solution3tep=1sub =1图67双层烟囱的温度分布dec 19 200200:11:0737.97853.4846