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

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

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

На рис. 7.12 представлена диаграмма направленности при d/λ = 0,5, β = π/2. Сравнивая рис. 7.12 и рис. 7.11, а, замечаем, что вследствие сдвига фаз между колебаниями монополей возникло сме- щение максимумов и минимумов характеристики направленности, что сопровождается изменением формы лепестков.

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

В конце этого параграфе сделаем интересное наблюдение. Форму- ла (7.36) определяет давление звукового поля синфазной пары моно- полей при любом расстоянии r. Исследуем, как изменится это выра-

жение, если точка наблюдения

находится в дальней зоне

(рис. 7.10, б). Здесь полагаем r r1 =

r = dcosθ, a A1 = A = ωρV0 /(4πr).

Поэтому, в рамках принятых приближений давление в дальней зоне можно определить так:

p = −iωρ

V0

exp(iωt + ikr ) 1 + exp(ikd cos θ) .

(7.54)

 

 

4πr

 

 

 

 

 

 

 

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

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

441

фазовым распределением скорости колебаний на поверхности излу- чателя.

7.7. Осциллирующая сфера

В параграфе 7.2 исследовали звуковое поле, которое соз- дается благодаря пульсирующему движению жесткой безмассовой сферы радиусом a. Сейчас рассмотрим ситуацию, когда центр сферы двигается вдоль оси Oz со скоростью υ = υ0exp(–iωt), т.е. сфера колеб- лется как одно целое вокруг своего положения равновесия z = 0 (рис. 7.13). Такой источник звука называют осциллирующей сферой, или

сферическим источником первого порядка, а волну, которую он соз-

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

Рис. 7.13. Пример осциллирующей сферы

Чтобы определить звуковое поле осциллирующей сферы, нужно найти такое решение p(r, θ, ψ, t) = p(r, θ, ψ)exp(–iωt)волнового уравнения, которое бы удовлетворяло граничным условиям на по- верхности сферы и, конечно, условию излучения. При таком движе- нии источника скорость частиц сферы и акустической среды будет иметь две составляющих: радиальную υr вдоль радиуса r и танген- циальную υθ перпендикулярную к радиусу r (см. рис. 7.13). Гранич- ное условие на поверхности сферы состоит в равенстве радиальной составляющей скорости частиц сферы υr = υ0cosθ exp(–iωt) и радиальной компоненты скорости частиц в звуковом поле при r = a, т.е.

1 p

 

= υ0 cos θ.

(7.55)

 

 

 

iωρ ∂r

 

r =a

 

 

 

Обратим внимание, что на тангенциальную компоненту скорости частиц сферы υ0sinθ exp(–iωt) не накладывается условие равенства

442

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

Поскольку в поставленной задаче имеет место осевая симметрия, то зависимость параметров звукового поля от сферической коорди- наты ψ отсутствует, и потому третье слагаемое в уравнении Гельм-

гольца (7.4) в сферических координатах равно нулю. Тогда, это урав- нение можно записать так:

p

 

 

r 2

 

+

 

r

r

 

1

p

+k2r 2 p = 0.

 

 

 

 

sin θ

 

(7.56)

 

 

sin θ ∂θ

∂θ

 

 

Найти решение уравнения (7.56) — сложная задача, но ее можно упростить, если в решении попробовать угадать зависимость от угла θ. Предположим, что решение имеет вид p(r, θ) = R(r)cosθ, т.е. зависи- мость от угла θ в звуковом поле такая же, как и для радиальной со- ставляющей скорости на поверхности сферы. Подставив это решение в (7.56), получим

d

2 dR

2 2

 

 

 

r

 

2R(r ) +k r

R(r ) = 0.

(7.57)

 

dr

dr

 

 

 

Таким образом, угловая зависимость решения в виде cosθ является правильной, и физически это означает, что звуковое поле по коорди- нате θ повторяет зависимость радиальной составляющей скорости на поверхности сферы. Следует отметить, что это не тривиальный ре- зультат. Например, если решение записать в виде p(r, θ) = R(r)cos(nθ) при n > 1, то сократить на cos(nθ) при подставке такого решения в уравнение (7.56) не удастся. Это свидетельствует о том, что поле не повторяет колебания сферы в случае зависимости радиальной со- ставляющей скорости типа cos(nθ), n > 1. Так исчезла привлекатель- ная возможность разложить сложное колебание сферы на гармоники (убедитесь в этом самостоятельно).

