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

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

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

Таким образом, пульсации подъемной силы на данной частоте со определяются возмущениями скорости с волновым параметром kx= ы/и0. Взаимная спектральная (плотность пульсаций подъем­ ной силы может быть определена из выражения

о

/

 

% 1 ,

AY (jfi, ш) Д Y* (Х2, w)

,

(3.33)

^Ayar(-yi,-y2,a ))= lim n ^ v

------

^

 

л

 

 

 

 

 

 

 

где АК (лгх, со)= —?—Г

ЬУ (x,t)Q-Mdt.

 

 

 

 

J

 

 

 

 

 

 

 

 

—т

 

 

 

 

 

 

Подставляя (3.32) в (3.33), получим

 

 

 

 

SbYAr (хи х 2, w)=(nQ0b)2 g (xv kx) g* (x2, kx) lim я

v

.

(3.34)

 

 

 

 

 

T—>-oo

1

 

 

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

(3.34)

представляет

собой

спектральную плотность энергии пульсаций скорости

 

 

 

ф™(*л-) = П т2л

V (kx) V*(kx)

 

(3. 35)

 

 

Ах

 

 

 

 

 

Ах-+оо

 

 

 

 

где &x=U0T и Т —юо.

 

 

 

 

 

 

Следовательно, уравнение (3.34) примет вид

 

 

 

5дГАг(АГ1,-«2,и))=

(лео^)2 -y-® ^(^)g(-«l.^)g*(-«2.^)-

(3.36)

Подставляя выражение

(3.36) в (3.26) и учитывая, что

 

 

 

.сог

.

 

пшг

d

 

 

 

 

 

sin2—

 

 

 

~ l— {Zi—Zi)

 

------— — = <Р

 

(3.37)

 

Ok

С о Г

dzl(jLz2=

 

 

 

 

 

 

/ сoz \2

 

 

 

 

 

 

 

\ сог )

 

 

 

в плоскости 2 = 0, получим

 

 

 

 

 

 

 

 

 

ф^ { к х)\0(К)\\

 

(3.38)

ь

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

где G(kx)= f g{x,kx)dx.

__b

2

Уравнение (3.38) можно представить в другом виде, характеризуя возмущения набегающего потока масштабом

 

оо

 

lz ^ ^ £> ,1 m

