传热学MATLAB温度分布大作业完整版.docx
《传热学MATLAB温度分布大作业完整版.docx》由会员分享，可在线阅读，更多相关《传热学MATLAB温度分布大作业完整版.docx（10页珍藏版）》请在冰豆网上搜索。
传热学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）/4;td=（2*c+200+h）/4
te=（100+a+f+i）/4;tf=（b+e+g+j）/4;
tg=（c+f+h+k）/4;th=（2*g+d+l）/4
ti={100+e-rm+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）/24;tp=（l+o+lOO）/12
三、源程序
【G-S迭代程序】
【方法一】
函数文件为：
function[yrn]=gauseidel（A,bzxQreps）
D=diag（diag（A））;
L=-tril（A,-1）;
U=-triu（A,1）;
G=（D-L）\U;
f=（D-L）\b;
y=G*xO+f;
n=l;
whilenorm（y-xO）>=eps
xO=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,0z0,-1,0,090z0f0z0f0,0,0f0;
0,—1,4,—1,0,0,—1,0z0z0,0f0f0,0,0,0;
0,0,-2,4,0,0,0,-1,0,0,0,0,0,0,0,0;
—1,090f0,4,—1,0r0,—1,0,0r0r0f0,0,0;
0,-1,0f0,-1,4,-1,0z0,-1,0,0z0z0/0,0;
0,0,-1,0,0,-1,4,-1,0z0,-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,0z0,-1,0,0,-1,4/-1/0r0,-1,0/0;
0,0,0,0/0z0,-1,0,0,0f0,-1,0;
0,0,0,0/0/0f0,-1,0/O7-2,4,0z090,-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/0z0,0z0f0f0,-2,0f0,-
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]';
[xzn]=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（xxzyy）;
Z=reshape（x,4,4）;
Z=Zf
contour（X,Y,Z,30）
Z=
139.6088150.3312153.0517153.5639
108.1040108.6641108.3119108.1523
84.142967.909663.379362.4214
20.155715.452114.874414.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;
fori=l:
10
t（2,2）=（300+t（3,2）+t（2Z3））/4;
t（3z2）=（200+t（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（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（5f5））/4;
t（2,5）=（2*t（2f4）+300+t（3,5））/24;
t（3z5）=（2*t（3z4）+t（2,5）+t（4r5）+200）/24;
t（4r5）=（2*t（4z4）+t（3,5）+55,5）+200）/24;
t（5,5）=（t（5,4）+t（4,5）+100）/12;
endcontour（t,,50）;
ans=
【Jacobi迭代程序】
函数文件为：
function[yzn]=jacobi（A,b,xOzeps）
D=diag（diag（A））;
L=-tril（A,-1）;
U=-triu（A,1）;
B=D\（L+U）;
f=D\b;
y=B*xO+f;
n=l;
whilenorm（y-xO）>=eps
xO=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,0f0,-1,0,0/0f0,090,0,0,0,0;
0,-1,4,-1,0,0,-1,0z0z0,0z0,0,0,0,0;
0,0,-2,4,0,0,0,-1,0,0,0,0,0,0,0,0;
-1,0,0f0,4,-1,0,0,-1,0,0/0,0r0,0,0;
0,-1,0f0,-1,4,-1,0/0,-1,0/0f0/0/0f0;
0r0,-1,0,0,—1,4,—1,0z0,—1,0z0/0z0,0;
0,0,0,-1,0,0,-2,4r0,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,0r0,-1,0,0;
0,0,0,0/0z0,-1,0,0,-0r0,-1,0;
0f0f0z0/0z0f0,-1,0/0z-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,0z0z0/0f0r0r-2/090,-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]';
[xrn]=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;
[XzY]=meshgrid（xx,yy）;
Z=reshape（x,4,4）;
Z=Z1
contour（X,Y,Z,30）n=97
Z=
139.6088
108.1040
84.1429
20.1557
150.3312
108.6641
67.9096
15.4521
153.0517
108.3119
63.3793
14.8744
153.5639
108.1523
62.4214
14.7746
四、不同初值时的收敛快慢
1、[方法1]在Gauss迭代和Jacobi迭代中，本程序应用的收敛条件均为norm（y-xO）>=eps,即使前后所求误差达到e的-6次方时，跳出循环得出结果。
将误差改为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=l:
10判断收敛，为迭代10次,若改为1:
20,则迭代20次。
2、在同样的误差要求下，误差控制在e的-6次方内，Gauss迭代用了49次达到要求，而Jacobi迭代用了97次，可见，在迭代中尽量采用最新值，可以大幅度的减少迭代次数，迭代过程收敛快一些。
在Gauss中，初值为100,迭代46次达到精确度l.Oe-6,初值为50时，迭代47次，初值为0时，迭代49次，初值为200时迭代50次，可见存在一个最佳初始值，是迭代最快。
这一点在jacobi迭代中表现的尤为明显。
五、上下边界的热流量：
上边界t=200℃,t=10℃,所以,X
200,
200-100.At200-t,4
热流量中1二、*[
木—+LAX+
△y2Ay
二1*（100/2+（200-139.6088）+（200-150.3312）+（200-153.0517）+（200-153.5639）/2）
=230.2264W
下边界热流量中2二人*[上5&什口1加・+三卫加・+口*4]-
AyAyAyAy2
Ay2AyAyAy△）】2
=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）I=-489.925W=489.25W
六、温度等值线
Gauss:
4一一一」1111
■T—：
r——一
3.5
―—二———―
-
3
--
2.5
-
■
2
—■■,一
15
'二二
1
——
一一
_
1
~——:
；；
1.522.53
3.5
4
Yacobi:
七、计算小结
导热问题进行有限差分数值计算的基本思想是把在时间、空间上连续的温度场用有限个离散点温度的集合来代替，即有限点代替无限点，通过求解根据傅里叶定律和能量守恒两大法则建立关于控制面内这些节点温度值的代数方程，获得各个离散点上的温度值。
要先划分查分网格，在建立差分代数方程组，用MATLAB或者其他软件编程求解。
高斯-赛德尔迭代法和雅克比迭代法区别在于使用新植和旧值进行下一次迭代，而采用新值迭代的高斯-赛德尔迭代收敛的更快些，但其求解代数方程是不一定得到收敛的解，其原因可能由于迭代方式不合适造成。
在计算热流量过程中，主要是正确利用傅里叶定律和牛顿冷却公式，本题中需要特别注意的一点是后边界是绝热的，因而左右方向上几乎不存在热量的传递，所以看似是二维稳态问题实际上是一维稳态的问题。
求解也比较简单。
程序运行出来的等温线结果也很好的说明了这一点，温度总体是从上向下递减，热量传递方向是自上而下。
升级会员 