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

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

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

Временной

множитель

exp(–iωt) писать не будем. Учитывая, что

kd << 1,

записываем

следующее

приближенное равенство:

exp(ikd cos θ) 1 ikd cos θ . Тогда

 

 

p = iωρ

V0d

 

ikr 1

cos θexp(ikr ).

 

4π r (r d cos θ)

 

 

 

Считая d/r << 1, получим следующее выражение для поля давления диполя:

p = iωρM0

ikr 1

cos θexp(ikr ),

(7.83)

 

 

4πr 2

 

где M0 = V0d. Величину M0 называют моментом диполя. Итак, поле

диполя определяется не объемной скоростью составляющих его моно- полей V0 и не расстоянием между монополями d отдельно, а произве- дением этих величин V0d. Одинаковое дипольное излучение можно получить с помощью разных пар противофазных монополей, подоб- рав объемные скорости монополей и расстояние между ними таким образом, чтобы моменты были одинаковыми.

Если расстояние d устремить к нулю, оставляя неизменными объ- емные скорости составляющих монополей, то момент диполя, а вме- сте с тем и поле диполя, также будет стремиться к нулю. Но если, уменьшая величину d, одновременно увеличивать объемные скоро- сти, так, чтобы момент диполя был неизменным M0 = const, то неизменным будет и поле диполя. Итак, приходим к понятию

точечного диполя.

Поле диполя можно также записать в виде:

p = iωρM0

exp(ikr )

 

 

 

 

.

(7.84)

 

 

 

z

4πr

 

 

Действительно, вследствие непосредственного вычисления по формуле (7.84) получаем

exp

(

ikr

)

 

=

exp

(

ikr

)

 

r

= ikr 1exp(ikr )

r

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

4πr

 

 

 

r

4πr

 

 

z

4πr 2

z

 

поскольку r/z = cosθ, то имеем результат, который совпадает с формулой (7.83). Соотношение (7.84) показывает, что можно ввести понятие диполя, продифференцировав поле монополя по координате источника z (рис. 7.17, а). Но первый вариант определения понятия диполя более приемлем с физической точки зрения.

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

451

ного излучателя является осциллирующая сфера. Интересным являет- ся случай осциллирующей сферы малого радиуса по сравнению с длиной излучающей волны (ka << 1). Рекомендуем читателю само- стоятельно осуществить преобразование формул (7.60), (7.77), пре- небрегая малыми величинами порядка (ka)2 по сравнению с едини- цей, и получить приближенное выражение для поля малой (ka << 1) осциллирующей сферы:

p = iωρ2πa3υ ikr 1cos θexp(ikr ).

(7.85)

0 4πr 2

 

Сравнивая выражения (7.83) и (7.85), можно сказать, что моментом малой осциллирующей сферы есть величина M0 = 2πa3υ0 . Поэтому,

если радиус сферы а и скорость υ0 изменяются так, что значение мо- мента М0 остается постоянным, то и поле излучения осциллирующей сферы будет неизменным.

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

Дальнее поле (kr >> 1) дипольного источника

p = −ωρk

M0

cos θexp(ikr )

(7.86)

 

 

4πr

 

имеет вид сферической волны с характеристикой направленности, которая определяется функцией косинуса (рис. 7.17, б).

Интересный результат можно получить при суперпозиции полей монополя и диполя. Итак, поместим в одну точку монополь с объем- ной скоростью V и диполь с моментом M. Совместная работа этих из- лучателей создает в дальней зоне (kr >> 1) поле

p = −iωρ 4Vπr exp(ikr )− ωρk 4Mπr cos θexp(ikr ) =

=

iωρ (V ikM cos θ)exp(ikr ).

(7.87)

 

4πr

 

Если выбрать V = ikM, то получим источник с характеристикой на- правленности в виде кардиоиды (рис. 7.18):

R(θ) =

1 + cos θ.

(7.88)

 

2

 

Максимум излучения имеем в направлении θ = 0 и нулевое излучение в противоположном направлении θ = π. Такой источник можно также получить с помощью одной сферы, которая одновременно выполняет

452

пульсации и осцилляции с соответственно подобранными амплитуда- ми.

