Математическая модель исследования нормальных акустических мод отсеков ракеты-носителя Текст научной статьи по специальности «Физика»

УДК 534.83+ 629.78

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ИССЛЕДОВАНИЯ НОРМАЛЬНЫХ АКУСТИЧЕСКИХ МОД ОТСЕКОВ РАКЕТЫ-НОСИТЕЛЯ

© 2012 П. А. Попов1, А. С. Белов1, А. Н. Крючков2

1ФГУП ГНПРКЦ «ЦСКБ-ПРОГРЕСС», г. Самара 2Самарский государственный аэрокосмический университет имени академика С. П. Королёва (национальный исследовательский университет)

Изложены метод, учитывающий влияние баков горючего и окислителя ракеты-носителя (РН) на распространение звука в межбаковом пространстве, пример нахождения резонансных акустических мод и частот отсека РН, анализ звукоизоляционных характеристик отсека.

Ракета-носитель, отсеки, акустическое давление, звукоизоляция, акустические и структурные моды, волновое уравнение, стоячая волна.

При проектировании изделий ракетно-космической техники (РКТ) одним из направлений работ инженеров является прогнозирование распространения акустических полей внутри отсеков ракет-носителей, так как акустическое поле, создаваемое струями двигательной установки, структурным шумом и турбулентным аэродинамическим потоком, является источником вибрации конструкции РН и бортовой аппаратуры при старте и полёте. Характеристики акустического поля внутри замкнутых цилиндрических отсеков РН неразрывно связаны с определением звукоизоляции их панелей. При этом, как показано в работах [1,2], явление снижения звукоизоляционной способности панелей обусловлено появлением структурных и акустических резонансов и их взаимодействием.

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

ние математической модели распространения акустических полей внутри отсека РН и анализ данных звукоизоляционных характеристик отсека на основе этой модели.

Пусть имеется цилиндрический отсек радиусом R и высотой Н, ограниченный на основаниях сферическими поверхностями с радиусами кривизны 5, при этом минимальное расстояние между сферическими поверхностями равно

А = H - 2(8 - 48^-Я^). Радиальные,

продольные и окружные моды акустического давления P в цилиндрическом отсеке описываются волновым уравнением в цилиндрической системе координат (г , р, г):

_ _ _ _ 1 Э2 Р

Э2 Р 1 ЭР 1 Э2 Р Э2 Р

~ЭТГ+Г ~ЭТ+Г р+~дгг С2 ~Э^

где с - скорость звука в среде, м/с.

0,

Из данного соотношения, используя

подстановку Р(г, р, г, /) = Р(г, р, г)е" лучим уравнение Гельмгольца:

Э 2Р 1 ЭР 1 Э 2Р Э 2Р ,2

по-

- + —

+ —г + Я2 Р = 0, (1)

Эх

Эт2 т Эт т2 Эр2

где Я = — - волновое число, рад /м; с

— -круговая частота, радТц.

Граничные условия для (1) выбираются следующим образом:

- на краях цилиндра и на сферических поверхностях производные по направлениям г и г от амплитуды акустических волн равны нулю;

- окружные моды внутри отсека повторяются с периодичностью 2^;

- на центральной линии отсека решение не является бесконечным.

Граничные условия для криволинейных границ днищ баков получены из геометрии, изображённой на рис. 1.

ё2 Я(т) 1 dR.iT) 1 ё 2Ф(р)

ёт2 т ёт + т 2 ёр 2 +

Жт) Ф(р) (4)

ё2 2 (2)

++ % = о.

2 (2)

Выражение (4) выполняется, когда ё2R(т) + 1 dR(т) 1 ё2Ф(р)

о + 2 7 2

ёт_______т ёт + Т р = _п 2 (5)

R(т) Ф(р) ~~т ' ()

Рис.1. Вспомогательная геометрия половины отсека РН:

1 - каркас отсека РН; 2 - дни-ще бака РН

Граничные условия имеют вид:

ЭР

— (т = R, р, 2) = 0,

Эт

ЭР

