МЕХАНИКА
Научная статья УДК 539.31
DOI: 10.34759/trd-2021-120-09
Нестационарное деформирование анизотропной круговой
цилиндрической оболочки
Наталья Александровна Локтева1, Дмитрий Олегович Сердюк213, Павел Дмитриевич Скопинцев3, Григорий Валерьевич Федотенков4
1,2,3,4Московский авиационный институт (национальный исследовательский университет), Москва, Россия
1,4НИИ механики МГУ имени М.В. Ломоносова, Москва, Россия
1nlok@rambler.ru
2d.serduk55@gmail.com3
3chgpashka@gmail.com
4greghome @mail. ru
Аннотация. Исследуется нестационарное деформирование тонкой неограниченной по длине круговой цилиндрической оболочки постоянной толщины при воздействии на ее боковую поверхность сосредоточенной и распределённой по произвольной области нагрузки с переменной во времени амплитудой. Материал оболочки полагается линейно упругим, анизотропным и обладающим симметрией относительно ее срединной поверхности. Движение оболочки рассматривается в цилиндрической системе координат, связанной с осью цилиндрической оболочки, и описывается при помощи гипотез Киргофа-Лява, а искомая функция нормального нестационарного перемещения строится на связи функции Грина с функцией
действующей нагрузки при помощи интегрального оператора типа свертки по пространственным переменным и времени. Функция Грина для анизотропной оболочки представляет собой решение специальной задачи о воздействии на оболочку сосредоточенной нагрузки, математически моделируемой дельта-функцией Дирака. Для построения функции Грина применяются разложения в экспоненциальные ряды Фурье, интегральное преобразование Лапласа по времени и интегральное преобразование Фурье по продольной координате. Обратное интегральное преобразование Лапласа выполняется аналитически, а оригинал интегрального преобразования Фурье находится с использованием численных методов интегрирования быстро осциллирующих функций. Интегралы сверки функции Грина с функцией нагрузки берутся при помощи квадратурных формул методом прямоугольников. В качестве численных примеров рассмотрены распространения нестационарных возмущений в неограниченной оболочке для нескольких вариантов симметрии упругой среды.
Ключевые слова: анизотропная цилиндрическая оболочка, нестационарная динамика, функция Грина, функция прогиба, обобщенные функции, интегральные преобразования, квадратурные формулы, нормальное перемещение, оболочка Кирхгофа-Лява
Финансирование: работа выполнена при финансовой поддержке РФФИ (проекты 2058-00023 и 19-08-00968)
Труды МАИ. 2021.Выпуск № 120 Trudy MAI. 2021.Issues no.120
Для цитирования: Локтева Н.А., Сердюк Д.О., Скопинцев П.Д., Федотенков Г.В. Нестационарное деформирование анизотропной круговой цилиндрической оболочки // Труды МАИ. 2021. № 120. DOI: 10.34759/trd-2021-120-09
MECHANICS
Original article
Unsteady deformation of anisotropic circular cylindrical shell
Natalia A. Lokteva1, Dmitry O. Serdyuk213, Pavel D. Skopincev3, Grigory V. Fedotenkov4
1,2,3,4Moscow Aviation Institute (National Research University), Moscow, Russia
1,4Lomonosov Research Institute of Mechanics, Moscow, Russia
1nlok@rambler.ru
2d.serduk55@gmail.comB
3chgpashka@gmail.com
4greghome @mail. ru
Abstract.The unsteady deformation of a thin, infinitely long circular cylindrical shell of constant thickness under the impact of a concentrated shock load distributed over an arbitrary region on its lateral surface is being studied. The shell material is assumed to be linearly elastic, anisotropic, and symmetrical about its middle surface. The Kirchhoff-Love model is employed to describe the shell motion. The motion of the shell is being considered in a cylindrical coordinate system associated with the axis of the cylindrical shell, and the sought-for function of the normal non-stationary deflection is constructed by connecting the
Труды МАИ. 2021.Выпуск № 120 Trudy MAI. 2021.Issues no.120
Green's function with the function of the operating load using an integral operator of the convolution type in spatial variables and time. The Green's function for an anisotropic shell is a solution to a special problem of the impact of an instantaneous concentrated load on the shell, mathematically modeled by the Dirac delta functions. Expansions in exponential Fourier series, integral Laplace transform in time and integral Fourier transform in longitudinal coordinate are being used to construct the Green's function. The inverse integral Laplace transform is being performed analytically, and the original integral Fourier transform is being found using numerical methods for integrating rapidly oscillating functions. The integrals of the convolution of the Gree n's function with the load function are being taken with quadrature formulas using the rectangle method. As an example, the unsteady dynamics of cylindrical shell was sstudied under the impact of arbitrarily time-dependent concentrated load and load distributed over the finite area belonging to the lateral surface of the shell. Several options of symmetry of the elastic medium (isotropic, orthotropic and anisotropic) herewith were analyzed, which demonstrates calculated solution versatility both in terms of influence nature and shell material. For the considered symmetry options, the study of the unsteady vibrations propagation character, which allowed evaluating the solution adequacy, was conducted.The presented approach to constructing the unsteady deflection function while transition to the dimensional values opens opportunities for the analysis of the stress-strain state of the extended cylindrical shells with account for various options of the material anisotropy and law of distribution of the unsteady loading along both coordinates and time.
Труды МАИ. 2021.Выпуск № 120 Trudy MAI. 2021.Issues no.120
Keywords: anisotropic cylindrical shell, unsteady dynamics, Green's function, deflection function, generalized functions, integral transformations, quadrature formulas, normal displacement, Kirchhoff-Love shell
Funding: the work was supported by RFBR projects 20-58-00023 and 19-08-00968 For citation: Lokteva N.A., Serdyuk D.O., Skopincev P.D., Fedotenkov G.V. Unsteady deformation of anisotropic circular cylindrical shell. Trudy MAI, 2021, no. 120. DOI: 10.34759/trd-2021-120-09
Введение
Оболочки применяются во многих отраслях промышленности, включая
ракетно-космическую, самолетостроительную, атомную, судостроительную,
нефтегазодобывающую и машиностроительную отрасли. Для усовершенствования и
создания новых перспективных конструкций, работающих в нестационарных
режимах, необходимо обладать знаниями о протекании процессов нестационарных
вынужденных колебаний в проектируемых оболочечных конструкциях.
Исследование нестационарных колебаний является важной задачей, так как при этом
существенно проявляется неоднородность искомого решения по времени и
координатам. Напряженно-деформированное состояние цилиндрических оболочек
при воздействии ударных нагрузок, имитируемых импульсными функциями,
представляет теоретический и прикладной интерес.
На данный момент, в рамках модели Кирхгофа-Лява, наиболее полно изучены
задачи по исследованию нестационарной динамики изотропных упругих пластин и
5
Труды МАИ. 2021.Выпуск № 120 Trudy MAI. 2021.Issues no.120
оболочек [1]. В работах [2,3] рассмотрен широкий круг проблем динамики ортотропных цилиндрических оболочек, их осесимметричное и неосесимметричное деформирование при продольном ударе. Так же рассмотрено неосесимметричное деформирование при нестационарном внешнем давлении. В труде [4] отражены проблемы деформирования подкрепленных цилиндрических оболочек при действии динамических сжимающих осевых нагрузок и внешнего давления. В работе [5] по изучению динамики оболочек автор рассматривал тонкостенные однородные цилиндрические оболочки из композиционных материалов при импульсном воздействии. За модель композиционного материала принималась ограниченная слоистая среда, составленная из ортотропных упругих слоев неизменной толщины. Нагрузкой являлся импульс внешнего или внутреннего давления. В работах [6-12] рассмотрено применение метода функций Грина для решения нестационарных задач теории упругости и теории оболочек. Исследованию нестационарных контактных задач для тонких сферических, цилиндрических оболочек и упругого полупространства посвящены труды [6-9]. Рассмотрены задачи для ограниченной балки Тимошенко под воздействием нестационарной нагрузки, проанализированы вопросы идентификации дефектов в упругих стержнях [10,11]. Исследовались случаи нестационарного воздействия жесткого индентора на упругую полуплоскость [12]. Исследованию проблем, связанных с динамикой, прочностью и устойчивостью оболочек, посвящены работы [13-15]. Исследования распространения нестационарных возмущениях в пластинах, в том числе и анизотропных, посвящены работы [16-28].
Труды МАИ. 2021.Выпуск № 120 Trudy MAI. 2021.Issues no.120
В свою очередь, задачи нестационарной динамики упругих анизотропных
оболочек остаются недостаточно хорошо изученными. В данной работе аналитически
исследуется процесс воздействия нестационарной нагрузки на упругую
анизотропную цилиндрическую оболочку. Строятся нестационарные функции
нормальных перемещений при действии сосредоточенной и распределенной нагрузки
для анизотропной оболочки, которые открывают возможности для обширных
прикладных и научных исследований. Найденные нестационарные функции
являются универсальными по отношению к свойствам материала оболочки, который
может быть изотропным, трансверсально-изотропным, ортотропным или
анизотропным. Представлены примеры расчетов.
1. Постановка задачи
Объектом исследования является тонкая неограниченная цилиндрическая оболочка (см. рис. 1). Оболочка имеет радиус R и толщину h. Материал цилиндрической оболочки принят упругим и анизотропным. При этом рассматривается случай анизотропии, при котором упругая среда имеет одну поверхность симметрии. В данном случае такой поверхностью является срединная поверхность оболочки. Для тонкой оболочки Кирхгофа-Лява рассматриваемый
материал имеет шесть независимых упругих постоянных:
_ 122 _ ^1112 _ ^2222 _f1222 _ ^1212
C11 ~ C , C12 ~ C , C16 ~ C , C22 ~ C , C26 ~ C , C66 ~ C ■
/Ma,z,11
z
Рис. 1. Цилиндрическая оболочка под влиянием нестационарной нагрузки В начальный момент времени оболочка находится в невозмущенном состоянии. Затем на нее воздействует нестационарная нагрузка Р(а,г,г) , распределенная по
произвольной области В. Движение оболочки рассматривается в цилиндрической системе координат ОЯаг.
Постановка задачи включает в себя уравнения движения упругой оболочки Кирхгофа-Лява, соответствующие геометрические и физические соотношения с учетом симметрии свойств материала исследуемой оболочки [1], [29], [30].
Уравнения движения в перемещениях анизотропной оболочки Кирхгофа-Лява имеют вид [1]:
где дифференциальные операторы определяются следующими соотношениями:
(1)
R2 да2 2 i
а + 2hci6 д2ua ,
R dadz
v
hc +Ic66
' 6 ^ R 2
2
д u
dz
2
K12(Uz )
hcl6 д u
R2 да2
+
hcn hc^ Ic,
+ ■
66
66
^ д 2u , д2u„
v R R R3 j
дaдz
■ + hc.
дz2
K0( w)
hcn ôw 2Ic66 д w
Ic6l д w
R2 да R2 дадz2 R3 дa2дz
+
hc61 Ic,
6i
R R3
дw Ic62 д w дz R дz3 '
K21(иа ) = K12(4x X
K 22 (Uz )
Ic66 hc
66
R R
д2u„ 2hc26 д2u
да2
+
R дaдz
+ hc
Ô\ 22 дz2
(3)
К23 ( w) =
I^ д3w 2Ic« Ô3W
+
66
, Ic62
R4 да3 R' да'ôz R2 дaдz
2
R дz
hc 61 Ic,
6i
R2
R4
ôw да
K31(Ua ) = -K13(4x X K32(u, ) = K 23 (uz X
K3 ( w) = - 4Ic62 д4 w
___4Icl6 д4w
R дaдz3 R3 да^
2Ic21 4Ic,
R2
+
66
R2
д 4w I^ д4 w
дa2дz2 R4 да4
j
Ic
д4 w 2Ic21 д2 w 4Ic16 д2 w 21^ д2 w
22 4
ôz 4
R2 ôz2
R3 дaдz R4 да2
rhÇiL + Iç^ R2 R4
w.
j
(4)
В (1) - (4): р- плотность материала оболочки, R - радиус оболочки, И - толщина оболочки, w - нормальное перемещение, иа, иг - компоненты вектора тангенциальных перемещений, Р - нормальная нагрузка, дг - тангенциальная нагрузка, I = И3/12 - погонный момент инерции.
Уравнения (1) - (4) совместно с начальными условиями
а It=0
дг
ôu„
= u_
z lt=0
t=0
дг
t=0
I ôw
= w^=0 = ~et
= 0,
(5)
t=0
образуют начальную задачу.
Труды МАИ. 2021.Выпуск № 120 Trudy MAI. 2021.Issues no.120
Цель исследования заключается в построении нестационарной функции
нормального перемещения w (a, z, t ) в ответ на воздействие нестационарной нагрузки P (a, z, t ).
2. Построение нестационарной функции нормального перемещения
Введем систему безразмерных величин и выполним обезразмеривание
соотношений (1) - (5):
г иа , иг , w , г , Я с.г иа = —, и2 = —, w = —, г = —, к = —, х =—,
а ь г ь ь ь ь ь
я;== М ч:=,У = рк,с. = Ь,ь=„, (6)
пси пси пси \ р
_ (-< ' _ (-< ' _ с22 (-< ' _ с26 (-< ' _ с66 с11 с11 с11 с11 с11
В (6): X - безразмерное время, ь - характерный размер, с - характерная скорость, к
- коэффициент отношения радиуса оболочки к ее толщине.
Уравнение движения и дифференциальные операторы в безразмерной форме (штрихи в безразмерных величинах опущены) примут вид:
д 2и
Эх2 Эй Эх Э2 w Эх2
= K1l(«a ) + K12^z ) + K13(W) + Ча,
2 = K2i(Ua ) + K22(uz ) + ^2з( w) + qz, (7)
= K31(Ua ) + K32(Uz ) + K33(W) + P-
^ / л 1 Э 2u
KH(Ua ) = F
/ N С Э 2u
K12("z) = С ôU +
^ , N 1 ôw С Э3w
K13(w) = TT^--./2 Л л 2
a + 2С2 ÔUa +
k ôaôz
С+
С
12k2 y
ô 2u
Cl + С5
v k
С
k 12k3
ô U
ôaôz
+С
С ô3w
k2 ôa 6k2 ôaôz2 12k3 ôa 2ôz
+
ôz2 '
ôU
ôz2 ôw
k 12k3
С ô3w
ôz 12k ôz
K 21(Ua) = K12(4aX
K22(Uz ) =
С С л
С5 + С 12k4 k2
ôU 2С ôU
ôa2
+
k ôaôz
+ С
3 ôz2
K23 ( w) =
С ô3w 12k24 ôa3
+
v
С2 + С2 k2 12k4
л
ôw С ô3w
+
+
С
ô w С ôw ■ + ■ 1
ôa 6k ôa ôz 12k ôaôz k ôz
(9)
K31(Ua) = - K13(UaX K32(Uz) = - K23(UzX
^ С4 ô 4 w С2 ô4 w
K33(w) = - 4 2
^ ^ л ô4w
3k ôaôz3 3k3 ôa3ôz
С С 6k2 3k2
ôa2ôz2 12k4 ôa
1 ô4 w
44
(10)
С ô4w С ô2w С ô2w 1 ô2w
12 ôz4 6k2 ôz2 3kj ôaôz 6k4 ôa
42
1
+
1
k2 12k4
w.
Начальные условия примут вид:
du
a 1х=0
ôx
ôu.
= u_
z lx=0
x=0
ôx
x=0
I ôw
= w x=0 = ôX
= 0.
x=0
(11)
Решение начальной задачи (7) - (11) может быть построено с помощью функции Грина ^ (а, г, х) для нормального перемещения [1]:
w(a, z, х) = G (a, z, x)*** p (a, z, x)
(12)
В (12) через * обозначены свертки по пространственным координатам а, г и
безразмерному времени х.
Труды МАИ. 2021.Выпуск № 120 Trudy MAI. 2021.Issues no.120
Обозначим за G (а,z,х) и G (а,z,х) функции Грина для тангенциальных
V ' uz \ '
перемещений. Тангенциальные давления qi в (7) примем равными нулю. Таким образом, функции Грина Gw (а, z,х), G„ (а,z,х) и GM (а,z,х) являются решениями следующей задачи:
д 2G.,
дх2 д 2G„
= KiiG ) + Ki2(Guz ) + Ki3(Gw ), 2 = K2G ) + K22G ) + K23(Gw),
дх
2
= Kзl(Gua ) + ) + K33G ) + 8(а, z )8(х).
G
= д- G
х=0 дх u
G
х=0
=д_ G
х=0 дх u
х=0
G I =— G
w U дх w
0,
(13)
х=0
где 8(*) - дельта-функция Дирака, а дифференциальные операторы K (G ) имеют
вид (8) - (10), где необходимо заменить соответствующие искомые функции на функции влияния.
Для решения поставленной начальной задачи применим к (13) разложения в экспоненциальные ряды Фурье по углу а , а также интегральные преобразования Фурье по координате г и Лапласа по времени т . В результате получим систему алгебраических уравнение относительно изображений & , & п, &п функций 1рина
в пространстве преобразований Фурье и Лапласа в коэффициентах рядов:
( *2+Qi ) gLF+Q2GLF+Q3GLF = 0
Q2GLF+( *2+Q4 ) GLF+Q5GLnF = 0 -Q3GF - QGLI +( *2+Q6 ) GLF = ±
(14)
Труды МАИ. 2021.Выпуск № 120 Trudy MAI. 2021.Issues no.120 где
Qi =
/
С
л
С5 +
V 5 12k у
q +
2пС2 ~k~
n
■q +
Q2 = C4q2 + С
пС5 + nC nC5
v k 12k3 у
q+
n C0
Q3=- с q
inC,
12k 6k
Q4 = C3q2 +
25 q2
n C2 C2 C2
V
2 2пС,
12k 3
2
k
k2 Л
i
у
iq
in
k
2
Q5 =
in С
42
12k2
q
k
С n2C5
. n С n С 4 q +—4 + 5
12k4 k
2 '
v k 6k у
тС2 тС2 - 1п3С2
iq
k2
12k4
„ С 4 пС 3
Q6 = 77: q +-;г q +
12
3k
n С - Q n С5
V
6k2
+
п С - пС2
3k3
n
3k2 n4 +1 1
q2+
ч у 6к4 12к4 к2' Здесь и далее верхний индекс « £ » у функции означает её преобразование по Лапласу, а « ^ » - её преобразование по Фурье, я - параметр преобразования Лапласа, ^ - параметр преобразования Фурье, п - номер коэффициента ряда.
Решая систему (14), получим изображения функций Грина для прогиба и тангенциальных перемещений:
1
GLZ ( q, * ) GZ ( q, * )=i С ( q, * )
*4 + *2 R + R
2л *6 + *4 R + *2 R + R '
1 - * Q + R
2л *6 + *4 R + *2 R + R '
1 * 2Q5 + R
(15)
2л *6 + *4R + *2R + R
Ri = Qi + Q4, R = Q1Q4 - Q22, R3 = Qi + Q4 + Q6,
R4 = 6164 + 6166 - а2 + 6э2 + QQ + 652, R5 = 616466 + - 6бб22 - 2626365 + 646Л R6 = 6265 - 6364, R = 6165 - 6263-
Оригинал коэффициентов (15) по Лапласу найден аналитически при помощи
таблиц [31] и разложения на простые дроби:
Gl (q,х) =i( Ash(*lх) + Bsh(*2х) + Сsh(*3х)), к
(16)
где:
л= *i4 + Ri*i2 + R2 *2 + Ri*22 + R2
4 , d „ 2
2*i (*i2 - *32 )(*i2 - *22 ) ' 2*2 (*i2 - *22 )(*22 - *32 )
С =
*34 + Ri*32 + R2 2*3 (*22 - *32 )(*i2 - *32 )
и * - корни знаменателя уравнений в (15):
*2 =
^6U1/3 (-2 RU1/3 + J ) pUi3(J2 -4R3Ui3 - J)
6U1/3 '
^-3U1/3 ( J2 + 4 RU1/3 + J )
*3 = U3 ■
J = U2/3 + 4R2 - i2R,
J = /V3u2/3 - 4/V3R2 + i2/V3R4
U = 36RR - i08R - 8R3 + i^i2RR3 - 3R2R32 - 54RRR + i2R43 + 8iR2
Труды МАИ. 2021.Выпуск № 120 Trudy MAI. 2021.Issues no.120
Для нахождения оригинала (16) по Фурье применим численный метод
интегрирования быстро осциллирующих функций [32]. Тогда оригиналы функции
Грина в коэффициентах рядов примет вид:
1 N- А
Gwn (т) = Т-2 XT
2л к=0 2
■í^z+ítz
Dfil ( qt, т)+АС ( Яы, т)
(17)
2A А ^ sin m , m cos m - sin m.
А = ^ m = ^ Di,2 =-f-2-^
N 2 mm
q = - A + к А, q+i = - A + (k + 1)А, к = 0..N -1.
где A, N - параметры численного интегрирования.
С учетом разложения в экспоненциальные ряды и соотношения (17) нестационарная функция Грина для нормального прогиба анизотропной неограниченной цилиндрической оболочки Кирхгофа-Лява примет вид:
i N-1 œ д
G, (a,z, т) = -L
2 k=0 n=-œ 2
( DGl ( ík-т)+АС ( т))
(18)
В случае воздействия внешнюю поверхность оболочки внезапно приложенной сосредоточенной нагрузки, изменяющейся во времени по закону Р (т), выражение
для нагрузки р(а,г,т) из (12) запишется так:
p (a, z, т) = P (т) H (т)5(а)5( z ),
(19)
где Н (т) - функция Хевисайда.
Тогда соотношение (12) с учетом (19) и свойств дельта функции Дирака преобразуется в следующий вид:
w(a,z,x) = Gw (a,z,x) * * * p(a,z,x) = JGw (a,z,x -1)P(t)dt.
(20)
Для вычисления интеграла в (20) используется метод прямоугольников [32]. Тогда приближенное выражение для искомой функции нестационарного прогиба примет вид:
^ X г
w(a, z, х)« G
i=1 n
V
Xi
a,z,x--
n
P
у
V n у
(21)
В случае воздействия на оболочку распределенной по области боковой поверхности нагрузки, выражение p (а, г, т) из (12) примет вид:
p (a, z, х) = P (х) H (х)
H Г Pi a + - - H f a - РЛ
V 2 у V 2 J
x
X
Г ь Л Г ь Л
H z + - -H z--
V 2 у V 2 у
(22)
что соответствует приложению к внешней поверхности оболочки нагрузки, распределенной по области В = |(а,2): , - — <г<—| и изменяющейся
во времени по закону Р (т) Н (т).
Безразмерные нормальные перемещения оболочки определяются по формуле (12) с учетом (22), в которой поверхностный интеграл с учетом геометрии области В заменяется повторным интегралом:
w (a, z, х) = G (a, z, х)*** p (a, z, х) =
P b x 2 2
= J dt J d ^J G (a-£, z -Ç, х-1 ) p (£, Ç, t ) d Ç.
0 P _b 2 2
(23)
0
Труды МАИ. 2021.Выпуск № 120 Trudy MAI. 2021.Issues no.120
Для вычисления интеграла в (23) используем квадратурную формулу метода
прямоугольников :
Соотношения (21) и (24) с учетом функции Грина (18) позволяют исследовать пространственно-временные распространения нестационарных возмущений в неограниченной оболочке Кирхгофа-Лява при воздействии сосредоточенной и распределенной по внешней поверхности нестационарной нагрузки.
Построенные функции при переходе к размерным величинам дают возможность исследовать нестационарное напряженно-деформированное состояние протяженных анизотропных цилиндрических оболочек.
Нестационарные функции нормального прогиба (21) и (24) являются универсальными по отношению к свойствам материала оболочки, который может быть изотропным, трансверсально-изотропным, ортотропным или анизотропным. Оценим характер распространения нестационарных возмущений в неограниченной оболочке для нескольких вариантов симметрии упругой среды - изотропной, ортотропной и анизотропной.
Для вычисления необходимых упругих постоянных Ск1, входящих в функцию прогиба через инженерные константы, воспользуемся связью матрицы упругих постоянных С с матрицей податливости О:
(24)
3. Примеры расчетов
где [33]:
С = D
D
1 V21 - V31
Ei e2 E3
V12 1 - V32
Ei E2 E3
-V13 V23 1
Ei E2 E3
0 0 0
0 0 0
K1.12 K2.12 K3.12
E1 E2
0 0 K 12.1 ^ G12
0 0 K 12.2 G12
0 0 K12.3 G12
1 Л 31.23 0
G23 G31
Л23.31 1 0
G23 G31
0 0 1 G12 y
v
ik
v
kl
k
kl
к
kl л
Л
ik .lm
Л
lm .ik
El Ek
Ei Gkl Gik Glm
(25)
В (25): Е - модули упругости первого рода, ^ - модули упругости второго рода, - коэффициенты Пуассона, кш, кш - коэффициенты взаимного влияния, цЛ1т - коэффициенты Ченцова.
Исследуем нестационарную динамику изотропной и ортотропной цилиндрической оболочки с соотношением радиуса к толщине к = 50 при воздействии внезапно приложенной нагрузки вида (22), распределенной по области боковой поверхности оболочки размерами 7 = 0.5, а = л /24 , с амплитудой, меняющейся по закону:
Р (т) = Н (7 -т).
Труды МАИ. 2021.Выпуск № 120 Trudy MAI. 2021.Issues no.120
Изотропная среда. Исследуем нестационарную динамику металлической
цилиндрической оболочке с модулем Юнга E = 198 ГПа и коэффициентом Пуассона
ц = 0.28. Коэффициенты взаимного влияния и коэффициенты Ченцова - нулевые.
Компоненты тензора упругих постоянных согласно (24) примут значения (в Па):
сп = 2.531-1011, с12 = 9.844 • 1010, с16 = 0, с22 = 2.531-1011, с66 = 7.734 -1010, с26 = 0.
Соответствующие безразмерные упругие константы согласно (6) таковы:
C = 0.389, C = 0, C = 1.000, C = 0, C = 0.306.
На рис. 2 представлены пространственные зависимости безразмерного
нестационарного нормального перемещения изотропной цилиндрической оболочки в
безразмерные моменты времени т = 3 и т = 6.
т = 3 х = 6
Рис. 2. Пространственные зависимости нестационарного прогиба изотропной
оболочки при воздействии распределенной нагрузки
Ортотропная среда. Исследуем нестационарную динамику полимерной
композитной оболочки с симметричной относительно срединной плоскости схемой
Труды МАИ. 2021.Выпуск № 120 Trudy MAI. 2021.Issues no.120
армирования. Приведенные характеристики пакета примем равными (модули упругости в Па):
E = 2.09 x 1011, E2 = 9.45 x 109, E3 = 9.45 x 109, G12 = 5.5 x 109, G23 = 3.9 x 109, G31 = 5.5 x 109, v12 = 0.27, v23 = 0.4, v31 = 0.27.
Коэффициенты взаимного влияния и коэффициенты Ченцова - нулевые. Компоненты тензора упругих постоянных примут значения (в Па):
сп = 2.113 • 1011, с12 = 4.300 • 109, с16 = 0, с22 = 1.134 • 1010, с66 = 5.500 • 109, с26 = 0.
Соответствующие безразмерные упругие константы из (6): с = 0.020, C2 = 0, C3 = 0.054, C4 = 0, Cs = 0.026. На рис. 3 представлены пространственные зависимости безразмерного нестационарного нормального перемещения ортотропной цилиндрической оболочки в безразмерные моменты времени т = 3 и т = 6.
т = 3 т = 6
Рис. 3. Пространственные зависимости нестационарного прогиба ортотропной
оболочки при воздействии распределенной нагрузки
Труды МАИ. 2021.Выпуск № 120 Trudy MAI. 2021.Issues no.120
Из рис. 2 и рис. 3 видно, что в случае изотропного материала распределения
нормальных перемещений носят концентрический характер. В случае ортотропного
материала различимы две плоскости симметрии распространения нестационарных
возмущений.
Анизотропная среда. Исследуем нестационарную динамику анизотропной цилиндрической оболочки с отношением радиуса к толщине к = 500 и к = 50 при воздействии сосредоточенной нагрузки вида (19) с амплитудой, меняющейся по закону:
P (т) = sin (т) sch ( 3т). Упругие константы примем согласно [34] с учетом (24) (в Па):
сп = 9.55 • 1010, с12 = 2.890 • 1010, с16 = 4.470 • 1010, с22 = 2.590 • 1010, с66 = 3.270 • 1010, с26 = 1.560 •Ю10,
что соответствует следующим безразмерным константам из (6):
с = 0.302, C2 = 0.468, C3 = 0.271, C4 = 0.163, Cs = 0.342.
На рис. 4 и 5 представлены пространственные зависимости безразмерного
нестационарного нормального перемещения анизотропной цилиндрической
оболочки в безразмерные моменты времени т = 3 и т = 6 с отношениями радиуса к
толщине к = 500 и к = 50 соответственно.
т = 3 т = 6
Рис. 4. Пространственные зависимости нестационарного прогиба анизотропной
оболочки при воздействии сосредоточенной нагрузки, к = 500
т = 3 т = 6
Рис. 5. Пространственные зависимости нестационарного прогиба анизотропной
оболочки при воздействии сосредоточенной нагрузки, к = 50 В представленных на рис. 4 и рис. 5 результатах решения, полученных для анизотропного материала, просматривается смещение возмущений относительно образующих цилиндра и, в целом, несимметричный характер прогиба оболочки. При этом чем меньше отношение радиуса оболочки к ее толщине тем сильнее проявляется данный эффект в оболочке Кирхгофа-Лява.
Труды МАИ. 2021.Выпуск № 120 Trudy MAI. 2021.Issues no.120
Представленные на рис. 2 - рис. 5 результаты демонстрируют обобщенность
построенных нестационарной функции (21) и (24) нормального перемещения для упругой анизотропной неограниченной оболочки Кирхгофа-Лява в вопросе частных случаев анизотропии материала оболочки и универсальность по характеру действующей нагрузки.
Выводы
В работе получены новые численно-аналитические решения динамической задачи теории упругости для анизотропной упругой тонкой неограниченной цилиндрической оболочки Кирхгофа-Лява при воздействии на нее нестационарной сосредоточенной и распределенной нагрузки с переменными во времени амплитудами с использованием нового фундаментального решения для динамических функций Грина. С применением функции Грина построены пространственно-временные зависимости нормальных перемещений неограниченной оболочки Кирхгофа-Лява. Установлена связь между геометрическими размерами анизотропной оболочки и несимметричным характером распространения возмущений. Для демонстрации универсальности построенных функций нормальных перемещений представлены примеры расчетов для нескольких вариантов симметрии упругой среды, характер распространения нестационарных возмущений в которых позволил дать оценку адекватности решения.
Представленный подход к построению нестационарной функции прогиба при
переходе к размерным величинам открывает возможности для анализа
Труды МАИ. 2021.Выпуск № 120 Trudy MAI. 2021.Issues no.120
нестационарного напряженно-деформированного состояния протяженных цилиндрических оболочек с учетом различных вариантов анизотропии материала и закона распределения нестационарной нагрузки, как по координатам, так и по времени. Кроме того, применение построенной функции Грина возможно при решении широкого круга нестационарных задач для анизотропных оболочек, например, нестационарных контактных и обратных задач, а также исследования напряженно-деформированного состояния анизотропных оболочек конечных размеров, например, с применением метода компенсирующих нагрузок.
Список источников
1. Горшков А.Г., Медведский А.Л., Рабинский Л.Н., Тарлаковский Д.В. Волны в сплошных средах. - М.: ФИЗМАТЛИТ, 2004. - 472 с.
2. Богданович А.Е. Деформирование и прочность цилиндрических композитных оболочек при динамических нагрузках. - Рига: 1985. -560 с.
3. Богданович А.Е. Нелинейные задачи динамики цилиндрических композитных оболочек. - Рига: Зинатне, 1987. - 295 с.
4. Кошкина Т.Б. Деформирование и прочность подкрепленных композитных цилиндрических оболочек при динамических сжимающих нагрузках. - Рига: Академия наук Латвийской ССР, 1984. - 180 с.
5. Сибиряков А.В. Динамика слоистых композиционных пластин и оболочек при импульсном нагружении: дисс... д.т.н. - М.: 2002. - 319 с.
Труды МАИ. 2021.Выпуск № 120 Trudy MAI. 2021.Issues no.120
6. Tarlakovskii D.V., Fedotenkov G.V. Two-Dimensional Nonstationary Contact of Elastic Cylindrical or Spherical Shells // Journal of Machinery Manufacture and Reliability, 2014, vol. 43, no. 2, pp. 145-152. DOI: 10.3103/S105261881401017
7. Tarlakovskii D.V., Fedotenkov G.V. Nonstationary 3D Motion of an elastic Spherical Shell // Mechanics of Solids, 2015, vol. 50, no. 2, pp. 208-217. DOI: 10.3103/S0025654415020107
8. Mikhailova E.Yu., Tarlakovskii D.V., Fedotenkov G.V. The impact of liquid filled concentric spherical shells with a rigid wall // Shell Structures: Theory and Applications, 2017, vol. 4, pp. 305-308. DOI: 10.1201/9781315166605-68
9. Mikhailova E.Yu., Tarlakovskii D.V., Fedotenkov G.V. Transient contact problem for liquid filled concentric spherical shells and a rigid barrier // Proceedings of the First International Conference on Theoretical, Applied and Experimental Mechanics, 2019, pp. 385-386. DOI: 10.1007/978-3-319-91989-8 92
10. Fedotenkov G.V., Tarlakovsky D.V., Vahterova Y.A. Identification of non-stationary load upon timoshenko beam // Lobachevskii Journal of Mathematics, 2019, vol. 40, no. 4, pp. 439-447. DOI:10.1134/S1995080219040061
11. Вахтерова Я.А., Федотенков Г.В. Нестационарная обратная задача по идентификации дефектов в упругом стержне // XII Всероссийский съезд по фундаментальным проблемам теоретической и прикладной механики: сборник трудов (Уфа, 19-24 августа 2019). - Уфа: Башкирский государственный университет, 2019. Т.3. С. 878-880.
Труды МАИ. 2021.Выпуск № 120 Trudy MAI. 2021.Issues no.120
12. Okonechnikov A.S., Tarlakovski D.V., Ul'yashina A.N., Fedotenkov G.V. Transient reaction of an elastic half-plane on a source of a concentrated boundary disturbance // IOP Conference Series: Materials Science and Engineering, 2016, vol. 158, no 1, pp. 012073. DOI:10.1088/1757-899X/158/1/012073
13. Иванов С.В., Могилевич Л.И., Попов В.С. Продольные волны в нелинейной цилиндрической оболочке, содержащей вязкую жидкость // Труды МАИ. 2019. №2 105. URL: http://trudymai.ru/published.php?ID=104003
14. Нуштаев Д.В., Жаворонок С.И., Клышников К.Ю., Овчаренко Е.А. Численно -экспериментальное исследование деформирования и устойчивости цилиндрической оболочки ячеистой структуры при осевом сжатии // Труды МАИ. 2015. № 82. URL: http://trudymai.ru/published.php?ID=58589
15. Карпов В.В., Семенов А.А., Холод Д.В. Исследование прочности пологих ортотропных оболочек из углепластика // Труды МАИ. 2014. № 76. URL: http://trudymai.ru/published.php?ID=49970
16. Жигалко Ю.П., Садыкова М.М. Динамика, тонкой круглой пластинки при нестационарном локальном нагружении // Исследования по теории пластин и оболочек. 1990. № 20. С. 184-191.
17. Моргачев К.С. Нестационарная динамика кольцевой пластины Тимошенко переменой толщины // Вестник Самарского государственного технического университета. 2007. T. 15. № 2. С. 162-164.
18. Дьяченко Ю.Г. Нестационарная задача динамики пластин переменного сечения
в уточненной постановке: Автореферат дисс.... к.ф-м.н. - Саратов: СГУ, 2008. - 19 с.
26
Труды МАИ. 2021.Выпуск № 120 Trudy MAI. 2021.Issues no.120
19. Шевченко В.П. Ветров О.С. Динамика ортотропной пластины под действием локальных внезапно приложенных нагрузок // Труды института прикладной математики и механики Национальной академии наук Украины. 2011. T. 22. С. 207215.
20. Nayfeh A.H., Chimenti D.E. Free Wave Propagation in Plates of General Anisotropic Media // Journal of applied mechanics-transactions of the ASME, 1989, vol. 56, no. 4, pp. 881 - 886. DOI: 10.1115/1.3176186
21. Wahab M.A., Jabbour T., Davies P. Prediction of impact damage in composite sandwich plates // Materiaux & Techniques, 2019, vol. 107, no. 2. DOI: 10.1051/mattech/2019006
22. Daros C.H. The dynamic fundamental solution and BEM formulation for laminated anisotropic Kirchhoff plates // Engineering analysis with boundary elements, 2015, vol. 54, no. 2, pp. 19 - 27. DOI: 10.1016/j.enganabound.2015.01.001
23. Igumnov L.A., Markov I.P. A boundary element approach for 3d transient dynamic problems of moderately thick multilayered anisotropic elastic composite plates // Materials physics and mechanics, 2018, vol. 37, no. 1, pp. 79-83. DOI: 10.18720/MPM.3712018 11
24. Sahli A., Boufeldja S., Kebdani S., Rahmani O. Failure analysis of anisotropic plates by the boundary element method // Journal of mechanics, 2014, vol. 30, no. 6, pp. 561-570. DOI: 10.1017/jmech.2014.65
25. Starovoitov E.I., Leonenko D.V. Deformation of an Elastoplastic Three-Layer Circular Plate in a Temperature Field // Mechanics of Composite Materials, 2019, vol. 55, no. 4, pp. 503-512. DOI: 10.1007/s11029-019-09829-6
Труды МАИ. 2021.Выпуск № 120 Trudy MAI. 2021.Issues no.120
26. Ryazantseva M.Y., Starovoitov E.I. Static and Dynamic Models of Bending for Elastic Sandwich Plates // Structural Integrity, 2019, vol. 8, pp. 294-297. DOI: 10.1007/9783-030-21894-2 54
27. Starovoitov E.I., Leonenko D.V., Tarlakovskii D.V. Thermoelastic Deformation of a Circular Sandwich Plate by Local Loads // Mechanics of Composite Materials, 2018, vol. 54, no. 3, pp. 299-312. DOI: 10.1007/s11029-018-9740-x
28. Starovoitov E.I., Leonenko D.V. Vibrations of circular composite plates on an elastic foundation under the action of local loads // Mechanics of Composite Materials, 2016, vol. 52, no. 5, pp. 665-672. DOI: 10.1007/s11029-016-9615-y
29. Тарлаковский Д.В., Федотенков Г.В. Общие соотношения и вариационные принципы математической теории упругости: Учебное пособие. - М.: Изд-во МАИ-ПРИНТ, 2009. - 112 с.
30. Михайлова Е.Ю., Тарлаковский Д.В., Федотенков Г.В. Упругие пластины и пологие оболочки: Учебное пособие. - М.: Изд-во МАИ, 2018. - 92 с.
31. Дёч Г. Руководство к практическому применению преобразования Лапласа и Z-преобразований. - М.: Изд-во Наука, 1971. - 288 с.
32. Бахвалов Н.С., Жидков Н.П., Кобельков Г.М. Численные методы. - М.: Наука, 1975. - 630 с.
33. Ашкенази Е.К. Анизотропия древесины и древесных материалов. - М.: Лесная промышленность, 1978. - 224 с.
Труды МАИ. 2021.Выпуск № 120 Trudy MAI. 2021.Issues no.120
34. Игумнов Л.А., Марков И.П., Пазин В.П. Гранично-элементное решение краевых задач трехмерной анизотропии теории упругости // Вестник Нижегородского университета им. Н.И. Лобачевского. 2013. № 1(3). С. 115 - 119.
References
1. Gorshkov A.G., Medvedskii A.L., Rabinskii L.N., Tarlakovskii D.V. Volny v sploshnykh sredakh (Waves in Continuous Media Study guide: for universities), Moscow, FIZMATLIT, 2004, 472 p.
2. Bogdanovich A.E. Deformirovanie i prochnost' tsilindricheskikh kompozitnykh obolochek pri dinamicheskikh nagruzkakh (Deformation and Strength for Cylindrical Composite Shells under Dinamic Loads), Riga, 1985, 560 p.
3. Bogdanovich A.E. Nelineinye zadachi dinamiki tsilindricheskikh kompozitnykh obolochek (Non-linear Problems of Dynamics of Cylindrical Composite Shells), Riga, Zinatne, 1987, 295 p.
4. Koshkina T.B. Deformirovanie i prochnost' podkreplennykh kompozitnykh tsilindricheskikh obolochek pri dinamicheskikh szhimayushchikh nagruzkakh (Deformation and Strength of Orthotropic Cylindrical Shells Subjected to Dynamic Compressive Loads), Riga, Akademiya nauk Latviiskoi SSR, 1984, 180 p.
6. Tarlakovskii D.V., Fedotenkov G.V. Two-Dimensional Nonstationary Contact of Elastic Cylindrical or Spherical Shells, Journal of Machinery Manufacture and Reliability, 2014, vol. 43, no. 2, pp. 145-152. DOI: 10.3103/S105261881401017
Труды МАИ. 2021.Выпуск № 120 Trudy MAI. 2021.Issues no.120
7. Tarlakovskii D.V., Fedotenkov G.V. Nonstationary 3D Motion of an elastic Spherical Shell, Mechanics of Solids, 2015, vol. 50, no. 2, pp. 208-217. DOI: 10.3103/S0025654415020107
8. Mikhailova E.Yu., Tarlakovskii D.V., Fedotenkov G.V. The impact of liquid filled concentric spherical shells with a rigid wall, Shell Structures: Theory and Applications, 2017, vol. 4, pp. 305-308. DOI: 10.1201/9781315166605-68
9. Mikhailova E.Yu., Tarlakovskii D.V., Fedotenkov G.V. Transient contact problem for liquid filled concentric spherical shells and a rigid barrier, Proceedings of the First International Conference on Theoretical, Applied and Experimental Mechanics, 2019, pp. 385-386. DOI: 10.1007/978-3-319-91989-8 92
10. Fedotenkov G.V., Tarlakovsky D.V., Vahterova Y.A. Identification of non-stationary load upon timoshenko beam, Lobachevskii Journal of Mathematics, 2019, vol. 40, no. 4, pp. 439-447. DOI: 10.1134/S1995080219040061
11. Vakhterova Ya.A., Fedotenkov G.V. XII Vserossiiskii s"ezd po fundamental'nym problemam teoreticheskoi i prikladnoi mekhaniki: sbornik trudov. Ufa, Bashkirskii gosudarstvennyi universitet, 2019, vol. 3, pp. 878-880.
12. Okonechnikov A.S., Tarlakovski D.V., Ul'yashina A.N., Fedotenkov G.V. Transient reaction of an elastic half-plane on a source of a concentrated boundary disturbance, IOP
Conference Series: Materials Science and Engineering, 2016, vol. 158, no 1, pp. 012073. DOI:10.1088/1757-899X/158/1/012073
13. Ivanov S.V., Mogilevich L.I., Popov V.S. Trudy MAI, 2019, no. 105. URL: http://trudymai.ru/eng/published.php?ID=104003
Труды МАИ. 2021.Выпуск № 120 Trudy MAI. 2021.Issues no.120
14. Nushtaev D.V., Zhavoronok S.I., Klyshnikov K.Yu., Ovcharenko E.A. Trudy MAI, 2015, no. 82. URL: http://trudymai.ru/eng/published.php?ID=58589
15. Karpov V.V., Semenov A.A., Kholod D.V. Trudy MAI, 2014, no. 76. URL: http://trudymai.ru/eng/published.php?ID=49970
16. Zhigalko Yu.P., Sadykova M.M. Issledovaniyapo teorii plastin i obolochek, 1990, no. 20, pp. 184-191.
17. Morgachev K.S. Vestnik Samarskogo gosudarstvennogo tekhnicheskogo universiteta, 2007, vol. 15, no. 2, pp. 162-164.
18. D'yachenko Yu.G. Nestatsionarnaya zadacha dinamikiplastin peremennogo secheniya v utochnennoi postanovke (Transient problem of dynamics of plates of variable cross-section in a refined formulation), author's abstract, Saratov, SGU, 2008, 19 p.
19. Shevchenko V.P. Vetrov O.S. Trudy instituta prikladnoi matematiki i mekhaniki Natsional'noi akademii nauk Ukrainy, 2011, vol. 22, pp. 207-215.
20. Nayfeh A.H., Chimenti D.E. Free Wave Propagation in Plates of General Anisotropic Media, Journal of applied mechanics-transactions of the ASME, 1989, vol. 56, no. 4, pp. 881 - 886. DOI: 10.1115/1.3176186
21. Wahab M.A., Jabbour T., Davies P. Prediction of impact damage in composite sandwich plates, Materiaux & Techniques, 2019, vol. 107, no. 2. DOI: 10.1051/mattech/2019006
22. Daros C.H. The dynamic fundamental solution and BEM formulation for laminated anisotropic Kirchhoff plates, Engineering analysis with boundary elements, 2015, vol. 54, no. 2, pp. 19 - 27. DOI: 10.1016/j.enganabound.2015.01.001
Труды МАИ. 2021.Выпуск № 120 Trudy MAI. 2021.Issues no.120
23. Igumnov L.A., Markov I.P. A boundary element approach for 3d transient dynamic problems of moderately thick multilayered anisotropic elastic composite plates, Materials physics and mechanics, 2018, vol. 37, no. 1, pp. 79-83. DOI: 10.18720/MPM.3712018 11
24. Sahli A., Boufeldja S., Kebdani S., Rahmani O. Failure analysis of anisotropic plates by the boundary element method, Journal of mechanics, 2014, vol. 30, no. 6, pp. 561-570. DOI: 10.1017/jmech.2014.65
25. Starovoitov E.I., Leonenko D.V. Deformation of an Elastoplastic Three-Layer Circular Plate in a Temperature Field, Mechanics of Composite Materials, 2019, vol. 55, no. 4, pp. 503-512. DOI: 10.1007/s11029-019-09829-6
26. Ryazantseva M.Y., Starovoitov E.I. Static and Dynamic Models of Bending for Elastic Sandwich Plates, Structural Integrity, 2019, vol. 8, pp. 294-297. DOI: 10.1007/978-3-03021894-2 54
27. Starovoitov E.I., Leonenko D.V., Tarlakovskii D.V. Thermoelastic Deformation of a Circular Sandwich Plate by Local Loads, Mechanics of Composite Materials, 2018, vol. 54, no. 3, pp. 299-312. DOI: 10.1007/s11029-018-9740-x
28. Starovoitov E.I., Leonenko D.V. Vibrations of circular composite plates on an elastic foundation under the action of local loads, Mechanics of Composite Materials, 2016, vol. 52, no. 5, pp. 665-672. DOI: 10.1007/s11029-016-9615-y
29. Tarlakovskii D.V., Fedotenkov G.V. Obshchie sootnosheniya i variatsionnye printsipy matematicheskoi teorii uprugosti (General ratios and variation principles of mathematical theory of elasticity) Moscow, Izd-vo MAI-PRINT, 2009, 112 p.
Труды МАИ. 2021.Выпуск № 120 Trudy MAI. 2021.Issues no.120
30. Mikhailova E.Yu., Tarlakovskii D.V., Fedotenkov G.V. Uprugie plastiny i pologie obolochki (Elastic Plates and Shallow Shells), Moscow, Izd-vo MAI, 2018, 92 p.
31. Dech G. Rukovodstvo k prakticheskomu primeneniyu preobrazovaniya Laplasa i Z-preobrazovanii (Guide to the Practical application of Laplace and Z-transforms), Moscow, Izd-vo Nauka, 1971, 288 p.
32. Bakhvalov N.S., Zhidkov N.P., Kobel'kov G.M. Chislennye metody (Numerical Methods), Moscow, Nauka, 1975, 630 p.
33. Ashkenazi E.K. Anizotropiya drevesiny i drevesnykh materialov (Wood and Wood Materials Anisotropy), Moscow, Lesnaya promyshlennost', 1978, 224 p.
34. Igumnov L.A., Markov I.P., Pazin V.P. Vestnik Nizhegorodskogo universiteta im. N.I. Lobachevskogo, 2013, no. 1(3), pp. 115 - 119.
Статья поступила в редакцию 23.09.2021; одобрена после рецензирования 30.09.2021; принята к публикации 22.10.2021.
The article was submitted 23.09.2021; approved alter reviewieng 30.09.2021; accepted for publication 22.10.2021.