Вычислительные технологии, 2022, том 27, № 4, с. 15-32. © ФИЦ ИВТ, 2022 Computational Technologies, 2022, vol. 27, no. 4, pp. 15-32. © FRC ICT, 2022
ISSN 1560-7534 elSSN 2313-691X
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ
DOI: 10.25743/ICT.2022.27.4.003
Численное моделирование собственных колебаний покоящейся на упругом основании цилиндрической оболочки, частично заполненной жидкостью
С. А. Бочкарёв
Институт механики сплошных сред УрО РАН, 614068, Пермь, Россия Контактный автор: Бочкарёв Сергей Аркадьевич, e-mail: bochkarev@icmm.ru Поступила 29 декабря 2021 г., доработана 10 июня 2022 г., принята в печать 25 июня 2022 г.
Представлены результаты исследований собственных колебаний круговых цилиндрических оболочек, полностью или частично заполненных неподвижной жидкостью и лежащих на упругом двухпараметрическом основании Пастернака. Поведение упругой конструкции и сжимаемой жидкости описывается в рамках классической теории оболочек и уравнений Эйлера. Уравнения движения оболочки, как и акустическое волновое уравнение, сводятся к системам обыкновенных дифференциальных уравнений. Решение сформулированной краевой задачи осуществляется методом ортогональной прогонки Годунова. Проанализированы зависимости минимальных частот колебаний от характеристик упругой среды с разными вариантами неоднородности по длине тела при различных уровнях заполнения оболочки жидкостью. Продемонстрировано, что с повышением уровня заполнения оболочек жидкостью влияние упругого основания на частотный спектр конструкции снижается.
Ключевые слова: классическая теория оболочек, сжимаемая жидкость, метод ортогональной прогонки Годунова, метод прямых, собственные колебания, упругая среда Пастернака.
Цитирование: Бочкарёв С.А. Численное моделирование собственных колебаний покоящейся на упругом основании цилиндрической оболочки, частично заполненной жидкостью. Вычислительные технологии. 2022; 27(1):1о 32. DOI: 10.25743/ICT.2022.27.4.003.
Введение
Тонкостенные цилиндрические оболочки, в том числе содержащие жидкость, как неотъемлемая часть многообразных индустриальных приложений находят применение в различных областях техники. В процессе эксплуатации они не только подвергаются разнообразным статическим и динамическим нагрузкам, но и контактируют с внешней окружающей средой. При моделировании этого взаимодействия широкое распространение получили приближенные модели (Винклера, Пастернака, Керра, упругого полупространства и т.п.), согласно которым действие окружающей среды заменяется упругим основанием [1-4]. При этом внешняя упругая среда оказывает существенное влияние на динамические характеристики взаимодействующей с ней конструкции. Этот факт послужил причиной появления значительного количества публикаций, где объектами исследований выступали тонкостенные тела, выполненные из разнообразных материалов,
находящиеся под действием различных температурных, механических, центробежных, аэро- и гидродинамических нагрузок, проанализированных с применением как аналитических, так и численных методов. Обзор публикаций, в которых осуществлен линейный анализ влияния упругого основания на спектр колебаний или границы устойчивости пустых круговых цилиндрических оболочек или прямолинейных труб, выполнен в [5].
В меньшей степени изучено воздействие упругого основания на динамические характеристики цилиндрических оболочек с жидкостью [6-16], Изотропные и функционально-градиентные (ФГ) оболочки, свободно опертые на обоих краях, полностью заполненные неподвижной жидкостью и целиком погруженные в двухпараметричеекую упругую среду, рассмотрены в [6, 7], Решение динамических уравнений с учетом гидродинамической нагрузки, определенной из акустического волнового уравнения, основано на методе распространяющихся волн. Вертикальные ФГ-оболочки, частично заполненные несжимаемой потенциальной жидкостью, анализируются в [8, 9], Здесь оболочки с различными комбинациями граничных условий также полностью размещаются на двухпараметричееком упругом основании, В первой работе с помощью метода Ритца, в котором полиномы Чебышёва используются для определения модальных функций в продольном направлении, оценено влияние коэффициентов основания на собственные частоты колебаний при различном уровне жидкости, В статье [9] численное решение строится на основе полуаналитического варианта метода конечных полос. Отмечается, что упругое основание Винклера оказывает большее воздействие на собственные частоты ФГ цилиндрической оболочки, частично заполненной жидкостью, чем основание Пастернака,
Влияние упругой среды на собственные частоты колебаний и границы гидроупругой устойчивости в случае текущей жидкости оценено в работах [10-12], Колебательные характеристики ФГ цилиндрической оболочки с жидкостью, погруженной в упругую среду (модель Винклера) и подвергающейся воздействию механических и температурных нагрузок, исследуются с помощью метода модальных расширений в [10], Продемонстрировано различное влияние упругого основания на частоты колебаний для пустых и заполненных текущей жидкостью оболочек, В [11] с применением метода обобщенных дифференциальных квадратур выявлено, что для углеродных нанотрубок с текущей внутри вязкой жидкостью, описываемой модифицированными соотношениями Навье — Стокса, повышение коэффициента постели в модели Винклера приводит к росту критических скоростей течения жидкости. Подземный нефтепровод моделируется в работе [12] в виде цилиндрической оболочки большого диаметра. Полученные аналитические соотношения для собственных частот конструкции учитывают влияние таких параметров, как скорость течения жидкости, глубина залегания и отклик однопарамет-рической упругой среды.
Оболочки с жидкостью, частично покоящиеся на упругом основании, изучены в статьях [13-16], В этих работах анализируются симметричные и антисимметричные моды колебаний изотропных и ФГ цилиндрических оболочек, покоящихся как на одно- [13, 14], так и двухпараметричееком [15, 16] упругом основании, неоднородном как в окружном, так и меридиональном направлениях. При этом, как правило, применяются подходы, предложенные этими же авторами для анализа пустых оболочек, В [14] используется разработанный в [17] гибридный конечный элемент, в котором точные функции перемещений определяются непосредственно из уравнений теорий оболочек Сандерса, В перечисленных публикациях влияние упругого основания, являющегося однородным в окружном направлении и неоднородным в меридиональном, на собственные частоты
колебаний оболочек вращения, полностью или частично заполненных жидкостью, не исследовалось. Выполнение такого анализа, является одной из целей настоящей работы.
Для решения сформулированной краевой задачи используется метод ортогональной прогонки Годунова [18]. В монографии [19] среди достоинств этого метода, как и других алгоритмов, сводящих краевые задачи к задачам Коши, выделяются такие, как численная устойчивость и высокая точность вычислений. Для задач гидроупругости этот метод применяется редко и традиционно [20, 21] используется в различных реализациях. В первой из них в уравнения движения включается аналитическое выражение для гидродинамического давления, полученное из решения волнового уравнения и записанное в членах функций Бесселя [22, 23]. В комбинированном методе волновое уравнение тем или иным способом сводится к системе обыкновенных дифференциальных уравнений, которая решается совместно с системой уравнений для упругого тела. В частности, в [19, 24] для анализа вынужденных колебаний коаксиальных оболочек, содержащих жидкость в кольцевом канале, предложено использовать процедуры сведения к замкнутой системе обыкновенных дифференциальных уравнений, основанные на методе прямых и кубической сплайн-аппроксимации. Насколько известно автору, комбинированный подход в рамках ортогональной прогонки Годунова в случае одиночных оболочек, полностью или частично заполненных жидкостью, ранее не использовался.
1. Постановка задачи
Рассматривается упругая вертикально расположенная цилиндрическая оболочка (рис. 1) длиной радиусом Я и толщиной ее внутренняя область У}- полностью или частично на высоту Н заполнена идеальной сжимаемой жидкостью. Оболочка целиком погружена в двухпараметрическое упругое основание с коэффициентами постели и Кр, отвечающими за растяжение — сжатие (модель Винклера) и сдвиг (модель Пастернака) соответственно. В меридиональном направлении упругое основание является как однородным, так и неоднородным. В последнем случае оно представляет собой чередование участков с его наличием или отсутствием в виде трех различных комбинаций (рис. 2). В первом варианте оболочка разбивается на два участка, один из которых (¿1) содержит упругое основание, а его размер изменяется от нуля до Ь. Без-
Рис. 1. Расчетная схема цилиндрической оболочки с жидкостью, расположенной на упругом основании
Fig. 1. Computational scheme for a cylindrical shell with fluid embedded in an elastic foundation
W-
W W
II
III
( I
h
1 '
1 I
h
1 1
-
Рис. 2. Различные варианты неоднороднсих) (но длине оболочки с жидкостью) упругого основания
Fig. 2. Different variants of elastic foundation heterogeneous along the length of the shell filled with
I
размерный коэффициент длины характеризующий область контакта упругой среды и поверхности оболочки, определяется в этом случае как £ = 1\/Ь. Во втором варианте оболочка разбивается па три участка, из которых два крайних содержат упругое основание и их размеры изменяются от нуля до Ь/2 (£ = [11 + 13}/Ь). Наконец, в третьем варианте оболочка также разбивается па три участка, из которых только центральный содержит упругое основание, а его размер равномерно увеличивается к обоим краям оболочки от ее середины (£ = 12/Ь). Необходимо проанализировать влияние упругого основания па собственные частоты колебаний оболочки с жидкостью при различных уровнях ее заполнения и разных вариантах граничных условий, задаваемых па краях конструкции.
2. Основные соотношения и метод решения
Дня классической теории оболочек, основанной на гипотезах Кирхгофа — Лява, компоненты вектора деформации Е^ в криволинейной системе ко ординат (х,9, г) могут быть записаны в виде |25|
Е11 = £ц + ZK.11, Е22 = £22 + ?К22, Е12 = £12 + 2^12,
где
^11
ди
дх'
«22
I
R дО
22
«12
1 idv
R\de +W
12
dv 1 ди
дХ + Rdd1
«11
d(h дх
1
R
idv _ dfhN V&X - ~дв)
---^ , 01 , 02 = - [v- —
dw дх
1
R
( dw\
(i)
Здесь и, V и т — меридиональная, окружная и нормальная составляющие вектора перемещений оболочки; — углы поворота недеформируемой нормали.
Физические соотношения, устанавливающие связь между вектором усилий и моментов Т = {Т11,Т22, Б, М11, М22, Н}т и вектором обобщенных деформаций
в матричном виде записываются как
£ = {£ 11, £221 £12) «11) «22, 2«12 }T,
Т = Бе
А В В С
Б
(2)
(а*,) = ! (1^г2)^ М = l,3
н
Я11 = Я22 = Е/[1 - А, Я12 = Я21 = »Я Я33 = Е/2/[1 + и].
Здесь Е, V — модуль упругости и коэффициент Пуассона изотропного материала оболочки.
Уравнения движения оболочки с учетом реакции упругой среды
дТп 1 дБ д2и дБ 1 дТ22 1 {п дН \ д2у _
^ + пдё - р0~дё = 0, д^ + + Я ^22 + -р0№ = 0,
дЯи + 1 дЯ22 1Т д2т К +К (д2т + 1 д2т\ + (3)
+ К^дГ - ЯТ22 - р0 ~д2+ Кр\-дх2 + Я Ж2) + * = ° (3) дМп 1дН дН 1 дМ22
+ -Я11 = 0, + К^Т-Я22 = 0,
где Ян ~ поперечные сплы; р0 = J рс1г, р — плотность материала; р — гидродинами-
н
ческое давление, которое в области V/ в цилиндрической системе координат (х,9, г) описывается акустическим волновым уравнением
у2 =1 +1. д>21 + = 1 д2^ ^
г дг \ дг) г2 дв2 дх2 с2 д12
Здесь с — скорость звука в жидкости. На смоченной (г = Я) и свободной (х = Н) поверхностях, оси вращения оболочки (г = 0) и нижнем крае (х = 0 давление р удовлетворяет условиям
д2 т
VP• Ч=д = -Р/ —, (5)
V • п|г=0 = 0, (6)
х = 0: ^ = 0, х = Н : р = 0, (7)
д х
где п — единичные внешние нормали к области жидкости У/', р/ — плотность жидкости,
х = Н
8/гее не перемещается, на ней отсутствуют динамическое давление и поверхностное натяжение.
Раскладывая все компоненты уравнений (1), (2) и (4) в ряды Фурье по окружной X(х, в) = ^Х,(х) совт, ¥(х, в) = ^Уз(х)зтав),
3=0 3=0
X = {и,т, 01, Ец, Е22, К11, К22, Тц, Т22, Мц, М22, Я11,р}, ¥ = {у, 62, Еи, Ки, Б, Н, Я22},
геометрические (1) и физические (2) соотношения, а также уравнения движения (3) сведем к системе восьми обыкновенных дифференциальных уравнений первого порядка
У1 = Тц, У2 = Б+
2 Н
К
Уз = М11, У4 = Яп + ЗкН, у5 = и, у6 = V, Ут = т, у8 = 01.
Здесь 3 — номер гармоники при разложении в ряд Фурье, зп = 3/К. С учетом этого и принимая во внимание у(£) = у ехр(иЛ;), искомая система может быть записана следующим образом:
^ = f 0>, у), (8)
где
дх
. (2 Н \ 2 г • гт Я22 2
11 = ^ I - У2) — Ш РоУб, 12 = ЗКТ 22 — PoУ6,
Т2
¡3 = УА — 2jRН,
к = К — ЗкЯ22 — "2роУт + Кшут + Кр (кц + ]2кут) — P,
(9)
¡5 = е^ ¡6 = £12 + 3 К У 5, ¡7 = - У8, ¡8
Кц.
Здесь ш — характеристический показатель; I2 = —1, Входящие в выражения (9) величины вычисляются по формулам
Н = Ьзз£12 + 2С33К12, я22 = —ЗкМ22,
к11
0-11 (Уз — Ьи£22 — С12К22) — Ьц {у1 — С112£22 — ^2^2)
а11Сц — Щ_1
£11
У1 — С12£22 — ЬцКц — ЬПК22
ап
£12
Т22 = С12 11 + С22 22 + 12 к11 + 22 к22,
У2 — 2 {ь-зз + 21) >я (| + ■»)
I. I
Ьзз + —
азз + 4-
Сзз
К
где
К
т 6
£22 = 3яУ6 + К К22 = ЗвЯ2, 02 = ^ — ^Ут,
к12
£12 + Зк У 5
К
— Зк У 8, М22 = Ь12611 + Ь22£22 + С12К11 + С22К22.
Для сведения уравнения (4) к системе обыкновенных дифференциальных уравнений используем метод прямых (дифференциально-разностный метод) [26], Для этого разделим область жидкости по радиусу па п равных частей и через внутренние точки деления проведем семейство параллелей г = Гг (г = 1, п — 1), На каждой такой прямой дифференциальное уравнение (4) заменим приближенным обыкновенным дифференциальным уравнением для функции р(х, Гг), предполагая, что она дважды непрерывно дифференцируема на рассматриваемом диапазоне. Для этого в уравнении (4) избавим-
дифференцирования
др дг
д 2р дг2
~ 7Г [P(х, П+1) — p(х, П-1)] , =п 2т
— [р(х, п+1) — 2р(х, п) + р(х, П-1)] 2
= К/ п
г=г
Для "крайних" прямых [27]: др дг
др дг
Вводя сокращенные обозначения
~ [4р(х, г 1) - р(х, Г2) - 3р(х, го)] , ~ — [3р(х, гп) - 4р(х, гп-г) + р(х, гп-2)}.
тз 2Т
г=К
, N / N дРп) ,, . д'2p(х, п) ,, p{х, Гг)= Pí{х), -—- = Pi(х), -дх2- = Р* (X),
подставляя формулы (10) в исходное уравнение (4), придем к системе (п - 1) уравнений в обыкновенных производных
р'! (х) = р\Рг(х) + /Зг2рг+1(х) + ргзРг-1(х), (12)
т 2 , Г пг т + 2гг 2гг -т где р\ = — + —2--о,Ро =--Рз =-^т3- система сводится к нормальному
Т2 Г2 С2 2г'г Т2 2г'г Т2
виду Коши следующим образом:
У9 = Р1(х), у'9 = У10, у'10 = р1У9 + р1У11 + р\ро, У11 = Р2 (х) , у'11 = У12, у'12 = Р1У11 + Р2У13 + РзУ9, .............................^..................................................(13)
У7+2(п-1) = Pn-1(х), У 7+2(п-1) = У 8+2(n-1),
У8+2(п-1) = Р1-1У7+2(п-1) + Р2-1Рп + Р3-1У5+2(п-1). Значения р0ъ рп в (13) определяются исходя из условий (5), (6) и уравнений для "крайних" прямых (11):
Ро = 1(4У9 - У11), Рп = 3 (2тр}ш2у7 + 4У7+2(п-1) - У5+2(п-1)) .
Примем во внимание, что в соотношениях (9) р = рп. Таким образом, совокупное число неизвестных, определяемых системами уравнений (8) и (13), составляет т = 8 + 2(п- 1) для заданной окружной гармоники.
На краях оболочки задаются однородные граничные условия
Уг1х=о 5г + Уг+41х=о (1 - й) = 0, 1=1,..., 4, (14)
Уг1х=Ь5г+4 + Уг+Ах=Ь (1 - $г+4)=0, i=1,..., 4, (15)
т
Уг1х=о = 0, г = 10, 12,..., т, (16)
2
т
Уг1х=н = 0, г = 9,11,...,~т - 1, (17)
где 5г = 0, если заданы кинематические, и ^ = 1, если заданы статические граничные условия.
Объединенная система уравнений (8), (13) с граничными условиями (14)—(17) решается методом ортогональной прогонки Годунова с численным интегрированием дифференциальных уравнений методом Рунге — Кутты четвертого порядка точности и орто-гонализацией векторов по методу Грамма —Шмидта, Для этих целей общее решение системы (8), (13) представляется в виде
т/2
У = У к,
к=1
где Ck — некоторые константы и yk — совокупность линейно независимых решений систем (8), (13), удовлетворяющих граничным условиям (14), (16), В результате интегрирования по заданному интервалу и удовлетворения граничных условий (15), (17) для определения постоянных Ck получаем следующую алгебраическую систему:
т/2
^Ckfik = 0, i=1. (18)
k=i
Решение задачи сводится к вычислению таких значений ш, при которых существует нетривиальное решение системы (18), Необходимым условием для этого является равенство нулю определителя матрицы | k(ш)| = 0, Для нахождения ш используется комбинация шагового метода, посредством которого вычисляются такие значения, при которых происходит смена знака определителя |f\k(ш)и метода, уточняющего ш в полученном диапазоне, В качестве последнего используется метод деления отрезка пополам. Численный алгоритм реализован в виде программного кода на языке Fortran,
3. Результаты численного моделирования
В численных примерах рассматривается цилиндрическая оболочка (L/R = 3, R/h = 301.5 v = 0.3), заполненная жидкостью (pf = 1000 кг/м3, с = 300 м/с) и имеющая следующие граничные условия: свободное опирание (v = w = Т11 = М11 = 0, SS) или жесткое закрепление (и = v = w = 91 = 0, СС) на обоих краях; консольное закрепление (Тп = 0, S + 2H/R = 0 Мп = 0,Q ii + jrH = 0, CF), Численными экспериментами установлено [28], что для указанной конфигурации 25 прямых обеспечивают вычисление собственных частот с подходящей точностью и приемлемой вычислительной эффективностью, Для отображения полученных результатов используются безразмерные параметры частоты Q, коэффициентов постели кр [29] и уровня жидкости q:
О / о^ ш2 Т^ Т^К2 Л М1 - П ЕЬ3 Н
В представленных ниже расчетах используется модель упругого основания Пастернака, при которой к = —— = —— и Ктах = Ктах = 3600,
' Г г ^тах ^тах ш р
к-ш кр
Для верификации разработанного алгоритма осуществлено сравнение собственных частот колебаний заполненных жидкостью оболочек с учетом и без учета (£ = 0) однородного упругого основания. Сопоставление низших собственных частот колебаний ш, определенных в рамках описанного алгоритма, с известными решениями для свободно опертой (Ь = 30 м, R = 25 м, К = 0.03 м, Е = 206 ГПа, и = 0.3, р = 7850 кг/м3) и консольно закрепленной (Ь = 0.231 м, R = 0.07725 м, К = 0.001^ м, Е = 205 ГПа, V = 0.3, р = 7800 кг/м3) цилиндрических оболочек, частично заполненных несжимаемой жидкостью, приведено в табл. 1 и 2 (£ = 0). Здесь М — количество полуволн в меридиональном направлении. На основании данных, представленных в этих таблицах, можно заключить, что значения частот, полученных с использованием различных методов, разных теорий оболочек и моделей, описывающих поведение жидкости, достаточно близки друг к другу. Отметим также, что второй столбец с данными из работы [30] в табл. 1 отражает результаты, полученные с учетом гравитационных эффектов
Таблица 1. Собственные частоты колебаний ш (рад/с) свободно опертой (SS) цилиндрической оболочки, частично заполненной жидкостью: j = 4, С = 0 £ = 0.72 Table 1. Natural frequencies of vibrations ш (rad/s) of a simply supported (SS) cylindrical shell partially filled with fluid: j = 4, £ = 0 = 0.72
м Amabili et al. [30] Amabili et al. [30] Shahbaztabar et al. [8] Kim et al. [31] Расчет
1 13.658 13.64 14.04 14.054 14.022
2 34.441 34.434 34.65 34.672 34.563
3 49.692 49.688 49.57 49.629 49.330
4 61.877 61.874 61.63 61.556 61.115
5 71.804 71.802 71.95 71.476 70.794
Таблица 2. Собственные частоты колебаний ш (Гц) консольно закрепленной (CF) цилиндрической оболочки, частично заполненной жидкостью: £ = 0 ^ = 0.697
Table 2. Natural frequencies of vibrations ш (Hz) of a cantilevered (CF) cylindrical shell partially filled with fluid: £ = 0 <i = 0.697
M 3 Maxuch et al., эксп. [32] Maxuch et al. [32] Бочкарёв и др. [33] Ergin et al. [34] Кашфутдинов и др. [35] Расчет
1 1 - 1407.4 - 1407.3 1419.6 1399.5
2 582.0 672.7 664.5 669.8 669.5 668.8
3 522.0 543.1 538.0 542.0 541.4 539.8
4 798.0 806.0 799.1 806.8 803.1 799.3
5 1196.0 1188.4 1178.8 1195.5 1181.9 1172.6
6 1679.7 - 1657.3 1693.6 1665.4 1648.5
2 3 1394.0 1425.3 1412.2 1434.0 1419.6 1417.7
4 1244.0 1253.2 1239.2 1261.4 1250.3 1246.5
5 1546.0 1553.8 1535.7 1557.7 1549.4 1544.4
6 - - 2092.0 2119.7 - 2092.7
на свободной поверхности жидкости (плескания). Минимальные расхождения с решением в более простой постановке подтверждает известный факт [22] о малости влияния гравитационных волн на свободной поверхности при осесимметричиых колебаниях,
В табл. 3 представлены результаты работы [7], где с помощью метода распространяющихся волн выполнен анализ цилиндрических оболочек, имеющих аналогичные параметры, за исключением длины (Ь = 0.41 м), но полностью погруженных в упругое основание, однородное как в окружном, так и меридиональном направлении (£ = 1), Здесь для шести окружных гармоник приведены низшие частоты ш (Гц), вычисленные как для пустой оболочки, так и оболочки с жидкостью (^ =1) при различных комбинациях коэффициентов постели (Кш, Кр). Из представленных данных следует, что получаемые в рамках разработанного алгоритма результаты идеально согласуются в случае пустых оболочек и имеют незначительные отличия для оболочек с жидкостью, которые снижаются с повышением номера окружной гармоники ].
На рис, 3, о приведены зависимости безразмерных минимальных частот колебаний О от безразмерного уровня заполнения оболочки жидкостью я, полученные для оболочек с различными граничными условиями при отсутствии упругого основания. Представленные результаты демонстрируют известную закономерность — с увеличением уровня заполнения вследствие возрастания присоединенной массы жидкости собственные
Таблица 3. Собственные частоты колебаний ш (Гц) свободно опертой (SS) пустой и полностью заполненной жидкостью цилиндрической оболочки, полностью лежащей на упругом основании с разными комбинациями коэффициентов постели (Kw, Кр): £ = 1, q = 1 Table 3. Natural frequencies и (Hz) of a simply supported (SS) empty or fluid-filled cylindrical shell completely resting on the elastic foundation at different combinations of foundation moduli (Kw,
кру. e = 1 ? = 1
(1.5 ■ 107, 0) (0,1.5 ■ 107)
3 Shah et al. [7] Расчет Shah et al. [7] Расчет
пустая с жидкостью пустая с жидкостью пустая с жидкостью пустая с жидкостью
8 354.16 145.74 354.16 149.25 6026.2 2456.2 6026.2 2537.2
9 360.25 155.86 360.25 158.63 6737.2 2887.5 6737.2 2962.3
10 384.24 173.65 384.24 175.90 7452.5 3336.9 7452.6 3404.7
11 421.95 198.17 421.95 199.96 8170.9 3802.6 8170.9 3862.4
12 470.61 228.72 470.61 230.05 8891.7 4282.9 8891.8 4333.5
13 528.41 264.85 528.41 265.66 9614.6 4776.6 9614.6 4816.7
n
6.5
3.5
0.5
ss \
\
4^=6 сс
4 cf ^5
0.25
0.5
0.75
Рис. 3. Зависимости частот колебаний О от уровня жидкости q, полученные для оболочек с различными вариантами граничных условий: a — без учета упругого основания; б — с учетом однородного упругого основания, к = 0.278, £ = 1
Fig. 3. Dependences of vibration frequencies О versus the fluid level q obtained for shells with various boundary conditions: a — without considering the elastic foundation; б — with consideration for homogeneous elastic foundation, к = 0.278, £ = 1
a
частоты колебаний конструкции снижаются. При этом происходит изменение номера обозначены места смены формы колебаний с минимальной частотой, а цифрами указаНа рис. 3, б представлены аналогичные зависимости, но полученные в случае погружения оболочки в однородное упругое основание с безразмерными коэффициентами постели к = 0.278. Наличие упругого основания приводит к значительному росту частот. При этом минимальная частота колебаний смещается в область низших гармоник, а поведение оболочки с жидкостью становится подобным балке. Несмотря на это, снижение частот с ростом уровня жидкости становится более существенным, а следовательно, конструкция с жидкостью, расположенная на упругом основании, в большей степени подвержена изменению частотного спектра.
СС
CF
SS
Рис. 4. Изоповерхности минимальной частоты колебаний О в зависимости от однородного упругого основания к и уровня жидкости ^ для оболочек с разными граничными условиями Fig. 4. Isosurfaces of minimum vibration frequency О as a function of homogeneous elastic foundation к and liquid level с for shells with different boundary conditions
a
Рис. 5. Изоповерхности минимальной частоты колебаний О в зависимости от длины упругого основания £ и уровня жидкости с для свободно опертой (SS) оболочки с разными вариантами неоднородности: a — к = 0.278; б — к = 0.556
Fig. 5. Isosurfaces of minimum frequency О as a function of elastic foundation length £ and fluid level с for simply supported shells (SS) with different variants of inhomogeneity: a — к = 0.278; б -к = 0.556
О
основания £ и уровня жидкости ^ для консольно закрепленной (CF) оболочки с разными вариантами неоднородности: a — к = 0.278; б — к = 0.556
Fig. 6. Isosurfaces of minimum vibration frequency О as a function of elastic foundation length £ and fluid level с for cantilevered (CF) shell with different variants of heterogeneity: a — к = 0.278; б -к = 0.556
a
Наглядно обрисовать зависимость собственных частот колебаний О от величины однородного упругого основания к и высоты жидкости ^ в оболочках с разными комбинациями граничных условий дают возможность данные, показанные на рис. 4. Приведенные здесь изоповерхности предоставляют более полную информацию о распределении частотного спектра и в зависимости от требований минимизации или максимизации частоты позволяют подобрать подходящую совокупность параметров. Исходя из качественной близости результатов, получаемых для свободно опертых и жестко закрепленных оболочек, дальнейшие вычисления будут выполнены только для граничных условий, соответствующих свободному опиранию и консольному закреплению.
Оценить влияние неоднородного в меридиональном направлении упругого основания на динамические характеристики оболочек, частично заполненных жидкостью, позволяют результаты, представленные на рис. 5 и б. Здесь приведены изоповерхности безразмерной минимальной частоты колебаний О, полученные при различной длине упругого основания £ и разных уровнях жидкости ^ в свободно опертой (рис. -5) и консольно закрепленной (рис. б) оболочках. Расчеты осуществлены для трех вариантов неоднородности упругой среды (см. рис. 2) и коэффициентов постели, равных к = 0.278 (рис. -5, а и 6, а) и к = 0.556 (рис. 5, б и 6,6). Из рисунков видно, что в зависимости
от граничных условий и размера, занимаемого упругим основанием, а также объема жидкости эволюция минимальных частот колебаний к значениям, соответствующим однородному случаю (£ = 1), существенно отличается для каждой из рассмотренных конфигураций. Если для некоторых из них наиболее значительные изменения в частотах происходят лишь тогда, когда уровень заполнения достигает половины длины оболочки, то для других требуется лишь незначительный объем жидкости. При этом рост обоих варьируемых параметров сопровождается интенсивным формоизменением, при котором минимальная частота колебаний, отвечающая различным окружным гармоникам, изменяется в широком интервале. Например, в случае свободно опертых краев окружная мода с минимальной частотой неоднократно чередуется в диапазоне от j = З до j = 9 для третьего варианта неоднородности.
Заключение
Представлены результаты численного исследования собственных частот колебаний вертикальных круговых цилиндрических оболочек, частично заполненных идеальной неподвижной жидкостью. Конструкция полностью окружена двухпараметричеекой упругой средой Пастернака, которая является неоднородной по длине и представляет собой различные комбинации участков, содержащих упругое основание, и участков, где оно отсутствует. Проанализировано влияние разных вариантов неоднородности на минимальные частоты колебаний оболочек, имеющих разнообразные граничные условия при различных значениях коэффициентов постели и уровнях жидкости. Численными экспериментами установлено, что с повышением уровня заполнения оболочек жидкостью влияние упругого основания на частотный спектр конструкции снижается. Продемонстрировано, что, изменяя местоположение, жесткость и размеры упругого основания, можно подобрать такую окружную форму гидроупругих колебаний конструкции, частично заполненной жидкостью, которая в большей степени будет соответствовать требуемым условиям ее эксплуатации.
Благодарности. Работа выполнена в рамках государственного задания (номер гос. регистрационной темы АААА-А19-119012290100-8),
Список литературы
[1] Пастернак П.Л. Основы нового метода расчета фундаментов на упругом основании при помощи двух коэффициентов постели. М.: Госстройиздат; 1954: 56.
[2] Власов В.З., Леонтьев H.H. Балки, плиты и оболочки на упругом основании. М.: Физ-матгиз; 1960: 491.
[3] Горбунов-Посадов М.И., Маликова Т.А. Расчет конструкций на упругом основании. М.: Стройиздат; 1973: 627.
[4] Kerr A.D. Elastic and viscoelastic foundation models. The Journal of Applied Mechanics. 1964; 31(3):491—498. DOLIO. 1115/1.3629667.
[5] Бочкарёв С.А. Собственные колебания цилиндрической оболочки, частично лежащей на упругом основании. Вычислительная механика сплошных сред. 2017; 10(4):406-415. DOÏ: 10.7242/1999-6691/2017.10.4.32.
[61 Shah A.G., Mahmood T., Naeem M.N., Arshad S.H. Vibrational study of fluid-filled functionally graded cylindrical shells resting on elastic foundations. ISRN Mechanical Engineering. 2011; (2011):892460. D01:10.5402/2011/892460.
[7] Shah A.G., Mahmood Т., Naeem M.N., Arshad S.H. Vibration characteristics of fluid-filled cylindrical shells based on elastic foundations. Acta Mechanica. 2011; 216(1—4):17—28. D01:10.1007/s00707-010-0346-l.
[8] Shahbaztabar A., Izadi A., Sadeghian M., Kazemi M. Free vibration analysis of FGM circular cylindrical shells resting on the Pasternak foundation and partially in contact with stationary fluid. Applied Acoustics. 2019; (153):87-101. D01:10.1016/j.apacoust.2019.04.012.
[9] Baghlani A., Khayat M., Dehghan S.M. Free vibration analysis of FGM cylindrical shells surrounded by Pasternak elastic foundation in thermal environment considering fluid-structure interaction. Applied Mathematical Modelling. 2020; (78):550-575. D01:10.1016/j.apm.2019.10.023.
[10] Sheng G.G., Wang X. Thermomechanical vibration analysis of a functionally graded shell with flowing fluid. European Journal of Mechanics — A/Solids. 2008; 27(6):1075-1087. D01:10.1016/j.euromechsol.2008.02.003.
[11] Mohammadi K., Barouti M.M., Safarpour H., Ghadiri M. Effect of distributed axial loading on dynamic stability and buckling analysis of a viscoelastic DWCNT conveying viscous fluid flow. The Journal of the Brazilian Society of Mechanical Sciences and Engineering. 2019; (41):93. DOI:10.1007/s40430-019-1591-4.
[12] Sokolov V., Razov I., Dmitriev A. Influence of the length parameter of an underground oil pipeline on the frequency of free oscillation. E3S Web of Conferences. 2020; (164):03024. D01:10.1051/e3sconf/202016403024.
[13] Amabili M., Garziera R. Vibrations of circular cylindrical shells with nonuniform constraints, elastic bed and added mass; Part I: Empty and fluid-filled shells. Journal of Fluids and Structures. 2000; 14(5):669-690. D01:10.1006/jfls.2000.0288.
[14] Tj H.G., Mikami Т., Kanie S., Sato M. Free vibrations of fluid-filled cylindrical shells on elastic foundations. Thin-Walled Structures. 2005; 43(11):1746-1762. DOI: 10.1016/j.tws.2005.07.005.
[15] Kim Y. Effect of partial elastic foundation on free vibration of fluid-filled functionally graded cylindrical shells. Acta Mechanica Sinica. 2015; (31):920-930. DOI:10.1007/sl0409-015-0442-5.
[16] Park K.J., Kim Y.W. Vibration characteristics of fluid-conveying FGM cylindrical shells resting on Pasternak elastic foundation with an oblique edge. Thin-Walled Structures. 2016; (106):407—419. D01:10.1016/j.tws.2016.05.011.
[17] Lakis A.A., Sinno M. Free vibration of axisvmmetric and beam-like cylindrical shells, partially filled with liquid. International Journal for Numerical Methods in Engineering. 1992; (33):235-268. D01:10.1002/nme.l620330203.
[18] Годунов С.К. Обыкновенные дифференциальные уравнения с постоянными коэффициентами. Т. I: Краевые задачи. Учебное пособие. Новосибирск: НГУ; 1994: 264.
[19] Юдин А.С., Сафроненко В.Г. Виброакустика структурно-неоднородных оболочек. Ростов н/Д: Изд-во ЮФУ; 2013: 424.
[20] Paidoussis М.Р. Fluid-structure interactions: slender structures and axial flow. Vol. 1. 2nd ed. London: Academic Press; 2014: 888.
[21] Paidoussis M.P. Fluid-structure interactions: slender structures and axial flow. Vol. 2. 2nd ed. London: Academic Press; 2016: 942.
[22] Горшков А.Г., Морозов В.И., Пономарев А.Т., Шклярчук Ф.Н. Аэрогидроупру-гость конструкций. М.: Физматлит; 2000: 591.
[23] Шмаков В.П. Избранные труды по гидроупругости и динамике упругих конструкций. М.: Изд-во МГТУ им. Н.Э. Баумана; 2011: 287.
Юдин А.С., Амбалова Н.М. Вынужденные колебания коаксиальных подкрепленных цилиндрических оболочек при взаимодействии с жидкостью. Прикладная механика. 1989; 25(12):63-68.
Кармишин А.В., Лясковец В.А., Мяченков В.И., Фролов А.Н. Статика и динамика тонкостенных оболочечных конструкций. М.: Машиностроение; 1975: 376. Березин И.С., Жидков Н.П. Методы вычислений. Т. 2. М.: Физматлит; 1962: 620. Демидович Б.П., Марон И.А., Шувалова Е.З. Численные методы анализа. М.: Наука; 1967: 368.
Bochkarev S.A. Natural vibrations of a cylindrical shell with fluid partly resting on a two-parameter elastic foundation. International Journal of Structural Stability and Dynamics. 2022; 22(6):2250071. D01:10.1142/S0219455422500717.
Khalifa M.A. Natural frequencies and mode shapes of variable thickness elastic cylindrical shells resting on a Pasternak foundation. Journal of Vibration and Control. 2016; 22(1):37—50. D01:10.1177/1077546314528229.
Amabili M., Paidoussis M.P., Lakis A.A. Vibrations of partially filled cylindrical tanks with ring-stiffeners and flexible bottom. Journal of Sound and Vibration. 1998; 213(2):259-299. DOI: 10.1006/jsvi. 1997.1481.
Kim Y.-W., Lee Y.-S., Ко S.-H. Coupled vibration of partially fluid-filled cylindrical shells with ring stiffeners. Journal of Sound and Vibration. 2004; 276(3-5):869-897. D01:10.1016/j.jsv.2003.08.008.
Maxuch Т., Horacek J., Trnka J., Vesely J. Natural modes and frequencies of a thin clamped-free steel cylindrical storage tank partially filled with water: FEM and measurement. Journal of Sound and Vibration. 1996; 193(3):669-690. D01:10.1006/jsvi.l996.0307. Бочкарёв С.А., Лекомцев С.В., Матвеенко В.П. Численное моделирование пространственных колебаний цилиндрических оболочек, частично заполненных жидкостью. Вычислительные технологии. 2013; 18(2):12—24.
Ergin A., Ugurlu В. Hvdroelastic analysis of fluid storage tanks by using a boundary integral equation method. Journal of Sound and Vibration. 2004; 275:489-513. DOI: 10.1016/j.jsv.2003.07.034.
Кашфутдинов Б.Д., Щеглов Г.А. Валидация свободного программного обеспечения Code^Aster применительно к задаче модального анализа цилиндрической оболочки с жидкостью. Наука и Образование: Научное издание. 2017; (6):101—117. DOI: 10.7463/0517.0001170.
Вычислительные технологии, 2022, том 27, № 4, с. 15-32. © ФИЦ ИВТ, 2022 ISSN 1560-7534
Computational Technologies, 2022, vol. 27, no. 4, pp. 15-32. © FRC ICT, 2022 elSSN 2313-691X
MATHEMATICAL MODELLING
DOI: 10.25743/ICT.2022.27.4.003 Numerical simulation of natural vibrations of a cylindrical shell partially filled with fluid and embedded in an elastic foundation
Bochkarev Sergey A.
Institute of Continuous Media Mechanics of the Ural Branch of Russian Academy of Science, 614068, Perm, Russia
Corresponding author: Bochkarev Sergey A., e-mail: bochkarevQicmm.ru Received December 29, 2021, revised June 10, 2022, accepted June 25, 2022.
Abstract
The paper presents the results of studies on natural vibrations of circular vertical cylindrical shells completely or partially filled with a stationary compressible fluid and embedded in a Pasternak two-parameter elastic foundation. In the meridional direction, the elastic medium is both homogeneous and inhomogeneous, and is an alternation of areas, in which this medium is present or absent. The behavior of the elastic structure and compressible fluid is described within the framework of classical shell theory based on the Kirchhoff— Love hypotheses and the Euler equations. The equations of motion of the shell combined with the corresponding geometric and physical relations are reduced to a system of ordinary differential equations for new unknowns. The acoustic wave equation is also reduced to a system of ordinary differential equations using the straight line method. The formulated boundary value problem is solved by the Godunov orthogonal sweep method involving the numerical integration of differential equations by the fourth order Runge^ Kutta method. The natural frequencies of vibrations are calculated using a combination of stepwise procedure and subsequent refinement by the half-division method. The validity of the obtained results is confirmed by comparing them with the known numerical-analytical solutions. The dependences of minimum vibration frequencies versus the level of fluid are analyzed for simply supported, clamped and cantilevered cylindrical shells of various types of inhomogeneitv along the length of the body and with different stiffness. It is demonstrated that with increase in the level of filling of shells with fluid, the influence of the elastic foundation on the frequency spectrum of the structure decreases.
Keywords: classical shell theory, compressible fluid, Godunov orthogonal sweep method, method of lines, natural vibrations, elastic Pasternak's medium.
Citation: Bochkarev S.A. Numerical simulation of natural vibrations of a cylindrical shell partially filled with fluid and embedded in an elastic foundation. Computational Technologies. 2022; 27(4):15—32. DOL10.25743/ICT.2022.27.4.003. (In Russ.)
Acknowledgements. The results were obtained within the framework of state task (state registration number of the topic AAAA-A19-119012290100-8).
References
1. Pasternak P.L. On a new method of analysis of an elastic foundation by means of two foundation constants fOsnovy novogo metoda rascheta fundamentov na uprugom osnovanii pri pomoshchi dvuh koeffitsientov postelij. Moscow: Gosstroyizdat; 1954: 56. (In Russ.)
2. Vlasov V.Z., Leont'ev N.N. Beams, plates, shells on elastic foundation [Balki, plity, obolochki na uprugom osnovanii]. Moscow: Fizmatgiz; 1960: 491. (In Russ.)
3. Gorbunov-Posadov M.I., Malikova T.A. Calculation of structures on elastic foundation [Raschet konstruktsiy na uprugom osnovanii]. Moscow: Stroyizdat; 1973: 627. (In Russ.)
4. Kerr A.D. Elastic and viscoelastic foundation models. The Journal of Applied Mechanics. 1964; 31(3):491-498. DOLlO.1115/1.3629667.
5. Bochkarev S.A. Free vibrations of a cylindrical shell partially resting on elastic foundation. Journal of Applied Mechanics and Technical Physics. 2018; 59(7):1242-1250. D01:10.1134/S0021894418070039.
6. Shah A.G., Mahmood T., Naeem M.N., Arshad S.H. Vibrational study of fluid-filled functionally graded cylindrical shells resting on elastic foundations. ISRN Mechanical Engineering. 2011; (2011):892460. DOI: 10.5402/2011/892460.
7. Shah A.G., Mahmood T., Naeem M.N., Arshad S.H. Vibration characteristics of fluid-filled cylindrical shells based on elastic foundations. Acta Mechanica. 2011; 21(¡(1 I): 17 2*. D01:10.1007/s00707-010-0346-l.
8. Shahbaztabar A., Izadi A., Sadeghian M., Kazemi M. Free vibration analysis of FGM circular cylindrical shells resting on the Pasternak foundation and partially in contact with stationary fluid. Applied Acoustics. 2019; (153):87-101. D01:10.1016/j.apacoust.2019.04.012.
9. Baghlani A., Khayat M., Dehghan S.M. Free vibration analysis of FGM cylindrical shells surrounded by Pasternak elastic foundation in thermal environment considering fluid-structure interaction. Applied Mathematical Modelling. 2020; (78):550-575. D01:10.1016/j.apm.2019.10.023.
10. Sheng G.G., Wang X. Thermomechanical vibration analysis of a functionally graded shell with flowing fluid. European Journal of Mechanics — A/Solids. 2008; 27(6): 1075-1087. DOI: 10.1016/j.euromechsol.2008.02.003.
11. Mohammadi K., Barouti M.M., Safarpour H., Ghadiri M. Effect of distributed axial loading on dynamic stability and buckling analysis of a viscoelastic DWCNT conveying viscous fluid flow. The Journal of the Brazilian Society of Mechanical Sciences and Engineering. 2019; (41):93. DOI:10.1007/s40430-019-1591-4.
12. Sokolov V., Razov I., Dmitriev A. Influence of the length parameter of an underground oil pipeline on the frequency of free oscillation. E3S Web of Conferences. 2020; (164):03024. D01:10.1051/e3sconf/202016403024.
13. Amabili M., Garziera R. Vibrations of circular cylindrical shells with nonuniform constraints, elastic bed and added mass; Part I: Empty and fluid-filled shells. Journal of Fluids and Structures. 2000; 14(5):669-690. D01:10.1006/jfls.2000.0288.
14. Tj H.G., Mikami T., Kanie S., Sato M. Free vibrations of fluid-filled cylindrical shells on elastic foundations. Thin-Walled Structures. 2005; 43(11):1746-1762. D01:10.1016/j.tws.2005.07.005.
15. Kim Y. Effect of partial elastic foundation on free vibration of fluid-filled functionally graded cylindrical shells. Acta Mechanica Sinica. 2015; (31):920-930. DC)I:10.1007/sl0409-015-0442-5.
16. Park K.J., Kim Y.W. Vibration characteristics of fluid-conveying FGM cylindrical shells resting on Pasternak elastic foundation with an oblique edge. Thin-Walled Structures. 2016; (106):407-419. DOI: 10.1016/j.tws.2016.05.011.
17. Lakis A.A., Sinno M. Free vibration of axisymmetric and beam-like cylindrical shells, partially filled with liquid. International Journal for Numerical Methods in Engineering. 1992; (33):235-268. DOI: 10.1002/nme. 1620330203.
18. Godunov S.K. Ordinary differential equations with constant coefficients. Providence: AMS; 1997: 283.
19. Yudin A.S., Safronenko V.G. Vibroakustika strukturno-neodnorodnykh obolochek [Vibroacoustics of structurally inhomogeneous shells]. Rostov-on-Don: Izd-vo YuFU; 2013: 424. (In Russ.)
20. Paidoussis M.P. Fluid-structure interactions: slender structures and axial flow. Vol. I. 2nd ed. London: Academic Press; 2014: 888.
21. Paidoussis M.P. Fluid-structure interactions: slender structures and axial flow. Vol. 2. 2nd ed. London: Academic Press; 2016: 942.
22. Gorshkov A.G., Morozov V.I., Ponomarev V.I., Shklyarchuk F.N. Aerogidrouprugost' konstruktsiy [Aerohydroelasticity of structures]. Moscow: Fizmatlit; 2000: 591. (In Russ.)
23. Shmakov V.P. Izbrannye trudy po gidrouprugosti i dinamike uprugikh konstruktsiy [Selected works on hydroelasticity and the dynamics of elastic structures]. Moscow: Izd-vo MGTU im. N.E. Baumana; 2011: 287. (In Russ.)
24. Yudin A.S., Ambalova N.M. Forced vibrations of coaxial reinforced cylindrical shells during interaction with a fluid. International Applied Mechanics. 1989; 25(12): 1222 1227. D01:10.1007/BF00887148.
25. Karmishin A.V., Lyaskovets V.A., Myachenkov V.I., Frolov A.N. Statika i dinamika tonkostennykh obolochechnykh konstruktsiy [The statics and dynamics of thin-walled shell structures]. Moscow: Mashinostroenie; 1975: 376. (In Russ.)
26. Berezin I.S., Zhidkov N.P. Computing methods. Vol. 2. N.Y.: Pergamon Press; 1965: 680.
27. Demidovich B.P., Maron I.A., Shuvalova E.Z. Chislennye metody analiza [Numerical methods of analysis. Approximation of functions, differential and integral equations]. Moscow: Nauka; 1967: 368. (In Russ.)
28. Bochkarev S. A. Natural vibrations of a cylindrical shell with fluid partly resting on a two-parameter elastic foundation. International Journal of Structural Stability and Dynamics. 2022; 22(6):2250071. D01:10.1142/S0219455422500717.
29. Khalifa M.A. Natural frequencies and mode shapes of variable thickness elastic cylindrical shells resting on a Pasternak foundation. Journal of Vibration and Control. 2016; 22(l):37-50. DOI: 10.1177/1077546314528229.
30. Amabili M., Paidoussis M.P., Lakis A. A. Vibrations of partially filled cylindrical tanks with ring-stiffeners and flexible bottom. Journal of Sound and Vibration. 1998; 213(2):259-299. DOI: 10.1006/jsvi. 1997.1481.
31. Kim Y.-W., Lee Y.-S., Ko S.-H. Coupled vibration of partially fluid-filled cylindrical shells with ring stiffeners. Journal of Sound and Vibration. 2004; 276(3-5) :869-897. D01:10.1016/j.jsv.2003.08.008.
32. Maxiich T., Horacek J., Trnka J., Vesely J. Natural modes and frequencies of a thin clamped-free steel cylindrical storage tank partially filled with water: FEM and measurement. Journal of Sound and Vibration. 1996; 193(3):669-690. DOI: 10.1006/jsvi. 1996.0307.
33. Bochkarev S.A., Lekomtsev S.V., Matveenko V.P. Numerical modeling of spatial vibrations of cylindrical shells, partially filled with fluid. Computational Technologies. 2013; 18(2):12-24. (In Russ.)
34. Ergin A., Ugurlu B. Hydroelastic analysis of fluid storage tanks by using a boundary integral equation method. Journal of Sound and Vibration. 2004; 275:489-513. D01:10.1016/j.jsv.2003.07.034.
35. Kashfutdinov B.D., Shcheglov G.A. Validation of the open source Code^Aster software used in the modal analysis of the fluid-filled cylindrical shell. Science and Education of the Bauman MSTU. 2017; (6):101-117. D01:10.7463/0617.0001170. (In Russ.)