do m=2,5 do n=2,11 tt(m,n)=0.25*(t(m-1,n)+t(m+1,n)+t(m,n-1)+t(m,n+1)) !内节点 end do end do do m=6,15 do n=8,11 tt(m,n)=0.25*(t(m-1,n)+t(m+1,n)+t(m,n-1)+t(m,n+1)) end do end do do m=2,5 tt(m,1)=0.25*(2*t(m,2)+t(m-1,1)+t(m+1,1)) !边界1 end do do n=8,11 tt(16,n)=0.25*(2*t(15,n)+t(16,n-1)+t(16,n+1)) end do do n=2,6 tt(6,n)=(2*t(5,n)+t(6,n+1)+t(6,n-1)+2*20*0.1*10/0.53)/(4+2*20*0.1/0.53) end do do m=7,15 tt(m,7)=(2*t(m,8)+t(m+1,7)+t(m-1,7)+2*20*0.1*10/0.53)/(4+2*20*0.1/0.53) end do do n=2,11 tt(1,n)=(2*t(2,n)+t(1,n+1)+t(1,n-1)+2*4*0.1*30/0.53)/(4+2*4*0.1/0.53) !边界3 end do do m=2,15 tt(m,12)=(2*t(m,11)+t(m+1,12)+t(m-1,12)+2*4*0.1*30/0.53)/(4+2*4*0.1/0.53) end do
tt(6,7)=(2*(t(5,7)+t(6,8))+t(7,7)+t(6,6)+2*20*0.1*10/0.53)/(6+2*20*0.1/0.53) !内角点 tt(1,12)=(t(2,12)+t(1,11)+2*4*0.1*30/0.53)/(2+2*4*0.1/0.53) !外角点 tt(1,1)=(t(2,1)+t(1,2)+4*0.1*30/0.53)/(2+4*0.1/0.53) tt(16,12)=(t(16,11)+t(15,12)+4*0.1*30/0.53)/(2+4*0.1/0.53) tt(6,1)=(t(5,1)+t(6,2)+20*0.1*10/0.53)/(2+20*0.1/0.53) tt(16,7)=(t(15,7)+t(16,8)+20*0.1*10/0.53)/(2+20*0.1/0.53) do m=1,16 do n=1,12 wc(m,n)=abs(t(m,n)-tt(m,n)) end do end do wcmax=maxval(wc(1:16,1:12)) do m=1,16 do n=1,12 t(m,n)=tt(m,n) end do
11
end do end do open(44,file=\ do m=1,16 do n=1,12 print*,\ write(44,*)m,n,t(m,n) end do end do tz=0 do n=2,7 tz=tz+t(6,n) end do do m=7,15 tz=tz+t(m,7) end do
qq=4*(0.1*20*(tz-15*10)+0.05*20*(t(6,1)-10)+0.05*20*(t(16,7)-10)) print*,'qq=',qq write(44,*)'qq=',qq end program question02
运行结果：
由于数据量较大所以各点温度见附录或者附件question02.txt
五．每米长度上通过壁面的冷量的计算
第一问：通过壁面的冷量既可以根据内表面计算也可以根据
12
外表面计算，这里根据内表面计算，由于第一问是恒温边界条件，所以只需求出内边界各个节点所获得的热量，并加和即可
Q总??qw*A?161.1212W
第二问：通过壁面的冷量既可以根据内表面计算也可以根据外表面计算，这里根据内表面计算，因为这种情况下是对流换热所以只需计算出所有内表面单个单元的传热量，求和即可
Q总??Qi??1?0.1?h1?(ti?tf1)?124.7063W
ii??等温图：可以根据各点数据，用origin软件画出等温图 问题一等温图：
问题二等温图：
13
附录：1.第一种情况各点温度 m n 温度
1 1 30.00000 1 2 30.00000 1 3 30.00000 1 4 30.00000 1 5 30.00000 1 6 30.00000 1 7 30.00000 1 8 30.00000 1 9 30.00000 1 10 30.00000 1 11 30.00000 1 12 30.00000 2 1 26.06844 2 2 26.08093 2 3 26.12242 2 4 26.20563 2 5 26.35332 2 6 26.59763 2 7 26.97121
14
2 8 27.48136 2 9 28.07841 2 10 28.71271 2 11 29.35634 2 12 30.00000 3 1 22.11191 3 2 22.13286 3 3 22.20311 3 4 22.34679 3 5 3 6 3 7 3 8 3 9 3 10 3 11 3 12 4 1 4 2 4 3 4 4 4 5 4 6 4 7 4 8 4 9 4 10 4 11 4 12 5 1 5 2 5 3 5 4 5 5 5 6 5 7 5 8 5 9 5 10 5 11 5 12 6 1 6 2 6 3
22.61001 23.06598 23.80584 24.87584 26.11955 27.41607 28.71267 30.00000 18.11347 18.13549 18.21038 18.36840 18.67396 19.25044 20.31035 22.09659 24.10791 26.11937 28.07825 30.00000 14.07099 14.08525 14.13453 14.24248 14.46698 14.95149 16.08852 19.09226 22.09612 24.87524 27.48096 30.00000 10.00000 10.00000 10.00000
15