我的for循环里有很多代码,而且有用,现帖一部分如下:
//中间节点的差分式
for i:=1 to (m-1)do
for j:=0 to (n-1)do
begin
if t0[i,j]>=Tl then
begin
lam:=lam_l;
cp:=cp_l;
rou:=rou_l;
end
else
if t0[i,j]<=Ts then
begin
lam:=lam_s;
cp:=cp_s;
rou:=rou_s;
end
else
begin
lam_sl:=lam_s+((lam_l-lam_s)/(Tl-Ts))*(t0[i,j]-Ts);
lam:=lam_sl;
cp_sl:=cp_s+lf/(Tl-Ts);
cp:=cp_sl;
rou:=rou_sl;
//判断是否满足初始条件
test_sl:=(lam_sl*dt)/(rou_sl*cp_sl*dx*dx);
if test_sl>(1/(4+(4*dx*h)/lam_sl)) then
begin
judge:=false;
exit;
end;
end;
t1[i,j]:=t0[i,j]+(lam*dt/(rou*cp*dx*dx))*(t0[i+1,j]+t0[i-1,j]+t0[i,j+1]
+t0[i,j-1]-4*t0[i,j]);
end;