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

UptoLike

ВУЗ: 

ТПУ | Томск

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

Рубрика: 

21
Ниже приведен листинг программы (явная схема) для решения
рассматриваемой задачи (на языке программирования Pascal)
uses crt;
const mf=500;
type
vector=array[0..mf] of real;
var {раздел описания переменных, которые мы будем использовать в
программе}
i, j, N : integer;
T, TT : vector;
a, lamda, ro, c, h, tau : real;
Tl, T0, Tr, L, t_end, time : real;
f, g : text;
begin
clrscr;
{с клавиатуры вводим все необходимые входные параметры}
Writeln('Введите количество узлов по пространственной координате,
N');
Readln(N);
Writeln('Введите окончание по времен
и, t_end');
Readln(t_end);
Writeln('Введите толщину пластины, L');
Readln(L);
Writeln('Введите коэффициент теплопроводности материала пластины,
lamda');
Readln(lamda);
Writeln('Введите плотность материала пластины, ro');
Readln(ro);
Writeln('Введите теплоемкость материала пластины, c');
Readln(c);
Writeln('Введите начальную температуру, T0');
Readln(T0);
Writeln('Введите температуру на границе х=0, Tl');
Readln(Tl);
Writeln('Введите температуру на границе х=L, Tr');
Readln(Tr);
{определяем коэффициент температуропроводности}
a:=lamda/(ro*c);
{определяем расчетный шаг сетки по пространственной координате}
h:=L/(N-1);

Страницы