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

книги из ГПНТБ / Кушнарев Д.М. Использование энергии взрыва в строительстве

.pdf
Скачиваний:
11
Добавлен:
27.10.2023
Размер:
19.55 Mб
Скачать

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

Чтобы определить форму истинного цилиндра выброса, нуж­ но принять более серьезные ограничения и отказаться от рас­ смотрения всего грунта в целом как несжимаемой идеальной жидкости. Первоначальную постановку задачи определяет, как известно, параметр Л, характеризующий действие взрыва на движение среды. Обычно в качестве параметра П выбирают ли­ бо энергию взрыва Е0, либо импульс /0 , полученный средой при взрыве. Параметр Е0, как правило, выбирают тогда, когда взрыв считается сильным. Такой случай подробно рассмотрен акад. Л. И. Седовым.

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

Так как скорость распространения ударной волны значи­ тельно превышает скорость расширения цилиндра выброса, то

условие несжимаемости

в непосредственной близости от

заря­

да действительно

будет

выполняться.

Воспользовавшись

ста­

ционарным уравнением

Бернулли в

критической точке

Я к р +

9

 

 

 

 

- \ — ^ — const,

можно

утверждать,

что предположение

о су­

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

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

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

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

20

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

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

Рассмотрим в такой постановке плоскую задачу об образо­ вании цилиндра выброса при взрыве горизонтального цилиндри­ ческого заряда радиуса R, заложенного на глубину h от поверх­ ности грунта (рис. 6). В силу симметрии достаточно рассмот­ реть только правую половину картины, которую будем назы­ вать областью Д.

Требуется найти гармоническую в области Д функцию српри следующих граничных условиях.

На границе круга Г действует постоянное вдоль границы им­ пульсное давление, поэтому на этой границе значение потенциа­

ла

постоянно:

 

 

 

 

Ф =

Фо =

. zer,

(1.39)

 

 

 

Р

 

где

Р — импульсное давление при взрыве,

равное по определе­

нию

 

т

 

 

 

 

 

 

P =

lim

\p(t)dt.

(1.39а)

На свободной поверхности давление Р, конечно, равно ат­ мосферному, поэтому импульсное давление равно нулю. Следо­ вательно,

(х,у = 0) = 0.

(1.40)

На неизвестной границе области Д выполняются два условия. Первое условие

да л = 0(п — внешняя нормаль к границе области) (1.41)

21

означает,

что граница цилиндра

выброса есть твердая

стенка,

а второе

условие

 

 

 

dq>

с

(1.42)

 

~dS

 

 

утверждает, что вдоль границы Л скорость скольжения по­ стоянна и равна критической С (С — единичный вектор, направ­ ленный по касательной к границе выемки). Кроме того, извест­

но, что на свободной

поверхности уфО и на

отрезке ,v=0, h—•

— / ? < / / < 0 скорости

направлены вверх,

на отрезке

х=0,

y>h-\-R— вниз, а на поверхности цилиндра Г — радналыю. Вследствие симметрии картины взрыва относительно оси у огра­ ничимся рассмотрением только правой половины области Д.

Если ввести функцию тока

то вместо условия (1.42) бу­

дем иметь:

 

ty(x,y)\A

= const,

или, выбирая на границе выемки константу, равную нулю, по­ лучим:

* | л = ° -

(І.42а)

 

Введем безразмерные переменные:

 

X = R

 

 

Ф = — . С =

cR_

(1.43)

 

R '

Фо

J

Фо

 

 

 

 

 

и комплексный потенциал

W(z) = ф + іяр.

Поэтому краевая задача формулируется следующим образом: найти границу неизвестной области Д, внутри которой опреде­ лена аналитическая функция W(z), граничные значения кото­ рой имеют следующий вид (черточка означает комплексную со­ пряженную величину):

R e r = l , arg

dW_

arg (z + ih)

при ггГ

x

 

 

 

 

dz

 

 

 

 

 

 

 

 

 

