二维稳态导热试验报告材料

1、实用标准文案矶家文冠乂学传热学二维导热物体温度场的数值模拟作 者： 陈振兴 学 号： 10037005学院（系）：化工学院专 业：过程装备与控制工程班 级：装备01指导教师：李增耀实验时间：2012-10精彩文档实用标准文案二维导热物体温度场的数值模拟、物理描述有一个用砖砌成的长方形截面的冷空气通道，其截面尺寸和示意图如图1-1所示，假设在垂直纸面方向上冷空气及砖墙的温度变化很小，可以近似地予以忽略。在以下情况下试计算：(1)砖墙横截面上的温度分布；(2)垂直于纸面方向的每米长度上通过砖墙的导热量。1、内外表面均为第三类边界条件，且已知：t：i =30 C,h1 =10.33 W/m2 C2_

2、t：2 =10C,h2 =3.93 W/m C砖墙的导热系数'=0.53W/m C2、内外壁分布均匀地维持在0及30 ；图1-1精彩文档实用标准文案:、数学描述该结构的导热问题可以作为二维问题处理，并且其截面如图1-1所示，由于对称性，仅研究其1/4部分即可。其网络节点划分如图2-1 ;上述问题为二维矩形域内的稳态、无内热源、常物性的导热问题，对于这样的物理问题，我们知道，描写其的微分方程即控制方程，就是导热微分方程：.：2t ；：2t .一 2. 2 = 0x 二 y第三类边界条件：内外表面均为第三类边界条件，且已知：2t：1；-30C,h1 =10.33 W/m C2t,2 -10

3、 C,h2 -3.93 W/m C砖墙的导热系数'=0.53W/m Cne图2-1精彩文档实用标准文案三：方程的离散如上图2-1所示，用一系列与坐标轴平行的网络线把求解区域划分成许多子区域， 以网格 线的交点作为需要确定温度值的空间位置，即节点，节点的位置已该点在两个方向上的标号 m n来表示。每一个节点都可以看成是以它为中心的小区域的代表，如上（ m, n）：对于（m, n）为内节点时：由级数展开法或热平衡法都可以得到，当Ax=Ay时:1 ,、t (t d t i t t - t i)，m,nm：；1,nm -1,nm,n 1 m,n J/4对于（mi, n）为边界节点时:位于平直边

4、界上的节点:1m,n.(tm 1,n 2tm,n t m,n 1) 4外部角点：如图2-1中a、b、d、e、f 点,内部角点：如图2-1中c点,由已知条件有，当t m,n - - ( tm 1 ,n 2tm,n 2tm,n 1 t m,n 1)6m=1或n=13时的节点的温度衡为tw1=30,当（m=6且n<9）和（n=8且6<m<17时的节点的温度为tw2=10 C四：编程思路及流程图开始输入已知参数说明边界条件取定初始试探值精彩文档TA （i,j） =0实用标准文案图3-1五、程序及运行结果第三类边界条件：1、实验程序（C语言）：/ l.cpp :定义控制台应用程序的入口

5、点/#include "stdafx.h"#include<stdio.h>精彩文档实用标准文案#include<math.h> int _tmain(int argc, _TCHAR* argv口)int i,j,l;float dt=1.0,dx=0.1,dy=0.1;float t1317,a1317;float q1=0,q2=0,q=0,e;float lmd=0.53,h1=10.33,h2=3.93,t1=30,t2=10,ep=1.0e-7;/*打印出题目*/printf("ttt二维稳态导热问题 tt");pri

6、ntf("ntttttt-陈振兴装备 n");printf("n 题目：二维导热物体温度场的电模拟实验n");printf("n 矩形区域，l1=2.2;l2=3;l3=2;l4=1.2,假设区域内无内热源，导热系数为常熟，内外表面均为第三类边界条件且已知t1=30;t2=10;h1=10.33;h2=3.93;LMD=0.53; 求该矩形区域内的温度分布及垂直于纸面方向的单位长度上通过墙体的导热量。n");/*各节点上的温度值*/for(j=0;j<17;j+)t0j=30.0;for(i=1;i<13;i+)ti0=3

