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

книги из ГПНТБ / Степанов И.Р. Элементы газовой динамики и теории ударных волн учеб. пособие

.pdf
Скачиваний:
12
Добавлен:
23.10.2023
Размер:
7.63 Mб
Скачать

Пренебрегая в уравнении (5-6") бесконечно малой величиной второго порядка, получаем уравнение пзоэнтропнческого процесса

dp _ rfp

Р ~ К~ '

Таким образом, изменение состояния в ударной волне беско­ нечно малой интенсивности является нзоэнтропнческим.

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

Замечательно также то, что все основные соотношения для удар­

ных волн

были установлены

из рассмотрения

уравнений

газовой

динамики

более 75 лет тому назад, когда не было

еще

никакого

опытного

материала до того,

как

ударная волна

была

изучена

экспериментально. По образному

выражению

Э. Жуге

«ударные

волны впервые появились на кончике пера теоретиков».

 

§ 5-6. Ударная волна в атмосфере

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

60

резкое и значительное сжатие среды, распространяющееся со сверх­ звуковой скоростью. Было установлено, что устойчивой формой движения в этом случае является такая, при которой имеется фронт ударной волны, на котором происходит скачкообразное изме­ нение всех параметров. Мы изучили основные особенности фронта ударной волны.

Рис. 5-6.

Рассмотрим теперь параметры за фронтом ударной волны. За фронтом ударной волны, распространяющейся в атмосфере, как

правило, снижаются давление, скорость частиц и температура

воз­

духа. Характер

изменения давления

в фиксированной точке

про­

странства при

прохождении

через

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

показан

на рис. 5-6. В волне наблюдается

область избыточного

давле­

ния Ар = р —ро>0 длительностью т

и область разрежения с мини­

мальным давлением Дрм < 0

длительностью т_. Величинах

назы­

вается временем действия избыточного давления, или просто вре­ менем действия ударной волны (временем фазы сжатия в ударной волне). Время действия пониженного давления т_ (длительность

фазы разрежения) обычно -в несколько раз более длительности фазы сжатия.

Пространственное распределение давления в ударной волне для некоторого момента времени имеет вид, аналогичный рис. 5-6, но на оси абсцисс отложено расстояние от центра взрыва (рис. 5-7). Точка А соответствует положению фронта ударной волны, а точ­ ка В — концу фазы избыточного давления. Отрезок AB называют длиной ударной волны и обозначают К. Расстояние от В до А фронт проходит за время т+, поэтому

A= DT+.

(5-26)

При взрывах обычных взрывчатых веществ величина т+ поряд­

ка 1

мсек и %— 0,5—1,5 м. Такие

волны называются короткими.

При

'взрывах ядерных зарядов

величина т+ получается поряд­

61

ка 1 сек и л — нескольких сотен метров. Такие ударные волны на­ зываются длинными.

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

Избыточное давление во фронте ударной волны при воздушном взрыве заряда взрывчатого вещества может быть определено по эмпирической формуле М. А. Садовского:

^ = 0 , 7 6 ^ - 2 , 4 7 ^ - , 6>5 ~ \кГ,см*\, (5-27)

где G — вес заряда тротила, кГ\

R — расстояние от центра взрыва, м.

Для взрыва других веществ за G принимается вес тротила, эквивалентный по выделяющейся при взрыве энергии рассматри­ ваемого заряда. Выражение (5-27) было получено в результате обработки опытных данных.

Время действия положительной фазы давления может оцени­ ваться по эмпирической формуле

Ч = К О Ѵ'Я-ІОГ3 \сек].

(5-28)

При взрыве в атмосфере вся энергия взрыва равномерно рас­ пределяется во всем пространстве. В случае взрыва на поверхности земли (наземный взрыв) почти вся энергия взрыва распростра­ няется в полупространстве. Поэтому взрыв на поверхности земли по ударной волне в воздухе оказывается эквивалентным взрыву удвоенного заряда в безграничной атмосфере и давление во фронте ударной волны при наземном взрыве может быть вы­

числено по уравнению

(5-20) с удвоенным зарядом или

по зави­

симости

 

 

 

ДР* = 0,96 ^

+ 3,9 ^

- 13-^ \кГ;см*\.

(5-29)

62

Важной характеристикой ударной волны, определяющей сило­ вое воздействие на поверхности, является импульс давления в фазе сжатия

"+

 

Л = \

Apd^.

(5-30)

о

 

 

Для его определения может быть использована эмпирическая формула .

І+ = А ^ - ,

(5-31)

где І +— импульс фазы Ожатия, кГ • сек/см2;

R — расстояние от центра взрыва, л/; G — вес заряда тротила, кГ;

А — постоянный коэффициент, для тротила равный 35.

'Подробные данные о параметрах ударной волны получены

втеории точечного взрыва.

§5-7. Некоторые сведения из теории точечного взрыва

ватмосфере

Постановка задачи о точечном взрыве состоит в

следующем.

В начальный момент времени при т = 0 в покоящейся

среде с за­

данными начальными параметрами в точке мгновенно выделяется конечная энергия ДоПри этом возникает резкий скачок парамет­ ров газа, вызывающий движение, распространяющееся в простран­ стве симметрично во все стороны. Требуется определить возникаю­ щее при этом движение, т. е. определить все параметры на любом расстоянии и в любой момент времени. Для идеального газа при условии адиабатичности движения за фронтом ударной волны за­ дача сводится к интегрированию системы уравнений в частных производных. В Г946 году Л. И. Седовым было получено решение задачи о сильном точечном взрыве без учета противодавления. Это

решение

хорошо согласуется

с результатами

экспериментов и

с успехом

используются для

описания ранней

стадии развития

взрыва, когда давление во фронте ударной волны еще достаточно

велико.

.

При

дальнейшем развитии взрыва давление в ущарной волне

падает и влияние давления невозмущенного воздуха перед фрон­ том ударной волны становится все более существенным. Учет цротнводавления сильно усложняет задачу. Ее решение в точной по­ становке приводит к необходимости интегрирования системы урав­ нений в частных производных на основе использования численных методов на ЭВМ. Решение задачи о сильном точечном взрыве с учетом противодавления дано в работах [4, 6, 10, 11, 12, 13 и др.]. Их результаты хорошо согласуются между собой. Результаты рас­ четов приводятся в этих работах в виде таблиц и графиков. Все

63

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

 

к — показатель адиабаты;

с =

R

 

—j—j— ■— относительное расстояеие от центра взрыва;

= —» - • 3-—j2-----z-------0

относительное время от момента взрыва.

Е

РоРо

 

Здесь R — абсолютное расстояние от центра взрыва;

т— абсолютное время от момента взрыва;

Е— энергия взрыва по ударной волне;

Ро

и

Р о

— начальные давление

и

плотность

невозмущеипого

 

 

 

воздуха.

Для

расчета было

приня­

 

 

 

то

 

 

 

 

Е = 10,03 -ІО12

кГм\

р 0 —

 

 

 

=

10321 к Г /.«*■ р0=О, 125кГ ■ се^/м'

 

 

 

и к= 1,4.

При этом

важно

отме­

 

 

 

тить, что при численном реше­

 

 

 

нии задачи о точечном взрыве

 

 

 

достаточно

произвести

 

расчет

 

 

 

только

для

одного

конкрет­

 

 

 

ного

случая

 

(для

определен­

 

 

 

ных

 

Ро,

Ро и к). Затем

пу­

 

 

 

тем

пересчета

можно

получить

 

 

 

зависимость

всех

величин

от

безразмерных параметров 5 и т для данного к, после чего для того же значения к можно определить характеристики дви­ жения, вызванного взрывом при

любых Ео, Ро и fv

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

2=Л_

- ^

 

 

Рис. 5-8.

Rj_

\Е 2

(5-32)

Я,

 

 

'Е, \~5

4

(5-33)

 

т2

 

 

 

6 4

Зависимости (5-32) и (5-33) выражают собой так называемый закон подобия при взрывах. Величины энергии Е могут быть заме­ нены весом зарядов G (для атомных взрывов, а также для взры­ вов различных взрывчатых веществ — тротиловым эквивалентом).

п

R

о

х

= т называются приведенными рас-

Величины -гг-?— =

Ухпр и

у G

 

V G

 

 

стоянием и временем. Часто в зависимости от них приводятся те или иные параметры ударных волц. На р.ис. 5-8 представлена зави­ симость давления на фронте ударной волны от приведенного рас­ стояния до центра взрыва. Кривая I соответствует воздушному взрыву, когда энергия взрыва распространяется равномерно вовсе стороны, пространства, а кривая 2 — наземному взрыву, т. е. взры­ ву на поверхности земли, когда почти вся энергия распространяется в полупространство над поверхностью земли.

По данным численных расчетов давление за фронтом ударной волны во времени изменяется так, как это показано на рис. 5-9. Кривые рис. 5-9 трудно поддаются аппроксимации простыми.анали­ тическими функциями. В работе [12] дается аппроксимация точных

5 Степанов И. Р.

65

численных решений, определяющих закон изменения давления во времени за фронтом ударных волн в виде зависимости

(5-34)

Д/>Ф

где а переменная величина, определяемая следующим образом:

при

Д/?ф < 1

0,5 + Арф,

 

 

 

 

 

 

 

 

 

(5-35)

при

1 < Арф< 3

 

 

 

 

 

 

а = 0,5 + Д/7Ф 1,1 — (0,13 + 0,2Д/?ф)

 

(5-36)

При 3 <

А/?ф

 

 

 

 

 

 

 

 

а = а Н------ ----- ,

 

 

(5-37)

 

 

 

 

\ + с —

 

 

 

 

 

 

 

'4-

 

 

 

где а =

\ -

0,231 +

0,388Д/?ф - 0,0332Д/?1

для

Д/;ф<

10;

|

0

 

 

для

Арф>

10;

 

11 Ар,,,

(0,88

0,072Д/рф)

 

для

Арф <

10;

 

 

 

для

Арф >

0;

 

дР ф (1,67 — 0,011 Дрф)

с - 8 . 7 1 :

' “4 1 0 .

 

 

 

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

данные о точечном взрыве можно получить

в работах [4, 6, 10, 11,-12]. В них не учитывается изменение зави­ симости показателя адиабаты от температуры, диссоциация и иони­ зация воздуха при высоких температурах, которые наблюдаются при взрывах. Численное решение задачи о точечном взрыве с при­ ближенным учетом свойств воздуха при высоких температурах было проведено Г. Броудом. Результаты этих расчетов, относя­ щиеся к зависимости максимального избыточного давления от рас­ стояния до центра взрыва приведены в работе [13]. Они показы­ вают, что решение задачи с учетом свойств воздуха при высоких температурах дает качественно такое же изменение давления во фронте ударной волны с ростом расстояния, как и для идеального газа. В количественном отношении учет процессов диссоциации и ионизации в предположении их равновесности ведет к некоторому, уменьшению избыточных давлений.

66

§ 5-8. Отражение плоской ударной волны от жесткой безграничной стенки

Пусть ударная волна с плоским фронтом, распространяясь в .пірост.раінопве, встречает .на своем пути абсолютно жесткую без­ граничную стенку, расположенную нормально к скорости движения фронта волны (рис. 5-10). Скорость частиц во фронте падающей волны .определяется фор­ мулой (5-5):

h — Я

и,

( ^ - £

)

 

Ро

І-'ф

В момент встречи ударной волны со стен­ кой .возникает отраженная ударная волна, распространяющаяся ів обратном направле­ нии с параметрами на фронте ртр, ротр. Рас­ пространение этой волны происходит в пото­ ке падающей волны с относительными ско­ ростями D0T,, и Uorр. В данном случае на­ чальными параметрами газа будут р ф и оф. Применяя формулу (5-6) для отраженной волны, будем иметь

1

Uотр, (Ротр Рф) рф Ротр (5-38)

'9

Ротр иотр

зл

отр V,

отр

Рис. 5-10.

Так как стенка по условию остается неподвижной, то в непо­ средственной близости от нее абсолютная скорость движения газа, определяемая как разность скоростей Ц, и U0TP, равна нулю:

U* - U0Tp= 0.

 

(5-39)

Подставив в уравнение (5-39) значения 1!ф

(5-5) и £/отр

(5-38),

получим

 

 

( Р ф - р.) (— - -) =( P w - P t )

(5-40)

Уравнение (5-40) устанавливает связь между давлением и плот­ ностью во фронте падающей и отраженной волн. Оно может быть записано в виде

Ротр

Рф

1

1

 

Ро

Рф

(5-40')

Р ф -

Po

J __L

 

4

 

Рф

Ротр

 

Прибавим к правой и левой части уравнения (5-40') единицу. Тогда

Ротр

Р0

Ротр

j

ро

(5-41)

Рф

Ро

Ротр

~

 

 

рф

67

Преобразуем уравнение (5-41):

 

Р отр

Рп_

5отр

Рф

- 1

 

Ра

Рф

 

рф

 

Р.1

 

 

X - * L

 

Ротр

 

- 1

На основании

р

ф

 

рф

 

 

 

 

уравнения

Гюгонпо имеем соотношения

 

Рф

 

К - р

1

р ф

 

',

,

 

 

К - -

1

 

 

 

К +

1

Рп

 

 

.

 

Ро

 

К

1

'.

 

Рф

 

 

 

 

 

 

 

Рп

 

Porp

 

К - г - 1

Рюр .

 

 

к - - 1

 

 

 

 

Рф

 

 

 

к

- Г

1

!

Р тр

 

рф

 

К --- 1

 

 

Рф

Подставим в уравнение

(5-42)

значения I отр и —

и (5-44); после преобразовании

 

 

 

 

 

 

 

3к

 

 

 

 

 

 

Р отр

 

____ х ' Р ± _ х

 

 

к —

' Ра

 

 

 

 

Рф

 

 

1 1

 

 

 

 

 

 

k — 1

 

 

Ра

(5-42)

(5-43)

(5-44)

(5-43)

(5-45)

Формула (5-45) определяет отношение давлений на фронте отра­ женной и падающей волн. Она получена впервые .Измайловым. Формула Измайлова может быть приведена к более удобному виду. Для этого перепишем ее в виде

(

Злг — 1 Рф

 

 

Р отр - Ра

К —

1

Ра

Рп ■( Р ф — Рп)-

(5-46)

Р ф ~ Р п

,

Рф

 

 

 

« +

1

 

 

 

Отношение Р отр Рп

К — I

 

 

 

 

^Potp

=п называется коэффициентом отра-

Рф — Рп

АРф

 

 

 

 

жения. Он показывает, во сколько раз увеличивается максималь­ ное избыточное давление в волне при отражении.

Произведя

преобразования в правой

части

уравнения

(5-46),

получим

к ■

 

к -Ь 1

 

 

 

 

п —2 +

к

у І Р ф Ро)

+

К

1I АР Ф

(5-47)

 

 

2к

 

Д Р ф +

2к

 

 

(Рф

- Р п ) +

■Ра

Р о

 

68

Йз формулы (5-47) следует, что для очень слабых волн, когда Д/?Ф=ь'0, «=2, т. е. избыточное давление при отражении удваивается. По этому закону, как известно из курса физики, происходит отра­ жение звуковых волн.

В другом предельном случае, когда падающая волна имеет бесконечно большое давление на фронте (Дрф -у оо),

при к 1,4 в этом случае п 8.

Таким образом, при к=1,4 всегда 2< «< 8 . Это означает, что при отражении ударной волны происходит очень большое увеличе­

ние

давления,

особенно

для

сильных

ударных

волн. Зави­

симость

п

от

 

 

для

ро=

= 1 tcFfcM2 и /с=1,4

приведена

на рис. 5-11.

 

 

 

давления

С

увеличением

Лрф

и

Др0Тр

увеличиваются

температуры

и

Готр.

Это

влечет за собой уменьшение к.

Последнее

приводит

к

еще

большему увеличению

п.

Так,

■например, при

к=.1,2 предель­

 

 

 

 

 

ное значение«

составит я = 2 +

 

 

 

 

 

1,2+ 1

= 13. Это означает,

 

что

учет зависимости

от темпера- •

1,2 - 1

 

 

 

 

 

 

 

 

 

туры дает увеличение.коэффициента отражения.

 

 

Преобразуем выражение для 1----. Для этого подставим в фор-

 

 

 

Р

 

 

рф

 

 

мулу (5-44) значение ----отр

(5-45). После преобразований

 

 

 

Рф

 

 

Рф

 

 

 

 

 

п

 

.

 

 

 

 

 

 

к

 

 

 

 

 

" о т р

 

ро

 

(5-48)

 

 

 

РФ

(к-

D — + 1

 

 

 

 

 

 

оРотр

 

 

Ра

 

 

В пределе

 

Рф

 

к

 

к

к

п ри ----- у с

 

к — 1 . Если к = 1,4,

=3,5.

 

 

 

Ра

- - - -- -

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для отношения

температур имеем зависимость

 

 

 

Рф

 

 

Рф

 

 

 

 

 

Т,

 

 

 

(5-49)

 

 

 

 

отр

Ротр'

 

 

 

 

 

 

 

Рф

Ротр

 

 

69

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