Аэродинамические источники шума
..pdfВ соответствии с преобразованием координат имеем
|
|
|
|
|
. |
, |
д |
. д |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
Sin ф------- cos Ф-—= — |
<?Ф’ |
|
|
|
|
|||||||
|
|
|
|
|
|
|
дх |
ду |
г |
|
|
|
|
||||
поэтому |
|
|
|
|
|
|
|
|
|
|
mnQ |
|
|
|
|
|
|
|
|
|
|
го 2л |
|
|
|
|
|
. |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
С о |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
rdrdty- |
(3.69) |
|||
|
PQ= - ~ |
|
\ \ |
в те тп^ ~ ^ — — I |
D |
|
|
||||||||||
|
|
|
|
ОО |
|
|
|
дф |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Вычислим производные под знаками интегралов |
в |
выражениях |
|||||||||||||||
<3.68) и (3.69) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
I |
_/»2* , |
\ |
|
|
_ /2 2 i. |
|
|
|
|
|
|
|
|
|
||
д |
I е |
С° |
|
I______ е |
с° |
/ imnQM0 |
, |
imtiQz |
. |
г |
\ |
||||||
дг |
' |
D |
|
. mnQ |
v |
|
в |
V |
с0Р2 |
|
|
№Ос0 |
' № |
}’ |
|||
|
|
( |
|
|
|
|
|
|
|
|
,тпЯ |
|
|
||||
|
|
|
'~1~7Гв\ |
|
х sin 4/ |
/ imnQ |
, |
|
р2 \ |
° |
|
||||||
|
|
|
е_______ 1 |
|
|
1 |
Св |
|
|||||||||
|
|
|
|
D |
|
) |
|
D2 |
{ |
с0 |
**” |
D ) |
|
|
|
|
|
Подставляя |
|
значения |
для производных |
в |
|
(3.68) |
и |
|
(3.69) и |
||||||||
складывая |
эти |
|
выражения, |
получим |
выражение |
для |
звукового |
||||||||||
давления |
|
|
|
|
|
~ Шп (+ + ,7 |
°)(Г Л „ (imnQMo |
imnQ_Z , £\"| |
|||||||||
Ртп — \Pi-\~PQ — |
|
|
|
|
|||||||||||||
|
4я |
) |
) |
|
|
|
Цд { |
СоР2 |
|
C0 p 2 D |
' D 2 jj + |
||||||
|
|
|
|
|
о |
о |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
rdrdb |
|
|
|
(3.70 |
Полученное уравнение позволяет определить звуковое давление в любой точке пространства при заданном распределении аэродина мической нагрузки, т. е. при известных Ат и Вт.
Для определения звукового давления в дальнем поле выраже ние (3.70) может быть упрощено. Для этого представим выражение для величины D в виде
D-- |
л 202*r cos Ф |
I |
P2r2 |
|
D\ |
1 |
Dl |
||
|
где Д)= Y z 2+ Р2-*2-
иметь
D. |
1 / |
2 ^ 2 Х Г COS Ф |
1 |
Di
Разлагая -подкоренное выражение в ряд по малому параметру (второй член под корнем) и пренебрегая величинами второго по рядка малости и выше, получим
p2jcr COS Ф
£ > « А . (3.71)
Da
Следовательно,
? |
M.QZ + Do |
xr cos ф |
(3 72) |
p2 Do
Подставим (3.71) и (3.72) в выражение (3.70). При этом будем считать, что D M Do, когда величина D определяет расстояние меж ду источником и точкой наблюдения. Однако в экспоненте, которая учитывает сдвиг фаз возмущений, приходящих в точку наблюдения от отдельных точек винта, будем пользоваться для D более точным выражением (3.71). В результате получим
|
г0 2л |
|
—I m nQ |
( M 0z + D |
x r cos ф\ |
|
|
£i(mnQt—mnty) g |
c0 |
V Pa |
D 0 ) |
|
|
Рт- ■ т Л \ |
|
|
|
|
|
|
|
о о |
Z \ 1 в тх sin |
|
|
|
|
/ imnQMg |
imnQz |
(imnQ | ft2 |
rdrdty. |
|||
xlt( c0p2 |
CQ$2JDQ |
1 |
л 2 |
|
|
|
|
|
|
Пренебрегая в амплитуде величинами второго порядка малости, будем иметь
|
mnQ |
im n Q |
|
---- {MQZ-\-DQ) го2к |
р(л*+£)+ |
|
P m n = i |
е |
|
СоРа |
|
||
|
|
|
||||
|
4яЛ0с0 |
|
|
0Ш0 |
|
|
|
|
|
|
|
|
|
|
|
£8шф] |
‘ ЩГ, Г-ГС08ф- 'тПф rdrdty. |
(3. 73) |
||
|
|
Do |
J |
|
|
|
Для вычисления аэродинамической нагрузки на лопасти рото ра, т. е. Ат и Вт, можно воспользоваться общими выражениями,, приведенными в подразд. 3.4.1. Однако здесь рассмотрим более простой случай, когда аэродинамическая нагрузка не зависит от угла ф. Используя свойства бесселевых функций, запишем
|
2— |
m nQ |
rx |
, . |
/ rn n O |
|
\ |
|
|
• |
l -=z----- |
COS Ф -./Ш П Ф |
|
||||
2* |
0 |
|
|
|
|
|
|
|
* |
, m nQ |
|
. |
2« r v |
|
I”"* |
||
^ sin |
^ |
‘ |
_ a a |
|
||||
|
|
|
Qrx |
|
mn I |
CQD Q |
Таким образом, выражение для звукового давления т-й гармоник» будет
__ ; т п + 1 |
mnQ е |
imnQ ( i |
AIQZATDQN |
' |
с0Ра ) X |
||
Р т п~ |
2DoCo |
|
|
х S [ * • Ы м°+ £ ) ■ - ® - 5 - 1 ( т £ г xr) rdr■ (3- 74)
Примем аэродинамическую нагрузку в виде сосредоточенной силы, приложенной, как это сделал Гутин, на радиусе гг«0,8/о и учтем выражения для тяги из (3.67а)
силы сопротивления вращению из (3.676) и (3.67в)
В результате выражение для определения максимального звуково го давления в дальнем поле примет вид
< 3 ' 7 5 )
В отсутствии скорости поступательного движения (Л4о=0) бу дем иметь
Р = 1, Д з=/?, Z / D Q= sin о, xfD 0= cose. Следовательно, выражение (3.75) примет вид
[ г й " “ - « 5 ; ] ( * ? * |
■> |
<3 - 76> |
Учитывая, что Q=MKV/rif видим полное совпадение выражения (3.76) с (3.66), полученного в подразд 3.4.1 для свободного рото ра без поступательной скорости.
3.5. ШУМ РОТОРА В КАНАЛЕ
Решение задачи по излучению шума ротора в канале вызвано необходимостью определения шума компрессора и вентилятора {57, 84, 101]. В общем виде компрессор можно представить состоя щим из ротора, расположенного в канале, перед ротором и за ним установлены решетки лопаток направляющего и спрямляющего ап парата — статора.
При решении задачи об излучении звука такой системой исполь зуем примененный в предыдущих разделах метод. В этом случае решение опять будет полностью определяться возможностью опре деления сил, с которыми лопатки ротора и статора воздействуют на воздушный поток в канале. Как и ранее, используем здесь прин цип суперпозиции, когда при наложении нескольких сил F{ звуко вое давление также суммируется. Последовательность импульсов сил, действующих в отрезок времени At вдоль траектории движе ния лопаток,, может быть различной.
Если в канале установлен ротор без статора, то здесь может »быть применено решение, приведенное в разд. 3.4. Если ротор с числом лопаток щ имеет впереди статор с числом лопаток я0, то в результате взаимодействия следов за лопатками статора с лопат ками ротора возбуждаются моды колебаний. Такие колебания соз даются круговым рядом источников, имеющих одинаковые частоты
и амплитуды, но разные фазы. Про изведя интегрирования звуковых давлений от всех элементарных ис точников в сечении потока, можно' получить общее решение для звуко вого поля в спектральной форме, используя разложение на гармони ческие составляющие.
Рассмотрим установленный в ци линдрическом канале осевой ком прессор (рис. 3.8), который сообща ет потоку поступательную скорость U0 и закручивает его с угловой ско
Рис. 3.8. Система координат для ростью Q [8 ]. Воспользуемся урав компрессора
нением Блохинцева (1.25), которое
вцилиндрических координатах г, ф,
хдля звукового давления р, распространяющегося в канале с од нородным и закрученным потоком, примет вид
д2р |
± |
( |
r |
Jd L \ |
+ - ! - ? E = |
( ± |
± + к |
± + |
± |
± ) 2 р . (3.77) |
± |
||||||||||
дх* |
г |
дг |
\ |
д г ) |
г2 дф2 |
\ с0 |
dt |
дх |
с0 |
дф / |
Решение волнового уравнения запишем в виде |
|
|
||||||||
|
|
|
|
|
р = р ( г )е цк*х+^ - ш‘), |
|
|
(3. 78) |
где kx — волновое число в направлении оси х. Подставив (3.78) в (3.77), получим
( 3 ' 7 9 )
где k = iо/с; M.u= QD/2c0 — окружное концевое число М рабочего колеса; D=R0—г<>; Ro — внешний радиус лопатки рабочего колеса; г0 — .радиус втулки. Обозначим
* , [ ( , _ M |
i |
— |
|
(3. 80> |
|
тогда уравнение (3.79) примет вид |
|
|
|||
d2£ + _L d£_ + A2o_ ^ |
о. |
||||
дг* |
' |
г |
дг Г\, |
г2 ) и |
|
Это уравнение Бесселя. Таким образом, общее решение урав |
|||||
нения (3.79) будет |
|
|
|
|
|
/ > = 2 |
2 |
v |
e' <v+s*“' ' ’ |
(3.81) |
|
|
q=>— <x а = 0
где Rg{Kqar) — линейная комбинация функции Бесселя Jg(x) и Неймана Nq(x) порядка q, удовлетворяющая граничным условиям.
в радиальном направлении; Аяп — постоянные коэффициенты, ко торые следует определить из граничных условий.
Для определения Aqa рассмотрим упрощенную модель компрес сора (сМ. рис. 3,8). Пусть в бесконечной цилиндрической кольце вой трубе вращается с угловой скоростью Q система вихревых ли ний, причем их число щ равно числу лопаток ротора и расположе ны они в плоскости х=0. Примем угловое положение ф=0 за начальное. Будем предполагать, что здесь находится какой-либо, источник возмущения, например лопатка неподвижного направляю щего аппарата. При прохождении каждого вихря через это поло жение происходит импульсное изменение окружной составляющей
.скорости, которое можно записать в виде
£/ф= / ( г ) 8 ( ф - ^ - ) , |
(3.82) |
где f(r) — радиальное распределение окружной скорости. Дельта функция принимает следующие значения:
|
( |
. 2лщ |
|
оо |
при ф = ----- - ; |
8 (♦ |
2nni |
по |
п о , О при ф ф — п— |
По
Используя уравнение (3.82), определим циркуляцию по радиусу лопатки
2тс |
|
■Г(г) = | U^rclty= f {г) г. |
(3.83) |
Рассмотрим взаимодействие п\ лопаток рабочего колеса с п0ло патками статора. Если лопатка статора с индексом «п» располо
жена под углом Фп ~ ~ — » гДе л = 0 , 1, 2 , (п0—1 ), то изменение
По
-окружной скорости можно записать в виде
(3.84)
где т — целое число.
Просуммируем это выражение с учетом всего числа лопаток гц статора. Используя свойство дельта-функции
2л г л = - ! - Y e " |
|
По J |
2л Q=>— оо |
запишем выражение (3.84) для п0 лопаток в виде
Я р— 1 |
со |
|
|
/ (г) |
exp ^iq ^ |
— ' ^ - j Jri2nnm — — i/ra/ZjQfj. |
(3. 85) |
п=0 |
|
|
|
Это выражение можно упростить, так как |
|
||
Т1о— 1 |
/г0, если /П/21— g целое число; |
|
|
S 6XP[f?f'(/ra/il_^]==~ |
^0 |
|
|
/яяг—а |
|
||
я=»О |
О, если — -—- ■ не равно целому числу. |
||
|
|
щ |
|
Обозначим ■— 1 ~~ q zPi, откуда |
|
||
m |
q = m n l —pinQ, |
(3, 8 6 ) |
|
|
|||
где р\ = 1, 2 , 3, |
— целые числа. |
|
Соотношение (3.86) между числами лопаток ротора и статора является весьма важным, так как оно определяет порядок мод ко лебаний. Впервые оно было приведено в работе [101]. После пре
образования уравнение (3.85) примет вид |
|
|||
U ^= |
по |
ехР |
~ imn^Qt) = -Г |
exp (iqty —imnfit). |
|
|
— оо |
|
тп=>— оо |
|
|
|
|
(3; 87) |
С другой стороны, окружную скорость можно определить из уравнения сохранения количества движения для закрученного по тока, которое в цилиндрической системе координат в направлении ф имеет вид
—Л.Ц ?U_JLQ;dlL —__ L AL . |
(3* 8 8 ) |
||
d t |
0 д х |
Qo |
|
Решение этого уравнения ищем в. форме |
|
|
|
|
U — U (г) е/(Л^ +<7ф_а)0. |
|
|
Подставив это выражение в (3 .8 8 ), |
получим |
после преобразо |
|
вания |
|
|
|
д_Р_:
и = |
—-------- ^ ------ . |
|
е0г <a.—LT0kx — qQ |
Производную др/д-ф определим из (3.81) |
|
% = i V ! 'S ) |
{\аГ) exp i { k xX + q*t — (tf).. |
Подставив это выражение в (3.89), получим
и 2 2 ^ ^ ° (Kqar) exp i (k x x + qA/— v>t)
____
бол ((о— U~0kx — q Q )
(3i 89)
(3t 90)
Приравняв (3.87) и (3.90) при х=0, будем иметь
а^О
Умножая это выражение на rRa(Xqar) и проинтегрируя от Го до Ro, получим для каждой qa моды
Ro |
|
Ro |
|
А<а \ R l ^ r d r ^ - U |
o k t - q Q ) ^1 |
^ Г {r)Ra{ \ ar)rdr. |
|
То |
|
Г0 |
|
Откуда |
|
|
|
= ^ |
% U\ t r ~ |
\ T ^ R^ a r ) r d r , |
(3. 91) |
|
|
То |
|
Ro
где lRa\ = § Rl(Kqar)rdr. Го
Подставив (3.91) в (3.81), определим звуковое поле компрессо ра в канале.
Г л а в а 4
РАСПРОСТРАНЕНИЕ ЗВУКА В КАНАЛЕ С ПОТОКОМ
Основное внимание в предыдущих разделах было уделено воп росам образования звука турбулентным потоком, считая эти зада чи наиболее актуальными в аэроакустике. Однако в случае, когда источник звука уже известен, возникает другая проблема, обуслов ленная необходимостью решения задач по распространению звука от источника до наблюдателя. Важность этой .проблемы для аэро акустики подчеркивается тем фактом, что в авиации постоянно имеем дело с движущимися объектами. Кроме того, изучение рас пространения, так же, как и образование звука, вызвано необхо димостью поиска методов снижения шума. Эта необходимость еще больше связывает два направления.
Воздействие шума может быть уменьшено как благодаря сни жению шума в самом источнике, так и по пути его распростране ния. Первый путь — наиболее действенный метод уменьшения шума — для источников шума турбулентного потока непосредст венно следует из рассмотрения физических основ образования зву ка, изложенных ранее в гл. 2 и 3. Несмотря на значительные ус пехи по уменьшению шума в источнике, в авиации пока не удалось создать достаточно малошумных силовых установок, обеспечиваю щих требуемый уровень шума на местности и в кабине самолетов. Поэтому для снижения шума в каналах силовой установки самоле та применяют звукопоглощающие конструкции. Следует отметить, что глушители шума устанавливаются также и в трубопроводах системы кондиционирования летательных аппаратов. Вследствие этого при рассмотрении распространения звука в движущейся сре де большое внимание уделяют задачам распространения звука в кан-алах.
4.1.УРАВНЕНИЕ РАСПРОСТРАНЕНИЯ ЗВУКА
ВКАНАЛЕ ПРИ НАЛИЧИИ СДВИГОВОГО ПОТОКА
Вгл. 1 были получены уравнения, описывающие распростране ние звука в потенциальном и изэнтропийном потоках. Однако в полной мере они не могут быть использованы при изучении распрос транения звука в каналах, поскольку воздушный поток, сущест вующий в каналах, нельзя считать потенциальным из-за наличия пограничного слоя. Завихренность, существующая в пограничном слое, в некоторых случаях оказывает существенное влияние на за-
туханйе звука. Поэтому представляется целесообразным получить уравнение такого же типа, как в разд. 1.2 , описывающее распрост ранений звука в завихренном потоке. Рассмотрим не произвольный завихрецный поток, а его частный и в то же время очень важный для акустики каналов случай — так называемый сдвиговый по ток, в котором скорость зависит только от поперечных координат. Это уравнение легко можно получить,из уравнений, описывающих распространение звука в неоднородной движущейся среде, рассмот ренных Д. И. Блохинцевым {5]. Эти уравнения, хотя и не учитыва ют влияние вязкости и теплопроводности, являются весьма общими за исключением двух приближений: они линейны иописывают толь ко адиабатические процессы. Тем не менее решения, полученные на основе этих уравнений, достаточно хорошо описывают реальную действительность, и поэтому они находят широкое применение при решении ряда практических задач [82]. В виду важности этих урав нений приведем кратко вывод их исходя из общих уравнений гид родинамики:
уравнение движения сжимаемой жидкости
уравнение неразрывности |
|
|
|
|
|
|
d t |
-|- v |
* |
+ Q ^ |
= 0 ; |
(4. 2) |
|
; |
д х j |
|
OX j |
|
|
|
уравнение для энтропии |
|
|
|
|
|
|
|
d s |
, |
dS |
=0 |
, |
(4.3) |
|
------г |
д х |
||||
|
d t |
1 |
|
|
|
где Vj — скорость; Q — плотность; р — давление; S — энтропия. Уравнение состояния удобно представить в следующем виде:
е = б (А S). |
(4.4) |
Дифференцируя его, получим
(4" 5)
Будем считать, что процесс адиабатический, следовательно,
- £ 1 = 0 . dS Iр
Обозначим
( * ) - А
где с — скорость звука.
Тогда из уравнения (4.5) будем иметь
-^ - = -L |
(4 . 6 ) |
d t |
& d t |
1 |
Представим скорость Vi в виде суммы двух величин: средней Uoi и дополнительной составляющей щ, (причем U i^U ^. Такие же предположения сделаем для давления, плотности и энтропии:
^/ = ^ 0/ + ^ л |
<<^£/0/; |
Q=eo+e'. |
б' <^ео; |
Р=Ро+Р'> |
Р’^Ро^ |
S = S 0+ S ', |
S ' ^ S 0. |
Подставим |
эти выражения в исходную систему уравнений |
|
(4.1) |
(4.3) |
и проведем линеаризацию этих уравнений. Для пер |
вого приближения получим систему линейных уравнений, которая описывает распространение звуковых возмущений в движущейся среде:
(4.7)
Уравнение состояния (4.4) представим в виде |
|
^ = C2 Q' + A5', |
(4.8) |
где
Система уравнений (4.7) и (4.8) содержит пять неизвестных, и она весьма сложна для решений. В соответствии со сказанным рассмотрим распространение звука в однородной, движущейся вод ном направлении среде, скорость которой зависит только от попе речных координат. В этом случае из системы уравнений (4.7) мож но достаточно просто получить волновое уравнение относительно давления. Итак, рассмотрим распространение звука в среде, кото рая движется в направлении оси Х \ = г со скоростью U0(x2= xt х3= =у). Относительно скорости UQпредполагаем, что она сравнима со скоростью звука, но меньше ее. Скорость потока в уравнениях (4.7) можно представить в следующем виде:
U0i= bu U0{x2, х 3), |
(4.9) |
где
Примем также
р0= |
const; |
|
Q0 = |
const; |
(4. 10) |
S0= const.