稳态热传导的有限元分析

在分析工程问题时,经常要了解工件内部的温度分布情况,例如发动机的工作温度、金属工件在热处理过程中的温度变化、流体温度分布等。物体内部的温度分布取决于物体内部的热量交换,以及物体与外部介质之间的热量交换,一般认为是与时间相关的。物体内部的热交换采用以下的热传导方程(Fourier方程)来描述,

        (6-1)

式中为密度,kg/m3为比热容,为导热系数,;T为温度,℃;t为时间,s;为内热源密度,w/m3

对于各向同性材料,不同方向上的导热系数相同,热传导方程可写为以下形式,

            (6-2)

除了热传导方程,计算物体内部的温度分布,还需要指定初始条件和边界条件。初始条件是指物体最初的温度分布情况,

        (6-3)

边界条件是指物体外表面与周围环境的热交换情况。在传热学中一般把边界条件分为三类。

  1. 给定物体边界上的温度,称为第一类边界条件。

物体表面上的温度或温度函数为已知,

或                (6-4)

  1. 给定物体边界上的热量输入或输出,称为第二类边界条件。

    已知物体表面上热流密度,

或                    (6-5)

  1. 给定对流换热条件,称为第三类边界条件。

    物体与其相接触的流体介质之间的对流换热系数和介质的温度为已知。

        (6-6)

其中h为换热系数,W/(m2 K);是物体表面的温度;是介质温度。

如果边界上的换热条件不随时间变化,物体内部的热源也不随时间变化,在经过一定时间的热交换后,物体内各点温度也将不随时间变化,即

这类问题称为稳态(Steady state)热传导问题。稳态热传导问题并不是温度场不随时间的变化,而是指温度分布稳定后的状态,我们不关心物体内部的温度场如何从初始状态过渡到最后的稳定温度场。随时间变化的瞬态(Transient)热传导方程就退化为稳态热传导方程,三维问题的稳态热传导方程为,

        (6-7)

对于各向同性的材料,可以得到以下的方程,称为Poisson方程,

        (6-8)

考虑物体不包含内热源的情况,各向同性材料中的温度场满足Laplace方程,

            (6-9)

在分析稳态热传导问题时,不需要考虑物体的初始温度分布对最后的稳定温度场的影响,因此不必考虑温度场的初始条件,而只需考虑换热边界条件。计算稳态温度场实际上是求解偏微分方程的边值问题。温度场是标量场,将物体离散成有限单元后,每个单元结点上只有一个温度未知数,比弹性力学问题要简单。进行温度场计算时有限单元的形函数与弹性力学问题计算时的完全一致,单元内部的温度分布用单元的形函数,由单元结点上的温度来确定。由于实际工程问题中的换热边界条件比较复杂,在许多场合下也很难进行测量,如何定义正确的换热边界条件是温度场计算的一个难点。

6.2稳态温度场分析的一般有限元列式

在前面我们已经介绍了有限元方法可以用来分析场问题,稳态温度场计算是一个典型的场问题。我们可以采用虚功方程建立弹性力学问题分析的有限元格式,推导出的单元刚度矩阵有明确的力学含义。在这里,介绍如何用加权余量法(Weighted Residual Method)建立稳态温度场分析的有限元列式。

微分方程的边值问题,可以一般地表示为未知函数u满足微分方程组,

        (在域内)        (6-10)

未知函数u还满足边界条件,

        (在边界上)        (6-11)

如果未知函数u是上述边值问题的精确解,则在域中的任一点上u都满足微分方程(6-10),在边界的任一点上都满足边界条件(6-11)。对于复杂的工程问题,这样的精确解往往很难找到,需要设法寻找近似解。所选取的近似解是一族带有待定参数的已知函数,一般表示为

        (6-12)

其中为待定系数,为已知函数,被称为试探函数。试探函数要取自完全的函数序列,是线性独立的。由于试探函数是完全的函数序列,任一函数都可以用这个序列来表示。

