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

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);M1=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:9for j=1:9if(i=1|i=9

2、|j=1|j=9)T(i,j)=Tw;% 邊界點溫度為100c elseT(i,j)=T0;endendendif(2*M1+2*M2<=1) %判斷是否滿足穩(wěn)定性條件while(t<tmax+dt)t=t+dt;k=k+1;for i=2:8for j=2:8T1(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); end end for 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=contour(x,y,T(i,j);cl

4、abel(C,H);axis square;xlabel('x','FontSize',14);ylabel('y','FontSize',14);title('1min 后模擬等溫線圖','FontSize',18)figure(2);xx=1:600;plot(xx,Tc,'k-','linewidth',2)xlabel('時 間 /s','FontSize',14);ylabel('溫度 / C ','F

5、ontSize',14);title('中心 點的冷卻曲線 ','FontSize',18)else disp('Error!) % 如果不滿足穩(wěn)定性條件，顯示“ Error!"end實驗結(jié)果:中心點的冷卻曲線1min后二維溫度場模擬圖100.011 30)91 01 0101 0 0100. 04100.0 21 00.01 .1min后模擬等溫線圖100.01100.06100.05100.04100.02100.0110095min后二維溫度場模擬圖5min后模擬等溫線圖105104102101100910min后模擬等溫線圖T 100.03T 1031021 01 ：-1QUO1 o O- 1 02 100-3100£2n OO1 InuDUI 11100干4MOO1006bloop -I I03o ()O1O5u.np40 QU1he041O.OU12U.DU1 1000305 100.03