传热学实验--墙角matlab导热问题

1、二维导热物体温度场的数值模拟姓名小明学号111111班级能动学院能动一、问题描述有一墙角模型，尺寸如图1所示，导热系数(mK),墙角内外壁为第一类边界条件。求解该模型的温度分布及导热量。图1F11L1二、计算原理根据热平衡法列出节点方程，各方向导入单元体的热量之和为零。内节点和绝热边界点(图1点划线上的点)的方程形式不同。图2图2所示的内节点和绝热边界节点方程如下内节点：+=九NSEW=0TOC o 1-5 h z(tt)Ax(tt)Ax(tt)/y(tt)/y*v*-*v*v*v*ij1I,J+ij1I,J+i+1jI,J+i1jI,JAyAy&绝热边界点：+=NSEW(tt)Ax(tt)A

2、x(tji,j+ji,j+0+i1Jt)Ayi,jAx=0三、计算过程用语言编写计算程序，初取网格步长Ax二Ay二0.1m5HTafunct:L=M;endform=6;16forn=6:12T皿nj=0:end.end1(16,2:5)=24,18,12.6;1(2:5,12)=24,18,12,6;vhile1for1=2:15forj=2:5；temp=I-ijj)；心jJ=(T(i-Lj)+T(1+LjHT(i,j-l)十心j+lD/4；eps=Iijj)-temp:end.endrwr?ifeps1e8bxsak:endendwhile1tenip=rr16j);1(16,j)-Q*

3、T(1.5,ji+TClG,j-D+Ttie,j+l)/4;eps=T16jj-temi);-endife-psle-8break;endendwhile1for1=2:5forj=2:11二:temp=T(ijj);龙jJ+rCi+UjJ+KLj-D+Tti,j+L)/4;epIijj)te皿p;end-endifepsle-8break;end-end-vhile1fori=2:5teHp=r(i.j12J；Ki,L2)=Q*T(l,llHKi-L12)+T(i+1,12)/4;eps=Tij12一七emp;?n.difepsleSbreak;end-endpcolorI1:shading

4、interpfigurejcontour(TJlu);运行结果：jO.O：jO.O：30.0：30.0：30.0：3C.0：3：.(03：.CD3：.?)3：.?)3：.(O30.0：jO.O：28.9:27.8=26.900526.9：26.02盯门2.79胆.邪24.13乩0524.0:jO.O：r茫25.2瓯：21.S72C.3；三门IE.79IE.3IE.22IE.0318.OEjO.O：26.223.4?20.217.3:1E.0：1-.6212.E6二労12.23二1)12.OEjO.O：药5221.16.612.3：E.6j7.216.59c.ZD5.15c.056.05jO.O

5、：24.4918.812.85F，用C.O：:.CD：.(O0.0：jO.O：24.2718.4512.466.2=：.(O:.CD:0.0：jO.O：14二18.2512.2c6.1:C.O：.(O:.CD：.(O0.0：jO.O：24.OE18.匸二6.0E：.(O:.CD2/0:）0.0：jO.O：24.0518.0?12.0?6.0:：:.CD：.：D2/0:）0.0：jO.O：24.0218.0412.06.02C.O：.(O:.CD：.(O0.0：jO.O：24.0118.012.026.01C.O：:.CD：.(O0.0：jO.O：24.0118.0112.016.01：.(O:.CD:0.0：jO.O：24.0：18.0112.016.0：C.O：.(O:.CD：.(O0.0：jO.O：24.0：18.0212.0：6.0：.(O:.CD2/0:）0.0：jO.O：24.0：18.0：12.0：6.0：:.CD：.：D2/0:）0.0：图1：各个点的温度数值图2：分层设色等温线分布图3：等温线分布（每两条线间隔为三度）四、小结本次数值模拟是运用matlab程序用于数值计算。小组成员共同讨论并复