传热学试验—-墙角matlab导热问题

1、二维导热物体温度场的数值模拟姓名小明学号 111111班级 能动学院能动问题描述有一墙角模型，尺寸如图1所示，导热系数（m K），墙角内外壁 为第一类边界条件。求解该模型的温度分布及导热量。二、计算原理根据热平衡法列出节点方程，各方向导入单元体的热量之和为零。内节点和绝热边界点（图1点划线上的点）的方程形式不同内节点维熱.边界图2所示的内节点和绝热边界节点方程如下:内节点:SEW?1?(ti,j1ti,j)X(Sti,j)x(tj1 lj) y(S,jti,j)yyyxx绝热边界点:?1?(lj i Ij) x(1门 Ij) x o (tj i,j Ij) y2 y2 yx三、计算过程用语言编

2、写计算程序，初取网格步长X y 0.1mfun亡tion I.for i=l:16Ki, 1)=30;endfot j=l:12T(I, j) = 30;endfor jn=6:16B for n=6:12I n) = 0;endendT (16, 2:5)= 24,18, 12,6;1(2:5, 12)= 24, 18, 12, 6: while 1 fot i=2:15 for j=2:5P3=1；+ *iap=T(i, j);I Ci, j)=TH- L j)+T (i+L j)+T (i, j-l)-hT (if j+l)/4;5ps=I(lj j)-tenp,k end卜endif

3、?ps+T(16, j+l)/4,epsKlSj j.)-t:endif &psle-8brak;endindnIHTUwhile Ifor i=2:5f cr j=2:11ep5-1:temp=T (ij j);I(i, j)= (T(i-lj j)+T(i十 1* j)十T(ij j-l)+l (i, j+1) )/4; eps= I (ij, j) -1 enp ,endendif ?ps ；I(i, 121 = (2*1(1, lD+TCi-L 12)+T(i+1,12) )/4 ；eps=T t.if 12)七ejiip ;*ndif epsle-8break;endendpcol&

4、r (IJ ； shading int erp-ligurej contour1 I , 10);运行结果:30.0030,0030.0030.0030030.0030,0031.530.0030.0030.0030. CO30. 0028,9327,8926.9006瓯h26.0225.3021792 A324.0624.0530. 002T.K525. 223. DZ1.S720.3919.3918.7913.2213.0918.0830. 0026. T223.4?20.29r ：7. ?&：5.0D13.6212.8612.6：2. 2312.1012. OS30, 0025.5921

5、,1316,64LI 30L637,21&596.306.156,066,0530- 00乩鸭18,3212,35rO.(CO.QD0.00O.OD0400.00O.OD30. 0024.27IS. 45r 12.46r 6.290.0：:i.(n:.：DD. P0.0030.0024.： 518.2512.25r 6.i&OXD0,00(HD0.000.000.000.U030.0024,0818, ：3:H6,080.0D0.000.000.0D0.000.Q0Q.0030- 0024.0518.0712,0?6,5OAC?(O:AO0.:0.0030. 0024.02IE. 0412.0

6、46.020.(1：0.0D0.00O.OD0.000.00O.OD30.0。24.0118.0212,02r 6,oi0,0D0,000,000.000.00CL 000,0030. 00201IE, 0112,01r 6,010.000.000.00O.OD0.000.00O.OD30. 00r?4.oor is.oir i2,oi6.0D0.00MO?(00：D0.0 J0.0030. 0024.00IE. 0012.00F 6.(1：0.(1：0.000.000.000.0u0.000.0030, O024.00LS.O012,00r 6.0D0,000,000,000.000.00CL 000,00图1各个点的温度数值图2:分层设色等温线分布I-20-15-10I24681012 U 16长度图3:等温线分布（每两条线间隔为三度）20406060100120140长度- 15四、小结本次数值模拟是运用 matlab 程序用于数值计算。小组成员共同