Ю.И. Димитриенко, И.Д. Димитриенко
14
1/2
1/2
1/2
1/2
1/2
11
12
(
) (
)
(
) (
)
2
j
j
j
j
n
n
n
n
j
j
L
B
n
j
j
j
j
L
j
B
t
Q
Q
S S
T T
V U V U
V U V U
N
.
n
j
t
P U
(25)
Шаг 3. Искусственная вязкость.
Для устранения эффекта нефизической осцилляции, возникающе-
го вследствие разности значений в четных и нечетных ячейках сетки,
введем в задачу искусственную вязкость:
1
1
n
n
n
j
j
j
U U L U
,
(26)
где
4
;
j
j
j
j
n
n
n
n
n
n
j
j
R
L
U B
L U U U U U U
α — коэффициент искусственной вязкости, который выбирался пе-
ременным в зависимости от направления фронта ударной волны.
Моделирование процесса динамического разрушения прегра-
ды из ГБКМ.
Представленные выше уравнения (17)—(21) позволяют
рассчитывать все характеристики ударника и мишени при непрерыв-
ных процессах деформирования, т. е. при отсутствии разрушения.
Однако сам процесс проникания ударников в преграду как раз и со-
стоит из стадий зарождения трещин, их распространения и образова-
ния отверстия в преграде. Вопросам численного моделирования рас-
пространения трещин в композитах посвящены работы [22, 23]. При
численных расчетах динамический процесс распространения трещин
моделируется в два этапа.
1. Этап зарождения трещин.
В процессе непрерывного дефор-
мирования в каждой расчетной точке ударника и мишени проверяют
условие отсутствия разрушения:
1,
где
— параметр поврежда-
емости материала, представляющий собой функции от инвариантов
тензора истинных напряжений Коши (силовой критерий):
( )
( )
O
I
T
.
Инварианты выбирают в соответствии с группой симметрии рас-
сматриваемого материала (изотропные или ортотропные). В этом
случае критерий прочности преграды и критерий пластичности по-
добны, в них пределы текучести заменяют на соответствующие пре-
делы прочности.
Для ортотропных ГБКМ, на основе тканевых структур будем
считать, что разрушение происходит только в случае, если осуществ-