Рис. 7.18. Характеристика направленности системы монополь-диполь

7.10. Коэффициент концентрации

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

Рассмотрим два излучателя ненаправленный и направленный. Пусть мощности излучения ненаправленного Р0 и направленного Рс излучателей одинаковы, т.е. Р0 = Рс. Ненаправленный излучатель на расстоянии r (в дальней зоне) создает звуковое поле интенсивностью I0, а интенсивность на том же расстоянии r вдоль некоторого направ- ления (θ,ψ) направленного излучателя Ic(θ,ψ). Коэффициентом кон-

центрации называют отношение

 

 

 

 

 

Ω(θ,ψ) =

Iс(θ,ψ)

 

 

(7.89)

 

 

.

I0

 

 

P0 =Pс

 

 

 

 

Пусть в направлении (θ0,ψ0) интенсивность направленного излучателя максимальная Iс(θ0,ψ0 ) = Iс max , тогда отношение

453

Ω(θ,ψ) =

Iс max

 

(7.90)

I0

 

 

P0 =Pс

 

 

называют коэффициентом осевой концентрации. Используя формулу

(7.65), записываем интенсивность направленного излучателя в виде

p 2 Iс = 2ρc .

Согласно определению характеристики направленности (7.33)) имеем p = p max Rp (θ,ψ), тогда (7.91) примет вид

Iс =

 

 

p

 

max

R2p (θ,ψ) = Iс maxR2p (θ,ψ).

 

 

 

 

 

 

 

 

 

 

 

 

 

2ρc

(7.91)