采用这种形式的近似解不能精确地满足微分方程和边界条件,所产生的误差就称为余量。

微分方程(6-10)的余量为,

            (6-13)

边界条件(6-11)的余量为,

            (6-14)

选择一族已知的函数,使余量的加权积分为零,强迫近似解所产生的余量在某种平均意义上等于零,

            (6-15)

称为权函数,通过公式(6-15)可以选择待定的参数

这种采用使余量的加权积分为零来求得微分方程近似解的方法称为加权余量法。对权函数的不同选择就得到了不同的加权余量法,常用的方法包括配点法、子域法、最小二乘法、力矩法和伽辽金法(Galerkin method)。在很多情况下,采用Galerkin法得到的方程组的系数矩阵是对称的,在这里也采用Galerkin法建立稳态温度场分析的一般有限元列式。在Galerkin法中,直接采用试探函数序列作为权函数,取

下面用求解二阶常微分方程为例,说明Galerkin法(参见,王勖成编著"有限元法基本原理和数值方法"的1.2.3节)。

例,求解二阶常微分方程

边界条件:当时,;当时,

取两项近似解:

,    

由公式(6-15)可以得到两个加权积分方程,

积分后可以得到一个二元一次方程组,解得,

近似解为,

该方程的精确解为,

近似解与精确解的结果比较见表6-1,

表6-1 近似解与精确解比较

 
x=0.25


x=0.5


x=0.75



0.04401


0.06975


0.06006



0.04408


0.06944


0.06008

假定单元的形函数为,

单元结点的温度为,

单元内部的温度分布为,

以二维问题为例,说明用Galerkin法建立稳态温度场的一般有限元格式的过程。二维问题的稳态热传导方程为,

            (6-16a)

第一类换热边界为

                                (6-16b)

第二类换热边界条件为,

                (6-16c)

第三类边界条件为,

        (6-16d)

在一个单元内的加权积分公式为,

            (6-17)

由分部积分得,

应用Green定理,一个单元内的加权积分公式写为,

            (6-18)

采用Galerkin方法,选择权函数为,

将单元内的温度分布函数和换热边界条件代入(6-18)式,单元的加权积分公式为,

        (6-19)

换热边界条件代入后,在(6-19)式内相应出现了第二类换热边界项,第三类换热边界项,但没有出现与第一类换热边界对应的项。这是因为,采用作为权函数,第一类换热边界被自动满足。写成矩阵形式有,

            (6-20)

公式(6-20)是n个联立的线性方程组,可以确定n个结点的温度。按有限元格式将(6-20)表示为,

            (6-21)

其中矩阵[K]e为单元的导热矩阵或称为温度刚度矩阵,{T}e为单元的结点温度向量,{P}e称为单元的温度载荷向量或热载荷向量(Thermal load vector)。对于某个特定单元,单元导热矩阵[K]e和温度载荷向量{P}e的元素分别为,

        (6-22)

            (6-23)

如果某个单元完全处于物体的内部,

在整个物体上的加权积分方程是单元积分方程的和,

        (6-24)

根据单元结点的局部编号与整体编号的关系,直接求和得到整体刚度矩阵,整体方程组为,

6.3三角形单元的有限元列式

图6-1 三角形单元

回顾第三章的内容可以发现,与计算弹性力学平面问题时所采用的方法一样,二维温度场问题计算中所采用的三角形单元可以使用相同的形函数,

单元内的温度分布用结点上的温度值表示为,

            (6-25)

在三角形单元上,采用Galerkin法可得,

        (6-26)

假定单元内的导热系数为常数,

    (6-27)

    (6-28)

单元的刚度矩阵为,

显然,单元的导热矩阵是对称的。

如果单元的内部热源为常数,由内部热源产生的温度载荷项为,

        (6-29)

由Green公式可得

        (6-30)

方便起见,把换热边界统一表示为第三类换热边界,

        (6-30)

如果在单元边上存在热交换,各条边上的边界换热条件在单元刚度矩阵中生成的附加项为,

        (6-31)

        (6-32)

        (6-33)

