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

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

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

В соответствии с преобразованием координат имеем

 

 

 

 

 

.

,

д

. д

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

 

 

 

 

 

)

)

 

 

 

Цд {

СоР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 , (п01 ), то изменение

По

-окружной скорости можно записать в виде

(3.84)

где т — целое число.

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

2л г л = - ! - Y e "

По J

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.