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

Аэродинамические источники шума

..pdf
Скачиваний:
58
Добавлен:
15.11.2022
Размер:
11.44 Mб
Скачать

цей, что в качестве

вместо

нулей

производной функции

Бессе­

ля следует понимать корни уравнения

 

 

 

 

 

J т{Со) ^ m(Q^o)

J пг (Q Со)^т(Со) — О

(4.62)

и вместо

(4.60)

для величины

Е необходимо

пользоваться

выра­

жением

 

 

да2

 

/

 

т 2

\

 

 

 

 

 

 

 

 

 

 

 

£= ■

‘ " “Со

"

Г

о

Г

Г

 

(4.63)

 

 

 

1

'

 

Q£o

'

 

 

 

 

 

+ (Х2/Q

 

 

 

 

где

 

 

t1= = j'm (C0)/^m(QC0)-

 

 

(4. 64)

При m = 0 для

низшей радиальной

моды

(0,1)

кольцевого

канала

формула

(4.63) упрощается

 

 

 

 

 

 

 

 

 

 

E = \ - Q .

 

 

 

(4.65)

Для плоского канала в формулах

(4.56)

(4.61) следует по­

ложить £ = 1 , £0 = яя/2а. Для

низшей

моды и в отсутствие

потока

имеем

 

 

С2=

— ika$.

 

 

(4.66)

 

 

 

 

 

Затухание на длине 2а будет

 

 

 

 

 

 

 

 

 

Д£0 =

8 , 6 8 1ш (V (2 to )2 + 4/Aap).

(4.67)

Предполагая, что выполняется также условие

|R ePK |A a — Imp|, для величины затухания из (4.67) будем иметь

4 A , - M 8 R . P |/ - = i = r

(4.68)

что совпадает с выражением' (4.49).

Возвратимся к цилиндрическому каналу и рассмотрим случай, когда частота звука намного превышает критическую, т. е. выпол­

нено условие £2 (^>*LI (1 — М2)в уравнении (4.26), что соответствует

Аа/(Со]/тГ=М'2) > 1 .

(4.69)

В этом случае формула

(4.45) значительно упрощается,

и выра­

жение для £2 в (4.54) примет вид

 

* 2

к2_____ \-ikCl

(4. 70)

 

 

~0 (1 + М)2£

Условие |pG| C l в данном случае означает

|Р|

1

(4.71)

(1 )2 ^

 

ч

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

М т п = 17,37

i m l / W - e L o - M 2)

(4.72)

 

1 —М2

 

Если подставить выражение (4.70) и в соответствии с условием (4.69) разложить корень в формуле (4.72), сохраняя члены, про­ порциональные р, то будем иметь

ДДтп= 17,37 __________ 0to) Re р____________

(4.73)

(1 + М )2 е У (& *)2 — (Со " )2 (1 — М 2)

 

где Со" —корень функции J'm{С0).

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

/

1+ |М|, \ 2

(4.74)

I

1-|М | ) '

 

Если частота близка к критической, т. е.

 

 

Ы

 

1 (Со Ф 0).

(4.75)

 

 

Со"“ V 1—М2

 

 

то формула

(4.54) для корней £ примет вид

 

 

; 2 = (

сГ )2- 2£’?

ка

' 0 ф руШ ),

(4.76)

 

(1 _ М 2 ) 2 £

 

а условие

|pG |< cl

в данном

случае может быть

переписано в

виде

 

 

 

 

 

 

 

1

!

(4.77)

 

 

(1 — М2)2

^

 

Затухание вблизи критической частоты может быть записано сле­ дующим образом:

ДД= 17,371/

