ВУЗ:
Составители:
146
Writeln('Введите теплоемкость материала пластины, c');
Readln(c);
Writeln('Введите коэффициент теплообмена на границах х,y = 0, kapa1');
Readln(kapa1);
Writeln('Введите коэффициент теплообмена на границах х = L, y = H,
kapa2');
Readln(kapa2);
Writeln('Введите температуру внешней среды относительно границ х,y
= 0, Te1');
Readln(Te1);
Writeln('Введите температуру внешней среды относительно границ х =
L, y = H, Te2');
Readln(Te2);
Writeln('Введите приведенную степень черноты, eps1');
Readln(eps1);
Writeln('Введите начальную температуру, T0');
Readln(T0);
{определяем расчетные шаги сетки по пространственным
координатам}
hx:=L/(Nx-1);
hy:=H/(Ny-1);
{определяем расчетный шаг сетки по времени}
tau:=t_end/1000.0;
{определяем поле температуры в начальный момент вр
емени}
for i:= 1 to Nx do
for j:= 1 to Ny do
T[i,j]:=T0;
{проводим интегрирование нестационарного уравнения
теплопроводности}
time:=0;
while time<t_end do {используем цикл с предусловием}
begin
{увеличиваем переменную времени на шаг τ}
time:=time+tau;
{запоминаем поле температуры на n-ом временном слое}
for i:=1 to Nx do
for j:=1 to Ny do
Tn[i,j]:=T[i,j];
{решаем СЛАУ в направлении оси Ох для определения поля
темпера
туры на промежуточном (n+1/2) временном слое}
for j:=1 to Ny do
Страницы
- « первая
- ‹ предыдущая
- …
- 144
- 145
- 146
- 147
- 148
- …
- следующая ›
- последняя »