Аэродинамические источники шума
..pdfцей, что в качестве |
вместо |
нулей |
производной функции |
Бессе |
||||||
ля следует понимать корни уравнения |
|
|
|
|
||||||
|
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 |
в данном |
случае может быть |
переписано в |
|
виде |
|
|
|
|
|
|
|
IР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) |