2 / 11
С.И. Климачков, В.Н. Охитин
2
Варианты расчета
Вариант
n
М
;
m
, кг
R
0
;
r
0
, м
Δ
r
, м
1
1
15
0,13
—
2
3
5
0,09
6
r
0
= 0,540
3
5
3
0,076
6
r
0
= 0,456
4
3
5
0,09
10
r
0
= 0,900
5
5
3
0,076
10
r
0
= 0,760
В численных расчетах использовали метод Лагранжа—Эйлера в
двумерной осесимметричной постановке. Для сокращения времени
счета учитывали симметрию размещения зарядов и рассматривали
лишь четверть области течения. Прямоугольная эйлерова сетка раз-
мером 6×3 м содержала 1,125·10
6
ячеек. Ось абсцисс совпадала с ли-
нией размещения зарядов, а ось ординат лежала в плоскости геомет-
рического центра разнесенных зарядов. На осях системы координат
ставили условия для жесткой стенки, а на остальных границах рас-
четной сетки — условия свободного протекания. Вдоль осей системы
координат размещали реперные точки, в которых регистрировали пара-
метры взрыва: по оси ординат 16 точек на расстоянии
r
= 5
R
0
и более, а
по оси абсцисс — 18 точек на расстоянии
r
=10
R
0
и более. Расчеты вы-
полняли с помощью программного комплекса AUTODYN.
Распространение детонации по зарядам взрывчатого вещества
(ВВ) рассчитывалось с использованием уравнения состояния JWL
для продуктов детонации [9]:
1
2
1
2
1
exp
1
exp
,
E
p A
RV B
R V
RV
R V
V
где
V
= ρ
ВВ
/ρ
0
— относительный удельный объем;
Е
— внутренняя
энергия на единицу объема;
А
,
В
,
R
1
,
R
2
, ω — эмпирические констан-
ты. Ниже приведены значения эмпирических констант, а также плот-
ности ВВ ρ
ВВ
, скорости детонации
D
ВВ
и давления в плоскости Чеп-
мена — Жуге
P
CJ
для тротила:
A
, Гпа
2026
0,3
B
, Гпа
3,23
BB
, кг/м
3
1630
R
1
4,15
D
ВВ
, м/с
6930
R
2
0,95
P
CJ
, Гпа
21,0
Для воды использовали ударную адиабату
D =
1500
+
2
u
, где
D
,
u
— волновая и массовая скорость в волне соответственно [10].
При взрыве одиночного заряда формируется сферическая ударная
волна (УВ). Расчетные значения максимального давления на фронте УВ
р
m
и удельного импульса положительной фазы давления
i
отличаются от