2020. 16(3). 219-225 Строительная механика инженерных конструкций и сооружений Structural Mechanics of Engineering Constructions and Buildings
HTTP://JOURNALS.RUDN.RU/STRUCTURAL-MECHANICS
DOI 10.22363/1815-5235-2020-16-3-219-225 научная статья
УДК 539.376
Устойчивость геометрически нелинейных пластинчатых систем под действием динамических нагрузок
С.П. Иванов1'2*, А.С. Иванова1, О.Г. Иванов1
1Поволжский государственный технологический университет, Российская Федерация, 424000, Йошкар-Ола, пл. Ленина, 3 2Марийский государственный университет, Российская Федерация, 424000, Йошкар-Ола, пл. Ленина, 1 *8р-геапоу@шаП.ги
История статьи:
Поступила в редакцию: 28 января 2020 г. Доработана: 04 марта 2020 г. Принята к публикации: 24 марта 2020 г.
Для цитирования
Иванов С.П., Иванова А.С., Иванов О.Г. Устойчивость геометрически нелинейных пластинчатых систем под действием динамических нагрузок // Строительная механика инженерных конструкций и сооружений. 2020. Т. 16. № 3. С. 219-225. http://dxdoi.org/10.22363/1815-5235-2020-16-3-219-225
Аннотация
Актуальность. Односвязные и многосвязные пластинчатые системы часто используют в строительстве, авиастроении, кораблестроении, машиностроении, приборостроении. Вследствие этого исследование устойчивости геометрически нелинейных пространственных пластинчатых систем является актуальной темой как при действии статических, так и динамических нагрузок. Но, несмотря на значительные достижения в этой области, имеется еще много нерешенных проблем. Так, запросы вышеуказанных областей применения тонкостенных пространственных систем требуют дальнейшего исследования статической и динамической устойчивости. Цель - разработка метода расчета на устойчивость геометрически нелинейных пластинчатых систем типа призматических оболочек под действием динамических сжимающих нагрузок. Методы. Рассматривается пластинчатая система, на которую в продольном направлении действуют динамические сжимающие нагрузки. Учитываются гипотезы Кирхгофа - Лява. Геометрическая нелинейность вводится через соотношения между деформациями и перемещениями. Диаграмма деформирования материала - линейная. Перемещение точек в нормальном направлении к срединной плоскости пластин определяется в виде разложения по Власову. Для решения задачи используются энергетический метод и вариационный метод Власова. Экстремальное значение полной энергии определяется с использованием уравнений Эйлера - Лагранжа. В результате получена система основных нелинейных дифференциальных уравнений для исследования потери устойчивости пластинчатой системы под действием динамических сжимающих нагрузок. Результатыи Разработанный метод применяется для расчета на устойчивость геометрически нелинейной призматической оболочки с замкнутым контуром поперечного сечения при центральном сжатии под действием динамической нагрузки. Края оболочки опираются на диафрагмы. Исследуется потеря устойчивости призматической оболочки в продольном направлении по одной и двум полуволнам синусоиды. Численное интегрирование нелинейных дифференциальных уравнений выполняется методом Рунге - Кутта. По результатам расчетов построены графики зависимости относительной величины прогиба от динамического коэффициента. Рассмотрено влияние скорости изменения сжимающего напряжения, начального несовершенства оболочки и других параметров на критерий динамической устойчивости пластинчатой системы.
Ключевые слова: динамическая устойчивость, геометрическая нелинейность, пластинчатая система, призматическая оболочка, вариационный метод Власова_
Иванов Сергей Павлович, доктор технических наук, профессор, заведующий кафедрой сопротивления материалов и прикладной механики ПГТУ; профессор кафедры электромеханики МарГУ; eLIBRARY SPIN-код: 5963-6739.
Иванова Анастасия Сергеевна, старший преподаватель, кафедра сопротивления материалов и прикладной механики; eLIBRARY SPIN-код: 9568-3451.
Иванов Олег Геннадьевич, кандидат технических наук, доцент, доцент кафедры сопротивления материалов и прикладной механики; eLIBRARY SPIN-код: 5052-9077.
© Иванов С.П., Иванова А.С., Иванов О.Г., 2020
ficci © I This work is licensed under a Creative Commons Attribution 4.0 International License ¿^HihI https://creativec0mm0ns.0rg/licenses/by/4.Q/
Введение
где
Для расчета на устойчивость пластин и оболочек используются различные аналитические, численно-аналитические и численные методы. К численным методам расчета относится метод конечных элементов (МКЭ), который совершенно непригоден для ручного счета. Вариационный метод В.З. Власова [1; 2] относится ко второму виду методов расчета, где в отличие от МКЭ конечные уравнения можно реализовать в замкнутом виде. Данный метод эффективно применяется для расчета пластин и оболочек на прочность как при статических, так и при динамических воздействиях [1]. Также он применяется и при исследованиях на устойчивость [1; 3], доказывающих, что метод В.З. Власова достаточно эффективен для изучения физически и геометрически нелинейных пластинчатых систем.
Тема, связанная с исследованием устойчивости пространственных пластинчатых систем, является одной из важнейших проблем механики деформируемого твердого тела. При движении таких систем в жидкости или воздушной среде (например, при воздействии ветровой нагрузки) может возникать динамическая потеря устойчивости [5; 6]. В настоящее время в области расчетов на устойчивость тонкостенных конструкций имеется достаточно большое число публикаций [6-15]. Анализируя указанные статьи, можно сделать общие выводы. Такого вида работы были выполнены ранее А.С. Воль-миром [5; 6] и П.А. Лукашем [13] только без применения метода МКЭ. Данная работа отличается от статей [14; 15] тем, что здесь не определяется верхняя и нижняя критические силы и не решается статическая задача. Критическое значение динамической силы определяется из условия бурного выпучивания. Вычисляется отношение динамической нагрузки к статической критической нагрузке.
1. Разработка математической модели
В отличие от работ [3; 4] настоящая статья посвящена исследованию устойчивости пластинчатых систем (типа призматических оболочек) под действием динамической нагрузки P(t) с учетом геометрической нелинейности (рис. 1).
Через u = u(x, s, t); v = v(x, s, t); w = w(x, s, t) (t - время) обозначим перемещения точки М срединной поверхности оболочки в направлении координатных осей x, s, z (рис. 1). Геометрическую нелинейность вводим через известные соотношения [1]:
0 du 0 dv 0 du dv 1 1 dw 2 i 1 ydw 2
ex =— ;es =— ;e„ =--1--;ex = — (-) ;es = — (-) ;
dx ds ds dx 2 dx 2 ds
s = e° + e1 - zx ;s = e° + e1
x x x r^x? s s s
-zx s ;
d2w
d2w
d2w
1 dw dw
e xs = ' X x = » 2 ' X s = -л 2 ' X xs = ^ ^
dx ds dx ds dxds
S xs eos + eis 2 ZX xs ,
(1)
Рис. 1. Общая схема многосвязной пластинчатой системы [Figure 1. The general diagram of a multi-connected plate system]
Определим полную энергию L системы:
L = П + К, (2)
состоящей из потенциальной П и кинетической К энергий:
1 (Мх Xx + Ms X s + 2M +j ^ W2 -I xs + P(tdxds;
2 l+Жл + Nses + 2Nxsexs J 2 J
К = P \ 2 ( + v2 + w2 )dt. (3)
g 2
Внутренние усилия (изгибающие моменты Мх, Ms, действующие соответственно в продольном х и поперечном s направлениях, крутящие моменты Mxs, нормальные Nx, Ns и сдвигающие Nxs силы) определяются по следующим формулам [1]:
Mx = -wxx + vwss); Nx = Ki (ux + 2 j;
Ms = -D( + vwxx);
Ns = Ki (vs + 2 w2 j; Mxs = Mx = -D(1 -v) ws;
NxS = G5(us + vx + WxWs), (4)
где в формулах (3) и (4) индексы при перемещениях u, v, w обозначают частные производные по переменным x, s, t:
е = е0 + е1; е = е0 + е1; е = е0 + е1 ;
X X s s s > xs xs xs>
D =
E53
-; K =
Eb
12(1 — V) 1 1 -V2
Здесь Е, О - модули упругости первого и второго рода; 5 - толщина пластин оболочки.
Все перемещения запишем в виде ряда [1; 2]:
и(х, ^ *) = 2 и (* )Ф. (х, я);
у( х, я, *) = ^Ук (* )у (х, я);
к
w( х, я, *) = 2^ (*) / (х, я);
а
(г = 1, 2, ... , т; к, а = 1, 2, ... , п), (5)
где Ц-(0, ЩО, - обобщенные перемещения, которые зависят только от времени и определяются из решения задачи; ф.(х,5), у^), /а(х,я) - функции распределения, которые задаются заранее [1; 3].
Число обобщенных перемещений можно сократить, используя условия совместности деформаций в узловых точках контура поперечного сечения оболочки, и при а = к [1]
ад = vk(t).
(6)
С учетом соотношения (6) определим минимум функционала (2) [1]:
8L d 8L 8L d 8L
■ = 0;----= 0.
8U, dt 8Uut
8Vk dt 8Vk,
(7)
В развернутой форме уравнения (7) принимают следующий вид:
Y^a- -ь,и -Ys* v -gGYßiUtt,=фУ i k gG i
Y
"Yie«
m.
+rhk ehk -a G
(8)
V +Yci U - gGYdk vhtt + Qk = % i g к
Tn«
(i, j = 1,2,3,..., m; k, h = 1,2,3,..., n),
где 71 =7 /(1 — V ), у = Е / О - отношение модуля
упругости Е к модулю сдвига О, О = Е /[2(1 + V)] ;
*
величина а - длина контура поперечного сечения оболочки, на который действует динамическая нагрузка Р(*); нагрузка Qк позволяет учитывать начальное несовершенство оболочки.
В системе уравнений (8) в функциях ип , ,
и и, У^ и индексы после запятой указывают на дифференцирование по времени *.
Коэффициенты уравнений (8) имеют вид [1; 3]:
а. =яф 1. хф., х5аяах;
х Б
ъР=\\ф 1 ф., ^ Ы^х;
х 3 х 3
сы = \\ук хФ,.Мах;
х ^
еНк =\\/ хх/к, хх^х
х 3
гш =я^, х у, х5аяах;
х 3
анк =\\(/хх/к, хх + у х 3
4 = Ц Л ххЛ,«^;
х £
пък = Ц А^Л^^;
j=а.
12
(9)
В выражениях (9) ф.- х = 8ф. / 8х, ф j,y = 8ф j / 8y,
fh,XX =82 fh /8x2, - • .
Правые части Ф, и Фh уравнений (8) имеют следующий вид:
Ф; = -Yi ^ YVkfk, x YVkfk ,xx9 jbds+,f YVkfk, x YVkfk ,x9,, sbds;
Ф h = -Yi[^ YU- ф.-, xx YYVkfk, xfhbds+f YU> Ф.-.X YYVkfk, fbds+2 .f (Vkfk, x )2 YVkfk,«fh] -
-f [Y U- ф.-.xs YVkfk,x+Yv^kx Y Vkfk,s + YVkfk, x Y Vkfk,s Y Vkfk,xs +Y Vkfk,xx (Y Vkfk,s )2 ]fhbds -
-f (]YVkfk,xx YVkfk, s+YVkfk,x YYVVkfk, xs )Vh bds. (10)
Динамика конструкций и сооружений 221
Получена общая система (m + n) дифференциальных уравнений (8). Данные уравнения пригодны для исследования устойчивости призматических оболочек под действием динамических нагрузок P(t), которые могут изменятся по различным законам [6]. А.С. Вольмир в работе [6] отмечает: «Говоря о непрерывном увеличении нагрузки, мы имеем при этом в виду, что важный для нас процесс прощелкива-ния оболочки происходит на восходящем участке диаграммы нагружения; дальнейшим поведением конструкции мы не интересуемся». Решающее влияние на бурное выпучивание оболочки оказывает скорость возрастания нагрузки. Далее это будет продемонстрировано на примере расчета призматической системы.
2. Алгоритм решения задачи
Уравнения (8) для решения конкретной задачи можно реализовать двумя способами.
Первый способ заключается в непосредственном интегрировании по времени нелинейных дифференциальных уравнений при заданной нагрузке P(t), изменяющейся во времени. При достижении критического значения времени t колебательный процесс нарушается, происходит резкое увеличение амплитуды колебаний, что соответствует бурному выпучиванию оболочки. Далее определяется соответствующее значение динамической нагрузки по формуле, которая задана для P(t).
Второй способ заключается в использовании новой переменной t* [6], которая связана с временной сжимающей нагрузкой P(t) и Эйлеровой критической нагрузкой Ркр следующей зависимостью t* = P(t)/PRp. А.С. Вольмир [6] отмечает: «Под величиной Кд понимается отношение динамической "критической" нагрузки к верхней статической критической нагрузке, вычисленной для идеальной оболочки или пластинки с теми же параметрами». Далее по преобразованным уравнениям строятся графики зависимости прогиба от t*.
Для численного интегрирования дифференциальных уравнений используем метод Рунге - Кут-та по программе, составленной на языке Фортран.
3. Пример расчета
Применим второй способ решения по полученным уравнениям (8) для исследования динамической устойчивости односвязной призматической оболочки при центральном сжатии (рис. 2) [1]. Оболочка опирается торцами на диафрагмы, которые считаются абсолютно жесткими в своей плоскости и абсолютно гибкими из плоскости.
Приняты следующие геометрические параметры замкнутой призматической оболочки: а = 2,4 м;
Ь = 1,6 м; 5 = 0,1 м; I = 31,4 м. Коэффициент Пуассона материала оболочки V = 0,2; объемный вес материала р = 20 кН/м3.
Положим, что динамическая нагрузка изменяется по линейному закону:
р(0 = к• г• а -5, (11)
где к - величина, которая характеризует скорость изменения сжимающего напряжения; а* - периметр контура поперечного сечения.
._т
*
Рис. 2. Односвязная призматическая оболочка: а - схема оболочки под действием динамической нагрузки; б - поперечное сечение оболочки [Figure 2. The single-connected prismatic shell: a - the diagram of shell under the action of dynamic load; б — the cross section of shell]
Для данной замкнутой оболочки при потере устойчивости от центрального сжатия по крутильной форме перемещения можно представить в виде [1]
u(x, s, t) = U1(t) Ф1 (x, s);
v( x, s, t) = V (t) Vi (x, s) + V (t) y2 (x, s);
w( x, s, t) = W (t )fi (x, s) + W2 (t )f2 (x, s), (12)
где Ui(t), Vi(t), V2(t), Wi(t), W2(t) - обобщенные перемещения в направлении осей x, s, z; 9i(x,s), Vi(x,s), V2(x,s), fi(x,s), f(x,s) - координатные функции, которые задаются при крутильной форме потери устойчивости в поперечном направлении (рис. 2, a, б).
Для данной оболочки в случае потери устойчивости по синусоиде в направлении оси х координатные функции можно записать следующим образом:
9i (x, s) = ф1 (s) ^s Xx; v1 (x, s) = v1 (s) sin Xx; V2 (x, s) = V2 (s) sin Xx; f (x, s) = f (s) sin Xx;
f2( x s) = f2( s) sin ^ (13)
где X = m1n/l, l - длина оболочки; m1 - число полуволн.
Функции распределения ф^), у^), ^(я), /¡(я), /.(я), зависящие от переменной я (рис. 2, б), представлены в работе [1]. Обход замкнутого контура ведется начиная с левого верхнего угла (рис. 2, б). Используя данные функции по формулам (9) с учетом (13) определяем коэффициенты уравнения (8).
Дифференциальные уравнения устойчивости (8) для данной призматической оболочки под действием динамической нагрузки принимают вид
(T1«11 -bu)Ui -cuVi -C12V2 +—aiiUl tt = ф
gG
2y b2 + a2 , r, +—— m„--P(t)
11 11 1 + v
ru +
2Y 1 + v
4
b2 + a2
V +
m11 -
P(t)
V2 +
+ C11U1 + gG (d11V1tt + d12V2tt ) = Ф2;
2y b2 - a2 r21 +—— m21--P(t)
'21 1 . '"21
1 + v
r +
22
2Y 1 + v
m22 - hn22
V +
b2 + a2
P (t)
V2 +
+ C21U1 + gG d'tt + d22V2,tt ) = Фз;
(14)
Здесь Ф1, Ф2, Ф3 учитывают геометрическую нелинейность системы и определяются по формулам (10) при г = 1; к, к = 1, 2.
Уравнения (14) можно решить первым способом, непосредственно интегрируя их по времени.
Конечную реализацию уравнений (14) проводим, используя второй способ, вводя следующие обозначения:
P(t) V V
P„ = 0,00134E; t* = ; z = —; Z2 = —;
P 5 2 5
- кр
S * =
52 (a* s)2
P • ^ • g • G'
(15)
где Ркр - величина статической критической нагрузки; ** - безразмерный параметр времени; ^ - относительная величина прогиба пластин оболочки; - параметр скорости изменения напряжения (нагрузки) [1].
Пренебрегаем депланацией оболочки, вводим новую переменную **, с учетом (15) после некоторых преобразований уравнения (14) приводим к следующему виду:
66,242 • (0,0549 -1* )Z1 + 25,478 • (0,0607 - t)Z2 + +S* (13,65Z1t.t. + 4,78Z2/,) + Q1 = Ф2(С1,С2); 25,478 • (0,0607 -1* )Z1 + 66,242 • (0,0549 -1% + +S * (4,78Z1t.t. + 13,65Z2/,) + Q2 = Фз(С1,С2), (16)
где с = с = Где C1,t*t* = ,2 , C2,t*t* =
dt
dt
2
По результатам численного интегрировании дифференциальных уравнений (16) с помощью метода Рунге - Кутта построены графики зависимости перемещения Z = v/5 узловой точки поперечного сечения, расположенного в середине пролета, от параметра времени t* = Р^)/Ркр (рис. 3).
На рис. 3 графики 7, 2, 3 построены соответственно для скоростей S* = 10, 50, 100 изменения сжимающей нагрузки при потере устойчивости по одной полуволне синусоиды m1 = 1 и параметру, учитывающему величину начального несовершенства оболочки Q1 = Q2 = 0,1. График 4 соответствует параметру S = 50 при m1 = 1, Q1 = Q2 = 0,01. График 5 построен для S = 50 при m1 = 2, Q1 = Q2 = 0,1.
Рис. 3. Графики зависимости перемещения Z от параметра времени t* [Figure 3. Graphs of the displacement Z against the time parameter t*]
Из анализа графиков следует, что с увеличением скорости действия нагрузки S величина динамического коэффициента Кд = t* возрастает значительно (см. графики 1, 2 и 3). При увеличении числа полуволн mi величина Кд также значительно возрастает (видно из сравнения графиков 2 и 5). При уменьшении начального несовершенства значение Кд увеличивается (см. графики 2 и 4).
Для сравнения выполнен расчет и первым способом. Оба способа дают близкие результаты.
Заключение
Разработана математическая модель расчета на устойчивость пластинчатых систем при действии динамической нагрузки с учетом геометрической
нелинейности. Получена система нелинейных дифференциальных уравнений для исследования динамической устойчивости призматических оболочек.
Предложен алгоритм конечной реализации уравнений двумя способами. Предпочтение отдано второму способу, так как его использование сразу выявляет, во сколько раз динамическая нагрузка превышает статическую критическую нагрузку (на основе величины
В качестве примера выполнен расчет на устойчивость оболочки замкнутого контура при действии динамической нагрузки, изменяющейся по линейному закону. Рассмотрено влияние скорости изменения сжимающего напряжения, начального несовершенства оболочки и числа полуволн на критерий динамической устойчивости оболочки.
Список литературы
1. Иванов С.П., Иванова А.С. Приложение вариационного метода В.З. Власова к решению нелинейных задач пластинчатых систем: монография. Йошкар-Ола: ПГТУ, 2015. 248 с.
2. Власов В.З. Тонкостенные пространственные системы. М.: Госстройиздат, 1958. 502 с.
3. Иванов С.П., Иванова А.С. Динамическая устойчивость физически нелинейных пластинчатых систем // Строительная механика инженерных конструкций и сооружений. 2014. № 4. С. 11-20.
4. Иванов С.П., Иванов О.Г., Иванова А.С. Динамическая устойчивость физически нелинейных пластинчатых систем при сжатии в двух направлениях // Строительная механика инженерных конструкций и сооружений. 2018. Т. 14. № 2. С. 132-141.
5. Вольмир А.С. Устойчивость деформируемых систем. М.: Наука, 1967. 984 с.
6. Вольмир А.С. Нелинейная динамика пластинок и оболочек. М.: Наука, 1972. 432 с.
7. Хамитов Т.К., Фатыхова Р.Р. Об устойчивости упругопластической цилиндрической оболочки при продольном ударе // Известия КГАСУ. 2016. № 4 (38). С. 490-496.
8. Трушин С.И., Журавлева ТА., Сысоева Е.В. Устойчивость нелинейно-деформируемых цилиндрических оболочек при действии неравномерных нагрузок // Строительная механика инженерных конструкций и сооружений. 2013. № 2. С. 3-10.
9. Трушин С.И., Журавлева ТА., Сысоева Е.В. Динамическая потеря устойчивости нелинейно-деформируемых сетчатых пластин из композиционного материала с различными конфигурациями решетки // Научное обозрение. 2016. № 4. С. 44-51.
10. Vescovini R., Dozio L. Exact refined buckling solutions for laminated plates under uniaxial and biaxial loads // Composite Structures. 2015. Vol. 127. Рр. 356-368.
11. Nazarimofrad E., Barkhordar A. Buckling analysis of orthotopic rectangular plate resting on Pasternak elastic foundation under biaxial in-plane loading // Mechanics of Advanced Materials and Structures. 2016. Vol. 23. No. 10. Рр. 1144-1148.
12. Ruocco E., Reddy J.N. A closed-form solution for buckling analysis of orthotopic Reddy plates and prismatic plate structures // Composites Part B: Engineering. 2019. Vol. 169. Pp. 258-273.
13. Лукаш ПА. Основы нелинейной строительной механики. M.: Стройиздат, 1978. 204 c.
14. Косицын С.Б., Акулич В.Ю. Определение критической нагрузки потери устойчивости стержневой и плоской моделей круговой цилиндрической оболочки, взаимодействующей с основанием // Строительная механика инженерных конструкций и сооружений. 2019. Т. 15. № 4. С. 291-298.
15. Мануйлов Г.А., Косицын С.Б., Грудцына И.Е. Численный анализ устойчивости подкрепленных пластин с некратными критическими нагрузками // Строительная механика инженерных конструкций и сооружений. 2020. Т. 16. № 1. С. 54-61.
RESEARCH PAPER
The stability of geometrically nonlinear plate systems under the action of dynamic loads
Sergey P. Ivanov1'2*, Anastasia S. Ivanova1, Oleg G. Ivanov1
1 Volga State University of Technology, 3 Lenin Sq, Yoshkar-Ola, 424000, Russian Federation 2Mari State University, 1 Lenin Sq, Yoshkar-Ola, 424000, Russian Federation *sp-ivanov@mail.ru
Article history: Abstract
Received: January 28, 2020 Relevance. Single-connected and multi-connected plate systems are widely
Revised: March 04, 2020 used in construction, aircraft, shipbuilding, mechanical engineering, instrument
Accepted: March 24, 2020 making. As a result, the study of the stability of geometrically nonlinear plate
Sergey P. Ivanov, Doctor of Science, Professor, Head of the Department of Strength of Materials and Applied Mechanics of VSUT; Professor of the Department of Electromechanics of MarSU; eLIBRARY SPIN-code: 5963-6739.
AnastasiaS. Ivanova, senior lecturer, Department of Strength of Materials and Applied Mechanics; eLIBRARY SPIN-code: 9568-3451.
Oleg G. Ivanov, Cand. Sc., Associate Professor, Department of Strength of Materials and Applied Mechanics; eLIBRARY SPIN-code: 5052-9077.
For citation
Ivanov S.P., Ivanova A.S., Ivanov O.G. The stability of geometrically nonlinear plate systems under the action of dynamic loads. Structural Mechanics of Engineering Constructions and Buildings. 2020;16(3):219-225. http://dx.doi.org/10.22363/1815-5235-2020-16-3-219-225 (In Russ.)
systems is an urgent topic. But, despite significant achievements in this area, there are still many unsolved problems. Thus, the requests of the above-mentioned areas of application of thin-walled spatial systems require further study of the issue of static and dynamic stability. The aim of the work - development of a method of the dynamic stability analysis of geometrically nonlinear plate systems such as prismatic shells under the action of dynamic compression loads. Methods. A plate system, which is subject to dynamic compression loads in the longitudinal direction, is considered. Kirchhoff - Love hypotheses are taken into account. The material stress-deformation diagram is linear. The displacement of points in the normal direction to the median plane of the plates is determined in the form of the Vlasov expansion. To derive the basic differential equations of stability, the energy method and the variational Vlasov method are used. The extreme value of the total energy is determined using the Euler - Lagrange equation. As a result, a set of basic nonlinear differential equations for studying the buckling of the plate system under the action of dynamic compression loads is obtained. Results. The developed method is used to stability analysis of a geometrically nonlinear prismatic shell with a closed contour of the cross section, under central compression under the action of dynamic loading. The edges of the shell rest on the diaphragm. The buckling of the prismatic shell in the longitudinal direction along one and two half-waves of a sinusoid is studied. The numerical integration of nonlinear differential equations is performed by the Runge - Kutta method. Based on the calculation results, graphs of the dependence of the relative deflection on the dynamic coefficient are constructed. The influence of the rate of change of compression stress, the initial imperfection of the system, and other parameters on the criteria for the dynamic stability of the plate system is investigated.
Keywords: dynamic stability, geometric nonlinearity, plate system, prismatic shells, variation method of Vlasov
References
1. Ivanov S.P., Ivanova A.S. Prilozheniye variacion-nogo metoda V.Z. Vlasov a k resheniyu nelinejnykh zadach plastinchatykh system [Application of V.Z. Vlasov's varia-tional method to solving nonlinear problems of plate systems]. Yoshkar-Ola: PGTU Publ.; 2015. (In Russ.)
2. Vlasov V.Z. Tonkostennye prostranstvennye sistemy [Thin-Walled spatial systems]. Moscow: Gosstrojizdat Publ.; 1958. (In Russ.)
3. Ivanov S.P., Ivanova A.S. The dynamic stability of physically nonlinear plate systems. Structural Mechanics of Engineering Constructions and Buildings. 2014;(4):11-20. (In Russ.)
4. Ivanov S.P., Ivanov O.G., Ivanova A.S. The dynamic stability of physically nonlinear plate systems under biaxial compression. Structural Mechanics of Engineering Constructions and Buildings. 2018;(2):132-141. (In Russ.)
5. Volmir A.S. Ustojchivost' deformiruemyh sistem [Stability of deformable systems]. Moscow: Nauka Publ.; 1967. (In Russ.)
6. Volmir A.S. Ustojchivost' deformiruemyh sistem [Nonlinear dynamic of plats and shells]. Moscow: Nauka Publ.; 1972. (In Russ.)
7. Khamitov T.K., Fatykhova R.R. On stability of elastic-plastic cylindrical shell under longitudinal impact. News of the KSUAE. 2016;(4):490-496. (In Russ.)
8. Trushin S.I., Sysoeva E.V., Zhuravleva T.A. The stability of nonlinear deformable cylindrical composite shells under non-uniform loads. Structural Mechanics of Engineering Constructions and Buildings. 2013;(2):3-10. (In Russ.)
9. Trushin S.I., Zhuravleva T.A., Sysoeva E.V. Dynamic buckling of nonlinearly deformable reticulate plates from composite material with different lattice configurations. Nauch-noe obozrenie [Scientific review]. 2016;(4):44-51. (In Russ.)
10. Vescovini R., Dozio L. Exact refined buckling solutions for laminated plates under uniaxial and biaxial loads. Composite Structures. 2015;(12):356-368.
11. Nazarimofrad E., Barkhordar A. Buckling analysis of orthotropic rectangular plate resting on Pasternak elastic foundation under biaxial in-plane loading. Mechanics of Advanced Materials and Structures. 2016;23(10):1144-1148.
12. Ruocco E., Reddy J.N. A closed-form solution for buckling analysis of orthotopic Reddy plates and prismatic plate structures. Composites Part B: Engineering. 2019;(169): 258-273.
13. Lukash, P.A. Osnovy nelinejnoj stroitel'noj mek-haniki [Fundamentals of nonlinear structural mechanics]. Moscow: Strojizdat Publ.; 1978. (In Russ.)
14. Kosytsyn S.B., Akulich V.Yu. The definition of the critical buckling load beam model and two-dimensional model of the round and two-dimensional model of the round cylindrical shell that interact with the soil. Structural Mechanics of Engineering Constructions and Buildings. 2019; 15(4):291-298. http://dx.doi.org/10.22363/1815-5235-2019-15-4-291-298 (In Russ.)
15. Manuylov G.A., Kositsyn S.B., Grudtsyna I.E. Numerical analysis of stability of the stiffened plates subjected aliquant critical loads. Structural Mechanics of Engineering Constructions and Buildings. 2020;16(1):54-61. http://dx. doi.org/10.22363/1815-5235-2020-16-1-54-61 (In Russ.)