Numerical modeling of natural vibrations of coaxial shells partially filled with fluid, taking into account the effects on the free surface

Abstract


This work is devoted to the numerical analysis of vertical elastic coaxial cylindrical shells completely or partially filled with a quiescent compressible fluid with account for sloshing of its free surface. The problem is solved in the axisymmetric formulation using a semi-analytical version of the finite element method. The fluid medium is described by the wave equation, which together with the conditions prescribed at its boundaries is reduced to the weak form by the Bubnov - Galerkin method. A mathematical formulation of the dynamic problem for thin-walled structures is developed using the variational principle of virtual displacements and linear theory of thin shells based on Kirchhoff - Love hypotheses. The fluid pressure on the walls of the structure is calculated according to Bernoulli's equation. Sloshing modes caused by gravitational effects on free surface of liquid medium are excluded from the resulting system of equations through the use of the iteration dynamic condensation method. The numerical model has been verified by comparing it with the known data for the case of a single shell partially filled with fluid. The influence of the fluid level on the lowest natural frequencies of the system vibrations at different variants of kinematic boundary conditions for shells (rigid clamping at both edges, cantilever support) and different widths of the annular gap between them has been evaluated. It has been found that for the examined configurations, the level of the fluid in the annular channel has a stronger influence on the frequency spectrum compared to the level of fluid in the cavity of the inner shell due to the fact that vibration frequencies change in a wider range.

Full Text

