(完整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(