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

UptoLike

ВУЗ: 

ТПУ | Томск

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

Рубрика: 

70
Readln(Th);
Writeln('Введите температуру на границе х = L области решения, Tc');
Readln(Tc);
Writeln('Введите начальную температуру, T0');
Readln(T0);
{определяем расчетные шаги сетки по пространственным
координатам}
hx:=L/(Nx-1);
hy:=H/(Ny-1);
{определяем коэффициент температуропроводности}
a:=lamda/(ro*c);
{определяем расчетный шаг сетки по времени}
tau:=t_end/100.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;
{решаем СЛАУ в направлении оси Ох для определения поля
температуры на промежуточном временном слое}
for j:=1 to Ny do
begin
{определяем начальные прогоночные коэффициенты на основе левого
граничного условия}
alfa[1]:=0.0;
beta[1]:=Th;
{цикл с параметром для определения прогоночных коэффициентов по
формуле (8)}
for i:= 2 to Nx-1 do
begin
{ai, bi, ci, fi – коэффициенты канонического представления СЛАУ с
трехдиагональной матрицей}
ai:=lamda/sqr(hx);
bi:=2.0*lamda/sqr(hx)+ro*c/tau;
ci:=lamda/sqr(hx);

Страницы