Из выражения (7.57) получаем уравнение для функции R(r):

d

dR

+ (k2r 2

2)R = 0.

 

r 2

 

 

dr

dr

 

 

Если kr >> 1, то (7.58) можно переписать в виде

d

2 dR

2 2

 

r

 

+k r R = 0,

 

dr

dr

 

(7.58)

(7.59)

что совпадает с уравнением (7.6) для пульсирующей сферы. Понятно, что в этом случае решение (7.59) совпадает с решением для пульси-

рующей сферы (7.8): R (r ) = Ar exp(ikr ). Решение Br exp(ikr ) отбрасы-

443

ваем, поскольку оно не удовлетворяет условию излучения. Таким об- разом, на большом волновом расстоянии от осциллирующей сферы (kr >> 1) зависимость давления звукового поля от расстояния r такая же, как и для пульсирующей сферы. (Вспомните вывод, сделанный относительно формулы (7.54).)

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

R(r ) =

A 1 +

i

exp(ikr ).

 

 

r

 

 

 

 

 

kr

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

p (r,θ,t ) =

A 1 +

i

cos θexp(iωt + ikr ),

(7.60)

 

 

r

 

 

 

 

 

 

kr

 

где постоянная А определяется из граничного условия (7.55). Если для давления пульсирующей сферы нет деления на ближнее и дальнее по- ле, то в случае осциллирующей сферы такое деление есть. При kr >> 1 формула (7.60) приобретает вид

p (r,θ,t ) =

A exp(iωt + ikr )cos θ.

(7.61)

 

r

 

Как и в случае пары монополей (см. (7.54)), в дальней зоне поле ос- циллирующей сферы представляет собой сферическую волну, ампли- туда давления которой на сферическом фронте определяется функ- цией, которая не зависит от r, а только от угла θ (в данном случае это cosθ). Так, при θ = 0 давление будет максимальное, а при θ = 90° равняется нулю. На рис. 7.14 приведена характеристика направлен- ности осциллирующей сферы по давлению Rp(θ) = |cosθ |. Отсюда следует важный вывод: направленное излучение звука получено за счет специального распределения колебательной скорости на поверх- ности источника.

Рис. 7.14. Характеристика направленности осциллирующей сферы

444

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

υ =

1

 

p =

1

A

exp(iωt + ikr )

ik

2i

 

2

cos θ,

 

 

 

 

 

 

 

r

iωρ ∂r

 

iωρ

 

 

 

 

r

 

kr 2 r

 

 

 

 

 

 

 

 

 

 

υθ =

1

 

1 p =

1

A

exp(iωt + ikr )

1 +

 

i

sinθ.

 

 

 

 

 

iωρ

 

 

 

 

 

 

iωρ r ∂θ

 

 

r 2

 

 

 

kr

 

(7.62)

(7.63)

Если kr >> 1, то υr 1/r, υθ 1/r 2, т.е. на большом волновом расстоя- нии от источника частицы среды двигаются практически радиально

(υr >> υθ).

Оказывается (и об этом мы еще поговорим), что при произвольном распределении скорости на поверхности источника, на большом вол- новом расстоянии (kr >> 1) от источника, звуковое поле имеет вид

p (r,θ,ψ,t ) = A

exp(iωt + ikr )

q (θ,ψ),

(7.64)

r

 

 

 

т.е. представлено в виде двух множителей: поля сферической волны нулевого порядка и функции q (θ, ψ), которая определяет направлен- ность излучения (для осциллирующей сферы q (θ, ψ) = сosθ). Другими словами, в дальней зоне звуковое поле произвольного излучателя представляет собой волну, в которой зависимость амплитуды давле- ния от расстояния r такая же, как в сферической волне; при этом распределение амплитуды давления на сферической поверхности оп- ределяется функцией |q(θ, ψ)| и не зависит от расстояния r.

Понятно, что представление поля в виде (7.64) позволяет вычис- лить интенсивность в дальней зоне произвольного излучателя по формуле (7.29) для интенсивности в сферической волне:

Ir =

 

p

 

2

,

(7.65)

 

 

 

 

 

 

2ρc

 

 

 

где p амплитуда давления в дальнем поле источника звука.

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