由边界换热条件生成的温度载荷向量为,

            (6-34)

            (6-35)

            (6-36)

6.4温度场分析举例

正方形截面的烟囱如图6-2所示,烟囱由混凝土建造,边长为60cm,通道的边长为20cm,混凝土的导热系数为。假定烟囱内表面的温度为100℃,烟囱外表面暴露在空气中,空气的温度为30℃,换热系数为。计算烟囱截面内的稳态温度场。(参见,Finite Element Method Theory and Application with ANSYS, p279)

    

图6-2 烟囱截面                        图6-3 有限元模型

图6-4 稳态温度分布

图6-5 热流量分布

稳态温度场分布与物体的初始状态无关,那么是否与材料的导热系数相关?我们把烟囱的模型做些修改,假定烟囱壁由两层材料构成。内层材料为混凝土,外表面的截面尺寸为,烟囱通道的尺寸不变,仍为。外层材料的导热系数为,外部表面的截面尺寸不变,内部表面的截面尺寸为。换热边界条件不变,双层烟囱的有限元模型如图6-6所示。

图6-6 双层烟囱的有限元模型

图6-7 双层烟囱的温度分布

图6-8 双层烟囱的热流量分布

对比两种不同结构烟囱的温度分布和热流量分布,稳态温度场分布与材料的换热系数相关。双层烟囱外层材料的导热系数比较小,接近保温材料,热量很快传进内层烟囱,但向外部环境传热慢了很多,所以内层烟囱的温度很高。比较热流量分布,保温材料的却能够有效的阻止热量散失。北方城市的供暖管道都包有一层隔热材料就是基于这个道理。

时间: 2024-11-05 02:24:51

稳态热传导的有限元分析的相关文章

SR latch D latch D filp-flop SR触发器 D触发器 D双稳态多谐震荡器 【数字电路】

SR  latch D latch D filp-flop SR Latch 典型的SR 触发器就长上面那样啦 不要以为简单,酒吧 S R Q Q' 之间的相对位置关系记错了! 组成的元素很简单--或非门,但是实现的功能却很重要 SR是两个输入端口,Q 和 Q' 是输出端口 不要慌,两个端口,就四种输入情况,一个个分析,就可以了解SR latch是怎么工作的 当 s == 1 , R == 0, 或门,只要是有1 ,或门处理结果就是1. 这里S端口输入的是1,那么Q' 得到的是0,这个0反馈到R

Neotec PIPEFLO v8.5.12 1CD(集输管网稳态软件)

