Проблемы машиностроения и надежности машин, 2022, № 2, стр. 68-79
Оценка изменения звукоизоляционных характеристик силовых панелей при их высокоинтенсивном акустическом нагружении
П. А. Попов 1, *, А. А. Иголкин 1, Е. В. Шахматов 1
1 Самарский национальный исследовательский университет им. С.П. Королева
Самара, Россия
* E-mail: banduir@rambler.ru
Поступила в редакцию 18.05.2021
После доработки 09.11.2021
Принята к публикации 20.12.2021
- EDN: TZEGTD
- DOI: 10.31857/S0235711922020122
Аннотация
В статье представлен подход, основанный на анализе акустических нагрузок, который позволяет по полученной из эксперимента звукоизоляции, найти ее параметры при более жестком случае нагружения. В результате получены данные, которые подтверждают изменение звукоизоляции конструкции разных силовых панелей при передаче внешнего акустического давления под экспериментальную панель, в то же время имеются общие закономерности параметров перепада акустического давления. Выявленные закономерности позволили провести анализ звукоизоляции трех экспериментальных образцов при их эксплуатационном акустическом нагружении.
В технической акустике возникают задачи оценки действия высокоинтенсивных пульсаций давления на различных этапах эксплуатации изделий. Примером могут служить аэроакустическое нагружение современных боевых самолетов, при котором амплитудные уровни нагрузок достигают 173 дБ [1], нагружения внешней части конструкции сборочно-защитных блоков ракет-носителей на старте и выведении, при которых суммарные среднеквадратические уровни акустического давления могут достигать величин свыше 150 дБ и 160 дБ соответственно [2]. Как показала практика, в широком диапазоне спектра от 18 Гц до 4470 Гц высокоинтенсивные уровни акустических нагрузок труднореализуемы в лабораторных условиях для некоторых акустических реверберационных стендов (АРС). В результате научно-информационного поиска открытых источников было обнаружено, что исследования в области определения реакции конструкций на широкополосные случайные воздействия акустического типа проводились в работах [3, 4], в которых решались задачи оценки звукоизоляции (ЗИ) и акустических вибраций летательных аппаратов расчетными и экспериментальными методами. Исходя из технических требований, выбирались те или иные способы защиты изделий внутри отсека от разрушающего действия шума путем размещения звукопоглощающих и вибродемпфирующих материалов. При этом коэффициенты безопасности конструкции при проектировании изделий, обоснованные вероятностно-статистическими методами, приведены, например, в работе [8]. В работах [5–7] представлены результаты исследований, содержащие обширную информацию о характеристиках обечаек и звукоизолирующих материалов, таких как индекс ЗИ, коэффициенты звукопоглощения, механических потерь, постоянных затухания и постоянных распространения для отсеков. Однако в этих работах предполагалось, что характеристики ЗИ являются инвариантными для различных внешних нагрузок. В работе [9] рассматривались вопросы использования метода наименьших квадратов (МНК) для уточнения звукоизоляции пластин из стали, дюралюминия и фанеры различных толщин, рассчитываемых по “закону массы”, однако представленные данные могут быть верны только для амплитудных характеристик внешнего акустического давления, которое использовалось в экспериментах.
В настоящее время анализ изменения звукоизоляционных характеристик силовых панелей при высокоинтенсивном диффузном нагружении в лабораторных условиях не проводился, а соответствующие параметры ЗИ вычислялись в предположении инвариантности звукоизоляции при высоких нагрузках или получались по результатам натурных измерений фактического акустического давления.
В настоящей статье для определения характеристик звукоизоляции предлагается использовать метод наименьших квадратов, с помощью которого получены линии положительного и отрицательного тренда. Данные линии тренда позволяют судить о характере изменения ЗИ в зависимости от угла их наклона к оси изменения акустической нагрузки, оценивать уровни звукоизоляционных характеристик при повышенных акустических нагрузках, которые невозможно реализовать в лабораторных условиях.
Объекты испытаний и оснастка. Экспериментальные исследования проводились в реверберационной камере объемом 1000 м3. Характеристики каждого образца приведены в табл. 1.
Таблица 1.
Параметры образца | Образец 1 | Образец 2 | Образец 3 |
---|---|---|---|
Поверхностная масса с ЗПМ, кг/м2 | 16.53 | 12.33 | 18.78 |
Поверхностная масса без ЗПМ, кг/м2 | 15.22 | 11.88 | 18.35 |
Масса оснастки с переходными элементами и испытуемым образцом, кг | 1000.9 | 911.2 | 998.5 |
Толщина ЗПМ, мм | 60 | 1.5 | 20 |
Плотность ЗПМ, кг/м3 | 35 | 360 | 35 |
Материал ЗПМ | пенополиуретан | ТТПСМ | пенополиуретан |
Из табл. 1 видно, что исследуемые образцы имеют достаточно высокую поверхностную массу за счет вклада силового набора, тогда как поверхностные массы несущих панелей образцов 1, 2, 3 составляют 5 кг/м2, 3.34 кг/м2, 3.34 кг/м2 соответственно, что составляет 30%, 27%, 18% от поверхностных масс образцов со звукопоглощающими материалами (ЗПМ). В силу особенностей эксплуатации поверхностная масса ЗПМ второго образца наносилась на внешнюю поверхность, тогда как у первого и третьего образцов на внутренние поверхности.
Испытаниям подвергались образцы от различных конструктивных зон сборочно-защитного блока (СЗБ), представляющие собой сегменты двух цилиндрических и одной конической оболочки (рис. 1).
Выбор толщин стенок оснастки проводился из следующих соображений. Предположим, что поверхность испытуемой панели площадью S1 излучает во внутреннее пространство акустическую энергию интенсивности I1, тогда соответствующую акустическую мощность N1 можно представить через поток акустической энергии проходящей через S1
Тоже самое справедливо для экспериментальной оснастки, если соответствующим параметрам присвоить индекс 2. Тогда суммарная акустическая мощность внутри оснастки запишется
(1)
${{N}_{\Sigma }} = \iint\limits_{{{S}_{1}}} {{{I}_{1}}ds} + \iint\limits_{{{S}_{2}}} {{{I}_{2}}ds}.$Предположим, что вследствие диффузности акустического поля, воздействие на все элементы конструкции оснастки и экспериментальной панели с внешней стороны одинаковые по амплитуде, поделим получившееся выражение (1) на пороговое значение мощности N0 = 10–12 Вт, раскрыв знаки интегралов, заменив величину I2 на I1/k, где k – параметр, подлежащий определению
(2)
$k = \frac{{{{S}_{2}}}}{{{{S}_{1}}}} \times {{10}^{{\left( {{{R}_{2}} - {{R}_{1}}} \right){\text{/}}10}}},$Переходя к логарифмическим величинам, получим для уровня акустической мощности, LΣ
(3)
${{L}_{\Sigma }} = 10\lg \frac{{{{I}_{1}}{{S}_{1}} + {{I}_{1}}{{S}_{2}}{\text{/}}k}}{{{{N}_{0}}}}.$Погрешность влияния оснастки, Δ определим по формуле
где ${{L}_{1}} = 10\lg \frac{{{{I}_{1}}{{S}_{1}}}}{{{{N}_{0}}}}$ – уровень акустической мощности, соответствующей излучению только экспериментальной панели.Чтобы исключить погрешность влияния оснастки, потребуем
(4)
$\Delta = 10\lg \left( {1 + \frac{{{{S}_{2}}}}{{k{{S}_{1}}}}} \right) \to 0,\quad {\text{при}}\quad k \to \infty .$Вычисления (4) с учетом (3) показывают, что погрешность эксперимента будет менее 1 дБ при изготовлении оснастки например из материала АМГ6 толщиной 30 мм.
Элементы конструкции монтировались в экспериментальную оснастку, в составе которой предусматривались необходимые переходные элементы для установки трех образцов, а также материальная часть в виде матов пенополиуретана толщиной 60 мм для звукоизоляции стенок оснастки. Глубина оснастки выбиралась исходя из расстояния от обечайки до солнечных батарей изделия, поперечные размеры ограничивались воротным проемом АРС.
Оценка характеристических параметров. Эффект положительного и отрицательного тренда. Для получения обобщенных характеристик звукоизоляции образцов для всех временных промежутков, выбранных в процессе испытаний (для образцов 1 и 2 их было 151, для образца 3 – 201), было проведено осреднение полученных уровней акустического давления снаружи Lн( f ) и внутри Lв( f ) экспериментальной оснастки c образцами по формулам согласно [9]
(5)
${{L}_{{\text{н}}}}\left( f \right) = 10\lg \left( {\sum\limits_{i = 1}^n {{{{10}}^{{{{L}_{{{{A}_{i}}}}}{\text{/}}10}}}} {\text{/}}n} \right),$(6)
${{L}_{{\text{в}}}}\left( f \right) = 10\lg \left( {\sum\limits_{i = 1}^n {{{{10}}^{{{{L}_{{{{A}_{i}}}}}{\text{/}}10}}}} {\text{/}}n} \right),$После этого для указанных временных промежутков ( j = 1…m, m = 151, 201) были получены характеристики ЗИ Rj( f ) и kj( f )
(7)
$\begin{gathered} {{R}_{j}}(f) = {{L}_{{{\text{н}}j}}}\left( f \right) - {{L}_{{{\text{в}}j}}}\left( f \right), \\ {{k}_{j}}\left( f \right) = {{10}^{{{{R}_{j}}(f){\text{/}}20}}}, \\ \end{gathered} $(8)
$\begin{gathered} {{k}_{{{\text{ср}}}}}\left( f \right) = \sqrt {\sum\limits_{j = 1}^m {\frac{{k_{j}^{2}\left( f \right)}}{m}} } , \\ {{R}_{{{\text{ср}}}}}(f) = 20\lg \left( {{{k}_{{ср}}}\left( f \right)} \right). \\ \end{gathered} $Осредненные характеристики ЗИ, вычисленные по формулам (5)–(8), показаны на рис. 2.
Экспериментальные результаты исследования звукоизоляционных характеристик свидетельствуют о следующей тенденции их изменения вдоль частотной оси:
1) падение звукоизоляции в полосе частот с центральной частотой 100 Гц обусловлено поперечным пространственным резонансом, который можно определить по формуле, приведенной в работе [11] $f = \frac{c}{{2l}}$, где с = 343 м/с – скорость звука в воздухе при нормальных условиях; l – максимальный поперечный размер воздушного пространства под оснасткой.
Перечень некоторых пространственных резонансов представлен в табл. 2;
Таблица 2.
Номер моды | Пространственные резонансные частоты, Гц | ||
---|---|---|---|
Образец 1 | Образец 2 | Образец 3 | |
1 | 98.1 | 99.8 | 93.2 |
2 | 100.8 | 103.1 | 100.8 |
3 | 137.4 | 139.0 | 137.1 |
4 | 174.8 | 173.5 | 186.3 |
5 | 192.4 | 195.4 | 191.2 |
6 | 206.5 | 206.0 | 211.6 |
7 | 213.9 | 214.3 | 212.6 |
2) падение звукоизоляции в полосах с центральными частотами 160 Гц, 200 Гц обусловлено нормальными к поверхности экспериментальной панели колебаниями. В табл. 3 представлены структурные резонансы оснастки с экспериментальными образцами;
Таблица 3.
Образец 1 | частота, Гц | 134.8 |
эффективная масса, кг | 180 | |
Образец 2 | частота, Гц | 218.6 |
эффективная масса, кг | 271.8 | |
Образец 3 | частота, Гц | 151 |
эффективная масса, кг | 123.3 |
3) после структурного резонанса, попадающего в полосы частот с центральными частотами 160 Гц и 200 Гц, возрастание звукоизоляции вдоль частотной оси обусловлено массовыми свойствами обечайки. На рис. 2 для образцов представлены данные расчетов по ЗМ. Влияние массы перекрестного силового набора на спектры звукоизоляции отсутствует, особенно на высоких частотах, где можно применить ЗМ (рис. 2).
Ежесекундная оценка уровней акустического давления при эксперименте (от меньшего уровня к большему уровню) позволило провести исследование изменения звукоизоляционных характеристик для рассматриваемых образцов. Для этого вычислялись характеристики изменения ЗИ, ∆R панели при изменении уровней внешнего акустического давления, ∆L относительно некоторого минимального уровня. Примеры полученных результатов для первого образца продемонстрированы на рис. 3.
Исходя из полученных результатов, для оценок изменения звукоизоляции были построены линии тренда по формулам [9]
Данные рис. 3 свидетельствуют о том, что параметр ΔR принимает как положительные, так и отрицательные значения при увеличении нагружения ΔL. В то же время наклон линии тренда показывает характер параметров ЗИ в исследуемой полосе частот и позволяет предположить дальнейшее изменение звукоизоляционной способности при увеличении внешних акустических нагрузок, которые не были реализованы в лабораторных условиях, но ожидаются при натурных испытаниях изделия.
На рис. 4 представлено изменение угла линий тренда от частоты при испытаниях образцов 1, 2, 3.
Наблюдается увеличение угла до значений соответствующих полосе с центральной частотой 63 Гц и далее происходит его спад вплоть до отрицательных значений, что говорит о снижении влияния эффекта положительного тренда на высоких частотах.
На рис. 5 изображено изменение ЗИ при увеличении уровня акустического давления на 20 дБ.
Данные рис. 4, 5 свидетельствуют об общих закономерностях изменения звукоизоляции для различных конструкций. Имеется ряд частотных диапазонов, в рамках которых наблюдается увеличение угла линии тренда в изменении ЗИ. Для образцов 1 и 3 таковыми являются диапазоны с центральными частотами до 63 Гц, для образца 2 – до 50 Гц и до 100 Гц. Далее наблюдается уменьшение угла линии тренда изменения ЗИ до значения, соответствующего диапазону с центральной частотой 200 Гц, после чего угол изменяется незначительно и даже может принимать отрицательные значения.
Полученный эффект можно объяснить с помощью данных табл. 1 и 3. Как можно видеть из табл. 1 массы сборок, включающих в себя оснастку, переходные элементы и испытуемую конструкцию, составляют 1000.9 кг, 911.2 кг и 998.5 кг соответственно для образцов 1, 2, 3, тогда как основные собственные частоты, представленные в табл. 3, составляют 134.8 Гц, 218.6 Гц, 151 Гц. Из этого можно сделать вывод, что сборка с образцом № 2 обладает большей жесткостью, нежели сборки с образцами 1 и 3. У данного образца (рис. 5) эффект положительного тренда проявлен в меньшей мере нежели у двух других именно в частотных полосах, где в большей степени звукоизоляция обусловлена жесткостью сборки [9], это подтверждается, в частности, данными рис. 3, на котором можно видеть снижение ЗИ в плоть до резонансных частот (табл. 3), попадающих в диапазон с центральными частотами 160 Гц и 200 Гц. На более высоких частотах, где звукоизолирующий эффект обусловлен массой конструкции, эффект положительного тренда проявлен в меньшей мере, можно сказать, что ЗИ на данных частотах в большей степени близка к линейным законам звукоизолирующих преград.
Пример оценки звукоизоляции при высокоинтенсивном акустическом нагружении. Проведем оценку звукоизоляции конструкции образца 1, при его акустическом нагружении с суммарным среднеквадратическим уровнем 162 дБ, отметим, что при лабораторных испытаниях по ряду 1/3-октавных диапазонов амплитудные значения не были достигнуты, но при анализе результатов были получены линии положительного и отрицательного тренда, позволяющие провести необходимую оценку. Уравнения линий тренда для каждой полосы частот, разность между целевым уровнем, Ltarg и минимальным получившимся уровнем акустического давления в процессе испытаний, Lmin представлены в табл. 4. Чтобы провести оценку неточностей испытаний вместо системы координат (ΔL, ΔR) с точками (рис. 3), была введена новая система координат (ΔL', ΔR'), ось абсцисс ΔL', которой совпадал с линией тренда, а ось ΔR' была перпендикулярна ΔL'.
Таблица 4.
Центральная частота 1/3-октавного спектра, Гц | Уравнение линии тренда, ΔR = aΔL | Ltarg – Lmin, дБ | σ, дБ |
---|---|---|---|
25 | ΔR = 0.12ΔL | 56.1 | 3.7 |
31.5 | ΔR = 0.33ΔL | 48.0 | 4.8 |
40 | ΔR = 0.22ΔL | 39.3 | 4.3 |
50 | ΔR = 0.65ΔL | 34.4 | 5.0 |
63 | ΔR = 0.79ΔL | 33.6 | 3.4 |
80 | ΔR = 0.73ΔL | 31.9 | 3.6 |
100 | ΔR = 0.35ΔL | 28.8 | 5.1 |
125 | ΔR = 0.29ΔL | 26.8 | 2.7 |
160 | ΔR = 0.11ΔL | 25.9 | 3.4 |
200 | ΔR = –0.02ΔL | 25.9 | 2.6 |
250 | ΔR = 0.11ΔL | 25.8 | 2.3 |
315 | ΔR = 0.05ΔL | 26.9 | 2.5 |
400 | ΔR = 0.04ΔL | 30.1 | 2.3 |
500 | ΔR = –0.05ΔL | 31.3 | 2.3 |
630 | ΔR = –0.02ΔL | 37.8 | 2.9 |
800 | ΔR = 0.07ΔL | 39.9 | 2.9 |
1000 | ΔR = 0.06ΔL | 41.7 | 2.8 |
1250 | ΔR = –0.05ΔL | 39.1 | 3.6 |
1600 | ΔR = –0.07ΔL | 44.3 | 4.5 |
2000 | ΔR = –0.12ΔL | 51.8 | 4.5 |
2500 | ΔR = –0.12ΔL | 61.8 | 4.4 |
3150 | ΔR = –0.14ΔL | 70.4 | 5.0 |
4000 | ΔR = –0.10ΔL | 77.1 | 5.1 |
Новые координаты точек вычислялись по формулам
После этого вычислялись среднеквадратические отклонения изменения звукоизоляции относительно оси линии тренда, σ по формуле
На рис. 6 представлены значения звукоизоляции при минимальном и максимальном уровнях внешнего акустического нагружения и найденные значения звукоизоляции при целевом уровне внешнего акустического нагружения.
Предложенная оценка (рис. 6), позволила предсказать рост ЗИ при целевом уровне акустического нагружения в диапазонах с центральными частотами 25–125 Гц до значений близких к ЗИ, соответствующей максимальному внешнему нагружению, для более высоких частот полученные значения ЗИ оказались близки к значениям ЗИ, соответствующим минимальным уровням нагружения.
Заключение. В рамках настоящей статьи были получены следующие результаты. 1. Выявлен эффект положительного и отрицательного тренда изменения звукоизоляционных характеристик при повышении акустических нагрузок, который позволяет судить о характере изменения звукоизоляции. 2. Получены общие закономерности изменения звукоизоляции при повышении нагрузки, показано, что в частотной области, где звукоизоляция обусловлена жесткостными характеристиками наблюдается ярко выраженный увеличивающийся положительный тренд исследуемой величины, после чего имеет место снижение влияния положительного тренда, на высоких частотах, где звукоизоляция обусловлена массовыми свойствами конструкции звукоизоляция близка к линейным закономерностям. 3. Полученный эффект положительного и отрицательного тренда звукоизоляционных характеристик позволяет оценить звукоизоляцию при повышенных акустических нагрузках, которые не возможно реализовать в лабораторных условиях, но которые могут наблюдаться при натурной эксплуатации.
В силу того факта, что эффект положительного и отрицательного тренда был выявлен случайно при обработке результатов экспериментов для элементов конструкции СЗБ, в рамках настоящей статьи не представлены предельные акустические нагрузки, при которых угол линии тренда к оси изменения нагрузки был бы равен нулю и звукоизоляция конструкции не изменялась. Поэтому дальнейшие исследования будут проведены для решения этого вопроса. Кроме этого, полученный эффект нуждается в математической интерпретации.
Список литературы
Абдрашитов Р.Г., Иванушкин Е.А., Логинов Г.М., Попов О.Ю., Шарунов А.В. Расчетно-экспериментальные исследования звукоизоляции подвесных контейнеров // Тезисы докладов пятой открытой Всероссийской (XVII научно-технической) конференции по аэроакустике (25–29 сентября 2017 г.). М.: Изд. отдел ЦАГИ, 2017. С. 213.
Карачун В.В., Мельник В.Н. Основные причины шума ракет-носителей // Авиационно-космическая техника и технология. 2007. № 9 (45). С. 7.
Popov P.A., Igolkin A.A., Shakhmatov E.V. Calculation of the Acoustic Field on the Surface of a Complex-Shaped Object within a Radiating Cylindrical Shell // Acoustical Physics. 2018. V. 64. Iss. 5. P. 637.
Яблоник Л.Р. Упрощенный метод расчета многослойной звукоизоляции, включающей слои волокнистого пористого материала // Акуст. журн. 2018. Т. 64. № 5. С. 639.
Liu Y., Liu X., Xu J., Hu X., Xia Z. Inverse identification of the acoustic porous parameters of double-layered poroelastic structures by acoustic rigidity approximation // J. Acoust. Soc. Am. 2017. V. 142. P. 72.
Chronopoulos D., Ichchou M., Troclet B., Bareille O. Predicting the broadband vibroacoustic response of systems subject to aeroacoustic loads by a Krylov subspace reduction // Applied Acoustics. 2013. № 74. P. 1394.
Белоус А.А., Корольков А.И., Шанин А.В. Экспериментальная оценка частотной зависимости коэффициента отражения звукопоглощающего материала при наклонном падении // Акуст. журн. 2018. Т. 64. № 2. С. 155.
Билык Н.А. Коэффициенты безопасности и запаса прочности механических конструкций // Проблемы машиностроения и надежности машин. 2019. № 5. С. 63.
Боголепов И.И. Промышленная звукоизоляция. Ленинград: Судостроение, 1986. 367 с.
Иванов Н. И. Инженерная акустика. Теория и практика борьбы с шумом. М.: Университетская книга, Логос, 2008. 424 с.
Шендеров Е.Л. Волновые задачи гидроакустики. Ленинград: Судостроение, 1972. 352 с.
Дополнительные материалы отсутствуют.
Инструменты
Проблемы машиностроения и надежности машин