(完整word版)西安交通大学——温度场数值模拟(matlab).docxVIP

温度场模拟matlab代码： clear,clc,clf L1=8;L2=8;N=9;M=9;% 边长为8cm的正方形划分为8*8的格子 T0=500;Tw=100; %初始和稳态温度 a=0.05; %导温系数 tmax=600;dt=0.2; % 时间限 IOmin 和时间步长 0.2s dx=Ll/(M-1) ;dy=L2/(N-l); M l=a*dt∕(dx^2) ;M2=a*dt/(dy 八 2); T=TO*ones(M,N); Tl=TO*ones(M,N); t=O;l=O;k=O; Tc=zeros( 1,600);%中心点温度，每一秒采集一个点 for i=l:9 for j=1:9 if(i==l∣i==9U==lU=9) T(i,j)=Tw;%边界点温度为100℃ else T(ij)=T0; end end end if(2*Ml+2*M2=l) %判断是否满足稳定性条件 while(ttmax+dt) t=t+dt; k=k+l; for i=2:8 for j=2:8 Tl(ij)=Ml*(T(i4j)+T(i+lj))+M2*(T(ij4)+T(i,j+l))+(l-2*Ml-2*M2)*T(i,j); end end for i=2:8 forj=2:8 T(i,jTl(i,j); end end if(k==5) 1=1+1; Tc(l)=T(5,5); k=0; end end i=l:9;j=l:9; [x,y]=meshgrid(i); figure(l); subplot( 1,2,1); mesh(x,y,T(i,j))% 画出 IOmin 后的温度场 axis tight; xlabel(x√FontSize,14)^label(,y7FoπtSize,,14)jzlabel(,T∕°C7FontSize,14) title( 1 min 后二维温度场模拟图VFontSize, 18) subplot(l,2,2); [C,H]=contour(x,y,T(i,j)); clabel(C,H);axis square; xlabel(x,,,FontSize,, 14);ylabel(y,FontSize, 14); title( 1 min 后模拟等温线图？FontSiZe, 18) figure(2); xx=l :600; plot(xx,Tc,,k-,,linewidth,2) xlabel(时 间 ∕s,,FontSize,, 14);ylabel(, 温度 / ℃ ,,FontSize1,14);title(中 心 点的冷却曲线 VFontSize, 18) else disp(E?τor!,) %如果不满足稳定性条件，显示“Error!” end 实验结果: 中心点的冷却曲线 时间/S Imin后二维温度场模拟图 Imin后模拟等温线图 5min后二维温度场模拟图 5min后模拟等温线图 I00.01 - TOO02 IOO03 IOo04 寸OOOL CM。。。L LOOOL a.O12oo?2 100.04 100.01 S s IOmin后二维温度场模拟图 IOmin后模拟等温线图 100.01 U 二 ∣ LOOOL 3,04 喘 Joo03 IOo.02 IOO01 ??β 80。L 寸。0° -? 3.00- 飞 IOO ?03 IoO≡ aoo? 100.02 100.01 Iomin后二维温度场模拟图（不满足稳定性条件） 出心点的冷却曲线（不满足稳定性条件）