Rp(θ,ψ) (см.

(7.92)

Интенсивность ненаправленного источника можно вычислить по формуле

I0

=

P0

.

(7.93)

4πr 2

 

 

 

 

Подставив (7.91) и (7.93) в (7.89), получим такое выражение для ко- эффициента концентрации:

Ω(θ,ψ) =

2π

 

pr

 

2

.

(7.94)

 

 

 

 

 

 

 

ρcP

 

 

 

 

 

 

 

 

0

 

 

 

 

Связь между коэффициентом концентрации в заданном направлении и коэффициентом осевой концентрации устанавливаем из их опреде- лений:

Ω(θ,ψ) =

I

с

(θ,ψ)

 

 

=

Iс maxR2p

(θ,ψ)

 

= Ω0R2p

(θ,ψ). (7.95)

 

 

 

 

 

 

 

 

I0

 

 

I0

 

 

 

 

 

 

P0

=Pс

 

 

P0 =Pс

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отметим, что по определению Ω0 > 1, тогда как коэффициент Ω(θ,ψ)

может быть меньше единицы и даже равным нулю, если для данного направления характеристика направленности R(θ,ψ) = 0.

Выведем формулу для коэффициента осевой концентрации через R(θ,ψ) направленного излучателя. Исходным будет равенство мощно- стей Р0 = Рс. Мощность излучателя определим, интегрируя его интен- сивность по сфере большого радиуса r (в дальней зоне). Для нена- правленного излучателя имеем

P0 = 4πr2I0.

(7.96)

454

Для направленного источника излучаемая мощность примет вид:

P =

I

с

(r,θ,ψ)dS = I

с max

r 2

π 2πR2

(θ,ψ)sinθdθdψ,

(7.97)

с

 

 

 

∫ ∫

p

 

 

 

(0)

 

 

 

 

 

0 0

 

 

 

где элемент сферической поверхности dS = r2sinθdθdψ. В соответст- вии с равенством Р0 = Рс приравняем выражения (7.96) и (7.97):

4πr 2I0 = Iс maxr 2

π 2π

(θ,ψ)sin θdθdψ ,

 

R2p

 

 

 

0 0

 

 

 

откуда

 

 

 

 

 

Ω0 =

 

 

4π

 

.

(7.98)

2π π

 

 

 

 

∫ ∫ R2p

(θ)sinθdθdψ

 

0 0

 

 

 

 

 

Если характеристика направленности имеет осевую симметрию отно- сительно θ = 0 и не зависит от угла ψ, то

Ω0 =

2

.

(7.99)

πR2p (θ)sinθdθ

 

 

 

 

0

 

 

Вычислим Ω0, например, для излучателя, у которого Rp (θ) = cos θ . Со-

гласно (7.99) имеем Ω0 = 3.

Равнозначным определению (7.90) будет следующее определение: коэффициент осевой концентрации это отношение мощностей ненаправленного и направленного излучателей при условии создания ими одинаковой интенсивности звука на некотором расстоянии r:

Ω0

=

P0

 

.

(7.100)

Pс

 

 

 

I0 =Iс max

 

 

 

 

 

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

455

7.11. Взаимодействие источников звука

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

Запишем поле давления синфазной (+ +) и противофазной (+ –) пары монополей в дальней зоне (рис. 7.10, б), временной множитель exp(–iωt) опускаем:

 

++

 

 

V0

 

exp(ikr )± (iωρ)

V0

exp(ik(r d cos θ))

 

 

 

p+− = −iωρ

 

 

 

 

4πr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4πr

 

 

 

 

 

 

 

 

или, иначе (сделайте преобразование самостоятельно):

 

 

 

для синфазной пары

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

++

= −iωρ

V

0

 

 

 

 

(ikr )2cos

kd

 

 

 

 

 

 

 

i

kd

 

 

 

 

 

 

 

exp

 

 

cos θ

exp

 

cos θ

,

(7.101)

 

4πr

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

для противофазной пары

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+−

 

V

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kd

 

 

 

 

 

 

 

 

kd

 

 

 

p

 

= −iωρ

 

 

 

exp

(ikr )2i sin

 

cos θ

exp

i

 

cos θ .

(7.102)

 

4πr

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Интенсивность отдельного монополя

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I

(0)

=

 

p

(0)

 

2

 

=

 

(ωρ)2

 

V02

,

 

 

 

 

 

(7.103)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2ρc

 

 

 

 

2ρc

 

(4πr )2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

по аналогии для пары монополей имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

I ++ =

 

p++

 

2

 

 

 

I (0)4cos2

kd

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

cos

θ ,

 

 

 

 

 

 

 

 

2ρc

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

I

+−

 

 

 

p+−

 

2

 

 

 

I

(0)

 

 

 

2

kd

 

 

 

 

 

 

(7.104)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

2ρc

 

 

=

 

 

 

4sin

 

cos θ .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

Определим мощность излучения исследуемых источников звука; ин- тегрирование интенсивности по поверхности сферы большого радиу-

456

са, которая окружает излучатель, можно выполнить аналогично как в параграфе 7.7 (см. рис. 7.15). Как следствие, получим

P ++ = P (0)2

1 +

sin(kd)

 

,

P +− = P (0)2

1

sin(kd)

 

, (7.105)

 

 

 

 

kd

 

 

 

 

kd

 

 

 

 

 

 

 

 

 

 

где P (0) = I(0)4πr 2 мощность единичного монополя. Очевидно, выра- жения, стоящие в квадратных скобках (7.105), характеризуют взаи- модействие монополей через создаваемое ими звуковое поле. Тогда мощность одного монополя, который работает в паре, равна

P ++ = P ++

= P (0) 1 +

sin(kd)

 

,

P +− = P +−

= P (0) 1

sin(kd )

.

(7.106)

 

 

1

2

 

kd

 

 

1

2

 

kd

 

 

 

 

 

 

 

 

 

 

Рис. 7.19. Нормированные зависимости P1++ (1) и P1+− (2) от волнового рас- стояния kd

На рис. 7.19 приведены зависимости P1++ /P (0) и P1+− /P (0) от вол-

нового расстояния между монополями kd. Эти зависимости позволя- ют сформулировать следующие выводы.

1.Мощность отдельного монополя зависит от волнового расстоя- ния между монополями.

2.Мощность одного монополя может быть как больше, так и меньше, чем мощность, которую он излучает, работая в одиночку. Ес- ли kd = nπ , n =1,2,3,..., то взаимодействие по полю отсутствует.

3.Два близко расположенных (kd << 1) синфазных монополя излу-

чают мощность в четыре раза большую, чем один монополь с такой

же производительностью V0.

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

поля равна υr = υ0 exp(–iωt). Тогда мощность, излучаемая монополем

457

P(0) =

υ02

Re Zи =

υ02 S Re ζ

 

S

=

υ02 S Re

p

 

 

 

=

V0

Re p

 

S

(7.107)

 

 

 

 

 

 

υ0

 

 

 

 

2

 

2

 

 

 

2

 

S

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

равняется объемной скорости V0 = υ0S, умноженной на действитель- ную часть звукового давления на поверхности монополя. Пусть моно- поль работает, находясь в поле давления другого монополя. Тогда давление на его поверхности состоит из собственного и дополнитель- ного давления, которое создается на его поверхности другим излуча- телем. Полагая радиусы монополей малыми по сравнению с расстоя- нием между ними, можно считать, что дополнительное давление на поверхности сферы распределено равномерно, как и собственное давление излучателя.

Запишем комплексную амплитуду давления монополя в виде ряда, разложив функцию exp(ikr) в ряд Тейлора по степеням величины kr :

 

V

0

exp(ikr ) = −iωρ

V

0

 

(kr )2

 

(kr )3

p = −iωρ

 

 

1 + ikr

 

i

 

4πr

 

 

2

6

 

 

4πr

 

 

 

 

 

 

 

 

 

 

 

 

V

0

 

 

i

+k + i

k2r

 

k3r 2

 

= ωρ

 

 

 

 

 

+... .

 

 

r

2

6

 

4π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+... =

(7.108)

При малой величине kr действительная часть давления Re p определя- ется вторым слагаемым в скобках и равняется ωρkV0/(4π) (мнимая часть соответствует первому слагаемому). Видим, что Re p не зависит от расстояния, пока членами высших порядков можно пренебрегать, т.е. пока величина kr достаточно мала по сравнению с единицей. Та- ким образом, для каждого из двух синфазных монополей действи- тельная компонента давления, которое создается одним из них на по- верхности другого, равняется действительной компоненте давления, создаваемого самым излучателем. Вследствие этого действительная составляющая давления удваивается, а, следовательно, в соответст- вие (7.107) удваивается и излучаемая мощность, причем независимо от расстояния между излучателями, пока это расстояние остается малым по сравнению с длиной волны.

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

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

458

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

7.12. Излучение звука источниками, которые имеют неодинаковую скорость движения частиц поверхности излучателя

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

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

Метод, который будет ниже использован для исследования, приго- ден лишь для излучателей, которые имеют вполне определенную форму поверхности (например, сфера, цилиндр, эллипсоид). Это свя- зано не с физическими, а с математическими соображениями. Для того чтобы вычислить излучаемую мощность или характеристику на- правленности, нужно определить звуковое поле излучателя по задан- ной скорости его поверхности, или, другими словами, найти решение уравнения Гельмгольца р + k2p = 0 при заданном значении скорости на поверхности излучателя. Полученное решение, кроме того, должно представлять собой волну, которая уходит на бесконечность от излу- чателя (это так называемое условие излучения).

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

459

чении колебаний мембраны. Рассмотренные примеры показали, что волновое уравнение в частных производных превращается в обычное дифференциальное уравнение. Существует всего 11 систем ортого- нальных криволинейных координат, которые допускают разделение переменных в уравнении Гельмгольца [ , с. 73—118].

Однако хорошо изучены лишь некоторые случаи; рассмотрим ци- линдрические и сферические координаты и соответствующие излуча- тели в форме цилиндра и сферы. Эти случаи не только самые про- стые, но и наиболее близкие ко многим реальным излучателям: на практике встречаются и цилиндрические, и сферические излучатели. Кроме того, эти примеры дают возможность принципиально разо- браться в особенностях излучения звука при неравномерном распре- делении скорости колебаний поверхности излучателя.

7.12.1. Цилиндрический излучатель

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

Рис. 7.20. Пример цилиндрического излучателя

Пусть частицы цилиндрической поверхности излучателя колеблют- ся по гармоническому закону с частотой ω. Для описания поля такого излучателя целесообразно использовать цилиндрические координаты (r,z,ψ). Уравнение Гельмгольца для комплексной амплитуды давления p(r,z,ψ) в цилиндрических координатах имеет такой вид [31, 52]:

Арфкен Г. Математические методы в физике. — М.: Атомиздат, 1970. — 712 с.

460