Введение Элементы машиностроительных конструкций в виде коаксиальных оболочек, кольцевой зазор между которыми и реже - внутренний объём заполнены текучей жидкостью, уже значительное время являются объектами разнообразных теоретических и экспериментальных исследований в связи с их широким распространением в различных областях техники [1]. Подробные библиографические обзоры, отражающие различные аспекты взаимодействия вложенных одна в другую соосных и имеющих эксцентриситет оболочек с неподвижной или текущей жидкостью, приведены в недавних статьях авторов [2-6]. Во введении данной работы будут процитированы только те немногочисленные публикации, в которых неподвижная жидкость располагается в обеих полостях, а упругая конструкция, как правило, имеет вертикальное расположение. В статье [7] с использованием метода Бубнова - Галёркина получены простые расчётные формулы для определения частот колебаний коаксиальных оболочек с разными радиусами и разной высотой жидкости. В численном примере оценивается влияние жёсткости оболочек на спектр колебаний конструкции, полностью заполненной жидкостью. Совокупность вложенных одна в другую оболочек, полностью заполненных несжимаемой жидкостью, анализируется в [8]. Отсутствие жидкости во внутренней полости рассматривается как частный случай. При задании условия плескания на верхней границе столба жидкости делается предположение, приводящее практически к тем же результатам, что и при его отсутствии. Бесконечно длинные оболочки исследуются в статьях [9; 10]. В первой из них предложен простой и точный метод вычисления собственных частот колебаний полностью заполненных коаксиальных оболочек. В работе [10] также для случая полного заполнения обеих полостей установлено существование синфазных (направление и количество меридиональных полуволн совпадает для обеих оболочек) и противофазных (направления противоположны) форм колебаний. Численный анализ линейных колебаний оболочек, частично заполненных тяжёлой жидкостью, выполнен в работе [11] с применением метода конечных разностей. Приведены результаты для низших частот колебаний, полученных для разных отношений радиусов оболочек, полностью заполненных жидкостью. В публикациях [12; 13] приводятся постановка и решение задачи, учитывающие наличие жидкости в обеих полостях, но численные примеры сведены к частному случаю, когда жидкость полностью или частично заполняет только кольцевой канал. В аналитическом исследовании [14] рассматривается несоосное расположение оболочек, полностью заполненных жидкостью. Решение волнового уравнения, описывающего поведение жидкости в обеих полостях, осуществляется с помощью функций Бесселя. С этой целью в рассмотрение вводятся смещённые координаты, преобразование к которым выполняется посредством теоремы Бельтрами. Оценено влияние эксцентриситета на частоты колебаний при наличии подвижной среды в обеих полостях или только в кольцевом зазоре. Конечно-элементное решение, основанное на использовании изопараметрических элементов, представлено в исследовании [15]. При фиксированной величине кольцевого зазора продемонстрировано влияние различия уровней жидкости внутри полостей на частоты колебаний слоистых оболочек при различных углах армирования. Влияние свойств вязкой несжимаемой жидкости, заполняющей коаксиальные оболочки, на амплитуды и скорости волн продольной деформации изучено в работе [16] с использованием физически нелинейных соотношений. В перечисленных публикациях фактически отсутствуют работы, где в полной мере исследовано влияние как различного уровня жидкостей в обеих полостях, так и их размеров на частотный спектр конструкции. Наличие такого анализа будет способствовать более качественной оценке восприимчивости проектируемых объектов к внешним нагрузкам различной природы. При частичном заполнении резервуаров воздействие гравитационного поля приводит к эффектам волнообразования (плескания) на свободных поверхностях жидкости [17]. В ряде вышеупомянутых работ [7; 8; 11; 15] соответствующее граничное условие учитывается в исходной постановке проблемы. Однако в данных публикациях не уточняется, что при этом решение модальной задачи осложняется появлением значительного количества плескательных собственных частот, величина которых или существенно уступает, или сопоставима с низшими связанными частотами колебаний конструкции. Данный факт значительно усложняет анализ полученных результатов [18]. Влияние гравитационных эффектов на динамические характеристики упругих конструкций, как свободно колеблющихся, так и подверженных внешним нагрузкам, более тщательно изучено для одиночных оболочек, частично заполненных жидкостью [19-36]. В некоторых работах предложены конструктивные решения (упругие или жёсткие мембраны) с целью предотвращения этого явления, представляющего опасность при определённых условиях эксплуатации. В работе [23] продемонстрировано, что при неблагоприятной комбинации геометрических и физико-механических параметров плескательные моды становятся сопоставимыми по величине со связанными модами гидроупругих колебаний и, следовательно, их учёт становится обязательным при моделировании. На основании вышесказанного целесообразным является разработка такого численного алгоритма, в котором учитывается плескание свободной поверхности, но соответствующие моды колебаний исключаются (редуцируются) из общего спектра частот с целью повышения вычислительной эффективности и удобства обработки получаемых результатов. Использование техники понижения (редукции) исходной размерности разрешающей системы уравнений в рамках дискретных численных методов имеет длительную историю [37; 38]. За прошедшее время разработано значительное количество методов редуцирования, которые классифицируются по подходам, применяемым к построению координатной трансформации. Помимо этих устоявшихся и широко используемых алгоритмов существуют различные варианты редукции. Некоторые из них учитывают, например, нестандартные методы решения (комбинация методов конечных и граничных элементов и т.п.) или специфические особенности поставленной краевой задачи и используются при изучении гидроупругих колебаний с учётом плескания жидкости [39-46]. Проведённые численные эксперименты показали, что наиболее простые и экономичные методы, связанные с редукцией физических координат (статическая конденсация Гайана, IRS-метод, метод динамической конденсации), не позволяют добиться или полного исключения плескательных мод, или приемлемой точности вычислений. Искомая цель достигается при использовании более затратного в вычислительном плане метода итерационной динамической конденсации [37] или терминологически корректнее - итерационной схемы решения для метода динамической конденсации. С его помощью в данной работе выполнено исследование фундаментальных частот колебаний коаксиальных оболочек, внутренние полости которых полностью или частично заполнены неподвижной сжимаемой жидкостью. 1. Математическая постановка задачи и основные соотношения Рассматривается система (рис. 1), состоящая из вертикально ориентированных коаксиальных цилиндрических оболочек длиной , радиусами и толщинами Кольцевой зазор между оболочками и пространство внутренней частично заполнены идеальной сжимаемой жидкостью на высоты и соответственно. Здесь и далее верхний индекс «(1)» обозначает принадлежность к внутренней оболочке и заполняющей её полость жидкости, а «(2)» - к внешней оболочке и жидкости в кольцевом канале. Необходимо исследовать влияние уровней заполнения полостей соосных оболочек на низшие собственные частоты колебаний системы при различных значениях кольцевого зазора и разных вариантах кинематических граничных условий для упругих тел. Рис. 1. Расчётная схема Fig. 1. Computational scheme Поведение идеальной сжимаемой среды описывается волновым уравнением относительной потенциала скорости которое в цилиндрических координатах принимает вид [47] , (1) где - оператор Лапласа, - частная производная, - время, - скорость звука в жидкой среде. Здесь и далее нижние и верхние индексы «f» и «(i)» () определяют жидкости, заполняющие полости оболочек, а индексы «s» и «(j)» () - упругие тела. На смоченных поверхностях оболочек задаётся условие непроницаемости где и - поверхности, ограничивающие объёмы жидкости и оболочек - вектора перемещений оболочек и единичной нормали к их поверхностям. На жёстком днище оболочек для потенциалов скоростей принимается граничное условие: (3) Волны на свободной поверхности жидкости описываются согласно допущению, что её фактическое положение находится на некоторой высоте относительно средней поверхности . (4) Для получения уравнения плескания на свободной поверхности воспользуемся обобщённым уравнение Бернулли (5) где - давление, - ускорение свободного падения, - плотность жидкостей. Уравнение (5) с учётом условия на свободной поверхности (6) принимает вид: (7) Дифференцируя обе части равенства (1) по времени и подставляя результат в (4) окончательно получим (8) С помощью метода Бубнова - Галёркина уравнение (1) с учётом граничных условий (2), (3) и (8) преобразуется к слабой форме [48] (9) где и - пробные решения для потенциалов скорости и векторов перемещений оболочек, и - базисные функции и их количество. Гидродинамическое давление, действующее со стороны жидкости на оболочки, вычисляется с использованием уравнения Бернулли (10) Компоненты вектора деформации определяются в цилиндрической системе координат с помощью классической линейной теории тонких оболочек, основанной на гипотезах Кирхгофа - Лява [49] (11) где , и - меридиональная, окружная и радиальная компоненты вектора перемещений, В матричном виде соотношения, устанавливающие связь между вектором усилий и моментов и вектором деформаций записываются как (12) Компоненты матрицы для изотропного материала известным образом определяется через модуль Юнга и коэффициент Пуассона Математическая постановка задачи динамики тонкостенных упругих тел выполнена с использованием вариационного принципа возможных перемещений, который с учётом (2), (10) и работы сил инерции может быть представлен в следующем виде (13) где - векторы поверхностных нагрузок оболочек; - плотность материалов оболочек. 2. Численная реализация Численное решение задачи осуществляется с использованием полуаналитического варианта метода конечных элементов (МКЭ), при котором искомые переменные раскладываются в ряды Фурье по окружной координате q [50]. В результате чего исходная задача сводится к совокупности двумерных задач для каждой окружной гармоники j. Стандартные процедуры МКЭ позволяют получить из соотношений (9) и (10) следующую связанную систему уравнений, описывающую взаимодействие соосных оболочек с жидкостями в обеих полостях (14) где типовые матрицы масс демпфирования и жёсткости для отдельных конечных элементов формируются следующим образом: (15) Здесь и - функции формы для потенциала скорости жидкости, вектора перемещений оболочек и его нормальной составляющей; и - количество конечных элементов, на которые разбиты объёмы жидкостей, оболочек, а также смоченные и свободные поверхности соответственно; - матрица градиентов, определяющая связь деформаций с узловыми перемещениями. Решение системы (14) ищется в виде (16) где - некоторые векторные функции координат, - характеристический показатель, - величина, отвечающая за демпфирование системы, - собственная частота колебаний, i - мнимая единица. С учётом (16) система уравнений (14) преобразуется к виду (17) Для нахождения квадратичная проблема на собственные значения (17) сводится к обобщённой для несимметричных матриц удвоенной размерности (18) где - единичная матрица. При численной реализации для оболочек используется конечный элемент в виде усечённого конуса с аппроксимацией меридиональной и окружной компонент вектора перемещений линейным полиномом, а нормальной компоненты - кубическим. Дискретизация объёма жидкости проводится с помощью прямоугольного конечного элемента с линейной аппроксимацией потенциала скорости. В случае заполнения обеих полостей приемлемая точность вычислений достигается при размерности нередуцированной системы, равной 3198 степеням свободы. Численный алгоритм реализован в среде MATLAB с использованием возможностей пакета ANSYS для создания конечно-элементной сетки. 3. Редуцирование системы уравнений Исключение плескательных мод колебаний, связанных с гравитационными эффектами на свободной поверхности жидкости, осуществляется методом итерационной динамической конденсации [37]. С этой целью перепишем уравнение (17), выделив главные (удерживаемые) и ведомые (редуцируемые) степени свободы, обозначив их нижними индексами «m» и «s» соответственно (19) Степени свободы, относящиеся к перемещениям упругих тел и потенциалу скорости жидкости, за исключением расположенных на свободной поверхности, обозначим за главные. В качестве ведомых примем те, которые относятся только к свободной поверхности жидкости. Представляя систему (19) в виде (20) (21) и выражая одни переменные через другие (22) получим редуцированную систему уравнений (23) Матрица преобразования является частотно зависимой, поэтому для решения редуцированной системы (23) требуется использование описанной ниже итерационной процедуры. Для каждого подлежащего определению k-го собственного значения вычисляемого с относительной погрешностью выполняются следующие шаги: 1. Задаётся начальное приближение 2. На каждой l-й итерации согласно (22) вычисляется матрица преобразования 3. В соответствии с (23) формируются матрицы редуцированной системы 4. Из редуцированной системы уравнений, приведённой к обобщённой проблеме на собственные значения для несимметричных матриц удвоенной размерности, определяется собственное значение 5. Итерационная процедура завершается при выполнении условия , где и - приближения, вычисленные на текущей и предыдущей итерациях. В противном случае последовательность действий повторяется, начиная с шага 2. 4. Верификация модели / численного алгоритма В связи с отсутствием подходящих для сравнения данных, полученных для коаксиальных оболочек, верификация численного алгоритма осуществляется на вертикальной цилиндрической оболочке, частично заполненной жидкостью ( кг/м3, кг/м3, и свободно опёртой с двух торцов (SS, Аналитическое решение, полученное методом Рэлея - Ритца для указанной конфигурации, представлено в работе [24]. В таблице приведено сравнение собственных частот колебаний w, вычисленных для полной и редуцированной постановок, с результатами статьи [24]. Здесь через m обозначено количество полуволн в меридиональном направлении. Данные таблицы демонстрируют хорошее соответствие результатов, определяемых в рамках описанного алгоритма, с известными сведениями. Применение метода итерационной динамической конденсации позволяет полностью исключить плескательные моды колебаний, а различие в связанных частотах, вычисляемых в полной и редуцированной постановках, не превышает %. 5. Результаты Представление полученных результатов расчётов осуществляется с помощью следующих безразмерных переменных: уровень заполнения жидкостью Сравнение собственных частот колебаний w (рад/с) свободно опёртой цилиндрической оболочки, частично заполненной жидкостью Comparison of natural vibration frequencies w (rad/s) of a simply supported cylindrical shell partially filled with fluid Моды Плескательные частоты Связанные частоты j m [24] Полная система Редуцированная система [24] Полная система Редуцированная система 4 1 1.4425 1.4431 - 13.658 14.061 14.061 2 1.9081 1.9111 - 34.441 34.614 34.614 3 2.2305 2.2382 - 49.692 49.532 49.532 4 2.5027 2.5176 - 61.877 61.641 61.641 5 2.7444 2.7695 - 71.804 71.793 71.793 а б Рис. 2. Зависимости собственных частот колебаний w (Гц) от безразмерного уровня заполнения внутренней оболочки h1 при различных уровнях заполнения кольцевого канала h2 = 0,25 (а) и h2 = 0,50 (б), полученные для различных значений окружной гармоники j: k = 1/100; CF Fig. 2. The dependencies of natural frequencies w (Hz) on the dimensionless fluid level in the inner shell h1 at different levels of filling of the annular channel h2 = 0,25 (a) and h2 = 0,50 (б), obtained for different values of the circular harmonic j: k = 1/100; CF величина кольцевого канала В численных примерах рассматривается система жёстко (CC, ) или консольно (CF, ) закреплённых оболочек со следующими геометрическими и физико-механическими параметрами: кг/м3, кг/м3, На рис. 2 отображены зависимости частот колебаний w (Гц) для нескольких окружных гармоник j от уровня заполнения внутренней оболочки h1 при различных уровнях заполнения кольцевого канала h2. На рисунках приводятся гармоники, соответствующие частоты колебаний которых близки к минимальному значению для каждой из рассмотренных конфигураций. Полученные графики позволяют наглядно оценить степень эволюции частот для каждого варианта заполнения. Для каких-то конфигураций различная чувствительность окружных гармоник к присоединенной массе жидкости в кольцевом канале обеспечивает смену форм с низшей частотой колебаний, тогда как для других этот фактор не имеет решающего значения. На рис. 3, 4 приведены зависимости низшей собственной частоты (Гц) от уровня заполнения одного из каналов при фиксированной высоте жидкости в другом. Проанализированы оболочки с разными вариантами граничных условий при величине кольцевого зазора между оболочками, равной следующим значениям: k = 1/2; k = 1/10; k = 1/100. Из представленных результатов следует, что увеличение уровня заполнения в любой из полостей приводит к снижению значений минимальных частот колебаний системы вследствие возрастания присоединённой массы жидкости. Характерные изломы на графиках, в основном наблюдаемые при h2 = 0,25, информируют об изменении гармоники в окружном направлении j с минимальной частотой колебаний, как это продемонстрировано на рис. 2, а (кривые j = 2 и j = 3). Для рассмотренных конфигураций уровень заполнения кольцевого канала оказывает большее влияние на низшие частоты, поскольку в этом случае их значения изменяются в более широком диапазоне. При узком кольцевом зазоре колебательные характеристики системы могут определяться практически только уровнем жидкости в нём, как это имеет место для жёстко закреплённых оболочек (рис. 3, д). Представленные на рисунке зависимости позволяют сделать вывод, что уровень жидкости в полостях коаксиальных оболочек можно использовать для подбора частотного спектра конструкции, соответствующего требуемым условиям эксплуатации. Возможность управления значением низшей частоты колебаний посредством изменения уровня заполнения в одной из полостей снижается с повышением уровня жидкости в другой для всех рассмотренных комбинаций граничных условий и кольцевого зазора. Особенно сильно такое поведение проявляется при варьировании уровня заполнения полости внутренней оболочки при фиксированном уровне жидкости в кольцевом канале, когда последний достигает некоторого порогового значения. Представленные на рисунках данные демонстрируют также известную закономерность [1], проявляющуюся в снижении частот при повышении размеров внутренней оболочки. Заключение В работе представлены математическая постановка задачи и алгоритм её численной реализации, предназначенные для изучения колебательных характеристик системы коаксиальных оболочек, внутренние полости а б в г д е Рис. 3. Зависимости низших собственных частот колебаний w (Гц) от безразмерных уровней заполнения внутренней оболочки h1 (а, в, д) и кольцевого зазора h2 (б, г, е), CC: (а, б) - k = 1/2; (в, г) - k = 1/10; (д, е) - k = 1/100 Fig. 3. Dependencies of the lowest natural vibration frequencies w (Hz) on the dimensionless filling level of the inner shell h1 (а, в, д) and the annular gap h2 (б, г, е), CC: (а, б) - k = 1/2; (в, г) - k = 1/10; (д, е) - k = 1/100 которых полностью или частично заполнены идеальной сжимаемой жидкостью. На свободной поверхности жидкости принимаются во внимание эффекты плескания, обусловленные гравитационными силами. С целью повышения вычислительной эффективности плескательные моды колебаний исключаются путём редуцирования системы уравнений с помощью метода итерационной динамической конденсации. Использование такого подхода позволяет полностью исключить связанные с плесканием частоты и не вносит погрешность в спектр связанных частот. С использованием разработанного алгоритма проанализировано влияние уровней заполнения внутренних полостей соосных оболочек на фундаментальные частоты колебаний системы при разных размерах внутренней оболочки и граничных условиях. Установлено, что наибольшее изменение в спектре частот достигается путём варьирования высоты жидкости в кольцевом канале между оболочками. При этом ни комбинации кинематических граничных условий, задаваемых на краях конструкции, ни размеры кольцевого зазора, не оказывают качественного влияния на выявленные закономерности. а б в г д е Рис. 4. Зависимости низших собственных частот колебаний w (Гц) от безразмерных уровней заполнения внутренней оболочки h1 (а, в, д) и кольцевого зазора h2 (б, г, е), CF: (а, б) - k = 1/2; (в, г) - k = 1/10; (д, е) - k = 1/100 Fig. 4. Dependencies of the lowest natural vibration frequencies w (Hz) on the dimensionless filling level of the inner shell h1 (а, в, д) and the annular gap h2 (б, г, е), CF: (а, б) - k = 1/2; (в, г) - k = 1/10; (д, е) - k = 1/100

About the authors

S. A Bochkarev

Institute of Continuous Media Mechanics UrB RAS

S. V Lekomtsev

Institute of Continuous Media Mechanics UrB RAS

A. N Senin

Institute of Continuous Media Mechanics UrB RAS

References

  1. Païdoussis M.P. Fluid-Structure Interactions. Slender Structures and Axial Flow. - London: Academic Press, 2016. - Vol. 2. - 942 p. doi: 10.1016/s1874-5652 (04) x8001-7
  2. Бочкарёв С.А., Лекомцев С.В., Сенин А.Н. Анализ пространственных колебаний коаксиальных цилиндрических оболочек, частично заполненных жидкостью // Вычислительная механика сплошных сред. - 2018. - Т. 11, № 4. - С. 448-462. doi: 10.7242/1999-6691/2018.11.4.35
  3. Bochkarev S.A., Lekomtsev S.V., Senin A.N. Analysis of spatial vibrations of piezoceramic eccentric cylindrical shells interacting with an annular fluid layer // Frat. ed Integrità Strutt. - 2019. - Vol. 49. - P. 814-830. doi: 10.3221/IGF-ESIS.49.15
  4. Бочкарев С.А., Лекомцев С.В. Гидроупругая устойчивость коаксиальных цилиндрических оболочек, выполненных из пьезо-электрического материала // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2019. - № 2. - С. 35-48. doi: 10.15593/perm.mech/2019.2.04
  5. Hydroelastic stability of partially filled coaxial cylindrical shells / S.A. Bochkarev, S.V. Lekomtsev, V.P. Matveenko, A.N. Senin // Acta Mech. - 2019. - Vol. 230, No. 11. - P. 3845-3860. doi: 10.1007/s00707-019-02453-4
  6. Бочкарёв С.А., Лекомцев С.В., Сенин А.Н. Анализ гидроупругой устойчивости несоосных цилиндрических оболочек, частично заполненных жидкостью // Вестник Самарского государственного технического университета. Серия: Физико-математические науки. - 2020. - Т. 24, № 1. - С. 95-115. doi: 10.14498/vsgtu1746
  7. Балакирев Ю.Г. К исследованию осесимметричных колебаний соосных цилиндрических систем оболочек с жидким заполнителем // Инженерный журнал. МТТ. - 1968. - № 3. - С. 133-140.
  8. Пшеничнов Г.И. Собственные колебания коаксиальных ортотропных цилиндрических оболочек заполненных жидкостью // Труды 8-й Всесоюзной конференции «Теория оболочек и пластин», Ростов-На-Дону, 1971. - М.: Наука, 1973. - С. 546-549.
  9. Krajcinovic D. Vibrations of two coaxial cylindrical shells containing fluid // Nucl. Eng. Des. - 1974. - Vol. 30, № 2. - P. 242-248. doi: 10.1016/0029-5493 (74) 90169-1
  10. Chen S.S., Rosenberg G.S. Dynamics of a coupled shell-fluid system // Nucl. Eng. Des. - 1975. - Vol. 32, № 3. - P. 302-310. doi: 10.1016/0029-5493(75)90101-6
  11. Антонов В.Н. Колебания соосных цилиндрических оболочек, частично заполненных сжимаемой жидкостью // Известия АН СССР. МТТ. - 1977. - № 3. - С. 118-124.
  12. Швец Р.Н., Марчук Р.А. Осесимметричные колебания заполненных жидкостью соосных ортотропных цилиндрических оболочек // Прикладная механика. - 1985. - Т. 21, № 8. - С. 46-51.
  13. Fujita K. Vibration analysis of fluid-coupled two coaxial axisymmetric shells containing liquid // Bulletin of JSME. - 1986. - Vol. 29, № 248. - P. 516-524. doi: 10.1299/jsme1958.29.516
  14. Jeong K.-H., Lee G.-M., Chang M.-H. Free vibration analysis of a cylindrical shell eccentrically coupled with a fluid-filled vessel // Comput. Struct. - 2001. - Vol. 79, № 16. - P. 1517-1524. doi: 10.1016/S0045-7949(01)00031-1
  15. Okazaki K., Tani J., Sugano M. Free vibrations of a laminated composite coaxial circular cylindrical shell partially filled with liquid // Trans. Jpn. Soc. Mech. Eng. Part C. - 2002. - Vol. 68. - P. 1942-1949. doi: 10.1299/kikaic.68.1942
  16. Mogilevich L., Ivanov S. Longitudinal waves in two coaxial elastic shells with hard cubic nonlinearity and filled with a viscous incompressible fluid // Recent Research in Control Engineering and Decision Making. ICIT 2020. Studies in Systems, Decision and Control / Dolinina O. et al. (eds). - Springer, Cham. - 2021. - Vol. 337. - P. 14-26. doi: 10.1007/978-3-030-65283-8_2
  17. Ibrahim R.A. Liquid Sloshing Dynamics: Theory and Applications. - Cambridge: Cambridge University Press, 2005. - 970 p.
  18. Bochkarev S.A., Lekomtsev S.V., Senin A.N. Natural vibrations and stability of loaded cylindrical shells partially filled with fluid, taking into account gravitational effects // Thin Wall. Struct. - 2021. - Vol. 164. - 107867. doi: 10.1016/j.tws.2021.107867
  19. Кулешов В.Б., Швейко Ю.Ю. Неосесимметричные колебания цилиндрических оболочек, частично заполненных жидкостью // Изв. АН СССР. МТТ. - 1971. - № 3. - С. 126-136.
  20. Free vibration analysis of cylindrical liquid storage tanks / T. Balendra, K.K. Ang, P. Paramasivam, S.L. Lee // Int. J. Mech. Sci. - 1982. - Vol. 24, № 1. - P. 47-59. doi: 10.1016/0020-7403 (82)90020-0
  21. Chiba M., Yamaki N., Tani J. Free vibration of a clamped-free circular cylindrical shell partially filled liquid - Part I: Theoretical analysis // Thin Wall. Struct. - 1984. - Vol. 2. - P. 265-284. doi: 10.1016/0263-8231(84)90022-3
  22. Chiba M., Yamaki N., Tani J. Free vibration of a clamped-free circular cylindrical shell partially filled with liquid - Part II: Numerical results // Thin Wall. Struct. - 1984. - Vol. 3. - P. 307-324. doi: 10.1016/0263-8231 (84) 90002-8
  23. Lakis A.A., Neagu S. Free surface effects on the dynamics of cylindrical shells partially filled with liquid //j. Sound Vib. - 1997. - Vol. 207, № 2. - P. 175-205. doi: 10.1006/jsvi.1997.1074
  24. Amabili M., Païdoussis M.P., Lakis A.A. Vibrations of partially filled cylindrical tanks with ring-stiffeners and flexible bottom //j. Sound Vib. - 1998. - Vol. 213, № 2. - P. 259-299. doi: 10.1006/jsvi.1997.1481
  25. Jeong K.-H., Lee S.-C. Hydroelastic vibration of a liquid-filled circular cylindrical shell // Comput. Struct. - 1998. - Vol. 66, № 2-3. - P. 173-185. doi: 10.1016/S0045-7949(97)00086-2
  26. Amabili M. Eigenvalue problems for vibrating structures coupled with quiescent fluids with free surface //j. Sound Vib. - 2000. - Vol. 231, № 1. - P. 79-97. doi: 10.1006/jsvi.1999.2678
  27. Biswal K.C., Bhattacharyya S.K., Sinha P.K. Dynamic response analysis of a liquid-filled cylindrical tank with annular baffle //j. Sound Vib. - 2004. - Vol. 274, № 1-2. - P. 13-37. doi: 10.1016/S0022-460X (03) 00568-6
  28. Kim Y.-W., Lee Y.-S., Ko S.-H. Coupled vibration of partially fluid-filled cylindrical shells with ring stiffeners //j. Sound Vib. - 2004. - Vol. 276, № 3-5. - P. 869-897. doi: 10.1016/j.jsv.2003.08.008
  29. Lakis A.A., Bursuc G., Toorani M.H. Sloshing effect on the dynamic behavior of horizontal cylindrical shells // Nucl. Eng. Des. - 2009. - Vol. 239, № 7. - P. 1193-1206. doi: 10.1016/j.nucengdes.2009.03.015
  30. Biswal K.C., Bhattacharyya S.K. Dynamic response of structure coupled with liquid sloshing in a laminated composite cylindrical tank with baffle // Finite Elem. Anal. Des. - 2010. - Vol. 46, № 11. - P. 966-981. doi: 10.1016/j.finel.2010.07.001
  31. Оболочки вращения, частично заполненные жидкостью, при сейсмических и импульсных нагрузках / В.И. Гнитько [и др.] // Проблемы вычислительной механики и прочности конструкций. - 2012. - № 19. - С. 65-72.
  32. Бочкарёв С.А., Лекомцев С.В. Собственные колебания частично заполненных жидкостью некруговых цилиндрических оболочек с учётом плескания свободной поверхности // Вычислительная механика сплошных сред. - 2014. - Т. 7, № 4. - С. 471-480. doi: 10.7242/1999-6691/2014.7.4.45
  33. Колебания жидкости в упругих оболочках с разным уровнем их заполнения / Д.В. Клименко [и др.] // Космическая техника. Ракетное вооружение. - 2015. - Вып. 3. - С. 24-29.
  34. Kashani B.K., Sani A.A. Free vibration analysis of horizontal cylindrical shells including sloshing effect utilizing polar finite elements // Eur. J. Mech. A/Solids. - 2016. - Vol. 58. - P. 187-201. doi: 10.1016/j.euromechsol.2016.02.002
  35. Bochkarev S.A., Lekomtsev S.V, Matveenko V.P. Dynamic analysis of partially filled non-circular cylindrical shells with liquid sloshing // Int. J. Appl. Mech. - 2016. - Vol. 8, № 3. - 1650027. doi: 10.1142/S1758825116500277
  36. Zarifian P., Ovesy H.R., Firouz-Abadi R.D. Sloshing effects on supersonic flutter characteristics of a circular cylindrical shell partially filled with liquid // Int. J. Numer. Methods Eng. - 2019. - Vol. 117, № 8. - P. 901-925. doi: 10.1002/nme.5984
  37. Qu Z.-Q. Model order reduction techniques with applications in finite element analysis. - London: Springer, 2004. - 369 p.
  38. Игнатьев В.А., Ромашкин В.Н. Алгебраическая проблема собственных векторов и собственных значений высокого порядка в задачах динамики и устойчивости конструкций (обзор) [Электронный ресурс] // Интернет-вестник ВолгГАСУ. - 2015. - Вып. 2 (38). - URL: http://vestnik.vgasu.ru/? source=4&articleno=1946) (дата обращения: 10.10.2021).
  39. Felippa C.A., Ohayon R. Mixed variational formulation of finite element analysis of acoustoelastic/slosh fluid-structure interaction //j. Fluids Struct. - 1990. - Vol. 4. - P. 35-57. doi: 10.1016/0889-9746(90)90036-5
  40. Andrianarison O., Ohayon R. Reduced models for modal analysis of fluid-structure systems taking into account compressibility and gravity effects // Comput. Methods Appl. Mech. Eng. - 2006. - Vol. 195, № 41-43. - P. 5656-5672. doi: 10.1016/j.cma.2005.11.013
  41. Noorian M.A., Firouz-Abadi R.D., Haddadpour H. A reduced order model for liquid sloshing in tanks with flexible baffles using boundary element method // Int. J. Numer. Methods. Eng. - 2012. - Vol. 89. - P. 1652-1664. doi: 10.1002/nme.3304
  42. Schotté J.-S., Ohayon R. Linearized formulation for fluid-structure interaction: Application to the linear dynamic response of a pressurized elastic structure containing a fluid with a free surface //j. Sound Vib. - 2013. - Vol. 332, № 10. - P. 2396-2414. doi: 10.1016/j.jsv.2012.07.036
  43. Григорьев В.Г., Григорьева Е.В. Редукция поверхностных волн в конечно-элементном анализе колебаний упругих тел, содержащих жидкость // Проблемы машиностроения и надежности машин. - 2014. - № 6. - С. 11-16.
  44. Ohayon R., Soize C. Vibration of structures containing compressible liquids with surface tension and sloshing effects. Reduced-order model // Comput. Mech. - 2015. - Vol. 55, № 6. - P. 1071-1078. doi: 10.1007/s00466-014-1091-4
  45. Ghalandari M., Bornassi S., Shamshirband S., Mosavi A., Chau K.W. Investigation of submerged structures’ flexibility on sloshing frequency using a boundary element method and finite element analysis // Eng. Appl.Comput. Fluid Mech. - 2019. - Vol. 13, № 1. - P. 519-528. doi: 10.1080/19942060.2019.1619197
  46. Sotoudehni E., Shahabian F., Aftabi Sani A. A dynamic order reduction method for fluid-structure systems // Appl. Math. Model. - 2021. - Vol. 89, part 1. - P. 136-153. doi: 10.1016/j.apm.2020.06.071
  47. Вольмир А.С. Оболочки в потоке жидкости и газа. Задачи гидроупругости. - М.: Наука, 1979. - 320 с.
  48. Бочкарев С.А., Матвеенко В.П. Численное исследование влияния граничных условий на динамику поведения цилиндрической оболочки с протекающей жидкостью // Известия РАН. МТТ. - 2008. - № 3. - С. 189-199.
  49. Бидерман В.Л. Механика тонкостенных конструкций. Статика. -М.: Машиностроение, 1977. - 488 с.
  50. Зенкевич О. Метод конечных элементов в технике. - М.: Мир, 1975. - 544 с.

Statistics

Views

Abstract - 14

PDF (Russian) - 9

Cited-By


PlumX


Copyright (c) 2022 Bochkarev S.A., Lekomtsev S.V., Senin A.N.

Creative Commons License
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies