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

книги из ГПНТБ / Хорошев Г.А. Шум судовых систем вентиляции и кондиционирования воздуха

.pdf
Скачиваний:
17
Добавлен:
23.10.2023
Размер:
7.65 Mб
Скачать

Используя эти значения в уравнении (48)

и принимая щ по­

стоянной по хорде, получаем

 

W = Эбяс8 і 1

(50)

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

§17. Влияние начальной турбулентности на входе

врабочее колесо на шум вентилятора

Одним из наиболее интенсивных источников аэродинамического шума лопаточных аппаратов является турбулентность потока на входе в решетку профилей. Как правило, судовые системы венти­ ляции и кондиционирования воздуха проектируют так, чтобы по­ ток на всасывании в рабочее колесо был равномерным, а степень его турбулентности мала. Но в некоторых случаях этого достиг­ нуть не удается (например, в случае, когда вентилятор располо­ жен после поворота всасывающего воздухопровода). Тогда турбу­ лентность потока на входе в колесо может иметь повышенные значения (5—10%), и уровень шума вентилятора соответственно также будет высоким. Возникновение интенсивного шума в этом случае обусловлено пульсациями давления на лопатке вследствие пульсаций подъемной силы [96].

Пульсации подъемной силы в свою очередь вызваны мгновен­ ным изменением угла атаки потока на лопатку. В самом деле, предположим, что к входной кромке лопатки пришел вихрь, ори­ ентированный так, что вектор мгновенного значения скорости сме­ стился и угол атаки стал отрицательным. В следующий момент по­ дойдет вихрь, вращающийся в противоположную сторону, и угол атаки станет положительным. Естественно, что эти изменения угла атаки приведут к изменениям подъемной силы профиля, пульса­ циям давления и, следовательно, к генерированию шума.

Для определения акустической мощности, излучаемой лопа­ точным аппаратом вентилятора вследствие начальной турбулент­ ности потока, воспользуемся зависимостью (34). Разницу пульсационного давления можно рассматривать как пульсации локаль­ ной подъемной силы, на единицу площади и выражать в виде коэффициента местной, подъемной силы:

