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

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

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Datum:2020
Hauptverfasser: Лимарченко, О.С., Нефедов, А.А., Семенович, Е.А.
Format: Artikel
Sprache:Russian
Veröffentlicht: Інститут механіки ім. С.П. Тимошенка НАН України 2020
Schriftenreihe:Прикладная механика
Online Zugang:http://dspace.nbuv.gov.ua/handle/123456789/188290
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Назва журналу:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Zitieren:Проявление вторичных резонансов при маятниковых колебаниях резервуаров с жидкостью / О.С. Лимарченко, А.А. Нефедов, Е.А. Семенович // Прикладная механика. — 2020. — Т. 56, № 6. — С. 105-115. — Бібліогр.: 21 назв. — рос.

Institution

Digital Library of Periodicals of National Academy of Sciences of Ukraine
id irk-123456789-188290
record_format dspace
spelling irk-123456789-1882902023-02-21T01:27:01Z Проявление вторичных резонансов при маятниковых колебаниях резервуаров с жидкостью Лимарченко, О.С. Нефедов, А.А. Семенович, Е.А. В настоящей работе показан характер изменения частот свободных колебаний резервуаров с жидкостью на основе модели совместного движения составляющих компонент системы для разных вариантов движения несущего тела и показано, что в системе при определенных вариантах маятникового закрепления несущего тела (которые исследуются значительно сложнее, чем случаи поступательного движения несущего тела) наблюдается проявление вторичных резонансов, под которыми понимаются резонансы, устанавливающиеся на основе нелинейных механизмов, которые не проявляются вообще в рамках линейного моделирования системы. Розглянуто задачу про коливання резервуара на маятниковому підвісі з рухомою точкою підвісу, частково заповненого ідеальною нестисливою рідиною. Задачу розглянуто в нелінійній сумісній постановці. Показано, що розгляд сумісності руху компонент системи призводить до суттєвої зміни величин частот и порядку розташування форм коливань системи при їх розміщенні в порядку зростання частот. Встановлено, що для коротких і середніх довжин маятникових підвісів перша антисиметрична форма, яка зазвичай є головною в розвитку резонансних процесів, тепер збуджується не на своїй резонансній частоті, а на частоті маятникових коливань і частотах вторинних (нелінійних) резонансів на основі нелінійного перерозподілу енергії. Показано відповідність розвитку динамічних процесів в системі основним закономірностям нелінійного хвилеутворення, встановленим в теоретичних і експериментальних роботах. A problem of oscillation of the reservoir on pendulum suspension with the movable suspension point and partially filled with the ideal incompressible liquid is considered. This problem is studied in the nonlinear joint statement. It is shown that consideration of the compatibility of the system components motion results in the considerable changingof values of frequencies and the order of arrangement the system normal modes of oscillationsfor their allocation according to ascending order of frequencies. It was ascertained thatfor the short and medium lengths of pendulum suspensions, the first anti-symmetric mode, which is usually the main in the development of resonant processes, is excited not on its resonant frequency, but the frequency of pendulum oscillations and frequencies of secondary (nonlinear) resonances on the bases of nonlinear redistribution of energy. A correspondence of development of the dynamical processes in the system with fundamental regularities of nonlinear wave generation, ascertained in еру theoretical and experimental publications, is shown. 2020 Article Проявление вторичных резонансов при маятниковых колебаниях резервуаров с жидкостью / О.С. Лимарченко, А.А. Нефедов, Е.А. Семенович // Прикладная механика. — 2020. — Т. 56, № 6. — С. 105-115. — Бібліогр.: 21 назв. — рос. 0032-8243 http://dspace.nbuv.gov.ua/handle/123456789/188290 ru Прикладная механика Інститут механіки ім. С.П. Тимошенка НАН України
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection DSpace DC
language Russian
description В настоящей работе показан характер изменения частот свободных колебаний резервуаров с жидкостью на основе модели совместного движения составляющих компонент системы для разных вариантов движения несущего тела и показано, что в системе при определенных вариантах маятникового закрепления несущего тела (которые исследуются значительно сложнее, чем случаи поступательного движения несущего тела) наблюдается проявление вторичных резонансов, под которыми понимаются резонансы, устанавливающиеся на основе нелинейных механизмов, которые не проявляются вообще в рамках линейного моделирования системы.
format Article
author Лимарченко, О.С.
Нефедов, А.А.
Семенович, Е.А.
spellingShingle Лимарченко, О.С.
Нефедов, А.А.
Семенович, Е.А.
Проявление вторичных резонансов при маятниковых колебаниях резервуаров с жидкостью
Прикладная механика
author_facet Лимарченко, О.С.
Нефедов, А.А.
Семенович, Е.А.
author_sort Лимарченко, О.С.
title Проявление вторичных резонансов при маятниковых колебаниях резервуаров с жидкостью
title_short Проявление вторичных резонансов при маятниковых колебаниях резервуаров с жидкостью
title_full Проявление вторичных резонансов при маятниковых колебаниях резервуаров с жидкостью
title_fullStr Проявление вторичных резонансов при маятниковых колебаниях резервуаров с жидкостью
title_full_unstemmed Проявление вторичных резонансов при маятниковых колебаниях резервуаров с жидкостью
title_sort проявление вторичных резонансов при маятниковых колебаниях резервуаров с жидкостью
publisher Інститут механіки ім. С.П. Тимошенка НАН України
publishDate 2020
url http://dspace.nbuv.gov.ua/handle/123456789/188290
citation_txt Проявление вторичных резонансов при маятниковых колебаниях резервуаров с жидкостью / О.С. Лимарченко, А.А. Нефедов, Е.А. Семенович // Прикладная механика. — 2020. — Т. 56, № 6. — С. 105-115. — Бібліогр.: 21 назв. — рос.
series Прикладная механика
work_keys_str_mv AT limarčenkoos proâvlenievtoričnyhrezonansovprimaâtnikovyhkolebaniâhrezervuarovsžidkostʹû
AT nefedovaa proâvlenievtoričnyhrezonansovprimaâtnikovyhkolebaniâhrezervuarovsžidkostʹû
AT semenovičea proâvlenievtoričnyhrezonansovprimaâtnikovyhkolebaniâhrezervuarovsžidkostʹû
first_indexed 2025-07-16T10:17:30Z
last_indexed 2025-07-16T10:17:30Z
_version_ 1837798314626514944
fulltext 2020 П Р И К Л А Д Н А Я М Е Х А Н И К А Том 56, № 6 ISSN0032–8243. Прикл. механика, 2020, 56, № 6 105 О . С . Л и м а р ч е н к о , А . А . Н е ф е д о в , Е . А . С е м е н о в и ч ПРОЯВЛЕНИЕ ВТОРИЧНЫХ РЕЗОНАНСОВ ПРИ МАЯТНИКОВЫХ КОЛЕБАНИЯХ РЕЗЕРВУАРОВ С ЖИДКОСТЬЮ Киевский национальный университет им. Т. Шевченко, просп. Глушкова, 4е, 01033, Киев, Украина; e-mail: olelim2010@yahoo.com Abstract. A problem of oscillation of the reservoir on pendulum suspension with the movable suspension point and partially filled with the ideal incompressible liquid is consid- ered. This problem is studied in the nonlinear joint statement. It is shown that consideration of the compatibility of the system components motion results in the considerable changing of values of frequencies and the order of arrangement the system normal modes of oscilla- tions for their allocation according to ascending order of frequencies. It was ascertained that for the short and medium lengths of pendulum suspensions, the first anti-symmetric mode, which is usually the main in the development of resonant processes, is excited not on its resonant frequency, but the frequency of pendulum oscillations and frequencies of second- ary (nonlinear) resonances on the bases of nonlinear redistribution of energy. A correspond- ence of development of the dynamical processes in the system with fundamental regularities of nonlinear wave generation, ascertained in еру theoretical and experimental publications, is shown. Key words: nonlinear dynamics, combined motion, pendulum suspension, secondary resonances. Введение. При исследовании задач динамики резервуаров с жидкостью со свободной по- верхностью в большинстве случаев предполагается, что движение резервуаров явля- ется заданным, а движение жидкости происходит без учета ее обратного влияния на движение несущего тела [5, 6, 16, 17]. Основная группа практически важных задач характеризуется тем, что масса жидкости значительно превосходит массу «сухой» несущей конструкции и это приводит к необходимости переходить к рассмотрению совместного движения резервуаров с жидкостью. В первую очередь, этот переход приводит к необходимости пересмотра собственных частот колебаний системы, что сразу ведет к существенным отличиям поведения системы в режиме совместных ко- лебаний от системы в режиме заданного движения несущего тела. Важным элементом такого изменения частот является не только изменение их величин, но и смена оче- редности расположения частот при расположении частот по степени их возрастания (а значит и приоритета в проявлении резонансных свойств соответствующих им форм колебаний). Проведенные исследования для случаев поступательного и вращательно- го движения несущего тела с жидкостью [7 – 14] показывают, что изменения частот могут быть существенными. Так, например, частоты колебаний по первой антисим- метричной форме колебаний для случая поступательного движения несущего тела могут увеличиться в 2-2,5 раза, а для случая углового движения несущего тела на ма- ятниковом подвесе – в 4-5 раз, что фундаментально меняет диапазон и форму проте- кания околорезонансных процессов. 106 В настоящей работе показан характер изменения частот свободных колебаний ре- зервуаров с жидкостью на основе модели совместного движения составляющих ком- понент системы для разных вариантов движения несущего тела и показано, что в си- стеме при определенных вариантах маятникового закрепления несущего тела (кото- рые исследуются значительно сложнее, чем случаи поступательного движения несу- щего тела) наблюдается проявление вторичных резонансов, под которыми понимают- ся резонансы, устанавливающиеся на основе нелинейных механизмов, которые не проявляются вообще в рамках линейного моделирования системы. 1. Определение частот собственных колебаний при разных способах закреп- ления несущего тела. Дальнейшие исследования проведены на основе модели резервуара как абсолют- но твердого тела, а жидкость рассматривается в рамках модели идеальной несжимае- мой жидкости. При этом движение жидкости предполагается безвихревым. В рамках таких предположений можно ввести потенциал скоростей жидкости, что позволяет определить дискретный спектр колебаний жидкости и далее представлять ее движе- ние в виде разложения по формам собственных колебаний. Такой подход к построе- нию нелинейной модели совместного движения системы был развит в работах [2] и на его основе было рассмотрено большой круг теоретических и прикладных задач [7 – 9, 13 – 15]. В дальнейших примерах мы ограничимся случаем вертикального кругового цилиндрического резервуара. Частотные параметры систем будем определять на основе уравнений движения системы для линейной модели [8], учитывающей только первую антисимметричную форму колебаний свободной поверхности жидкости, движение которой характеризу- ется амплитудой 1a . Поступательное движение системы вдоль оси Ox характеризует- ся амплитудой отклонения от начального состояния x , а угловое движение вокруг оси Oy характеризуется углом 1 . Система уравнений совместного движения в рам- ках таких предположений примет вид 1 1* 2 1 1 1 11 1 1 1 1 1 1 0; 2x xv v a B E a          1 1 1 1 1 2 0;x x r l r l a B F M M M M          (1)    1* 1 11 2 1 11 2 1 11 1 1 2 2 2 / 4 / 2 0.x res r x l x g a E F J E M R H M R H                      Здесь остальные обозначения представляют собой коэффициенты уравнений движе- ния, определяемые через квадратуры от форм колебаний свободной поверхности жидкости и потенциала Стокса – Жуковского. Также обозначено  – плотность жид- кости; 11 resJ – момент инерции резервуара относительно оси Oy, проходящей через центр невозмущенной свободной поверхности жидкости; rM и lM – массы резерву- ара и жидкости; H – глубина заполнения резервуара; xR – длина маятникового под- веса, отсчитываемая от точки подвеса к центру невозмущенной свободной поверхно- сти жидкости. Если предположить, что в системе движение происходит без учета совместного характера движения составляющих компонент, то можно получить следующие парци- альные частоты 1 1tanhp a g g R R    – для жидкости; 0p a  – для поступательного движения несущего тела и 11 p r r l l a res M h M h J    – для углового (маятникового) движе- 107 ния несущего тела (здесь 1 – волновое число для первой антисимметричной формы колебаний свободной поверхности жидкости, g – ускорение свободного падения, R – радиус цилиндрического резервуара, rh и lh положения центров масс резервуара и жидкости относительно центра невозмущенной свободной поверхности жидкости). Однако для определения собственных частот колебаний системы в совместном дви- жении резервуара с жидкостью следует воспользоваться таким характеристическим уравнением для определения частот  , вытекающим из системы (1)     2 2 1 2 1* 2 1 1 11 1 1 1 2 2 1 2 1 2 1* 2 1 2 11 2 2 11 2 11 1 1 2 0. 1 2 2 2 / 4 / 2 xv v x r l r l res r x l x B E B F M M M M g E F J E M R H M R H                                Приведем некоторые численные результаты определения частот совместных ко- лебаний, полученные для случая 0,1r lM M для разных длин маятникового подвеса резервуара и 1R  м. В таблице в первой колонке приведены значения длин подвесов резервуара, во второй колонке приведена частота парциальных колебаний по первой антисимметричной форме, в третьей колонке приведена парциальная частота маятни- ковых колебаний резервуара для случая неподвижной точки подвеса, в четвертой ко- лонке приведены частоты совместных колебаний резервуара и жидкости для случая точки подвеса маятника, движущейся по заданному закону (верхнее значение соответ- ствует движению близкому к маятниковым колебаниям, нижнее – близкому к колеба- ниям свободной поверхности жидкости), в пятой колонке представлены результаты для полностью связанной задачи с подвижной точкой подвеса, движение которой наперед неизвестно (верхнее значение соответствует движению близкому к маятнико- вым колебаниям, нижнее – близкому к колебаниям свободной поверхности жидкости). 1 2 3 4 5 1l R 4,14431 2,47842 2,24312 2,26729 6,69986 16,63085 2l R 4,14431 1,96189 1,85304 1,86112 6,10821 12,45135 8l R 4,14431 1,07442 1,05895 1,05934 5,52733 9,75099 Как видно из результатов табличных данных частоты совместных колебаний си- стемы для маятниковых колебаний имеют тенденцию к незначительному убыванию, а частоты колебаний жидкости сильно возрастают (от 23% до 400%), что соответствует результатам теоремы об изменении частот в механических системах при переходе к связанным колебаниям. Такое существенное изменений частот колебаний жидкости в режиме связанных колебаний в системе показывает, что игнорирование учета фактора совместности колебаний компонент системы является недопустимым. Известно [3, 4, 8, 13, 14], что при совместном движении резервуара и жидкости связанное движение устанавливается только между параметрами движения резервуа- ра и теми формами, которые приводят к изменению положения центра масс объема жидкости в ее возмущенном движении. Для кругового цилиндрического резервуара это будут только антисимметричные формы (окружной номер 1),m  для которых и проявляется существенное изменение частот по сравнению с их парциальными значе- 108 ниями. В тоже время частоты колебаний жидкости для форм с окружным номером, отличным от 1, не будут меняться при переходе к режиму движения системы с учетом совместности движения ее компонент. Поэтому существенное изменение частот сов- местных колебаний антисимметричных форм может привести к изменению порядка следования форм колебаний в режиме совместных колебаний по сравнению со случа- ем заданного движения несущего тела, если располагать формы колебаний в порядке возрастания частот. Проиллюстрируем это на рис. 1, где для случая заданного движе- ния точки подвеса и длины маятникового подвеса l R представлена тенденция из- менения частот от парциальных к совместным с указанием их расположения на чис- ловой оси. Рис. 1 На рисунке парциальные частоты, соответствующие маятниковым колебаниям ре- зервуара и колебаниям свободной поверхности жидкости по первой антисимметрич- ной форме обозначены звездочкой. При переходе к моделированию в рамках модели совместного движения новые значения частот становятся 2,24 вместо 2,47 (уменьше- ние минимальной частоты) и 6,69 вместо 4,14 (увеличение максимальной частоты). Как видно из рисунка, теперь частота колебаний свободной поверхности жидкости по первой антисимметричной форме становится уже не второй по порядку возрастания частот (после маятниковой частоты), а четвертой, превосходя частоты колебаний сво- бодной поверхности жидкости по формам колебаний с окружными номерами 2 и 0. Естественно ожидать ослабления проявления ее резонансных свойств по сравнению со случаем заданного движения несущего тела, который соответствует частоте 4,14. Приведенный рисунок соответствует длине подвеса l R . Для больших длин подвеса (как это следует из табл. 1) увеличение частоты колебаний свободной поверхности по антисимметричной форме будет проявляться слабее. Соответственно, в зависимости от длины маятникового подвеса поведение системы «резервуар – жидкость со сво- бодной поверхностью» можно разбить на такие три диапазона. Первый диапазон (большие длины подвеса), когда частота антисимметричной формы колебаний будет меньше 5,45 и будет занимать второе место после частоты маятниковых колебаний. Второй диапазон (средние длины подвеса), когда частота антисимметричной формы колебаний будет находиться в интервале от 5,45 до 6,12, то есть будет на третьем ме- сте. И третий диапазон (короткие длины подвеса), когда частота антисимметричной формы колебаний будет больше 6,12 (четвертое место в порядке возрастания частот). Исходя из известных свойств проявления резонансных свойств системы в зависимо- сти от положения резонансных частот в общем ряду собственных частот колебаний следует ожидать разные тенденции в поведении системы для разных длин подвесов резервуара. Заметим при этом, что в рамках линейной теории резонансы для форм колебаний с окружными номерами 2m  и 0m  вообще не проявляются. Поэтому остается открытым вопрос о поведении системы в рамках нелинейной модели совмест- ного движения. 2. Нелинейная математическая модель совместного движения системы. Моделирование резонансных свойств системы «резервуар – жидкость» будем производить на основе математической модели, выведенной в [3], прошедшей много- стороннюю апробацию, включая сравнение с качественными результатами экспери- ментов [2, 7, 8, 9, 13 – 15]. За независимые параметры, характеризующие движение системы, приняты амплитуда возбуждения форм колебаний свободной поверхности жидкости ia , параметры поступательного i и вращательного i движения несущей конструкции. Согласно теореме о том, что безвихревое движение идеальной однород- 109 ной жидкости полностью определяется движением ее границ, выбранные параметры полностью характеризуют динамику системы, а именно, ia определяют движение свободной поверхности, а i и i – движение твердых границ области, занимаемой жидкостью, При этом по этим параметрам полностью восстанавливаются характери- стики движения свободной поверхности жидкости  , поле скоростей жидкости и уг- ловая скорость движения несущего тела i . Математическая модель системы «резервуар – жидкость» представляет собой не- линейную систему обыкновенных дифференциальных уравнений второго порядка, в которой вторые производные неизвестных входят линейно, что создает предпосылки для аналитического приведения рассматриваемой системы уравнений к форме Коши, удобной для последующего численного интегрирования. Система уравнений движе- ния с обозначениями, приведенными в [2] преобразована к форме, соответствующей движению точки маятникового подвеса по заданному закону )(t , имеет такой вид: 3 3 3 4 1* 2* 3* , 1 1 , 1 2 p i ir j rij j k rijk s pr i pri i j prijv i j j k s p i i jr s a a A a a A E a E a a E                                   3 4 1 2 3 4 , , , , , , 1 i j ijr i j k ijkr r i ri i j rij i j k rijkv i j i j k i i j i j kr a a C a a a C B a B a a B a a a B                               3 2* 2* 3* 3* 3* 3* 1 , 1 2 p i pir pri i j pijr pirj prij prjiv p i i jr a E E a a E E E E                     3 3 2 3 3 ( ) 1* 2* 3* , 1 1 , 1 1 2 2 k p s psr i psir psri p pr i pri i j prijv v p s i p i i jr r E a E E E a E a a E                           3 2 3 4 2 3 3 1 1 2i ir i j ijr i j k ijkr p pr i pir priv i i i p ir a D a a D a a a D F a F F                                      4 4 4 1 2 3 1 3 , cos sin cos sin sin c r i j pijr pirj prij v i j r a a F F F g                     1 2 1 3 1 2 3cos cos sin cos cos sin sin s r r rv v r r N g a g             ; 3 1* 2* 3* 1 , p i pi j pij j k pijk i p j j kr a E a E a a E                    3 3 2 2 3 1 , 1 , 1 2 p p ps resn ps i psi i j psij n p s i i jr r A a E a a EJ                              3 * ( ) 2 2 3 , , 1 , 1 2 pk ps resp r s p ps i psi i j psij p s i i jr A a E a a EJ                         3 1 2 3 4 1 , 2 p p i pi i j pij i j k pijk p i i j ir F a F a a F a a a F                        110                        3 1 , 321* , 3 1 , *3 , *2*1* , 2 p i ji pijjipiiprp p i ji pijkkji ji pijjipiirp FaaFaFEaaaEaaEa                              3 1 , *3*2 3 1, , 32 222 p ji pijkkji i pijji r p sp ji psijji i psiis r p EaaaEaaEaaEa              3 2 3 1 2 1 , 2 2 2 cos cosp i pi i j pij T T F F p i i jr r g a F a a F M h M h                              1 2 3 1 32 cos sin cos sin sin c c i i ir g a Hl                   1 2 3 1 3cos sin sin sin cos s s i i i a Hl               . (2) Система уравнений движения содержит две группы связанных уравнений: N урав- нений для определения обобщенных амплитудных параметров движения свободной поверхности жидкости ia ; три уравнения для определения параметров углового дви- жения конструкции. В целом система содержит 3N  уравнений, которые описывают совместное движение нелинейной системы «конструкция – жидкость». В систему уравнений (2) входят индексные выражения, характеризующие дина- мические связи в системе, инерциальные и статические характеристики конструкции, а также внешние моменты, действующие на конструкцию. При этом все эти коэффи- циенты определяются в квадратурах от форм колебаний свободной поверхности жид- кости и потенциалов Стокса – Жуковского [2]. Далее для численных примеров рас- сматривается модель движения жидкости, учитывающая 12N  форм колебаний сво- бодной поверхности жидкости. 3. Особенности проявления вторичных резонансов при колебаниях резервуа- ра с жидкостью на маятниковом подвесе. В системе «резервуар – жидкость со свободной поверхностью на маятниковом подвесе с точкой подвеса», движущейся по заданному закону, основным резонансом является резонанс на частоте близкой к частоте колебаний маятникового подвеса без учета подвижности жидкости. В режиме совместных колебаний собственной формой колебаний будет вид движения, включающий в себя не только маятниковые колеба- ния резервуара, но и противофазные колебания свободной поверхности жидкости по антисимметричной форме. Такой резонанс наиболее чувствителен к кинематическим возмущениям движения точки подвеса. Для коротких и средних длин подвесов следующей по порядку идет частота коле- баний 5,45, соответствующая колебаниям по форме с окружным номером 2.m  Та- кой резонанс не проявляется в рамках модели, не учитывающей совместность движе- ния и в рамках линейной теории. Рассмотрим колебания системы при возмущении движения точки подвеса маятника по закону sinx A t  . Для длины подвеса l R , 0,0022A  и 5,4571  получим следующую картину развития процессов в системе. На рис. 2 представлены законы изменения возмущения свободной поверхности жид- кости на стенке резервуара  и угла отклонения маятника  во времени. Обращает на себя внимание существенный рост амплитуд колебаний жидкости и угловых возмущений в системе, а также проявление модуляции колебаний. Законы изменения этих параметров не выходят на установившиеся значения с ростом време- ни. Анализ изменения внутренних параметров системы показывает, что рост ампли- туд колебаний жидкости на стенке резервуара в основном происходит из-за роста ам- 111 плитуд антисимметричных форм колебаний ( 1),m  хотя сами колебания происходят на частоте колебаний формы 2m  . При этом амплитуды колебаний формы 2m  изменяются незначительно (до 0,025R). Фактически форма колебаний 2m  играет в этом случае роль своеобразного катализатора проводника энергии от поступательного движения точки подвеса к форме колебаний 1,m  а сама возбуждается незначительно. Рис. 2 Анализ развития вторичного резонанса на частоте формы колебаний 2m  пока- зал, что с увеличением длины подвеса его проявление ослабевает, а для поступатель- ного движения резервуара он не проявляется вообще. Исследования показали, что наблюдаемый в таком резонансе энергообмен полностью определяется коэффициен- тами уравнений движения 1* piE , 2* pijE , 3* pijkE , которые отвечают за связь волнового дви- жения жидкости и углового движения резервуара. В частности, этими коэффициента- ми определяется также действие сил Кориолиса на частицы жидкости. Естественно, вклад таких членов при поступательном движении несущего тела не происходит, что и объясняет отсутствие такого резонанса для этого варианта движения несущего тела. Рассмотрим теперь случай проявления внутреннего резонанса системы. Подберем длину маятникового подвеса таким образом, чтобы частоты форм колебаний для 1m  и 2m  совпали. Это происходит для длины маятникового подвеса 12l R . При таких же параметрах движения точки маятникового подвеса получим аналогич- ный закон изменения возмущений свободной поверхности жидкости на стенке резер- вуара во времени (рис. 3). Рис. 3 Анализ изменения амплитуд отдельных форм колебаний показывает, что основ- ной вклад в возмещения свободной поверхности жидкости на стенках резервуара дает форма колебаний 1,m  а возмущение других форм на порядок меньше. Заметим, 112 однако, что колебания по осесимметричной форме 3a ( 0)m  развиваются таким об- разом, что амплитуда этой формы 3a всегда больше нуля (рис. 4), что обеспечивает несимметричность профилей волн, когда высота гребня волны выше глубины впади- ны, подтверждаемую экспериментально. Рис. 4 В целом настройка системы на условия внутреннего резонанса 12l R не внесла существенных изменений в развитие процессов по сравнению со случаем l R . Таким образом, показано, что в системе «резервуар – жидкость на маятниковом подвесе» в широком диапазоне изменения длин подвеса проявляется резонанс по форме 2,m  который в рамках линейной модели или модели, не учитывающий сов- местность движения составляющих компонент, не проявляется вообще. Важно также отметить, что в случае поступательного движения несущего тела такого типа резонанс для формы колебаний 2m  не проявляется вообще. Рассмотрим теперь специфику проявления резонансных свойств формы колеба- ний m=0, которая для коротких подвесов является третьей частотой при их размеще- нии в порядке возрастания значений частот. Анализ численных примеров показывает, что практически для всех значений коротких и средних длин подвесов этот резонанс не проявляется. Однако в режиме настройки на внутренний резонанс между формами колебаний, соответствующих 0m  и 1m  удалось обнаружить проявление резо- нанса. Для значений 1,95l R (резонансная длина подвеса), амплитуды колебаний точки подвеса маятника 0,01A  и частоты 6,1261  получена такая картина воз- мущения волн на стенке резервуара (рис. 5). Рис. 5 113 Причем анализ вклада амплитуд форм колебаний в формирование этого возмуще- ния свободной поверхности жидкости на стенке показывает, что примерно до 200 с это возмущение формируется в основном за счет первой антисимметричной формы колебаний. Однако далее происходит вовлечение в процесс колебаний первой осе- симметричной формы, что видно из графика (рис. 6). Рис. 6 Наличие такого резонанса подтверждается также спектральной характеристикой системы, построенной за 300 с колебаний (рис. 7). Из этой диаграммы видно наличие двух ярко выраженных пиков для частот 5,45 (форма колебаний 2)m  и 6,12 (сов- мещенная частота, соответствующая формам 1m  и 0).m  Причем из диаграммы следует, что резонанс на частоте 6,12 энергетически выражен сильнее. Рис. 7 Таким образом, резонансные свойства третьей по порядку формы колебаний (осе- симметричная форма), которая в рамках линейного моделирования или в рамках мо- дели, не учитывающей совместный характер движения компонент, не проявляются вообще, выражены значительно слабее, чем маятниковые колебаний системы и коле- бания в окрестности резонанса по форме 2.m  Только настройка системы на внут- ренний резонанс с формой колебаний 1m  создает условия для проявления этого нелинейного резонанса. Наконец, для формы колебаний 1m  (первая антисимметричная форма колеба- ний) для коротких подвесов, которая является четвертой в последовательности раз- мещения форм колебаний по степеням роста их собственных частот, обнаружить про- явление резонансных свойств практически не удается. В совокупности получается, что применение коротких и средних маятниковых подвесов позволяет избежать резо- 114 нансов на частоте первой антисимметричной формы, которая в большинстве других случаев является основной при резонансных колебаниях. Важно отметить общее свойство для всех типов резонансов, проанализированных в данной работе, а также для резонансов по маятниковой форме колебаний [13, 14]: для всех режимов нет асимптотического выхода на постоянное значение амплитуды колебаний свободной поверхности жидкости, однако сильно проявляется модуляция колебаний. Это отличие от результатов работ [5, 16] в последнее время получило ши- рокое экспериментальное подтверждение [5, 19 – 21] и качественное обоснование на основе парадокса Ферми – Паста – Улама [1], в частности, с применением к многоча- стотным волновым процессам на поверхности идеальной несжимаемой жидкости [18], где указано, что в многомодовых волновых процессах с трансцендентным соот- ношением частот волновой процесс никогда не выйдет на периодический режим ко- лебаний. Именно многочастотный характер развивающихся процессов приводит к необходимости исследования нелинейных задач динамики резервуаров с жидкостью на основе многомодовых моделей с учетом совместного движения компонент систе- мы, что в свою очередь может сильно преобразовать частотные характеристики си- стемы. Заключение. На основе теоретического моделирования резонансных процессов исследованы закономерности развития угловых колебаний резервуара и частично заполняющей его жидкости в рамках нелинейной модели, учитывающей совместный характер движе- ния составляющих компонент системы. Показано, что учет совместности движения компонент привел к существенному изменению величин частот и порядка следования форм колебаний системы при их размещении в порядке возрастания частот. Установ- лено, что для коротких и средних длин маятниковых подвесов первая антисиммет- ричная форма, которая обычно является главной в развитии резонансных процессов, теперь возбуждается не на своей резонансной частоте, а на частоте маятниковых ко- лебаний системы и частотах вторичных (нелинейных) резонансов на основе нелиней- ного перераспределения энергии. Показано соответствие развития динамических про- цессов в системе основным закономерностям нелинейного волнообразования, уста- новленным в теоретических и экспериментальных работах. Р Е З Ю М Е . Розглянуто задачу про коливання резервуара на маятниковому підвісі з рухомою точкою підвісу, частково заповненого ідеальною нестисливою рідиною. Задачу розглянуто в неліній- ній сумісній постановці. Показано, що розгляд сумісності руху компонент системи призводить до суттєвої зміни величин частот и порядку розташування форм коливань системи при їх розміщенні в порядку зростання частот. Встановлено, що для коротких і середніх довжин маятникових підвісів перша антисиметрична форма, яка зазвичай є головною в розвитку резонансних процесів, тепер збу- джується не на своїй резонансній частоті, а на частоті маятникових коливань і частотах вторинних (нелінійних) резонансів на основі нелінійного перерозподілу енергії. Показано відповідність розвит- ку динамічних процесів в системі основним закономірностям нелінійного хвилеутворення, встанов- леним в теоретичних і експериментальних роботах. 1. Заславский Г.М., Сагдеев Р.З. Введение в нелинейную физику: От маятника до турбулентности и хаоса. – Москва: Наука, 1988. – 368 с. 2. Лимарченко О.С., Ясинский В.В. Нелинейная динамика конструкций с жидкостью. – Киев: НТТУ КПИ, 1997. – 338 с. 3. Микишев Г.Н. Экспериментальные методы в динамике космических аппаратов. – Москва: Маши- ностроение, 1978. – 247 с. 4. Abramson H. The dynamic behavior of liquid in moving container. – Washington: NASA SP–106, 1966. – 467 p. 5. Faltinsen O.M., Rognebakke O.M., Timokha A.N. Transient and steady-state amplitudes of resonant three- dimensional sloshing in a square base tank with a finite fluid depth // Physics of Fluids. – 2006. – 1, N 18. – P. 1 – 14. 6. Ibrahim R.A. Liquid sloshing dynamics: theory and applications. – Cambridge: Cambridge University Press, 2005. – 950 p. 115 7. Konstantinov A.V., Limarchenko O.S. Effect of the Viscosity and Capillarity of Fluid on the Nonlinear Dynamics of a Tank Partially Filled with a Fluid // Int. Appl. Mech. – 2017. – 53, N 2. – P. 130 – 138. 8. Konstantinov A.V., Limarchenko O.S., Lukianchuk V.V., Nefedov A.A. Dynamic Methods of Damping the Oscillation in Structure-Free-Surface Fluid System // Int. Appl. Mech. – 2019. – 55, N 1. – P. 58 – 67. 9. Konstantinov A.V., Limarchenko O.S., Mel’nik V.N., Semenova I.Yu. Generalizing the Faraday Problem of the Parametric Oscillations of a Cylindrical Tank Partially Filled with a Fluid // Int. Appl. Mech. – 2017. – 53, N 1. – P. 59 – 66. 10. Kubenko V.D., Koval’chuk P.S. Modeling the Nonlinear Interaction of Standing and Traveling Bending Waves in Fluid-Filled Cylindrical Shells Subject to Internal Resonances // Int. Appl. Mech. – 2014. – 50, N 4. – P. 353 – 364. 11. Kubenko V.D., Koval’chuk P.S. Stability and Nonlinear Vibrations of Closed Cylindrical Shells Interact- ing with a Fluid Flow (Review) // Int. Appl. Mech. – 2015. – 51, N 1. – P. 12 – 63. 12. Kubenko V.D., Koval’chuk P.S. Nonlinear Vibrations of Fluid-Filled Cylindrical Shells with Initial Im- perfections // Int. Appl. Mech. – 2017. – 53, N 2. – P. 121 – 129. 13. Limarchenko O., Nefedov A. Resonant modes of the motion of a cylindrical reservoir on a movable pen- dulum suspension with a free-surface liquid // Mathematical Modeling and Computing. – 2018. – 5, N 2. – P. 178 – 183. 14. Limarchenko O.S., Semenovich K.O. Energy redistribution between the reservoir and liquid with free surface for angular motions of the system // J. of Mathematical Sciences. – 2017. – 222, N 3. – P. 296 – 303. 15. Limarchenko O.S., Tkachenko R.V. Influence of Spring Attachment on the Dynamics of a Fluid-Filled Cylindrical Tank on a Moving Platform // Int. Appl. Mech. – 2015. – 50, N 3. – P. 289 – 294. 16. Lukovsky I. A. Nonlinear Dynamics. Mathematical Models for Rigid Bodies with a Liquid. – Berlin: De Gruyter, 2015. – 410 p. 17. Miles J. Resonantly forced surface waves in a circular cylinder // J. Fluid Mech. – 1984. – 149. – P. 15 – 31. 18. Onorato M., Vozella L., Proment D., Lvov V. Route to thermalization in the α-Fermi – Pasta – Ulam system // Proc. Nat. Acad. Scie. of USA. – 2015. – 112, N 14. – P. 4208 – 4213. 19. Pal P. Sloshing of liquid in partially filled container – an experimental study // Int. J. of Recent Trends in Engineering. – 2009. – 1, N 6. – P. 1 – 5. 20. Shaoa W., Yanga J., Hu Z, Tao L. Coupled analysis of nonlinear sloshing and ship motions // Appl. Ocean Research. – 2015. – 47. – P. 85 – 97. 21. Zhang C., Li Y., Meng Q. Fully nonlinear analysis of second-order sloshing resonance in a three- dimensional tank // Computers and Fluids. – 2015. – 116. – P. 88 – 104. Поступила 09.04.2019 Утверждена в печать 09.07.2020