《传热学实验—-墙角matlab导热问题》由会员分享，可在线阅读，更多相关《传热学实验—-墙角matlab导热问题（12页珍藏版）》请在金锄头文库上搜索。
1、二维导热物体温度场的数值模拟姓名小明学号111111班级能动学院能动问题描述有一墙角模型，尺寸如图1所示，导热系数（mK）,墙角内外壁为第一类边界条件。求解该模型的温度分布及导热量。二、计算原理根据热平衡法列出节点方程，各方向导入单元体的热量之和为零。内节点和绝热边界点（图1点划线上的点）的方程形式不同。内节点绝热边界图2所示的内节点和绝热边界节点方程如下:内节点：储1t,1X(ti,j1t,)X(tiljtpy(ti1,jt,)yNSEW?1?yyxx绝热边界点:?1?(ti,j1ti,j)X(ti,j15X0(ti1,j5y2y2yx三、计算过程用语言编写计算程序，初取网格步长y0.1mf
2、unctionI.fori=l:出1)=30；endfotj=l:12T(I,j)=30;卜endforjn=6:16Bforn=6:12加n)=0endend1(16,2:5)=24,18,12,6：1(2:5,12)=24,18,12,6;while1中fori=2:15forj=2:5P3=1；+*inp=T(i,j.;ICi,j)=TH-Lj)+T(i+Lj)+T(i,(ifj+l)/4;eps=Ij)-tenp,卜endendif?p+T(16,j+l)/4,eps=l(16;endifeps1e-8break;nTHTEendindwhileIfari=2:5-forj=2:11e
3、ps-1；teinp=T(i,j);ICi,j)=(T(i-Lj)+T(itl,j)+1(1,j-D+Ki,j+l)/4;eps=I(i5j)-teup;endendif?psle-Sbrealt;*ndendwhil*1臼fori=2:5lentp=r(ij12);I(i,12)=(2*1(1,1D+TC1-L12)+T(i+1,12)/4；eps=Tt.if12JtejiLjj;*ndifeps00.000.0D0.00QJQ0.0030-0024.0518.0712,0?6,5口：二的:1内0.：0：的0.0：30.0024.021E,0412.046.020.(1：0.0D0.000.
4、0D0.000.000.0030.0024.0118.0212,02r6,010,0D0,000,000.000.00CLOO0,(030.001E,0112,01r6,010.000.000.000.000.000.000.0030.00r?4.oor18.01r12,oi6.0DCODMO二10：.：00.：00.0J0.0：30.0024.00IE.0012.00C6.(1：0.(1：0.000.000.00O.0IJ0.000.0030,0024.00LS.O012.00r6.0D0,000,30,000.000.00CLOO0,(0图1:各个点的温度数值121110337554321246610121415图2：分层设色等温线分布蛔图3:等温线分布（每两条线间隔为三度）252015W20406060100120140长度15长度四、小结本次数值模拟是运用matlab程序用于数值计算。小组成员共同讨论并复习了热传导问题的数学描述和热平衡法；从模拟过程中练习了不同节点迭代方程的建立；并简单学习了matlab语言的使用。这次大作业对于我们以后的学习和可能的研究来说是一个很好的锻炼机会。
