Разностные методы решения задач теплопроводности. Кузнецов Г.В - 54 стр.

UptoLike

ВУЗ: 

ТПУ | Томск

Составители: 

Рубрика: 

54
Readln(c1);
Writeln('Введите теплоемкость материала второй части пластины, c2');
Readln(c2);
Writeln('Введите температуру на границе х = 0, Tl');
Readln(Tl);
Writeln('Введите температуру на границе х = L, Tr');
Readln(Tr);
Writeln('Введите начальную температуру, T0');
Readln(T0);
{определяем общее число узлов в пластине}
N:=N1+N2+1;
{определяем расчетный шаг сетки по пространственной координате}
h:=L/(N-1);
{определяем коэффициенты температуропроводности}
a1:=lamda1/(ro1*c1);
a2:=lamda2/(ro2*c2);
{определяем расчетный шаг сетки по времени}
tau:=t_end/100.0;
{определяем поле температуры в начальный момент вр
емени}
for i:= 1 to N do
T[i]:=T0;
{проводим интегрирование нестационарного уравнения
теплопроводности}
time:=0;
while time<t_end do {используем цикл с предусловием}
begin
{увеличиваем переменную времени на шаг τ}
time:=time+tau;
{определяем начальные прогоночные коэффициенты на основе левого
граничного условия}
alfa[1]:=0.0;
beta[1]:=Tl;
{цикл с параметром для определения прогоночных коэффициентов по
формуле (8) в первой части пластины}
for i:= 2 to N1 do
b
egin
{ai, bi, ci, fi – коэффициенты канонического представления СЛАУ с
трехдиагональной матрицей}
ai:=lamda1/sqr(h);
bi:=2.0*lamda1/sqr(h)+ro1*c1/tau;
ci:=lamda1/sqr(h);

Страницы