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

Основы_акустики_Гринченко_Вовк

.pdf
Скачиваний:
62
Добавлен:
10.04.2023
Размер:
16.81 Mб
Скачать

Рассмотрим случай, когда напряжение на поверхности описывает- ся одной пространственной гармоникой exp(iβx1) с круговой про- странственной частотой β, и определим волны, которые порождаются

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

Итак, пусть на границе x3 = 0

(рис. 6.20) заданны механические

напряжения

σ33 = exp(iβx1).

 

σ31 = σ32 = 0,

(6.178)

Круговая пространственная частота β связана с пространственным периодом L соотношением β = 2π/L .

Волны, возникающие в среде, описываются векторной функцией смещения u (x1,x3 ), которую можно представить в виде суммы

u = ue + ut , где ue и ut - решения волновых уравнений (6.116). Век- торы: ue = e1ue(1) + e3ue(3), где ue(1) и ue(3) проекции вектора ue на оси координат Ox1 и Ox3 соответственно, ut = e1ut(1) + e3ut(3), где ut(1) и ut(3) проекции вектора ut на оси координат Ox1 и Ox3. Амплитудно- фазовые соотношения между векторами смещений ue и ut должны

быть такими, чтобы удовлетворялись граничные условия (6.178). Рассмотрим волновые уравнения относительно проекций векторов

ue и ut . Понятно, что форма этих уравнений одинакова:

 

u u(n )

или u u(n ) ,

u +k2u = 0 ,

 

 

(6.179)

где

n =1,3 ; k = ω/c

e

или k = ω/c

t

соответствен-

 

e

t

 

 

 

но. Будем искать решение уравнения (6.179) в виде

 

 

 

 

u (x1,x3 )exp(iωt ) = g (x3 )exp(iωt + iβx1).

(6.180)

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

ругих волн и зависимость от x1 будет определяться exp(iβx1). Если это

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

361

ординаты x3 ничего определенного сказать нельзя, кроме одного: если

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

Предлагаем читателю подставить решение (6.180) в уравнение (6.179) и получить такие выражения для компонент векторов смеще- ний ue и ut :

u(n )

= A(n )

exp i (βx

− η

x

3

)

,

(6.181)

e

e

 

1

e

 

 

 

 

u(n )

= A(n )

exp i (βx

− η x

3

)

,

(6.182)

t

t

 

1

t

 

 

 

 

где n =1,3 , ηe = ke2 − β2 , ηt =

kt2 − β2 ,

ke

= ω/ce , kt = ω/ct ,

Ae(n ) и

At(n ) - неизвестные пока комплексные амплитуды.

При условиях β ≤ ke и β ≤ kt эти выражения определяют продоль-

ную (P-волну) (6.181) и поперечную (SV-волну) (6.182) однородные пло- ские волны. Поскольку волновые векторы ke и kt ориентированы в

направлении распространения соответствующей волны (см. рис. 6.21), то угол между этим направлением и осью Ox3 равен

θ = arcsin

β

,

θ = arcsin

β

.

(6.183)

 

 

e

ke

 

t

kt

 

 

 

 

 

Рис. 6.21. Иллюстрация к определению направлений распространения продольной и поперечной волн

Учитывая, что β = 2π/L , ke = 2π/λe , kt = 2π/λt , где λe и λt - длины продольной и поперечной волн, определим углы θe и θt через соотно-

шение пространственного периода возбуждения L с длинами волн λe и λt :

362

θ = arcsin λe

,

θ = arcsin λt .

(6.184)

e

L

 

t

L

 

 

 

 

 

Формулы (6.184) можно получить согласно наглядным геометриче- ским построениям. На рис. 6.22 синусоида отображает график про-

странственного возмущения σ33 (x1,t ) в некоторый момент времени t .

Параллельные прямые линии представляют собой фронты продольной волны, которые находятся на расстоянии λe один от другого, а прямая

