第三章非稳态导热分析解法.docx
- 文档编号:8838556
- 上传时间:2023-05-15
- 格式:DOCX
- 页数:22
- 大小:193.04KB
第三章非稳态导热分析解法.docx
《第三章非稳态导热分析解法.docx》由会员分享,可在线阅读,更多相关《第三章非稳态导热分析解法.docx(22页珍藏版)》请在冰点文库上搜索。
第三章非稳态导热分析解法
第三章非稳态导热分析解法
本章主要要求:
1、重点内容:
①非稳态导热的基本概念及特点;
2集总参数法的基本原理及应用;
3一维及二维非稳态导热问题。
2、掌握内容:
①确定瞬时温度场的方法;
②确定在一时间间隔内物体所传导热量的计算方法。
3、了解内容:
无限大物体非稳态导热的基本特点。
许多工程问题需要确定:
物体内部温度场随时间的变化,或确定其内部温度达某一极限值所需的时间。
如:
机器启动、变动工况时,急剧的温度变化会使部件因热应力而破坏。
因此,应确定其内部的瞬时温度场。
钢制工件的热处理是一个典型的非稳态导热过程,掌握工件中温度变化的速率是控制工件热处理质量的重要因素;金属在加热炉内加热时,要确定它在炉内停留的时间,以保证达到规定的中心温度。
§3-1非稳态导热的基本概念
一.非稳态导热
1、定义:
物体的温度随时间而变化的导热过程称非稳态导热。
2、分类:
根据物体内温度随时间而变化的特征不同分:
=8曲
图3・1非稳态导热
过程中的温度分布
1)物体的温度随时间的推移逐渐趋于恒定值,即:
2)物体的温度随时间而作周期性变化
如图3-1所示,设一平壁,初值温度t0,令其左侧的表面温度突然升高到勺并保持不变,而右侧仍与温度为s的空气接触,试分析物体的温度场的变化过程。
首先,物体与高温表面靠近部分的温度很快上升,而其余部分仍保持原来的t0o
如图中曲线HBD,随时间的推移,山于物体导热温度变化波及范用扩大,到某一时间后,右侧表面温度也逐渐升高,如图中曲线HCD、
HE、HFo
最后,当时间达到一定值后,温度分布保持恒定,如图中曲线HG
(若入二const,则HG是直线)o
山此可见,上述非稳态导热过程中,存在着右侧面参与换热与不参与换热的两个不同阶段。
(1)第一阶段(右侧面不参与换热)
温度分布显现出部分为非稳态导热规律控制区和部分为初始温度区的混合分布,即:
在此阶段物体温度分布受t分布的影响较大,此阶段称非正规状况阶段。
(2)第二阶段,(右侧面参与换热)
当右侧面参与换热以后,物体中的温度分布不受to影响,主要取决于边界条件及物性,此时,非稳态导热过程进入到正规状况阶段。
正规状况阶段的温度变化规律是本章讨论的重点。
2)二类非稳态导热的区别:
前者存在着有区别的两个不同阶段,而后者不存在。
3、特点;
非稳态导热过程中,在与热流量方向相垂直的不同截面上热流量不相等,这是非稳态导热区别于稳态导热的一个特点。
原因:
山于在热量传递的路径上,物体各处温度的变化要积聚或消耗能量,所以,在热流量传递的方向上羊8心。
二、非稳态导热的数学模型
1、数学模型
非稳态导热问题的求解=规定的{初始条件,边界条件}下,求解导热微分方程。
2、讨论物体处于恒温介质中的第三类边界条件问题
在第三类边界条件下,确定非稳态导热物体中的温度变化特征与边界条件参数的关系。
已知:
平板厚2X初温to、表面传热系数h、平板导热系数兄,将其突然置于温度为①的流体中冷却。
图3-3毕廉数对平板温度场变化的影响
试分析在以下三种情况:
/h、8">>l/h、"丸二1/h时,平板中温度
场的变化。
1)1/h«"入
因为1/h可忽略,当平板突然被冷却时,其表面温度就被冷却到s随着时间
的延长,平板内各点
al
tfT如图3-3(a)o
2)l/h»"入
因为"丸忽略不计,即平板内导热的流量接近于无穷大,所以任意时刻平板中各点温度接近均匀,随着时间的延长,平板内各点t一心,而且整体温度下降如图3-3(b)o
3)l/h=§5
平板中的温度分布介于二者之间,如图3-3(c)o
山此可见,表面对流换热热阻1/h与导热热阻"兑的相对大小对物体中非稳态导热的温度场的分布有重要影响,因此,引入表征二者比值的无量纲数,毕渥数。
3、毕渥数
E._%_尿
1)定义式:
XA
(3-1)
毕渥数属特征数(准则数)。
2)Bi物理意义:
Bi的大小反映了物体在非稳态条件下内部温度场的分布规律。
3)特征数(准则数):
表征某一物理现象或过程特征的无量纲数。
4)特征长度:
是指特征数定义式中的儿何尺度。
§3-2集总参数法的简化分析
一、集总参数法
1、定义:
当固体内的时,固体内的温度趋于一致,此时可认为整个固体在同一瞬间均处于同一温度下,这时需求解的温度仅是时间的一元函数,而与坐标无关,好象该固体原来连续分布的质量与热容量汇总到一点上,而只有一个温度值那样。
这种忽略物体内部导热热阻的简化分析方法称为集总参数法。
2、集总参数法的计算
已知:
有一任意形状的物体,其体积为V,面积为A,初始温度为tO,在初始时刻,突然将其置于温度恒为九的流体中,且to>S固体与流体间的表面传热系数h,固体的物性参数均保持常数。
试根据集总参数法确定物体温度随时间的依变关系。
解:
①建立非稳态导热数学模型
方法一:
据非稳态有内热源的导热微分方程:
•••物体内部导热热阻很小,忽略不计。
•••物体温度在同一瞬间各点温度基本相等,即t仅是t的一元函数,二与坐标x、y、z无关,即
dhd2tdhdx2十dy2十dz1
dt①
一
dupc
◎可视为广义热源,而且热交换的边界不是计算边界(零维无任何边
界)。
:
.界面上交换的热量应折算成整个物体的体积热源,即:
一①产=A.h Vt>I物体被冷却,•••①应为负值 a), b)式得: "狰% 由 (3-2)这就是瞬时时刻导热微分方程式。 方法二: 根据能量守恒原理,建立物体的热平衡方程,即 物体与环境的对流邱热量=物体内能的减少量乎=-Ah(t—匕)=①r 则有: ②物体温度随时间的依变关系引入过余温度: £=—9pcV—^-AhO dr &〔0)=勺一匕= 则上式表示成: 其初始条件为: 空二辽=_旦 将E~dr~“'分离变量求解微分方程,&Q卩对时间匸从0TT积分,则: 臥8尤pcV 9hA- =—£ In旳k 即: (3-3) kA qF77) kA沁? 和774) ——t=——r= .5-=BzjrPoy 其中: °7肚a戸a0/丄)2 其中: V/A是具有长度的量纲,记为打%__恥__毕渥数;呀--S--傅立叶数;而V说明Fov、Biv中的特征长度为V/A 0 “=e^pC-BtyFay) 故得: &0£0一,8(3~4) 山此可见,采用集总参数法分析时,物体内的过余温度随时间成指数曲线 1、叭 1)定义: 表征固体内部单位导热面积上的导热热阻与单位面积上的换热热阻(即外部热阻)之比。 —% 戲卩越小,表示内热阻越小,外部热阻越大。 此时采用集总参数法求解更为合适。 2)物理意义: &扩的大小反映了物体在非稳态导热条件下,物体内温度场的分布规律。 2、Fc,v 1)定义: 肥卩表征两个时间间隔相比所得的无量纲时间。 分子T是从边界上开始并生热扰动的时刻起到所计时刻为止的时间间 隔。 分母可视为边界上发生的有限大小的热扰动穿过一定厚度的固体层扩散到尸的面积上所需的时间。 2)物理意义: 表示非稳态导热过程进行的程度,兀卩越大,热扰动就越深入地传播到物体内部,因而物体内各点的温度越接近周围介质的温度。 §3—3—维非稳态导热的分析解 本节介绍第三类边界条件下: 无限大平板、无限长圆柱、球的分析解及应用。 如何理解无限大物体,如: 当一块平板的长度、宽度>>厚度时,平板的长度和宽度的边缘向四周的散热对平板内的温度分布影响很少,以至于可以把平板内各点的温度看作仅是厚度的函数时,该平板就是一块“无限大”平板。 若平板的长度、宽度、厚度相差较小,但平板四周绝热良好,则热量交换仅发生在平板两侧面,从传热的角度分析,可简化成一维导热问题。 一、无限大平板的分析解 已知: 厚度2別勺无限大平板,初温S初始瞬间将其放于温度为人的流体中,而且九〉S流体与板面间的表面传热系数为一常数。 试确定在非稳态过程中板内的温度分布。 解: 如图3-5所示,平板两面对称受热,所以其内温度分布以其中心截面为对称面。 对于X>0的半块平板,其导热微分方程: 乔丽(0 定解条件: t(x,0)= 5(0 冼(兀Q hh=0(边界条件) 妣(玄£)_心]—-火酬兀T)|x* 八"」&"(边界条件) 引入过余温度: E 7="兀巧一丄8 ded20 则 drJ去彳(0 o(x,0)=2(0 条件) 阳⑴叭ft-o 去|/=°(边界条件) 肋°T)--严GEIX-& 去(边界条件) 对偏微分方程乔萌分离变量求解得: &gT)「9狞Sin(炼6口[(忌6卡] 殆2扫色占十3啦炸占)COS(禺£) (3-10) t3R(^yy=,A2=*1,2, 其中离散值8”是下列超越方程的根,称为特征值。 色方…… (3-11) 其中Bi是以特征长度为$的毕渥数。 山此可见: 平板中的无量纲过余温度%。 与三个无量纲数有关: 以平板厚度一半,/竺◎=心」1_乞=Bi,-} 0为特征长度的傅立叶数、毕渥数及力即: s切一匕& (3-12) 2.非稳态导热的正规状况阶段 1、平板中任一点的过余温度与平板中心的过余温度的关系 前述得到的分析解是一个无穷级数,计算工作量大,但对比计算表明,当Fo>0.2时,釆用该级数的第一项与采用完整的级数计•算平板中心温度的误差小于1%, 因此,当Fo>0.2时,采用以下简化结果: 对圆柱体、 &(") 且及 歹⑺可表示为: 琴空=人吨(-圧(血Q (3-14) 此式反映了非稳态导热过程中一种很重要的物理现象: 即当Fo>0.2以后,虽然&(x,T)与叭T)各自均与T有关,但其比值则与T无关,而仅取决于儿何位置及边界条件(Bi)o也就是说,初始条件的影响已经消失, &(x,r)/ 无论初始条件分布如何,只要Fo>0.2,沧⑺之值是一个常数,也就是无量 纲的温度分布是一样的。 山此可见,当Fo>0.2时,非稳态导热过程进入正规状况阶段。 2、在一个时间间隔内非稳态导热过程中传递的热量 1)从物体初始时刻平板与周围介质处于热平衡,这一过程中传递的热量: Qq=pdQo一上8) (3-15) 此值为非稳态导热过程中传递的最大热量。 2)从初始时刻到某一时间T,这段时间内所传递的热量 00—t-z(x.巧灯? (3-16) 0O->r 3)q之比: ftlTr_3卜[力_2(兀,~0彩7_1『血一匕)_(—8)护Co-垢)_卩"巾-匕 =1-丄Lla诃=1-— (3-17) 其中: 歹=0⑴是时刻T物体的平均过余温度,一討以心o _对于无限大平板,当Fo>0.2,将式(3-13)代入©的定义式,可得: 毛®/=丄[一匕甘严—2sin“i©-3: 凤)Sind /&0ViVS_如旳+sin口\cos^Ai( 8) %=吓 球体"一臣>0.2时,无穷级数的解也可用第一项近似代替,并 =Auxp(-肿 (3-30) 其中: n为无量纲儿何位置,对平板…心、对柱体及球体"=小虑,r为外表面半径,系数A、B及函数心小的表达式取决于儿何形状,见教材表3-2所示。 3.正规阶段状况的实用计算方法 当Fo>0.2时,可采用上述计算公式求得非稳态导热物体的温度场及交换的热量,也可釆用简化的拟合公式和诺模图求得。 1、诺模图: 工程技术中,为便于计算,采用按分析解的级数第一项绘制的一些图线,叫诺模图。 2、海斯勒图: 诺模图中用以确定温度分布的图线,称海斯勒图。 首先根据(3—13)式给出第何随Fo及Bi变化的曲线(此时x/5二0),然后根据(3-14)式确定9/險的值,于是平板中任意一点的却佝值? 耳旦 便为: 现外%(3-21) 同样,从初始时刻到时刻t物体与环境间所交换的热量,可釆用(3—15)、 (3—17)作岀曲线。 3、诺模图法评述 优点: 简洁方便。 缺点: 准确度有限,误差较大。 U前,随着计算技术的发展,直接应用分析解及简化拟合公式计算的方法受到重视。 4.分析解应用范围的推广及讨论 1、推广范围 1)对物体被冷却的情况也适用; 2)也适于一侧绝热,另一侧为第三类边界条件的厚为5的平板; 3)当固体表面与流体间的表面传热系数hts时,即表面换热热阻T0时,所以欣-8时分析解就是固体表面温度发生一突然变化然后保持不变时的解,即第一类边界条件的解。 2、讨论Bi与Fo对温度场的影响: 1)傅立叶数Fo: 由(3-10)、(3-13)式及诺模图可知: 物体中各点的过余温度随时间t的 增加而减小;而Fo与成正比,所以物体中各点过余温度亦随Fo的增大而减小。 2)毕渥数Bi Bi对温度的影响从以下两方面分析: —方面,从教材图3—6可知,Fo相同时,Bi越大,&妙/%越小。 因为,Bi越大,意味着固体表面的换热条件越强,导致物体的中心温度越迅速地接近周围介质的温度;当BiTo□时,意味着在过程开始瞬间物体表面温度就达到介质温度,物体中心温度变化最快,所以在诺模图中1/Bi=0时的线就是壁面温度保持恒定的第一类边界条件的解。 另一方面Bi的大小决定于物体内部温度的扯平程度。 如: 对于平板,从诺模图3—7中可知: 丄 当焉>10(即Bi<0.1)时,截面上的过余温度差小于5% 当Bi下限一直推到0・01时,其分析解与集总参数法的解相差极微。 综上可得如下结论: 介质温度恒定的第三类边界条件下的分析解;当BiTeo时,转化为第一类边界条件下的解,Bit0时,则与集总参数法的解相同。 §3-4二维及三维非稳态导热问题的求解 1.求解方法 、对于典型的儿何形状的物体,可利用一维非稳态导热问题分析解的组合求 得。 如图3-9所示: 无限长方柱体的非稳态导热问题,属二维导热问题。 截面 尺寸为: 2"广2旳的方柱体可视为两块厚度分别为2用及2勺的无限大平板垂直相交所截出的物体。 讨论的U的: 找出二维温度场与两块无限大平板的温度场之间的关系。 已知: 方柱体初温为to,初始时放于t流体中,表面传热系数为h・试求: 温度场分布。 解: 如图3-9所示,建立坐标系,山于其对称性,只研究其U4截面的温度分布, (3-2 截面上的温度分布山下列导热微分方程和定解条件确定: —“升=0 必dya (II) (III) (IV) e=心尹力-匕=& 式中: 1_殆一爲为无量纲过余温度。 如果无量纲过余温度与©,3力分别是处于与方柱体同样定解条件 下的厚度分别为2勺及2®的无限大平板的分析解,则它们必须满足各自的导热微分方程及定解条件,即: 饥ax2 (3-23) (VI) ©2力+学竺倍 hOx (VII) 及 (3-24) ◎30)=1 Cy,T) 旧页飞― 只要证明: 两块无限大平板分析解的乘积就是上述无限长方柱体的分析解,即: =®x(,2C.T)®yO•巧 (3-25) 证明: 首先证明式(3-25)满足导热微分方程(3-22),为此将式(3-25)代入式(3-22)的左右两端得: 旳_冷丿)_◎氾丿十⑥◎ 左端: drdrxQtydr .・・证明职(“)满足微分方程 其次证明: 满足初始条件 根据及即3弋)的初始条件®x(^o)=1和630)=1得: ®x(3®>>0.0)=1 ・•・证明⑧心.诈2©满足初始条件*.”0)=1 最后证明: 耳(")即3仍=满足边界条件 将式(3-25)代入边界条件(I),并注意到式(VII)的关系得 兄短X(x,£) Ad®x(X,T) 二OyCy,T)[◎”內工)+》Z]二07^r)*O=O 同样可以证明它也满足式(3-34)o再将式(3-25)代入边界条件(III),并注意到式(VI)的关系 得 39(x,t) -Fl^-°.巧二0*^y(y.t)=0 同理可证明它也满足式(IV)。 .・・证明职(“)满足边界条件 综上可知: 职穴生矩八严⑺是上述无限长方柱体导热微分方程的解。 结论: 此方法是多维非稳态导热的求乘积解法,此法适用于第一类边界条件,且£。 =const时o 同理,对长圆柱体,矩六柱体等二维,三维非稳态导热问题,可以用相应的二个或三个一维问题的解的乘积来表示其温度分布。 一—表示无限大平板的解 ——表示无限长圆柱体的解 则: (a): t)]A[®o,r)]A (b): ㊀0・兀巧=巧]pgg巧l (c).&(x.y.s.T)=[®(^r)iA[®o,r)]A 二、乘积解法的适用条件 •初始温度为常数,zo=const; •第一类边界条件,”=const; •第三类边界条件,心=const、h=const; •线性微分方程,且定解条件均为齐次,即乘积解中温度必须以过余温度或无量纲过余温度的形式表示。 说明: 对于形状复杂或边界条件复杂,分析解法无能为力,应借助其它的求解的方法, 如①数值解法;②实验模拟法。 §3-5半无限大物体的非稳态导热 一.半无限大物体的概念 儿何上是指从x二0的界面开始可以向正的x方向及其他两个坐标(x,y)方向无限延伸的物体,称半无限大物体。 实际中不存在该物体,但研究物体中非稳态导热的初始阶段,可把实物看为该物体处理。 如: 有限厚度的平板,起初有均匀温度,后其侧表面突然受到热扰动,如 •壁温突然升高到一定值并保持不变; •壁面突然受到恒定的热流量密度加热; •壁面受到温度恒定的流体的加热或冷却。 当扰动的影响只局限在表面附近,而尚未进入平板内部时,就可视该平板为,“半无限大”物体。 二.第一类边界条件下半无限大物体非稳态导热温度场的分析解 如图3-11所示: 已知半无限大物体初始温度均匀为(to儿当存0时,X二0侧表面温度突然升高到",并保持不变,试确定物体内温度随时间的变化和在时间间隔[0,T]内的热流量。 解: 1、物体内的温度分布 根据半无限大物体的定义,得出其导热微分方程: dT~adx2(3-26) 初始条件为: T二o时,心Q=S 边界条件为: x二0时,"6巧=5 X—2时,=/0 JI入过余温度: &=“兀,巧则冇: =a—— 無/(3-27) T=0时,。 (厂0)=內=8—S x=0时,°(°Q=S—5=° XfB时9巩8>巧=f0 将微分方程訂F分离变量并求解得分析解为: =erf(2-7^)=erfn(3-28) X 其中: 无量纲变量Q=玮;erfQ称为误差函数,它随Q的变化而变化,由附录表可知: 0X 当n二2时,6=0.9953,就是说当n>2即2后22时,该处x的温度仍认为等于to(无量纲过余温度的变化小于5%),由此得到以下两个重 要参数: 1从儿何位置上说,若心而,则时刻t时x处的温度可认为未发生变化。 所以,对^=0且厚为28的平板,当其一侧温度突然变化到另一恒定温度时,若8"后,则在t时刻之前•该平板中瞬时温度场的汁算可采用半无限大物体模型处理。 2从时间上看,如果「応,则此时x处的温度可认为完全不变,所以把石视 兀2 为tfi性时间,即X时x处的温度可认为仍等于叽 2、表面上的瞬时热流密度及在[0,t]时间间隔内放出或吸收的热量: 物体中任意一点的热流密度®"畤八恥f報如 =兄空二匕汀,久也巧 (3-29) 则,表面上的热流密度为 曲此可见: ①半无限大物体在第一类边界条件影响下被加热或冷却时,界面上的瞬时热流量与时间的平方根成反比; ②在时间[0,J内交换的总热量则正比于何及时间的平方根。 其中: 何称为吸热系数,表示物体与其接触的高温物体吸热的能力。 三、半无限大物体概念的适用范围 只适于物体非稳态导热的初始阶段,当物体表面上的热扰动已深入传递到物体内部时,就不再适用,则应釆用前述分析方法。 补充: 非稳态导热问题求解思路 •解一维非稳态导热问题的基本思路 1)首先,用Bi检验是否满足集总参数法的条件,若性质属于h或5未知,可先假设,然后校核; 2)若不能用集总参数法,可采用分析解法{诺模图法和近似公式法}: 3)若2),1)方法均不能求解,则釆用数值解法。 •多维非稳态导热问题的求解方法 1)是否满足乘积解法的条件; 2)合理将一个多维问题分析成儿个一维问题。 二.物体内速度变化的规律 •温度变化最慢的点位于物体的体心或形心; •温度变化最快的点位于离物体的体心或形心最远处。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 第三 稳态 导热 分析 解法