ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 15.10.2024
Просмотров: 94
Скачиваний: 1
РАСЧЕТ ПОДВОДНОГО ВЗРЫВА С УЧЕТОМ СЖИМАЕМОСТИ 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 бар. |