ab определяет луч, вдоль которого распространяется продольная вол- на. Как видно, справедливость первой формулы (6.184) становится очевидной. Аналогичные построения можно провести и для попереч- ной волны.

Рис. 6.22. Геометрическая картинка соотношения между пространственным периодом возбуждения L с длиной продольной волны λe

Говоря образно, к пространственной гармонике возбуждения до- страиваютсядве волны: продольная и поперечная, которые распро- страняются под разными углами. Формулы (6.184) можно переписать еще в таком виде:

θ = arcsin

2πce

,

θ = arcsin

2πct

.

(6.185)

 

 

e

ωL

 

t

ωL

 

 

 

 

 

Как видим, при постоянном пространственном периоде возбуждения L с изменением частоты ω происходит поворот лучей продольной и поперечной волн (при уменьшении частоты ω углы θe и θt увеличи-

ваются). Тот же эффект имеет место если зафиксировать частоту и изменять пространственный период L . Если L и ω таковы, что λe = L

(или λt = L ), то луч, который определяет направление движения волны,

располагается параллельно поверхности. Если же при уменьшении L или ω получаем λe > L , то имеем β > ke . Аналогично при λt > L имеем

β > kt . При таких обстоятельствах, как мы уже знаем, образовываются неоднородные продольная и поперечная волны.

363

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

векторы смещения в этих волнах в виде

 

)

 

 

 

 

 

 

 

 

 

 

u

= e u(1) + e u(3) =