f Qw (k*'z ) dz>

(3- 39)

Rvv (kx* 0)

J

 

 

0

 

где Rvv(kx, z) — пространственно-волновая функция корреляции турбулентности в потоке; Rvv(kx, 0) — автокорреляционная функ­ ция.

Используя выражения для взаимной спектральной пульсаций скорости в виде

оо

^vv (^х) ^ Rvv (^л*» ^0

О

и спект|ральной плотности

(“>) = Rvv (kx, 0)/Uо,

получим

Qof/\и>ЧЫ

/(г , ш)

32лс§г2

К? ( * * ) 1 Ч М $ * ( “ )•

 

Проведем размерный анализ полученных выражений.

мощности

(3.40)

-— — bkx ^>— \ высокочастотном приближении функ­

ция G(kx) может быть представлена следующей зависимостью [55]:

 

 

Q{kx)-

(3.41)

Мощность

W,

излучаемая профилем в полосе частот Да, пропор­

циональна

/ ( а ) Да. Вводя характерные величины для скорости V,

размеров L и частоты a~V /L и учитывая (3.38) и

Ф(АЖ)— LV2, |Ох

х ( ^ ) |~ у = р

получим

 

W ~ QoVlL2 = Q0V3L2М2, С0

де М = U 0lcQ.

Б. При

( j

 

низк°частотном приближении функ­

ция G(kx)

может быть

представлена

следующей

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

[56]:

 

 

 

 

 

 

 

(т=*!*-) Г/0

\(1 - М2)

\ _ и

/-М »М

 

V|1 -

М2| / [

/

\ (1 -

М2)

 

 

 

 

 

 

(3.42)

где S ( —kxb) — функция Сирса [98]; / 0, / 1 — функции Бесселя, пер­ вого рода нулевого и первого порядка соответственно. Поскольку здесь G(kx) безразмерна, то

W

= Q0V3L2№

С0

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

На рис.

3.4

приведен *

спектр

* <!-

о*8 «1^,

I о-Эксперимент'^

шума профиля в 1/3-октавных поло­

 

 

1 ®

Расчет

 

сах частот

(200... 10 000 Гц)

в точ­

 

 

 

°1го. 1

 

 

ке,

расположенной

прямо

над

про­

 

 

 

 

 

филем на расстоянии

у = 2 м. Про­

 

 

 

!

°! о в .

 

 

 

 

!

о

• •

филь представлял собой часть ло­

 

 

 

 

 

 

j

1

° Ч

*

пасти модели несущего винта вер­

 

 

 

толета с хордой

й =

180 мм и отно­

 

 

 

!

i

 

 

сительной толщиной

12%.

Профиль

 

 

 

 

 

устанавливался

в области

потенци­

100

200

•500

1000

2000

5000 f j и.

ального ядра струи, истекающей из

 

 

 

 

 

 

 

сопла диаметром 400 мм, на рассто­

Рис. 3.4. Спектр шума

обтека­

янии 200 мм от среза сопла. Здесь

 

 

ния

лопасти

 

 

же

приведен спектр,

полученный

 

 

(3.40).

Ввиду

того,

расчетным

путем на основании

формулы

 

что

профиль

обтекался

неоднородным

 

потоком,

вклад в

суммарный уровень

шума

от

различных

сечений профиля

был

различным. Поэтому на основании гипотезы плоских сечений

про­

филь разбивался на участки с шагом

10 мм, для которых все па­

раметры набегающего потока U0, Svv(со), lz(со) считались постоян­ ными по размаху и определялись экспериментальным путем. За­ тем проводилось сложение интенсивности шума от каждого элемен­

та

профиля. В рассматриваемом диапазоне частот справедливо

-

}> — , поэтому при расчете использовалась функция реак-

ции профиля в форме (3.41). Сравнение результатов расчета с экс­ периментом показывает хорошее совпадение в диапазоне частот f —200 1000 Гц, а при более высоких частотах наблюдается рас­ хождение. Несовпадение на высоких частотах объясняется тем, что выражение (3.41) было получено для тонкой пластины, оно не учитывает толщину и кривизну профиля.

3.4. ШУМ СВОБОДНОГО РОГОРА

Основным элементом применяемых в авиации лопаточных ма­ шин таких, как самолетные винты и винтовентиляторы, несущие и рулевые винты вертолетов, вентиляторы, компрессоры и турбины реактивных двигателей, является ротор, т. е. набор профилирован­ ных лопастей или лопаток. Обтекание потоком вращающегося ро­ тора приводит к интенсивному излучению звука [89, 107]. Если ро­ тор помещен в канал, диаметр которого -существенно больше дли­ ны волны излучаемого звука, что имеет место для вентиляторов, компрессоров и турбин реактивных двигателей, то влиянием кана­ ла на процессы генерации звука можно пренебречь. Наличие кана­ ла, в котором вращается ротор, сказывается, в основном, на усло­ виях распространения и излучения звука из канала. Таким обра­ зом, процессы шумообразования в лопаточных машинах можно рассматривать с единых позиций. При. этом в качестве элементар-

* Измерения выполнены А. В. Топоровым и П. С. Антохиным.

необходимое для распространения звука ла расстояние D со скоро­ стью с0.

 

Поскольку

 

 

 

 

 

 

 

 

 

 

x — R cos a cos б; г/=/? coso sin б; z= ./?sino,

 

то

D — У х 2-\-у2

z2-\-г2 —2г(у sin ф-|-.х: cos ф)=

 

 

 

— У R2 + r2 — 2Rr cos а cos (ф — б).

(3.45)

 

Запишем выражение (3.44) в декартовых координатах

 

р {х, у,

=

eimnSl J

гд Р

 

p

 

__p

dD \ у

 

д х

 

 

 

 

 

 

 

 

т

д у

mz

d z )

 

 

х

д ( 1

 

 

)} ,5 .

 

(3.46)

 

 

 

dD \ D

 

 

 

 

 

 

 

Воспользовавшись выражением (3.45), получим

 

 

 

d D

 

R

Q

г

,

 

 

 

 

 

----=

— cos a cos 0 -------cos Ф;

 

 

 

 

 

д х

 

D

 

D

т

 

 

 

 

 

d D

 

R

л

г

• .

 

 

 

 

 

----=

— cos

а sm 0 -------sm Ф;

 

 

 

 

 

д у

 

D

 

D

 

 

 

 

 

 

d D

 

R .

 

 

 

 

 

 

 

 

----= — sm а .

 

 

 

 

 

 

 

d z

 

D

 

 

 

 

 

 

 

Следовательно, уравнение

(3.46) примет вид

 

 

 

/

1

im n 9 ( t ~ )

 

 

 

 

 

 

p{x-y-Z J ) = M

e

 

Я Ч 1cos о cos 0 --- D cos

 

~\~Fmy

cos о sin 6

sin ф^ + Т7,

^

sin о j

X

 

 

 

x

 

e-WdS.

 

 

(3.47)

Таким образом, задача сводится к определению сил Fmx, Fmv, Fmz. Рассмотрим случай равномерного распределения сил вдоль хор­ ды элемента лопасти. Временная зависимость импульса силы на элементе лопасти rdrdty приведена на рис. 3.6, а. Сдвинув начало координат в середину хорды (рис. 3.6, б), разложим выражение

для силы в ряд Фурье

 

■ Р (г,ф ,^ )= 2 [ап cos (тпШ')-\-Ьтsin (tnriQt')],

(3.48)

 

m-1

 

 

TZ/ПЯ

 

где

ат = ^ - \

F{r, ф, t') cos (mnQt')dt';

(3.49)

 

Ьл

J

 

 

—«/л2

 

 

Tz/nQ

 

 

t6m = - ^ - \

/г(г,ф, t') sin (mnQt')dt'.

(3.50)

 

Ьл J

 

 

 

—%/nQ

 

f(!i

f(t)

A +

 

JL 9

V b

. V

t

b

0

Я + Я r

Z°+ H

* * *

- 2 Я 7

19 r

 

a)

 

 

S)

Рис. 3.6. Схема импульсов, действующих на ротор

Учитывая, что

 

 

 

 

 

 

 

О

при

nQ.

^

^

----

 

 

F

°2r'пQ ' ’

F { r , ^ n =

А

при

__ 1 _ t ' t '

b

2rQ

 

 

 

2rQ

 

 

 

О

При

<

t' <

——

 

 

v

2rQ

^

^

nS.

заменяя в (3.49)

и

(3.50) пределы интегрирования, получим

 

 

а„

2А

( m nb

\

 

 

 

----sin »------ >*

 

 

 

т

wL

(—

)•

 

 

 

 

m nb

 

 

 

 

ьт= о.

 

 

 

Следовательно, .выражение (3.48) будет

 

 

 

 

оо

 

 

 

 

F(r, ф,/') =

—— sin [

\ 2r

cos (тпШ'),

(3. 51)

 

 

m=>1 m n b

J

 

где Л = /(г,ф).

 

 

 

 

 

 

Используя выражения (3.43) и (3.51), получим

 

2

, ,

. ч . I тпЬ\

(тл2/ - 9ф -А )

 

f ‘= W

/ ( r ’*) s m (—

) е

<3-52)

Запаздывание по времени (/ — появилось вследствие перено-

са начала отсчета в середину хорды Ь. Составляющие силы Fi будут

F y {r,bt)

2

/(r,+ )sin

cos фе

 

(3. 536)

nmb

 

 

 

 

 

 

 

 

 

 

 

_ . .

 

2

l4 .

( m n b \ « (mn2t-qC/— )

(3 .5 3 B)

F z ( r ^ J ) = ~ b f { r ^ ) s ™ \ ^ — \ f i \

2 r >

Таким образом, выражение

(3.47) будет

 

 

X f (/">40 [cos о sin (<|>— 0) -f- sin o] |-j—_|_г ^?Н| rdrdty. (3.54)

l D Co J

Итак, для определения звукового давления необходимо знание распределения аэродинамической нагрузки / (г, г|э) по радиусу и уг­ лу. С этой целью рассмотрим силы, действующие на элемент ло­ пасти. Они могут 'быть представлены как тяга Т и сила сопротив­ ления вращению Q. Интегрируя эти силы по радиусу (размаху) ло­ пасти при каждом значении г|э, а затем по тр, получаем

2г. г0

\ / i ( r > fldrdf,

11.

}

(3-55)

г ^ \ \ М г Л ) 4 Ы Ь

 

Оо

 

 

где Го — радиус лопасти.

 

 

Функции /2(г, г|?) и fi(r, \J?) представим в виде

 

/а(Л 40= W

( г ) Н (ф); |

(3. 56)

 

 

f\ { r ,^ ) = Q aH{r)H{^). j

Для вертолета на режиме висения, а для воздушного винта, ра­ ботающего в статических условиях, распределение тяги и сопро­ тивления вращению по радиусу достигают максимального значения в каком-то сечении лопасти и падают до нуля у конца лопасти (при г= /о). Зависимость для Н(г) представим следующим образом:

(3. 57)

На рис. 3.7 дается зависимость распределения Н(г) от тр Следует отметить, что для случая нагрузки, сосредоточенной на радиусе г= = 0,8/‘о, величина ri = 8,47.

Функция Н (ф) может быть представлена в виде

 

Н (4») = [1 + ki sin Щ+ kj cos уф].

(3. 58)

Первый член — единица соответствует постоянной составляю­ щей нагрузки неизменной при вращении лопасти. Величины Qa и То определяются путем использования уравнений (3.55) и (3.56):

H(r)

Рис. 3.7. Функция Я (г)

Q o = ~

П

f J Н (г)Л(ф) drdii

о о

(3.59)

т -Z.

1

о——

2%^^

 

"

I f J Н (г) М (<Юdrdty

 

 

[1

 

 

о

(3.60)

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

Рт&, у, г, t)-- =_ R _

r e' " ‘ ( ' - t ) ( l

s i n f u l

X

2n2D2mb

I

 

JJ

\ 2 r j

 

 

imnQ

v J

-ШЧ

 

x e

со

■}H

X

D

 

 

 

 

 

x [1 + k t sin ty + kj cos m

A +

T’oj ( ^ ) 2- ( - ^ ) 71] X

X [1+^,- sin /ф+AyCos jty\'E2}rdrdty,

(3.61)

где £i = cos crsin (ф—0); £ 2=sin <т, так как тяга действует только по оси Z, а сила сопротивления вращению (крутящего момента) — по осям X и У.

Если пренебречь зависимостью нагрузки от азимута, тогда вы­ ражения (3.55) принимают вид

Т = п j f 2{r)dr\

О

Q=n f fi{r)dr.

о

Рассмотрим частный случай, когда вся нагрузка сосредоточена на эффективном радиусе г*, равном 0,8 радиуса лопасти. Кроме то­

го, заменим

sin [

^ - j на его аргумент

.

Точку

наблюдения

возьмем при

0=0

и в дальнем поле, т. е. r<^D; —---- ►0. В этом

случае при определении амплитуды примем

D ~ R f а

при учете

фазы D ~ R —г cos a cos ф. Используя эти соотношения, запишем уравнение (3.54) в следующем виде:

 

 

mnQ

lmn2

(t —) 2те

 

 

 

;p=i

--------- e

 

'

c°/

^ (Q sin ф cos о -f-

 

 

 

2jt2Rc0

 

 

 

о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I m n Q r t

cos a

созф —Iqty

 

 

 

 

 

 

 

_i_____

 

 

 

 

+ rs in o )e

 

dty.

 

(3. 62)

Используем соотношения

 

 

 

 

 

 

 

 

j

е1хсоз* - (ч* dty=2niq J q{x)\

 

(3.63)

 

2it

 

 

 

 

 

 

 

 

 

^ e/xcos+- /9+ sjn w

=

_ _ £ 2я*% (*).

 

(3.64)

В нашем случае

 

x = rnnZricosg

 

 

 

 

 

 

 

 

 

 

 

 

Со

 

 

 

 

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

 

 

.<7?+1i mnQ

,tmnSax I(l‘ - *

) г

.

п соЧ Л г (ntnQri п

\

Р т п,-= 1

------ е

'

с« 1 \ Т

Sin

а—О - - ■ - ] / „ { ------- ‘ COS Я .

 

2nRc0

 

 

L

 

Qrimn] q \

с0

)

 

 

 

 

 

 

 

 

 

(3. 65)

Учитывая, что в нашем случае q=mn и вводя обозначение для крутящего момента Мкр= Qn, получим выражение для максималь­ ного звукового давления /теп-й гармоники в дальнем звуковом поле на расстоянии R от оси ротора

Ртп

RCQ

£о_1

COSO

(3.66)

 

Q r]\

«о

 

Это уравнение впервые было получено в 1936 г. Л. Я. Гутйным[11]. Представим выражение (3.65) в другом виде, вводя обозначе­ ние Mi=Qri/c0 и выразив тягу и крутящий момент через подъем­

ную силу

T = LTCOS р, M KP= L Trt sin р.

В результате будем иметь

Рт= 2лЛсо (C°S ^ Sln ° — Sin

J ^ mtlMi C0Sa^ (3*66')

3.4.2. Ротор при поступательном движении

Рассмотрим практически весьма важный случай звукового поля ротора, движущегося в направлении оси вращения со скоростью У0 (самолетный винт в установившемся полете). Точку наблюде­ ния возьмем в плоскости X, Z, т. е. У=0 и 0=0 (см. рис. 3.5). Для «определения звукового давления при действии на среду от эле-

ментов лопасти сосредоточенной силы Fi воспользуемся выражени­ ем [71]

imnQ

Fmi^___ rdrdty,

D

о о

где ^ = ^ « + ^ + ^ « 1 D = Y Z 2-j-p2(x22xr cos.Ф+ r2) - рас­ стояние между точками наблюдения (л:, z) и излучения с учетом

скорости полета p = -j/ 1 — Mo; а= М° ^ + D; М0= —

Приведенное йыражение удобно представить в .виде двух сла­ гаемых

г0 2я

imriS i t — —

 

\ с0 *)

rdrd%

ОО

То2тс

ОО

D

rdrdtyy

где /?т и pQ — величины звукового давления, обусловленные соот­ ветственно тягой и силой сопротивления вращению. Составляющие силы Fi аналогично (3.43) можно записать

dFmz= Ате тп

rdrdf,

(3. 67а)

dFmx= Bmsin фе/т*<2/--Ф> rdrdty,

(3. 676}

dFту— —Втcos

rdrdty,

(3. 67в)

где Ат, Вт — функции, определяющие распределение аэродинами­ ческих сил по лопасти.

Полное звуковое давление благодаря тяге лопасти винта с уче­ том выражения (3.67а) будет

. гпп2\

 

'~ ~ Б ~ ) ЫпЩ.

(3.68)

Переходя к полярным координатам в. плоскости вращения вин­ та x = r cos\J), y = r sin гр и учитывая выражения (3.67 6) и (3.67 в), получим выражение для полного звукового давления от действия на среду силы сопротивления вращению

То2*

.

тп2

~ 1

а

pQ= —

\ \ S

е'ш (2/" ф) (sin ф ——

cos ф — \ —----------- rdrdф.

J J

\

д х

д у I

D

 

О*0