р(х, у, т') = у-ри? (X, у) Су (х,

у, т'),

(51)

где і'і — местная средняя скорость потока,

параллельного

поверх­

ности.

. . .

 

61

Если предположить, что любая пульсационная составляющая скорость меньше осредненной скорости, то

 

(X, у , т') = - L 9ѵ\ (X,

у)

(X, у , т').

(52)

 

от

2

 

от

 

С учетом

последнего

формулу

(34)

можно записать

в виде

W

Р

 

 

 

(53)

48лс3

Данное интегрирование относится только к одной поверхности лопатки.

В турбулентном потоке корреляционная площадь S зависит главным образом от размера вихрей набегающего турбулентного потока. Она будет изменяться по поверхности только с изменением масштаба турбулентности. Если хорда пластины достаточно мала, то при допущении постоянства скорости Ѵ\ по хорде выражение (53) преобразуется к виду

где

(х, т')1 \

определяет средний интервал по у.

 

дт

J / у

 

Учитывая нормальные составляющие пульсационнои скорости по отношению к поверхности и принимая ее постоянной по хорде, получаем

< с и > и = Ф-

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

Тогда

 

 

(55)

С учетом (55) зависимость (54) примет вид

 

W ■

Ф*Ьѵ\ (х) S X, дѵ'

(56)

48лс3

дт

 

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

62

часть энергии спектра пульсации скорости турбулентного потока содержится в полосе центральной частоты f.

Среднеквадратичное значение пульсации скорости определится зависимостью

= CÜ2(ü')a,

Если обозначить корреляционные радиусы в направлении по­ тока и поперек потока через т\ и г2, то корреляционная площадь должна быть порядка ягіг2, т. е. должна быть равна площади эл-

Рис. 31. Уровни шума, излучаемого плоской пластинкой.

/ — шум, измеренный при */£>=8;

2 — расчет шума от нате­

кающей турбулентности при */£)=8;

3 — шум,

измеренный

при x / D 1; 4 — расчет шума от

срыва

вихрей

при x / D —\\

5 —расчет шума от натекающей

турбулентности

при x,'D= 1.

липса с полуосями гх и г2 . Исследования, проведенные Е. В. Вла­ совым и другими авторами, показали, что в турбулентных струях продольный размер вихря примерно в три раза больше попереч­ ного. С учетом этого можно принять, что S составляет величину порядка /у2.

Характерная частота конвекционной турбулентности может быть выражена как ѵ/Іи где 1\ — типичная «длина волны» вихрей. Исследования турбулентных струй свидетельствуют о том, что 1\ больше, чем гх. Как показали исследования, эффективное расстоя­ ние между вихрями одинаковой ориентации, обладающими боль­ шой энергией, оказывается примерно в 2я раз больше корреляци­ онного в направлении потока радиуса вихря. В этом случае Sco2 приобретет значение порядка щ2, и акустическая мощность, излу­ чаемая решеткой z профилей, согласно уравнению (56) будет равна

W = - ä - Г ФѢѵ\ (v')2 dx..

(57)

48яс3 "

63

Экспериментальная проверка этой зависимости была осущест­ влена И. Шарландом [96] для двух случаев: изолированной пла­ стины, обтекаемой потоком воздуха, вытекающим из круглого от­ верстия, и ротора осевого вентилятора.

Пластину размером 25 мм по хорде и 63 мм в поперечном на­ правлении помещали в воздушную струю, выходящую из форсунки диаметром 25 мм. На рис. 31 показаны уровни воздушного шума, возникающего при обтекании данной пластины. Кривая 3 соответ­ ствует случаю, когда пластина размещалась около выхода струи, так что большая ее часть находилась в среде относительно спокой-

Рис. 32. Зависимость уровня шума осевого вен­ тилятора от угла закрытия жалюзи.

ного потока. Кривая 1 соответствует случаю, когда пластина уста­ навливалась в такой точке, где поток по всей поверхности пластины был турбулентным. Приведенные кривые наглядно показывают, что излучение шума значительно (на 10—12 дБ) увеличивается при наличии турбулентности в потоке впереди пластины. Оценка шума, создаваемого пластиной в полностью турбулентном потоке

(x/D = 8),

произведенная Шарландом [96] на основании уравне­

ния (57),

показывает вполне удовлетворительное совпадение с эк­

спериментом. Турбулентность в вентиляторе создавали путем по­ мещения непосредственно в раструбе заборной трубы изогнутого в кольцо круглого стержня диаметром 0,9 см. Диаметр кольца был равен осредненному диаметру рабочего колеса и примерно в 20 раз больше диаметра самого стержня. В этом случае расчетное и эк­ спериментальное значения уровней шума вентилятора также имели удовлетворительное совпадение.

В работе [80] описано влияние жалюзи на шум осевого венти­ лятора. При изменении угла наклона жалюзи изменялась турбу­ лентность потока на входе в рабочее колесо вентилятора. На

64

рис. 32 показано изменение общего уровня воздушного шума вен­ тилятора в зависимости от угла наклона жалюзи. Как видно из рисунка, в этом случае не наблюдается такого резкого увеличе­ ния уровня шума, хотя характеристики турбулентности должны были бы изменяться значительно (уровень шума увеличивается только на 5—8 дБ). Это, очевидно, объясняется изменением импе­ данса выходного сечения вентилятора при уменьшении живого се­ чения на входе в вентилятор. Однако качественная картина изме­ нения уровня шума вентилятора от турбулентности натекающего потока полностью подтверждается и в этом случае.

§ 18. Неоднородность структуры потока за рабочим колесом

Рассмотрим, как влияет на возникновение шума в судовых вен­ тиляторах обтекание неоднородным потоком лопаточного аппа­ рата или какого-либо тела, находящегося за лопаточным аппара­ том (например, языка улитки центробежного вентилятора, растяжки двигателя осевого прямоточного вентилятора). Под неод­ нородностью потока понимается различие в скоростном и стати­ ческом давлении в аэродинамическом следе за каждой лопаткой и в области основного потока. Шум, возникающий вследствие воз­ действия неоднородного потока на обтекаемое тело, проявляется, как было показано в § 2 и 3, на частотах

е __

km

' г

6(Г'

При выводе зависимости интенсивности звукового давления на частоте f примем следующие допущения:

— форма входной кромки обтекаемого тела — прямоугольная и симметричная относительно вектора скорости набегающего по­ тока;

звуковой импульс, возникающий в районе обтекаемого тела, распространяется без искажения через межлопаточные каналы рабочего колеса и всасывающий воздухопровод;

корпус вентилятора абсолютно жесткий;

переменная аэродинамическая сила, действующая на тело, изменяется по такому же закону, как и скорость потока на выходе из решетки лопаток.

При выводе уравнения для спектральной составляющей звуко­ вого давления рассмотрим для простоты случай, когда решетка профилей неподвижна, а обтекаемое тело, находящееся за решет­ кой, движется вдоль оси решетки со скоростью и, равной скорости движения решетки (рис. 33). Когда движущееся тело находится в середине межлопаточного канала, на него действует аэродина­

мическая сила, пропорциональная полному давлению потока Н =

 

О

Рс

рѵі . При попадании тела в зону аэродинамического следа

3 Г. А. Хорошев и другие

65

эта сила за короткий промежуток времени изменяется, в резуль­ тате чего тело будет воздействовать на окружающую среду с пе­ ременной силой, которая по своим свойствам является акустиче­ ским источником звука осциллирующего типа [72]. Если обтекае-

Рис. 33. Схема обтекания тела неоднородным потоком.

мым телом является лопатка рабочего колеса осевого вентиля­ тора, выполненного по схеме НА + РК, наряду с появлением пульсационной аэродинамической силы на входной кромке лопатки

 

появляется

еще

пульсационная

 

аэродинамическая

 

сила

вслед­

 

ствие

флуктуаций

 

подъемной

 

силы профиля.

 

 

 

 

 

Генерирование шума от неод­

 

нородности

потока

из-за пульса­

 

ции

подъемной

силы

профиля

 

в данном

параграфе не рассмат­

 

ривается. Этот вопрос хорошо

 

освещен Е. Я- Юдиным

[72] и

 

другими авторами

[93].

Кроме

Рис. 34. Проекция входной кромки

того,

также не рассматривается

картина возникновения вихревого

обтекаемого тела, наклоненного на

шума

вследствие-

обтекания ка­

угол Ѳ в системе координат хуг.

кого-либо тела неоднородным по­

 

не отличается от физической

током, так как она принципиально

картины

генерирования

вихревого

шума при обтекании лопаток турбулентным потоком.

 

 

 

Итак, рассмотрим возникновение шума

вследствие

обтекания

неоднородным потоком какого-либо тела [73]. Совместим проек­ цию входной кромки этого тела на плоскость, нормальную вектору скорости набегающего потока с плоскостью хоу (рис. 34). Поло­ жительное направление оси z примем совпадающим с направле­ нием, обратным движению потока. Обозначим через AF (%) пуль­ сирующую аэродинамическую силу, действующую на единичный

66

элемент входной

кромки и изменяющуюся во времени по закону

£ = £(т). Так как

аэродинамический след изменяется с увеличе­

нием расстояния от выходных кромок лопаток, то эта сила также зависит от расстояния As обтекаемого тела до выходных кромок лопаток. Положим, что эта сила в зависимости от расстояния As изменяется по закону X(As). Тогда ее можно представить в виде

AF (т) = AF% (т) X (As).

(58)

В свою очередь этот элемент будет воздействовать на окружаю­ щую среду с такой же по величине, но противоположно направлен­ ной силой. В элементе dS на среду действует сила AFl(x)X(As)dS. Очевидно, что

J AFI (т) X(As) dS = F (т) X (As).

(59)

Для элемента dS потенциал скоростей, обусловленный сосредото­ ченной силой, будет равен [12]

Ф =

і

Zi

_д_ e-iks J_

(60)

4яр£]С

dz

 

 

 

где Zi = AF%(x)X(As)dS, ki — волновое число;

г — расстояние до точки замера; р — плотность окружающей среды.

Потенциал скоростей, обусловленный всей системой сил, дей­

ствующих на обтекаемое тело, выражается уравнением

(61)

С учетом того,

что

e-iks

 

 

 

 

 

 

д_

 

d

 

(62)

 

 

 

dz

г

I

dr

 

 

 

 

 

 

где

у — угол

между

осью

z и

радиусом-вектором г,

выражение

(61)

можно представить в виде

 

 

 

 

Ф =

і

fj AFUx)X(As) _a_

cos ydS.

(63)

 

 

 

Anpkxcs

 

dr

 

 

Дифференцируя функцию, стоящую под интегралом, и произ­

водя некоторое упрощение, получаем

 

 

 

д_

 

 

 

 

 

 

(64)

 

dr

 

 

 

 

 

 

 

Звуковое

давление, обусловленное

потенциалом

скоростей,

' равно

 

 

 

 

 

 

 

Р = Р

дер

(65)

 

д т

 

 

 

з*

67

Подставив сюда выражение для потенциала (63) с учетом

(64)

и (59), получим для амплитуды звукового давления

 

P==l b " k lF (t) M As)] cos;Y,

(66)

где /г( т ) — аэродинамическая сила, действующая на обтекаемое тело, находящееся на некотором расстоянии от выходных кромок лопаток.

В зависимости от конструкции рабочих колес вентилятора сила F (х) может быть периодической (при постоянном шаге лопаток) и

Рис. 35. Трапецеидальная эпюра скорости и положение наклонного тела в различные моменты времени.

непериодической (при переменном шаге лопаток).

В случае равно-

бо

(z — число ло­

мерного шага сила F(x) имеет период Т ==----

паток колеса) и может быть разложена в ряд Фурье:

 

 

00

 

 

 

 

 

 

F (т) = N0-f- 2

A/к cos ( k m —ф),

(67)

где

 

k = \

_____

 

 

 

 

 

 

 

 

 

 

NK= ~\fЯк+ tf;

 

 

 

 

t g

^ ;

 

 

 

T

 

 

г

 

2

г

F (t ) cos kmdx\

bK=

2

c

(68)

aK— J

— J F (x) sin kmdx.

T o

 

 

T o

 

Тогда выражение для звукового давления может быть запи­ сано в виде

4пег дх

NKcos ( k m —ty) cos у-

(69)

k = \

 

68

Как видно из выражения (69), для расчета звукового давления необходимо знать закономерность изменения силы F(т) или фор­ му аэродинамического следа по шагу на выходе из решеток про­ филей.

Рис. 36. Эпюра скорости на выходе из рабочего колеса по шагу лопаток.

/ — действительная

эпюра скорости;

2 — схематизированная

(прямоугольная) эпюра скорости.

Экспериментальное

исследование

формы аэродинамического

следа на' выходе из рабочих колес судовых вентиляторов с по­

мощью термоанемометра

фирмы

Диза

(см. § 24)

показало, что

все эпюры скорости можно в основном

 

 

свести к двум видам:

эпюра

ско­

 

 

— трапецеидальная

 

 

рости (рис. 35), наблюдаемая у

цен­

 

 

тробежных вентиляторов с радиаль­

 

 

ными лопатками и у осевых вентиля­

 

 

торов на расчетном режиме, а также

 

 

на режимах, меньше номинальных;

 

 

 

— прямоугольная эпюра

скорости

Рис. 37. Переменная сила, дей­

(рис. 36), встречающаяся

у центро­

ствующая на

обтекаемое тело.

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

при производительности номинальной и,больше номинальной. Определим коэффициенты Фурье для случая трапецеидальной

формы эпюры скорости на выходе из решетки лопаток, наиболее часто встречающегося у судовых вентиляторов. Развернем на пло­ скость эпюру скорости за рабочим колесом вентилятора (см. рис. 35). Там же показано взаимное расположение обтекаемого тела, наклоненного на угол Ѳ относительно эпюры в различные моменты времени хп- Для наглядности и простоты дальнейших выводов будем считать, что форма входной кромки обтекаемого тела — прямоугольная. При перемещении обтекаемого тела из по­ ложения 0 в положение 3 сила, действующая на него, изменяется приближенно по тому же закону, что и Дп (рис. 37). Величины хп определяются следующими зависимостями:

69

4 =

А Р A B — d/cos6

_

Ахг -f- Дл:3 + А.х3b tg Ѳ — d/cos Ѳ .

 

 

aR

 

2aR

 

 

 

 

 

т2— %\ —

DE -f- d/c o s Ѳ

b tg Ѳ + d /c o s Ѳ — Ax3

(70)

aR

 

aR

 

 

 

 

 

Т —т2 -Axt -f- Ax2Ax3

b tg Ѳ — d /c o s 0 .

 

 

 

 

2aR

 

 

где а — угловая скорость вращения рабочего колеса.

Сила F (х) изменяется во времени по следующему закону:

при

 

о

F(x) = F ± i

 

 

і

 

 

 

 

 

Т1

при

 

т2—

F (т) = F = const;

при

 

Т — т2

F(x) = p l = l .

 

 

 

 

 

Т — т3

Подставив (71)

в (68), найдем

 

 

 

г

 

 

 

ак =

j* F (т) cos kmdx =

F — cos kmdx +

 

 

о

 

 

Ti

 

 

 

 

 

+ J F cos kmdx -f- J F

— — coskaxdx

 

 

 

 

r,

2

 

2F

(COS k(£)Xi -j- COS &CBT2—cos kwT — 1).

 

TkWxt

Аналогично для bKполучим

 

 

,

 

2F

(sin kcoXi + sin Ы х г—sin ka>T).

 

Tk2^

 

 

 

 

Тогда

 

 

 

 

 

\j

1 f

2

, «2

2F

. 6шт, . k(äX3

 

 

 

+

 

 

(71)

(72)

(73)

(74)

Подставляя выражение (74) в формулу для амплитуды звуко­

вого давления (69) и дифференцируя,

получаем

выражение для

звукового давления на частоте fz:

 

/

 

р = 4,15- 101—7 F -sin

k(äX1 sin- &С0Т,

(75)

rkx1

2

2

 

В зависимости (75) неизвестной величиной остается пульсаци онная аэродинамическая сила

F = Fi - F 2,

где F1 — полная аэродинамическая сила, действующая на обте­ каемое тело, когда оно находится в центре межлопаточного канала;

70

Соседние файлы в папке книги из ГПНТБ