
- •Отчет по лабораторной работе №1
- •Постановка задачи:
- •Аналитическое решение
- •Метод конечных разностей (явная схема)
- •Метод конечных разностей (неявная схема)
- •Регулярный режим охлаждения
- •Результаты расчетов Сравнение численных решений с аналитическим
- •Зависимости температуры по сечению
- •Темп охлаждения
- •1 Метод.
- •2 Метод.
- •3 Метод (графический).
- •Кривые охлаждения
- •Аналитический метод
- •Мкр явная схема
- •Мкр неявная схема
Мкр явная схема
program yavnaya;
const
alpha=22; lamda=0.16; r0=0.021; p=1200; c=1380; n=21; tau=3200;
a=lamda/(p*c);hi=r0/n; hj=hi*hi/(2*a); m=round(tau/hj+1); t0=95; tcp=15;
var
i,j:integer; r:real;
t:array[0..n,0..m] of real;
begin
for i:=0 to n do
t[i,0]:=t0;
for j:=0 to m-1 do
begin
for i:=1 to n-1 do
begin
r:=i*hi;
t[i,j+1]:=a*hj/hi/hi*(t[i+1,j]-2*t[i,j]+t[i-1,j])+a*hj/2/hi/r*(t[i+1,j]-t[i-1,j])+t[i,j];
end;
t[0,j+1]:=t[1,j+1];
t[n,j+1]:=(tcp*alpha*hi/lamda+t[n-1,j+1])/(1+alpha*hi/lamda);
end;
for i:=0 to n do
begin
write (' | ',t[i,0]:2:3,' |');
write (' | ',t[i,round(m/100)]:2:3,' |');
write (' | ',t[i,round(m/2)]:2:3,' |');
writeln (' | ',t[i,m]:2:3,' |');
end;
end.
Мкр неявная схема
program neyavnay;
const
r0=0.021; t0=95; tcp=15; lambda=0.16; ro=1380; c=1200; alfa=22; a=lambda/ro/c;
tau=3200; n=21; m=619; hi=r0/n; hj=tau/m;
var
i,j,tg:integer;
r,X,Y,Z:real;
t:array [0..n,0..m] of Real;
alpha,betta:array[1..30] of real;
begin
for i:=0 to n do
t[i,0]:=t0;
for j:=0 to m-1 do
begin
for i:=1 to n-1 do
begin
r:=i*hi;
X:=a*hj/hi*(1/hi-1/2/r);
Y:=1+2*a*hj/sqr(hi);
Z:=a*hj/hi*(1/hi+1/2/r);
alpha[1]:=1;
betta[1]:=0;
alpha[i+1]:=Z/(Y-X*alpha[i]);
betta[i+1]:=(X*betta[i]+t[i,j])/(Y-X*alpha[i]);
end;
t[n,j+1]:=(tcp*alfa*hi/lambda+betta[n])/((1+alfa*hi/lambda)*(1-alpha[n]/(1+alfa*hi/lambda)));
for i:=n-1 downto 0 do
begin
t[i,j+1]:=alpha[i+1]*t[i+1,j+1]+betta[i+1];
t[0,j+1]:=t[1,j+1];
end;
end;
for i:=0 to n do
begin
write (' | ',t[i,0]:2:3,' |');
write (' | ',t[i,round(m/100)]:2:3,' |');
write (' | ',t[i,round(m/2)]:2:3,' |');
writeln(' | ',t[i,m]:2:3,' |');
end;
end.