Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Подводные и подземные взрывы сб. ст

.pdf
Скачиваний:
14
Добавлен:
22.10.2023
Размер:
15.49 Mб
Скачать

РАСЧЕТ ПОДВОДНОГО ВЗРЫВА С УЧЕТОМ СЖИМАЕМОСТИ 61

ской волне для случаев плоской и сферической расхо­ дящейся волн. Для эксперимента с плоским слоем эф­ фективная плотность при Ч. — Ж.-состоянии составляет

Рис. 1. Рассчитанные и PHERMEX экспериментальные данные по плотности для волн Тейлора. Взрывчатое вещество В-3 с плотно­

стью

1,73 г/см3.

/ —плоская волна, распространение

на 8,8 см; 2—сферическая расходящаяся

волна,

радиус 6,4 см.

-------рассчитанные кривые*,------------- экспериментальные кривые.

2,4 ± 0,05 г/см3, а экспериментальная максимальная плотность в расходящемся потоке равна 2,2 ± 0,05 г/см3. На основе у-закона или уравнения состояния БКВ было получено отличное совпадение между вычисленным нэк-

62

Ч. Л. МЕЙДЕР

спериментальным распределениями плотности в тейло­ ровской волне, что и показано на рис. 1 для передней четверти волны. Очевидно, что тейлоровское автомо-

Р и с. 2. Вычисленные и экспериментальные значения координаты ударного фронта в воде как функция времени для взрыва в воде 4,5-дюймовой сферы из вещества 9205.

Относительные единицы взяты такими же, как п в работе [9],

где

относитель­

ное время

есть

(/0 + a/D)/(a/D).

относительное

расстояние

есть

(йо + а)/а;

время

от

момента выхода

волны

детонации

в

воду, а —радиус

сферы.

D — скорость

детонации, а йо—расстояние фронта

ударной волны от

поверх­

 

 

 

 

ности сферы.

 

 

 

 

 

 

-------

рассчитанная кривая;

О экспериментальные точки.

 

дельное решение является неточным и что численное ре­ шение для расходящейся волны адекватно воспроизво­ дит действительное течение.

Другая расходящаяся конфигурация, для которой доступны экспериментальные данные, состоит из сферы взрывчатого вещества 9205 с плотностью 1,69 г/см3 при диаметре 4,5 дюйма, детонация которого осуществляется из центра при помощи 0,25-дюймового в диаметре дето­ натора. Этот заряд окружен слоем воды толщиной 13 см.

РАСЧЕТ ПОДВОДНОГО ВЗРЫВА С УЧЕТОМ СЖИМАЕМОСТИ 63

Данная конфигурация была изучена Хантелем и Дэви­ сом [9]. Давление Ч. — Ж. по уравнению состояния Б КВ для 9205 равно 280 кбар, а вычисленное максимальное давление составило 230 кбар. Рассчитанные и экспери­ ментальные координаты фронта ударной волны как функции времени даны на рис. 2. Они совпадают с точ­ ностью ошибки эксперимента.

Теоретического анализа расходящихся детонацион­ ных волн нет. Однако мы воспользовались такой число­ вой моделью, которая, по-видимому, вполне разумна, поскольку она воспроизводит имеющиеся эксперимен­ тальные данные по одномерному расходящемуся те­ чению.

IV. ЧИСЛЕННЫЙ МЕТОД

Метод SIN расчета одномерных течений в лагранжевых переменных был использован при исследовании те­ чения [4], возникающего от детонации сферического за­ ряда радиусом 3,27 см из тетрила в воде при различных начальных давлениях и плотностях среды. Для того чтобы сохранить (машинное) время вычислений в рас­ четах, использовалось только 1000 ячеек, причем вна­ чале взрывчатое вещество разбивалось через 0,1 см, а первые 10 см воды — через 0,25 см. Для вычислений при низких начальных давлениях ширина водяных ячеек была увеличена для того, чтобы иметь возможность про­ следить за распространением ударной волны в течение интересующего нас интервала времени. Это разбиение было таковым, чтобы обеспечить независимость основ­ ных характеристик течения жидкости от расчетной сет­ ки и от (искусственной) вязкости.

В этой методике используются более мелкие времен­ ные шаги и затрачивается больше машинного времени, чем это было бы необходимо, если бы применялись но­ вые неявные схемы для расчета сжимаемых течений, например такие, как метод ICE [10]. После того как неявные схемы будут получше разработаны и испы­ таны, дальнейший анализ вопросов рассматриваемого типа следовало бы проводить с использованием таких схем.

64

Ч. Л. МЕЙДЕР

 

V. РЕЗУЛЬТАТЫ

Численные расчеты подводных взрывов были выпол­

нены многими

исследователями, а работа Стернберга

и Уолкера [11] является выдающимся примером1) про­ гресса в работе в течение недавних лет. Подробное рас­ смотрение уравнения состояния продуктов ВВ, исследо­ вание расходящейся детонации и схлопыванмя пузыря для низких гидростатических давлений — в этом состоит главное отличие данного исследования от тех, которые предшествовали ему. Наши результаты отличаются в некоторых деталях (например, начальное давление в ударной волне в воде является более низким), од­ нако главные черты течения остаются такими же, какие были у Стернберга и Уолкера.

