西安交通大學(xué)——溫度場數(shù)值模擬matlab.doc

1、溫度場模擬 matlab 代碼：clear,clc,clfL1=8;L2=8;N=9;M=9;% 邊長為 8cm 的正方形劃分為 8*8 的格子 T0=500;Tw=100; % 初始和穩(wěn)態(tài)溫度a=0.05; % 導(dǎo)溫系數(shù) tmax=600;dt=0.2; % 時間限 10min 和時間步長 0.2s dx=L1/(M-1);dy=L2/(N-1);M仁 a*dt/(dxA2);M2=a*dt/(dyA2);T=T0*ones(M,N); T1=T0*ones(M,N);t=0;l=0;k=0;Tc=zeros(1,600);% 中心點溫度，每一秒采集一個點for i=1:9 for j=1:

2、9if(i=1|i=9|j=1|j=9)T(i,j)=Tw;% 邊界點溫度為100CelseT(i,j)=T0;endendendif(2*M1+2*M2=1) % 判斷是否滿足穩(wěn)定性條件while(ttmax+dt)t=t+dt;k=k+1;for i=2:8for j=2:8 T1(i,j)=M1*(T(i-1,j)+T(i+1,j)+M2*(T(i,j-1)+T(i,j+1)+(1-2*M1-2*M2)*T(i,j);endendfor i=2:8for j=2:8T(i,j)=T1(i,j);endendif(k=5)l=l+1;Tc(l)=T(5,5);k=0;endendi=1:9

3、;j=1:9;x,y=meshgrid(i);figure(1);subplot(1,2,1);mesh(x,y,T(i,j)% 畫出 10min 后的溫度場axis tight;xlabel(x,FontSize,14);ylabel(y,FontSize,14);zlabel(T/ C ,FontSize,14)title(1min 后二維溫度場模擬圖,FontSize,18)subplot(1,2,2);C,H=co ntour(x,y,T(i,j);clabel(C,H);axis square;xlabel(x,Fo ntSize,14);ylabel(y,Fo ntSize,14)

4、;title(1min 后模擬等溫線圖,FontSize,18)figure(2);xx=1:600;plot(xx,Tc,k-,li newidth,2)xlabel(時 間 /s,FontSize,14);ylabel(溫度 / C ,FontSize,14);title(中心 點的冷卻曲線,Fo ntSize,18)else disp(Error!) % 如果不滿足穩(wěn)定性條件，顯示 Error!”end實驗結(jié)果:中心點的冷卻曲線106105100.06100.05100.04100.03100.02100.0110095min后二維溫度場模擬圖10410310min后二維溫度場模擬圖1min后模擬等溫線圖II100.01100.06100.05100.04100.03100.02100.01100 .951100.02876y 54321123456789X4unu1 -mu .so 9 1OU-2UODI10001 1QU100-061 001少10 0.040510 0.0 15min后模擬等溫線圖910103050410210110UO1，05876y 54321123456789X10min后模擬等溫線圖9876y 54321100.0 1(6 10 0.100.0 21 00 .01 -10,%1 00. 03100.0403 100.100.4 031000