Optcalc v2001 1CDLED.Tool.v5.0 1CDTechWiz LCD 3D v15.0.10.1202-ISO 1CD(液晶模拟系统)TFCalc.v3.5.6(著名的光学薄膜设计软件)TFCalc 材料库TFCalc 用户手册(中文)薄膜光学——理论与实践薄膜光学的基础理论(中文)镀膜工艺与镀膜系统配置(中文)光学镀膜技术(中文)Readiris.Pro.11 1CD(中文版.領先的光學文字識別-OCR軟件)Specman.Pro.2006 1CD(对1D和2D数据进行N

[家里蹲大学数学杂志]第033期稳态可压Navier-Stokes方程弱解的存在性

1. 方程  考虑 R3 中有界区域 Ω 上如下的稳态流动: {div(?u)=0,div(?u?u)?μ△u?(λ+μ)?divu+??γ=?f+g.(1) 2. 假设  先作一些初步的假设: 2.1. γ>32 ---保证对流项 div(?u?u) 可看成 (1)2 的扰动; 2.2. μ>0 , λ+23μ≥0 ---Navier-Stokes 假设; 2.3. 考虑非滑动边界---u=0 在 ?Ω 上. 这时可将 Navier-Stokes 假设放宽为 μ>0,λ+43μ>

homeostasis稳态

目录: 1.反馈系统 2.正反馈调节 3.负反馈调节 4.behavioral response 行为反应,physiological response生理反应 1.反馈控制系统(feedback control system) 是一种“闭环”系统,即控制部分发出信号,指示受控部分活动,而受控部分的活动可被一 定的感受装置感受,感受装置再将受控部分的活动情况作为反馈信号送回到控制部分,控制部分可以根据反馈信号来改变自己的活动,调整对受控部分的指令,因而能对受控部分的活动进行调节.可见,在这样的控

2017-11-4—稳态和暂态/瞬态(对运放积分电路的思考)[待仿真]

先直接截图了,暂态或者说瞬态都是暂时的状态,是从一个稳定态到另一个稳定态的过程. 之所以要了解这个概念是因为对于使用运放搭建的模拟PID有很多的疑惑,比如负反馈没有电阻满不满足"虚短"原则? 网上的答案大多都是按照已经满足"虚短""虚断"原则下进行分析的,包括积分电路也是如此. 积分电路 如果Vi保持正值,Vo输出为0,这个电路就是一个比较器,电容两端电压固定.这么分析的话跟积分半毛钱关系都没有啊... 积分的过程只能是发生在瞬态/暂态(充放电的

解析双稳态肖特基二极管的设计

肖特基二极管是将P型半导体或者N型半导体与金属表面衬体交互作用,并最终形成稳定结构,可分为物理吸附类型与化学吸附类型.肖特基二极管由于其独特的微观结构而具有整流效应,被广泛应用于传感器.分子开关.集成电路等诸多电子器件中.现如今,如何降低肖特基高度而且同时又能保持分子的完整性,成为设计高效肖特基二极管的难点.因此目前急需设计开发出将物理吸附和化学吸附结合在一起,使其拥有双重优点的肖特基二极管.ADT分子吸附在铜表面的双稳态在实验上得到了印证,而且被应用于分子开关中,并取得了很好的效应.与传统有机

从随机过程的熵率和马尔科夫稳态过程引出的一些思考 - 人生逃不过一场马尔科夫稳态

1. 引言 0x1:人生就是一个马尔科夫稳态 每一秒我们都在做各种各样的选择,要吃青菜还是红烧肉.看电影还是看书.跑步还是睡觉,咋一看起来,每一个选择都是随机的,而人生又是由无数个这样的随机选择组成的结果.从这个前提往下推导,似乎可以得出一个结论,即人生是无常的,未来是不可预测的.但事实真的是如此吗? 以前的老人流行说一句话,三岁看小,七岁看老.这似乎是一句充满迷信主义色彩的俗语,但其实其中暗含了非常质朴而经典的理论依据,即随机过程不管其转移概率分布如何,随着时序的增大,最终会收敛在某个稳态上.

稳态视觉诱发电位 SSVEP 简介

目录 稳态视觉诱发电位(SSVEP) SSVEP信号原理与基础 获取SSVEP信号要注意的点 稳态视觉诱发电位(SSVEP) 诱发电位主要有视觉诱发电位.听觉诱发电位以及触觉诱发电位三种.视觉诱发电位由于简单.方便被广泛应用于脑电信号的研究中.当视觉收到光或者图形闪烁等刺激时,脑电信号的电位会发生变化,这些电位的变化就是视觉诱发电位(Visual Evoked Potential, VEP). 视觉诱发电位可以大致分为三类: 1.瞬时视觉诱发电位(Transient Visual EvokedP

正弦稳态电路的分析

RLC串联电路如图(a)所示,其中R=15Ω,L=12mH,C=5μF,端电压uS=100(2^(1/2))cos(5000t)V. 试求: (1)电路中的电流i(瞬时表达式)和各元件的电压相量: (2)电路的等效导纳和并联等效电路. 解(1)计算各部分阻抗为: ZR=15Ω ZL=jωL=j60Ω ZC=-j(1/ωC)=-j40Ω Zeq=ZR+ZL+ZC=(15+j20)Ω =25∠53.13°Ω 电流相量为 I·=U·S/Zeq=(100∠0°)/(25∠53.13°)A=4∠-53.1