Мы выбрали для изучения сферический заряд тетри­ ла весом 0,55 фунт на различных глубинах, потому что эта взрывная система была изучена подробно экспери­ ментально и описана Коулом [12]. Давление Ч,—Ж. по уравнению состояния БКВ для тетрила с плотностью 1,70 г/см3 равно 251 кбар, а вычисленное максимальное давление в детонационной волне равно 200 кбар. Ре­

зультаты

нашего исследования

суммированы в табл.

И.

 

Суммирование результатов

Таблица

II

 

 

 

Гидростатическое

Глубина,

Период,

Максимальный

давление, бар

фут

М К С

радиус, см

1

4660

~ 156000

200

6,3

 

2

462

— 15 500

1 225

12,5

 

3

74,6

~ 2 500

5 400

23,5

 

4

9,91

~ 3 0 0

25 500

46,5

 

Рассчитанные зависимости давления и радиуса газового пузыря от времени для взрыва в воде сферического заряда радиусом 3,27 см при гидростатических давле­ ниях, равных 4660, 462, 74,6 и 9,91 бар показаны на рис. 3—10. На рис. 10 также дана экспериментальная)*

*) См. также стр. 121 настоящего сборника. Прим, ред,

РАСЧЕТ ПОДВОДНОГО ВЗРЫВА С УЧЕТОМ СЖИМАЕМОСТИ 65

кривая изменения радиуса пузыря со временем, опубли­ кованная Коулом [12]. Оказалось, что результаты рас­ чета воспроизводят данные опыта. На рис. 11 представ­ лены расчетные и экспериментальные периоды пузыря как функции от глубины в воде, причем согласие между кривыми является удовлетворительным.

На рис. 12 и 13 даны графики, соответствующие на­ чальной стадии процесса для некоторых из рассмотрен­ ных случаев. Интересно отметить, что системы с более высокими гидростатическими давлениями имеют почти совпадающие графики давление— время с аналогич­ ными графиками для систем с более низкими гидроста­ тическими давлениями вплоть до моментов времени, ко­ гда давление на границе раздела между продуктами взрыва тетрила и водой падает почти до гидростатиче­ ского давления. Градиент давления за ударной волной в воде становится менее крутым в случае систем с более высокими гидростатическими давлениями, и по измене­ нию давления можно определить, когда будет охлопы­ ваться пузырь. При более низких гидростатических давлениях имеют место многократные отражения в про­ цессе расширения и сжатия пузыря, и зависимости дав­ ление— время, давление — расстояние становятся очень сложными.

Как отметили Притчетт [13] и др. [12], согласно про­ стой модели несжимаемой жидкости, можно сделать

вывод, что

максимум радиуса пузыря обратно пропор-

 

 

 

 

Таблица III

Расчеты

Отношение

Отношение

Отношение

Отношение

максимальных

величин

периодов

величин

 

радиусов

р '/.

 

р5/е

Vi

7,38

7,78

127,5

168,6

3,72

3,60

20,8

24,6

Vs

1,98

1,96

4,7

5,4

3/ 2

1,88

1,84

4,4

4,6

V,

3,73

3,97

27,0

31,3

2/ |

1,98

2,16

6,1

6,8

3 Зак. 741

Р и с.

3. Давление

на фронте ударной

волны

в воде

(/) и

на

по­

верхности раздела

тетрил — вода

(2)

как

функция

времени

для

случая

сферы из

тетрила радиусом

3,27

см

при взрыве

в

воде

 

 

с давлением

4660 бар.

 

 

 

 

Р и с. 4. Радиус пузыря как функция времени для сферического заряда из тетрила радиусом 3,27 см при взрыве в воде с давлением

4660 бар,

Р и с. 5.

Давление

на фронте ударной

волны

в воде

(/) и

на

по­

верхности раздела

тетрил — вода (2)

как

функция

времени

для

случая

сферы из

тетрила радиусом 3,27

см

при взрыве

в

воде

 

 

с давлением 462

бар.

 

 

 

 

Рис. 6. Радиус пузыря как функция времени для сферического заряда из тетрила радиусом 3,27 см при взрыве в воде с давле­ нием 462 бар.

Р и с. 7. Давление на фронте ударной волны в воде (/) и на по­ верхности раздела тетрил — вода (2) как функция времени для слу­

чая сферы из тетрила радиусом 3,27 см при взрыве в воде с дав­ лением 74,6 бар.

Р и с . 8. Радиус пузыря как функция времени для сферического заряда из тетрила радиусом 3,27 см при взрыве в воде с давле­ нием 74,6 бар.

Р u с. 9. Давление на фронте ударной волны в воде (/) и на по­ верхности раздела тетрил — вода (2) как функция времени для случая сферы из тетрила радиусом 3,27 см при взрыве в воде с давлением 9,91 бар.

Рис. 10. Радиус пузыря как функция времени для сферического заряда из тетрила радиусом 3,27 см при взрыве в воде с давле­ нием 9,91 бар. Показаны также экспериментальные данные

Н — -)■

Рис. 11. Вычисленный (+ ) и экспериментальный (°) периоды как функция глубины воды.

Рис. 12. Давление на фронте ударной волны в воде {1) и на гра­

нице раздела тетрил — вода (2) как функция времени для ранних моментов времени.

••■давление 4660бдр* —»— давление 462 бар; -------

давление 9,91 бар.

Соседние файлы в папке книги из ГПНТБ