(

e A(1) + e

A(3)

exp i (βx

− η

x

3

)

,

(6.186)

e

1

e

3

e

1 e

 

3

e

 

1

e

 

 

 

 

 

u

= e u

(1)

+ e u

(3) =

(

e A(1)

+ e A(3)

)

exp i (βx

− η x

3

)

,

(6.187)

t

1 t

3 t

1 t

3

t

 

1

t

 

 

 

 

 

где (следите за проекциями маленьких стрелочек на оси координат

Ox1,Ox3 на рис. 6.21)

 

A(1)

= sin θ A

=

 

β

A ,

 

 

(6.188)

 

 

 

 

 

 

 

e

 

 

e

e

 

 

 

 

 

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

ke

 

 

 

A(3) = −cos θ

A

= −

ηe

A

 

= −

ke2 − β2

A ,

(6.189)

 

 

 

e

e

 

e

ke

e

 

 

 

 

ke

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A(1)

= cos θ

A

= ηt

A

 

=

 

kt2 − β2

A ,

(6.190)

 

 

 

 

t

 

t

t

kt

 

t

 

 

 

 

 

kt

t

 

 

 

 

 

 

 

 

 

 

 

β

 

 

 

 

 

 

A(3) = sinθ

A =

 

A .

 

 

(6.191)

 

 

 

 

 

 

t

 

 

 

t

 

t

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

kt

 

 

 

Здесь Ae и At еще неизвестные амплитудные множители.

Понятно,

что формулы (6.186) и (6.187) одинаково справедливы как для описа- ния однородных, так и неоднородных волн. В последнем случае вели-

чины ηe

и ηt

будут чисто мнимыми.

 

 

 

 

 

 

 

 

 

Амплитудные множители Ae и At

найдем из граничных условий

(6.178), а именно

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u3 +

u1

 

 

 

x3 = 0 ,

 

 

 

 

 

 

 

 

 

 

 

σ31

= μ

 

= 0

при

 

 

 

(6.192)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x1

x3

 

 

 

 

 

 

 

 

 

 

 

σ

33

= (λ + 2μ)u3 + λ

u1 = exp(iβx )

при x

3

= 0 .

(6.193)

 

 

 

 

 

x3

x1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Третье условие σ32 = 0 при x3 = 0

выполняется автоматически. Под-

ставим

в условия

(6.192),

(6.193)

выражения

u

 

= u(1)

+ u(1)

и

 

 

= u(3) + u(3)

, где u(1),u(1) и u(3),u(3)

 

 

1

e

t

 

u

3

определяются согласно форму-

 

e

t

 

 

 

e

t

e

t

 

 

 

 

 

 

 

лам (6.186)-(6.191). После несложных преобразований, которые реко- мендуем читателю провести самостоятельно, получим систему из двух уравнений относительно искомых амплитудных множителей Ae и At :

A

+

kt2 2β2

 

ke

A

= 0,

(6.194)

 

 

e

 

 

 

 

t

 

 

 

 

2βηe kt

 

364

 

 

 

 

 

 

 

A

 

2βηt

 

 

ke

A

= −

i

 

 

 

ke

.

(6.195)

 

 

 

 

 

 

μ k2

2β2

e

k2

2β2 kt

 

 

t

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

Решая ее, находим

 

 

 

 

 

 

 

 

 

 

 

kt2 2β2

 

 

 

 

 

 

 

A

= −i

ke

 

,

 

 

(6.196)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

 

 

 

 

μ

 

 

 

 

D (β)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

= i

kt

 

2βηe

,

 

 

 

(6.197)

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

μ D (β)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где определитель

D (β) = (kt2 2β2 )2 + 4β2ηeηt .

 

 

 

 

(6.198)

Перепишем формулы (6.196), (6.197) в виде удобном для проведе- ния численных расчетов. Для этого разделим числитель и знаменатель

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

β

1

− β

2

,

указанных формул на kt

и введем обозначение: β =

kt

, ηe =

q2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ηt′ = 1− β′2 , Ae′ = μkt Ae ,

At

= μkt At , где, согласно формуле (6.140), ве-

личина q =

ce

=

kt

= 2

 

1− ν

 

 

, ν

- коэффициент Пуассона. Тогда фор-

 

 

1

 

 

 

 

 

c

t

k

2ν

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

мулы (6.196), (6.197) примут вид

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ae

 

= −

 

 

 

 

 

12β

 

 

 

 

 

,

 

 

(6.199)

 

 

 

 

 

 

q

 

 

 

2 2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(12β

)

+ 4β ηe

ηt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

At

= i

 

 

 

 

 

2β ηe

 

 

 

 

.

 

 

 

(6.200)

 

 

 

 

 

 

 

 

 

 

2 2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

′ ′

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(12β

) +

4β ηeηt

 

 

 

 

 

 

 

 

 

Здесь может быть полезной такая аналогия. Если рассматривать пространственную гармонику возбуждения в качестве входного сигна- ла, саму среду как линейную систему, продольную и поперечную вол-

ны как выходной эффект, то зависимости Ae (β) и At (β) сыграют роль частотной характеристики данной системы.

365

Рис. 6.23. Зависимости модулей амплитудных коэффициентов Ae(β′)

(кривая 1) и At(β′) (кривая 2) от нормированной пространственной частоты

β′ ; ν = 0,3

На рис. 6.23 представлены зависимости модулей амплитудных ко-

эффициентов

 

 

A

(β′)

 

и

 

A(β′)

 

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

 

 

 

 

 

 

 

 

e

 

 

 

 

t

 

 

 

β

 

 

 

λ

 

 

 

 

 

 

частоты β′ =

 

 

=

Lt ,

 

коэффициент Пуассона ν = 0,3 . Анализируя эти

k

 

 

 

t

 

 

 

 

 

 

 

 

 

 

зависимости, нужно обратить особое внимание на возможность равен- ства нулю определителя D (β) . Как мы уже знаем (см. формулы (6.165)-

(6.167)), определитель D (β) имеет единственный действительный ко- рень β = kR . При β = kR амплитудные коэффициенты Ae (β) и At (β)

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

ет волновому процессу, который может существовать в твердом полу- пространстве со свободной границей без внешнего влияния. Этот вол- новой процесс называется волна Релея, которую мы исследовали в па- раграфе 6.11.

Отметим, что, согласно формулам (6.196), (6.197), отношение ам- плитудных множителей At /Ae при β = kR совпадает с формулой

(6.169), которая определяет связь между амплитудами неоднородных волн продольного и поперечного типа в волне Релея. Конечно, так и должно быть (убедитесь в этом самостоятельно).

Глядя на рис. 6.23 отметим также следующее:

1) если β = 0 , т.е. пространственный период L → ∞ , то имеем по-

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