---—----Im{e- ^ pvaj.

(4.78)

у

Е { 1 — М2)3

r

v

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

iPQ|»i*

Если попытаться в этом случае применить процедуру, подоб­ ную формулам (4.51) (4.54), то в результате получим алгебраи­ ческое уравнение относительно Z? более высокой ступени, которое не решается в радикалах, и поэтому поступим иначе.

Если импеданс равен нулю, т. е. р->оо, то корни уравнения

(4.39) совпадают с корнями функции Бесселя У(^ш)= 0. Когда параметр pG велик, естественно искать корни в виде

рСг

где

—корень функции J m{$n).

Подставляя (4.79) в уравнение (4.39) и разлагая функция Бес­ селя в ряд, при сравнении одинаковых степеней р<3 получим

6о“ , следовательно,

? = ( С ? - 2 ( С ? ~ = ( С у

2 Z W 1)*

(4.80)

О(С2) ’

 

рО

 

 

где 1 = р-1

безразмерный удельный

 

акустический

импеданс.

С уравнением

(4.80) можно проделать

итерационную

процедуру.

Внулевом положении С2 =(Ео,я)2-

Впервом приближении

c2= a D 2

2Z ( ^ у

(4.81)

а а о)

 

 

В следующем приближении

 

 

 

2Z(5?",)2

с2 = (С п)2—

 

(4.82)

 

2 Z ( f f ”)2

 

^(5Г)2-

(C )

 

0

При условии (4.69), т. е. для частот, значительно превышающих критическую для данной моды, из (4.81) и (4.45) для корней £2 получим формулу

С2=(5отя) - 2 / - ^ (С")2(1+ М)2,

(4.83)

k a

 

которая справедлива при условии

 

В _ (1+ М )2« 1 .

(4.84)'

k a

 

Затухание отдельной моды, отвечающее формуле (4.83), полу­ чим из уравнения (4.72)

Д/,= 17,37 ReZ (SH2(< + М)2

(4.85)

(1 — M 2 )(fta )2

 

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

[

1 +

|М |

\ 2

I

1 —

|М |

)

Вблизи критической частоты, т. е. при условии (4.75), из фор­ мул (4.81) и (4.45) следует

/ е т п )2

С учетом только старшего члена разложения для затухания от­ дельной моды из формулы (4.86) получаем

Д 1= 17,37 Im J/^

(4.87)

Вблизи критической частоты, как и ранее, затухание не зависит от направления распространения звуковых волн и воздушного потока.

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

4.3.РАСПРОСТРАНЕНИЕ ЗВУКА В КАНАЛЕ

СПОТОКОМ ПРИ НАЛИЧИИ ПОГРАНИЧНОГО СЛОЯ НА СТЕНКЕ

При движении воздуха в канале на его стенках образуется по­ граничный слой. Наличие пограничного слоя оказывает влияние на распространение и затухание звука, поскольку фронт звуковой вол­ ны может искажаться в пограничном слое. Для простоты рассмот­ рим направленный вдоль оси z плоский канал высотой а, у которо­ го одна стенка жесткая, а другая имеет импеданс Z и на ней су­ ществует пограничный слой с постоянной толщиной б (см. рис. 4.2). Поток является однородным по сечению канала и спадает по ли­ нейному закону до нуля на стенке при у = О, т. е.

I

иоУ

0

< г / < 8;

u { y ) = \

5

 

(4.8 8 )

I

и 0

8

< г/ < а .

Уравнение, описывающее распространение звука в плоском ка­

нале с осевым

потоком,

зависящим

только от у, имеет вид

 

1 д

, „ д

.2

р

^

4 “

-2г.

(4.89)

~ ~ + М 1 7 j р

ж '

ф

д у

д г

 

Решение уравнения для р ищем в виде

 

 

 

 

p= F {y)exp [i{kz z —со/)].

 

(4.90)

Предполагая, что v имеет такую же

зависимость,

воспользуемся

линеаризованным уравнением движения в форме

 

 

 

d v

, ,г

d v __

1

д р

 

 

 

d t

0

d z

е0

д у

 

 

В результате получим

F'y exp [t (kz z

taQ]

 

 

 

 

 

 

 

 

iQQCk

M 0-

k

 

 

 

где М= и/с; M0= U 0/c\

v

 

 

 

M= M0 y/8.

 

 

 

 

Подставив (4.90) и (4.91) в уравнение (4.89), получим

 

2 kzF' d

М

 

 

 

 

 

 

F" +

 

dy

 

- F [ { \ - t A k zf - k z] = Q.

 

k — M к г

 

 

 

 

 

 

 

 

 

Сделав подстановку

 

 

 

 

 

 

 

 

 

| =

1 —М —

= 1 —М0-^ -

 

(4.93)

 

 

 

 

k

 

 

k b

 

 

уравнение (4.92) можно представить в виде

 

 

dp £

dk _|_02 (\?2_

^&\ /Iг = 0 ,

(4.94)

где a = k4/MQkz.

 

 

 

 

 

 

 

 

 

Введя подстановку

 

 

 

 

 

 

 

 

/*■(£)= ехр

 

i ^ U

(ri),

r i = — ia?,

(4.95)

получим для U(t])

вырожденное

гипергеометрическое уравнение

dXJ ,

(--М

 

dU

 

~ ) U = Q.

(4.96)

т]— - 4 -

 

 

 

 

d-ц2

V 2

 

7

d-ц

\

4

 

 

Частное решение этого уравнения будет

 

 

 

и = Ф \

------ ;

 

1

 

 

 

-------;

aip

 

 

 

 

 

 

4

 

2

 

 

Удобно ввести функцию р (у) — переменный адмитанс, опреде­

ляемую следующим образом:

 

 

 

 

 

т

__

QCV

 

 

 

 

 

(4.97)

=

-

 

 

 

 

 

 

 

 

 

 

 

 

“ ' ( '

- • •

т т )

 

Если выразим

р через

\ из соотношения

(4.93), то выражение

(4.97) можно переписать

 

 

 

 

 

 

 

 

 

 

 

1

Р’ (6)

 

 

(4.98)

 

 

 

Р0 0

= <°£ F (£) *

 

 

Выражение (4.98) определяет нормальный адмитанс в сдвиговом потоке через давление и его производную.

Из уравнений (4.94) и (4.98) следует, что функция р(£) удов­ летворяет следующему уравнению:

 

P - —Z-

& =

Ь&р2 + i a -------—

£

£

которое эквивалентно уравнению

(4.94).

Поскольку

уравнение

(4.94) имеет частное решение вида

 

 

 

 

*2

 

 

 

 

ia ■ft2

1

1

2

(4.100)

/?(?) = е ''2 е,ф \ -

4

4

Y ’

~ 1а?

 

то уравнение (4.99) имеет частное решение следующего вида:

 

dR

 

Pi(5)=

dk

(4. 101)

Riai

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

ш = -

^ <е>

■ ш

(4.102)

 

;2

dS

 

 

с, + у

 

Л2(С)

Константу Ci определяем из граничных условий. При у = 0, т. е. при | = 1, имеем р = р0>где р0 — адмитанс стенки. Определив кон­ станту Си получаем

 

Р(6)= -

_________S//?2(£)____________

*j(6)

(4.103)

 

 

 

io&

 

 

 

 

 

 

 

 

 

Л2 ( 0

Ро i<sR(\)

Л2 dС

 

 

 

 

I

 

 

Полагая

у = 8, т. е.

£ = 1 — М0

получим

выражение для

адмитанса пограничного слоя

 

 

 

Р

 

 

6i/#2 (Si)

/CT£I/?(5х)

(4.104)

 

 

 

 

 

 

 

 

 

 

 

 

JoC2

 

 

У~5

Л2 (1)[р0 - Л ; ( 1)//оЛ(1)

r J Л»(O'

 

 

 

 

5

1

 

 

Поскольку в основной части канала dN[/dy=0y то решение урав­ нения (4.90) можно записать в виде

F (у) — A cos (а 0 + 0).

Подставляя это выражение в формулу (4.97), получим выражение для адмитанса при у=8

F'

ci sin а (5 — а)

ik{ \ — Щ ^ -^F

ik ^1 — Mo- — j cos а (В — а)

__

c c t g ( z ( a — 5)

ik

Приравнивая (4.104) и (4.105), получим

a a' tg a a '= — ika' р

где а' = а — Ъ.

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

Л?2 (5)= 1 + a 2S2^ 1 . Подставляя эти выражения в (4.104), будем иметь

( ' - М о ^ г Ь - щ р - м ^ - А )

(

кг

(4. 107)

2 *! \

1 -? о Ш ^ 1 -М о -^ + М0 — )

Как и следовало ожидать^ адмитанс зависит от угла падения зву­ ковой волны. Действительно, если ф — угол падения звуковой вол­ ны между направлением распространения и нормалью к поверх­ ности, то выражение (4.107) примет вид

 

sin ср

 

 

S i n

ср

 

 

sin2 ср

М0

) &>-»«( 1 -Мп

Мр sin

ср

(1 +

М0 sin ср)2 J

1 +

Mo sin ср

 

1 +

1 -

ЭоШ 1 -

Мо

S i n ср

М

 

 

sin 2 ср

 

Mo sin ср

з

(1

+

Mo sin

ср)2

 

 

1 - f

 

 

 

 

 

 

 

 

(4.108)

При нормальном падении

(ср = 0) будем иметь

 

 

 

 

 

п

__

Эо — №Ъ

 

 

 

 

(4. 109)

 

 

у ь

1— шро

 

 

 

 

 

что соответствует формуле для проводимости слоя малой толщины с нагрузкой р0.

Рассмотрим частный случай — распространение нулевой моды низкой частоты, для которой можно получить аналитическое реше­ ние. Если |Aa|3|<Cl, то* в знаменателе в правой части формулы (4.107) вторым членом, малым по сравнению с единицей, можно пренебречь. В результате уравнение (4.106) примет вид.

a0a tgала = -ika% 0{ \ — М0- ~ ')

2

— k2ab

z 0

1 — М0

 

где а! жа.

Решая приближенно это уравнение, получим

I + J t )

2а

(4.111)

Выделяя мнимую часть уравнения (4.111), получим следующее выражение для определения затухания звука нулевой моды:

AL = 8,68

Re $ !У Ъ

+

 

(I ± М0 / б )2 +

M§rf2

 

462

 

1 + Ь

±

зм0(1 ± м0уъ —6)

 

 

 

V I

 

 

1 ± Мр Уь_________ I

(4.

112)

 

(I ± м0/б)3 +

(1 ±

4М0 У ь ) 2 j

 

 

 

 

46

(1 ±

Мо УЪ)2 .1

 

 

 

Rep

 

Ь +

 

 

где

Ь =

 

Imp

 

 

 

 

 

k a

 

 

ka

 

 

 

 

Верхний злак относится к распространению по потоку, а ниж­

ний — против потока.

 

 

 

 

 

 

 

Анализ формулы (4.112)

показывает,

что при распространении

звуковой волны

по потоку

наличие

средней

скорости приводит к

уменьшению затухания, а наличие пограничного слоя — к увеличе­ нию затухания. При распространении звуковой волны против по­ тока наблюдается обратное явление: средний поток увеличивает затухание, а пограничный слой уменьшает, причем влияние погра­ ничного слоя в этом случае более существенно. Таким образом, на­ личие пограничного слоя на стенке канала ослабляет влияние по­ тока на затухание звука.

Формула (4.107) для импеданса в случае плоского канала при достаточно малой толщине пограничного слоя по сравнению с ра­ диусом может применяться и для цилиндрического канала.

Действительно, в случае цилиндрического канала радиусом а, уравнение, аналогичное уравнению (4.92), имеет вид

 

d

M

d t F .

2 k z

d r

d F

dr% d r

( k — M k z )

Так как r = a —yb, у — (а — г)1Ъ, то

d F

___ 1_

d r

5

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

dF_

d W

1

d*F

d y

drl

52

rfi/2

 

 

 

ud ivм1 \

 

 

 

 

 

 

25 ----k z

 

 

 

d^F

d F

5

_ -^-— \^-Fb2\{k-bAkzf - k 2z------- ^ — 1 = 0;

d y t

d y

\ a yb ^

k

— М Л*/

L

(a — 1/5)2 J

 

Пренебрегая b2/a2, получим

 

 

 

 

 

 

d M

 

 

 

d^F

 

d F (

25—— kz

 

 

 

 

 

d r

[(ft_ M * J 2 - £

l - - ^ l = 0 .

 

d y 2

 

----1 ----- г ------------- +F& 2

 

 

d y \ a T Л - М Л , J

V

°-2 J

 

 

 

 

 

 

 

(4.

114)

Решение этого уравнения сводится к аналогичному решению для плоского канала при условии

 

 

 

d М

 

 

2 5 -------k z

 

 

sd

 

 

а

k — М kz

или, если

d М

М0

 

d r

, ТО

 

 

 

 

 

 

в

2Mpkz

a^ k — M k z

4.4.ОПТИМИЗАЦИЯ ЗАТУХАНИЯ ЗВУКА В КАНАЛЕ

Сточки зрения снижения шума исследование распространения звука в каналах .направлено прежде всего на выбор звукопогло­

щающей конструкции (ЗГ1К), которая обеспечивала бы максималь­ ное затухание. Возникает вопрос, какими критериями руководст­ воваться при выборе импеданса стенок канала, который должен быть реализован с помощью конкретной ЗПК. Одним из таких мо­ жет быть критерий, основанный на оптимизации затухания звуко­ вых волн в канале.

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

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

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

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

В цилиндрическом и кольцевом каналах моды различаются ази­ мутальными и радиальными числами, причем при фиксированном азимутальном числе т могут быть слиты любые две соседние ра­ диальные моды (п, п+ 1). Это могут быть как наименее ослабляе­ мая, так и наиболее возбуждаемая моды. Таким образом, под оп­ тимизацией понимаем слияние мод, т. е. образование двойных или даже тройных мод, которое происходит при определенных дискрет­

ных значениях импеданса Z„pf, называемых оптимальными.

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

импеданса

из множества оптимальных (а следовательно, и вы­

бор ЗПК)

в случае каналов авиационных двигателей непосредст­

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

4.4.1. Двойные моды

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

/ М ) = 0 ,

(4.115)