




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、东南大学能源与环境学院课程作业报告作业名称:传热学大作业一一利用 matlab程序解决热传导问题院系:能源与环境学院专业:建筑环境与设备工程学号:姓名:2014年11月9日12一、题目及要求1. 原始题目及要求2. 各节点的离散化的代数方程3. 源程序4. 不同初值时的收敛快慢5. 上下边界的热流量(入=1W/(m C)6. 计算结果的等温线图7. 计算小结题目:已知条件如下图所示:100'C绝然Tf=lg X=lW/(m X?)二、各节点的离散化的代数方程各温度节点的代数方程ta=(300+b+e)/4 ; tb=(200+a+c+f)/4; tc=(200+b+d+g)/4; td
2、=(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)/24; tp=(l+o+100)/12三、源程序【G-S迭代程序】【方法一】函数文件为:function y,n=gauseidel(A,b,x0,eps)D=diag(diag(A);L=-tri
3、l(A,-1);U=-triu(A,1);G=(D-L)U;f=(D-L)b;y=G*x0+f;n=1;while norm(y-x0)>=epsx0=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,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
4、,-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,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;
5、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.0e-6) xx=1:1:4;yy=xx;X,Y=meshgrid(xx,yy);Z=reshape(x,4,4);Z=Z'contour(X,Y,Z,30)Z =139.6088 150.3312 153.0517 1
6、53.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(1,3)=100;t(1,4)=100;t(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
7、,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)/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,
8、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,5)+200)/24;t(5,5)=(t(5,4)+t(4,5)+100)/12;t'endcontour(t',50);ans =100.0000 200.0000 200.0000 200.0000 200.0000100.0000 136.8905 146.967
9、4 149.8587 150.7444100.0000 102.3012 103.2880 103.8632 104.3496100.0000 70.6264 61.9465 59.8018 59.6008100.0000 19.0033 14.8903 14.5393 14.5117Jacobi迭代程序】函数文件为: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)>=epsx0=y;y=B*x0+f;n=
10、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,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
11、,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,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
12、,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,yy);Z=reshape(x,4,4);Z=Z'contour(X,Y,Z,30)n =97Z =139.6088 150.3312 153.0517 153.5639108.1040 108.6641 108.3119 108.152384.1429 67.9096 63.3793 62.421420.1557 15.4521 1
13、4.8744 14.7746四、不同初值时的收敛快慢1、 方法1在Gauss迭代和Jacobi迭代中,本程序应用的收敛条件均为norm(y-x0)>=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=1:10判断收敛,为迭代10次,若改为1:20,则迭代20次。2、在同样的误差要求下,误差控制在 e的-
14、6次方内,Gauss迭代用了 49次达到要求,而Jacobi 迭代用了 97次,可见,在迭代中尽量采用最新值,可以大幅度的减少迭代次数,迭代过程收敛快 一些。在Gauss中,初值为100,迭代46次达到精确度1.0e-6,初值为50时,迭代47次,初值为0时,迭代49次,初值为200时迭代50次,可见存在一个最佳初始值,是迭代最快。这一点在jacobi迭代中表现的尤为明显。五、上下边界的热流量:上边界t=200 C, t =10 C,所以,热流量 1=入* 22*M220dAAx + 2Ax+22Ax+A* 竺2.:y:y:y.:y2=1* ( 100/2+(200- 139.6088)+(2
15、00-150.3312)+(200-153.0517)+(200-153.5639)/2) =230.2264W下边界热流量2=|入*100-10*7心m . :x+ InE .x+里上 .侦 + 址E * 快- yyyy 2k* .:x+k* .次 + 1* *+*-)|yyyy 2=|1*( 84.1429-20.1557)+(67.9096-15.4521)+(63.3793-14.8744)+(62.4214-10)+(14.7746-10)/214.7746 )/2)-10*(90/2+(20.1557-10)+(15.4521-10)+(14.8744)| = |-489.925|W =489.25W六、温度等值线Gauss:Yacobi:七、计算小结导热问题进行有限差分数值计算的基本思想是把在时间、空间上连续的温度场用有限个离散点温度的集合来代替,即有限点代替无限点,通过求解根据傅里叶定律和能量守恒两大法则建立关于控制 面内这些节点温度值的代数方程,获得各个离散点上的温度值。要先划分查分网格,在建立差分代数方程组,用MATLAB者其他软件编程求
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 农村合伙购房合同范例
- 4000美元供气合同范例
- 假期转租合同范例
- 兑诊所合同范例
- 公路钢筋安装施工合同范本
- 两人合租协议合同范例
- 个人房屋赠与合同范例
- 供热改造维修合同范例
- 个人签订物业服务合同范例
- 内控外包合同范例
- 人教版四年级上册语文《一单元》测试卷【及答案】
- 复数算符在图像处理中的应用
- 百融云创风险决策引擎V5产品操作手册
- GB 15979-2024一次性使用卫生用品卫生要求
- 2024年合肥市轨道交通集团有限公司招聘笔试冲刺题(带答案解析)
- CJJT8-2011 城市测量规范
- 故事绘本后羿射日
- 产前筛查标准技术操作规程
- ISO27001:2022信息安全管理手册+全套程序文件+表单
- 2024年广州市高三一模高考物理试卷试题答案(精校打印)
- 国测省测四年级劳动质量检测试卷
评论
0/150
提交评论