366

увеличении пространственной частоты (β > 0 ), т.е. уменьшении вели-

чины L , в среде присутствуют две волны продольная и поперечная; 2) продольная волна является однородной в диапазоне изменения пространственной частоты 0 ≤ β < ke = kt /q , а поперечная волна есть

однородная при 0 ≤ β < kt . Для всех реальных материалов величина q >1, так при ν = 0,3 имеем q 1,87 . Тогда при β′ >1/q 0,535 (см. рис. 6.23) продольная волна будет неоднородной, при β′ >1 неоднород-

ной будет и поперечная волна. При дальнейшем увеличении величины β′ амплитуды возбуждения обеих неоднородных волн остаются посто-

янными;

0,707 в среде рас-

3) согласно (6.199) и рис. 6.23, при β = 0,5

пространяется только однородная поперечная волна.

6.13. Методы возбуждения объемных и поверхностных волн

Сформулируем перед собой такую задачу: рассмотреть мо- дель источника, которая описывает практическую возможность воз- буждения продольных, поперечных и поверхностных волн (волн Релея). Ключевым моментом для решения этой задачи является вывод, полу- ченный согласно результатам предыдущего параграфа: для того чтобы создать в полупространстве ту или иную волну, нужно возбудить по- верхность полупространства пространственной гармоникой с подоб- ранной должным образом пространственной частотой β (β = 2π/L , L -

пространственный период). Чтобы получить в полупространстве про- дольную волну под углом θe к оси Ox3 , необходимо выбрать период

пространственной гармоники L = λe /sinθe , где λe - длина продольной

волны (см. рис. 6.22). Одновременно с продольной волной будет рас- пространяться и поперечная волна под углом θt , таким что

sin θt = λt /L . Поскольку λt < λe , угол θt < θe . Если пространственный период λt < L < λe , то в среде однородной остается только поперечная волна. Когда величина L < λt в полупространстве имеем две неодно-

родных волны.

Возможно ли, получить только одну однородную продольную волну, которая распространяется под некоторым углом θe ? Теоретически -

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

367

Таким образом, используя тот или другой излучатель, необходимо стремиться создать пространственный период возбуждения с пра- вильно подобранным пространственным периодом L . Поскольку лю- бой излучатель имеет конечные размеры, естественно, речь может ид- ти о периодичности возбуждения только на участке, занятом излучате- лем. Влияние конечных размеров источника рассмотрим в следующем параграфе, а пока, рассматривая методы возбуждения, будем игнори- ровать это обстоятельство.

Существует несколько вариантов технической реализации перио- дического возмущения поверхности полупространства. Рассмотрим три метода (которые возможно читатель уже изобрел”).

Рис. 6.24. Возбуждение волн в твердом полупространстве:

а- метод клина, 1 -пьезопластина, 2 - клинообразная накладка;

биммерсионный метод, 1 - пьезопластина, 2 - жидкость;

в- метод гребенки, 1 - пьезопластина, 2 – гребенчатая накладка

Метод клина (рис. 6.24, а). В методе клина используется пьезопла- стина 1, которая совершает толщинные колебания и возбуждает тем самым в клинообразной накладке 2 продольную волну. Эта волна соз- дает на поверхности полупространства возбуждения с пространствен-

368

ным периодом L = λe(кл)/sin θкл , где λe(кл) - длина продольной волны в материале клинообразной накладки, θкл - угол клина.

Иммерсионный метод (рис. 6.24, б) отличается от метода клина тем, что продольная волна излучается в жидкую среду 2, которая кон- тактирует с полупространством. Этот метод удобен в случаях, когда деталь двигается относительно источника (например, лист проката).

