книги из ГПНТБ / Кушнарев Д.М. Использование энергии взрыва в строительстве
.pdfХотя сама форма цилиндра выброса при таком подходе не определяется, выемку можно условно принять как внутреннюю область, ограниченную тон линией тока, которая в месте пере сечения с поверхностью грунта определяет край выемки. Полу чаемый таким образом цилиндр выброса будем называть услов ным, и, как вытекает из энергетических соображений, он будет меньше фактического.
Чтобы определить форму истинного цилиндра выброса, нуж но принять более серьезные ограничения и отказаться от рас смотрения всего грунта в целом как несжимаемой идеальной жидкости. Первоначальную постановку задачи определяет, как известно, параметр Л, характеризующий действие взрыва на движение среды. Обычно в качестве параметра П выбирают ли бо энергию взрыва Е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. Определяем поле скоростей:
_ |
2Р |
|
1 |
( Л \ |
2Ур |
|
dx |
яр |
1 + |
( J L f |
х |
я (.ѵ= + у-) Р |
|
|
|
|
|
|
||
d<p |
|
|
|
|
2хР |
(1.56) |
|
|
|
|
|
|
Уяр(.ѵ2 + г/=) '
1+
Вчастности, на свободной поверхности скорости вылетаю щих частиц равны:dij яр
Ѵ , І І Р = О = 0, |
ѵ ^ = о |
2Р |
(1.57) |
|
лхр |
||||
|
|
|
и координата точки В', являющейся краем условного цилиндра выброса, определяется из условия:
|
D' |
(1.58) |
|
|
2 |
||
|
|
||
D' = |
4P |
(І.58а) |
|
ярС |
|||
|
|
||
На рис. 10 показаны поле скоростей, эпюра скоростей |
вылетаю |
щих частиц и условный цилиндр выброса при расположении
ленточного заряда перпендикулярно |
открытой поверхности |
||||
(см. рис. 2, а). |
Скорости |
вылетающих |
частиц меняются |
по ги |
|
перболическому |
закону. |
|
|
|
|
Абсолютная величина |
скорости в произвольной точке |
|
|||
|
|
V = |
2Р |
|
|
|
|
рлг |
|
|
|
|
где г |
|
|
у , |
(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