%%詳細推導求解過程請移步?微信公眾號?CFD入門指南?%%
%%二維穩態導熱問題第一類邊界%%
clc
clear

L1=0.6;?????%板長
L2=0.4;?????%板寬
T1=100;?????%邊界溫度
T2=60;
lan=200;????%導熱系數
m=61;???????%長度的節點數量
n=81;???????%寬度的節點數量
dx=L1/（m-1）;????%長度的網格尺寸
dy=L2/（n-1）;????%寬度的網格尺寸

v1=zeros（mn）;

for?i=1:m????????%邊界條件
????v1（i1）=T2;
????v1（in）=T1;
end
for?j=1:n
????v1（1j）=T2;
????v1（mj）=T2;????
end

dt=1;
v2=v1;
k=0;

while?dt>1e-3????????%兩次迭代誤差小于10^-3時計算結束
????k=k+1;
????dt=0;
????for?i=2:m-1??????%兩端的溫度值在初始化時已定義，所以不需再計算
????????for?j=2:n-1
????????????v2（ij）=（v1（ij+1）+v1（ij-1））*（dx^2/2/（dx^2+dy^2））+（v1（i+1j）+v1（i-1j））*（dy^2/2/（dx^2+dy^2））;
????????????dt=max（（v2（ij）-v1（ij））dt）;
????????end
????end
????v1=v2;
end

%%詳細推導求解過程請移步?微信公眾號?CFD入門指南?%%

figure
surf（v2）;
a=xlabel（‘長度‘）;
b=ylabel（‘寬度‘）;
c=zlabel（‘溫度‘）;
set（a‘FontSize‘15‘）;
set（b‘FontSize‘15‘）;
set（c‘FontSize‘15‘）;
k
%%詳細推導求解過程請移步?微信公眾號?CFD入門指南?%%

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1094??2020-03-11?16:18??Steady_2D\Steady_2D_B1.m

?????文件???????1130??2020-03-11?16:18??Steady_2D\Steady_2D_B2.m

?????文件???????1444??2020-03-11?16:18??Steady_2D\Steady_2D_B3.m

?????目錄??????????0??2020-03-11?16:26??Steady_2D

-----------?---------??----------?-----??----

?????????????????3668????????????????????4