Math-Net.Ru
А. П. Янковский, Моделирование упругопла-стического поведения гибких пространственно-армированных пластин в рамках уточненной теории изгиба, Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2019, том 23, номер 1, 90-112
001: https://doi.org/10.14498/vsgtu1644
Использование Общероссийского математического портала Math-Net.Ru подразумевает, что вы прочитали и согласны с пользовательским соглашением
http://www.mathnet.ru/rus/agreement
Параметры загрузки:
IP: 109.252.26.93
29 ноября 2019 г., 20:53:20
Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2019. Т. 23, № 1. С. 90—112 ISSN: 2310-7081 (online), 1991-8615 (print) d https://doi.org/10
УДК 539.4
Моделирование упругопластического поведения гибких пространственно-армированных пластин в рамках уточненной теории изгиба
А. П. Янковский
Институт теоретической и прикладной механики им. С. А. Христиановича СО РАН,
Россия, 630090, Новосибирск, ул. Институтская, 4/1.
Аннотация
На базе алгоритма шагов по времени построена структурная модель упругопластического деформирования изгибаемых пластин с пространственными структурами армирования. Неупругое поведение материалов фаз композиции описывается определяющими уравнениями теории пластического течения с изотропным упрочнением. Возможное ослабленное сопротивление армированных пластин поперечному сдвигу учитывается на основе уточненной теории, из которой в первом приближении получаются соотношения теории Редди. Геометрическая нелинейность задачи рассматривается в приближении Кармана. Решение сформулированных начально-краевых задач строится по явной численной схеме типа «крест». Исследовано динамическое неупругое деформирование пространственно- и плоско-перекрестно армированных металло-композитных и стеклопластиковых гибких пластин разной относительной толщины под действием нагрузки, вызванной воздушной взрывной волной. Продемонстрировано, что для относительно толстых стеклопла-стиковых пластин замена плоско-перекрестной структуры армирования на пространственную структуру с сохранением общего расхода волокон приводит к уменьшению податливости конструкции в поперечном направлении почти в 1.5 раза, а также к уменьшению максимума интенсивности деформаций в связующем материале в два раза. Для относительно тонких как стеклопластиковых, так и металлокомпозитных пластин замена плоско-перекрестной структуры 2Б-армирования на пространственные структуры 3D- и 4Б-армирования не приводит к заметному уменьшению их прогибов, но позволяет уменьшить интенсивность деформаций в связующем на 10 % и более. Показано, что широко используемая неклассическая теория Редди не позволяет получать надежные результаты расчетов упругопластического динамического поведения изгибаемых пластин как с плоскими, так и пространственными структурами армирования даже при малой относительной толщине конструкций и слабой анизотропии композиции.
Научная статья
3 ©® Контент публикуется на условиях лицензии Creative Commons Attribution 4.0 International (https://creativecommons.org/licenses/by/4.0/deed.ru) Образец для цитирования
ЯнковскийА. П. Моделирование упругопластического поведения гибких пространственно-армированных пластин в рамках уточненной теории изгиба // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2019. Т. 23, № 1. С. 90-112. doi: 10.14498/vsgtu1644. Сведения об авторе
Андрей Петрович Янковский & https://orcid.org/0000- 0002- 2602- 8357
доктор физико-математических наук; ведущий научный сотрудник; лаб. физики быстро-
протекающих процессов; e-mail: lab4nemir@rambler.ru
.14498/vsgtu1644
Ключевые слова: композитные пластины, пространственно-перекрестное армирование, плоско-перекрестное армирование, теория Редци, уточненная теория изгиба, упругопластическое деформирование, геометрическая нелинейность, взрывные нагрузки, схема типа «крест».
Получение: 18 сентября 2018 г. / Исправление: 11 февраля 2019 г. / Принятие: 4 марта 2019 г. / Публикация онлайн: 20 марта 2019 г.
Введение. Конструкции из композиционных материалов (КМ) находят широкое применение в современной инженерной практике [1-8,10]. При этом все более активно используются пространственно-армированные материалы [6-8,10], так как такой тип армирования упраздняет серьезный недостаток плоско-армированных сред — расслоение в силу слабого сопротивления поперечным сдвигам и отрыву, а также позволяет локализовать трещины в пределах нескольких ячеек периодичности. Упругое деформирование волокнистых сред с пространственными структурами армирования моделировалось в [9,11,12]. Однако КМ-изделия часто испытывают высокоинтенсивное внешнее воздействие [1-3, 6, 13]. При этом КМ-конструкции могут деформироваться упругопластически. В связи с этим в [14] была построена структурная модель упругопластического деформирования пластин, перекрестно армированных в своей плоскости. Структурная же модель неупругого поведения КМ-среды с пространственной структурой армирования, использующая определяющие уравнения теории пластического течения, пока не построена.
Так как при пространственном ортогональном армировании [6-10] тонкостенные элементы КМ-конструкций типа пластин также могут обладать ослабленным сопротивлением поперечным сдвигам, актуальной является проблема не только математического моделирования упругопластического деформирования КМ-сред, но и адекватного описания механического поведения КМ-пластин с такими структурами армирования при изгибе с учетом их возможного ослабленного сопротивления поперечным сдвигам. Последнее обстоятельство традиционно учитывается в рамках теорий Рейсснера [13,15-18] и Редди [14,19-22]. В [19,23] продемонстрировано, что при линейно-упругом изгибном деформировании КМ-пластин уточнение теории Редди, которая в первом приближении учитывает искривление поперечной нормали, не требуется. Однако, согласно [23], при упругопластическом изгибном деформировании пластин, перекрестно армированных в своей плоскости, теория Редди не гарантирует получения приемлемых результатов расчетов их динамического поведения. Для адекватного описания ослабленного сопротивления таких КМ-пластин поперечным сдвигам следует использовать кинематические гипотезы более высокого порядка [23,24], чем это принято в теории Редди.
Согласно всему вышеизложенному, данное исследование посвящено математическому моделированию упругопластического поведения пространственно-армированных гибких пластин при использовании более точных, чем в теории Редди, кинематических гипотез. Для численного интегрирования возникающих при этом начально-краевых задач предполагается использовать явную пошаговую схему типа «крест».
1. Моделирование упругопластического деформирования пространственно-армированной гибкой пластины с учетом ослабленного сопротивления поперечным сдвигам. Рассматриваем изгибное деформирование КМ-пластины толщиной 2к (рис. 1), с которой свяжем прямоугольную декартову систему координат Хг так, чтобы отсчетная плоскость Х\Х2 (хэ = 0) совпадала со срединной плоскостью (|жэ| ^ к). Конструкция плоско-или пространственно усилена К семействами волокон с относительным их объемным содержанием в репрезентативной ячейке , 1 ^ к ^ К (интенсивности армирования). На рис. 1, а изображено ортогональное плоское 2Б-армирование при К = 2; на рис. 1, Ь — ортогональное пространственное 3Б-армирование при К = 3 [7,9,10]; на рис. 1, с — неортогональное пространственное 4Б-армирование при К = 4 [6]. В поперечном направлении хэ структура армирования считается квазиоднородной. Относительное объемное содержание связующего материала в представительном элементе композиции определяется так [6,10]:
к
= 1 Шк. (1)
к= 1
Для учета возможного (например в случаях ортогональных структур армирования, изображенных на рис. 1, а и Ь) ослабленного сопротивления КМ-пластины поперечным сдвигам используем кинематические соотношения более высокого порядка, чем принято в теории Редди. При этом предполагаем, что на лицевых поверхностях хэ = ±к отсутствуют внешние касательные
с
Рис. 1. Элементы КМ-пластин: a) с плоской ортогональной структурой 2Б-армиро-вания; b) с пространственной ортогональной структурой 3Б-армирования; с) с пространственной неортогональной структурой 4Б-армирования
[Figure 1. Elements of composite plates: a) with a flat-orthogonal 2D reinforcement structure; b) with a spatial-orthogonal 3D reinforcement structure; с) with a spatial non-orthogonal 4D reinforcement structure]
нагрузки. Тогда, согласно [23], средние деформации в композиции е^ и перемещения точек КМ-пластины ^ при учете геометрической нелинейности в приближении Кармана аппроксимируются так:
£гЗ &, г) = ^РгЩ + д^Уа) - Х^д^-Ш +
^ Хт+1 / Л2 Х2 \ / N , N 1
+£ V (^ - (*#>++(2)
т=0 (2)
^, г) = , *),
т=0
X = {Х1 ,Х2}, г = {Х1,Х2,Х3}, 3 = 1, 2;
^ Хт+1 / Л2 Х2 \ ( N иг(1, Г) = ^, Х) - Х3дг» + 2 £ ХЛг ( ^П - тХ+з) £*)(*, Ж), (3)
т=0
, г)=»(г, х), X еС, |хз| < Л, ¿0, г = 1,2,
где » — прогиб; у — тангенциальные перемещения точек отсчетной плоскости (хз = 0) в направлениях х%; ¿0 — начальный момент времени ¿; д% — оператор дифференцирования по хг (г = 1, 2); М — целое число, задающее количество слагаемых, которые удерживаются в степенных разложениях по Хз; С — область, занимаемая КМ-пластиной в плане. Последнее равенство (3) соответствует традиционной для тонких пластин кинематической гипотезе, согласно которой изменяемостью поперечного смещения £з(£, г) по переменной х3 можно пренебречь [13-23]. При М = 0 из равенств (2), (3) следуют кинематические допущения теории Редди [14,19-22].
Таким образом, в соотношениях (2), (3) неизвестными являются функции
», у и е(т) (г = 1, 2, 0 ^ т ^ М), которые зависят от двух пространственных переменных Х1 , Х2 и времени .
В силу того что установить действительное распределение полей деформаций и напряжений в КМ, в котором связующее содержит многочисленные армирующие включения достаточно произвольной ориентации (см. рис. 1), весьма затруднительно [21] (особенно при упругопластическом деформировании компонентов композиции), для построения практически приемлемых в инженерных приложениях определяющих уравнений неупругого деформирования КМ-среды используем исходные предпосылки, аналогичные принятым в [14,21].
1. В пределах представительного элемента КМ можно рассматривать как макроскопически анизотропное квазиоднородное тело. (При равномерном и достаточно густом наполнении связующего тонкими армирующими элементами это допущение вполне приемлемо [21].)
2. Между волокнами и связующим реализуется идеальный механический контакт.
3. В пределах репрезентативной ячейки, выделенной из КМ-среды на ми-ниуровне, напряжения и деформации во всех компонентах и в композиции в целом кусочно-постоянны. Эффекты более высоких порядков, характеризующие изменение напряжений и деформаций на микроуровне в окрестностях границ контакта волокон и связующего, не учитываем.
4. Поля деформаций и напряжений в КМ-среде осредняем по объему представительного элемента, т. е. согласно гипотезе 3 пропорционально Шк, 0 < к ^К (см. (1)).
5. Компоненты композиции однородны и изотропны, а их упругопластиче-ское деформирование описывается соотношениями теории Прандтля— Рейсса—Хилла (теории течения с изотропным упрочнением) [25,26]:
&к = Вк ¿к (В* = Ак - Рк), к = 0,1, 2,..., К, (4)
где
^ = га№ а(*0 а№\т_1а(к) а(к) а(к) а(к) *(*) л^т
^ = (а1 ,а2 ,л3 ,а4 ,а5 ,а6 ) = \а11 ,а22 ,лээ ,Л2Э ,ЛЭ1 , а12 ) , х ={+(к) .(к) .(к) .(к) .(к) -(к)Хт= г.(к) .(к) .(к) 2 .(к) 2 . (к) 2 . (ЛЬТ ьк= \ь 1 , ь2 , ьэ , ь 4 , ь 5 , ь 6 / = 1ь 11 , ь 22 , ь ээ , 2Ь2Э , 2ЬЭ1 , 2Ь12/ ;
1 2 э 4 5 6 11 22 ээ 2э э1 12 (5)
а,
(к) (к)
Ч
,
Ч
компоненты тензоров напряжений и деформаций в й-той фа-
зе композиции; В& = (Ь^), Ак = (а(^), Р& = (р(^) —симметричные
(к).
-г]
гз
г]
6х6-матрицы.
Ненулевые компоненты матриц Ак, Рк в (4) вычисляются по формулам [26]:
,(к)
= 2ц№ + \(к), а^ = а[к) = \(к),
(к)
А(к) =
г] зг_____
О = ЬЬЗ = 1, 3, г = 4, 6, 1,п = 1, 6); п(к)с(к) Е (к)
Ц С (1 -«(*)), 2ц№ = Е
ц(к), Р^ = А(к)зП) 4к)
X
(к)
1 + и(кУ
\(к) =
„(к) Е(к)
(1 + и(к))(1 - 2г/№)' э6
¿к)
' =
Ц
(к)
Ц
(кУ
X
(к)
„(к)
=
2
0 при
1 при
ж(й) =
г=1
г=1
X(к) < X
т № = т № Х2 = Х2*
г=4
(й) Л) Л)
или Х2 = ,
(&) (й)
4) = а)
э г=1
ж(А;) < 0, ж(й) > 0,
,(*) _ „со
(6)
'г
'г
ь з
еЛ ' =е\ '
ЭС0 _ „(*)
г
=
г=1
4? = шах{х£),х&)}, ¿ = М,
I = 4, 6;
Е(й) — коэффициент Пуассона и модуль Юнга й-того компонента компо-
касательный модуль на диаграмме чистого сдвига для того же (;г) —функция переключения, задающая активное упругопласти-
Лк)
зиции; ц материала;
ческое нагружение или разгрузку элемента /г-той фазы композиции; Х2р
- т(к)
значение второго инварианта тензора-девиатора напряжений , при котором й-той компонент композиции начинает впервые деформироваться пла-
стически; J^ — наибольшее значение , которое достигнуто за предыдущую историю нагружения элемента k-той фазы композиции; точка — дифференцирование по времени. Равенства (5) задают соответствия между шестью
элементами /(к) (i = 1, 6) некоторого вектор-столбца fk и компонентами симметричного тензора второго ранга /j? (j,l = 1, 3, 0 ^ k ^ К). В соответствии
с (6) s(к) и е(к) (i = 1, 6) — компоненты (возможно, удвоенные) девиаторов напряжений и деформаций в k-той фазе композиции.
С каждым k-тым семейством волокон свяжем локальную ортогональную
систему координат х(к) таким образом, чтобы х?к) совпадала с траекторией
армирования, а оси х(к), х3к) были перпендикулярны этим траекториям. При этом направление волокна k-того семейства можно задать углами сферической системы координат и p? (рис. 2). Направляющие косинусы ¿(k) между
х(к) и Xj (i,j = 1, 3) в этом случае определяются так: 4? = sin вк cos^fc, ® = sin вк sinpk, ® = cos вк,
= — sin^fc, = cos^fc, 4з) = 0, (7)
4? = — cos вк cos pk, tí? = — cos вк sin pk, 4? = sin вк, 1 ^ k ^ К.
При переходе от глобальной х^ к локальной (г,] = 1, 3) системе координат для векторов-столбцов, аналогичных (5), выполняются преобразования
& к = Gk&k £к = Qk^k
а,
( к) _
,(к) (к)
j=i
6
= Е
j=i
»j j
г = 1, 6,
где G? = (Щ»j)) и Qk = (g^-) — 6х6-матрицы, элементы которых вычисляют-
Hj
Ажз
Рис. 2. Локальная система координат, связанная с арматурой fc-того семейства [Figure 2. Local coordinate system associated with fc-th family of armature]
ся по формулам (см. (7))
(к) _ (к) _ ,(к),(к) (к) _ ( к) _ ,(к),( к)
5*11 _ <?11 _ ¿11 А11 , 5*12 _ ^12 _ '12 П2 , • • • ,
„(к) _ 9_(к) _ 2,(к),(к) 2 (к) _ (к) _ 2,(к),(к)
У16 _ 2416 _ 2' 12 '11 , • • • , 2У61 _ 461 _ 2' 21 '11
(9)
(к ) _ (к ) _ ,(к ),(к ) . ,(к ),(к ) 1 . , . К 5*66 _ <?66 _ П1 ¿22 + П2 (21 , 1 ^ К ^ К
Остальные элементы матриц О к и Qк, не выписанные в (9), приведены в [21, табл. 21.40 и 21.44]. Чертой сверху в (8) обозначены величины, заданные
в системе координат ж(к) (г _ 1, 3).
Согласно условиям сопряжения полей напряжений и перемещений на границах контакта связующего с волокнами, при учете второго и третьего допущений имеем (см. (8))
4к) _¿<ЙЧ0), £$Чк) _£$Ч0), ^1 < *<(10)
3 = 1 3 = 1 3 = 1 3 = 1
Здесь первое соотношение характеризует равенство удлинений арматуры того семейства и связующего материала вдоль волокна; остальные равенства— условия контакта между арматурой и связующей матрицей в напряжениях на боковой поверхности волокна.
Исходные предположения 1-5 аналогичны допущениям, которые ранее были приняты в [14], поэтому, повторяя рассуждения из [14] и учитывая соотношения (1), (4), (6) и (10), окончательно в матричном виде получим следующее определяющее уравнение упругопластического деформирования КМ-среды с пространственной структурой армирования:
& _ Бё, (11)
где
Б _ ( ШоБо + £ шкБкЕк ) Н, к =1
/ к л-1 (12)
Н _ шо1 +£шкЕк , 1 ;
к =1
Ек _ Б-1Ск, 1 О ;
&, £ — шестикомпонентные векторы-столбцы скоростей средних напряжений ац и деформаций в КМ, по структуре подобные (5); I — единичная 6x6-матрица; Б, Н, Е к, С к, Ок — 6х6-матрицы; Ю-1 — 6х6-матрица, обратная Ок.
Элементы с(к) и матриц Ск и Ок с учетом (4), (6), (7), (9), (10) имеют следующие выражения:
66 (к ) ,(к ) (к ) ( к) ^ ( кк(0) ,( к) ^ ( к) 7 (к) . 7• Т~а
1=1 1=1
При выводе соотношений (11) и (12), как и в [14], попутно получаются дополнительные матричные равенства
¿о = Не, ¿к = Еке0, 1 4 к 4К. (13)
Первое соотношение (13) определяет скорости деформаций в связующей матрице ¿о через скорости осредненных деформаций композиции ¿, а второе равенство выражает скорости деформаций волокон к-того семейства ¿к через скорости деформаций связующего материала ¿о, т. е. в конечном итоге опять же через ¿.
В настоящем исследовании рассматривается механическое поведение КМ-пластины как гибкой тонкостенной системы, поэтому напряжение азз(£,г) с приемлемой для инженерных приложений точностью можно линейно аппроксимировать по толщине [18]:
азз(1 ,г) =аз(1, г) = 33 ( , 33 ( , ) Хз + 33 ( , ) 2 33 ( , ), (14)
где а^^, х) = азз(£, х, ±К) —нормальные напряжения на нижней (—) и верхней (+) лицевых поверхностях конструкции, известные из соответствующих граничных условий.
Матричное соотношение (11) —система шести алгебраических уравнений. Из третьего равенства этой системы при учете условий соответствия, аналогичных (5), можем выразить скорость поперечной линейной деформации КМ-пластины:
¿зз = ¿з = 5^ (аз - £ (1 - ¿3г)Ъзг¿Л , (15)
г=1
где Ьц — компоненты матрицы В в (11); ¿зг (я = 1, 6) — символ Кронекера; аз определяется путем дифференцирования (14) по времени. В правой части (15) скорости деформаций ¿^ можно получить за счет дифференцирования по Ь соотношений (2), т. е. можно выразить их через двумерные функции ш,
-ш, щ, (1 = 1, 2, 0 4 т 4М).
Как и в [23], для понижения размерности уравнений движения гибкой КМ-пластины применим метод взвешенных невязок [27]. В качестве весовых функций используем полиномы Х3. Тогда с учетом (14) и последнего равенства (3) будем иметь [23]
2^ рш = ^ дЛм™ + £ М^9гш) + а33} - а<- + х30), =1 =1
= 2
рй[1) = £ д3 (м« - мЦдгш) - 1м1~1] + 1М&~ 1)дгш- (16) ,=1
-(О
'33
-Ь1 [а(3) - (-1)'а(3)] дгш + Х(
1 = 1, 2, 0 4 М + 1, х еС, г0,
где
к . . г Н
" ГН
р = роШо + £ Рк Шк, Х()(£, ж) = / Х^,
к=1 ■)-Н
/Н гН
, г)хг3йх3, и®(£, ж) = из(Ь, г)хг3йх3, -н }-н
[Н
(£, ж) = I / СТ33(£, г)х31йх33 _ ■)-Н
(Н гн
мУу, ж) = аф, г)х3^Х3, , ж) ^
, ^ '-Н ,н -1-Н (17)
■1-Н
_ | [(*3+ + а3-))_^- (-1)') + у+у (*3+ - а3-))(1 + (-1)')
г,./_Т73, в_1,2, 0 ^/^М + 1;
Ро, Рк — объемные плотности материалов связующего и арматуры й-того семейства; Хг — объемные нагрузки, действующие на КМ и определяемые по правилу смеси аналогично р. Из третьего равенства (17) следует, что М(о) = ^, М(1) = М^ — тангенциальные силы и изгибающие и крутящие
моменты в пластине; М^ = ^¿3 (я, ,7 _ 1, 2) —перерезывающие силы; остальные силовые факторы в (16) — моменты высших порядков [27].
Для однозначного интегрирования исследуемой задачи нужно задать граничные и начальные условия. Для понижения размерности этих условий вновь применим метод взвешенных невязок, используя те же весовые функции. Тогда на части кромки пластины (обозначим ее Гст), где заданы силовые граничные условия, получим следующие равенства в силовых факторах (17) [23]:
- М^^ад) + ^(М^ - М23)9гад) _ Р{()(1, ж),
£%(м<3о) + _р3о)(;,*), (18)
3=1 ^ 1=1 '
г _ 1, 2, 0 ^/^М + 1, х е Га, ¿о,
а на остальной части кромки (обозначим ее Гад), где заданы перемещения, получим граничные условия в кинематических переменных (см. (3) и (17)) [23]:
г _ 1, 2, 0 ^/^М + 1, х е Ги, ¿о,
, , ж), и()(£, ж)_и*/(£, ж),
- и
где
¡■Н Г-Н
Р,()(£, ж) = / р,-(£, г)хг3ЙХ3, , ж) = / , г)хг3ЙХ3
¿-к З-Н
(] _ 1, 3, г _ 1, 2), п1 _ео8 7, п2 _8Ш7;
р^- _ 1, 3) — приложенные к торцевой поверхности пластины внешние нагрузки; — известный на прогиб; и*г — заданные на торцевой поверх-
ности пластины тангенциальные перемещения; 7 — угол, определяющий направление функции.
правление внешней нормали к кромке Гст; (г _ 1, 2)—известные на
В момент времени ¿о нужно использовать начальные условия (см. (3), (17)) [23]:
■ш^о, х) = Щз(х), уи^о, х) = Уоз(х), у^^о, х) = ио01(х), у^^ о, х)= ьЦЦ(х), X е С, г = 1,2, 0 4 I 4 М + 1,
(21)
где
и.
0í
рп рп
(х) = ! ио^г)х13(1хз, уо}(х) = J Уо^г)х13(1хз, (22)
'-п J-n
х = 1, 2, 0 4 I 4 м + 1;
иоз, Уоу а = 1, 3) —заданные при ¿о перемещения и скорости точек КМ-пластины; и^, у^^ (г = 1, 2) — известные в области С функции.
Согласно выражениям (20), Р-о) и Р—известные на кромке Гст усилия,
приложенные в направлениях Xí и хз, а Р-1) (г = 1, 2)—известные механические моменты; остальные величины, стоящие в правых частях (18),— известные моменты высших порядков от внешних сил, заданных на торцевой поверхности.
Таким образом, чтобы однозначно проинтегрировать исследуемую задачу, в каждой точке области С при £ = ¿о нужно задать начальные условия (21) с учетом выражений (22), а на кромках пластины — силовые (18) или кинематические (19) граничные условия с учетом обозначений (20). Допускается задание и смешанных из (18), (19) граничных условий.
Проинтегрируем первые соотношения (3) по толщине пластины с весами х\ (0 4 I 4 М + 1), тогда при учете четвертого соотношения (17) получим матричные равенства
Се^ = и1 + -шд^, г = 1,2, (23)
где * = (24)
е- = 1У е{о) е{1) е{М-1) £(м)\т-
т = — (М + 2)-компонентный вектор-столбец, С = (с^) — матрица размерности (М + 2)х(М + 2), элементы которых определяются так:
= ^ - (-1)), н1+1
= 1+1^ + (-1)1), (25)
С1+1,т+2 =2к1+т+2(1 - (-1)1+т) 0 4 I 4 М + 1, 0 4 т 4 М.
(т + 1)(1 + т + 2) (т + 3)(1 + т + 4)
Из равенств (25) видно, что компоненты матрицы С и вектора-столбца 'ш рассчитываются один раз, поэтому уравнения (23) целесообразно преобразовать к виду
е- = С-1щ + Нод№, г = 1, 2, (26)
где
w = C-1w; (27)
C-1 — матрица, обратная C.
Если из каких-либо соображений в текущий момент времени известны значения функций w и uf^ (0 ^ I ^ М + 1), то из (26) с учетом (24), (27)
получаем значения функций ui и е^ (0 ^ т ^ М), определяющие осреднен-ные деформации композиции (2) и тангенциальные перемещения U (г = 1, 2) точек КМ-пластины (см. (3)).
2. Численный метод расчета. Для построения численного решения исследуемой задачи используем алгоритм шагов по времени [14,17,23,26,28], т. е. будем определять значения неизвестных функций в дискретные моменты времени tn+1 = tn + т (п = 0,1, 2,...), где т = const > 0 — заданный шаг по времени. Предполагаем, что в два предшествующих момента времени ts уже известны значения следующих функций:
W(ж) = w(is, ж), uf )(ж) = uf) (is, ж),
(33)(ж) = CT33)(iS, x), <73^)(Ж) = <3±)(is, Ж),
Gij(r) = aij(ts, Г), <iJfc)(r) = <iJfc)(is, r^ (28)
s _
X,-(r) = X(is, r), ¿ = 1, 2, j = 1, 3,
s = n - 1,п, 0 ^/^М + 1, 0 < k < К, ж eG, |жз| < k,
поэтому, согласно (17), в момент времени ¿п известны и все силовые факторы М®, которые входят в (18) и в правые части уравнений (16).
Чтобы получить явную численную схему, аппроксимируем встречающиеся производные по времени их центральными конечными разностями на трехточечных шаблонах [14,23]. Согласно этому, конечно-разностные аналоги уравнений движения (16) при учете обозначений (28) примут вид
п ,
^ ("Й1 -2 й + V) = Е % (М,? + Е М«?а й) + »<+» - з<-> + . <»>,
А ("+1<'> - 2^«'»+V«) = Еа (м§> - М>Й)- (29)
- МГ> +1М <^1>а«1Й-Л' [г <+> - (-1)' <"!,;>] а« й +.'>,
г = 1, 2, 0 ^/^М + 1, п = 1, 2, 3,...
Правые части в (29) при £ = Ьп известны, поэтому из этих уравнений, учи/ \ / \ "+1 "+1('> тывая соответствующие граничные условия (18) и (19), вычисляем а им «
п+1
в следующий момент времени ¿п+ь На основании (26) по уже известным а и "+1«'> (0 < I < М + 1) определяем "+1« и "+1(™> (г = 1,2, 0 < т < М),
п_1
а затем — деформации е ^ (см. (2)). При £ = ¿п_ 1 деформации е^ предполагаются уже известными (см. (2) при учете (28)), поэтому, используя (15)
и формулы численного дифференцирования по времени, можем вычислить
п
и скорости средних деформаций После чего по формулам (13) последо-
п
вательно определяем скорости деформаций в составляющих композиции ,
п
а из (4) —скорости напряжений в компонентах композиции ак (0 4 к 4 К). Дальнейшее решение рассматриваемой задачи строится совершенно так же, как это подробно описано в [14,23].
Согласно структуре левых частей уравнений (29), для начала проведения расчетов по разработанной явной численной схеме необходимо использовать
значения функций гг и известные из начальных условий (21) при учете (22), и требуется предварительно определить значения функций гг и и*}) (см. (29) при п = 1). Эти величины можно вычислить по формуле Тейлора, учитывая начальные условия (21) и уравнения движения (16) при £ = ¿о [23]:
1 о 0 г2 0
го(х) = го(х) + тг(х) + —ги(х) + 0(т3) « 0,
^) (х) = и?) (х) + ТЦ1) (х) + ^Ц1) (х) + 0(т3) « 0, X <ЕС, г = 1,2, 0 414М + 1.
Здесь приближенные равенства выполняются с точностью порядка г3, причем нулевые значения получаются, когда в начальный момент времени пластина покоится в естественном состоянии (см. (21) и (22) при иог = 0, Уо% = 0), а внешние нагрузки отсутствуют (см. (14), (16) и (17) при ^3±)(¿о, х) = 0,
Хг(1 о, г) = 0, 1 = 13).
В случае прямоугольной области С, заменяя в равенствах (18) и (29) производные дг(-) их конечно-разностными аналогами, в конечном итоге получим явную схему «крест» [14,17,23]. Если область С неканоническая, то дискретизация соотношений (18), (29) по пространственным переменным хг (г = 1, 2) может быть осуществлена путем применения вариационно-разностного подхода [17]. Необходимые условия устойчивости явной схемы «крест» вытекают из условий Куранта—Фридрихса—Леви [17] и для КМ-пластины характеризуются неравенствами (60) из [14].
3. Обсуждение результатов расчетов. В качестве примеров рассмотрим динамическое упругопластическое деформирование КМ-пластин толщиной 21ъ = 2 см, имеющих прямоугольную форму в плане ( С : |х 11 4 а, |х21 4 Ь, а = 4Ь). По всем кромкам пластины жестко закреплены (см. (19) и (20) при и*г = 0), а в начальный момент времени ¿о = 0 находятся в естественном состоянии покоя (см. (21) и (22) при иог = 0 и Уог = 0). Объемными нагрузками пренебрегаем (см. (16), (17) и (29) при Хг = 0, г = 1, 3). Конструкции нагружаются со стороны нижней лицевой поверхности давлением, соответствующим воздушной взрывной волне [28] (см. (14), (16), (17), (28) и (29)):
4+^) ^ 0,
0 4 ^ 4 ímax, -а(г - ¿тах)], t> ь тах)
^хо={ £тах
где
а = - ln(0.01)/(imin - ¿max) > 0, imin » U
(30)
tтах — момент времени, при котором нагрузка р(£) > 0 достигает наибольшего значения ртах; ¿т1п — момент времени, при котором р(£) становится пренебрежимо малым по сравнению с ртах (так, (30) соответствует предположению р(¿т1п) = 0.01 ртах). Согласно экспериментальным данным [28], в расчетах принято ¿тах = 0.1 мс и ¿т1п = 2 мс. Следовательно, при £ > 2 мс колебания исследуемых ниже пластин можно рассматривать как свободные.
Пластины изготовлены либо из эпоксидной смолы и армированы стекловолокнами марки 8-994 [29,30] (стеклопластики), либо из алюминиевого сплава АДМ и усилены стальной проволокой У8А [29] (металлокомпозиции). Упругопластическое поведение компонентов этих композиций на стадии их активного нагружения описывается идеализированной диаграммой с линейным упрочнением:
а
Е (fc)e,
sign(e)aifc) + Е^к\е - sign(e)eifc)),
|e| < 4fc) ^aifc)/E(fc), 0 < k < K,
^ (fc) > ,
где а, e— осевые напряжение и деформация при растяжении и сжатии; ст^,
Ее — условный предел текучести и модуль линеиного упрочнения материала k-той фазы композиции. Физико-механические характеристики компонентов указанных композиций приведены в таблице, где а = \JЕ/р — скорость звука в соответствующем материале.
Физико-механические характеристики компонентов композиций пластин [29,30] [Physico-mechanical characteristics of the components of composite plates [29,30]]
Materials p, kg/m3 V ae, MPa E, HPa Ee, HPa a, m/s
Epoxy 1210 0.33 20 2.8 1.114 1521.2
Fiberglass S-994 2520 0.25 4500 86.8 6.230 5868.9
ADM (Aluminum alloy) 2710 0.30 30 71.0 0.143 5118.5
U8A steel wire 7800 0.31 3968 210.0 6.973 5188.7
Структуры армирования считаются однородными и прямолинейными: = const, = const и = const, 1 ^ k ^ K (см. (1) и (7)). Исследуется динамическое поведение КМ-пластин со следующими структурами:
1) ортогональное плоское 2Б-армирование (рис. 1, а), когда два семейства волокон (K = 2) укладываются в направлениях Ж1 и с интенсивно-стями армирования wi = 0.266 и W2 = 0.324 соответственно;
2) ортогональное пространственное 3Б-армирование (рис. 1, b), когда три семейства волокон (K = 3) укладываются по направлениям Ж1, Ж2 и
с интенсивностями w1 = 0.235, w2 = 0.324 и w3 = 0.031 [9];
3) неортогональное пространственное 4Б-армирование (рис. 1, c), когда первое и второе семейства укладываются по направлениям Ж1 и Ж2, а третье и четвертое семейства — по траекториям, задаваемым углами (см. рис. 2 и (7)): 03 = ^/4, 04 = 3^/4, = = ^/2 (т. е. на рис. 1, c 0 = ^/4).
В последнем случае интенсивности армирования имеют значения: w = 0.126, w2 = 0.324 и w3 = w4 = 0.07. При всех структурах общий расход арматуры один и тот же.
Для проведения численных расчетов по направлениям х\ и Х2 введем равномерную сетку с шагом Дж1 = Дж2 = 26/100, а шаг по времени т зададим равным 0.25 мкс. Далее рассматриваются пластины с наименьшим размером в плане 2 b =1 ми 2 b = 20 см, для которых получаем Дж1 /т = 40 км/с и Дж1 /т = 8 км/с соответственно. Эти отношения значительно больше значений а, приведенных в таблице для всех компонентов композиций. Следовательно, необходимые условия устойчивости схемы типа «крест» выполняются для каждого компонента композиции (см. неравенства (60) в [14]), а значит, и для самих композиций, причем с запасом [14].
На рис. 3 изображены поперечные колебания центральных точек металло-композитных пластин с разными структурами армирования (wo(t) = w(t, 0, 0)), полученные при pmax = 16 МПа и b = 50 см. Кривые 1 рассчитаны по уточненной теории при М = 7 (см. (2), (3)) для конструкции с 3Б-армированием, а кривые 2 — либо по теории Редди (М = 0) для пластин с той же структурой 3Б-армирования (рис. 3, а), либо по уточненной теории (М = 7) для конструкции с плоской структурой 2Б-армирования (рис. 3, b). Сравнение кривых 1 и 2 на рис. 3, а показывает, что для 3Б-армированных пластин при t > 60 мс прогиб, определенный по теории Редди, начинает отличаться от прогиба, рассчитанного по уточненной теории; с увеличением расчетного интервала времени это различие возрастает. В еще большей степени проявляется соответствующее различие при сравнении деформированных состояний
а
Q Time in ms
b
Рис. 3. Осцилляции центральных точек тонких металлокомпозитных пластин с разными структурами: a) 3Б-армирование; расчет по уточненной теории и по теории Редди; b) 2D-и 3D-армирование; расчет по уточненной теории [Figure 3. Oscillations of the central points of thin metal-composite flexible plates with
different structures: a) spatial reinforcement; line 1—calculation by refined theory, line 2—calculation Reddy theory; b) line 1—spatial reinforcement, line 2—flat-cross 2D reinforcement; calculation by refined theory]
в компонентах композиции, рассчитанных по этим теориям. Так, на рассматриваемом интервале времени 0 4 t 4 100 мс максимальное (достигаемое на первой осцилляции) значение интенсивности деформаций в связующем, определенное по теории Редди, на 15.1 % меньше аналогичной величины, рассчитанной по уточненной теории.
Сопоставление кривых 1 и 2 на рис. 3, Ь свидетельствует о том, что замена 2Б-структуры армирования на 3Б-структуру с тем же расходом волокон в случае относительно тонких (Н/Ь = 1/50) металлокомпозитных пластин практически не приводит к уменьшению максимального значения прогиба, который достигается на первой осцилляции. Однако такая замена структуры армирования приводит к уменьшению наибольшего значения интенсивности деформаций в связующем на 7.4 %.
Согласно таблице, АДМ-У8А-композиция имеет слабую анизотропию (Е(к^/Е(0) « 3, 1 4 к 4 К). Традиционно считается, что расчеты упругого изгибного поведения КМ-пластин со слабой анизотропией вполне достаточно проводить по классической теории [19,21]. Поведение же кривых на рис. 3, а демонстрирует, что в случае динамического упругопластического изгибно-го деформирования 3Б-армированных пластин со слабой анизотропией даже неклассическая теория Редди при расчетных временах порядка 0.1 с и более не гарантирует получения надежных результатов определения прогибов. А значит, необходимо использовать более точные теории изгиба.
На рис. 4 приведены зависимости 'Шо(Ь) для стеклопластиковых пластин с разными структурами армирования, рассчитанные при ртах = 5 МПа и Ь = 50 см. Кривые 1 определены по уточной теории (М = 7) для конструкции с 4Б-армированием, а линии 2 — по теории Редди для пластины с той же структурой 4Б-армирования (рис. 4, а) или по уточненной теории для пластины с плоской 2Б-структурой (рис. 4, Ь).
Анизотропия стеклопластиковой композиции является сильно выраженной (Е(к)/Е(0) = 31, 1 4 к 4 К). Сравнение кривых 1 и 2 на рис. 4, а показывает, что и в случае 4Б-армирования относительно тонкой пластины зависимость -Ш0(£), рассчитанная по теории Редди (кривая 2), при £ > 110 мс существенно отличается от 'Шо(Ь), определенной по уточненной теории (кривая 1).
Сопоставление кривых на рис. 4, Ь свидетельствует о том, что замена плоской структуры 2Б-армирования на пространственную 4Б-структуру в случае относительно тонких стеклопластиковых пластин также практически не приводит к уменьшению максимального значения прогиба, который по-прежнему достигается на первой осцилляции. Кроме того, такая замена структуры армирования не позволяет заметно уменьшить и максимальное значение интенсивности деформаций в связующем материале.
Выше исследовались относительно тонкие КМ-пластины (Н/Ь = 1/50).
На рис. 5, 6 изображены зависимости -ш0(1) и ^т?^) = тах£*0)(£, г), х е С,
г
|жэ| 4 Н, где е*0) —интенсивность деформаций в связующем, рассчитанные при ртах = 6 МПа и Ь = 10 см (относительно толстые пластины из стеклопластика: Н/Ь = 1/10). Кривые с номерами 1 получены для конструкций с 4Б-армированием, а кривые с номерами 2 — с плоским 2Б-армированием. Сплошные кривые, номера которых помечены одним штрихом, определены
а
О Time in ms
b
Рис. 4. Осцилляции центральных точек тонких стеклопластиковых пластин с разными структурами: а) 4Б-армирование; расчет по уточненной теории и по теории Редди; b) 2D-и 4D-армирование; расчет по уточненной теории [Figure 4. Oscillations of the central points of thin fiberglass flexible plates with different structures: a) spatial 4D reinforcement; line 1—calculation by refined theory, line 2—calculation Reddy theory; b) line 1—spatial 4D reinforcement, line 2—flat-cross 2D reinforcement;
calculation by refined theory]
по уточненной теории ( M = 7), а штриховые кривые, номера которых помечены двумя штрихами, — по теории Редди (M = 0). Сравнение кривых 2 и 2' на рис. 5, а показывает, что в случае плоского 2Б-армирования относительно толстых стеклопластиковых пластин расчеты прогибов по уточненной теории и по теории Редди начинают заметно различаться уже в окрестности начального момента времени (при t > 2 мс). Поведение же кривых 1 ' и 1'' на рис. 5, а свидетельствует о том, что в случае пространственного 4Б-армиро-вания это различие не столь заметно, но с увеличение расчетного интервала времени также начинает проявляться и к моменту времени t = 100 мс это различие становится существенным (см. рис. 5, b).
В еще большей степени это различие проявляется при сравнении зависимостей e 2!)(i ). Действительно, из сравнения кривых 2 и 2' на рис. 6 видно, что в случае 2Б-армирования теория Редди завышает расчетное значение
maxeii = maxein^i ) (кривая 2'') по сравнению с аналогичной величиной,
то
определенной по уточненной теории (кривая 2'), на 18.4 %. Наоборот, из сравнения кривых 1' и 1'' на рис. 6 следует, что в случае 4Б-армирования теория
Редди занижает величину max ei на 10.8 % по сравнению с уточненным расчетом.
Сопоставление кривых 1 и 2 на рис. 5, а свидетельствует о том, что в случае относительно толстой стеклопластиковой пластины замена плос-
а
Рис. 5. Зависимости от времени поперечных колебаний центральных точек относительно толстых стеклопластиковых пластин с 2D- и 4D-структурами армирования, рассчитанные по разным теориям в окрестности начального момента времени (a) и t = 100 мс (b)
[Figure 5. Oscillations of the central points of thick fiberglass plates with different structures calculated by different theories: a) data for initial moment of time, b) data for moment of time t = 100 ms; line 1'—calculation by refined theory for spatial 4D reinforcement, line 1"—calculation by Reddy theory for spatial 4D reinforcement, line 2—calculation by refined theory for flat-cross 2D reinforcement, line 2''—calculation by Reddy theory
for flat-cross 2D reinforcement]
Time in ms с
Рис. 6. Зависимости от времени наибольших значений интенсивности деформаций в связующем
[Figure 6. Oscillations of the greatest strain intensity of the binding structure for thick
fiberglass plates; line 1 —calculation by refined theory for spatial 4D reinforcement, line 1 —calculation by Reddy theory for spatial 4D reinforcement, line 2 —calculation by refined theory for flat-cross 2D reinforcement, line 2 —calculation by Reddy theory
for flat-cross 2D reinforcement]
кой структуры 2Б-армирования на пространственную структуру 4Б-армиро-вания приводит к уменьшению максимального значения прогиба на 47.3%, а сравнение кривых 1 и 2 на рис. 6 показывает, что такая замена структуры
(о)
армирования приводит к уменьшению величины max m в два раза.
При упругопластическом деформировании КМ-пластины величину (о)
max m можно рассматривать как меру поврежденности материала связующей матрицы, поэтому чем больше значение max е 2°), тем больше повреждено связующее. Согласно поведению кривых на рис. 6, существенно меньшую по-врежденность связующего имеет относительно толстая пластина с пространственной структурой 4Б-армирования, причем она же имеет и наименьшую податливость в поперечном направлении (см. рис. 5, а).
Заключение. Сравнительный анализ изгибного динамического упруго-пластического деформирования пластин с пространственными и плоскими структурами армирования продемонстрировал, что в случаях относительно толстых КМ-пластин с сильно выраженной анизотропией композиции (стеклопластики) замена плоской структуры 2Б-армирования (рис. 1, а) на пространственную структуру 4Б-армирования (рис. 1, с) с сохранением общего расхода волокон позволяет уменьшить максимальные значения прогиба такой конструкции почти в 1.5 раза, а максимальные значения интенсивности деформаций в связующем — в два раза. В случаях относительно тонких пластин замена плоской 2Б-структуры армирования на пространственные структуры 3D- (рис. 1, b) и 4Б-армирования практически не приводит к уменьшению их податливости в поперечном направлении как при использовании композиций со слабой анизотропией (металлокомпозиты), так и при использовании композиций с сильно выраженной анизотропией (стеклопластики). Однако такая замена структуры армирования все-таки может приводить к уменьшению максимальных значений интенсивности деформаций в связующем материале на десятки процентов.
Проведенные расчеты показали, что при исследовании неупругого динамического поведения гибких КМ-пластин как с пространственными, так и плоскими структурами армирования традиционно используемая неклассическая теория Редди может приводить к результатам расчетов, которые существенно отличаются от полученных на основе уточненных теорий изгиба как для относительно тонких, так и для относительно толстых конструкций, причем как в случаях сильно-, так и в случаях слабовыраженной анизотропии композиции. С увеличением расчетных интервалов времени это различие значительно возрастает, поэтому для получения адекватных результатов расчетов неупругого динамического поведения гибких КМ-пластин целесообразно использовать более точные по сравнению с теорией Редди модели изгибного деформирования таких элементов конструкций.
Конкурирующие интересы. У меня нет конфликта интересов в авторстве и публикации этой статьи.
Авторская ответственность. Я несу полную ответственность за предоставление окончательной версии рукописи в печать. Окончательная версия рукописи мною одобрена.
Финансирование. Работа выполнена в рамках Программы фундаментальных научных исследований государственных академий наук на 2017-2020 годы (проект
№ 23.4.1 «Механика деформирования и разрушения материалов, сред при механических нагрузках, воздействии физических полей и химически активных сред»).
Библиографический список
1. Qatu M. S., Sullivan R. W., Wang W. Recent research advances on the dynamic analysis of composite shells: 2000-2009 // Composite Structures, 2010. vol. 93, no. 1. pp. 14-31. doi: 10. 1016/j.compstruct.2010.05.014.
2. Zhu F., Wang Z., Lu G., Nurick G. Some theoretical considerations on dynamic response of sandwich structures under impulsive loading// Int. J. Impact Eng., 2010. vol.37, no. 6. pp. 625-637. doi:10.1016/j.ijimpeng.2009.11.003.
3. Kazanci Z. Dynamic response of composite sandwich plates subjected to time-dependent pressure pulses// Int. J. Nonlin. Mech., 2011. vol.46, no. 5. pp. 807-817. doi: 10.1016/j. ijnonlinmec.2011.03.011.
4. Gill S. K., Gupta M., Satsangi P. Prediction of cutting forces in machining of unidirectional glass-fiber-reinforced plastic composites// Front. Mech. Eng., 2013. vol.8, no. 2. pp. 187200. doi: 10.1007/s11465-013-0262-x.
5. Gibson R. F. Principles of composite material mechanics. New York: Taylor & Francis Group, 2015. 815 pp.
6. Тарнопольский Ю. М., Жигун И. Г., Поляков В. А. Пространственно-армированные композиционные материалы: Справочник. М.: Машиностроение, 1987. 224 с.
7. Schuster J., Heider D., Sharp K., Glowania M. Measuring and modeling the thermal conductivities of three-dimensionally woven fabric composites// Mech. Compos. Mater., 2009. vol.45, no. 2. pp. 241-254. doi: 10.1007/s11029-009-9072-y.
8. Mohamed M. H., Bogdanovich A. E., Dickinson L. C., Singletary J. N., Lienhart R. R. A new generation of 3D woven fabric performs and composites// SAMPE J., 2001. vol.37, no. 3. pp. 3-17.
9. Тарнопольский Ю. М., Поляков В. А., Жигун И. Г. Композиционные материалы, армированные системой прямых взаимно ортогональных волокон. 1. Расчет упругих характеристик// Механика полимеров, 1973. №5. С. 853-860.
10. Жигун И. Г., Душин М. И., Поляков В. А., Якушин В. А. Композиционные материалы, армированные системой прямых взаимно ортогональных волокон. 2. Экспериментальное изучение// Механика полимеров, 1973. №6. С. 1011-1018.
11. Крегерс А. Ф., Тетерс Г. А. Структурная модель деформирования анизотропных, пространственно армированных композитов // Механика композитных материалов, 1982. №1. С. 14-22.
12. Янковский А. П. Определение термоупругих характеристик пространственно армированных волокнистых сред при общей анизотропии материалов компонент композиции. 1. Структурная модель // Механика композитных материалов, 2010. Т. 46, №5. С. 663678.
13. Соломонов Ю. С., Георгиевский В. П., Недбай А. Я., Андрюшин В. А. Прикладные задачи механики композитных цилиндрических оболочек. М.: Физматлит, 2014. 408 с.
14. Янковский А. П. Применение явного по времени метода центральных разностей для численного моделирования динамического поведения упругопластических гибких армированных пластин// Вычислительная механика сплошных сред, 2016. Т. 9, №3. С. 279-297. doi:10.7242/1999-6691/2016.9.3.24.
15. Reissner E. The effect of transverse shear deformations on the bending of elastic plates // J. Appl. Mech., 1945. vol.12, no. 2. pp. 69-77.
16. Vasiliev V. V., Morozov E. Advanced Mechanics of Composite Materials and Structural Elements. Amsterdam: Elsevier, 2013. xii+412 pp.
17. Абросимов Н. А., Баженов В. Г. Нелинейные задачи динамики композитных конструкций. Н. Новгород: Изд-во ННГУ, 2002. 400 с.
18. Богданович А. Е. Нелинейные задачи динамики цилиндрических композитных оболочек. Рига: Зинатне, 1987. 295 с.
19. Амбарцумян С. А. Теория анизотропных пластин. Прочность, устойчивость и колебания. М.: Наука, 1987. 360 с.
20. Reddy J. N. Mechanics of laminated composite plates. Theory and analysis. Boca Raton, FL: CRC Press, 2004. xxiii+831 pp.
21. Малмейстер А. К., Тамуж В. П., Тетерс Г. А. Сопротивление жестких полимерных материалов. Рига: Зинатне, 1972. 500 с.
22. Belkaid K., Tati A., Boumaraf R. A simple finite element with five degrees of freedom based on Reddy's third-order shear deformation theory // Mech. Compos. Mater., 2016. vol. 52, no. 2. pp. 367-384. doi: 10.1007/s11029-016-9578-z.
23. Янковский А. П. Построение уточненной модели упругопластического поведения гибких армированных пластин при динамическом нагружении // Механика композитных материалов и конструкций, 2017. Т. 23, №2. С. 283-304. doi: 10.25590/mkmk.ras.2017. 23.02.283_304.10.
24. Whitney J., Sun C. A higher order theory for extensional motion of laminated composites // J. Sound Vibration, 1973. vol.30, no. 1. pp. 85-97. doi: 10.1016/S0022-460X(73)80052-5.
25. Зубчанинов В. Г. Механика процессов пластических сред. М.: Физматлит, 2010. 352 с.
26. Иванов Г. В., Волчков Ю. М., Богульский И. О., Анисимов С. А., Кургузов В. Д. Численное 'решение динамических задач упругопластического деформирования твердых тел. Новосибирск: Сиб. унив. изд-во, 2002. 352 с.
27. Векуа И. Н. Некоторые общие методы построения различных вариантов теории оболочек. М.: Наука, 1982. 268 с.
28. Houlston R., DesRochers C. G. Nonlinear structural response of ship panels subjected to air blast loading// Computers & Structures, 1987. vol. 26, no. 1-2. pp. 1-15. doi: 10.1016/ 0045-7949(87)90232-X.
29. Композиционные материалы: Справочник / ред. Д. М. Карпинос. Киев: Наукова думка, 1985. 592 с.
30. Handbook of Composites/ ed. G. Lubin. New York: Springer, 1982. xi+786 pp. doi: 10. 1007/978-1-4615-7139-1.
Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki
[J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2019, vol. 23, no. 1, pp. 90-112
d https://doi.org/10.14498/vsgtu1644
ISSN: 2310-7081 (online), 1991-8615 (print)
MSC: 74K20
Modeling of elastoplastic behavior of flexible spatially reinforced plates under refined theory of bending
A. P. Yankovskii
Khristianovich Institute of Theoretical and Applied Mechanics,
Siberian Branch of the Russian Academy of Sciences,
4/1, Institutskaya st., Novosibirsk, 630090, Russian Federation.
Abstract
On the basis of the time-steps algorithm the structural model is constructed for elastic-plastic deformation of bended plates with spatial reinforcement structures. The inelastic behavior of the composition phase materials is described by equations of the theory of plastic flow with isotropic hardening. The possible weakened resistance of the reinforced plates to the transverse shear is taken into account on the basis of the refined theory, from which the relations of the Reddy theory are obtained in the first approximation. The geometric nonlinearity of the problem is considered in the Karman approximation. The solution of the formulated initial boundary value problems is based on an explicit numerical "cross" scheme. The dynamic inelastic deformation of spatially- and flat-cross-reinforced metal-composite and fiberglass flexible plates of different relative thickness is investigated in the case of the load caused by an air blast wave. It is demonstrated that for relatively thick fiberglass plates, the replacement of the flat-cross reinforcement structure by the spatial structure with the preservation of the total fiber consumption leads to a decrease in the structural flexibility in the transverse direction by almost 1.5 times, as well as to a decrease of the maximum of intensity of deformation in the binder by half. For relatively thin both fiberglass and metal-composite plates, the replacement of flat-cross 2D reinforcement structure with 3D and 4D spatial structures does not lead to a noticeable decrease in their deflections, but allows to reduce the intensity of deformations in the binder by 10 % or more. It is shown that the widely used non-classical Reddy theory does not allow obtaining reliable results of calculations of the elastic-plastic dynamic behavior of the bended plates, both with plane and spatial reinforcement structures, even with a small relative thickness of the structures and weak anisotropy of the composition.
Keywords: composite plate, spatial-cross reinforcement, flat-cross reinforcement, Reddy theory, refined theory of bending, elastic-plastic deformation, geometric nonlinearity, explosive load, "cross" scheme.
Research Article
Q ©® The content is published under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/) Please cite this article in press as:
Yankovskii A. P. Modeling of elastoplastic behavior of flexible spatially reinforced plates under refined theory of bending, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2019, vol. 23, no. 1, pp. 90-112. doi: 10.14498/vsgtu1644 (In Russian). Author's Details:
Andrei P. Yankovskii A https://orcid.org/0000-0002-2602-8357
Dr. Phys. & Math. Sci.; Leading Research Scientist; Lab. of Fast Processes Physics;
e-mail: lab4nemir@rambler .ru
Modeling of elastoplastic behavior of flexible spatially reinforced plates...
Received: 18th September, 2018 / Revised: 11th February, 2019 / Accepted: 4th March, 2019 / First online: 20th March, 2019
Competing interests. I hereby declare that I have no conflicts of interest in the authorship and publication of this article.
Author's Responsibilities. I take full responsibility for submitting the final manuscript in print. I approved the final version of the manuscript.
Funding. The research was (partly) carried out within the framework Program of the fundamental scientific research of the state academies of sciences for the years 2017-2020 (project no. 23.4.1 "Mechanics of deformation and destruction of materials, media, under mechanical loads, the influence of physical fields and chemically active media").
References
1. Qatu M. S., Sullivan R. W., Wang W. Recent research advances on the dynamic analysis of composite shells: 2000-2009, Composite Structures, 2010, vol.93, no. 1, pp. 14-31. doi: 10. 1016/j.compstruct.2010.05.014.
2. Zhu F., Wang Z., Lu G., Nurick G. Some theoretical considerations on dynamic response of sandwich structures under impulsive loading, Int. J. Impact Eng., 2010, vol.37, no. 6, pp. 625-637. doi:10.1016/j.ijimpeng.2009.11.003.
3. Kazanci Z. Dynamic response of composite sandwich plates subjected to time-dependent pressure pulses, Int. J. Nonlin. Mech., 2011, vol.46, no. 5, pp. 807-817. doi: 10.1016/j. ijnonlinmec.2011.03.011.
4. Gill S. K., Gupta M., Satsangi P. Prediction of cutting forces in machining of unidirectional glass-fiber-reinforced plastic composites, Front. Mech. Eng., 2013, vol. 8, no. 2, pp. 187-200. doi:10.1007/s11465-013-0262-x.
5. Gibson R. F. Principles of composite material mechanics. New York, Taylor & Francis Group, 2015, 815 pp.
6. Tarnopol'skii Yu.. M., Zhigun I. G., Polyakov V. A. Prostranstvenno-armirovannye kompoz-itsionnye materialy [Spatially Reinforced Composite Materials]. Moscow, Mashinostroenie, 1987, 224 pp. (In Russian)
7. Schuster J., Heider D., Sharp K., Glowania M. Measuring and modeling the thermal conductivities of three-dimensionally woven fabric composites, Mech. Compos. Mater., 2009, vol.45, no. 2, pp. 241-254. doi: 10.1007/s11029-009-9072-y.
8. Mohamed M. H., Bogdanovich A. E., Dickinson L. C., Singletary J. N., Lienhart R. R. A new generation of 3D woven fabric performs and composites, SAMPE J., 2001, vol. 37, no. 3, pp. 3-17.
9. Composite materials reinforced with a system of three straight, mutually orthogonal fibers. I. Calculation of the elastic characteristics, Polymer Mechanics, 1973, vol. 9, no. 5, pp. 754-759. doi: 10.1007/BF00856271.
10. Composites reinforced with a system of three straight mutually orthogonal fibers. 2. Experimental study, Polymer Mechanics, 1973, vol.9, no. 6, pp. 895-900. doi: 10.1007/ BF00856974.
11. Kregers A. F., Teters G. A. Structural model of deformation of anisotropic three-dimensionally reinforced composites, Mech. Compos. Mater., 1982, vol. 18, no. 1, pp. 10-17. doi: 10.1007/BF00605085.
12. Yankovskii A. P. Determination of the thermoelastic characteristics of spatially reinforced fibrous media in the case of general anisotropy of their components. 1. Structural model, Mech. Compos. Mater., 2010, vol. 46, no. 5, pp. 451-460. doi: 10.1007/s11029-010-9162-x.
13. Solomonov Yu. S., Georgievskii V. P., Nedbai A. Ya., Andryushin V. A. Prikladnye zadachi mekhaniki kompozitnykh tsilindricheskikh obolochek [Applied Problems of Mechanics of Composite Cylindrical Shells]. Moscow, Fizmatlit, 2014, 408 pp. (In Russian)
14. Yankovskii A. P. Using of explicit time-central difference method for numerical simulation of dynamic behavior of elasto-plastic flexible reinforced plates, Computational Continuum Mechanics, 2016, vol.9, no. 3, pp. 279-297 (In Russian). doi: 10.7242/1999-6691/2016.9. 3.24.
15. Reissner E. The effect of transverse shear deformations on the bending of elastic plates, J. Appl. Mech., 1945, vol. 12, no. 2, pp. 69-77.
16. Vasiliev V. V., Morozov E. Advanced Mechanics of Composite Materials and Structural Elements. Amsterdam, Elsevier, 2013, xii+412 pp.
17. Abrosimov N. A., Bazhenov V. G. Nelineinye zadachi dinamiki kompozitnykh konstruktsii [Nonlinear problems of the dynamics of composite structures]. Nizhni Novgorod, Nizhni Novgorod State Univ., 2002, 400 pp. (In Russian)
18. Bogdanovich A. E. Nelineinye zadachi dinamiki tsilindricheskikh kompozitnykh obolochek [Nonlinear Problems in the Dynamics of Cylindrical Composite Shells]. Riga, Zinatne, 1987, 295 pp. (In Russian)
19. Ambartsumyan S. A. Theory of Anisotropic Plates: ¡Strength, Stability, and Vibration. Stamford, Conn., Technomic Pub., 1970, viii+248 pp.
20. Reddy J. N. Mechanics of laminated composite plates. Theory and analysis. Boca Raton, FL, CRC Press, 2004, xxiii+831 pp.
21. Malmeister A. K., Tamuzh V. P., Teters G. A. Soprotivlenie zhestkikh polimernykh mate-rialov [Strength of Polymer and Composite Materials]. Riga, Zinatne, 1972, 500 pp. (In Russian)
22. Belkaid K., Tati A., Boumaraf R. A simple finite element with five degrees of freedom based on Reddy's third-order shear deformation theory, Mech. Compos. Mater., 2016, vol. 52, no. 2, pp. 367-384. doi: 10.1007/s11029-016-9578-z.
23. Yankovskii A. P. Construction of refined model of elastic-plastic behavior of flexible reinforced plates under dynamic loading, Mekhanika kompozitsionnykh materialov i konstruktsii, 2017, vol.23, no. 2, pp. 283-304 (In Russian). doi: 10.25590/mkmk.ras.2017.23.02.283_ 304.10.
24. Whitney J., Sun C. A higher order theory for extensional motion of laminated composites, J. Sound Vibration, 1973, vol.30, no. 1, pp. 85-97. doi: 10.1016/S0022-460X(73)80052-5.
25. Zubchaninov V. G. Mekhanika protsessov plasticheskikh sred [Mechanics Processes of Plastic Media]. Moscow, Fizmatlit, 2010, 352 pp. (In Russian)
26. Ivanov G. V., Volchkov Yu. M., Bogul'skii I. O., Anisimov S. A., Kurguzov V. D. Chislennoe reshenie dinamicheskikh zadach uprugoplasticheskogo deformirovaniia tverdykh tel [Solving Numerically Dynamic Problems of Elastoplastic Deformation of Solids]. Novosibirsk, Sib. Univ. Izd-vo, 2002, 352 pp. (In Russian)
27. Vekua I. N. Shell theory: general methods of construction, Monographs, Advanced Texts and Surveys in Pure and Applied Mathematics, vol. 25. New York, John Wiley & Sons. Inc., 1985, xvi+287 pp.
28. Houlston R., DesRochers C. G. Nonlinear structural response of ship panels subjected to air blast loading, Computers & Structures, 1987, vol.26, no. 1-2, pp. 1-15. doi: 10.1016/ 0045-7949(87)90232-X.
29. Kompozitsionnye materialy [Composite Materials. Handbook], ed. D. M. Karpinos. Kiev, Naukova Dumka, 1985, 592 pp. (In Russian)
30. Handbook of Composites, ed. G. Lubin. New York, Springer, 1982, xi+786 pp. doi: 10. 1007/978-1-4615-7139-1.