мулой (7.16а). При этом необходимо вычислить интеграл ∫∫ pυ dS . По-

S

скольку звуковое поле обладает осевой симметрией относительно оси Oz, то в качестве элемента интегрирования можно взять поясок между близкими параллелями, приняв за переменную интегриро- вания угол θ (рис. 7.15). Площадь такого пояска равна 2πa2sinθdθ.

445

Итак, ∫∫ pυ dS = 2πa2

π

(pυr )

 

r =a sin θdθ. Подставив сюда выражения

 

 

S

0

 

 

 

 

 

 

(7.60) и (7.62) для давления p и колебательной скорости υr , после ин-

тегрирования получим

 

 

 

 

 

 

 

 

 

 

 

 

 

2 2

 

 

A

 

2

 

 

 

i

 

 

 

 

 

∫∫ pυ dS = 2πa

 

 

 

 

 

 

i 1

+

 

 

 

ik +

3

 

 

ωρ

 

S

 

 

 

 

ka

 

2i

2

 

(7.66)

 

 

 

ka2

a

 

 

 

Рис. 7.15. Пример элемента интегрирования

Импеданс осциллирующей сферы, согласно определению (7.16а),

 

Zи =

 

1

 

 

∫∫ pυ dS

(7.67)

 

(r = a,θ = 0)

 

2

 

υ

 

S

 

 

r

 

 

 

 

 

Учитывая (7.66) и (7.62) (при r = a, θ = 0), после несложных преобразо- ваний получим

 

4

πa2

 

(ka )

4

i (ka )

3

+ 2ka

 

ρc = ρc S

(Rи iXи), (7.68)

Zи =

 

 

 

 

 

 

 

 

 

 

3

 

4 + (ka )4

4 + (ka )4

 

3

 

где S = 4πa2 площадь поверхности сферы;

Rи =

 

(ka)4

; Xи = (ka)3 + 2ka .

(7.69)

 

+ (ka)4

4

4 + (ka )4

 

Как видим, Zи имеет комплексный характер, т.е. имеется синфазная

составляющая между р и υr, а также составляющая, для которой при- сущ сдвиг фаз на 90° между р и υr . Как и для пульсирующей сферы,

446

мнимая часть импеданса излучения осциллирующей сферы имеет ха- рактер массы, т.е. iIm[ Zи ] = iωMпр, где Мпр присоединенная масса:

 

 

 

Mпр =

4

πa

3 (ka )2

+ 2

ρ.

(7.70)

 

 

 

3

 

(ka )4

+ 4

 

 

 

 

 

 

 

 

 

 

Если волновой размер сферы мал, то

 

 

 

Mпp =

1

 

4

πa3

 

ρ

при условии ka << 1,

(7.71)

2

 

3

 

 

 

 

 

 

 

 

 

 

 

 

т.е. равняется массе среды в половине объема сферы.

Соотношение между реактивной и активной частями импеданса излучения зависит от соотношения между радиусом сферы и длиной звуковой волны, т.е. от величины ka. На рис. 7.16 представлены за- висимости Rи и Xи от ka. Для предельного случая ka >> 1 реактив-

ная часть импеданса излучения стремится к нулю, а активная к

ρcS/3.

Рис. 7.16. Зависимости коэффициентов Rи и Xи от ka для осциллирующей сферы

На практике важен противоположный предельный случай малый радиус сферы по сравнению с длиной волны (ka << 1). Поскольку ак- тивная составляющая Zи связана с распространением энергии в сре-

де, а реактивная составляющая Zи с локальным перераспределе-

нием энергии вблизи излучателя, то эффективность излучателя мож- но характеризовать отношением Re[Zи]Im[Zи] . При условии ka << 1

имеем:

для пульсирующей сферы

Re Zи(0)

(ka ) ,

(7.72)

(0)

 

 

Im Zи

 

 

для осциллирующей сферы

447

Re Zи(1)

(ka )3 .

(7.73)

(1)

Im Zи

 

 

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

7.8. Энергетические характеристики осциллирующей сферы

Поскольку среда не имеет потерь, то мощность, излучае- мую осциллирующей сферой, можно рассчитать как поток мощности, который проходит через замкнутую поверхность, окружающую ис- точник, например, сферу большого радиуса r, описанную вокруг ис- точника как из центра. При kr >> 1 давление и радиальная скорость стремятся к величинам (см. (7.61) и (7.62)):

