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

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

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

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

71

цнонален гидростатическому давлению в степени

1/3,

а период обратно пропорционален гидростатическому давлению в степени 5/6. Наши численные результаты находятся в хорошем согласии с этими выводами для

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

——-----

давление 74,6 бар; --------

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

максимума радиуса пузыря, но находятся в плохом согласии с ними для периодов, как показано в табл. 111. Периоды, полученные по расчетам, обратно пропор­ циональны гидростатическим давлениям в степени

0,877.

72 Ч. Л. МЕПДЕР

VI. ВЫВОДЫ

Согласие между наблюдениями за поведением под­ водного взрыва в течение длительного времени и по­ дробными расчетами одномерных движений сжимаемой жидкости с использованием нашего лучшего уравнения состояния для продуктов детонации взрывчатого веще­ ства и описания расходящейся детонационной волны дает основание предполагать, что рассчитанное распре­

деление

энергии

взрыва между продуктами детонации

и водой

является

достаточно точным для того, чтобы

его можно было использовать в будущих исследованиях многомерных задач о механизмах образования волн при подводных взрывах.

Хотя такое согласие и является необходимым, оно, конечно, не достаточно для исключения возможности того, что вычисленное распределение энергии является неточным.

В настоящее время усиленно изучается детонация при расходящейся геометрии и с выделением зоны реак­ ции, и это должно обеспечить лучшее понимание теоре­ тических вопросов расходящихся детонационных волн.

Список литературы

1. Mader Ch. L., Detonation properties of condensed explosives com­ puted using the Becker — Kistiakowsky — Wilson equation of sta­ te, Los Alamos Scientific Laboratory Report LA-2900, 1963.

2.Rice M. H., Walsh J. M., Equation of state of water to 250 kilobars, J. Client. Phys., 26, 824 (1957).

3.Walsh J. M„ Christian R. H., Equation of state of metals from shock wave measurements, Phys. Rev., 97, 1554 (1955).

4.Mader Ch. L., Gage W. R., FORTRAN SIN. A one-dimensional hydrodynamic code for problems which include chemical reactions, elastic-plastic flow, spalling, and phase transitions, Los Alamos Scientific Laboratory report LA-3720, 1967.

5Mader Ch. L., Comment on paper of Lee and Hornig, Twelfth Symposium (Internationa!) on Combustion, The Willliams and Wilkins Company, Baltimore, Md., 1968, p. 499.

6Курант P., Фридрихе К., Сверхзвуковое течение и ударные вол­ ны, ИЛ, М., 1950.

7. Mader Ch. L., Oneand two-dimensional flow calculations of the reaction zones of ideal gas, nilromethane and liquid TNT de tonations. Twelfth Symposium (International) on Combustion, The Williams and Wilkins Company, Baltimore, Md., 1968, p. 701.

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

73

8. Variable D.,

частное сообщение, 1968; Rivard W. С.,

Venable

D.,

Picket! W.,

Davis W. C, Flash

X-ray

observations

of marked

mass points

in explosive products,

Fifth

International

Symposium

on Detonation, ONR report DR-163, 1970.

9.Mantel L. W„ Davis \V. C., Spherical explosions in water, Fifth Symposium (International) on Detonation ONR DR-163, 1970.

10.Harlow F. H„ Arnsden A. A., Numerical calculations of almost incompressible flow, J. Comp. Phys., 3, 80 (1968).

11.Sternberg H. M., Walker W. A., Artificial viscosity method cal­ culation of an underwater detonation. Fifth Symposium (Interna­ tional) on Detonation, ONR Report DR-163, 1970.

12.Коул P., Подводные взрывы, ИЛ, M., 1950.

13.Pritchett J. W., Incompressible calculations of underwater explo­ sion phenomena, Second International Conference on Numerical Methods in Fluid Dynamics, 1970; см. русский перевод на стр. 44

настоящего сборника.

ВЗРЫВЫ ВБЛИЗИ ПОВЕРХНОСТИ ВОДЫ')

Ч. Л . М ейдер

