Numerical modeling of eccentric cylindrical shells partially filled with a fluid

Full Text

Введение. Содержащие жидкость или газ коаксиальные цилиндрические оболочки являются частью многих индустриальных приложений и используются в различных областях техники. Они уже на протяжении значительного времени выступают в качестве объекта разнообразных теоретических исследований. Обширная библиография, посвящённая анализу коаксиальных оболочек, взаимодействующих как с неподвижной, так и с текущей жидкостью, представлена в монографии [1]. В недавних публикациях авторов [2, 3] также упомянуты редко цитируемые, как правило, пионерские работы отечественных исследователей. Ниже приводится обзор статей, в которых рассмотрены различные аспекты динамического поведения взаимодействующих с жидкостью соосных и несоосных конструкций в линейной постановке. Бесконечно длинные коаксиальные оболочки, содержащие текущую как в кольцевом канале, так и во внутренней оболочке жидкость, исследованы в [4]. В аналитических моделях [5–6] движение оболочек конечной длины описывается с помощью теории тонких оболочек Флюгге. Обобщённые гидродинамические силы определяются в рамках потенциальной теории и вычисляются с помощью преобразования Фурье. Решение задачи осуществляется методом Галёркина. Полученные в [7] экспериментальные данные демонстрируют, что при жёсткой внешней оболочке потеря устойчивости осуществляется на значительно меньших скоростях течения, чем предсказывает аналитическая модель. Предполагается, что причиной таких расхождений являются несовершенства формы конструкции. Приближённая теория для оценки динамического поведения системы коаксиальных оболочек, взаимодействующих с несжимаемой жидкостью, предложена в [8]. Показано, что собственные частоты колебаний в случае неподвижной жидкости и критические скорости потери устойчивости хорошо согласуются с имеющимися численными результатами. Анализ устойчивости коаксиальных оболочек, взаимодействующих с двумя потоками идеальной сжимаемой жидкости при различных комбинациях граничных условий и свойствах материалов оболочек, выполнен в [9–10] с применением метода конечных элементов (МКЭ). В [11] решение трёхмерных линеаризованных уравнений Навье—Стокса ищется в виде суммы скалярного и векторного потенциалов, а системы коаксиальных оболочек, содержащих жидкость, — в классе бегущих волн. Свободные и вынужденные колебания коаксиальных оболочек с вязкой жидко96 стью в кольцевом зазоре, которая описывается линеаризованными двумерными уравнениями Навье—Стокса, исследуются аналитически в работе [12]. Продемонстрировано, что размер кольцевого зазора и вязкость жидкости оказывают значительное влияние на демпфирование колебаний. Аналогичный подход используется в [13] для текущей жидкости. Здесь анализируются сложности в реализации условия прилипания, задаваемого на стенках свободно опёртых оболочек. Показано, что влияние нестационарных вязкостных сил возрастает с уменьшением ширины кольцевого канала. В [14] для жёстко закреплённых и в [6] для консольных оболочек вводятся в рассмотрение стационарные силы вязкого сопротивления. Показано, что они оказывают существенное влияние на критические скорости течения жидкости. В рамках аналогичной модели в [15] выполнено исследование влияния на устойчивость ряда системных параметров при кольцевом течении жидкости. В [16] представлена модель, учитывающая как стационарные, так и нестационарные силы вязкого сопротивления. Они определяются из решения линеаризованных уравнений Навье—Стокса с использованием численной процедуры, основанной на конечно-разностном методе. Продемонстрировано, что эта модель лучше согласуется с экспериментальными данными [7, 17], чем модель, учитывающая только стационарные силы вязкого сопротивления. Исследование колебаний оболочек с неоднородными ограничениями в кольцевом потоке как невязкой, так и вязкой жидкости осуществлено в [18–19]. Для этих целей использован метод Релея—Ритца, в котором в качестве допустимых функций были применены формы колебаний свободно опёртых оболочек в вакууме. В работах [20–22] в ходе численного решения задачи методом конечных элементов обнаружено существенное расхождение с известными численоаналитическими решениями для тех случаев, когда потеря устойчивости осуществляется на высоких модах колебаний. В статьях [23–24] представлено исследование коаксиальных оболочек, содержащих текущую только в кольцевом канале идеальную [23] или вязкую [24] жидкость, в том числе с учётом влияния температурных эффектов. Анализ амплитудно-частотных характеристик подкреплённых конструкций, взаимодействующих с кольцевым слоем вязкой несжимаемой жидкости, осуществлён в [25]. Отметим, что в перечисленных выше работах рассматриваются оболочки, в которых оси вращения совпадают. При этом жидкая среда полностью заполняет как пространство между двумя оболочками, так и, как правило, объём внутренней оболочки. В этом случае, как и при частичном заполнении вертикально ориентированных оболочек, задача может быть рассмотрена в осесимметричной постановке. Частичное заполнение жидкостью горизонтально расположенных конструкций, как и несовпадение их осей вращения при любой ориентации, нарушает симметричность по окружной координате. Это приводит к необходимости использования более сложных пространственных моделей. В качестве примера можно привести работы [2, 3], в которых представлены трёхмерные алгоритмы, предназначенные для анализа собственных колебаний и гидроупругой устойчивости соосных горизонтально ориентированных оболочек, частично заполненных идеальной [3] или вязкой [2] жидкостью. Библиографический список исследований, в которых рассматриваются несоосные оболочки, существенно ограничен. В [26] анализируется влияние 97 соосности оболочек на динамическое поведение системы с текущей жидкостью. Показано, что смещение осей оказывает воздействие на колебательный процесс только в случае узкого кольцевого зазора между оболочками. Влияние эксцентриситета на скорости потери устойчивости не оценивалось. Двумерные исследования бесконечных оболочек с неподвижной и текущей жидкостью представлены в [27, 28]. В первой из работ отмечается, что с ростом эксцентриситета происходит возрастание присоединённой массы жидкости и, следовательно, снижение собственных частот колебаний. В следующих работах рассматриваются собственные колебания несоосных оболочек, расположенных вертикально. В аналитических исследованиях [29, 30] решение волнового уравнения, описывающего поведение жидкости в обеих полостях, осуществляется с помощью функций Бесселя. С этой целью в рассмотрение вводятся смещённые координаты, преобразование к которым выполняется посредством теоремы Бельтрами. Достоверность полученных результатов не проверяется. В работе [31] анализ пьезокерамических оболочек с кольцевым слоем жидкости осуществляется в трёхмерной постановке с использованием МКЭ. Представлено исследование низших собственных частот и форм колебаний при различных вариантах граничных условий на краях оболочек, уровне заполнения жидкостью и смещении внутренней оболочки, в том числе с учётом электроупругих свойств материала. Из приведённого библиографического обзора можно заключить, что в настоящий момент в должной степени не исследовано влияние несоосности вложенных друг в друга горизонтальных цилиндрических оболочек на критические скорости потока жидкости, текущей между ними. Выполнение такого анализа является целью настоящей работы. 1. Постановка задачи и основные соотношения. Рассматриваются горизонтально ориентированные упругие цилиндрические оболочки длиной , имеющие толщины (1) и (2) и радиусы (1) и (2) , пространство между которыми на высоту заполнено текущей со скоростью идеальной сжимаемой жидкостью (рис. 1). Здесь и далее верхние индексы «(1)» и «(2)» характеризуют внутреннюю и внешнюю оболочки соответственно. На высоту заполнения накладывается ограничение (2) (1) 6 6 2 (2) , гарантирующее смоченность обеих оболочек. Ось вращения внутренней оболочки смещена в поперечном сечении относительно наружной на величину в направлении угла , отсчитываемого от оси вращения внешней оболочки против часовой стрелки в направлении, обратном оси аппликат. Смещение подчинено условию | | < (2) (1) , обеспечивающему отсутствие контакта между двумя оболочками. Необходимо исследовать влияние отклонения оболочек от соосного расположения на границы гидроупругой устойчивости при различных значениях кольцевого зазора и уровнях его заполнения жидкостью. Определяющие поведение упругих оболочек соотношения записываются в предположении, что криволинейная поверхность конструкции может быть представлена в виде совокупности плоских сегментов [32]. Деформации в каждом из них вычисляются в рамках классической теории тонких пластин [33] в декартовых координатах ( , , ), связанных с боковой поверхностью тела: 98 Рис. 1. Сечение некоаксиальных цилиндрических оболочек, кольцевой зазор между которыми частично заполнен текущей жидкостью [Figure 1. Section of eccentric cylindrical shells with the annular gap partially filled with a flowing fluid] { ( ) ( ) ( ) } { ( ) ( ) ( ) ( ) } ( ) = , , = , , + + { 2 ( ) 2 ( ) 2 ( ) } + ( ) , , 2 . (1) 2 2 Здесь и далее прямой чертой сверху обозначены величины, записанные в координатах ( , , ), ( ) , ( ) и ( ) — компоненты перемещений срединной поверхности в направлении этих осей, = 1, 2. Физические соотношения, связывающие вектор обобщённых усилий и моментов T ( ) и вектор обобщённых деформаций ( ) , записываются следующим образом: { ( ) ( ) ( ) ( ) ( ) ( ) } ( ) . T ( ) = , , , , , = D( ) Для изотропного материала коэффициенты, входящие в матрицу жёсткостей D( ) , определяются известным образом [33]. Математическая формулировка задачи динамики упругих тел основана на вариационном принципе возможных перемещений, который в матричном виде записывается как ( ( ) ) ( ) ( ( ) ) ( ) ( ) ( ( ) ) ( ) ( ) d d + u P d = 0, (2) D u u ( ) ( ) ( ) { ( ) ( ) ( ) ( ) } где — плотность материалов оболочек; u ( ) = ( ) , ( ) , ( ) , , , — векторы перемещений и углов поворота внутренней и внешней оболочек; { } P( ) = 0, 0, ( ) , 0, 0, 0 — векторы поверхностных нагрузок оболочек; ( ) — гидродинамическое давление. Движение идеальной сжимаемой жидкости описывается потенциальной теорией, волновое уравнение которой для потенциала возмущения скорости 99 в декартовых координатах ( , , ) имеет вид [34] 2 = 1 2 2 2 2 2 + 2 , + 2 2 2 2 (3) где — скорость звука в жидкой среде. Вообще говоря, сжимаемость жидкости в случае дозвукового течения оказывает крайне незначительное влияние на критические скорости потери устойчивости, как это показано в [35]. Однако реализуемый в работе алгоритм является универсальным и может быть использован для задач аэроупругости. Предполагается, что свободная поверхность жидкости free не перемещается, на ней отсутствуют динамическое давление и поверхностное натяжение. Соответствующее граничное условие имеет вид [36] = 0. (4) Потенциал возмущения скорости на входе и выходе из кольцевого канала между оболочками подчиняется следующим граничным условиям: =0: = 0, ( ) = : / = 0. (5) ( ) На смоченных поверхностях = задаются условия непроницаемости ( ( ) ( ) ) =± + , (6) ( ) где и — поверхности, ограничивающие объёмы жидкости и оболочек ( ) . Здесь и далее знак перед формулой зависит от направления нормалей к внешним поверхностям оболочек. Гидродинамическое давление , действующее со стороны жидкости на оболочки, вычисляется из уравнения Бернулли ( ) + , ( ) = ± где — плотность жидкости. Уравнение (3) вместе с граничными условиями (4), (5) и (6) преобразуется к слабой форме с помощью метода Бубнова—Галёркина [35]: ^ d + ( 1 2 ^ 2 2 ^ 2 2 ^ ) + + d + 2 2 2 2 2 ( ^ (1) ^ (1) ) + + d (1) ( ^ (2) ^ (2) ) + d = 0, = 1, , (7) (2) где ^ и ^ ( ) — аппроксимация потенциала возмущения скорости и нормальной компоненты вектора перемещений оболочек; и — базисные функции и их количество. 100 Численная реализация. Численное решение задачи выполнено методом ^ базисные функконечных элементов [32]. Потенциал возмущения скорости , ( ) ( ) ции и мембранные перемещения оболочек ( , ) аппроксимируются линейными функциями формы, а изгибные перемещения оболочек ( ) — неконформными функциями формы Эрмита. Дискретизация областей жидкости и оболочек осуществляется пространственными призматическими и плоскими четырёхугольными конечными элементами соответственно. Связанная система уравнений, описывающая взаимодействие упругих оболочек (2) и жидкости (7), формулируется в координатах ( , , ) и записывается в матричном виде: { } { } { } + C u (1) , u (2) , + (K + A) u(1) , u(2) , = 0, M u (1) , u (2) , (8) где { } { } (2) (2) K = diag K(1) M = diag M(1) , K , K , , M , M , (1) (1) 0 0 C 0 0 A (2) (2) C= 0 0 C , A = 0 0 A , (1) (2) (1) (2) C C C A A A ( ( ) ) ( ) ( ) ( ( ) ) ( ) K ( ) B D B d , M ( ) = ( ) N N d , = ( ) ( ) ( ) F ( ( ) ) ( ) ( ) d , N Fd , A = N( ) C = ( ) ( ) ( 1 F F F F F F ) + + d , M = F F d , K = 2 2 F 2 F F C = F d , A = d , 2 2 ( ) ( ) ( ) N F d . C = F N( ) d , A = ( ) ( ) Здесь u( ) и — обобщённые векторы перемещений и углов поворота внутренней и наружной оболочек и потенциала возмущения скорости; B( ) — матрица градиентов, связывающая вектор деформаций с вектором узловых перемеще( ) ний оболочечного конечного элемента; F, N( ) и N — функции формы для потенциала возмущения скорости жидкости, обобщённого вектора узловых перемещений оболочек и его нормальной составляющей. В определяющих соотношениях (1) отсутствует уравнение для поворота вокруг оси . Если все элементы, имеющие общий узел, будут компланарны, то жёсткость в этом направлении станет нулевой. В этом случае любое возмущение, способствующее повороту, существенно повлияет на корректность решения. С целью устранения этой проблемы вводится фиктивный момент , который добавляется в ( ) ( ) матрицу жёсткости конечного элемента оболочки [32]. Матрицы K и M формируются в координатах ( , , ), связанных с боковой поверхностью конструкции. Преобразование узловых перемещений к глобальным декартовым 101 координатам ( , , ) осуществляется для каждого элемента с помощью матрицы направляющих косинусов следующим образом: u( ) = L u ( ) , где ( ) ( ) ( ) ( ) L = diag { , , , , , , , }. Тогда K = L K L, M = L M L. Ана( ) ( ) ( ) ( ) логично формируются матрицы C , A , C и A , возникающие в результате гидроупругого взаимодействия. С учётом представления возмущённого движения оболочек и жидкости ( ) ( ) exp( ) исходная система уравнений (8) в виде u(1) , u(2) , = u (1) , u (2) , сводится к обобщённой задаче на собственные значения, формулируемой как [ M 0 0 I ] { x x } [ + C K+A I 0 ] { x x } = 0, (9) — векторные функции координат, = + i — характеристигде u ( ) и v ческий показатель (i = 1), в котором полагается, что — собственная частота колебаний, а — величина, отвечающая за демпфирование системы, { } . I — единичная матрица, x = u (1) , u (2) , Вычисление комплексных собственных значений системы (9) осуществляется посредством алгоритма, основанного на неявно перезапускаемом методе Арнольди [37]. Оценка устойчивости базируется на анализе характеристических чисел , получаемых при последовательно возрастающей скорости течения жидкости. Результаты расчётов. В представленных ниже примерах рассмотрена устойчивость системы горизонтально ориентированных жёстко закреплённых на обоих краях ( = = = = = = 0) цилиндрических оболочек ( = 1 м, (2) = 0.1 м, (1) = (2) = = 5 · 10 4 м, модули Юнга (1) = (2) = = 2 · 1011 Па, коэффициенты Пуассона (1) = (2) = (1) (2) = = 0.3, = = = 7800 кг/м3 ), в кольцевом канале между которыми содержится поток идеальной сжимаемой жидкости ( = 1000 кг/м3 , = 1500 м/с). Расчёты выполнялись при различных значениях уровня заполнения кольцевого канала и радиуса оболочки (1) , определяемого внутренней (2) (1) (1) из соотношения = ( ) . Здесь параметр характеризует безразмерную величину кольцевого зазора для случая соосного расположения оболочек. Для представления полученных результатов используются безразмерные уровень заполнения жидкостью , смещение внутренней оболочки (эксцентриситет) и критические скорости потери устойчивости = (2 (2) ) 1 , = ( (2) (1) ) 1 , = [ (1 2 )/ ]0.5 · 103 . В предыдущих работах авторов [2, 3] верификация разработанной модели была осуществлена посредством сравнения собственных частот колебаний и критических скоростей течения с известными решениями, полученными для случая неподвижной [2] и текущей жидкости [3], полностью заполняющей кольцевой канал соосной конструкции. Поскольку для несоосных оболочек отсутствуют публикации с доказанной достоверностью приводимых результатов, подтверждение корректности решений в рамках описанного выше 102 конечно-элементного алгоритма осуществляется путём сравнения с данными, полученными в пакете ANSYS. В последнем при моделировании тонкостенных тел используются элементы SHELL181 (теория Рейсснера—Миндлина), а при тоделировании жидкости — FLUID30 (уравнения акустической среды) [38]. Количество используемых элементов в обеих моделях одинаково и равно 4 800 элементам для каждой из оболочек (по 40 и 120 в меридиональном и окружном направлениях) и 96 000 элементам для жидкости (по 40, 120 и 20 в меридиональном, окружном и радиальном направлениях соответственно). В таблице приведены собственные частоты колебаний (Гц) для системы с упругой внутренней и абсолютно жёсткой наружной оболочками ( = 0.3 м, (1) = 0.1 м, (2) = 1.3 м, (1) = 2 · 10 3 м, (1) = 6.9 · 1010 Па, (1) (1) = 0.3, = 2700 кг/м3 ), кольцевой зазор между которыми полностью заполнен неподвижной жидкостью ( = 1000 кг/м3 , = 1483 м/с). Здесь через и обозначено количество полуволн в окружном и меридиональном направлениях. Представленные данные, полученные при различной величине смещения внутренней оболочки , демонстрируют хорошее соответствие (погрешность не превышает 0.9%), а незначительные расхождения могут быть объяснены разными подходами к описанию оболочек. На рис. 2 представлены графики зависимостей безразмерных скоростей потери устойчивости от безразмерного смещения внутренней оболочки в направлении угла , полученные при различных значениях кольцевого зазора и уровнях его заполнения . При полном заполнении кольцевого канала жидкостью и отсутствии массовых сил угол , в направлении которого происходит смещение внутренней оболочки, не оказывает влияния на границы устойчивости, тогда как отклонение от соосного положения (увеличение по абсолютному значению) приводит к снижению критических скоростей. В случае частичного заполнения картина значительно усложняется. При смещении внутренней оболочки вдоль свободной поверхности происходит перераспределение присоединённой массы жидкости. Это приводит к падению критических скоростей независимо от направления смещения и уровня заполнения жидкостью. При = 90 положительные значения эксцентриситета соответствуют подъёму внутренней оболочки из жидкости, а отрицательные — её погружению. В результате изменения площади смачиваемой Сравнение собственных частот колебаний (Гц) при различном эксцентриситете [Comparison of natural vibration frequencies (Hz) at different values of eccentricity ] = 0.0 = 0.2 = 0.4 = 0.6 = 0.8 ANSYS Calc. val. ANSYS Calc. val. ANSYS Calc. val. ANSYS Calc. val. ANSYS Calc. val. 3 4 1 2 5 6 4 5 3 1 1 1 1 1 1 2 2 2 469.88 473.91 532.20 534.82 616.48 890.39 924.65 944.73 1023.8 469.28 472.80 531.99 534.37 613.29 882.54 921.50 940.29 1021.1 456.45 478.70 514.48 551.66 614.50 887.87 910.39 947.01 1016.4 455.70 477.75 514.17 551.31 611.32 880.05 907.08 942.77 1013.7 429.05 473.90 500.67 567.11 608.07 879.38 870.32 942.45 978.25 428.37 473.11 500.40 566.74 604.94 871.63 867.15 938.43 976.24 383.75 453.44 495.87 579.89 595.45 861.32 788.88 881.21 936.57 383.17 452.71 495.25 579.51 592.42 853.73 786.20 879.10 933.12 310.61 417.70 490.81 590.27 571.89 823.22 631.70 766.01 893.17 310.12 416.96 490.11 589.87 569.07 815.95 629.72 763.95 889.79 103 Рис. 2. Зависимость безразмерных критических скоростей от смещения внутренней оболочки в направлениях = 0 (слева) и = 90 (справа) при различных уровнях заполнения кольцевого канала для = 1/2, 1/10 и 1/100 [Figure 2. Dependence of dimensionless critical velocities on the eccentricity in the directions = 0 (left) and = 90 (right) at different levels of filling of the annular channel for = 1/2, 1/10, and 1/100] 104 поверхности и перераспределения гидродинамического давления происходит, соответственно, повышение или снижение критических скоростей течения. Данные, приведённые на рис. 2, также позволяют сделать вывод о том, что размер кольцевого зазора не оказывает качественного влияния на зависимости критических скоростей от смещения внутренней оболочки. Уменьшение значения приводит к снижению границ устойчивости, что установлено ранее для соосных цилиндрических оболочек при полном или частичном заполнении кольцевого канала [1, 3]. На рис. 3 слева представлены зависимости критических скоростей потери устойчивости от смещения внутренней оболочки при различных уровнях заполнения и направлении смещения . Из представленных результатов видно, что c повышением уровня заполнения происходит снижение влияния угла смещения на границу устойчивости до тех пор, пока при полном заполнении он не перестанет оказывать какое-либо влияние. Минимальное значение скорости всегда достигается при одной и той же величине эксцентриситета ( 1) вне зависимости от угла и уровня заполнения . Максимальное же значение может достигаться при различных значениях смещения на полуинтервале [0; 1) в зависимости от выбора и , что объясняется изменением площади смачиваемой поверхности внутренней оболочки и величины присоединённой массы жидкости, как было отмечено ранее. Высказанные утверждения справедливы при положительных значениях угла , а отрицательные величины смещения соответствуют положительным в направлении отрицательных углов . Данные, представленные на рис. 3 слева, позволяют также говорить о возможности существования таких конфигураций, в которых при надлежащем подборе всех параметров обеспечивается более высокий порог гидроупругой устойчивости по сравнению с соосным расположением оболочек. Зависимости критической скорости от угла при различном уровне заполнения и смещении внутренней оболочки представлены на рис. 3 справа. При уменьшении угла имеет место такая область его значений, где критическая скорость практически остаётся неизменной. Размер данной области возрастает с ростом уровня жидкости до тех пор, пока при полном заполнении угол смещения не перестанет оказывать влияние (прямая параллельна оси абсцисс) на границы устойчивости. При этом с ростом эксцентриситета критическая скорость меняется в больших пределах, предоставляя более широкие возможности управления динамическим поведением. На рис. 4, 5 приведены собственные формы колебаний оболочек для = 1/10 при различных вариантах заполнения кольцевого канала, величины и направления смещения внутренней оболочки. При построении мод скорость течения жидкости для каждой конфигурации задавалась близкой к критической скорости потери устойчивости cr . На изображённых сечениях пунктиром показаны оболочки в недеформированном состоянии, а сплошными линиями — в деформированном; уровень заполнения показан серым цветом. В пространственных формах колебаний перемещения отмасштабированы для наглядности представления результатов. Реальные значения, полученные из решения спектральной задачи (9), приведены на цветовой шкале, которая является общей для обеих оболочек. Здесь красным цветом обозначено перемещение в направлении внешней нормали к поверхности оболочки, а синим — 105 Рис. 3. Зависимость безразмерных критических скоростей для = 1/10 при = 0.25, 0.5, 0.75 от смещения внутренней оболочки при различных значениях угла (слева) и угла при различных значениях смещения (справа) [Figure 3. Dependence of dimensionless critical velocities on the eccentricity at different values of the angle (left) and the angle at different values of the eccentricity (right) for = 1/10 at = 0.25, 0.5, and 0.75] 106 = 1.0, = 0.0 = 1.0, = 0.5 = 0.5, = 0.0 = 0.5, = 0.5, = 0 = 0.5, = 0.5, = 90 = 0.5, = 0.5, = 270 = 0.25, = 0.0 = 0.25, = 0.5, = 0 = 0.25, = 0.5, = 90 = 0.25, = 0.5, = 270 Рис. 4. Собственные формы колебаний взаимодействующих с жидкостью несоосных оболочек в поперечном ( = /2) и продольном ( = 0) сечениях для = 1/10 при = 0.25, 0.5 и 1.0, cr и различных значениях смещения и угла [Figure 4. The transverse ( = /2) and longitudinal ( = 0) cross-sections of mode shapes of eccentric cylindrical shells interacting with the fluid for = 1/10 at = 0.25, 0.5, and 1.0, cr and different values of the eccentricity and the angle ] 107 Рис. 5. Собственные формы колебаний взаимодействующих с жидкостью несоосных оболочек (поперечное сечение = /2, внутренняя и наружная оболочки) при смещении внутренней оболочки на = 0.95 в направлениях = 0 (сверху), = 90 (по центру) и = 270 (снизу); = 1/10, = 0.25, cr [Figure 5. Mode shapes of eccentric shells interacting with the fluid (the cross section at = /2, the inner and outer shells) in the case of the eccentricity = 0.95 in the directions = 0 (top), = 90 (center), and = 270 (bottom); = 1/10, = 0.25, cr ] 108 в противоположном ей. В случае полностью заполненных ( = 1.0) соосных оболочек ( = 0.0) окружные полуволны имеют одинаковую высоту, размер которой различен для внутренней и наружной оболочек. При наличии эксцентриситета ( = 0.0) высота полуволн в пределах одной оболочки будет разной. Максимальные перемещения возникают на участках боковых поверхностей, соответствующих минимальному расстоянию между оболочками. При изменении уровня жидкости ( < 1.0) также происходит образование полуволн разной высоты, а максимальные перемещения формируются на тех частях боковых поверхностей оболочек, которые взаимодействуют с жидкостью. В случае частичного заполнения несоосных конструкций наибольший размер полуволн наблюдается на смоченных поверхностях. На основании представленных данных можно заключить, что уровень заполнения кольцевого зазора жидкостью оказывает большее влияние на перемещения оболочек в процессе колебаний, чем эксцентриситет. Заключение. Представлено численное исследование гидроупругой устойчивости тонкостенных горизонтально ориентированных несоосных цилиндрических оболочек, взаимодействующих с потоком жидкости, текущим в кольцевом канале между ними. Моделирование осуществлено в трёхмерной постановке с использованием метода конечных элементов. Проанализировано влияние отклонения оболочек от соосного положения на критические скорости потери устойчивости при различном уровне заполнения кольцевого канала. Установлено, что, как и для случая соосных оболочек, снижение уровня заполнения приводит к повышению границы гидроупругой устойчивости при прочих равных параметрах системы. Продемонстрировано, что с понижением уровня заполнения кольцевого зазора жидкостью возрастает влияние величины отклонения внутренней оболочки от соосного положения, стабилизирующий или дестабилизирующий характер которого определяется направлением смещения. За счёт подбора комбинации «смещение–угол» для определённых конфигураций можно добиться повышения границ устойчивости по сравнению с соосным положением оболочек. При этом наибольшее изменение критической скорости имеет место при смещении внутренней оболочки в направлении ортогональном свободной поверхности жидкости.

About the authors

Sergey Arkadievich Bochkarev

Institute of Continuous Media Mechanics UB RAS

Email: bochkarev@icmm.ru
Candidate of physico-mathematical sciences, no status

Sergey Vladimirovich Lekomtsev

Institute of Continuous Media Mechanics UB RAS

Email: lekomtsev@icmm.ru
Candidate of physico-mathematical sciences, no status

Aleksandr N. Senin

Institute of Continuous Media Mechanics UB RAS

References

Supplementary files