X{z:|z +ih\

=R\

;

 

 

 

 

ReW=0,

a r g ^

= - | при y = 0,

0 <

x <

-|-

 

I m W = 0 ,

d\V

= С на неизвестной границе Л.

(1.44)

 

dz

 

 

 

 

 

 

 

 

< p ' < R e l F < l t

dW

 

 

при

x = О,

 

a r g — =

 

 

 

-H<y<~(h

dz

+ R);

 

 

 

 

 

 

 

 

 

 

О <

Re И7< 1, arg

dW

я

при x

= О,

 

 

= ~

 

 

 

,

arg

dz

ni

 

 

 

 

 

 

 

&

2

1

 

 

 

 

 

 

•(h

 

R)Ky<0.

 

 

 

 

22

Таким образом, задача сводится к отысканию некоторой ана­ литической в области Д функции, действительная и мнимая ча­ сти которой принимают на границе области заданные значения,

т.е. к смешанной краевой задаче.

В1937 г. М. В. Келдыш и Л. И. Седов [14] вывели форму­ лу, которая позволила эффективно решать смешанную краевую задачу для гармонических функций, благодаря чему она широ­ ко применяется в практических расчетах.

Смешанная задача решена М. В. Келдышем и Л. И. Седо­ вым только для односвязной области. Ограничимся случаем, когда область Д представляет собой верхнюю полуплоскость; к нему с помощью конформного отображения сводится, очевид­ но, случай произвольной односвязной области.

Пытаясь

применить формулу Келдыша — Седова для реше­

ния задачи

о взрыве мелкозаглубленного горизонтального ци­

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

мо

обобщить

формулу Келдыша — Седова на случай двусвяз-

ных

областей.

Чтобы обойти указанную трудность, воспользу­

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

что этот

метод, несмотря на свою эффективность, приводит

к весьма

громоздким вычислениям.

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

Предположим, что изучается взрыв на поверхности грунта. Решить такую задачу можно, так как здесь применима формула Келдыша — Седова. В результате находим форму выемки и, в частности, ее край (точка В, см. рис. 6). Затем, считая грунт идеальной несжимаемой жидкостью, определяем условный ци-

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

D

*

пами, можно показать, что отношение а=-^-

всегда будетооль-

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

23

В самом деле, при взрыве зарядов выброса в грунте основ­ ная часть энергии расходуется на выброс и только ничтожная ее доля — на деформацию и другие потерн. Если взрыв проис­ ходит на поверхности жидкости, то энергия распределяется та­ ким образом, что только часть ее расходуется на образование условного цилиндра, а остальная ее часть — на сообщение ки­ нетической энергии жидкости, находящейся за пределами ус­ ловной воронки. Поэтому всегда D'<CD и а > 1 .

Далее, очевидно, что в первом приближении D и D' обратно пропорциональны С и прямо пропорциональны энергии взрыва. Поэтому физическая характеристика грунта С и параметр, ха­ рактеризующий действие взрыва на грунт Р, выпадают из отно-

шения

D

у,

a ~ j j r -

Наконец, как был о показано выше, в точках

с координатами D, D'> R конформное отображение, при по­ мощи которого некоторая форма заряда сводится к стандарт­ ной, является тождественным и форма заряда влияет только на вид течения вблизи самого заряда, а ие на место расположения

края выемки. Поэтому отношение а = ~^г будет лишь незначи­ тельно зависеть от указанных выше характеристик. Строгое до­ казательство высказанного утверждения будет дано ниже, где в результате конкретных вычислений будет также найдена ве­

личина а.

Смысл величины а можно выяснить следующим образом. Пусть при взрыве на поверхности грунта и на поверхности не­ сжимаемой жидкости выделяется одно и то же количество энер­ гии Е0. Объем цилиндра выброса, рассчитанный на единицу дли­ ны в продольном направлении, будет в первом случае пропор­

ционален

D2,

а во втором случае

~(D')2.

Однако только у-я

часть энергии

Е0

израсходовалась

па

образование условного

цилиндра

выброса,

поэтому

 

 

 

 

 

 

 

 

ш . = а

, = & = у <

1 .

( І . 4 5 )

Итак, а=

У

у, т - е

-

квадрат

отношения

а дает ту долю

энергии,

которая расходуется на образование условной выемки по отно­ шению к полной энергии взрыва.

 

Как видно из определения,

величина

а

имеет

геометричес­

кий смысл и должна зависеть

только от

геометрических

пара­

метров. За исключением D и D', мы имеем следующие парамет­

ры

длины: R, Іг и Я

(см. рис.

6). Глубина

условной

выемки h'

по

ее определению

совпадает

с глубиной

заряда

h.

Так

как

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

A.

JL

А.

H

' H '

h '

24

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

/ h

а = а —

U

Мы будем считать заряд мелкозаглубленным, т.е. k малым

по сравнению с Я, и заменим

а его приближенным

значением:

а =

а(0).

(1.46)

Но даже и в тех случаях, когда H не намного превышает h, а не должна сильно меняться, так как по смыслу эта величи­ на— вялая функция от — . Фактически а есть универсальная

H

константа, которую мы и определим в дальнейшем, а сейчас предположим ее известной величиной. Тогда по формуле

 

 

 

 

D(h)

= ^ а-

 

 

 

 

 

(1.47)

определяем край цилиндра выброса при взрыве

мелкозаглуб-

ленного заряда в грунте на глубине

h

и находим

координату

точки

В (см. рис. 6). Затем,

пользуясь

решением

задачи

о по­

верхностном взрыве в грунте, определим, какая должна

быть

величина поверхностного

заряда, чтобы

получить

край

выемки

в нужном

месте,

т. е. в

конкретной

точке В,

находящейся на

расстоянии

D(h) — а

от начальной точки,

и по уже

извест-

ным

формулам

рассчитаем

форму

цилиндра

выброса.

 

Здесь

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

5. ОПРЕДЕЛЕНИЕ УСЛОВНОГО ЦИЛИНДРА ВЫБРОСА

 

ПРИ ВЗРЫВЕ

ЗАРЯДОВ РАЗЛИЧНЫХ

КОНСТРУКЦИЙ

 

Рассмотрим сначала

заряд в виде

тонкой

ленты

бесконечной

ширины (/->- оо), расположенной вертикально

(см. рис. 2). Если

рассматривать

грунт

как несжимаемую идеальную

жидкость,

то в результате

взрыва начинают

двигаться

все частицы грун­

та и задача сводится к частному случаю, рассмотренному ра­

нее при С->0. При этом краевые условия

(1.44) принимают вид:

Re Q = 1, arg —

= 0 при х =

0,

_ о о < г , < 0 ;

dz

 

 

(І.44а)

 

 

 

Re Q = 0,

arg - ^ - = у

при

у = О,

25

где Q — величина, комплексно сопряженная по отношению к

комплексному потенциалу Q=W.

В плоскости комплексного потенциала ÇÏ = u-\-iv область те­ чения занимает бесконечную вертикальную полосу шириной Im (рис. 7).

Будем искать в этой области

/

г

,dQ

In

dz

Re/(Q) = In dQ I у - d

— , Im/(Q) = arg dz I

J?=U *IV

V,'j»f(sï)=o

У

'/

(1.48)

Q

dz

7

Рис. 7. Область течения для бесконеч­ ной вертикальной полосы

Рис. 8. Отображение точек на ниж­ нюю полосу

Л = Л Х « - %

77777777777?77777777777777777777Я7777777777777Ъ

Для функции /(Q) на границах полосы выполняются следу­ ющие граничные условия:

Im / (Q) = — при

Ф = О,

 

 

 

2

 

(І.48а)

Im/(Q)

= 0 при

ф = I . J

 

Функция h(Q)=emP~

переводит

данную полосу в

верхнюю

полуплоскость, причем

Q =

0 переходит в A = l , Q = l — в

Іі = — I

(рис. 8), и на положительной полуоси lmf(Q) — ~ , а на отри-

26

цательной I m / ( Q ) = 0 . Наконец, при помощи дробно-линейного преобразования переводим верхнюю полуплоскость снова в

верхнюю полуплоскость, но с нужным

соответствием точек

(рис. 9):

 

^ = т т \ и л и

( L 4 9 )

Вводя вместо f(Q) функцию

 

F (Q) =—if

(1.50)

приходим к следующей задаче: найти аналитическую в верхней полуплоскости функцию F[Q{'Q)] = Ф ( £ ) , действительная часть

Рис.

9. Соответствие

точек

ß(- оо)

D

2 f

Re<P=Q

при

переходе

отображаю­

 

 

щей функции в верхнюю по­

 

 

 

 

 

 

луплоскость

 

 

 

 

 

 

которой на

вещественной оси принимает

следующие

 

значения:

 

R e Ф ( 0 =

0 при о о < £ < — 1 ,

1 < £ < о о ; ]

(1.51)

 

Re Ф (£) =

-£- при — 1 < £ < 1 .

 

 

 

 

 

 

 

Теперь

воспользуемся

формулой Келдыша — Седова [14],

которая сводится в данном случае

просто к интегралу

Шварца.

В самом деле, мнимая часть искомой функции на действитель­

ной оси равна нулю, а ее действительная

часть отлична

от нуля

только на отрезке от —1 до + 1 . Функция g(£,) =

h поэтому

 

1

 

 

 

 

Ф (О =

= — I n ъ

+ Ф(оо),

(1.52)

Так как при £ > 1 , К е Ф ( £ ) = 0 ,

то константу

Ф(оо)

следует

 

 

 

 

— Qj и воз­

вращаясь к функции

f(й),получим:

 

 

 

 

 

і ig ( — Q ) — 1

 

 

/ ( Q ) = - L l n

ï—L

.

 

(1.53)

 

 

+ 1

 

 

Отсюда, пользуясь

(1.48), имеем:

 

 

 

 

 

=

іе

 

 

27

Интегрируя это выражение, получаем решение задачи в сле­ дующем виде:

• П+

К.

 

(1.54)

Из условий (1.54) находим, что постоянная

/ ( = 0 .

Таким об­

разом, получаем окончательно:

 

 

 

Q = J _ l n / _ J L 2

ф

" I ' ,

(1.55)

tit

 

 

 

откуда, переходя к размерным переменным, имеем:

Ф— arg 2 - f я = — arctg H

Л

( л

яр

д:

р

(І.55а)

4P ,

2

 

 

 

ІП

 

 

 

 

яр

I 2

 

 

 

 

т. е. течение представляет собой плоский вихрь с центром в на­ чале координат 0. Определяем поле скоростей:

_

 

1

( Л \

2Ур

 

dx

яр

1 +

( J L f

х

я (.ѵ= + у-) Р

 

 

 

 

 

 

d<p

 

 

 

 

2хР

(1.56)

 

 

 

 

 

 

Уяр(.ѵ2 + г/=) '

1+

Вчастности, на свободной поверхности скорости вылетаю­ щих частиц равны:dij яр

Ѵ , І І Р = О = 0,

ѵ ^ = о

(1.57)

лхр

 

 

 

и координата точки В', являющейся краем условного цилиндра выброса, определяется из условия:

 

D'

(1.58)

 

2

 

 

D' =

4P

(І.58а)

ярС

 

 

На рис. 10 показаны поле скоростей, эпюра скоростей

вылетаю­

щих частиц и условный цилиндр выброса при расположении

ленточного заряда перпендикулярно

открытой поверхности

(см. рис. 2, а).

Скорости

вылетающих

частиц меняются

по ги­

перболическому

закону.

 

 

 

 

Абсолютная величина

скорости в произвольной точке

 

 

 

V =

 

 

 

 

рлг

 

 

 

где г

 

 

у ,

(1.59)

28

также меняется по гиперболическому закону. Кинетическая энергия, заключенная в кольцевом слое толщиной dr, очевидно, равна:

dEr г = 2 шаг,

где р — плотность энергии.

Энергия, приходящаяся на полукруг радиуса г с центром в точке О, Е,- пропорциональна

и растет логарифмически, а плотность энергии убывает как

Е = -ML. ~

l£L f

(1.60)

 

яг-

г2

 

т. е. несколько медленнее,

чем по закону обратных квадратов

В этом наиболее простом

случае

условный цилиндр

выброса

имеет форму круговой выемки (см. рис. 10).

 

Рассмотрим случай, изображенный на рис. 2,6. Задача фор­ мулируется следующим образом. Требуется найти аналитиче­

скую в нижней полуплоскости z функцию W(z),

которая на гра­

нице области принимает следующие значения:

 

Re W (z) = 1 при у = 0, — 1 < X <

1;

ReW{z)

= 0

при у = 0, \х\>

1.

(1.61)

Решение дается

интегралом

Шварца

 

(1.52), который прини­

мает в данном случае следующий вид:

 

 

 

 

 

1

 

 

 

 

W ( z ) = —

f -Л—

= ±-ln-^î-.

(1.62)

 

я'

J t — z

711

z -f- 1

 

 

 

 

 

 

 

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

29

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