Описываются результаты экспериментального п теоретического исследования сжимаемого течения, возникающего в результате взрыва сферического заряда взрывчатого вещества 9404, причем взрыв инициируется в центре сферы, а сама сфера частично погру­ жена в воду. Экспериментально такое течение было исследовано при помощи рентгенографической и фотографической техники. Гидродинамические расчеты с учетом сжимаемости были выполнены при помощи эйлерова численного метода для многокомпонентных сред.

1. ВВЕДЕНИЕ

Теоретическое определение волн в воде, создаваемых взрывами большой мощности, обычно основывается на экстраполяции эмпирических корреляционных зависи­ мостей, построенных по экспериментальным данным для взрывов малой мощности. Точность такого определения остается неизвестной, поэтому имеется необходимость в подробном описании механизма, который создает волны при взрывах. Надо также лучше понять влияние глубины заложения заряда. «Верхняя критическая глу­ бина», т. е. глубина, при которой имеет место пик высоты гребня волны, еще не изучена. Ранее было уде­ лено лишь небольшое внимание теоретическому и экс­ периментальному исследованию в области верхней кри­ тической глубины таких эффектов в ближней зоне взрыва (где течение является существенно сжимаемым), как водяной вал, распространение ударной волны, от­ рыв ударной волны от пузыря, а также волновой кар­ тине вблизи взрыва.

‘) Mader Ch. L., Detonations near the water surface, Los Alamos Scientific Laboratory of the University of California, LA-4958, UC-34, Los Alamos, New Mexico, 87544, June 1972.

ВЗРЫВЫ ВБЛИЗИ ПОВЕРХНОСТИ ВОДЫ

75

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

Чтобы выяснить целесообразность проведения гидро­ динамических расчетов с учетом сжимаемости при по­ мощи методов, развитых в течение последних 15 лет для описания явления детонации и уравнения состоя­ ния, мы рассчитали раннюю стадию поведения расходя­ щейся волны при детонации и сравнили найденные ре­ зультаты с экспериментальными данными, полученными рентгенографическим методом. Для ударной волны в воде, которая возникает при расходящейся детонации, результаты расчетов сравнивались также с детальными данными измерений при близких временах и с данными измерений при далеких временах, когда взаимодействие продуктов детонации и воды происходило после по край­ ней мере одного полного колебания пузыря. Обнару­ женное согласие наблюдаемого поведения течения при подводной детонации для близких и далеких времен с результатами детального одномерного гидродинами­ ческого расчета с учетом сжимаемости позволяет пред­ полагать, что рассчитанное распределение энергии меж­ ду продуктами детонации и водой является достаточно точным для того, чтобы его можно было использовать при многомерном анализе механизмов образования волн.

2. ЭКСПЕРИМЕНТАЛЬНЫЕ РЕЗУЛЬТАТЫ, ПОЛУЧЕННЫЕ НА УСТАНОВКЕ PHERMEX

В Лос-Аламосской лаборатории имеется рентгено­ графическая установка PHERMEX (импульсная рентге­ нографическая машина с высокой энергией, испускаю­ щая рентгеновские лучи), описание которой дано в ра­ боте [2]. Она представляет собой установку, в которой импульс рентгеновских лучей создается за счет удара

овольфрамовую мишень диаметром 0,7 мм электронов

сэнергией 27 МэВ, образующихся в линейном ускорителе

78

Ч. Л. МЕЙДЕР

детонации. Детонационная волна достигала наружной поверхности взрывчатого вещества через 1,5 мне. На

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

мени для

случая детонации сферического заряда радиусом 1,27 см

в воде при давлении 1 бар. Здесь

также приведены положения

этих

границ, определенные при

анализе рентгенограмм.

а —радиус ударной волны в воде; Ф —данные, полученные на установке PHERMEX; 6 —радиус пузыря; данные вдоль горизонтальной (верхняя черточка) и вертикальной (нижняя черточка*) осей рентгенограмм, полученные на уста­ новке PHERMEX.

рентгенограмме (рис. 3), снятой через 61,3 мкс после начала детонации, сферический заряд погружен в воду ка две трети диаметра. По рентгенограммам можно оп­ ределить положение границы раздела между продукта­ ми детонации и водой, всплесковую волну и ударную