— (т, (р, 2 = ±И) = 0,

Эх

Р(т,р, 2) = Р(т,р + 2р, 2),

Р(т = 0, р, 2) < ¥,

И = ^ + 8 -482—Ё2.

где

2

Задача решается методом разделения переменных:

Р(т, р, 2, і) = R(т)Ф(р) 2 (2)Т (і).

(3)

Подставим (3) в (1) и поделим получившееся на (3):

ё2 2 (2)

ё2 2 =-п 2

2 (2) ' ’

(Я2 = т2 + л2).

Первое из уравнений (5) умножим на т2 и после преобразования получим

2 ё2R(r) dR(т)

т----+ т Л 2

—+т 2т2 +-р- = 0. (6) R(т) м Ф(р) V ^

Выражение (6) выполняется, когда

2 ё2Чт) + тёй(т)

ёт

R(т)

+ т 2т2 = а2,

(2) Ф(р)

2Ф(р)

ёр 2 2

г = -а .

(7)

Второе из уравнений (7) описывает окружные моды внутри отсека. Рассмотрим следующую задачу нахождения собственных чисел а{ (задача Штурма-Лиувилля):

Ф(2р) + а 2Ф(р) = 0,

ёр

Ф(р) = Ф(р + 2р).

Из (8) окружные моды примут вид

Фг (р) = Аг соъ(агр) + Бг ), Таблица 1. Перечеть собственных чисел

(9) и соответствующих им резонансных частот

(аг = I € 7).

Первое из выражений (7) является уравнением Бесселя, из (9) можно сделать вывод, что его порядок - целое число. Решения представляют собой сумму двух линейно независимых функции Бесселя и Неймана:

Ri (г) = С1 (т • г)+Б?г (т • г),

где ^ (т • г) - функция Бесселя г - го порядка;

Уг (т • г) - функция Неймана г - го порядка.

Известно, что функции Неймана сингулярны в нуле, поэтому в соответствии с последним из условий (2) можно сделать вывод, что радиальные моды будут представлены лишь функциями Бесселя:

R (г) = С1 (ткг), (10)

где к - номер корня функции Бесселя порядка г.

Исходя из первого из условий (2), найдём собственные числа тк (их перечень и резонансные частоты, соответствующие радиальным модам, представлены в табл. 1: 1 - порядок г функции (ткК); 2 - номер корня, к; 3 - корень,

mkR; 4 - радиальная резонансная частота

Были рассмотрены радиальные и круговые моды акустического поля. Для дальнейшего анализа рассмотрим продольные моды и найдём акустические резонансы. Для этого воспользуемся вторым уравнением из (5) с граничными условиями вида

1 2 3 4

0 1 3,84 152,4

2 7,01 278,2

3 10,18 404,0

4 13,3 527,9

5 16,48 654,1

6 19,61 778,3

7 22,77 903,7

1 1 1,8 71,4

2 5,33 211,5

3 8,53 338,5

4 11,7 464,4

5 14,86 589,8

6 18,02 715,2

7 21,16 839,8

2 1 3,05 121,1

2 6,7 265,9

3 9,96 395,3

4 13,17 522,7

5 16,34 648,5

6 19,51 774,1

3 1 4,2 166,7

2 8,01 317,9

3 11,34 450,1

4 14,58 578,7

5 17,78 705,7

6 20.97 166,7

4 1 5,31 210,7

2 9,28 368,3

3 12,68 503,3

4 15,96 633,4

5 19,19 761,6

5 1 6,41 254,4

6 1 7,5 297,7

7 1 8,57 340,1

8 1 9,64 382,6

9 1 10,71 425,1

Условие (11) даст:

Л

- рл з Э1п(л з(у+8-^8 2 - т 2)) + + Ел соъ(лзА + 8 -л/82 -т2)) = 0

з2 .А

з Э1п(лз(— + 8-^82 -т2)) +

+ Ел . соъ(лз(А + 8 -л/82 -т2)) = 0.

(12)

2

