УДК: 534.121.1 OECD: 1.03.AA
Акустическое излучение пластины в поле действия турбулентных
пульсаций давления
Светлов В.В.
К.т.н., ФАУ «ЦАГИ», ул. Радио 17, 105005, г. Москва, РФ
Аннотация
В статье описан алгоритм расчета акустического излучения пластины, находящейся в поле действия пристеночных турбулентных пульсаций давления применительно к излучению панелей фюзеляжа в салон самолета. Используемый подход основан на разложении колебаний пластины по базисным функциям. В качестве воздействующей силы используется нагрузка, которая задается моделью, предложенной Ефимцовым Б.М. Вибрационные поля рассчитываются на основе дифференциальных уравнений колебаний пластины. Интенсивность прошедшего за пластину звука определяется посредством численного интегрирования в пределах акустической области волновых чисел. Приведено сравнение полученных результатов с известными асимптотическими соотношениями.
Ключевые слова: неподкрепленная пластина, турбулентный пограничный слой, виброускорение, интенсивность акустического излучения.
Noise generation of the plate in the field of turbulent pressure pulsations
Svetlov V.V.
PhD, Central Aerodynamics Institute, Moscow Branch, 17, Radio str., 105005, Moscow, Russia
Abstract
The article describes an algorithm for calculating the acoustic radiation of a plate located in the wall pressure structuration field beneath a turbulent boundary layer in relation to the radiation of sound into the aircraft cabin. The approach used is based on the decomposition of plate vibrations by basic functions. The load is given by the Efimtsov B.M. model. Vibration fields are calculated on the basis of differential equations of plate vibrations. The intensity of the sound passing behind the plate is determined by numerical integration within the acoustic region of the wave numbers. The obtained results are compared with the known asymptomatic relations.
Keywords: unsupported plate, turbulent boundary layer, vibration acceleration, acoustic radiation
intensity.
Введение
Среди общих проблем защиты от негативного акустического воздействия важное место занимают вопросы, связанные с авиационным шумом. Авиационный шум отличается от других источников, таких как автомобильный и железнодорожный транспорт, инженерно-технологическое оборудование предприятий, проведение
E-mail: svetlov-valeriy@yandex.ru (Светлов В.В.)
строительных работ тем, что имеет высокие уровни звука. Основными источниками шума, создаваемого воздушными судами (ВС) на местности, являются силовая установка и планер с элементами его механизации, а в салоне - силовая установка, источники, связанные с обтеканием поверхности фюзеляжа ВС и системы внутреннего бортового оборудования.
Отличительной чертой основных источников шума в салоне ВС является их аэродинамическая природа возникновения. Они имеют различные пути и механизм воздействия, что приводит к необходимости использования разных подходов к снижению негативного акустического воздействия. В зависимости от вида силовой установки, в салоне ВС может формироваться поле с разным спектральным составом. В случае двухконтурных турбореактивных двигателей преимущественно формируется поле со сплошным спектром в широком диапазоне частот, при использовании турбовинтовой силовой установки - сплошным спектром с ярко выраженными узкополосными составляющими. Механическое воздействие силовой установки на пилон сопровождается появлением структурного шума. Внутреннее бортовое оборудование, такое как система вентиляции и кондиционирования воздуха - носит аэродинамический характер воздействия. Исследования показывают, что при высоких дозвуковых скоростях (числах Маха более 0,4) значимым, по суммарному вкладу в общий уровень, источником становится шум, вызванный колебаниями упругой поверхности фюзеляжа под действием турбулентного пограничного слоя (ТПС) [1].
Существующие методы расчета вибрации и акустического излучения пластин фюзеляжа ВС под воздействием турбулентного пограничного слоя [2, 3] основаны на различных подходах, каждый из которых имеет свои сильные и слабые стороны. К недостаткам имеющихся методов можно отнести: ограниченность частотного диапазона расчета, использование в качестве геометрического параметра только толщины пластины, пренебрегая ее габаритными размерами, отсутствие точных исходных данных, которые зачастую получаются экспериментальным путем и только для конкретных конструкций и пр. Использование же численных методов, в частности метода конечных элементов, позволяет с высокой точностью определить искомые виброакустические параметры, однако требует такого числа конечных элементов, что практическая реализация требует существенных временных затрат, что становится серьезной проблемой при его использовании.
1. Перемещения неподкреплённой пластины
Рассмотрим подход, основанный на разложении колебаний пластины по базисным функциям [4]. В общем случае, уравнение, описывающее колебания тонкой пластины (рисунок 1) для изгибных волн, имеет вид:
d 2w
D(1 + zn)A2w + ph—- = Q(x,y,t) (1)
ot2
где D - цилиндрическая жесткость пластины, n - коэффициент потерь материала пластины, А - оператор Лапласа в декартовой системе координат, w - нормальное смещение пластины, p - плотность материала пластины, h - толщина пластины, t -время, Q - обобщенная сила, которая в данном случае описывает некоторое внешнее воздействие пристеночных турбулентных пульсаций давления, оказываемое на пластину.
Применив преобразование Фурье к уравнению (1) и разложив его по базисным функциям [5] fmim2 = sin ("¿р) sin ("¿J^) получим выражение для виброперемещений свободно опертой пластины:
[D(1 ++ k2 ) phuj ]Wmim2 /mim2 ^ ] Qmim2 fmim2
mim2 mi m2
где к\ = , к2 = , т1 и т2 - целые числа (1, 2, 3...), ш - круговая частота.
Рис. 1. Расчетная схема тонкой неподкрепленной пластины под действием ТПС
Правая часть в (2) является коэффициентами разложения обобщенной силы Q(x,y,ш) по базисным функциям. Используя ортогональность базисных функций, получим явное выражение для этих коэффициентов:
1 Г^1 Г^2
Qmlm2 = "Г / Q(x,У,ш)fml „2 йхйу (3)
V2 ./о
Нормировочная константа С = находится из выражения:
СJ (^„2 )2 йхс1у = j j sin2(klx)sin2(k2y)dxdy = —— (4)
После подстановки полученных выражений в (2) получим формулу общего решения для нахождения перемещений неподкреплённой пластины:
w ( ) _ \ л W f _ \ л _Qmim2 fmim2_
W (x,y,u)_ Wmim2/mim2 _ D(1 + in)(k2 + k|)2 - phuj2
mi m2 mi m2
2. Спектр виброперемещений пластины под действием ТПС
В уравнении (5) обобщенная сила имеет размерность давления, которое представляет собой силу на единицу поверхности. В случае воздействия на пластину пульсаций давления ТПС, нагрузка носит случайную природу, поэтому, в данном случае, невозможно найти такую детерминированную величину как виброперемещение. Корректно говорить о спектральной плотности этой величины.
Для нахождения спектральной плотности виброперемещений |Ш |2 необходимо Ш в уравнении (5) умножить на комплексно-сопряженную величину Ш*. Это действие вполне обоснованно при рассмотрении колебаний панели обшивки ВС, вызываемых безотрывным ТПС:
,о „ ^^ Qm^m2 Q„1 „ з1и2 (к^) вШ2 (к2у)
Ш 2 = ШШ* = > -т1„2-^^-V- (6)
1 1 „т [Я(1 + ¿п)(к2 + к2)2 - р^2][£(1 - ¿п)(к2 + к2)2 - рки2] (6)
Числитель в (6), в соответствие с положениями работы [4], также разложим по базисным функциям и представим в виде:
2 2
пЬ 1 пЬ2 1*1*1 гЬ2 (7)
= С2 / / / / ф(х,у,х',у',ш) вШ^^х) вт2(к2у) 81п2(к1Ж/) вт2(к2у')йхйуйх'йу' ио ио ио ио
где ф(х,у,х',у',ш) - спектр пространственных корреляций случайного поля турбулентных пульсаций давления, определяющий интенсивность и взаимную корреляцию спектральных составляющих в двух произвольных точках наблюдения (х,у) и (х',у').
Заменив для краткости записи знаменатель в (6) на ^т1т2, получим формулу для нахождения квадратичного значения спектра виброперемещений в любой точке пластины:
те ф
Бщ = ^ 7т1т2 81п2(к1х) в1п2(к2у) (8)
т1т2 1 2
Для нахождения среднеквадратичного значения виброперемещений воспользуемся теоремой Лагранжа о среднем значении:
_ 1 сЬ1 гЬ2 Ф 1 те ф
БЩ = ^ 81п2(к1х) 81п2(к2у)йхйу £ ф^ = 1 £ ф^ (9)
Ь1Ь2^ ^ тт ^т1т2 4 т1т2 ^т1т2
3. Модель ТПС
Существует большое число моделей частотно-волнового спектра ТПС, которые позволяют определить спектральную плотность пространственных корреляций, необходимую для вычисления спектра виброперемещений пластины. В качестве расчетной модели в данной работе выберем модель спектра пространственных корреляций случайного поля турбулентных пульсаций давления, предложенную Ефимцовым Б.М. [6].
В выбранной модели учитывается зависимость пространственной корреляции пульсаций давления от толщины пограничного слоя (8), что является удобным в рассматриваемой задаче. Согласно этой модели, спектральная плотность давления в двух различных точках (х,у) и (х',у') задается в следующем виде:
. . р2 (— \х — х'|\ /-|у — у'\\ ( — гш(х — х')\ , . ф(х,у,х ,у ,ш) = — ех^ех^) ех^ ^ ) (10)
где р2 — среднеквадратичное значение турбулентных пульсаций давления, Аш — ширина полосы круговой частоты, Л1 и Л2 — масштабы корреляции, ирь, — фазовая скорость.
Пространственные масштабы корреляции по соответствующим осям координат определяются соотношениями из [6].
Характеристики ТПС по модели Ефимцова Б.М. могут быть вычислены методом, изложенным в [7]. Этот метод позволяет по известному числу Маха (Мте), высоте полета и расстоянию от носа ВС вычислить спектр пульсаций давления в третьоктавной полосе частот.
Интеграл в выражении (7) для краткости записи заменим на обозначение 11. Данный интеграл рассчитывался численным методом.
После подстановки замены в (10) получаем итоговую формулу для нахождения коэффициентов разложения обобщенной силы (7) получим:
фт1т2 = С2 ■ ■ 11 (И)
Аш
4. Спектр виброускорения и интенсивности акустического излучения пластины под действием ТПС
Найдем связь между волновыми числами в распространяющейся звуковой волне, используя волновое уравнение для давления в системе координат, представленной на рисунке 1.
Обозначим скорость колебательных движений в воздухе за пластиной, сопровождающуюся переносом звуковой энергии через Уа. Учитывая малую амплитуду колебаний пластины, можно считать, что колебательная скорость (V) направлена по нормали к ее поверхности и равна колебательной скорости частиц воздуха, движение которых сопровождается излучением акустической энергии, т.е.
Vа(k х,ку|^=0 — V (к ку ,и) (12)
Используя линеаризованное уравнение Эйлера можно получить связь между колебательной скоростью частиц среды и звуковым давлением:
1 dp iup0 dz
Va(kx,ky ,u)\z=q = -г-— \z=0 (13)
где р0 - плотность воздуха за пластиной.
Представив звуковое давление в поле излучения пластины согласно [3] получим связь между волновыми числами:
к2 — ко2 - (кХ + кУ) (14)
где к0 - волновое число для воздуха.
Подставляя в (13) связь между волновыми числами и, учитывая, что др — гкгр(кх,ку) получаем:
п , ч ро^(кх,ку)ко
Р(кх,ку) —--, = (15)
\Jk 0 — (kX + Щ)
Выразим интенсивность акустического излучения (мощность на единицу поверхности) с использованием соотношения из [8]. Таким образом для спектральных величин получим окончательное соотношение:
о _ [ [ Р0иБУ(кхку) „ „
о/ — , ■ акхйку
«Ц<0 Чк2 - (кх + к2> (16)
Для того, чтобы воспользоваться формулой (16), необходимо представить Бу в волновом пространстве. Для этого воспользуемся теоремой Парсеваля и выполним преобразование Фурье, учитывая связь спектра виброперемещения (Бщ) со спектром виброскорости (Бу) в частотном пространстве:
те
1y
= T-y£ \W(x
4п2
mim2
pLi 1-L2
/ sin(kix)e-ikxXdx / sin(k2y)e-iky ydy 00
Вычислим соотношения из (17):
2
Вычислив интеграл в (17) и раскрыв модуль получим: Б = ^ Фт1ТО2 2к2[1 - (-1)т1 008(^1)] 2^2 [1 - (-1)т2 0С8(куЬ2)] (18)
Т 4п2 ¿2 ^ ■ к - кХ]2 ■ [к22 - к2]2 (18)
Подставив взаимосвязь между спектром виброскорости со спектром виброперемещения Бу = в (16), получим формулу для расчета спектра интенсивности:
1 = РР1Г Ц /; 2 ^2 ■ ^ ^^ (19)
k2 i k2<k2
(kX + kj/)
Для того, чтобы вычислить (19), удобно выполнить замену переменных в интеграле, чтобы интегрирование проходило по прямоугольной области. Замена переменных через kx = ho cos ф и ky = ho sin ф и вычисление якобиана такого преобразования позволяют записать выражение для интенсивности акустического излучения в следующем виде:
г р0ш3 1 fk° [2ж Sw(h0 cos ф,ho sin ф) I = —— -г—; -, ^йфй^ =
2 4П Л ^ у/к2 - (kX + k2)
= ро^! fk0 ^ Ф^т 2(mf)2[1 - (-1)m1 cos(ho сов(ф)Ьх)^ (20)
8П2 -/О Уо ¿m2 ' ykf-^0 ■ [(mf )2 - (ho COs^))2]2 X
2(mL^f )2[1 - (-l)™2 cos(ho сов(ф)Ь2)], ,.„
X [(^)2 - (ho sin^))2]2
5. Асимптотические соотношения для спектра виброускорения и интенсивности акустического излучения пластины под действием ТПС
Для проверки расчета, выполненного по формуле (9), воспользуемся асимптотическим соотношением для расчета спектра виброускорения упругой конструкции, находящейся под действием ТПС из [2]:
Ф .. (f) = Фд(f)F(в)n(f) (21)
(f) = 2(ph)2n (21)
где Фд - спектральная плотность пульсаций давления, n(f) - отношение плотности собственных частот упругой конструкции к плотности собственных частот пластины (в случае плоской пластины принимается равным 1), р - плотность материала конструкции, П - коэффициент потерь, F(в) - функция безразмерных параметров, определяемая из [2].
Для проверки расчета, выполненного по формуле (20), воспользуемся асимптотическим соотношением для расчета спектра интенсивности акустического излучения упругой конструкции, находящейся под действием ТПС из [3,9]:
Ф(.ч = Фд(f)F(в)n(f) Пи (22)
1 (f) 2phu n ( )
где все множители соответствуют (21), а Пи/n - отношение коэффициента потерь на излучение к коэффициенту потерь в пластине определятся экспериментально и может быть принят из [3].
6. Вычисление спектра виброускорения и интенсивности акустического излучения пластины под действием ТПС
Получив выражения (9) и (20) произведем вычисление спектров виброускорения и интенсивности акустического излучения пластины.
Первоначально необходимо задаться параметрами, характеризующими саму пластину. Примем длину пластины L = 1 м, ширину пластины L2 = 0.75 м, толщину h = 10-3 м, коэффициент потерь п = 0.01, материал - алюминий. Средой, находящейся за пластиной, примем воздух с параметрами c0 = 343 м/с и р0 = 1.2 кг/м3.
Расчеты для модели, предложенной Ефимцовым Б.М., произведем с помощью метода из [7], для этого примем, что пластина находится на расстоянии 30 м от носа ВС, полет происходит на высоте 9000 м с числом Маха равным 0.8.
Сравнение результатов расчета уровня виброускорения и уровня интенсивности акустического излучения пластины под действием ТПС для ш\ = 100 и m2 = 100 представлено на рисунках 2-3.
Частота. Г и
Рис. 2. Сравнение расчета 1/3 октавных уровней виброускорения тонкой неподкрепленной пластины под действием ТПС по формуле (9) с асимптотическим
соотношением по формуле (21)
-105 -'-'-'-'-'-'-'-'-'-
0 1000 2000 3000 4000 5000 6000 7000 SOOO 9000 10000
Частота, Гц
Рис. 3. Сравнение расчета 1/3 октавных уровней интенсивности акустического излучения тонкой неподкрепленной пластины под действием ТПС по формуле (20) с асимптотическим соотношением по формуле (22)
Влияние числа мод на спектры виброускорения и интенсивности акустического излучения пластины представлено на рисунках 4-5.
Рис. 4. Влияние числа мод на 1/3 октавные уровни виброускорения пластины,
формула (9)
-so
Ш
4 -BO
2 -90 Q
_J
О
^-100
л
ь
о
ш -110
£ О X
03
* -120 -130 -140
Рис. 5. Влияние числа мод на 1/3 октавные уровни интенсивности акустического
излучения пластины формула (20)
Заключение
По результатам сравнения данных по предложенному алгоритму расчета с асимптотическими выражениями для виброускорения и интенсивности акустического излучения тонкой неподкрепленной пластины, находящейся в поле действия пристеночных пульсаций давления турбулентного пограничного слоя, получена удовлетворительная сходимость.
В отличие от асимптотических соотношений, предложенный алгоритм учитывает влияние собственных форм колебаний пластины в области низких частот, а также влияние явления волнового совпадения в области высоких частот. Алгоритм обладает высокой скоростью расчета по сравнению с численными методами и позволяет при сравнительно небольшом числе мод, получать значения во всем слышимом диапазоне 1/3 октавных полос частот для наиболее характерных размеров пластин, применяемых в авиастроении.
Примененный подход может быть использован для решения задач по определению виброакустических полей как при других видах нагрузки, так и при иных геометрических формах пластины.
Работа выполнена при финансовой поддержки РНФ (грант 21-71-30016).
Список литературы
1. Ефимцов Б.М. Колебания обшивки фюзеляжа самолета, вызываемые турбулентным пограничным слоем // Учен. зап. ЦАГИ. 1987. Т. XVIII, №1 стр. 90-97.
2. Ефимцов Б.М., Колебания цилиндрической панели в поле турбулентных пульсаций давления // Акуст. журн. 1986. Т. 32. № 4. с. 536-538.
3. Авиационная акустика. Часть 2 / под ред. А.Г. Мунина. М.: Машиностроение, 1986. 264 с.
4. Martin N.C., Leehey P., Low wavenumber wall pressure measurements using a rectangular membrane as a spatial filter // Journal of Sound and Vibration. 1977, V. 52, N. 1. pp. 95-120.
5. Вибрации в технике. Том 1 / под ред. В.В. Болотина. М.: Машиностроение, 1978. 352 с.
6. Ефимцов Б.М., Характеристики поля пристеночных турбулентных пульсаций давления при больших числах Рейнольдса // Акуст. журн. 1982. Т. 28. № 4. с. 491-497.
7. David Alan Bies, A review of flight and wind tunnel measurements of boundary layer pressure fluctuations and induced structural response". NASA CR-626, 1966. p. 94.
8. Заборов В.И. Теория звукоизоляции ограждающих конструкций. Издание 2-е. М.: Стройиздат, 1969. 184 с.
9. Ефимцов Б.М. Применение энергетического статистического метода для оценки акустического излучения пластин при псевдозвуковом нагружении // Тр. ЦАГИ. 1978. Вып. 1902 с. 3-8.
References
1. Efimtsov B.M. Vibration of the aircraft fuselage skin caused by a turbulent boundary layer // Scientific notes of TsAGI. 1987. V. XVIII. № 1. P. 90-97.
2. Efimtsov B.M. Cylindrical panel vibrations in a turbulent pressure pulsation field // Akustichesky zhurnal, 1986. V. 32. № 4. P. 536-538.
3. Aviation acoustics. Path 2 / edited by A.G. Munin. M.:Mashinostroenie, 1986. P.
264.
4. Martin N.C., Leehey P., Low wavenumber wall pressure measurements using a rectangular membrane as a spatial filter // Journal of Sound and Vibration. 1977, V. 52, N. 1. P. 95-120.
5. Vibrations in technic. Volume 1 / edited by V.V. Bolotin. M.:Mashinostroenie, 1978. P. 352.
6. Efimtsov B.M. Characteristics of the wall turbulent pressure fluctuation field at high Reynolds numbers // Akustichesky zhurnal, 1982. V. 28. № 4. P. 491-497.
7. David Alan Bies, A review of flight and wind tunnel measurements of boundary layer pressure fluctuations and induced structural response". NASA CR-626, 1966. P. 94.
8. Zaborov V.I. Theory of sound isolation of enclosing constructions. Edition 2. M:Stroyizdat. 1969. P. 184.
9. Efimtsov B.M. Application of the energetic statistical method for evaluation of the acoustic radiation of plates under pseudosound loading // Proceedings of TsAGI. 1978. V. 1902. P. 3-8.