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

UptoLike

94
Rewrite(g);
for i:=1 to N do
writeln(g,' ',h*(i-1):6:3,' ',T[i]-273:8:5);
close(g);
end.
Получены следующие распределения температуры:
Рис. 22. Распределения температуры по толщине пластины в различные
моменты времени
Рассмотрим чисто неявную схему.
()
λ
λ=
λ
+
+
+
++
+
+
+
h
TT
h
TT
hx
T
T
x
n
i
n
i
n
i
n
i
n
i
n
i
1
1
1
1
21
11
1
1
21
1
,
где
.
2
,
2
11
1
1
21
1
1
1
1
21
++
+
+
+
+
+
+
λ+λ
=λ
λ+λ
=λ
n
i
n
i
n
i
n
i
n
i
n
i
Таким образом, в результате аппроксимации частных
производных соответствующими конечными разностями получаем
следующую систему линейных алгебраических уравнений:
,0 ,1,,2
,
1
1
1
1
1
21
11
1
1
21
1
=
λ
λ=
τ
ρ
+
+
+
++
+
+
+
+
nNi
h
TT
h
TT
h
TT
с
n
i
n
i
n
i
n
i
n
i
n
i
n
i
n
i
K
(41)