Диполь в твердом теле может быть двух видов: это может быть сфера, вложенная в полость того же радиуса так, что среда может свободно скользить по поверхности сферы, и это может быть сфера, «вмороженная» в среду и увлекающая за собой среду. Хотя граничное условие в первом случае похоже на граничное условие для жидкой среды, поля в жидкости и в твердом теле будут совершенно различны: помимо потенциального поля смеще ний продольной волны (похожего по форме на поле в жидкости), в твердом теле в обоих случаях появится еще и вихревое поле сдвиговых волн, для которого никакой аналогии в жидкой среде нет.
Дело в том, что для диполя в жидкой среде на границе сферы имеется только одно граничное условие — совпадение нормальных
смещений (или |
скоростей) сферы и среды: иг = и cos 0 (угол |
отсчитывается |
от оси |
осцилляций). Это же требование остается |
и для осциллирующей |
сферы в твердом теле, но к нему добав |
ляется еще условие обращения в нуль касательного напряжения на поверхности сферы г — а: огѳ = 0 (первый случай), либо условие равенства касательных смещений щ — — и sin Ѳ (второй случай). Поэтому, взяв скалярный потенциал в таком же виде, как и для жидкости,
Ф — М - 1 ё ЪіГcos Ѳ, ( 1 52.1)
мы получим желательное угловое распределение (по косинусу) нормальных скоростей частиц, но второе условие окажется не удовлетворенным. Поэтому в среде возникает еще и вихревое поле. Его можно описать векторным потенциалом, который давал бы такое же угловое распределение для нормальных скоростей. Это — вектор с единственной не равной нулю компонентой, направ ленной по параллели и равной
ф = AM ik‘4r ~_ 1 е‘V sin Ѳ• /ф, |
(152.2) |
где /<р — орт параллели.
Теперь можно распорядиться двумя остающимися пока неопре деленными величинами М и А так, чтобы удовлетворить обоим граничным условиям.
Начало расчета одинаково для обоих видов второго граничного условия. Смещения частиц будут происходить в меридиональных плоскостях с полярной осью, совпадающей с направлением осцил ляций сферы. Компоненты смещений выражаются через скалярный потенциал ф и не равную нулю компоненту ф векторного потен циала следующим образом:
|
_ дф |
+ |
дф |
|
|
|
дг |
дѲ + — ф ctg Ѳ, |
|
0 |
= _ L i2 ._ ilË ._ _ L ,h |
(152.3) |
г |
дѲ |
dr |
г |
Нормальные и касательные напряжения в меридиональной пло скости равны в данном случае агг = Киаа -(- 2цигг; or0 = 2\uirQ. Компоненты тензора деформации равны
игг |
д и г |
|
1 / д и а |
1 |
, |
1 д и г \ |
|
* |
дг > |
Ura — - г г * ( |
— |
---------------------- ~ иѳ Л |
T ~ w ) ' |
м а а |
— Д ф - |
|
|
гѲ~ |
2 V |
дг |
|
|
|
|
|
Подставляя |
сюда |
|
% |
и пользуясь уравнениями движения: |
(152.3) |
|
|
|
|
|
|
а2<р |
7T C tg0-||- = |
— fe?<p, |
|
|
|
|
|
г2 аѳ2 + |
Дф = |
А (0, 0, ф) = /ф |
а2ф I _2_ |
|
а2Ф |
|
|
Ö/-2 |
' |
г дг |
та- |
аѳ2 1 |
|
|
|
|
|
+ -jsr ctg Ѳ-Ц- — -і- (cosec2 Ѳ)ф] = |
— /г?ф/ф. |
приведем выражения для компонент тензора напряжений к виду
г . 2 . |
4 аф I |
2 а / . а аф \ |
агг — — 11 [^ ф + |
Т ~дг + |
г2sin ѳ Ж |
\sin ѳ аѳ ) “ |
|
|
|
(152.4)- |
°гѲ= - |
- 2ч { ч ^ ) |
+ 2ч [ - т ч ("М]}' |
Подставляя значения ф и ф из (152.1) и (152.2), можно переписать, формулы (152.3), (152.4) в виде
Ur = ~ tL cose [ ч Г + А ~Т~ R‘) ’ |
|
“e = |
- - H - s,n 0 [ - r ^ |
+ |
i4- r 4 - ( r ^ ) ] * |
= |
- |
ѵ |
~ |
ш |
cosѳ [ $ Rl + |
4 М |
ч |
Ri ) - |
|
|
|
|
|
|
|
|
|
(152.5); |
|
|
|
|
|
- ^ |
а 4 ? { ч |
М . г |
^ ѳ= |
- |
^ |
- |
^ |
5іпѲ{2 ^ ( - Г ^ ) |
+ |
|
|
+ А [А/Я/+ 217 ( " Г ^ г № > )]}* |
где введены обозначения |
|
|
|
|
J . |
|
|
|
|
|
/?, - |
|
|
|
еік‘г и |
Я, - |
‘ѴгГ -' е'*'г. . |
Дифференцирование дает:
dR t __ 2 — 2ik t r — { k i r f j i t , г |
|
d r |
~ |
Г3 |
е |
’ |
d (rRi) _ |
1— iktr —(k/r)“ j k , r |
|
dr |
~ |
r2 |
’ |
|
|
= |
|
e . |
V , |
_d_ |
— 3 + 3tk ,r + |
2 (V)'- — i (fefr)3 /i r |
rfr
Соответственные выражения производных для і?; получатся просто заменой kt на k[ в этих формулах. Подставляя найденные выражения в (152.5), найдем:
|
|
|
|
|
|
Щ = |
|
cosѲfіI2 - 2jV |
- (Ѵ )2К 'Ѵ + |
|
|
|
+ |
2А (iktr— 1)eiktr), |
|
иѳ = |
~ д а sin Ѳ (ОѴ — 1) еік‘г |
|
|
|
|
+ Л [ 1— і£,г — (ktrf] еік‘г), |
|
° г г = |
— Ц |
4^7 C0SѲ fKV)2(»V—1) + |
(152.6) |
|
+ |
4 [3 — 3tV — (Ѵ )211 |
_ |
|
|
|
— 4Л [3 — ЗіѴ — (V )2] е'Ѵ), |
|
<Дѳ = —Ц 4 ^ г sin Ѳ {2 [3 — 3/£,r — {ktr)2\ eik‘r—
— Л [6 — 6iktr — 3 (ktrf -j- i (ktrf] e‘k‘r).
Характеристики направленности для ur и |
arr — обычные |
дипольные характеристики. Направленность для |
и агѳ дается |
такой же характеристикой, как и для крутоля. Обе характеристики в пространстве можно получить, вращая вокруг оси х систему четырех одинаковых окружностей, касающихся осей х и у в начале координат. При этом ось х соответствует направлению осцилляций.
Теперь найдем выражение для силы F, с которой сфера дей ствует на среду. Эту силу можно выразить следующим интегралом по всей поверхности сферы:
F = — J (оггcos Ѳ— агg sin Ѳ) dS.
Ввиду осевой симметрии всей картины можем принять за эле мент поверхности идущую по параллели полоску ширины а d0; тогда dS = 2іш2 sin 0 dQ и интеграл примет вид
Л
F — — 2яа2 j (агг cos Ѳsin 0 — агѲsin2 Q)dQ.
о
Подставляя сюда значения orr и агд для г — а из (152.6), получим после простых преобразований
F = - v ^ ß 2l(1~ ia)eCa + 2А 0 - #)*'*], |
(152.7) |
где введены обозначения: а — k[ü\ ß = kta.
Случаи скользящей и вмороженной сферы рассмотрим отдельно. Для скользящей сферы одним из граничных условий является равенство нулю касательного напряжения на ее поверхности.
Полагая в последней формуле (152.6) о>е = 0, |
найдем А: |
А = |
6 |
— біа— 2а2 |
(152.8) |
6 — |
6iß—3ß2+ iß3еі (а—ß) |
Полученные выше формулы точные; они справедливы для лю бого размера сферы (или для любой частоты). Теперь дадим при ближенные величины для сферы, малой по сравнению с длиной сдвиговой волны (и подавно по сравнению с длиной продольной волны в среде): будем считать малым ß (и подавно а) и будем пре небрегать величинами а 2 и ß2 по сравнению с единицей. Подставляя (152.8) в (152.7), получим при этом условии
F = - l i - ^ = -p(o*M . |
(152.9) |
Далее, полагая г = а в (152.6), выразим М через амплитуду колебаний сферы и (пренебрегая, как и в предыдущей формуле, малыми второго порядка по а и ß). Для этого подставим в первое
.уравнение (152.6) величину (152.8); вначале придется сохранять все члены в числителе и в знаменателе, потому что старшие члены сокращаются. Произведя вычисления, найдем, пользуясь соотно шением иг = и cos 0:
м _ _ |
12ла3и |
/ 1 _ |
• 2ß3+ а3\ |
|
3ß2+ a2 |
\ |
1 3ß2+ а2 у ’ |
Теперь и силу, действующую со стороны сферы на среду, сможем выразить через смещение сферы:
F = |
рсо2 |
12яа3и |
. 2ß3+ а3 \ |
3ß2+ а2 |
1 3ß2+ aa ) * |
Если плотность сферы равна плотности среды р, то к F следует еще добавить инерциальную силу, соответствующую массе сферы,
равную — рсо2 -|- па3щ однако этой величиной можно прене
бречь по сравнению с упругой реакцией среды. Таким образом,. F можно считать силой диполя в твердом теле. Очевидно,
|
рсо2 |
12яа3I и I |
\ F\ |
3ß2+ а2 ' |
Отсюда видно, что коэффициент упругости для реакции среды равен
_ |F| |
12рсо2яа3 , 0 |
|х(А.4-2ц) |
_ |
12яар. |
Х —7 “Т ~ |
3ß2+ а2 — |
Яй ЗД,+ 7ц |
— з + [(.і/(А,-1-2р)] ' |
Здесь учтена только реактивная часть силы упругости, соответ ствующая данному смещению сферы.
Полученных данных достаточно для того, чтобы найти энер гию, излучаемую осциллирующей сферой в виде продольных и поперечных волн в окружающую среду. В самом деле, мощность излучения равна половине произведения амплитуды скорости колебания сферы на компоненту силы диполя, находящейся в фазе с этой скоростью.-Скорость равна —/соц, а активная, компо нента силы есть (с той же точностью, что и выше)
2ß3+ а3
— і I F I 3ß2 + а2
Излучаемую мощность найдем в виде
J — |
1Im F 11 соц I = |
I F Г- (2ß3+ а3) |
|
|
|
24ярсоа3 |
|
откуда после простых преобразований получим |
|
J' = ^sr [2 (Lp ) 2Pc' + (xtfi^)2Pc'] - |
(152-10> |
Сравнивая с (104.2), найдем, что мощность излучения про дольных волн в жидкость и в твердую среду (с той же скоростью, продольных волн) одинакова при равных силах диполя, прило женных к среде со стороны. Если же сравнивать мощность излу чения в твердую среду и в жидкость при одинаковых смещениях сфер одинаковых (малых) радиусов, то результаты получатся разные: излучение продольной волны в твердой среде оказывается
в |
.. |
36 |
- |
- |
|
,-^г раз больше, |
чем в жидкости: при данной ско- |
[3 |
(й,а)2 + |
(ft/а )2]2 г |
|
рости движения требуемая сила в твердой среде гораздо больше,
чем в жидкости.
Аналогичный расчет можно произвести и для «вмороженной» сферы. При этом условия на поверхности сферы имеют вид и, = = ц cos Ѳ, UQ — — и sin Ѳ,' что дает, согласно первым двум фор мулам (152.6),
(2 — 2іа — а 2) еіа + 2А(iß — 1 )е® = (іа — l)ela-\-A(l — iß — ß2) e‘p,
откуда точное значение величины А получается в виде
А = |
3 |
— Зіа |
—а2 |
(а-Р) |
(152.11) |
3 |
—3iß |
— ß2е |
Переходя ’ к приближенным |
формулам для |
ß -С 1, найдем |
для силы, действующей со стороны сферы на среду, снова ту же формулу-(152.9). Подставляя во второе уравнение (152.6) величину (152.11), найдем и величину М\
|
|
|
|
|
|
|
|
|
М = — |
12яа3и |
/. |
. 2ß3+ а3\ |
|
|
2ß2+ а2 V1 |
1 2ß2+ ct2) • |
Следовательно, |
|
12ла3м |
|
2ß3+ а3\ |
|
|
F, = роз,2 |
|
’ |
|
2ß2+ а2 |
|
2ß2+ а2) |
что дает коэффициент |
упругости, |
равный |
|
|
|F| _ |
12рш2яа3 _ |
, 0 |
ц(Я.+ 2(х) |
_ |
12яац |
1«| |
2ß2+ a2 |
п а |
2Я,+.5р |
2 + [|і/(Я,+ 2(і)]- |
Наконец, находя тем же способом, что и выше, излучаемую мощность, обнаружим, что она выразится той же формулой (152.10). Но смещения сферы оказываются теперь другими: отно шение смещений для скользящей и «вмороженной» сфер больше
единицы и равно |
|
|
|
< |
3ß2 + |
a 2 |
|
ЗА, + |
7|х |
2ß2 + |
a2 |
~ |
2Х + |
5Ң. ‘ |
Первый и второй члены в скобках выражения для излучаемой мощности (152.10) соответствуют излученным сдвиговым и' про дольным волнам соответственно. Распределение мощности между этими типами волн определяется отношением / СДв^прод = = 2 (ct/c() независимо от граничных условий на сфере. Большая часть энергии всегда идет в сдвиговые волны. В водоподобной среде (et/С[ -С 1) отношение смещений для двух типов граничных условий равно приблизительно 3/2. Эффективные коэффициенты упругости для водоподобной среды оказываются для двух слу чаев равными 4лл|.і и бяац.
Михаи л Александрович Исакович
ОБЩ А Я АКУСТИКА
М ., 1973 г ., 496 стр. с нлл.
Редакторы В. А . Григорова. Н. А . Райская
Техн. редактор С. Я ■ Шкляр
Корректор Н . Б. Румянцева
Сдано в набор 22/Ш 1973 г. Под. к печати 31/Х 1973 г. Бумага 60х90*/ц . Тип. № 2. Фнз. печ. л . 31. Уел. печ. л . 31. У ч .-нзд . л . 32,24.Тираж 17000 экз. Т-16939. Цена книги 1 р. 23 к . Заказ № 140
Издательство «Наука» Главная редакция физико-математической
литературы 117071, Москва, В-71, Ленинский проспект, 15
Ленинградская типография ЛГ« 6 Союзполнграфпрома Іірн Государственном комитете Совета Министров СССР по делам издательств, полиграфии и книжной торговли 193144, г. Ленинград, ул. Моисеенко, 10