ВУЗ:
Составители:
137
Writeln('Введите начальную температуру, T0');
Readln(T0);
Writeln('Введите температуру на границе х = 0, Tc');
Readln(Tc);
Writeln('Введите температуру фазового перехода, Tfr');
Readln(Tfr);
Writeln('Введите теплоту фазового перехода, Qfr');
Readln(Qfr);
Writeln('Введите влагосодержание грунта, w');
Readln(w);
{определяем расчетный шаг сетки по пространственной координате}
h:=L/(N-1);
{определяем коэффициенты температуропроводности}
a1:=lamda1/(ro1*c1);
a2:=lamda2/(ro2*c2);
{определяем поле температуры в начальный момент времени}
time:=0;
for i:= 1 to N do
T[i]:=T0;
{определяем файл, содержащий шаги по времени}
Assign(f1,'ti
me.txt');
Rewrite(f1);
{определяем положение границы фазового перехода}
k:=1;
{проводим интегрирование нестационарного уравнения
теплопроводности}
while time<t_end do {используем цикл с предусловием}
begin
{запоминаем поле температуры на предыдущем временном слое}
for i:= 1 to N do
Tn[i]:=T[i];
{граница фазового перехода на каждом соответствующем временном
шаге смещается на пространственный шаг вправо}
inc(k);
{цикл с постусловием, позволяющий итерационно вычислять поле
температуры, всл
едствие наличия нелинейности в граничном условии
(61)}
repeat
{запоминаем поле температуры на предыдущей итерации}
for i:= 1 to N do
Ts[i]:=T[i];
Страницы
- « первая
- ‹ предыдущая
- …
- 135
- 136
- 137
- 138
- 139
- …
- следующая ›
- последняя »