传热学MATLAB温度分布大作业完整版.docx

1、传热学MATLAB温度分布大作业完整版东南大学能源与环境学院课程作业报告作业名称：传热学大作业利用matlab程序解决热传导问院系：能源与环境学院专业：建筑环境与设备工程学号：姓名：2014年11月9日一、题目及要求1.原始题目及要求2.各节点的离散化的代数方程3.源程序4.不同初值时的收敛快慢5.上下边界的热流量(入=lW/(mC)6.计算结果的等温线图7.计算小结题目：已知条件如下图所示：Tf=10, X=lW/(nr4C)h=10W/ (m2 )二、各节点的离散化的代数方程各温度节点的代数方程ta=(300+b+e)/4 ; tb=(200+a+c+f)/4;tc=(200+b+d+g)

2、/4; td=(2*c+200+h)/4te=(100+a+f+i)/4; tf=(b+e+g+j)/4;tg=(c+f+h+k)/4 ; th=(2*g+d+l)/4ti= 100+e-rm+j ) /4; tj= (f+i+k+n) /4;tk=(g+j+l+o)/4; tl=(2*k+h+q)/4tm=(2*i+300+n)/24; tn=(2*j+m+p+200)/24; to=(2*k+p+n+200)/24; tp=(l+o+lOO)/12三、源程序【G-S迭代程序】【方法一】函数文件为：function yr n=gauseidel(A,bz xQreps)D=diag(diag

3、(A);L=-tril(A,-1);U=-triu(A,1);G=(D-L)U;f=(D-L)b;y=G*xO+f;n=l;while norm(y-xO)=epsxO=y;y=G*xO+f;n=n+l;end命令文件为：公=4,-1,0,0,-1,0,0,0,0,0,0,0,0,0,0,0;-1,4,-1,0z 0,-1,0,09 0z 0f 0z 0f 0,0,0f 0;0,1,4,1,0,0,1,0z 0z 0,0f 0f 0,0,0,0;0,0,-2,4,0,0,0,-1,0,0,0,0,0,0,0,0;1,09 0f0,4,1,0r 0,1,0,0r 0r 0f0,0,0;0,-1,

4、0f 0,-1,4,-1,0z 0,-1,0,0z 0z 0/0,0;0,0,-1,0,0,-1,4,-1,0z 0,-1,0/0/0,0,0;0,0,0,-1,0,0,-2,4,0,0,0,-1,0,0,0,0;0,0,0,0,-0,4,0,0,0,-1,0,0,0;0,0,0,0z 0,-1,0,0,-1,4/-1/0r 0,-1,0/0;0,0,0,0/0z 0,-1,0,0,0f 0,-1,0;0,0,0,0/0/0f 0,-1,0/O7-2,4,0z 09 0,-1;0,0,0,0,0,0,0,0,-2,0,0,0,24,-1,0,0;0/0,0,0,0,0,0,0,0,-2,0,0

5、,-1,24,-1,0;0,0,0,0/0z 0,0z 0f 0f 0,-2,0f 0,-0,0,0,0,0,0,0,0,0,0,0,-1,0,0,-1,12;b=300,200,200,200,10 0,0,0, 0, 100,0, 0, 0,300, 200,200,100;xz n=gauseidel(A,b, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0 ,，l.0e-6) xx=l:1:4;yy=xx;XrY=meshgrid(xxz yy);Z=reshape(x,4,4);Z=Zfcontour (X, Y, Z, 30)Z =139.6088 150.3312

6、 153.0517 153.5639108.1040 108.6641 108.3119 108.152384.1429 67.9096 63.3793 62.421420.1557 15.4521 14.8744 14.7746【方法2】 t=zeros (5,5);t(1,1)=100;t(1,2)=100;t(I,3)=100;t(lz4)=100;t(1,5)=100;t(2,1)=200;t(3,1)=200;t(4zl)=200;t(5,l)=200;for i=l:10t (2,2) = (300+t (3, 2) +t (2Z 3) ) /4 ;t(3z2) = (200+t(

7、2,2)+t(4z2)+t (3,3)/4;t(4,2) = (200+t (3,2)+t (5,2)+t (4,3) )/4;t(5z2)=(2*t(4,2)+200+1(5,3)/4;t(2,3) = (100+t(2,2)+t (3,3)+t (2,4)/4;t(3,3) = (t(3,2)+t(2,3)+t(4,3)+t(3,4) )/4;t(4,3) = (t(4,2)+t(3,3)+t(5,3)+t(4,4)/4;t(5,3) = (2*t (4,3)+t (5,2)+t (5,4) )/4;t(2,4) = (100+t (2,3)+t (2,5)+t (3,4) ) /4;t(

8、3,4) = (t(3,3)+t(2,4)+t(4,4)+t(3,5) )/4;t(4,4) = (t(4,3)+t(4,5)+t(3,4)+t(5,4) )/4;t(5z4) = (2*t(4z4)+t(5z3)+t(5f 5)/4;t(2,5) = (2*t(2f4)+300+t (3,5)/24;t (3z5) = (2*t (3z4)+t(2,5)+t (4r 5)+200) /24;t(4r 5) = (2*t (4z4)+t (3, 5)+55,5)+200)/24;t(5,5)=(t(5,4)+t(4,5)+100)/12;end contour (t,50);ans =【Jac

9、obi迭代程序】函数文件为：function yz n=jacobi(A,b,xOz eps)D=diag(diag(A);L=-tril(A,-1);U=-triu(A,1);B=D (L+U);f=Db;y=B*xO+f;n=l;while norm(y-xO)=epsxO=y;y=B*xO+f;n=n+l;end命令文件为：A=4,-l,0,0,-l,0,0,0,0,0,0,0,0,0,0,0;-1,4,-1,0f 0,-1,0,0/0f 0,09 0,0,0,0,0;0,-1,4,-1,0,0,-1,0z 0z 0,0z 0,0,0,0,0;0,0,-2,4,0,0,0,-1,0,0,

10、0,0,0,0,0,0;-1,0,0f 0,4,-1,0,0,-1,0,0/0,0r 0,0,0;0, -1, 0f 0, -1, 4, -1, 0/ 0, -1,0/ 0f 0/ 0/ 0f 0;0r0,-1,0,0,1,4,1,0z 0,1,0z 0/0z 0,0;0,0,0, -1,0,0,-2,4r 0,0,0,-1,0,0,0,0;0,0,0,0,0,4,0,0,0,-1,0,0,0;0/0/0/0/0,-1,0,0,-1,4,-1,0r 0,-1,0,0;0,0,0,0/0z 0,-1,0,0,-0r 0,-1,0;0f0f0z 0/0z 0f 0,-1,0/0z-2/4/0/0,

11、0,-1;0,0,0,0,0,0,0,0,-2,0,0,0,24,-1,0,0;0,0,0,0z 0z 0/0f 0r0r-2/09 0,-1,24,-1,0;0, 0, 0, 0,0, 0, 0, 0,0,0, -2, 0,0, -1,24, -1;0,0,0,0,0,0,0,0,0,0,0,-1,0,0,-1,12;b= 300,200,200,200, 10 0,0,0, 0, 100,0, 0, 0,300,200,2 00,100;xr n=jacobi (A,b, 0, 0,0, 0,0, 0,0,0,0, 0,0,0,0,0,0,0 11 . Oe-6);xx=l:1:4;Xz

12、Y=meshgrid(xx,yy);Z=reshape(x,4,4);Z=Z 1contour (X, Y, Z, 30) n =97Z =139.6088108.104084.142920.1557150.3312108.664167.909615.4521153.0517108.311963.379314.8744153.5639108.152362.421414.7746四、不同初值时的收敛快慢1、方法1在Gauss迭代和Jacobi迭代中，本程序应用的收敛条件均为norm (y-xO) =eps, 即使前后所求误差达到e的-6次方时，跳出循环得出结果。将误差改为0.01时，只需迭代25

13、次，如下x, n =gauseidel (A,b, 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.01)运行结果为将误差改为0.1时，需迭代20次，可见随着迭代次数增加，误差减小，变化速度也在减小。方法2通过i=l:10判断收敛，为迭代10次,若改为1:20,则迭代20次。2、在同样的误差要求下，误差控制在e的-6次方内，Gauss迭代用了49次达到要求，而Jacobi 迭代用了 97次，可见，在迭代中尽量采用最新值，可以大幅度的减少迭代次数，迭代过程收敛快 一些。在Gauss中，初值为100,迭代46次达到精确度l.Oe-6,初值为50时，迭代47次，初值为0 时，迭代

14、49次，初值为200时迭代50次，可见存在一个最佳初始值，是迭代最快。这一点在jacobi 迭代中表现的尤为明显。五、上下边界的热流量：上边界t=200,t =10,所以, X200,200-100 .At 200-t, 4热流量中1二、* 木+ LAX +y 2 Ay二 1* (100/2+(200-139.60 88) + (200-150.3312) + (200-153.0517) + (200-153.5639)/2)=230.2264W下边界 热流量中2二人*上5&什口1加+三卫加+口*4-Ay Ay Ay Ay 2Ay 2 Ay Ay Ay )】 2=1*(84.1429-20.

15、1557)+(67.9096-15.4521)+(63.3793-14.8744)+(62.4214-14.774 6)/2)-10*(90/2+(20.1557-10) + (15.4521-10) + (14.8744-10) + (14.7746-10) /2)I = -489.925 W =489. 25W六、温度等值线Gauss:4 一一一 1 1 1 1 T： r一 3.5二 -3 -2.5-2 , 一1 5二二 1 一 一 _ 1: ； ； 1.5 2 2.5 33.54Yacobi:七、计算小结导热问题进行有限差分数值计算的基本思想是把在时间、空间上连续的温度场用有限个离散点温

16、度 的集合来代替，即有限点代替无限点，通过求解根据傅里叶定律和能量守恒两大法则建立关于控制 面内这些节点温度值的代数方程，获得各个离散点上的温度值。要先划分查分网格，在建立差分代数方程组，用MATLAB或者其他软件编程求解。高斯-赛德尔迭代法和雅克比迭代法区别在于使用新植和旧值进行下一次迭代，而采用新值迭代的 高斯-赛德尔迭代收敛的更快些，但其求解代数方程是不一定得到收敛的解，其原因可能由于迭代 方式不合适造成。在计算热流量过程中，主要是正确利用傅里叶定律和牛顿冷却公式，本题中需要特别注意的一点是 后边界是绝热的，因而左右方向上几乎不存在热量的传递，所以看似是二维稳态问题实际上是一维 稳态的问题。求解也比较简单。程序运行出来的等温线结果也很好的说明了这一点，温度总体是从 上向下递减，热量传递方向是自上而下。