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

《传热学MATLAB温度分布大作业完整版.pdf》由会员分享，可在线阅读，更多相关《传热学MATLAB温度分布大作业完整版.pdf（11页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、东南大学能源与环境学院东南大学能源与环境学院课程作业报告课程作业报告作业名称作业名称：传热学大作业利用传热学大作业利用 matlabmatlab 程序解决热传导问题程序解决热传导问题院系：能源与环境学院院系：能源与环境学院专业：建筑环境与设备工程专业：建筑环境与设备工程学号：学号：姓名：姓名：20142014 年年 1111月月 9 9 日日请浏览后下载，资料供参考，期待您的好评与关注！一、题目及要求一、题目及要求1.2.3.4.5.6.7.原始题目及要求各节点的离散化的代数方程源程序不同初值时的收敛快慢上下边界的热流量（=1W/(m）计算结果的等温线图计算小结题目：已知条件如下图所示：二、各

2、节点的离散化的代数方程二、各节点的离散化的代数方程各温度节点的代数方程ta=(300+b+e)/4;tb=(200+a+c+f)/4;tc=(200+b+d+g)/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+m+j)/4;tj=(f+i+k+n)/4;tk=(g+j+l+o)/4;tl=(2*k+h+q)/4请浏览后下载，资料供参考，期待您的好评与关注！tm=(2*i+300+n)/24;tn=(2*j+m+p+200)/24;to=(2*k+p+n+200)/

