西安交通大學材料制備與成型實驗——溫度場數值模擬,matlab編程

溫度場模擬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; % 時間限10min和時間步長0.2s dx=L1/(M-1);dy=L2/(N-1); M1=a*dt/(dx^2);M2=a*dt/(dy^2);

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:9

if(i==1|i==9|j==1|j==9)

T(i,j)=Tw;% 邊界點溫度為100℃ else

T(i,j)=T0;

end end end

if(2*M1+2*M2<=1) % 判斷是否滿足穩定性條件 while(t

for 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); end end

for i=2:8

for j=2:8 T(i,j)=T1(i,j); end end if(k==5) l=l+1; Tc(l)=T(5,5); k=0; end

總結

以上是生活随笔為你收集整理的温度 数值模拟 matlab,西安交通大学——温度场数值模拟(matlab)的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。