П соъ(п (- + 5-л/52 - г2 )>12,

П ЪШ(П (А + 5^5 2 - г2 )) 21,

(13)

П j С08(^ (- +5^52 - г2 )) 22 = 0.

Уравнение (13) преобразуется к виду А

ъ1п(2лз (— + 8 - л/82 - т2)) = 0.

(14)

Имеется набор собственных чисел, полученных из (14):

р • з

л

А + 2(8 -V82 -т2)'

Продольные моды имеют вид

р • 7

2(2) = рзС0э(--------:—и 2 2 2) +

А + 2(8-V8 -т2)

+ Е, ъ1п(

р • з

А + 2(8 -л/82 -т2)

2).

Найдём величину Я*:

= т: + (

р • з

А + 2(8 -V52 - т2)

)2, (17)

тІ - находится из табл. 1.

Подставим в (4) выражения (10), (11), (16):

Как известно, уравнение (12) имеет решение при условии:

ёе^-'Пj яп(пДу + 5-л/52 -г2))п

Р = Цф)(4 СОЪф) +

+ът^фУ)^. соъ( Р з

+Ез ъ1п(-

■/™А+2(8-V82 -т2) Р • з

2) +

А+2(8-л/ 82 - т2)

2)).

(18)

На рис. 2-4 показаны радиальные акустические моды, возникающие в меж-баковом пространстве на частоте 71,4 Гц, на различных высотах от срединной плоскости между баками (Н=0 м) до высоты

Н = А (днище бака).

(15) Рис. 2. Акустическая мода, распространяющаяся

по срединной плоскости отсека,

Н=0 м (функция Бесселя нулевого порядка)

(16)

Из (16) следует, что вдоль центральной линии (г=0), между днищами баков межбакового отсека имеется семейство стоячих (продольных) волн:

Zj(г) = ^ соз(—а^" г) + Е} эш(^ г).

Рис. 3. Акустическая мода,

распространяющаяся на высоте н = А

4

Рис. 4. Акустическая мода,

распространяющаяся на высоте н = А

2

Как видно из рис. 2, вдоль центральной плоскости цилиндрического отсека со сферическими поверхностями на основаниях распространение акустических волн происходит аналогично цилиндру с плоскими основаниями. Однако при изменении высоты Д влиянием баков горючего и окислителя РН пренебречь нельзя (см. рис. 3,4) и необходимо использовать модель, предложенную в данной статье.

Для дальнейшего анализа найдём структурные резонансы обечайки отсека

РН. Уравнение динамики при малых смещениях точек тел имеет следующий вид:

[М ]х и +[б]х и + [К ]х и = Р (і), (19)

где [М] - матрица масс;

и - вектор ускорений в узловых точках;

[В] - матрица демпфирования;

и - вектор скоростей в узловых точках;

[К] - матрица жёсткости структуры;

и - вектор перемещений узлов;

Р (і) - вектор узловых сил как функция времени.

Принимая равными нулю воздействующие силы, уравнение (19) преобразуем к виду

[М ]х и +Б]х и + [к ]х и = 0. (20)

Решение уравнения (20) показало, что собственные оболочечные колебания обечайки отсека РН на первом тоне составляют 270 Гц. Форма колебания представлена на рис. 5.

Десять оболочечных тонов представлены в табл. 2.

б

Рис. 5. Формы собственных колебаний каркаса отсека на частотах: а - 270 Гц; б - 273 Гц Таблица 2. Оболочечные тона

Номер тона 1 2 3 4 5 6 7 8 9 10

Собственная частота, Гц 270 273 283 287 289,5 294,5 322 349 357 370

а

Как отмечалось ранее, звукоизоляционные характеристики зависят от распределения акустических и структурных резонансов внутри отсека и на его поверхности. Выше такие распределения были получены. Эти данные могут объяснить падения звукоизоляции при пусках РН типа "Союз" в области частот 160-600 Гц. Действительно, как показано в табл.1, акустические резонансы начинаются с 71,4 Гц, а основная их доля попадает в области частот 1/3-октавного спектра с центральными частотами 160, 200, 250, 315, 625 Гц. Оболочечные резонансы (табл. 2, рис. 6) начинаются с 270 Гц. С учётом их высокой плотности следует вывод, что могут возникать на частотной оси поддиапазоны, в которых структурные резонансы совпадают с акустическими резонансами, а основные минимумы звукоизоляции попадают в 1/3 октавные полосы частот с центральными частотами 160, 250, 315, 625 Гц.

В табл. 3 показаны соотношения между акустическими и структурными резонансами и центральные частоты в 1/3октавном спектре, на которых наблюдается падение звукоизоляционных характеристик панели отсека (рис. 6): 1 - структурные резонансы подкреплённой обо-