3、24;tp=(l+o+100)/12三、源程序三、源程序【G-S迭代程序】【方法一】函数文件为：function y,n=gauseidel(A,b,x0,eps)D=diag(diag(A);L=-tril(A,-1);U=-triu(A,1);G=(D-L)U;f=(D-L)b;y=G*x0+f;n=1;while norm(y-x0)=eps x0=y;y=G*x0+f;n=n+1;end命令文件为：A=4,-1,0,0,-1,0,0,0,0,0,0,0,0,0,0,0;-1,4,-1,0,0,-1,0,0,0,0,0,0,0,0,0,0;0,-1,4,-1,0,0,-1,0,0,0,0

4、,0,0,0,0,0;请浏览后下载，资料供参考，期待您的好评与关注！0,0,-2,4,0,0,0,-1,0,0,0,0,0,0,0,0;-1,0,0,0,4,-1,0,0,-1,0,0,0,0,0,0,0;0,-1,0,0,-1,4,-1,0,0,-1,0,0,0,0,0,0;0,0,-1,0,0,-1,4,-1,0,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,-1,0,-1,0,4,0,0,0,-1,0,0,0;0,0,0,0,0,-1,0,0,-1,4,-1,0,0,-1,0,0;0,0,0,0,0,0,-1,0,0,

5、-1,4,-1,0,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,-2,0,0,0,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,-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,100,0,0,0,100,0,0,0,300,200,200,100;x,n=gauseidel(A,b,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1.

6、0e-6)xx=1:1:4;yy=xx;X,Y=meshgrid(xx,yy);Z=reshape(x,4,4);Z=Zcontour(X,Y,Z,30)Z=139.6088 150.3312 153.0517 153.5639请浏览后下载，资料供参考，期待您的好评与关注！108.1040 108.6641 108.3119 108.1523 84.1429 67.9096 63.3793 62.4214 20.1557 15.4521 14.8744 14.7746【方法2】t=zeros(5,5);t(1,1)=100;t(1,2)=100;t(1,3)=100;t(1,4)=100;t(

7、1,5)=100;t(2,1)=200;t(3,1)=200;t(4,1)=200;t(5,1)=200;for i=1:10t(2,2)=(300+t(3,2)+t(2,3)/4;t(3,2)=(200+t(2,2)+t(4,2)+t(3,3)/4;t(4,2)=(200+t(3,2)+t(5,2)+t(4,3)/4;t(5,2)=(2*t(4,2)+200+t(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)/

8、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(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(5,4)=(2*t(4,4)+t(5,3)+t(5,5)/4;t(2,5)=(2*t(2,4)+300+t(3,5)/24;t(3,5)=(2*t(3,4)+t(2,5)+t(4,5)+200)/24;t(4,5)=(2*t(4,4)+t(3,5)+t(5,

9、5)+200)/24;t(5,5)=(t(5,4)+t(4,5)+100)/12;tendcontour(t,50);ans=100.0000 200.0000 200.0000 200.0000 200.0000 100.0000 136.8905 146.9674 149.8587 150.7444 100.0000 102.3012 103.2880 103.8632 104.3496 100.0000 70.6264 61.9465 59.8018 59.6008 100.0000 19.0033 14.8903 14.5393 14.5117请浏览后下载，资料供参考，期待您的好评与关

10、注！【Jacobi迭代程序】函数文件为：function y,n=jacobi(A,b,x0,eps)D=diag(diag(A);L=-tril(A,-1);U=-triu(A,1);B=D(L+U);f=Db;y=B*x0+f;n=1;while norm(y-x0)=eps x0=y;y=B*x0+f;n=n+1;end命令文件为：A=4,-1,0,0,-1,0,0,0,0,0,0,0,0,0,0,0;-1,4,-1,0,0,-1,0,0,0,0,0,0,0,0,0,0;0,-1,4,-1,0,0,-1,0,0,0,0,0,0,0,0,0;0,0,-2,4,0,0,0,-1,0,0,0,

11、0,0,0,0,0;-1,0,0,0,4,-1,0,0,-1,0,0,0,0,0,0,0;0,-1,0,0,-1,4,-1,0,0,-1,0,0,0,0,0,0;0,0,-1,0,0,-1,4,-1,0,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,-1,0,-1,0,4,0,0,0,-1,0,0,0;0,0,0,0,0,-1,0,0,-1,4,-1,0,0,-1,0,0;0,0,0,0,0,0,-1,0,0,-1,4,-1,0,0,-1,0;0,0,0,0,0,0,0,-1,0

12、,0,-2,4,0,0,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,-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,100,0,0,0,100,0,0,0,300,200,200,100;x,n=jacobi(A,b,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1.0e-6);xx=1:1:4;yy=xx;X,Y=meshgrid(xx,y

13、y);Z=reshape(x,4,4);Z=Zcontour(X,Y,Z,30)n =97Z=139.6088 150.3312 153.0517 153.5639 108.1040 108.6641 108.3119 108.1523 84.1429 67.9096 63.3793 62.4214 20.1557 15.4521 14.8744 14.7746请浏览后下载，资料供参考，期待您的好评与关注！四、不同初值时的收敛快慢四、不同初值时的收敛快慢1、方法1在Gauss迭代和Jacobi迭代中，本程序应用的收敛条件均为norm(y-x0)=eps，即使前后所求误差达到e的-6次方时，跳出

14、循环得出结果。将误差改为0.01时，只需迭代25次，如下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=1:10判断收敛，为迭代10次，若改为1:20，则迭代20次。2、在同样的误差要求下，误差控制在e的-6次方内，Gauss迭代用了49次达到要求，而Jacobi迭代用了97次，可见，在迭代中尽量采用最新值，可以大幅度的减少迭代次数，迭代过程收敛快一些。在Gauss中，初值为100，迭代46次达到精确度1.0e-6,初值为5

15、0时，迭代47次，初值为0时，迭代49次，初值为200时迭代50次，可见存在一个最佳初始值，是迭代最快。这一点在jacobi迭代中表现的尤为明显。五、上下边界的热流量：五、上下边界的热流量：上边界t=200,t=10，所以，热流量1=*200-tb200-tc200-tdx200-100 x200-tax+x+x+*+yyyy2y2=1*（100/2+(200-139.6088)+(200-150.3312)+(200-153.0517)+(200-153.5639)/2)=230.2264W下边界tl-tqxtj-totk-tpti-tmx+热流量2=|*-*x+x+yy2yyh*(t-tx

16、t-t100-10 xtn-tt-t*x+o*x+m*x+p*+)|y2yyyy2请浏览后下载，资料供参考，期待您的好评与关注！=|1*(84.1429-20.1557)+(67.9096-15.4521)+(63.3793-14.8744)+(62.4214-14.7746)/2)-10*(90/2+(20.1557-10)+(15.4521-10)+(14.8744-10)+(14.7746-10)/2)|=|-489.925|W=489.25W六、温度等值线六、温度等值线Gauss:Yacobi:请浏览后下载，资料供参考，期待您的好评与关注！七、计算小结七、计算小结导热问题进行有限差分数

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