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

UptoLike

ВУЗ: 

ТПУ | Томск

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

Рубрика: 

91
где
,
2
,
2
1
21
1
21
n
i
n
i
n
i
n
i
n
i
n
i
λ+λ
=λ
λ+λ
=λ
−
−
+
+
при этом
n
i
λ вычисляются по
формуле (37), например,
(
)
3
10
10942.0
560
5500
n
i
n
i
n
i
T
T
−
⋅+
+
=λ . Добавляя к
системе (40) конечно-разностные аналоги краевых условий:
.0 ,
;0 ,
;1,2 ,
1
0
0
>=
>=
−==
nTT
nTT
NiTT
c
n
N
h
n
i
получим замкнутую разностную задачу.
Полученную систему можно свести к наиболее общему виду:
i
n
ii
n
ii
n
ii
FTCTBTA =⋅+⋅−⋅
+
−
++
+
1
1
11
1
,
где
n
ii
n
i
n
i
i
n
i
ii
T
c
F
c
h
B
h
CA
τ
ρ
−=
τ
ρ
+
λ+λ
=
λ
==
−++
, ,
2
2121
2
21
,
,
2
,
2
1
21
1
21
n
i
n
i
n
i
n
i
n
i
n
i
λ+λ
=λ
λ+λ
=λ
−
−
+
+
()
3
10
10942.0
560
5500
n
i
n
i
n
i
T
T
−
⋅+
+
=λ .
Прогоночные коэффициенты находятся по формулам (8). Далее
неизвестное поле температуры определяется по выражению (7).
Блок-схема аналогична представленной в пункте 2.1.
Ниже приведен листинг программы для решения рассматриваемой
задачи (на языке программирования Pascal)
uses crt;
const mf=500;
type
vector=array[1..mf] of real;
var {раздел описания переменных, которые мы будем использовать в
программе}
i, j, N : integer;
T, alfa, beta : vector;
ai, bi, ci, fi : real;
ro, c, h, tau : real;
Th, T0, Tc, L, t_end, time : real;
f, g : text;

Страницы