p =

A cos θexp(iωt + ikr ),

υ

=

1

A cos θexp(iωt + ikr ). (7.74)

 

 

r

r

 

ρс r

Тогда интенсивность Ir = Re (pυr* )/2 будет определяться по формуле

Ir

=

 

 

A

 

2

cos2 θ.

(7.75)

 

 

 

 

 

 

 

2ρсr 2

 

 

 

 

Мощность находим, интегрируя интенсивность Ir по всей поверхности

 

 

 

 

 

 

 

 

 

 

π

 

 

сферы радиусом r (см. рис. 7.15):

 

P (1) = Ir 2πr 2 sin

θdθ.

С учетом

 

 

 

 

 

 

 

 

 

 

0

 

 

 

(7.75) находим мощность осциллирующей сферы:

 

 

P

(1)

=

2

π

 

A

 

2

.

 

 

(7.76)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

ρс

 

 

 

 

 

 

 

 

 

 

 

Постоянную А определим из граничного условия (7.55):

 

 

 

 

υ ka3 exp(ika )

 

 

A = −iωρ

 

0

 

 

 

 

 

 

 

.

 

(7.77)

 

 

 

 

 

 

 

 

2

 

 

 

2ka + 2i i (ka )

 

 

Итак, мощность осциллирующей сферы определяется по формуле

448

P (1) =

2

π

(ωρ)2

υ02k2a6

=

υ02

Re Zи(1),

(7.78)

3

ρс

4 + (ka )4

 

 

 

2

 

 

где Zи(1) сопротивление излучения осциллирующей сферы (см. (7.68)).

Как и следовало ожидать, мощность пропорциональна активной части сопротивления излучения.

Рассмотрим важную для практики ситуацию, когда волновой раз- мер источника мал, т.е. ka << 1 (в электроакустике все излучатели низких частот находятся в таких условиях). В этом случае (7.26) и (7.78) принимают такой вид:

для пульсирующей сферы

P (0) =

υ02

ρcS (ka )2 ,

(7.79)

 

 

2

 

 

 

 

для осциллирующей сферы

 

 

 

 

 

 

 

P (I) =

υ02

ρc S (ka )4

,

(7.80)

 

 

2

 

3

4

 

 

откуда

 

 

 

 

 

 

 

 

P (I)

=

(ka )2

 

 

 

 

 

 

12 .

 

(7.81)

 

P (0)

 

Как видим, при условии ka << 1 излучение осциллирующей сферы бу- дет менее эффективным, чем пульсирующей сферы. Этот факт мож- но объяснить тем, что при осцилляциях сферы разности давлений пе- ред и за малой сферой успевают выровняться местными паразитны- ми потоками” — перетеканием жидкости из мест с большим давлени- ем в места с малым давлением. Как следствие, эффективность излу- чения резко снижается вследствие значительного уменьшения актив- ной части сопротивления излучения и большого значения реактивной составляющей. Данная ситуация имеет образное название: “акусти- ческое короткое замыкание”. При пульсациях сферы одинаковое аку- стическое давление создается сразу на всей поверхности сферы жидкости перетекать некуда. При высоких частотах осцилляций сфе- ры, когда длина волн намного меньше размера источника, изменение давления происходит настолько быстро, что растяжения и сжатия не успевают взаимно компенсироваться.

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

449

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

7.9. Диполь

В параграфе 7.5 была построена математическая модель точечного источника (монополя), которая широко используется при анализе звуковых полей. Другая модель источника, которая также широко используется в теоретических исследованиях, представляет собой пару противофазных монополей (рис. 7.17, а), размещенных на расстоянии d, малом по сравнению с длиной волны. Такой источник называют диполем.

Рис. 7.17. Диполь (а) и его характеристика направленности (б)

Пара противофазных монополей создает поле давления

p = iωρ

V0

 

exp(ikr )iωρ

 

 

V0

 

 

 

exp(ikr ikd cos θ) =

 

4πr

4π(r d cos θ)

 

 

 

 

 

 

 

 

 

 

 

 

 

V0

 

 

exp

(

ikd cos θ

)

 

 

 

 

= iωρ

1

 

 

 

exp(ikr ).

(7.82)

 

 

 

 

 

 

 

 

 

 

 

 

4π r

 

 

r d cos θ

 

 

 

450