Этап 1. Предиктор
U
n
+1/2
j
=
U
n
j
−
−
Δ
t
2
ˉ
K
n
j
F
(
U
n
R
j
)
−
F
(
U
n
j
)
S
R
j
−
S
j
Q
21
j
+
F
(
U
n
U
j
)
−
F
(
U
n
j
)
T
U
j
−
T
j
Q
22
j
!
−
−
Δ
t
2
ˉ
S
n
i, j
G
(
U
n
R
j
)
−
G
(
U
n
j
)
S
R
j
−
S
j
Q
21
j
+
G
(
U
n
U
j
)
−
G
(
U
n
j
)
T
U
j
−
T
j
Q
22
j
!
−
−
Δ
t
2
ˉ
L
n
i, j
W
(
U
n
R
j
)
−
W
(
U
n
j
)
S
R
j
−
S
j
Q
11
j
+
W
(
U
n
U j
)
−
W
(
U
n
j
)
T
U j
−
T
j
Q
12
j
!
−
−
Δ
t
2
ˉ
N
n
i, j
V
(
U
n
R
j
)
−
V
(
U
n
j
)
S
R
j
−
S
j
Q
11
j
+
V
(
U
n
U
j
)
−
V
(
U
n
j
)
T
U
j
−
T
j
Q
12
j
!
.
(32)
Этап 2. Корректор
˜
U
n
+1
j
=
1
2
U
n
+1
/
2
j
+
U
n
j
−
−
Δ
t
2
ˉ
K
n
j
F
(
U
n
j
)
−
F
(
U
n
L
j
)
S
j
−
S
L
j
Q
21
j
+
F
(
U
n
j
)
−
F
(
U
n
B
j
)
T
j
−
T
B
j
Q
22
j
!
−
−
Δ
t
2
ˉ
S
n
i, j
G
(
U
n
j
)
−
G
(
U
n
L
j
)
S
j
−
S
L
j
Q
21
j
+
G
(
U
n
j
)
−
G
(
U
n
B
j
)
T
j
−
T
B
j
Q
22
j
!
−
−
Δ
t
2
ˉ
L
n
i, j
W
(
U
n
j
)
−
W
(
U
n
L
j
)
S
j
−
S
L
j
Q
11
j
+
W
(
U
n
j
)
−
W
(
U
n
B
j
)
T
j
−
T
B
j
Q
12
j
!
−
−
Δ
t
2
N
n
i, j
V
(
U
n
j
)
−
V
(
U
n
L
j
)
S
j
−
S
L
j
Q
11
j
+
V
(
U
n
j
)
−
V
(
U
n
B
j
)
T
j
−
T
B
j
Q
12
j
!
.
(33)
Этап3. Учет правой части
˜
U
n
+1
j
= ˜
U
n
+1
j
+ Δ
t P
(
U
n
j
)
.
(34)
Этап 4. Сглаживание
U
n
+1
i
= ˜
U
n
+1
i
+
α
(
U
n
R
j
+
U
n
L
j
−
4
U
n
j
+
U
n
U
j
+
U
n
B
j
)
,
(35)
где
α
— коэффициент искусственной вязкости.
Моделирование процесса динамического разрушения и прони-
кания ударника.
Приведенные разностные соотношения (31)–(34) по-
зволяют рассчитывать все характеристики ударника и мишени при
108
ISSN 1812-3368. Вестник МГТУ им. Н.Э. Баумана. Сер. “Естественные науки”. 2008. № 4