ВЗРЫВЫ ВБЛИЗИ ПОВЕРХНОСТИ в о д ы

79

волну в воде. На рис. 1—3 показаны статические и ди­ намические рентгенограммы и схемы, изображающие характерные детали этих рентгенограмм.

Р и с. 5. Давление

на ударной волне в воде (кривая /) и на гра­

нице раздела

между продуктами детонации 9404 и водой (кривая 2)

как функция

времени для

случая

взрыва сферического заряда ра­

диуса

1,27 см

в воде

при давлении 1 бар.

Мы провели также одномерные расчеты по схеме SIN (аналогично тому, как это было описано в работе [1]) для случаев взрыва сферического заряда в воде при давлении 1 атм и в воздухе при атмосферном давлении в Лос-Аламосе. В уравнении состояния для воды были взяты те же значения параметров, которые были использованы в работе [1]. Применявшиеся значения

80

Ч. Л. МЕПДЕР

Таблица 1

Постоянные параметры в уравнении состояния типа НОМ

Параметр

ТЭН

 

9-104

 

А

-3,10639

+00

—2,88303

+00

В

—2,25210

+00

—2,25910

+00

С

+ 1,93865

-0 1

+ 2,09836

-0 1

D

+ 1,06761

—02

— 1,62402

—02

Е

-5,71317

- 0 5

+4,14247

-0 4

К

—М38"'0

+00

-1,27244

+ 00

L

+4,17(30

-0 1

+4,27159

-0 1

М

+4,43146

- 0 2

+4,61539

- 0 2

N

+2,43302

—03

+2,54544

- 0 3

О

+5,15057

- 0 5

+5,31474

- 0 5

0

+8,10009

+00

+8,24707

+00

R

—3,67433

—01

—4,89534

—01

S

+ 1,38196

- 0 3

+6.12169

- 0 2

Т

+8,14028

- 0 3

—3,22067

- 0 3

и—7,34294 —04 —5,13495 —06

с '

+0,5

 

+0,5

 

Z

+0,1

 

+0,1

 

Ро

+ 1,55 г/см3

+

1,844

 

^чж

+0,231

мбар

 

0,363

 

d

+0,740

см/мкс

 

0,8880

 

Тч ж

+3369

К

 

2460

 

Ртах

+0,40

 

+0,40

 

Pmln

+ 1,0

- 0 8

+

1,0

- 0 8

Давление

0,40

 

 

0,40

 

изэнтропы

 

 

 

 

 

Воздух при давлении 0,76 бар

—2,36733

+00

-1,23356

+00

+2,15170

- 0 2

-2,95528

- 0 3

+ 1,22549

-0 4

—5,53376

-0 1

+2,44880

- 0 3

+ 1,80516

- 0 2

-1,21968

- 0 3

—2,53726

- 0 5

+9,88588

+00

-2,35014

-0 1

+3,36987

—02

-4,21156

- 0 3

+ 1,63045

- 0 4

+0,5

 

+0,1

+0,0010757

+

1,0

- 0 3

+

1,0

- 0 8

 

5,0

- 0 4

параметров в уравнениях состояния типа НОМ для

взрывчатых веществ

9404,

ТЭН

и

воздуха даны

в

табл.

I 1). Радиус ударной волны

в

воде

и радиус

пу­

зыря,

определенные

в этих

расчетах,

построены

на

рис. 4 как функции времени. Здесь приведены также со­ ответствующие данные, полученные при анализе рент­

генограмм, изображенных на

рис. 1—3. Рассчитанное*)

*) Ч и сл ов ы е

 

зн ач ен и я

бол ь ш ей

части эт и х

п а р а м ет р о в

п р е д ­

ставлен ы

зд е с ь

в н о р м а л и зо в а н н о й

д еся ти ч н ой

ф ор м е: сл ева

д а н а

м ан т и сса

чи сла,

а

сп р ава —

д еся ти ч н ы й п ор я д ок .

У р ав н ен и е с о с т о я ­

ния типа

Н О М

п

р и в од и т ся

в р а б о т е [4]. — Прим, иерее.

 

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