CC BY
16
Строительная механика инженерных конструкций и сооружений Structural Mechanics of Engineering Constructions and Buildings
2018. 14 (2). 132-141
http://journals.rudn.ru/ strUctural-mechanics
Расчеты на устойчивость
УДК 539.376
DOI: 10.22363/1815-5235-2018-14-2-132-141
ДИНАМИЧЕСКАЯ УСТОЙЧИВОСТЬ ФИЗИЧЕСКИ НЕЛИНЕЙНЫХ ПЛАСТИНЧАТЫХ СИСТЕМ ПРИ СЖАТИИ В ДВУХ НАПРАВЛЕНИЯХ
СП. ИВАНОВ*;**, О.Г. ИВАНОВ*;**, АС. ИВАНОВА*
*Поволжский государственный технологический университет
пл. Ленина, 3, Йошкар-Ола, Россия, 424000 ** Марийский государственный университет пл. Ленина, 1, Йошкар-Ола, Россия, 424000
(поступила в редакцию: 10 сентября 2017 г.; принята к публикации: 05 марта 2018 г.)
В работе представлен метод расчета на динамическую устойчивость пластинчатых систем с несмещающимися ребрами. Рассмотрена пластинчатая система, на которую в двух взаимно перпендикулярных направлениях действуют динамические сжимающие нагрузки. В основе расчетов учитываются гипотезы Кирхгофа - Лява, гипотеза о нелинейно-упругом теле. Материал пластинчатой системы принимается физически нелинейным, диаграмма деформирования аппроксимируется в виде кубического полинома. Перемещение точек в нормальном направлении к срединной плоскости пластин представлено в виде разложения по Власову. Для вывода основных дифференциальных уравнений устойчивости используется энергетический метод и вариационный метод Власова. Экстремальное значение полной энергии системы определяется с использованием уравнения Эйлера - Лагранжа, после раскрытия которого получена система основных нелинейных дифференциальных уравнений для исследования потери устойчивости пластинчатой системы с несмещающимися ребрами под действием динамических сжимающих нагрузок. В качестве примера выполнен расчет на устойчивость физически нелинейной Т-образной пластинчатой системы, края которой закреплены шарнирно по контуру. Потеря устойчивости пластинчатой системы в продольном направлении происходит по одной полуволне синусоиды. При решении задачи в первом приближении выведено нелинейное дифференциальное уравнение, численное интегрирование которого проводилось методом Рунге -Кутта. По результатам расчетов построены графики зависимости относительной величины прогиба от динамического коэффициента. Исследовано влияние на динамический критерий потери устойчивости пластинчатой системы степени физической нелинейности материала, скорости изменения динамической сжимающей нагрузки и других параметров.
Ключевые слова: динамическая устойчивость, физическая нелинейность, пластина, пластинчатая система, сжимающая нагрузка, энергетический метод, вариационный метод Власова
Исследование устойчивости пространственных пластинчатых систем (тонкостенных конструкций) - одна из важных проблем строительной механики. К настоящему моменту в области расчетов на устойчивость тонкостенных конструкций накоплен значительный опыт: построено большое число математических моделей деформирования пластин и оболочек, разработаны эффективные аналитические и численные методы решения задач. Но, несмотря на достиже-
ния в теории и практике, имеется еще много нерешенных вопросов. Проблеме расчета на устойчивость пластин и оболочек при статических и динамических воздействиях посвящены работы многих авторов [1-7].
Цель настоящей работы заключается в разработке метода расчета на устойчивость пластинчатых систем с несмещающимися ребрами под действием динамических сжимающих нагрузок в двух взаимно перпендикулярных направлениях (рассматривается общий случай действия сжимающих нагрузок) с учетом физической нелинейности материала.
Рассмотрим пластинчатую систему (рис. 1), на которую в двух взаимно перпендикулярных направлениях действуют сжимающие нагрузки P(t) и aP(t). Динамическая нагрузка P(t) прикладывается к ребрам пластинчатой системы в поперечном направлении, aP(t) - в продольном направлении, а = 1, 2,..., п.
Принимаем, что сжимающая нагрузка быстро изменяется во времени t по следующему закону:
P(t) = s ■ t , (1)
где s - величина, характеризующая скорость изменения динамической нагрузки.
Обозначаем: w = w(x, y, t) - перемещение точек в нормальном направлении к срединной плоскости пластин.
Рис. 1. Общая схема пластинчатой системы с действующими нагрузками [Fig. 1. General diagram of plate system under the action of loads]
Результаты испытаний различных конструкционных материалов (композитов, сплавов, цветных металлов) показывают, что для таких материалов характерна физическая нелинейность, т.е. они имеют нелинейную диаграмму деформирования. Зависимость между интенсивностями напряжений а, и интенсивно-стями деформаций e, можно принять в виде кубического полинома:
at = E■ei-El e3, (2)
где Е - начальный модуль упругости материала; Ej - постоянная, учитывающая степень физической нелинейности материала [8].
В основе расчетов используем гипотезы Кирхгофа - Лява и гипотезу о нелинейно-упругом теле.
Принимаем соотношения между деформациями и перемещениями в виде
sx ~ —x; Sy _ —y; Sxy ~ —2zXxy, (3)
где X* =
ô2w
X ôx2'
X y
ô2w
У Я.2'
ôy
X xy
ô2w ôxôy
Определяем интенсивность деформаций е^ и объемную деформацию 0 (с учетом гипотез Кирхгофа - Лява и сжимаемости материала (стг = 0, еХ2 = 0, еу2 = 0)):
Л
2 2 2 3 2
2(1 + vW(s * "Sy) + (s у -s z ) + (s z -s *) + 2 S ХУ
где деформация
1 - 2v
0=--(s x + s y ),
1-v
(s x +s y )•
1-v
(4)
(5)
(6)
Используем энергетический метод. Записываем выражение полной энергии системы:
Ь = П + Т,
где потенциальную энергию П определяем по формуле
п = Д
л - 2 pU )
2 ^ ô w
vôx2 у
—aP(t ) 2
2 ^ ôw
2
ôy
+ q • w
V "s ;
dxdy,
кинетическая энергия равна:
Т=1 Г*
2" g
ôw ~dt
dxdy.
(7)
(8)
(9)
В выражениях (8), (9) введены следующие обозначения: р - объемный вес материала; 5 - толщина элементов, составляющих пластинчатую систему; g -ускорение свободного падения; Р(?) и аР(?) - динамические сжимающие нагрузки; q - интенсивность нагрузки, которая действует в нормальном направлении к срединной плоскости пластин и позволяет учитывать начальное несовершенство системы; А - работа внутренних сил, отнесенная к единице площади поверхности пластинчатой системы, определяющаяся так:
5/2
А = {Фdz. (10)
-5/2
Здесь Ф - удельная энергия изменения объема и формы [9]:
, e;
1
2
Ф = -К-02 + - J(1 + v)• а,- • det
3
(11)
0
где К = Е/[3(1 - 2у)] - модуль объемного сжатия; V - коэффициент Пуассона. Перемещение V представим в виде разложений по В.З. Власову [10]:
Vх,О = £Жк (?)/к (х,^), (к = 1,2,3,...,и), (12)
к
где обобщенные перемещения Жк(() - функции, которые определяются из решения задачи и зависят от переменной ? (? - время). Выбор координатных функций /к(х,у) осуществляется по виду деформированного состояния системы.
v
2
2
Учитывая соотношение (12), определяем минимум функционала (7), составляя для него уравнение Лагранжа:
cL d CL
W dt dWk,
■ = 0,
(13)
где индексы после запятой обозначают частные производные от обобщенных перемещений по соответствующим переменным.
После раскрытия уравнения (13) приходим к системе нелинейных дифференциальных уравнений для исследования устойчивости пластинчатых систем (под действием динамических сжимающих нагрузок и оР(0):
I
k
fl aP(f)) 2b + fl Р(1) aik \ 1--^ I -2bik + cik\ 1"
* „ a D
*
l D
W + IdikWk,tt -Gi = ФГ, (14) g •D k
где а - длина контура пластинчатой системы, на который действует нагрузка аР(0; l - суммарная длина ребер пластинчатой системы, на которые действует нагрузка Р^); D - цилиндрическая жесткость пластин системы; Gi - величина нагрузки, которая позволяет учитывать начальное несовершенство системы.
Величины D и Gi определяются по формулам:
Е5 ? ^=яqfdxdy. (15)
D =■
12(1 -v2)
x у
(16)
Правая часть уравнений (14) Ф у учитывает физическую нелинейность материала и имеет вид:
Фнел. = Ц^1,хх / dxdy +| 1^2 /у,уу dxdy -{{N3,х/ dxdy, ху ху ху
где = /дх; ^^ = д2N /дх2,...;. у = 1, 2,3.
Величины N1, N2, N3 в выражении (16) определяются по формулам: N1 = Д • с • (У1Х х + 0,5У2Х у);
N2 = Д • с • (У1Х у + 0,5У2Х х );
N3 = Д •с ху ;
(17)
где Д = 0,9E2 •б2
е2 = E •(1 v2; c=v1(xi+i2y)+v2XxX y
E •(! + v)
V1 =-
1 3
(1 + v)2
■+1
v2 =-2 3
2v
(1 -v)2
--1
Коэффициенты уравнения (14) записываются в виде:
aik = J \fk,хх fi,хх dx dy; bik = J \fkxxyfi
yx
cik JJ fk, yyfi
yx
УУ
ik ¡JJ k ,xyJi, xy у x
dik = JJfkfi dxdy.
(18)
yx
1
1
v
В качестве примера выполним расчет на устойчивость Т-образной пластинчатой системы под действием динамических сжимающих нагрузок P(t) и aP(t) в двух взаимно перпендикулярных направлениях (рис. 2). Края пластинчатой системы закреплены шарнирно по контуру. Геометрические параметры: толщина пластин 5 = 0,1а; длина l = a.
Рис. 2. Схема Т-образной пластинчатой системы [Fig. 2. The diagram of T-shaped plate system]
При решении задачи в первом приближении, согласно граничным условиям, координатная функция (форма потери устойчивости в поперечном направлении х) (рис. 3) имеет вид
ft \ x 1т 3x 1 x
f1(x) =-| -1 +-----:
a I 2 a 2 a
(19)
Рис. 3. Вид координатной функции [Fig. 3. Formula for a coordinate function]
При потере устойчивости в продольном направлении у по одной полуволне синусоиды прогибы определяются по формуле
w(x,y,t) = Wi(t)fi(x,y) = Wi(t)—| -1 ^ Isi^.
3 x 1 x2 | • ny
2 a 2 a
l
(20)
Тогда дифференциальное уравнение (14) принимает следующий вид
fi aP(t ) J + fi P(t)
Ч1 -OD J-2b11 +M1 - JD
•W1 + duWUt - G1 = 0^. (21) g • D
a
Вводим обозначение для параметра времени:
P(t) st
*
t =■
PP
1 кр 1 кр
(22)
Для пластинчатой системы, составленной из квадратных пластин (I = а), величина погонной статической критической нагрузки составляет [11]
P
Ркр.
кр.
= 23,5-
D
.2 '
(23)
Параметр времени t принимаем за динамический коэффициент Кд. После преобразований дифференциальное уравнение (21) запишем с новой переменной t в виде
1 - К1(1 + а) • t
1 3
•W +—г-W * * = Q + K•Wï3, 1 s * l,t t 1 '
(24)
где приняты следующие обозначения:
(a11 + c11) • Ркр .
К =
а • (ац - 2Ьц + Сц) • D '
K =
Фнел. 1
S* = (а11 - 2Ь11 + С11) • g • D • Ркр .
р- d
Q =
11
а11- 2Ь11 + С11 ' ^ .
а11 - 2Ь11 + Сп'
(25)
Величину Фнел. определяем, используя формулы (16) и (17), а коэффициенты уравнения (21) вычисляем по формулам (18).
Интегрирование дифференциального уравнения (24) выполнено численным методом Рунге - Кутта на ПЭВМ.
По результатам расчетов построены графики зависимости относительной величины прогиба = Щ/ 5 от динамического коэффициента Кд = {* (при различных значениях параметров Е/Е, 5*, а) (рис. 4, табл. 1). Для всех графиков параметр Q = 0,0001.
Рис. 4. Графики зависимости относительной величины прогиба W от динамического коэффициента Кд [Fig. 4. Graphs of the relative magnitude of deflection W* against the dynamic coefficient Кд]
a
Таблица 1
Номер графика [Graph number] Степень физической нелинейности [Degree of physical nonlinearity] E1/E Коэффициент [Coefficient] a Параметр [Parameter] S" Динамический коэффициент Кд при W = 0,1 [Dynamic coefficient Кд at W" = 0,1]
1 5-106 2 10 « 2,7
2 106 1 10 « 3,0
3 103 1 1,0 « 3,4
4 106 1 0,1 « 6,5
5 106 0 0,1 « 15,0
Рассмотрим влияние на критерии динамической устойчивости пластинчатой системы степени физической нелинейности материала, параметра S (зависит от скорости изменения сжимающей нагрузки), величины коэффициента а.
С увеличением степени физической нелинейности материала график зависимости W от Кд смещается влево, т.е. динамический коэффициент Кд уменьшается (см. графики 3, 2, 1 на рис. 4). Таким образом, чем больше степень физической нелинейности материала, тем меньше динамическая «критическая» нагрузка, т.е. при меньшем значении динамической нагрузки P(t) происходит бурное выпучивание пластин системы.
При уменьшении параметра S от 10 до 0,1 (что соответствует увеличению скорости изменения динамической нагрузки в 10 раз при прочих равных условиях) динамический коэффициент Кд увеличивается от 3,0 до 6,5 при W = 0,1 (сравниваем графики 2 и 4 на рис. 4).
Если сжимающая динамическая нагрузка F(t) действует только в одном направлении (при а = 0, см. график 5), то динамический коэффициент принимает в несколько раз большее значение (Кд ~ 15 при W = 0,1), чем при действии F(t) и aF(t) в двух направлениях (при прочих равных условиях, см. график 4).
Выводы. Разработанный метод позволяет рассчитывать на устойчивость пластинчатой системы с несмещающимися ребрами под действием динамических сжимающих нагрузок в двух направлениях, учитывая физическую нелинейность материала. Из результатов проведенных исследований видно, что на критерии динамической устойчивости пластинчатой системы значительное влияние оказывают скорость изменения динамической сжимающей нагрузки и степень физической нелинейности материала.
© Иванов С.П., Иванов О.Г., Иванова А. С., 2018
Список литературы
1. ВольмирА.С. Устойчивость деформируемых систем. М.: Наука, 1967. 984 с.
2. Трушин С.И., Журавлева Т.А., Сысоева Е.В. Динамическая потеря устойчивости нелинейно деформируемых сетчатых пластин из композиционного материала с различными конфигурациями решетки // Научное обозрение. 2016. № 4. С. 44-51.
3. Vescovini R., Dozio L. Exact refined buckling solutions for laminated plates under uniaxial and biaxial loads // Composite Structures. 2015. Vol. 127. P. 356-368.
4. 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. № 10. P. 1144-1148.
5. Хамитов Т.К., Фатыхова Р.Р. Об устойчивости упругопластической цилиндрической оболочки при продольном ударе // Известия Казанского государственного архитектурно-строительного университета. 2016. № 4 (38). С. 490-496.
6. Cao G., Chen Z., Yang L., Fan H., Zhou F. Analytical Study on the Buckling of Cylindrical Shells With Arbitrary Thickness Imperfections Under Axial Compression // Journal of Pressure Vessel Technology Transactions of the ASME. 2014. Vol 137. № 1. DOI: 10.1115/ 1.4027179.
7. Трушин С.И., Журавлева Т.А., Сысоева Е.В. Устойчивость нелинейно деформируемых цилиндрических оболочек из композиционного материала при действии неравномерных нагрузок // Строительная механика инженерных конструкций и сооружений. 2013. № 2. С. 3-10.
8. Иванов С.П., Иванова А.С. Приложение вариационного метода В.З. Власова к решению нелинейных задач пластинчатых систем: монография. Йошкар-Ола: ПГТУ, 2015. 248 с.
9. ЛукашПА. Основы нелинейной строительной механики. М.: Стройиздат, 1978. 204 c.
10. Власов В.З. Тонкостенные пространственные системы. M.: Госстройиздат, 1958. 502 c.
11. Иванов С.П., Иванова А.С. Статическая устойчивость пластин и пластинчатых систем с несмещающимися ребрами при сжатии в двух направлениях // Строительная механика инженерных конструкций и сооружений. 2015. № 4. С. 63-67.
Об авторах:
Иванов Сергей Павлович - доктор технических наук, профессор, заведующий кафедрой сопротивления материалов и прикладной механики, Поволжский государственный технологический университет; профессор кафедры электромеханики, Марийский государственный университет. Опубликовал 147 научных статей, 2 монографии, 4 учебника, 20 наименований учебно-методической литературы. Область научных интересов: расчеты на прочность, устойчивость и колебания физически и геометрически нелинейных стержней, пластин и пластинчатых систем. Контактная информация: e-mail - [email protected]
Иванов Олег Геннадьевич - кандидат технических наук, доцент, доцент кафедры сопротивления материалов и прикладной механики, Поволжский государственный технологический университет. Опубликовал 35 научных статей, 1 монографию и 6 наименований учебно-методической литературы. Область научных интересов: расчеты на прочность и устойчивость физически нелинейных пластин и пластинчатых систем, контактирующих с упругой средой. Контактная информация: e-mail - [email protected]
Иванова Анастасия Сергеевна - аспирант, старший преподаватель кафедры сопротивления материалов и прикладной механики, Поволжский государственный технологический университет. Научный руководитель - д.т.н., проф. С.П. Иванов, Поволжский государственный технологический университет. В настоящее время работает над кандидатской диссертацией «Динамическая устойчивость физически нелинейных стержней, пластин и пластинчатых систем» по специальности 05.23.17 - Строительная механика. Опубликовала 15 научных статей и 1 учебное пособие. Область научных интересов: расчеты на устойчивость физически нелинейных стержней, пластин и пластинчатых систем. Контактная информация: e-mail - [email protected]
Для цитирования:
Иванов С.П., Иванов О.Г., Иванова А.С. Динамическая устойчивость физически нелинейных пластинчатых систем при сжатии в двух направлениях // Строительная механика инженерных конструкций и сооружений. 2018. Т. 14. № 2. С. 132-141. DOI: 10.22363/1815-5235-2018-14-2-132-141.
References
1. Volmir, A.S. (1967). Stability of deformable systems, Moscow: Science, 984 p. (In Russ.)
2. Trushin, S.I., Zhuravleva, T.A., Sysoeva, E.V. (2016). Dynamic buckling of nonlinear-ly deformable reticulate plates from composite material with different lattice configurations, Nauchnoe Obozrenie, No 4, p. 44-51. (In Russ.)
РАСЧЕТЫ НА УСТОЙЧИВОСТЬ
139
3. Vescovini, R., Dozio, L. (2015). Exact refined buckling solutions for laminated plates under uniaxial and biaxial loads, Composite Structures, Vol. 127, p. 356-368.
4. Nazarimofrad, E., Barkhordar, A. (2016). Buckling analysis of orthotopic rectangular plate resting on Pasternak elastic foundation under biaxial in-plane loading, Mechanics of Advanced Materials and Structures, Vol. 23, No 10, p. 1144-1148.
5. Khamitov, T.K., Fatykhova, R.R. (2016). On stability of elastic-plastic cylindrical shell under longitudinal impact, News of the Kazan State University of Architecture and Engineering, No 4 (38), p. 490-496. (In Russ.)
6. Cao, G., Chen, Z., Yang, L., Fan, H., Zhou, F. (2015). Analytical study on the buckling of cylindrical shells with arbitrary thickness imperfections under axial compression, Journal of Pressure Vessel Technology Transactions of the ASME, Vol. 137, No 1.
7. Trushin, S.I., Zhuravleva, T.A., Sysoeva, E.V. (2013). The stability of nonlinear deformable cylindrical composite shells under non-uniform loads, Structural Mechanics of Engineering Constructions and Buildings, № 2, p. 3-10. (In Russ.)
8. Ivanov, S.P., Ivanova, A.S. (2015). Prilozheniye variacionnogo metoda V.Z. Vlasova k resheniyu nelinejnykh zadach plastinchatykh system: Monographiya, Yoshkar-Ola: PGTU, 248 p. (In Russ.)
9. Lukash, P.A. (1978). Osnovy nelinejnoj stroitel'noj mekhaniki, Moscow: Strojizdat, 204 p. (In Russ.)
10. Vlasov, V.Z. (1958). Tonkostennye prostranstvennye sistemy, Moscow: Gosstroj-izdat, 502 p. (In Russ.)
11. Ivanov, S.P., Ivanova, A.S. (2015). The static stability of plates and plate systems with nonshifting ribs under biaxial compression, Structural Mechanics of Engineering Constructions and Buildings, No 4, p. 63-67. (In Russ.)
THE DYNAMIC STABILITY OF PHYSICALLY NONLINEAR PLATE SYSTEMS UNDER BIAXIAL COMPRESSION
S.P. IVANOV*;**, O.G. IVANOV*;**, A.S. IVANOVA*
*Volga State University of Technology
3 Lenin Sq., Yoshkar-Ola, 424000, Russian Federation
**Mari State University
1 Lenin Sq., Yoshkar-Ola, 424000, Russian Federation
(received: September 10, 2017; accepted: March 05, 2018)
The article presents the method of dynamic stability analysis of plate systems with non-shifting ribs. A plate system under the biaxial dynamic compression loads is considered. The Kirchhoff - Love hypotheses, the nonlinear-elastic body hypothesis are considered the basis of the calculations. The material of the plate system is assumed to be physically nonlinear, stress-deformation diagram is approximated in the form of a cubic polynomial. The displacement of points in normal direction to middle plane of plates is presented in the form of Vlasov expansion. To derive the basic differential equations of stability, the strain-energy method and Vlasov's variation method are used. The extreme value of total energy of the system is defined using Euler - Lagrange equation, after solving of which the set of basic nonlinear differential equations of buckling of the plate system with non-shifting ribs under dynamic compression loads is given. As an example, the stability calculation of physically nonlinear T-shaped plate system hinge-supported along the contour is carried out. Buckling of the plate system occurs longitudinally on one half-wave of sinusoid. At the solution of a task in the first approximation, a nonlinear differential equation is derived, the numerical integration of which was carried out by the Runge - Kutta method. Based on the results of the calcu-
lations, graphs of the relative magnitude of deflection against the dynamic coefficient are plotted. The influence of the degree of physical nonlinearity of the material, the rate of change of the dynamic compressive load on the dynamic criterion of buckling of the plate system was studied.
Keywords: dynamic stability, physical nonlinearity, plate, plate system, compression load, strain-energy method, variation method of Vlasov
About the authors:
Ivanov Sergey Pavlovich - Doctor of Science, Professor, Head of Department of Strength of Materials and Applied Mechanics, the Volga State University of Technology; Professor of Department of Electro-Mechanics, the Mari State University. He is the author of 147 scientific articles, 2 monographs, 4 textbooks, 20 names of educational literature. General research interests: strength, stability and vibrations analyses of the physically and geometric nonlinear rods, plates and plate systems. Contact information: 3 Lenin Sq., Yoshkar-Ola, 424000, the Republic of Mari El, Russia, e-mail - [email protected]
Ivanov Oleg Gennadevich - Cand. Sc, Assistant Professor, Associate Professor of the Department of Strength of Materials and Applied Mechanics, the Volga State University of Technology. He is the author of 35 scientific articles, 1 monograph, 6 names of educational literature. General research interests: strength and stability analyses of the physically nonlinear plates and plate systems resting on elastic foundation. Contact information: 3 Lenin Sq., Yoshkar-Ola, 424000, the Republic of Mari El, Russia, e-mail - [email protected]
Ivanova Anastasia Sergeevna - Aspirant, Senior Lecturer of the Department of Strength of Materials and Applied Mechanics, the Volga State University of Technology. At the present time she works on the Candidate's dissertation «The dynamic stability of physically nonlinear rods, plates and plate systems» in the specialty 05.23.17 - Structural Mechanics. She is the author of 15 scientific articles, 1 name of educational literature. General research interests: stability analyses of the physically nonlinear rods, plates and plate systems. Contact information: e-mail - [email protected]
For citation:
Ivanov, S.P., Ivanov, O.G., Ivanova, A.S. (2018). The dynamic stability of physically nonlinear plate systems under biaxial compression. Structural Mechanics of Engineering Constructions and Buildings. 14 (2), 132-141. DOI: 10.22363/1815-5235-2018-14-2-132-141.
