Файл: Подводные и подземные взрывы сб. ст.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 15.10.2024

Просмотров: 96

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

340 ДЖ. Р. МЁРФИ

Подстановка (27) в

(20) дает горизонтальную компо-

.ненту смещения

J A0(q) е 1квм q h (1 - ф ) ~ Ъ dq, (28)

;' ,<Л. = е*W4) (k/2nrf

 

-Do

где В (q) = rq -f- 2 Я (1 — q2)'1*, и мы использовали равен­ ство е-Дзл/4) _ _ е( (я/4)# в выражении для В мы теперь

заменим переменную q на р, используя (26), и напишем

В (р) = г [q\ — 2 (1 — <7jj)1/s ре-' (я/4) + /р2],/а +

+ [1 - q\ + 2 (1 - <7*)V ре-' <я'4>- <p8]V*. (29)

Разлагая В(р) в ряд около р = 0, получаем ikB (р) = ikB (<7о) +

+ ^ p [ 2 ^ - r ( l - p ^ / p 0] e - 'W ) - ^ Р 2+ .... (30)

где последующими членами этого выражения можно

пренебречь

для

больших удалений;

при qo — sin е =

= r/RR, где

Rr =

(г2+ АН2)к (рис.

1),

имеем

B(q0) = RR,

2 H - r ( l - q $ b l q 0 = 2 H-(R*R-r*)'b = o. (31)

Таким образом, экспоненциальный член можно пред­

ставить в виде

 

 

 

 

ikB (р) ~ ikRR -

р \

(32)

Выразим теперь подинтегральное выражение в (28) че­

рез переменную р, используя

(32) и тот факт,

что

<7(1- q2)''h dq =

- e ~ ‘im dp.

(33)

Тогда получим (помня, что пределами интегрирования вдоль D0 являются оо и — оо)

00 •

-‘ ' и \ = Я/ к**(к12пгр0)'1> | А0(р) exp (— krp2l2ql')dp. (34)

— ОО

Из (34) ясно, что подинтегральное выражение будет быстро уменьшаться, когда |р | увеличивается, и, сле­ довательно, только значения .р, близкие к нулю (т. е.



РАСЧЕТНЫЕ ВСТУПЛЕНИЯ ВОЛН СЖАТИЯ

341

q qо), будут давать значительный вклад. В большин­ стве приложений Яо(<7о) будет слабо изменяющейся функцией <7о, за исключением окрестности qo = n. Сле­ довательно, если <7о не близко к п, мы можем А0(р) вы­ нести за знак интеграла, как и Aofpo), и использовать тот факт, что

Jоо exp ( -krp2/2ql) dp = (2пф^кг)\

(35)

чтобы получить

_

lk R R

 

Ч . = ? , Л М

\ - .

(36)

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

 

 

_

i&Rft

 

 

 

 

R

(37)

 

 

 

 

Для

окрестности qo = п мы используем (23), чтобы пе­

реписать (34)

в форме

JP (p) exp ( - krp2/2q§ dp

 

и°н, =

q / kRR {k/2nrq0)'k j

 

 

 

Vo o

 

 

 

-

J Q(p)(n2 — q2)'k exp(-krp°-/2ql)dp\. (38)

Для большинства приложений P(qo) и Q{qo) будут

достаточно слабо изменяющимися функциями qo, так что можно положить Р(р) ж Р(0) =P(q0), Q(p) » « <3(0) = Q{qo) в (38) и вынести эти члены за знак

интеграла. Тогда

 

ikR R (

U%, = %

( Р Ы ~ RrQЫ (k/2nrq0)'laX

 

ОО

 

X J («2 — q2)'h exp ( - krp2/2q$)dpJ . (39)


342

ДЖ. Р. МЁРФИ

Вводя переменнуюг/°=(/гг/2^)/з [(1 — n2)h—(^—ql)h] и используя соответствующую аппроксимацию для квад­

ратного корня (/г2 — р2)'/г в уравнении (39), Червени [3] получил следующее выражение для поля отраженной

Р и с. 3. Схема пути пробега головной волны Р„

волны сжатия, которое справедливо для значении qo, близких к п (и включающих п):

где

( Д. =■V

К to) - ( т Г « w о <Л ]. (40)

 

 

 

1 =

 

 

 

 

“ЗГ*

 

 

 

g (п ) =

n h ( I -

n i h Q (n)

 

 

(4я )1

(41)

 

 

 

 

 

 

G (y°) = -fir

[ [г/° — p exp ( - ш/4)]’АX

 

 

К Я

*J

 

 

X exp ( - p’2)dp — г'2у' / p°.

РАСЧЕТНЫЕ ВСТУПЛЕНИЯ ВОЛН СЖАТИЯ

343

Сравнение уравнении (40) и (36) показывает,

что

член в скобках в (40) можно интерпретировать как

ко­

эффициент отражения сферической волны, причем до­ полнительное слагаемое определяет зависящую от ча­ стоты поправку для коэффициента отражения плоских волн. Благодаря уравнению (40), справедливому вблизи q0 =zii, и (36), справедливому всюду в остальной обла­ сти, мы имеем полное решение для смещения в дальней зоне, связанного с волнами типа РР от гармонического источника. Однако для q0 > п к общему полю отражен­ ной волны добавляется вклад, связанный с интегрирова­ нием вдоль D* и соответствующий преломленной волне Рп в критической области. Для значений qo, близких к п (но больших п), поле волны Рп будет интерферировать с полем РР, и наблюдаемое смещение будет результатом суперпозиции смещений двух волн. Для q0 п времена пробега этих двух волн будут существенно различаться, так что интерференция для нестационарного источника уже не будет иметь места.

Вводя переменную у" = L \kn (1 — п2)/2г]'Гг, где L — расстояние, проходимое волной Р„ параллельно границе (рис. 3), Червени [3] получил следующее приближенное решение для интеграла вдоль D*, которое справедливо для значений q0, близких к п:

 

 

г ; . = 7

( I - ) ' ' ^

f t

(42)

где

 

 

 

 

 

 

gi (!/*)=

 

 

 

 

 

I

 

 

 

со

 

 

I

s. pi (7л/8) .

 

 

I

= 2 1

V

r -

| / р е х р [ - р 2— /2 д (1 + г )г /* Ц э .

| (43)

 

я

■’

 

 

 

В (п) = nL +

(Z +

Z0) sec е*,

sine’ = H.

 

Когда L велико, можно перейти к асимптотическому со­ отношению

£/* =

in2Q (п)

oikB in)

(44)

kYTL'4\-n*)'h

 

 

которое часто использовалось для расчетов [8, 16].


344

ДЖ. Р. МЕРФИ

В окрестности qo — п волновые поля РР и Рп будут интерферировать. Следовательно, чтобы получить сум­ марное поле в этой области, необходимо сложить обе компоненты (40) и (42). Тогда получим

 

i k R —

/ 5

\V

 

 

и\= ?о-V - [Ло Ы -(-Jfj'g W G^

+

 

 

,

я*

(1-«*)‘'7

я \ ч,п в м

, ,ч

’ ]• (45)

 

+

-

— Г ~ ( т )

 

где у = [kRR kB (/z)]v\

В окрестности

критической

точки

имеем

у ~

у0 ж у*

для qo ^

п и i/s; —у0 для

q0< n .

Следовательно,

суммарное

поле

отраженной

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

в виде [3]

IkR

 

- (4Г sC ‘)0 to)].

 

 

 

 

UR, = % п г 1 К Ы

(46)

 

'R

 

 

 

 

{ - [ k R

km 2Hk (1 — n2)'/i]',\

q0< n,

(47)

\ +

[kR -

km - 2Hk{ \ - rt2)'/j]v’,

<7o > n .

 

Наконец,

так как

 

 

 

G (у) = V2 exp [/ (я/8 -

уЩ] D4t [у (i - 1)] -

Й*А/ у ,

(48)

где Dу, — протабулированная функция Вебера поряд­ ка '/г, уравнение (46) представляется в легко реализуе­ мой расчетной форме.

Таким образом, уравнение (46) вместе с уравнения­ ми (44) и (36) дает полное описание основных волн сжатия, образуемых при отражении сферической волны сжатия от плоской границы, разделяющей два твердых полупространства. Из (36) можно видеть, что для qo «С п отраженная волна будет иметь ту же самую временную зависимость, что и падающая волна, поскольку величи­

на A0(qo) не зависит от

частоты и

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

для

qo < п. Таким образом,

например,

горизонтальная

со­

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