Рис. 6. Перепад звукового давления через обечайку отсека во время полёта РН при различных пусках

лочки отсека, /рез, Гц; 2 - радиальная акустическая частота отсека, /рез, Гц; 3 -продольная акустическая частота (между полусферами днищ баков), /рез, Гц; 4 -места "провалов" в звукоизоляционном спектре, /-У ,Гц.

&

Таблица 3. Соотношения между акустическими и структурными резонансами и частотами, на которых наблюдается падение коэффициента звукоизоляции панели

1 2 3 4

71,44

121,1

152,4 160

166,7

210,7 200

270 254 250

273

287

289 297,7

322 340,1 315

382,6

425,1

637,5 625

Таким образом, разработана математическая модель распространения акустических мод внутри замкнутого отсека, ограниченного двумя сферическими поверхностями. Исследованы акустические и упругие резонансные частоты, значения которых представлены в табл. 1, 2. Проведённый анализ объясняет снижение звукоизоляционных способностей отсека в 1/3-октавных полосах частот с центральными частотами 160, 250, 315, 625 Гц.

Полученная математическая модель полезна при прогнозировании акустического нагружения аппаратуры и конструкции РН с большими диаметрами отсеков (например, РН "Союз-2-3-в") и может быть модифицирована для случаев распространения акустических полей между

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

Библиографический список

1. Исследования по повышению звукоизоляции самолетных конструкций (по материалам открытой зарубежной печати) [Текст]. - ЦАГИ, 1990. - 101 с.

2. Экспериментальное определение звукоизоляции межбакового отсека блока III ступени РН "Союз" [Текст] / НТО, инв.№5535. - ЦАГИ, 1999. - 38 с.

3. Владимиров, В . С. Уравнения математической физики [Текст] / В. С. Владимиров - М.: Наука, 1976. - 392 с.

MATHEMATICAL MODEL OF INVESTIGATING NORMAL ACOUSTIC MODES

OF LAUNCH VEHICLE MODULES

© 2012 P. A. Popov1, A. S. Belov1, A. N. Kryuchkov2

1State Research and Production Space Rocket Center «TsSKB-Progress»

Samara State Aerospace University named after academician S. P. Korolyov (National Research University)

The article presents a method which considers the influence of fuel and oxidizer tanks of a launch vehicle on sound propagation in the inter-tank space. An example of finding resonant acoustic modes and frequencies of a launch vehicle module and the analysis of sound-insulation characteristics of the module are also given.

Launch vehicle, rocket modules, acoustic pressure, sound insulation, acoustic and structural modes, wave equation, standing wave.

Информация об авторах

Попов Павел Александрович, начальник группы, ФГУП ГНПРКЦ «ЦСКБ-ПРОГРЕСС». E-mail: banduir@rambler.ru. Область научных интересов: звукоизоляция и акустические нагрузки.

Белов Анатолий Сергеевич, инженер-конструктор, ФГУП ГНПРКЦ «ЦСКБ-ПРОГРЕСС». E-mail: AnatoliyBelov@Gmail.com. Область научных интересов: исследование нагружения оборудования КА и РН от вибрационных и акустических воздействий.

Крючков Александр Николаевич, доктор технических наук, профессор кафедры автоматических систем энергетических установок, Самарский государственный аэрокосмический университет имени академика С. П. Королёва (национальный исследовательский университет). E-mail: kan@ssau.ru. Область научных интересов: виброакустика и динамика машин.

Popov Pavel Aleksandrovich, chief of group, Space Rocket Center «TsSKB-Progress». E-mail: banduir@rambler.ru. Аrea of research: sound insulation and acoustic loads.

Belov Anatoly Sergeevich, design engineer, Space Rocket Center «TsSKB-Progress». E-mail: AnatoliyBelov@Gmail.com. Area of research: space vehicle and launch vehicle equipment loading caused by vibration and acoustic influence.

Kryuchkov Alexander Nikolaevich, doctor of technical science, professor, the department of automatic systems of power installations, Samara State Space University named after academician S. P. Korolyov (National Research University). E-mail: kan@ssau.ru. Аrea of research: vibration of acoustics and dynamics of machines.