7、0.0;for(i=7;i<13;i+)ti7=10.0;for(j=8;j<17;j+)t7j=10.0;for(i=1;i<7;i+)for(j=1;j<17;j+) tij=20,aij=0;for(i=7;i<13;i+)for(j=1;j<7;j+)精彩文档实用标准文案tij=20,aij=0;)while(dt>=ep)(for(i=1;i<7;i+)for(j=1;j<17;j+)aij=tij;)(for(i=7;i<13;i+)for(j=1;j<7;j+)aij=tij;)(for(i=6;i<12;i

8、+)for(j=2;j<6;j+)tij=(ti-1j+ti+1j+tij-1+tij+1)/4;)(for(i=2;i<6;i+)for(j=2;j<16;j+)tij=(ti-1j+ti+1j+tij-1+tij+1)/4;)(for(j=2;j<6;j+)t12j=(t12j-1+t12j+1+2*t11j)/4;)(for(i=2;i<6;i+)ti16=(ti-116+ti+116+2*ti15)/4;)(for(i=2;i<12;i+)ti1=(dx*h1*t1/lmd+(ti+11+ti-11)/2+ti2)/(2+dx*h1/lmd);)(f

9、or(j=2;j<16;j+)t1j=(dy*h1*t1/lmd+(t1j+1+t1j-1)/2+t2j)/(2+dy*h1/lmd);)(for(i=7;i<12;i+)精彩文档实用标准文案ti6=(dx*h2*t2/lmd+(ti+16+ti-16)/2+ti5)/(2+dx*h2/lmd);for(j=7;j<16;j+)t6j=(dy*h2*t2/lmd+(t6j+1+t6j-1)/2+t5j)/(2+dy*h2/lmd);t11=(h1*dx*t1+lmd*(t21+t12)/2)/(lmd+h1*dy);t116=(h1*dx*t1+lmd*(t115+t216)

10、/2)/(lmd+h1*dy);t616=(h2*dx*t2+lmd*(t615+t516)/2)/(lmd+h2*dy);t121=(h1*dx*t1+lmd*(t122+t111)/2)/(lmd+h1*dy);t126=(h2*dx*t2+lmd*(t125+t116)/2)/(lmd+h2*dy);t66=(h2*dy*t2+lmd*(t56+t65+t76/2+t67/2)/(3*lmd+h2*dx);for(i=1;i<7;i+)for(j=1;j<17;j+)dt=dt+abs(ti叶aij);dt=dt/(6*16);for(i=7;i<13;i+)for(j

11、=1;j<7;j+)dt=dt+abs(tij-aij);dt=dt/(6*6);printf(" 温度分布为：ttttttttt");l=0;for(i=1;i<7;i+)for(j=1;j<17;j+)printf("%3.1f ",tij);l=l+1;if(l=16)printf("n");l=0;精彩文档实用标准文案l=0; for(i=7;i<13;i+)for(j=1;j<7;j+)printf("%3.1f "阳j);l=l+1;if(l=6) printf("

12、;n"); l=0; for(j=2;j<16;j+) q1=q1+(30-t1j)*h1*dx;for(i=2;i<12;i+)q1=q1+(30-ti1)*h1*dy;for(j=7;j<17;j+)q2=q2+(t6叶10)*h2*dx;for(i=7;i<12;i+)q2=q2+(ti6-10)*h2*dy;q1=q1+h1*(dx/2*(30-t116)+dy/2*(30-t121)+dx*(30-t11);q2=q2+h2*(dx/2*(t616-10)+dy/2*(t126-10)+dx*(t77-10);q=(q1+q2)/2;e=abs(q2

13、-q1)/q);printf(" 单位长度上1/4墙体的导热量为：4.2fW偏差为：3.2f",q,e); getchar();getchar();return 0;运行结果图：精彩文档实用标准文案c:UsersAd mi nistratorD«skto p - 1D ebu gl. exf二维稳态导热题目：二维导热物体温度场的电模拟实能域二亲?.区,分2,11 =2 2 12=3； 1?Z12假设区域内无内热源，导热系鹦为常熟，内夕 类边界条件且已知11 =3小七2 =1匹hl=10.33;h2=3.93;LMD=0.53;求该矩形区X 幸直于纸面方向的单位长度

14、上通过墙体的导热量,为：7 29.5 29,3 29.1 29,0 28,R 2R.7 28,? 28,6 28.6 28.5 28.S 2心5 28,t22.? 22.? 22.7 22*?21.9 26.Q 18.717.517,2 17,1 17.014.614.314-2 14.1 14.013.?1312-917.916.929.1 27.424.828.7 25.9 23.0 20.0 16.9 13.823.1 20,Z 17,2 14.128.7 26.2 23.6 20,9 17,9 14.629.1 25.9 23.0 19.9 16,7 13.Q三位长度上工/4墙体三一导

15、热量为：26.73皿偏差为：0.能16.6 16,4图3-2实验算得导热量为97.62皿 与数值模才H的偏差为(26.73*4W-97.62W)/ (26.73*4 ) W*100%=8.7%2、数值模拟程序(matlab ):z=29.9 29.7 29.5 29.3 29.1 29.0 28.8 28.7 28.7 28.6 28.6 28.5 28.5 28.5 28.5 28.229.7 29.1 28.5 28.027.426.926.5 26.2 26.025.9 25.8 25.7 25.6 25.625.625.529.5 28.5 27.6 26.625.724.824.1

16、23.6 23.323.1 22.9 22.8 22.8 22.722.722.729.3 28.0 26.6 25.223.822.521.6 20.9 20.520.2 20.0 19.9 19.9 19.819.819.829.1 27.4 25.7 23.821.919.918.7 17.9 17.517.2 17.1 17.0 16.9 16.916.916.929.0 26.9 24.8 22.519.916.715.2 14.6 14.314.2 14.1 14.0 14.0 14.014.013.928.8 26.5 24.1 21.518.615.20 0 0 0 0 0 0

17、 00 0;28.7 26.2 23.5 20.817.814.60 0 0 0 0 0 0 00 0;28.6 25.9 23.2 20.317.314.20 0 0 0 0 0 0 00 0;28.5 25.7 22.8 20.017.014.00 0 0 0 0 0 0 00 0;28.3 25.4 22.6 19.716.713.60 0 0 0 0 0 0 00 0;精彩文档实用标准文案27.2 25.1 22.5 19.5 16.3 12.5 0 0 0 0 0 0 0 0 0 0;v=18 22 26;xx,yy=meshgrid(y,x);surf(xx,yy,z);color

18、bar;xlabel( 'x' );ylabel( 'y' );zlabel( 'z' );az = 0;el =- 90;view(az, el);shading interp;axis tight;figure,contour(xx,yy,z,v);grid on数值模拟图：Figure No. 2File Edit Tock Window Help75 口昌 it A Z Z母月C精彩文档123456789012-11 11LLi i 1X1r*1" 1ii1 1 1111 /| 1r1 111 i一4-11IIIii 1 - H

19、i -h -*II1-T- 1 1 I 1i 1 1 i 曲I 1 1 i 1 i|i 1 1 1 1 III 1 1 1 I 1 | I 1 1- 1i 1 1246J.10121416图3-3实用标准文案图3-4第一类边界条件1、实验程序(C语言)：#include "stdafx.h"#include<stdio.h>#include<math.h>int _tmain(int argc, _TCHAR* argv口)精彩文档实用标准文案int i,j,l；float dt=1.0,dx=0.1,dy=0.1;float t1216,a1216

20、;float q1=0,q2=0,q=0,e;float lmd=0.53,t1=30,t2=0,ep=1.0e-7;/* 打印出题目*/printf("ttt二维稳态导热问题 tt");printf("ntttttt-陈振兴装备 n");printf("n 题目：二维导热物体温度场的电模拟实验n");printf("n矩形区域，l1=2.2;l2=3;l3=2;l4=1.2,假设区域内无内热源，导热系数为常熟，内外表面均为第一类边界条件且已知t1=30;t2=0;LMD=0.53;求该矩形区域内的温度分 布及垂直于纸面方向

21、的单位长度上通过墙体的导热量。n");/*各节点上的温度值*/for(j=0;j<16;j+)t0j=30.0;for(i=1;i<12;i+)ti0=30.0;for(i=5;i<12;i+)ti5=0.0;for(j=6;j<16;j+)t5j=0.0;for(i=1;i<5;i+)for(j=1;j<15;j+) tij=20,aij=0;for(i=5;i<11;i+)for(j=1;j<5;j+) tij=20,aij=0;while(dt>=ep)精彩文档实用标准文案for(i=1;i<5;i+)for(j=1;

22、j<15;j+) aij=tij； )(for(i=5;i<11;i+)for(j=1;j<5;j+) aij=tij;)(for(i=5;i<11;i+)for(j=1;j<5;j+)tij=(ti-1j+ti+1j+tij-1+tij+1)/4;)(for(i=1;i<5;i+)for(j=1;j<15;j+)tij=(ti-1j+ti+1j+tij-1+tij+1)/4;)(for(j=1;j<5;j+)t11j=(t11j-1+t11j+1+2*t10j)/4;)(for(i=1;i<5;i+)ti15=(ti-115+ti+115

23、+2*ti14)/4;)(for(i=1;i<5;i+)(for(j=1;j<15;j+)dt=dt+abs(tij-aij);dt=dt/(4*14);)for(i=5;i<11;i+)(for(j=1;j<5;j+)dt=dt+abs(tij-aij);dt=dt/(4*6);)精彩文档实用标准文案)printf("温度分布为：ttttttttt");1=0;(for(i=0;i<6;i+)(for(j=0;j<16;j+)printf("%4.1f ",tij);1=1+1;if(l=16)(printf(&qu

24、ot;n");l=0;)l=0;for(i=6;i<12;i+)(for(j=0;j<6;j+) printf("%4.1f ",tij);l=l+1;if(l=6)printf("n");l=0;)for(j=1;j<15;j+)q1=q1+(30-t1j)*lmd;for(i=1;i<11;i+)q1=q1+(30-ti1)*lmd;for(j=6;j<15;j+)q2=q2+t4j*lmd;for(i=6;i<11;i+)q2=q2+ti4*lmd;q1=q1+lmd*(30-t115)/2+(30-t

25、111)/2);q2=q2+lmd*(t415/2+t114/2+t44);精彩文档实用标准文案q=(q1+q2)/2;e=abs(q2-q1)/q);printf("单位长度上墙体的导热量为：4.2fW偏差为：3.2f",4*q,e);getchar();getchar();return 0;运行结果图：回目：二维导热物体温度场的电模拟实验懒黑n.区第的分3,112.2；12斗门32门4nl.2,假设区域哑 类边界条件且已知11 =3小七2 =0 ; LMD=0.5力求百 住长度上通过墙体的导热量“ 为:0 30,0 30,0 30.0 30.Q 30.0 30.0 38

26、,0 30一一陈振兴装备29.0 28,1 Q.0 27.1 0.0 262 Q.0 25.S 0.0 24.9 Q.0 24.5 0.Q 24.3 8.Q 24.22«.1 27.126.1 24.224.2 21.222.3 18.120.7 15-S19,& 13.918.9 L3.B18.5 12.618.3 12-326.2 25.522.3 20.?18.1 15.513.6 9.19,1 0.B7.4 0.0 6-7 6.06-4 0.0 fi.224.9 24.519.6 18.913.9 13.a?.4 6.7B.0 0.024.3 24-21B.5 18.

27、312.S 12-36.4 6*2 0.024.1 24.18.2 LB.1 12 .2 12.16.1 E,1 0.0 0.024.0 24.018.a 18.B12.0 12.06.0 6.00.0 0.024.0 24.01B.B 18.B12.B 12.B 6*0 6.00.Bc:UsertAcl mi nistratorDeskto p" 1D e bu gl. ex?B.U 24.1 1«.2 12-2 60.Q三位长度上。4墙体引导热量为：5"20凡偏差为：0.042、数值模拟程序(matlab ):z=30.0 30.0 30.0 30.0 30.

28、0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0 30.0;30.0 29.0 28.1 27.126.2 25.524.9 24.524.3 24.224.1 24.124.0 24.024.024.0;30.0 28.1 26.1 24.222.3 20.719.6 18.918.5 18.318.2 18.118.0 18.018.018.0;30.0 27.1 24.2 21.218.1 15.513.9 13.012.5 12.312.2 12.112.0 12.012.012.0;30.0 26.2 22.3 18.113.6

29、9.1 7.4 6.7 6.46.2 6.1 6.1 6.0 6.0 6.0 6.0;30.0 25.5 20.715.59.10.0 0.0 0.0 0.00.00.00.00.0 0.0 0.0 0.0;30.0 24.9 19.613.97.40.0 0.0 0.0 0.00.00.00.00.0 0.0 0.0 0.0;30.0 24.5 18.913.06.70.0 0.0 0.0 0.00.00.00.00.0 0.0 0.0 0.0;精彩文档实用标准文案30.0 24.318.512.66.40.0 0.00.0 0.00.0 0.00.00.00.00.0 0.0;30.0 24.218.312.36.20.0 0.00.0 0.00.