Аналитические методы построения конечно-разностных сеток для расчета аэротермодинамики спускаемых космических аппаратов - page 18

где
M
1
число узлов на первом участке
,
включая первый узел на оси
симметрии
(
точка
A
)
и последний узел
(
точка
B
),
являющийся одновре
-
менно первым узлом на втором участке
;
z
0
осевая координата центра
кривизны носового сферического сегмента
.
Второй участок
(
BC
):
z
i
=
z
B
+(
i
1)
h
2
cos
ϕ
1
;
r
i
=
r
B
+(
i
1)
h
2
sin
ϕ
1
, i
= 2
,
3
, . . . , M
2
;
h
2
=
L
BC
M
2
1
;
L
BC
=
r
C
r
B
sin
ϕ
1
;
r
C
=
R
b
R
C
(1
cos
ϕ
1
) ;
r
B
=
R
n
cos
ϕ
1
,
где
M
2
число узлов на втором участке
,
включая точки
B
и
C
.
Заме
-
тим
,
что осевая координата точки
C
рассчитывается также по формуле
z
C
=
z
0
R
n
sin
ϕ
1
+
L
BC
cos
ϕ
1
.
Третий участок
(
CD
).
Сначала рассчитываются координаты точки
U
,
отвечающей миделеву сечению КА
:
z
u
=
z
C
+
R
c
cos
ϕ
1
=
z
0
R
n
sin
ϕ
1
+
L
BC
cos
ϕ
1
+
R
c
cos
ϕ
1
;
r
u
=
R
b
.
Тогда радиальная координата центра симметрии кругового сегмента
CD
равна
r
V
=
R
b
R
c
.
Координаты узлов сетки на поверхности этого сегмента рассчиты
-
ваются по формулам
:
z
i
=
z
u
+
R
c
cos
α
i
;
r
i
=
r
V
+
R
c
sin
α
i
;
α
i
=
h
3
(
i
1) +
ϕ
1
, i
= 2
,
3
, . . . , M
3
;
h
3
=
ϕ
1
+
ϕ
2
M
3
1
,
где
M
3
число узлов на третьем участке
,
включая граничные точки
C
и
D
.
Координаты точки
D
могут быть также рассчитаны по формулам
:
z
D
=
z
U
+
R
c
sin
ϕ
2
;
r
D
=
r
V
+
R
c
cos
ϕ
2
.
Четвертый участок
(
DE
).
Осевая координата точки
E
находится
по формуле
z
E
=
z
D
+
r
D
r
E
tg
ϕ
2
=
z
D
+
r
D
R
a
tg
ϕ
2
,
так как
r
E
=
R
a
.
ISSN 0236-3941.
Вестник МГТУ им
.
Н
.
Э
.
Баумана
.
Сер
. “
Машиностроение
”. 2004.
2 41
1...,8,9,10,11,12,13,14,15,16,17 19,20,21,22,23,24,25,26,27
Powered by FlippingBook