Метод гребенки (рис. 6.24, в). Толщинные колебания пьезопластины 1 передаются в упругую среду через гребенчатую накладку 2, расстоя- ние между зубцами которой равняется L . В двух предыдущих методах можно без колебаний указать, в какую сторону от источника распро- страняется волна в твердой среде вправо. При использовании метода гребенки волны распространяются как влево, так и вправо. Это обу- словленное тем, что пространственное возбуждение в твердой среде

имеет характер стоящей волны cos (βx1)exp(iωt ). Его можно предста- вить в виде суммы бегущих навстречупространственных гармоник: exp(iβx1 iωt ) + exp(iβx1 iωt ) , каждая из которых возбуждает волну

соответствующего направления.

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

6.14. Волны в твердом полупространстве, возбуждаемые источником конечных размеров

Вызывает практический интерес вопрос, какие волны воз- буждает в полупространстве излучатель конченых размеров, который создает на конченом участке поверхности твердого полупространства некоторое периодическое пространственное возмущение. Рассмотрим такую задачу. Пусть компоненты тензора напряжения на поверхности твердого полупространства определены так:

σ

33

(x ,t )

= f (x

)exp

i (ωt − β

x

)

,

 

1

1

 

 

0 1

 

 

σ31

= σ32 = 0,

при

x3 = 0,

 

 

 

где f (x1) - некоторая функция. В частном случае

 

1,

x

[a /2, a /2],

f (x

 

1

 

) =

x

[a /2, a /2].

1

0,

 

 

1

 

(6.201)

(6.202)

369

Представим σ33 как функцию x1 в виде интегралу Фурье:

 

σ33 (x1,t ) = exp(iωt ) S (β)exp(iβx1)dβ.

(6.203)

−∞

 

Такая форма записи разрешает воспользоваться результатами па- раграфу 6.12. Действительно, выражение (6.203) задает напряжение на поверхности полупространства в виде бесконечного количества пространственных гармоник с пространственными частотами β. Каж-

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

и неоднородными. В соответствии с формулами (6.186), (6.187) их ам- плитуды обозначены Ae(n ) и At(n ) , где n =1,3 . Суммируя все эти волны, получаем звуковое поле, которое отвечает граничному условию (6.201):

u = u1 + u3 = e1u1 + e3u3 = e1 (ue(1) + ut(1) )+ e3 (ue(3) + ut(3) ), (6.204)

где

u

= u(n ) + u(n )

=

A(n )

(

β

S

(

β

)

exp i

(

βx

n

e

t

 

e

)

 

 

 

1

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

(n )

( )

S

(

β

)

 

(

1

 

 

 

t

 

 

 

 

A

 

 

β

 

 

exp i

 

βx

−∞

− η

x

3

) dβ +

e

 

 

 

 

 

 

 

 

(6.205)

− η x

3

) dβ,

t

 

 

 

здесь n =1,3 ; u1,u3 - проекции вектора смещения частичек среды;

ηe = ke2 − β2 , ηt = kt2 − β2 , ke = ω/ce , kt = ω/ct ; временной множи- тель exp(iωt ) в дальнейшем не пишем.

Для того чтобы получить результат интерференции бесконечного количества плоских волн необходимо вычислить интегралы (6.205). Математика дает возможность провести их корректное приближенное вычисление. Но мы воспользуемся более простым приемом на основе

физических соображений. Обратимся к функции S (β), которая пред- ставляет собой преобразование Фурье от σ33 по координате x1 . Из

(6.201) и (6.203) имеем

 

(

 

)

 

1

 

(

 

1)

 

 

 

(

0 ) 1

 

1

 

 

 

 

 

 

 

 

 

 

S

 

β

 

=

2π −∞

f

 

x

exp

i

 

β − β

x

dx .

(6.206)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В частном случае, когда f (x1)

определена формулой (6.202) получаем

 

 

 

 

 

 

 

 

 

 

) = a

sin (β − β

)a /2

 

S (β) F (β − β

0

 

 

0

 

 

 

.

(6.207)

 

(β − β0 )a /2

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

370