哈工程传热大作业

传热学大作业班级：20121515学号：2012151531姓名：张永宽第一题：如图所示，一个无限长矩形柱体，其横截面的边长分别为L1和L2，常物性。该问题可视为二维稳态导热问题，边界条件如图中所示，

传热学大作业 20121515 班级： 2012151531 学号： 姓名：张永宽 第一题： LL 如图所示，一个无限长矩形柱体，其横截面的边长分别为和，常物性。 12 L =0 该问题可视为二维稳态导热问题，边界条件如图中所示，其中。 1 LTTλ 6m,=0.4m =60,=20=200W/m·K) ，（ ，℃℃。 2w1w2 1) . （编写程序求解二维导热方程 xLyL (2 =/2=/2 ）绘制和处的温度场，并与解析解进行比较。已知矩形内 12 的温度场的解析解为。 (1 ） 164b)Tmn=Tm+1n+Tm-1n+Tmn+1+Tm,n-1/4; 1cm. 根据课本页公式（，（，，，）取步长为 编出以下程序迭代求解内部个点温度。 a=zeros4161; %4160 （，）生成＊的矩阵。 k=060; ： a(41,:=20*sinpik/60) Tw2=sin(pi ）（。＊；％矩形上边温度满足＊ x/L1. ） a=a+60; ％使四周都为给定的边界条件。 for x=1:10000 %10000( 迭代次估计能满足要求精度）。 for i=240 ： for j=260 ： ai,j=(ai—1j)+a(ij-1)+a(i+1,j)+aij+1)/4; （）（，，（，） %. 内部每一个点都为周围四个点温度和的四分之一 end end end mesh(a ） title'’,'Fontsize’,18) （第一题（张永宽作请勿抄袭） xlabel'xcm’'Fontsize',14) （轴张永宽作请勿抄袭，单位， ylabel’y,cm'’Fontsize',14) （轴单位， zlabel(’t℃',’Fontsize’,14) 轴，单位