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

UptoLike

ВУЗ: 

ТПУ | Томск

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

Рубрика: 

98
Ниже приведен листинг программы для решения рассматриваемой
задачи (на языке программирования Pascal)
uses crt;
const mf=500;
eps=1e-5;
type
vector=array[1..mf] of real;
var {раздел описания переменных, которые мы будем использовать в
программе}
i, j, N, s : integer;
T, Ts, Tn, alfa, beta : vector;
ai, bi, ci, fi, max1, max2 : real;
ro, c, h, tau : real;
Th, T0, Tc, L, t_end, time : real;
f, g, f1, f2 : text;
function lamda(x:real):real;
{функция вычисления коэффициента теплопроводности по формуле
(31)}
begin
lamda:=5500/(560+x)+0.942*(1e-10)*x*sqr(x);
end;
begin
clrscr;
{с клавиатуры вводим все необходимые входные параметры}
Writeln('Введите количество пространственных узлов, N');
Readln(N);
Writeln('Введите окончание по времени, t_end');
Readln(t_end);
Writeln('Введите толщину пластины, L');
Readln(L);
Writeln('Введите плотность материала пластины, ro');
Readln(ro);
Writeln('Введите теплоемкость материала пластины, c');
Readln(c);
Writeln('Введите начальную температуру в К, T0');
Readln(T0);
Writeln('Введите температуру в К на границе х=0, Th');
Readln(Th);
Writeln('Введите температуру в К на границе х=L, Tc');
Readln(Tc);

Страницы