CC BY
76
Бочкарев С.А., Лекомцев С.В. Собственные колебания нагретых функционально-градиентных цилиндрических оболочек, содержащих жидкость // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2015. - № 4. - С. 19-35. DOI: 10.15593/perm.mech/2015.4.02
Bochkarev S.A., Lekomtsev S.V. Natural vibrations of heated functionally graded cylindrical shells with fluid. PNRPU Mechanics Bulletin. 2015. No. 4. Рр. 19-35. DOI: 10.15593/perm.mech/2015.4.02
ВЕСТНИК ПНИПУ. МЕХАНИКА № 4,2015
PNRPU MECHANICS BULLETIN
http://vestnik.pstu.ru/mechanics/about/inf/
DOI: 10.15593/perm.mech/2015.4.02 УДК 539.3
СОБСТВЕННЫЕ КОЛЕБАНИЯ НАГРЕТЫХ ФУНКЦИОНАЛЬНО-ГРАДИЕНТНЫХ ЦИЛИНДРИЧЕСКИХ ОБОЛОЧЕК, СОДЕРЖАЩИХ ЖИДКОСТЬ
С.А. Бочкарев, С.В. Лекомцев
Институт механики сплошных сред УрО РАН, Пермь, Россия
О СТАТЬЕ
АННОТАЦИЯ
Получена: 11 июля 2015 г. Принята: 29 октября 2015 г. Опубликована: 25 декабря 2015 г.
Ключевые слова: цилиндрическая оболочка, функционально-градиентный материал,температурное нагружение, потенциальная жидкость, метод конечных элементов, собственные колебания
Представлены результаты исследований собственных колебаний нагретых круговых цилиндрических оболочек, содержащих неподвижную идеальную жидкость и выполненных из функционально-градиентных материалов. Температурнозависимые эффективные свойства материала, состоящего из смеси оксида циркония и сплава титана, изменяются по толщине оболочки согласно степенному закону. Распределение температуры по радиальной координате определяется из решения квазилинейного одномерного уравнения теплопроводности. В качестве математической формулировки задачи динамики оболочки используется классическая теория оболочек и принцип возможных перемещений. Поведение жидкости описывается в рамках потенциальной теории. Соответствующее волновое уравнение совместно с условием непроницаемости и граничными условиями преобразуются к системе уравнений с использованием метода Бубнова-Галеркина. В результате решение задачи, осуществляемое с помощью полуаналитического варианта метода конечных элементов, сводится к вычислению комплексных собственных значений связанной системы уравнений. Достоверность результатов, получаемых в рамках разработанного алгоритма, оценена путем сравнения с известными численно-аналитическими решениями. Для круговых цилиндрических оболочек с различными вариантами граничных условий продемонстрированы зависимости минимальных частот колебаний от температурной нагрузки при разных консистенциях функционально-градиентного материала. Определены критические значения температур при разных условиях нагрева и геометрических размерах. Проанализировано различие, оказываемое нагревом на динамические свойства пустых оболочек и оболочек, взаимодействующих с жидкостью. Показано, что наличие жидкости внутри оболочки проявляет наиболее существенное влияние на колебательные характеристики нагретых консольно закрепленных оболочек.
©ПНИПУ
© Бочкарев Сергей Аркадьевич - кандидат физико-математических наук, старший научный сотрудник, e-mail: [email protected]
Лекомцев Сергей Владимирович - кандидат физико-математических наук, научный сотрудник, e-mail: [email protected]
Sergey A. Bochkarev - PhD in Physical and Mathematical Sciences, Senior Researcher, e-mail: [email protected] Sergey V. Lekomtsev - PhD in Physical and Mathematical Sciences, Researcher, e-mail: [email protected]
19
Bochkarev S.A., Lekomtsev S.V. /PNRPUMechanics Bulletin 4 (2015) 19-35
NATURAL VIBRATIONS OF HEATED FUNCTIONALLY GRADED CYLINDRICAL SHELLS WITH FLUID
S.A. Bochkarev, S.V. Lekomtsev
Institute of Continuous Media Mechanics RAS, Perm, Russian Federation
ARTICLE INFO
ABSTRACT
Received: 11 July 2015 Accepted: 29 October 2015 Published: 25 December 2015
Keywords:
cylindrical shell, functionally graded material, thermal loading, potential fluid, finite element method, natural vibrations
The paper presents the results of studying the natural vibrations of heated cylindrical circular FGM shells containing a quiescent ideal liquid. The temperature dependent effective properties of the material representing a mixture of zirconium oxide and titanium alloy vary through the thickness of the shell according to the power law. The distribution of temperature along the radial coordinate is determined by solving a quasi-linear onedimensional heat conduction equation. A mathematical formulation of the problem is based on the classical theory of shells and the principle of virtual displacements. The behavior of the liquid is described in the framework of the potential theory. The corresponding wave equation together with the impermeability condition and boundary conditions are transformed into a system of equations using the Bubnov-Galerkin method. As a result, the solution of the problem, which is sought with the use of semi-analytical version of the finite element method, reduces to calculations of complex eigenvalues of the coupled system of equations. The reliability of the results, obtained by application of the developed algorithm, is verified through a comparison with the known numerical-analytical solutions. The data obtained for circular cylindrical shells with different boundary conditions have revealed the dependence of the minimal vibration frequency on temperature at different volume fractions of FGM. The critical values of temperature have been determined for different heating regimes and geometrical dimensions. The difference between the dynamic properties of empty and liquid-containing shells caused by heating has been analyzed. It has been shown that in the case of cantilevered shells the presence of liquid inside the shell exerts the most notable effect on the vibrational behavior of the system.
©PNRPU
Введение
При всем многообразии численно-аналитических и численных работ, посвященных исследованию конструкций, выполненных из функционально-градиентных (ФГ) материалов [1], только в незначительном числе статей затрагиваются вопросы их взаимодействия с газообразной [2-6], а также неподвижной [7-10] или текущей [11-13] жидкой средой. Такой ограниченный список источников не является очевидным, поскольку высокие жаропрочные свойства, обусловленные использованием керамики в качестве одного из составляющих, подразумевают применение ФГ-материалов в авиационной и ракетной технике, где упругие конструкции взаимодействуют с нагретым потоком жидкости или газа. Только в некоторых из вышеперечисленных работ учитывается температурное нагружение, а свойства материала принимаются зависимыми от температуры. Однако для упрощения используется постановка, в которой коэффициенты теплопроводности не зависят от температуры. Такое допущение позволяет использовать аналитическое выражение для нахождения распределения температуры по толщине оболочки. Более точная формулировка требует непосредственного численного решения квазилинейного уравнения теплопроводности. Несмотря на то, что получаемые результаты достовернее описывают поведение реальных машиностроительных конструкций, этот подход распространен в меньшей степени [14-17].
Ни в одной из указанных работ, где рассматриваются ФГ-оболочки, содержащие неподвижную жидкость, не принимаются во внимание температурные нагрузки. Влияние
20
Бочкарев С.А., Лекомцев С.В. /Вестник ПНИПУ. Механика 4 (2015) 19-35
различных параметров на собственные частоты колебаний свободно опертых ФГ-обо-лочек с жидкостью исследовалось в [7] при помощи трехмерных уравнений теории упругости. В работе [8] оценено влияние различных комбинаций граничных условий на собственные частоты колебаний оболочки, заполненной жидкостью и выполненной из стали и никеля. Аналогичный подход используется в [9] для анализа расположенной на упругом основании оболочки с жидкостью. Представленные в этих статьях численно-аналитические решения применимы лишь к цилиндрическим оболочкам, содержащим несжимаемую жидкость, потому что гидродинамическое давление вычисляется с помощью аналитического выражения, записанного относительно функций Бесселя. Более общая формулировка, основанная на использовании метода конечных элементов, предложена в работе [10]. Здесь исследовалось влияние различных комбинаций граничных условий на собственные частоты колебаний оболочек с различными геометрическими размерами. Таким образом, влияние температурной нагрузки на спектр собственных частот цилиндрических оболочек, содержащих неподвижную жидкость, остается не изученным. В связи с этим нашей целью является проведение такого анализа при различных комбинациях кинематических граничных условий. Моделирование жидкости осуществляется в рамках потенциальной теории, согласно которой единственной искомой величиной является потенциал скорости. В результате решение задачи сводится к вычислению комплексных собственных значений несимметричной матрицы. Отметим, что незначительная модификация алгоритма позволяет учитывать как течение, так и вращение невязкой сжимаемой жидкости или газа [18-19].
1. Постановка задачи и основные соотношения
Рассматривается выполненная из функционально-градиентного материала упругая цилиндрическая оболочка (рис. 1) длиной L и радиусом R, внутри которой находится идеальная сжимаемая жидкость плотностью р/. Одна поверхность оболочки, например наружная, нагревается до заданной температуры To, тогда как температура другой всегда остается равной T = 300 K. Целью работы является исследование влияния свойств ФГ-материала на собственные колебания нагретой оболочки при различных вариантах кинематических граничных условий и линейных размерах конструкции.
z \
Рис. 1. Содержащая жидкость цилиндрическая оболочка, выполненная из функционально-градиентного материала
21
Bochkarev S.A., Lekomtsev S.V. /PNRPUMechanics Bulletin 4 (2015) 19-35
Для описания неподвижной жидкости в области Vf вводится в рассмотрение потенциал скорости ф, который в цилиндрической системе координат (r, 0, x) в случае малых возмущений описывается волновым уравнением [20]
V2„ д2ф , 1 S2ф, д2ф , 1 5ф V ф —----;■ + ;---;■ +-;■ +-----
dr2 r2 д02 dx2 r dr
1 д2ф
2 д.2 ,
c dt
(1)
где с - скорость звука в жидкости. Давление жидкости Ру на упругую конструкцию вычисляется по линеаризованной формуле Бернулли
Р^ — —ру дф на поверхности Sa — Sy О Ss.
(2)
Здесь s - меридиональная координата оболочки; Sу, Ss - площади, ограничивающие объемы жидкости и оболочки соответственно. На поверхности раздела оболочка-жидкость Sa задается условие непроницаемости
0ф —Sw, (3)
dn dt
где n - нормаль к поверхности; w - нормальная составляющая вектора перемещений оболочки. Потенциал скорости подчиняется следующим граничным условиям:
x — 0: ф — 0, x — L : dф/dx — 0. (4)
Для дальнейшей численной реализации задачи на основе полуаналитического варианта метода конечных элементов (МКЭ) дифференциальное уравнение для потенциала скорости (1) с граничными условиями (3) и (4) преобразуются с помощью метода Бубно-ва-Галеркина [18].
Компоненты вектора деформации в криволинейной системе координат (s, 0, z) определяются с помощью классической теории оболочек, основанной на гипотезах Кирхгофа-Лява. Для эквидистантной поверхности они могут быть записаны следующим образом [21]:
£11 — Е11 +zk11, ^22 — Е22 + zk22, £12 — Е12 +zk12,
где
и
Е11 — е1 + 1/20J2, k11 —k1, (1 о2), Е12 — е12 +0102, k12 —2т
^1 —
du
ds
S 2 —
1 С dv л
----+ w
\d0 j
R
dv 1 du dw
Sl2 I , 01 ■
12 ds R d0 1 ds
dw ^ , d2 w k — 1 dv d2 w ^ 1 dv 2 d w
d0 j , k1 — Я 2 , ds 1 2 — R2 [d0 d02 j , T— R yds dsd0 у
(5)
02 —-
2 R
Здесь u, v - меридиональная и окружная составляющие вектора перемещений; 0г- - углы поворота недеформируемой нормали.
Компоненты деформаций срединной поверхности (5) в матричной форме можно представить в виде
22
Бочкарев С.А., Лекомцев С.В. /Вестник ПНИПУ. Механика 4 (2015) 19-35
£ — £* +1/2 Ее,
где £ = {Е11,Е22,Е12,k11,к22,k12}T, £* = {{,е2,е12,k1,k2,2x}T - линейная часть деформации,
e = {O,O,O,O,01,02,O,O,O}T; E - матрица линейных множителей, которая имеет ненулевые
компоненты Е15 = 01, Е\6 = 02, Е35 = 02, Е36 = 01-
Физические соотношения, устанавливающие связь между вектором усилий и моментов T — {1, 722, 712 ’ M11 ’ M22 ’ M12 } и вектором деформаций е, запишем в матричном виде:
T — De- Tth
A B If Tth 1 b cje [MthJ,
(6)
где Tth — {7h, 72th ,O} и Mth = {м^, M22,O} - векторы температурных усилий и моментов. Коэффициенты, входящие в матрицу жесткостей D, и температурные компоненты определяются как
(,bij,cij) —Я1,z,z2)Qijdz, (г',j —1,2,3), (7)
h
(t*. Mth ) = J(1>z)gsaArdz, a = {a,# (z,r),a,f (z, 7),O}T, (8)
h
где
Q11 — Q22 — Eeff (z,7)/ [1 - Vef (z.7)] . Q12 — Q21 — Vef (z, 7)Q11'
Q33 — Eff (z,7)/2/[1 + Vf (z,7)].
Здесь A 7 - перепад температуры относительно значения, не приводящего к появлению начальных напряжений (3OO K); h - толщина оболочки; Ef (z, 7), veff (z, 7), aef (z, 7) -
эффективные модуль упругости, коэффициенты Пуассона и температурного расширения материала. Для ФГ-материала эти характеристики определяются свойствами составляющих его материалов P и их объемными долями V . В случае функционально-градиентного материала, представляющего собой смесь металла m и керамики c [22],
Peff — PmVm + PcVc-
В случае оболочки постоянной толщины объемные доли являются функциями радиальной координаты z, отсчитываемой от срединной поверхности, и изменяются согласно степенному закону
V J 2 z+h
\N
2h
V = 1 - V
r m 1 r c’
где N - показатель объемной доли, изменяющийся от нуля до бесконечности. С учетом последнего представления эффективные характеристики, под которыми, помимо перечисленных выше, понимаются также плотность материала peff (z,7) и коэффициент теплопроводности kf (z, 7), вычисляются следующим образом:
23
Bochkarev S.A., Lekomtsev S.V. /PNRPUMechanics Bulletin 4 (2015) 19-35
Pff (z, T) = Pm (T) + [Pc (T) - Pm ( T)] (^ .
Температурно-зависимые свойства материала определяются как
P(T) = P0 (p-1t 1 +1 + P1T + P2T2 + P3T3).
Здесь P0,P-1,P,P2,P3 - коэффициенты температуры T, уникальные для каждого материала. Распределение температуры по толщине оболочки находится из численного решения квазилинейного уравнения установившейся теплопроводности с соответствующими граничными условиями
= 0, T \z=-h/2= Ti, T |z=h/2= To . (9)
Для математической формулировки задачи используется принцип возможных перемещений. С дополнительным слагаемым, учитывающим предварительное нагружение, его можно представить в матричной форме [23]:
j SeT De„ dS + jSd T p0ddS + jSe T о °edS - jSd T PdS = 0. (10)
S S S S
Здесь d и P = |0 0 P^ | - векторы обобщенных перемещений и поверхностных нагрузок, p0 =jpeff (z,T)dz . Матрица о0, элементы которой находятся из условия
ET D£° = о0е, имеет следующие ненулевые компоненты: а°5 = T11, а^6 = T22. Здесь верхний индекс «0» означает начальное равновесное состояние, а усилия Tu вычисляются по формуле (6).
2. Численная реализация
_d_
dz
kf (z,T )dT
Численное решение задачи осуществляется с использованием полуаналитического варианта МКЭ [24], основанного на представлении решения в виде ряда Фурье по окружной координате 0:
р, "W ф) = Х”=0 (uj, wj, фj )cos j0, v = SJ=o vj sin j0, (11)
где j - номер гармоники.
Выражая в (11) искомые переменные через их узловые значения, получим для произвольных конечных элементов оболочки и жидкости
U = {, v, W = Nde, (12)
ф=ГТ, (13)
где N и F - матрицы функций формы оболочечного конечного элемента и потенциала скорости; de и fe - векторы узловых значений. Для оболочек используется конечный элемент в виде усеченного конуса с аппроксимацией меридиональной и окружной компонент вектора перемещений линейным полиномом, а нормальной компоненты - кубическим полиномом. Для жидкости применяется треугольный конечный элемент с линейной аппроксимацией потенциала скорости. С учетом (12)
24
Бочкарев С.А., Лекомцев С.В. /Вестник ПНИПУ. Механика 4 (2015) 19-35
£*= Bdе, e = Gdе
(14)
где B и G - матрицы связи деформаций £* и e с узловыми значениями.
Применяя стандартные процедуры МКЭ, получим из преобразованного волнового уравнения для потенциала скорости и принципа возможных перемещений (10) с учетом (12)—(14) связанную систему двух уравнений. Ее можно записать в матричном виде, представляя возмущенное движение оболочки и жидкости как (d, ф) = (q, f) exp(i*wt).
(К -ю2М + ГюС )q, f )T = 0,
K = К , + K * 0 ' , M = "M , 0 ' , c = " 0 Cf"
0 K f _ 0 M f _ 1 о 'ft, и 1
K s = Xj BT DBdS, M s = Xj NT p0NdS, M f = Xj 1FT FdV,
ms Ss
K f = U
m ss
dFT dF 1 dFT dF dFT dFЛ
mfVf
- + -
■ + ■
m V y dr dr r2 50 50 dx dx j
f Vf
dV, C/s=-X|fT NdS,
ms Sn
Cf =Pf X j NTFdS, K* = X j GTa0GdS.
ms S„
ms Ss
(15)
Здесь mf и ms — число конечных элементов, на которые разбиваются области жидкости
и оболочки; N — матрица функций формы нормальной составляющей вектора перемещения оболочки; ю — характеристический показатель, в общем случае являющийся комплексным; i* = V—1 .
Таким образом, задача о свободных колебаниях нагретой цилиндрической ФГ-оболочки, содержащей неподвижную жидкость, сводится к выполнению следующей последовательности шагов:
• из решения квазилинейного уравнения (9) определяется распределение температуры по толщине оболочки AT(z) при заданных граничных условиях на внутренней Ti и наружной To поверхностях. Соответствующий алгоритм подробно описан в [14];
• на основе полученного распределения температуры по толщине оболочки вычис-
v* rwith
ляется вектор температурных усилий и моментов T ;
• из решения осесимметричной статической задачи Ksd = P0 , где P0e = l B TTthdS ,
*’Ss
по формуле (6) формируются компоненты матрицы начальных усилий и моментов a0;
• в результате применения итерационного алгоритма на основе метода Мюллера [25], из решения системы (15) определяются и анализируются собственные частоты ю.
3. Результаты расчетов
В численных примерах рассматривается цилиндрическая функционально-градиентная оболочка (h = 1х10—3 м), выполненная из оксида циркония и сплава титана. Физикомеханические характеристики составляющих материалов приведены в табл. 1. Расчеты выполнялись для нескольких конфигураций (табл. 2), в которых чередуются базовые материалы и условия нагрева на внутренней и наружной поверхностях оболочки. Вычисле-
25
Bochkarev S.A., Lekomtsev S.V. /PNRPUMechanics Bulletin 4 (2015) 19-35
ния осуществлялись для оболочек, один край которых всегда жестко закреплен (u = v = w = dw/ds = 0, C). На другом краю задаются следующие граничные условия: свободное опирание (v = w = 0, S), жесткое закрепление, свободный край (F). Численными экспериментами было установлено, что для достижения приемлемой точности вычислений достаточно 40 оболочечных конечных элементов и 1000 элементов для жидкости.
Таблица 1
Физико-механические свойства составляющих материалов и их зависимость
от температуры
Свойства P0 P-1 P1 P 2 P3
Оксид циркония, ZrO2 [26]
р, кг/м3 5,7000х103 0 0 0 0
E, Па 2,4427x10" 0 —1,3707х 10-3 1,2139х10-6 -3,6814х 10-10
V 0,2882 0 1,1335 х 10-4 0 0
а, 1/К 1,2766х10-5 0 -1,4900х10-3 1,0х10-6 -0,6775х10-11
k, Вт/мК 1,7 0 1,2760х10-4 6,6485х10-6 0
Сплав титана, Ti-6Al-4V [26]
р, кг/м3 4,4290х103 0 0 0 0
E, Па 1,2256x10" 0 -4,5864х10-4 0 0
V 0,28838 0 1,1214х10-4 0 0
а, 1/К 7,5788х10-6 0 6,5000х10-4 3,1467х10-7 0
k, Вт/мК 1,2095 0 1,3938х10-2 0 0
Сталь, SUS304 [8]
р, кг/м3 8166 0 0 0 0
E, Па 2,0104x10" 0 3,079х10-4 -6,534х10-7 0
V 0,3262 0 -2,002х10-4 3,797х10-7 0
Никель, [8]
р, кг/м3 8900 0 0 0 0
E, Па 2,2395x10" 0 -2,794х10-4 -3,998х10-9 0
V 0,31 0 0 0 0
Таблица 2
Варианты расчетов
Конфигурация I Конфигурация II Конфигурация III Конфигурация IV
ZrO2 / Ti-6Al-4V нагрев снизу ZrO2 / Ti-6Al-4V нагрев сверху Ti-6Al-4V / ZrO2 нагрев снизу Ti-6Al-4V / ZrO2 нагрев сверху
3.1. Тестирование алгоритма
Тестирование разработанного алгоритма осуществлено путем сравнения с результатами работы [8], в которой собственные частоты колебаний заполненной жидкостью ФГ-оболочки (L/R = 20, R/h = 500) получены волновым методом. Оболочка выполнена из стали (внутренняя поверхность) и никеля (наружная поверхность). Свойства материалов не зависят от температуры и определяются из табл. 1 при 300 K. Сравнение двух низших частот колебаний (Гц) для двух гармоник в случае пустой и заполненной жидкостью цилиндрической оболочки для различных значений показателя объемной доли N и граничных условий представлено в табл. 3. Здесь через m обозначено число полуволн в мери-
26
Бочкарев С.А., Лекомцев С.В. /Вестник ПНИПУ. Механика 4 (2015) 19-35
диональном направлении. В верхней строке приведены данные из [8], а в нижней строке -полученные в настоящей работе. Результаты демонстрируют хорошее совпадение с известными данными, особенно в случае свободно опертых оболочек, где различие не превышает 2 %.
Таблица 3
Сравнение собственных частот пустой и содержащей жидкость ФГ-оболочки при различных граничных условиях
m ] N = 0,3 N = 1 N = 5
Пустая С жидкостью Пустая С жидкостью Пустая С жидкостью
Граничные условия SS
1 2 4,4160 0,9031 4,4737 0,9049 4,5498 0,9076
4,4159 0,9137 4,4736 0,9155 4,5497 0,9188
3 4.0940 0,9599 4,1478 0,9620 4,2183 0,9650
4,0940 0,9656 4,1478 0,9678 4,2183 0,9708
2 2 16,628 3,3988 16,848 3,4062 17,132 3,4158
16,628 3,4676 16,848 3,4751 17,132 3,4848
3 8,7166 2,0419 8,8299 2,0462 8,9807 2,0527
8,7165 2,0929 8,8298 2,0973 8,9805 2,1038
Граничные условия CC
1 2 9,5233 1,9473 9,6487 1,9514 9,8121 1,9569
9,3979 1,9511 9,5219 1,9553 9,6824 1,9607
3 5,7422 1,3458 5,8166 1,3486 5,9161 1,3529
5,7325 1,3682 5,8080 1,3713 5,9065 1,3755
2 2 25,549 5,2189 25,888 5,2303 26,325 5,2448
24,674 5,1588 24,999 5,1699 25,420 5,1840
3 12,847 3,0073 13,015 3,0138 13,236 3,0232
12,653 3,0567 12,818 3,0635 13,036 3,0728
Граничные условия CS
1 2 6,7063 1,3714 6,7943 1,3742 6,9096 1,3782
6,6459 1,3778 6,7334 1,3808 6,8474 1,3847
3 4,7475 1,1129 4,8093 1,1153 4,8914 1,1188
4,7372 1,1238 4,7994 1,1263 4,8810 1,1298
2 2 20,873 4,2651 21,149 4,2744 21,506 4,2864
20,542 4,2895 20,813 4,2988 21,164 4,3106
3 10,648 2,4935 10,787 2,4988 10,971 2,5067
10,578 2,5484 10,716 2,5539 10,898 2,5617
3.2. Собственные колебания нагретой оболочки без жидкости
Для оценки влияния жидкости на динамические характеристики ФГ-оболочки целесообразно проанализировать влияние температурной нагрузки на собственные частоты колебаний пустой оболочки. Внешнее воздействие повышенной температуры, осуществляемое, например, в результате контакта с горячим потоком жидкости или газа, оказывает дестабилизирующее влияние на упругое тело и может привести к потере устойчивости. В табл. 4 приведены значения минимальной температуры нагрева T (K), при которой низшая частота колебаний оболочки ю становится равной нулю и, следовательно, осуществляется потеря устойчивости дивергентного типа. Здесь представлены данные, соответствующие четырем конфигурациям (см. табл. 3), полученные для двух вариантов граничных условий. Показанные в табл. 4 результаты демонстрируют известные
27
Bochkarev S.A., Lekomtsev S.V. /PNRPUMechanics Bulletin 4 (2015) 19-35
закономерности, имеющие место при исследовании выполненных из ФГ-материалов конструкций. С увеличением показателя объемной доли N критические значения температуры возрастают или, наоборот, уменьшаются, в зависимости от того, на какой поверхности располагается более жесткий материал. В некоторых случаях критические температуры, соответствующие промежуточным значениям N, ограниченны значениями, полученными для чистых составляющих материалов (т. е. для N = 0 и N = да).
Таблица 4
Критические значения температуры T (K) в зависимости от объемной доли N для пустых оболочек с различными граничными условиями
Конфигурация N = 0 N = 0,3 N = 1 N = 5 N = 10 N = 30 N = да
Граничные условия CC
I 767,76 850,61 966,23 1079,1 1089,9 1091,6 1089,6
II 771,07 798,34 825,47 956,63 1013,29 1065,5 1095,5
III 1089,3 925,97 821,32 773,14 771,37 770,24 767,84
IV 1096,5 1046,5 972,39 858,28 823,57 791,49 771,77
Граничные условия CF
I 1714,7 1814,0 1874,2 1871,3 1858,5 1853,6 1858,2
II 978,23 986,15 1001,4 1134,8 1220,0 1302,1 1346,1
III 1858,6 1798,4 1722,6 1698,9 1716,9 1725,3 1715,4
IV 1347,5 1350,6 1340,4 1216,1 1133,3 1040,9 980,34
Анализ данных табл. 4 позволяет выявить некоторые специфические особенности, обусловленные или комбинацией рассматриваемых материалов, или заданным вариантом граничных условий. Например, для некоторых конфигураций имеется определенный диапазон показателя объемной доли N, на котором критические значения температур превышают значения для чистого материала с наибольшей жесткостью (т.е. керамики). Аналогичные зависимости были получены в [15] при исследовании нагретых пустых конических оболочек в случае нескольких комбинаций материалов, в частности, сталь-нитрид кремния. Возможного объяснения такого нехарактерного поведения предложено не было. Хотя в работе [15], как и в настоящем исследовании, рассматривались материалы, свойства которых зависят от температуры, нетипичные зависимости были получены авторами настоящего исследования и в том случае, когда свойства материалов не зависят от температуры. Более того, подобные особенности наблюдались, если вместо нелинейного распределения температуры по толщине оболочки задавалось линейное или однородное распределения. Дополнительные расчеты показали, что в решениях задачи теплопроводности (9), статической задачи, и в распределении температурных усилий и моментов Tth по толщине оболочки отсутствуют закономерности, которые могли бы однозначно характеризовать ту или иную зависимость. Было установлено, что нехарактерное поведение проявляется в том случае, когда материалы имеют существенно различные по величине коэффициенты температурного расширения (как для двух первых материалов в табл. 2). В случае материалов с близкими значениями такая зависимость отсутствует для любой из возможных конфигураций. Немаловажную роль играют также и непосредственные физико-механические свойства некоторых материалов: значительная жесткость приводит к повышению порога потери устойчивости, а увеличение коэффициента температурного расширения, наоборот, способствует его снижению [15]. В результате объединения этих двух факторов возникают предпосылки для появления некоторых специфических зависи-
28
Бочкарев С.А., Лекомцев С.В. /Вестник ПНИПУ. Механика 4 (2015) 19-35
мостей в изменении критических температур от консистенции ФГ-материала, которые повышают эксплуатационные свойства выполненных из них конструкций.
Другой особенностью, отраженной в табл. 4, являются результаты, полученные для консольных оболочек. Для этой комбинации граничных условий изменение стороны нагрева приводит к значительному росту критических температур потери устойчивости. Для рассматриваемых в работе идеализированных условий нагрева, при которых температура на одной из поверхностей остается неизменной, смена нагреваемой поверхности приводит к изменению знака температурных моментов Mth. Если в случае жесткого закрепления это ведет к незначительному изменению критических температур, то для консольного закрепления изменение знака способствует большему раскрытию свободного края оболочки и как результат - существенному снижению порога устойчивости.
Из табл. 4 можно также увидеть, что для одного и того же материала в зависимости от того, на какой поверхности он располагается, имеет место незначительное расхождение в критических температурах. Это различие обусловлено тем, что в случае N = да рассматривается конечное значение показателя объемной доли. Поэтому при вычислении жесткостных характеристик, усилий и моментов по формулам (7)-(8) при z = ± h/2 принимаются во внимание свойства другого материала.
3.3. Собственные колебания нагретой оболочки с жидкостью
На рис. 2, а для жестко закрепленной ФГ-оболочки, внутри которой содержится неподвижная жидкость (р/ = 1000 кг/м3, c = 1500 м/с), представлены зависимости низших частот колебаний ю (Гц) от номера гармоники j, полученные при температуре T = 300 K и различных значениях объемных долей N. Здесь показана уже упомянутая ранее типичная картина, имеющая место при исследовании функционально-градиентного материала, а именно: вибрационные свойства конструкции практически всегда ограничены предельными значениями, полученными для чистых материалов. Кривые частотного спектра, вычисленные с произвольным показателем объемной доли N, располагаются между кривыми, полученными для керамики и металла (N = 0). Таким образом, содержащаяся внутри оболочки жидкость не оказывает качественного влияния на характер динамического поведения конструкций, выполненных из ФГ-материала.
Как и в случае пустой оболочки, повышение температуры окружающей среды приводит к снижению собственных частот колебаний конструкций, содержащих жидкость. Это демонстрируется данными, приведенными на рис. 2, б. Здесь для жестко закрепленной ФГ-оболочки (N = 1) показаны зависимости минимальных частот колебаний ю (Гц) от номера гармоники j, полученные при различных значениях температуры T (K) внутренней поверхности. Из представленных данных видно, что последовательное ее повышение приводит к постепенному снижению собственных частот ю до тех пор, пока они резко не обращаются в ноль при определенных гармониках. Более подробно это явление демонстрируется на рис. 3, а. Здесь показано изменение собственной частоты колебаний ю от температуры T (K) для гармоник j = 5, 13 и 14. Из рисунка видно, что частота пятой гармоники является минимальной на достаточно широком диапазоне температур, пока собственные частоты более высоких гармоник не достигнут меньших значений. Это происходит благодаря тому, что они сильнее реагируют на возрастание температуры. Такая зависимость в целом определяет характер изменения собственных частот колебаний от
29
Bochkarev S.A., Lekomtsev S.V. /PNRPUMechanics Bulletin 4 (2015) 19-35
температуры для оболочек с некоторыми комбинациями граничных условий, что демонстрируется результатами, показанными на рис. 3, б и 4. На этих рисунках приведены зависимости минимальных частот колебаний ю (Гц) от температуры T (К), полученные для оболочек с жидкостью, при разных вариантах граничных условий. Результаты для оболочек с граничными условиями CF (см. рис. 3, б) и CC (рис. 4, а), вычисленные при различных значениях показателя объемной доли N, демонстрируют резкое снижение частоты при достижении некоторого температурного порога, обусловленного, как показано ранее, изменением гармоники с минимальной частотой колебаний. Отметим, что вне зависимости от граничных условий потеря устойчивости осуществляется при j = 14, тогда как в случае пустых оболочек - при j = 13. Для оболочек с граничными условиями CS (рис. 4, б) рост температуры приводит лишь к незначительному уменьшению частоты и потери устойчивости не происходит. Это обусловлено отсутствием закрепления по осевой координате на правом краю оболочки. Тело имеет возможность свободно расширяться в этом направлении, как и в случае граничных условий CF, что препятствует появлению значительных напряжений, возникающих в случае жесткого закрепления обоих краев.
Рис. 2. Зависимости минимальных частот колебаний ю от номера гармоники j для жестко закрепленной ФГ-оболочки, содержащей жидкость, при различных температурах T (K); конфигурация I: а - T = 300 K; б - N = 1
Более подробные результаты исследований ФГ-оболочек для двух вариантов граничных условий и различных конфигураций сведены в табл. 5. Представленные здесь критические значения температур T (K), по сравнению с полученными для пустой оболочки данными (табл. 4), позволяют сделать вывод о том, что наличие жидкости внутри конструкции, помимо снижения частотного спектра, может оказывать влияние и на характер зависимости динамических свойств конструкций из ФГ-материала от температуры. В частности, наличие жидкости для некоторых конфигураций приводит к расширению диапазона показателя объемных долей, на котором проявляется нетипичная температурная зависимость (CC, конфигурация I), или даже снижает границы устойчивости, определенные значениями для чистых материалов (CC, III).
30
Бочкарев С.А., Лекомцев С.В. /Вестник ПНИПУ. Механика 4 (2015) 19-35
Рис. 3. Зависимости частот колебаний ю от температуры T для консольной ФГ-оболочки при N = 0,3 (а) и разных значениях показателя объемной доли N (б); конфигурация I
Рис. 4. Зависимости минимальных частот колебаний ю от температуры T при разных значениях показателя объемной доли N для ФГ-оболочек с различными граничными условиями; Конфигурация I: а - CC; б - CS
Наличие жидкости внутри оболочки оказывает наиболее сильное влияние в случае консольного закрепления. Если для жестко закрепленных оболочек нагрев внешней поверхности приводит к возрастанию критических температур, как для пустых оболочек, так и оболочек с жидкостью, то для консольных оболочек нагрев на разных поверхностях по-разному влияет на критические температуры пустых и наполненных жидкостью оболочек. Для более наглядной демонстрации этих различий на рис. 5 приведены зависимости минимальных частот колебаний ю (Гц) от температуры T (К), полученные для разных конфигураций пустых и содержащих жидкость консольных оболочек. Из представленных данных следует, что в случае пустых оболочек нагрев внешней поверхности приводит к снижению критических температур потери устойчивости. Однако для оболочек, содержащих жидкость, имеет место, наоборот, стабилизирующее воздействие. В качестве возможного объяснения этого явления могут быть приведены следующие рассуждения. Выше было отмечено, что изменение поверхности нагрева оказывает влияние на величину
31
Bochkarev S.A., Lekomtsev S.V. /PNRPUMechanics Bulletin 4 (2015) 19-35
нормальных перемещений свободного края оболочки. С другой стороны, значение гидродинамического давления, действующего по нормали к упругой поверхности, имеет прямо пропорциональную зависимость от величины нормальной составляющей вектора перемещений оболочки. Следовательно, раскрытие свободного края оболочки, отрицательно влияя на ее температурную устойчивость, тем не менее способствует уменьшению гидродинамического давления, оказывая стабилизирующее воздействие.
Таблица 5
Критические значения температуры T (K) в зависимости от объемной доли N для оболочек с различными граничными условиями, содержащих жидкость
ГУ R/h N = 0 N = 0,3 N = 1 N = 5 N = 10 N = 30 N = да
Конфигурация I
CC 200 702,01 786,31 899,97 997,97 1000,4 994,35 987,88
500 495,53 532,22 580,15 609,19 604,72 598,24 594,11
CF 200 1031,4 1238,9 1402,6 1414,1 1398,4 1392,9 1397,5
500 671,83 814,73 1013,9 1050,5 1034,7 1032,6 1040,8
Конфигурация II
CC 200 727,91 762,54 803,02 919,38 965,93 1008,5 1033,2
500 501,15 514,62 525,58 562,04 581,04 602,35 617,29
CF 200 1637,9 1648,6 1655,3 1735,8 1765,2 1781,1 1784,8
500 1055,8 1072,0 1090,9 1221,2 1296,7 1365,1 1400,2
Конфигурация III
CC 200 987,43 814,19 720,15 689,73 693,98 699,32 701,91
500 593,85 527,12 494,77 488,42 491,48 494,31 495,48
CF 200 1397,9 1193,5 1044,3 1022,7 1034,9 1037,8 1031,8
500 1041,5 774,36 675,65 668,65 674,66 675,19 672,01
Конфигурация IV
CC 200 1034,1 972,26 890,21 790,47 764,74 742,03 728,40
500 617,89 596,22 564,94 528,36 517,81 507,72 501,38
CF 200 1784,8 1790,5 1789,9 1739,6 1707,7 1669,9 1639,1
500 1401,3 1401,0 1388,6 1277,9 1203,8 1117,6 1057,9
Рис. 5. Зависимости минимальных частот колебаний ю от температуры T для консольной ФГ-оболочки при разных конфигурациях: N = 1; арабские цифры - пустая оболочка; римские цифры - оболочка с жидкостью; сплошные линии - нагрев снизу; пунктирные линии - нагрев
сверху
32
Бочкарев С.А., Лекомцев С.В. /Вестник ПНИПУ. Механика 4 (2015) 19-35
В табл. 5 также отражены критические температуры, полученные при других геометрических размерах (h = 4х10-4 м, R/h = 500). В этом случае качественных изменений не наблюдается. Таким образом, для некоторых конфигураций нетипичная зависимость от температуры сохраняется и для более тонких оболочек.
Работа выполнена при финансовой поддержке РФФИ (грант № 13-01-96049).
Библиографический список
1. Birman V., Byrd L.W. Modeling and analysis of functionally graded materials and structures // Appl. Mech. Rev. - 2007. - Vol. 60. - No. 5. - P. 195-216.
2. Haddadpour H., Mahmoudkhani S., Navazi H.M. Supersonic flutter prediction of functionally graded cylindrical shells // Compos. Struct. - 2008. - Vol. 83. - No. 4. - P. 391-398.
3. Mahmoudkhani S., Haddadpour H., Navazi H.M. Supersonic flutter prediction of functionally graded conical shells // Compos. Struct. - 2010. - Vol. 92. - No. 2. - P. 377-386.
4. Sabri F., Lakis A.A. Efficient hybrid finite element method for flutter prediction of functionally graded cylindrical shells // J. Vib. Acoust. - 2014. - Vol. 136. - No. 1. - 011002.
5. Бочкарев C.A., Лекомцев С.В Исследование панельного флаттера круговых цилиндрических оболочек, выполненных из функционально-градиентного материала // Вестник ПНИПУ. Механика. - 2014. - № 1. - С. 57-75.
6. Бочкарев C.A., Лекомцев С.В. Анализ панельного флаттера нагруженных круговых цилиндрических оболочек, выполненных из функционально-градиентного материала // Неравновесные процессы в соплах и струях: материалы X Междунар. конф. - М: Изд-во МАИ, 2014. - C. 335-337.
7. Chen W.Q., Bian Z.G., Ding H.J. Three-dimensional vibration analysis of fluid-filled orthotropic FGM cylindrical shells // Int. J. Mech. Sci. - 2004. - Vol. 46. - No. 1. - P. 159-171.
8. Vibration characteristics of FGM circular cylindrical shells filled with fluid using wave propagation approach / Z. Iqbal [et al.] // Appl. Math. Mech. - 2009. - Vol. 30. - No. 11. - P. 1393-1404.
9. Vibrational study of fluid-filled functionally graded cylindrical shells resting on elastic foundations / A.G. Shah [et al.] // ISRN Mechanical Engineering. - 2011. - Vol. 2011. - 892460.
10. Бочкарев C.A., Лекомцев С.В. Собственные колебания цилиндрических оболочек, содержащих неподвижную жидкость и выполненных из функционально-градиентных материалов // Современные проблемы математики, механики, информатики: материалы междунар. науч. конф. -Тула: Изд-во Тул. гос. ун-та, 2013. - С. 182-189.
11. Sheng G.G., Wang X. Thermomechanical vibration analysis of a functionally graded shell with flowing fluid // Eur. J. Mech. A-Solid. - 2008. - Vol. 27. - No. 6. - P. 1075-1087.
12. Sheng G.G., Wang X. Dynamic characteristics of fluid-conveying functionally graded cylindrical shells under mechanical and thermal loads // Compos. Struct. - 2010. - Vol. 93. - No. 1. - P. 162-170.
13. Бочкарев С.А., Лекомцев С.В. Гидроупругая устойчивость функционально-градиентных цилиндрических оболочек, содержащих жидкость // Современные проблемы математики, механики, информатики: материалы междунар. науч. конф. - Тула: Изд-во Тул. гос. ун-та, 2014. - С. 133140.
14. Kadoli R., Ganesan N. Buckling and free vibration analysis of functionally graded cylindrical shells subjected to a temperature-specified boundary condition // J. Sound Vib. - 2006. - Vol. 289. -No. 3. - P. 450-480.
15. Bhangale R.K., Ganesan N., Padmanabhan C. Linear thermoelastic buckling and free vibration behavior of functionally graded truncated conical shells // J. Sound Vib. - 2006. - Vol. 292. - No. 1-2. -P. 341-371.
16. Malekzadeh P., Fiouz A.R., Sobhrouyan M. Three-dimensional free vibration of functionally graded truncated conical shells subjected to thermal environment // Int. J. Pres. Ves. Pip. - 2012. -Vol. 89. - P. 210-221.
33
Bochkarev S.A., Lekomtsev S.V. /PNRPUMechanics Bulletin 4 (2015) 19-35
17. Akbari M., Kiani Y., Eslami M.R. Thermal buckling of temperature-dependent FGM conical shells with arbitrary edge supports // Acta Mech. - 2015. - Vol. 226. - No. 3. - P. 897-915.
18. Бочкарев C.A., Матвеенко В.П. Численное исследование влияния граничных условий на динамику поведения цилиндрической оболочки с протекающей жидкостью // Изв. РАН. МТТ. -2008. - № 3. - C. 189-199.
19. Бочкарев C.A. Собственные колебания вращающейся круговой цилиндрической оболочки с жидкостью // Вычислительная механика сплошных сред. - 2010. - № 2. - С. 24-33.
20. Вольмир А. С. Оболочки в потоке жидкости и газа. Задачи гидроупругости. - М.: Наука, 1979. - 320 c.
21. Алфутов Н.А., Зиновьев П.А., Попов В.Г. Расчет многослойных пластин и оболочек из композиционных материалов. - М.: Машиностроение, 1984. - 264 с.
22. Reddy J.N., Chin C.D. Thermomechanical analysis of functionally graded cylinders and plates // J. Therm. Stresses. - 1998. - Vol. 21. - No. 6. - P. 593-626.
23. Бочкарев С.А., Матвеенко В.П. Численное моделирование устойчивости нагруженных оболочек вращения при внутреннем течении жидкости // ПМТФ. - 2008. - № 2. - С. 185-195.
24. Зенкевич О. Метод конечных элементов в технике. - М.: Мир, 1975. - 544 с.
25. Матвеенко В.П. Об одном алгоритме решения задачи о собственных колебаниях упругих тел методом конечных элементов // Краевые задачи теории упругости и вязкоупругости. - Свердловск, 1980. - С. 20-24.
26. Praveen G., Chin C., Reddy J. Thermoelastic analysis of functionally graded ceramic-metal cylinder // J. Eng. Mech. - 1999. - Vol. 125. - No. 11. - P. 1259-1267.
References
1. Birman V., Byrd L.W. Modeling and analysis of functionally graded materials and structures. Appl. Mech. Rev., 2007, vol. 60, no. 5, pp. 195-216. DOI: 10.1115/1.2777164
2. Haddadpour H., Mahmoudkhani S., Navazi H.M. Supersonic flutter prediction of functionally graded cylindrical shells. Compos. Struct., 2008, vol. 83, no. 4, pp. 391-398. DOI: 10.1016/j.compstruct.2007.05.011
3. Mahmoudkhani S., Haddadpour H., Navazi H.M. Supersonic flutter prediction of functionally graded conical shells. Compos. Struct., 2010, vol. 92, no. 2, pp. 377-386. DOI: 10.1016/j.compstruct.2009.08.018
4. Sabri F., Lakis A.A. Efficient hybrid finite element method for flutter prediction of functionally graded cylindrical shells. J. Vib. Acoust., 2014, vol. 136, no. 1, 011002. DOI: 10.1115/1.4025397
5. Bochkarev S.A., Lekomtsev S.V. Investigation of panel flutter of functionally graded circular cylindrical shells. PNRPU Mechanics Bulletin, 2014, no. 1, pp. 57-75.
6. Bochkarev S.A., Lekomtsev S.V. Analiz panel'nogo flattera nagruzhennykh krugovykh tsilin-dricheskikh obolochek, vypolnennykh iz funktsional'no-gradientnogo materiala [Analysis of panel flutter of loaded circular cylindrical shells made of functionally graded material]. Materialy X Mezhdunarodnoi konferentsii “Neravnovesnyeprotsessy v soplakh i struiakh”. Moscow, 2014, pp. 335-337.
7. Chen W.Q., Bian Z.G., Ding H.J. Three-dimensional vibration analysis of fluid-filled orthotropic FGM cylindrical shells. Int. J. Mech. Sci., 2004, vol. 46, no. 1, pp. 159-171. DOI: 10.1016/j.ijmecsci.2003.12.005
8. Iqbal Z., Naeem M.N., Sultana N., Arshad S.H., Shah A.G. Vibration characteristics of FGM circular cylindrical shells filled with fluid using wave propagation approach. Appl. Math. Mech, 2009, vol. 30, no. 11, pp. 1393-1404. DOI: 10.1007/s10483-009-1105-x
9. 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, vol. 2011, 892460. DOI: 10.5402/2011/892460
10. Bochkarev S.A., Lekomtsev S.V. Sobstvennye kolebaniia tsilindricheskikh obolochek, soderzhashchikh nepodvizhnuiu zhidkost' i vypolnennykh iz funktsional'no-gradientnykh materialov
34
Бочкарев С.А., Лекомцев С.В. /Вестник ПНИПУ. Механика 4 (2015) 19-35
[Natural vibrations of cylindrical shells containing quiescent fluid and made of functionally graded materials]. Materialy mezhdunarodnoi nauchnoi konferentsii “Sovremennye problemy matematiki, mekhaniki, informatiki”. Tulskii gosudarstvennyi universitet, 2013, pp. 182-189.
11. Sheng G.G., Wang X. Thermomechanical vibration analysis of a functionally graded shell with flowing fluid. Eur. J. Mech. A-Solid, 2008, vol. 27, no. 6, pp. 1075-1087. DOI: 10.1016/j.euromechsol.2008.02.003
12. Sheng G.G., Wang X. Dynamic characteristics of fluid-conveying functionally graded cylindrical shells under mechanical and thermal loads. Compos. Struct, 2010, vol. 93, no. 1, pp. 162-170. DOI: 10.1016/j.compstruct.2010.06.004
13. Bochkarev S.A., Lekomtsev S.V. Gidrouprugaia ustoichivost' funktsional'no-gradientnykh tsilindricheskikh obolochek, soderzhashchikh zhidkost' [Hydroelasticity stability of functionally graded cylindrical shells containing fluid]. Materialy mezhdunarodnoi nauchnoi konferentsii “Sovremennye problemy matematiki, mekhaniki, informatiki”. Tulskii gosudarstvennyi universitet, 2014, pp. 133-140.
14. Kadoli R., Ganesan N. Buckling and free vibration analysis of functionally graded cylindrical shells subjected to a temperature-specified boundary condition. J. Sound Vib., 2006, vol. 289, no. 3, pp. 450-480. DOI: 10.1016/j.jsv.2005.02.034
15. Bhangale R.K., Ganesan N., Padmanabhan C. Linear thermoelastic buckling and free vibration behavior of functionally graded truncated conical shells. J. Sound Vib., 2006, vol. 292, no. 1-2, pp. 341371. DOI: 10.1016/j.jsv.2005.07.039
16. Malekzadeh P., Fiouz A.R., Sobhrouyan M. Three-dimensional free vibration of functionally graded truncated conical shells subjected to thermal environment. Int. J. Pres. Ves. Pip., 2012, vol. 89, pp. 210-221. DOI: 10.1016/j.ijpvp.2011.11.005
17. Akbari M., Kiani Y., Eslami M.R. Thermal buckling of temperature-dependent FGM conical shells with arbitrary edge supports. Acta Mech, 2015, vol. 226, no. 3, pp. 897-915. DOI: 10.1007/s00707-014-1168-3
18. Bochkarev S.A., Matveenko V.P. Numerical study of the influence of boundary conditions on the dynamic behavior of a cylindrical shell conveying a fluid. Mech. Solids, 2008, vol. 43, no. 3, pp. 477486. DOI: 10.3103/S0025654408030187
19. Bochkarev S.A. Sobstvennye kolebaniia vrashchaiushcheisia krugovoi tsilindricheskoi obolochki s zhidkost'iu [Natural vibrations of a rotating circular cylindrical shell containing fluid]. Computational continuum mechanics, 2010, vol. 3, no. 2, pp. 24-33. DOI: 10.7242/1999-6691/2010.3.2.14
20. Volmir A.S. Shells in liquid and gas flow: Hydroelastisity problems. Moscow: Nauka, 1979. 320 p.
21. Alfutov N.A., Zinov’ev P.A., Popov B.G. Raschet mnogosloinykh plastin i obolochek iz kompozitsionnykh materialov [Analysis of multilayer plates and shells of composite materials]. Moscow: Mashinostroenie, 1984. 264 p.
22. Reddy J.N., Chin C.D. Thermomechanical analysis of functionally graded cylinders and plates. J. Therm. Stresses, 1998, vol. 21, no. 6, pp. 593-626. DOI: 10.1080/01495739808956165
23. Bochkarev S.A., Matveenko V.P. Numerical modelling of the stability of loaded shells of revolution containing fluid flows. J. Appl. Mech. Tech. Phys., 2008, vol. 49, no. 2, pp. 313-322. DOI: 10.1007/s10808-008-0043-1
24. Zienkiewicz O.C. The finite element method in engineering science. London: McGraw Hill, 1971. 521 p.
25. Matveenko V.P. Ob odnom algoritme resheniia zadachi o sobstvennykh kolebaniiakh uprugikh tel metodom konechnykh elementov [On an algorithm of solving the problem on natural vibrations of elastic bodies by the finite element method]. Kraevye zadachi teorii uprugosti i viazkouprugosti. Sverdlovsk, 1980, pp. 20-24.
26. Praveen G., Chin C., Reddy J. Thermoelastic analysis of functionally graded ceramic-metal cylinder. J. Eng. Mech., 1999, vol. 125, no. 11, pp. 1259-1267. DOI: 10.1061/(ASCE)0733-9399(1999)125:11(1259)
35
