Аэродинамические источники шума
..pdfТаким образом, пульсации подъемной силы на данной частоте со определяются возмущениями скорости с волновым параметром kx= —ы/и0. Взаимная спектральная (плотность пульсаций подъем ной силы может быть определена из выражения
о |
/ |
|
% 1 , |
AY (jfi, ш) Д Y* (Х2, w) |
, |
(3.33) |
||
^Ayar(-yi,-y2,a ))= lim n ^ v |
------ |
^ |
||||||
|
л |
|
|
|
|
|
|
|
где АК (лгх, со)= —?—Г |
ЬУ (x,t)Q-Mdt. |
|
|
|
|
|||
2л |
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. Функция Я (г)
2я
Q o = ~
П
f J Н (г)Л(ф) drdii
о о
(3.59)
т -Z. |
2я |
|
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 от оси ротора
Ртп |
2я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(x2— 2xr 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ф. |
||
4я |
J J |
\ |
д х |
д у I |
D |
|
О*0 |
|
|
|
|