4 / 14
И.Ф. Кобылкин, А.В. Петюков
4
где
— массовая доля ВВ;
0
,
— начальная и текущая плотность
в расчетной ячейке;
max 1
1max 2
, , , ,
, , , , ,
, , , , ,
ig
G
I b a x
G c d y
G e g z
2min
G
— эмпирические коэффициенты, приведенные в табл. 2. Перво-
начально эти коэффициенты были выбраны такими же, как для состава
C
-4 — аналога ПВВ-12М [5]. Затем их корректировали, исходя из луч-
шего приближения к экспериментальным результатам работы [6].
Таблица 2
Коэффициенты кинетики Ли —Тарвера для ПВВ-12М
I
, мкс
–1
1,4·10
11
G
1
, Мбар
–
y
· мкс
–1
250
G
2
, Мбар
–
y
· мкс
–1
0
b
0,667
с
0,667
е
0,667
а
0,0367
d
0,333
g
0,667
х
7,0
y
2,0
z
3,0
λ
ig
max
0,022
λ
G
1max
1,0
λ
G
2max
0
Для описания сжимаемости стальных пластин использовали УРС
Ми — Грюнайзена. Механическое поведение материала пластин опи-
сывалось моделью Мизеса [4].
На рис. 2 и 3 приведены результаты расчета в виде распределе-
ний поля давления и массовой доли ПД в слое ВВ в различные мо-
менты времени. Толщина пластин составляла 7 мм.
а
б
в
Рис. 2.
Распределения поля давления в слое ВВ и пластинах, полученные с помо-
щью программы
LS-DYNA
:
а
—
t
= 3,5 мкс;
б
—
t
= 4,5 мкс;
в
—
t
= 6 мкс