Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2021. Т. 25, № 2. С. 343-364 ISSN: 2310-7081 (online), 1991-8615 (print) d https : //doi
УДК 539.4
Уточненная модель вязкоупругопластнческого деформирования гибких пространственно-армированных цилиндрических оболочек
©А.П. Янковский
Институт теоретической и прикладной механики им. С. А. Христиановича СО РАН,
Россия, 630090, Новосибирск, ул. Институтская, 4/1.
Аннотация
Разработана модель вязкоупругопластнческого деформирования гибких круговых цилиндрических оболочек с пространственными структурами армирования. Мгновенное пластическое поведение материалов композиции определяется теорией течения с изотропным упрочнением. Вязкоупругое деформирование компонентов композиции описывается уравнениями модели тела Максвелла—Больцмана. Геометрическая нелинейность задачи учитывается в приближении Кармана. Используемые соотношения позволяют с разной степенью точности рассчитывать остаточные перемещения точек конструкции и остаточное деформированное состояние компонентов композиции. При этом моделируется возможное слабое сопротивление армированной оболочки поперечному сдвигу. В первом приближении используемые уравнения, начальные и граничные условия редуцируются в соотношения неклассической теории Амбарцумяна.
Численное решение сформулированной начально-краевой задачи строится по явной схеме «крест». Исследовано упругопластическое и вяз-коупругопластическое динамическое деформирование тонких стеклопла-стиковых оболочек под действием внутреннего давления взрывного типа. Рассматриваются две структуры армирования:
1) ортогональное армирование в продольном и окружном направлениях;
2) пространственное армирование в четырех направлениях.
Показано, что даже для относительно тонких композитных оболочек
теория Амбарцумяна неприемлема для получения адекватных результатов расчетов их вязкоупругопластнческого динамического деформирования. Продемонстрировано, что расчет по теории упругоплаетиче-ского деформирования армированных конструкций не позволяет даже приближенно оценить остаточные состояния композитных оболочек после их динамического нагружения. Показано, что даже для относительно тонкой и длинной цилиндрической оболочки замена традиционной
Научная статья
@® Контент публикуется на условиях лицензии Creative Commons Attribution 4.0 International (https : //creativecommons. org/licenses/by/4. 0/deed ,ru) Образец для цитирования
Янковский А. П. Уточненная модель вязкоупругонластического деформирования гибких пространственно-армированных цилиндрических оболочек // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2021. Т. 25, № 2. С. 343-364. https://doi.org/10.14498/vsgtul838.
Сведения об авторе
Андрей Петрович Янковский А https://отcid. org/0000- 0002- 2602- 8357
доктор физико-математических наук; ведущий научный сотрудник; лаб. физики быстро-
протекающих процессов; e-mail: lab4nemir@rambler .ru
© Самарский государственный технический университет
343
«шгосю»-перекрестной структуры армирования на пространственную структуру позволяет существенно уменьшить интенсивность остаточных деформаций связующего материала. В случаях относительно толстых и, особенно, коротких оболочек положительный эффект от такой замены структур армирования проявляется в значительно большей степени.
Ключевые слова: цилиндрическая оболочка, пространственное армирование, армирование по эквидистантным поверхностям, вязкоупруго-пластическое деформирование, нагрузки взрывного типа, уточненные теории изгиба, теория Амбарцумяна, геометрическая нелинейность, явная численная схема «крест».
Получение: 8 декабря 2020 г. / Исправление: 12 марта 2021 г. / Принятие: 11 мая 2021 г. / Публикация онлайн: 5 июня 2021 г.
Введение. Оболочки из композиционных материалов (КМ) находят широкое применение в инженерной практике [1—7], поэтому актуальна проблема моделирования их механического поведения в особенности при действии высокоинтенсивных нагрузок, характерных для современных КМ-изделий [7-9].
Механическое поведение анизотропных и слоистых оболочек в рамках классической теории изгиба моделировалось в [10,11], однако при этом не принимается во внимание плохое сопротивление таких конструкций поперечным сдвигам. Традиционно это негативное свойство тонкостенных КМ-изделий учитывают в рамках теорий Тимошенко—Рейсснера [5,7,12-14], Амбарцумяна [11,15,16] и Редди—Немировского [2,17]; реже применяются более точные теории [5,14,18—20]. При этом большинство публикаций посвящено рассмотрению линейно-упругого деформирования КМ-конструкций [2,5,10-14,17-19]. Однако при нагрузках высокой интенсивности компоненты композиции изделия могут деформироваться пластически [7,8,21,22]. Проблема моделирования неупругого поведения КМ-конструкций сейчас находится на стадии начального становления. Так, в работе [14] исследовалось упругопластическое деформирование слоистых конструкций с изотропными слоями. В [15,16,20] моделировалось вязкоупругопластическое поведение армированных гибких пластин, пологих и цилиндрических оболочек, причем в [15,20] рассматривались только традиционные «плоско»-перекрестные структуры армирования. Известно, что тонкостенные элементы конструкций с таким армированием плохо сопротивляются поперечному сдвигу и отрыву [9]. Для ликвидации этого серьезного недостатка КМ-изделий используют пространственные структуры армирования [9,20,23,24]. В работе [16] была разработана структурная модель вязкоупругопластического деформирования пространственно-ар-мированных сред, примененная для исследования механического поведения гибких пологих оболочек.
В [17] было показано, что в случае линейно-упругого деформирования компонентов композиции армированных оболочек уточнения расчетов, выполненных по теории Редди—Немировского, не требуется. (Напомним, что в теориях Редци—Немировского и Амбарцумяна [11] используются одинаковые кинематические соотношения, а различие заключается в уравнениях силового и моментного балансов.) Однако в [20] продемонстрировано, что
при исследовании динамического вязкоупругопластпческого поведения гибких цилиндрических оболочек, традиционно армированных в окружном и продольном направлениях, целесообразно использовать кинематические соотношения более высоких порядков, чем в теориях Редди—Немировского и Амбарцумяна. При этом открытым остался вопрос: нужно ли использовать теории изгиба КМ-оболочек высоких порядков точности при пространственных структурах армирования таких конструкций?
Решения нелинейных задач динамики тонкостенных конструкций, как правило, строят с использованием явных численных схем [14-16,20], например схемы типа «крест», или неявных методов Ньюмарка [8,25].
В связи с вышеизложенным данная работа посвящена моделированию вязкоупругопластпческого деформирования пространственно-армированных цилиндрических оболочек в рамках уточненной теории изгиба. Численные решения соответствующих начально-краевых задач предполагается при этом строить по явной схеме «крест».
1. Моделирование вязкоупругопластического деформирования пространственно-армированной цилиндрической оболочки. Рассматривается тонкая замкнутая круговая цилиндрическая оболочка радиуса Я, длиной Ь и толщиной 2/?, <С мп(1, Я) (рис. 1). С оболочкой связана цилиндрическая система координат Хг так, что х\ — осевая координата (0 ^ Х\ ^ Ь), Х2—угловая координата в окружном направлении (0 ^ х-2 < 27т), х$ — радиальная координата . Я /> Я ■ /г). Конструкция усилена «плоско» -перекрестно (традиционное армирование по эквидистантным поверхностям) или пространственно-перекрестно N семействами волокон с плотностями армирования ¿1,- (1 к ^ И). Структура армирования по толщине оболочки квазиоднородна. На рис. 2 изображены представительные элементы таких композиций: «плоская» ортогональная 2Б-структура армирования (рис. 2, а) и неортогональная пространственная 4Б-структура (рис. 2, Ь) [9].
(к)
С волокнами к-то семейства свяжем ортогональную систему координат х\ так, как показано на рис. 3. Направление армирования при этом задается двумя углами сферической системы координат вд; и ¡Дц, а направляющие косинусы между осями х^ и ;Г; = 1,3, 1 "С к ^ ./V) вычисляются по формулам (27) из [16].
Рис, 1. Цилиндрические КМ-оболочки с жестким закреплением левой кромки (а) и с жестким закреплением обеих кромок (b) [Figure 1. Cylindrical CM-shells with rigid fastening of the left edge (a) and with rigid fastening of both edges (b)]
a
b
ухз
^yowvyyyvv.
рооооооооло
роооооооооо
¿/щшшшшшшшяшш Х1 Щшшшпппп
_Х2
Рис. 2. Репрезентативная ячейка композиции с ортогональным 2В-армированием (a.) и с пространственным 40-армированием (b) [Figure 2. A representative cell of a composition with orthogonal 2D-reinforcement (a) and with spatial 4D-reinforcement (b)]
Рис. 3. Локальная система координат, связанная с волокнами k-го семейства
[Figure 3. Local coordinate system associated with A^th family of fibers]
Рассматривается частный, но практически важный случай нагружения оболочки, когда на ее лицевых поверхностях (жз Я =Ь /г) можно пренебречь распределенными внешними касательными силами (рис. 1, а). При этом предполагается: пространственные структуры армирования таковы, что если арматура некоторого к-то семейства имеет наклонную траекторию (0 < % < < 7г/2), то обязательно найдется 1-е семейство наклонных волокон, изготовленных из того же материала и имеющих параметры армирования 0/ = п—О^, Ш = <Рк, = Шь, (см. рис. 3). Структуры армирования с та-
ким свойствами, как правило, и встречаются на практике [9,24]. В частности, к ним относится структура, изображенная на рис. 2, Ь, а также ортогональная ЗБ-структура армирования [24]. (При нарушении этих двух условий не удается разработать явную численную схему типа «крест».)
Для учета возможного плохого сопротивления армированной оболочки поперечным сдвигам [2,5,7,11,13,14,17,20] используем более точные кинематические соотношения, чем в теориях Редци—Немировского [2,17] и Амбарцу-мяна [11]. При этом, согласно [20], осредненные деформации композиции ец и перемещения точек гибкой цилиндрической оболочки V% в рамках уточненной теории аппроксимируются так (геометрическая нелинейность задачи учитывается в приближении Кармана):
4|i(i, 1-) = dim - zdfw + f[m\z)dх4зг) + ^(diw)2,
m=О
£22 (t, г) =
1
R + z 1
R + z R
м
d2u2 + diw +W+J2 ft\z)d2еЩ +1(
Ж
m=О
(m),
1 / 0¿tt! \2
2\R + Z,
2e12(t, r) = R ^ d2m - zdid2w + ^ f[m\z)d2e
Ám) -13
m=0
+^dlU2 + ft&t» + £ ft\z)d+
м
(m).
+
(1)
R
Г) =
/г2 - z2 h2
м
m=0
E(í)m4m)(tx), ¿ = 1,2;
d\wd2w,
m=0
M
R + Z
U2(t, r) = —p«2(í,x) + + £ /ГЧ*)4з
|/l(t, r) = Ul(í, x) - zftu, + £ /Ím)(«^(í, x),
m=0
(m) f л ^(m).
M
(2)
Д
m=O
где
f%(í,r) = w(t,x); x = € O, O = {x : 0 ^ ®i ^ L, 0 ^ x2 < 2тг},
r = Щ), x-3 = R + Z, 4 h, t to,
Úm\¿) - p - , < = 1, 2,
1 v ; hm(l + m) 2 hm J Фm(z) = L l--Rm4R + z)-,
(3)
1=0
m — l
r — вектор места; x — вектор точки срединной поверхности (z = 0); z — введенная для удобства новая радиальная переменная; w — прогиб; щ, ^—перемещения точек срединной поверхности в осевом и окружном направлениях; t-o — начальный момент времени £; Q — область, которую занимает срединная поверхность оболочки в координатах Х\, х2\ <9¿ —оператор частного дифференцирования по Xi\ М — целое число, задающее количество слагаемых, которые удерживаются в степенных разложениях деформаций поперечного сдвига e¿3- При М = 0 из (1)—(3) получаются кинематические соотношения теорий Редди—Немировского [2,17] и Амбарцумяна [11]. В равенствах (1) и (2) неизвестны функции w, Ui и е¿з^ (i = 1, 2, 0 ^ m М), зависящие от времени t и двух пространственных координат х2.
Материал k-то компонента композиции предполагается изотропным, а его вязкоупругопластическое деформирование описывается определяющими уравнениями [15,16,20]:
rí? = - BWctW + SÍÍDWCTW -
а,
v
гз
-11 ^ Uij иц
-¿(4fc>eM. m=IA o (4)
где
<"' хт__^_
1+9(Ю' 9 0(4- л (1 + ,/(Ч)(1-2|/«>>)'
-(А)
2 1 + г/Ш' 3 1 - 21/(*0 '
0, Г, < тЯ{Х«) иди Тк = т?\ < (т^{Х«))2/Л
1, п = 4*\х{к)1 > (г^оЛ )7Л
По
(к) (к) , (к) (к) 1 (А=)
=4 -Мо > =д4 ;
0 — , е^- —компоненты тензоров напряжений и малых деформаций; р\- —
компоненты тензора пластических деформаций (р^ = 0); — мгно-
венные модуль Юнга и коэффициент Пуассона; —коэффициенты
линейной вязкости при чистом сдвиге и объемном расширении-сжатии; ть- — функция упрочнения при чистом сдвиге, которая зависит от параметра Од~ квиста @ —касательный модуль на мгновенной диаграмме чистого
сдвига т'^ ~ 7^ — пластическая составляющая полной угловой деформации при сдвиге; т^—касательное напряжение при сдвиге; —символ Кронекера; индекс к — номер компонента композиции (к = 0 — связующее,
1 ^ к ^ ./V —арматура к-го семейства); точка — дифференцирование по времени. Соотношения (4) при учете (5) получены в [15] в предположении, что
(к)
деформации е- можно разложить на сумму вязкоупругих и пластических
(к)
составляющих. При этом пластическое поведение материала описывается уравнениями теории течения с изотропным упрочнением, а вязкоупругое деформирование — уравнениями модели тела Максвелла—Больцмана.
Как и в работах [15, 16, 20], решение исследуемой задачи будем строить численно, используя пошаговый алгоритм [8,25]. Согласно этому, значения искомых функций будем вычислять в дискретные моменты времени Iп~\ = ¿/г + А, п = 0,1,2,..., где А —шаг по времени. В [16] была разработана структурная модель вязкоупругопластического деформирования КМ с пространственным армированием, механическое поведение материалов композиции которого описывается соотношениями (4) при учете (5). При этом в рассматриваемый момент времени Ьп определяющее уравнение для композиции записывается в матричной форме
п &
п п п В е + Р,
п = 0,1,2,...
(6)
где
a = (<7i a-2 a-s 04 a5 сг6)Т = (сц <т22 033 0"23 C31 ci2)T,
£ = (с | £2 £3 ^б)Т = (til £22 ^33 2tT23 2„31 2ti2)T, (7)
P = (p 1 P2 Pi P4 P5 Pe)T;
n n n n
&{r) = &(tn, r), e(r) = ¿{tn, r), P(r) = p(tn, r), B(r) = B(tn, r); (8)
iTijt €ij — компоненты тензоров усредненных напряжений и деформаций в КМ; В — матрица размером 6x6; индекс т — операция транспонирования.
Согласно [16], соотношение (6) линеаризовано по методу, аналогичному методу переменных параметров упругости [26], причем матрица В и вектор р с учетом (7) вычисляются по матричным формулам (24) из [16]. На данной итерации указанного метода при t = tn элементы матрицы В и вектора-столбца р уже известны и определяются механическими характеристиками материалов композиции, напряженным состоянием в них (которое рассчитано на предыдущей итерации) и структурой армирования: углами 9k, <fik (см-рис. 3) и плотностями Wfc (1 к N) армирования.
Далее моделируется механическое поведение КМ-оболочки как гибкой тонкостенной системы, поэтому нормальное напряжение 033 (t, г) с приемлемой для инженерных приложений точностью можно линейно аппроксимировать по координате z [13]:
<733(f, Г) ЕЕ (73 i, г) = 33 ' 33 JZ +
21)
+ 33 v ; 33 v , ^ х € |„| ^ h? t ^ ^ ^
где х) = (7зз(^х, ±/г) — напряжения на внешней (+) и внутренней (—)
лицевых поверхностях, которые известны из силовых граничных условий.
Матричное равенство (6) является системой шести алгебраических уравнений. В силу условий соответствия (7) из третьего равенства этой системы выразим скорость линейной поперечной деформации конструкции
е
11 п п 1 /п х - п п П \
¿33 = ¿з = (баз) [¿а ~ С1 " М Ъм-Рз), (10)
г=1
п __п
где Ьзг (?' = 1,6)—элементы матрицы В в уравнении (6); производная <тз определяется путем дифференцирования по времени выражения (9). Скорости деформаций ¿г в правой части (10) получаются дифференцированием
по I равенств (1), т. е. зависят от функций го, го, щ, егд (0 ^ т ^ М) и их частных производных по щ (I = 1,2).
Для окончательной формулировки начально-краевой задачи вязкоупруго-иластического изгибного деформирования тонких армированных цилиндрических оболочек в рамках уточненной теории к равенствам (1), (2), (6), (10) следует добавить двумерные уравнения движения и соответствующие начальные и граничные условия (см. [20], соотношения (7), (17), (19) и (21)).
2. Численный метод расчета. В разделе 1 уже отмечалось, что для численного интегрирования сформулированной начально-краевой задачи будем использовать метод шагов по времени [8,14—16,20,25,26]. Согласно этому, предполагаем, что в дискретные моменты времени ^ уже известны следующие функции:
т/ \ ^ т(1)/ \ (О/Л т(±)с (±)/.1 \
•ш(х) ее щ (х) = Щ (£т,х), СГззЧх) = (т^з'(¿т.х),
= 4з ау(г) ее (7у(£т,г), г = 1,2, 3 = 1,3, (И)
т = II - |,п, 0 < / «С Л / • 1. х е О, < /г,
где
Л
(£, X) = / г)х1йх% г = 1,2, 0 < I < М + 1, (12)
■1-1),
— введенные новые кинематические неизвестные переменные [20].
За исключением необходимых преобразований уравнений (4) (см. [16], формулы (19) и (20)), обеспечивающих устойчивость разрабатываемой численной схемы [15], производные по времени далее заменяем центральными конечными разностями на трехточечном шаблоне {¿п-ь • Это позво-
ляет построить явную схему численного интегрирования рассматриваемой задачи. После аппроксимации в уравнениях движения вторых производных по времени от кинематических переменных го и (см. (12)) их конечно-разностными аналогами, учитывая обозначения, аналогичные (8) и (11), получим [20]
2Нр /п+1 Л п п- к „ / п (о) п(0)~ п —( го -2го + го ) = (м^ + М„91 го + и>) Ф
ТЬ тъ ^ ^
+ а2(мГз} + (1з)
- г«« ) = эх(м® -л^&й) -
- I +1М^дг & -Ь1 - (-1)г 5Н] 01 й,
- 2^ + = д1 (Й2 - + п ^ ть
+ 92 (м§ - ^23 го) - I + 1 Мзз1]с>2 & ~
- И1 [(Д + /г)"1 ?<+> - (-1)г(й - /г)"1 +
0 % I < М • 1. хеО, п = 1,2,3,...,
где (см. выражение (9))
N А
р = роШо + Е рксок, х) = / т)г1йя,
к= 1
¡■Н
= 1 / (ТззС^г)
г1 =
¥
(4з} + 4з0(1 - (-1)') + - 0(1 + (-1)0]. (м
Д")
= I£ = ¿(4з} + - ф|(-л» +
N
ООО = 1 - ^ и^, г, э = ТТЗ, 0 < I < М + 1;
а-=1
Ро, Ра- — объемная плотность материала связующего и волокон к-то семейства; функции Ф/(л) определены в (3). В уравнениях (13) массовые нагрузки не учитываются.
На основании формул (14) с учетом (11) при t = можем определить все силовые
факторы йЦ\ М® и внешние силы с!^, которые входят
в правые части равенств (13). Учитывая при этом необходимые граничные условия (см. (17) и (19) в [20]), из (13) определяем по явной схеме значения
функций гю , // : .0 С / х? М 1. и следующий момент времени После этого по формуле (23) из [20], которая вытекает из (2) и (12) при учете (3),
вычисляем функции и| и П£ г = 1,2, 0 < т ^ М. Затем на основании
(1) получаем осредненные деформации композиции В предшествующий
момент времени деформации уже известны (см. (1) при учете (11)), поэтому, используя численное дифференцирование по времени и соотношение
п
(10), определяем и скорости деформаций ¿^ в текущий момент времени 1П.
Дальнейшее численное решение исследуемой задачи строится совершенно так
же, как и в работах [16,20].
Структура левых частей равенств (13) свидетельствует о том, что для
начала проведения расчетов по разработанной явной схеме нужно знать зна-0 0 п)
чения функций ио и Щ , известные из начальных условий (см. (21) в [20]),
а также ни и и^ (см. (13) при п = 1). Значения последних функций вычисляются по формуле Тейлора при учете начальных условий и уравнений движения, рассматриваемых при t = I о (см. последние формулы раздела 2 в [20]).
Для замкнутой круговой цилиндрической оболочки область Г2 (см. (2)) является прямоугольной. Следовательно, аппроксимируя в соотношениях (13) производные с^(-) их конечно-разностными аналогами от уже известных (см. (14) при учете (11)) в текущий момент времени Ьп сеточных функций, окончательно получим явную численную схему «крест» [14-16,20]. Необходимые условия устойчивости этой схемы для однородных тонкостенных конструкций приведены в [14]. Если эти условия выполняются для каждого компонента композиции, то они с запасом выполняются и для КМ-оболочки.
3. Обсуждение результатов расчетов. Исследуем неупругое динамическое изгибное поведение относительно тонкой замкнутой цилиндрической КМ-оболочки толщиной 2/г = 2 см, радиуса Я = 1 м и длиной Ь = 4 м (2/г/гшп(Ь, Я) = 1/50), которая жестко закреплена на кромках (го = 0,
и^Р = 0 при х\ — О, Ь и £ > ¿о [20]) и до начального момента времени I = ¿о = 0
покоится в естественном состоянии (го = 0, п^ = 0, г = 1,2, 0 ^ I ^ М + 1, х £ О, £ < [20]). Конструкция нагружается изнутри равномерным давлением р(1) по закону (см. (9), (13), (14)) [25]:
(+) _
33
= 0. -а.
(-) 33
=p(t) =
\ Pmaxt/ tri
о ^ i < L,
\ртахбхр[ Ql(t imax)]i t ^maxi
a = - ln(O-Ol)/(imin - imax) > 0, imin > tmax, p(imln) = 0.01j9max,
(15)
где смысл параметров нагрузки ртах, ¿тах, ¿тш и а вполне очевиден и подробно описан в [16,20]. Используя экспериментальные данные [25], в расчетах примем ртах = 8 МПа, ¿тах = 0.1 мс и £т;п = 2 мс. Эти данные условно характеризуют давление, вызванное воздушной взрывной волной.
Оболочка изготовлена из эпоксидной смолы и армирована стекловолокнами. Мгновенное упругопластическое поведение к-го материала композиции при активном нагружении описывается диаграммой растяжения-сжатия с линейным упрочнением
а =
sign(fc)(rifc) + Е^к)(е - sign(e)eifc))
и ^ 4к} ее ¿к)/е(к\ -{к) 0 € к < N,
t > £в
где а, е — осевое напряжение и линейная деформация; Е$ —модуль упроч-
(к)
нения; —предел текучести. Физико-механические характеристики материалов композиции приведены в таблице, причем объемная вязкость не учитывается: —> оо (см. (5)).
Physico-mechanical characteristics of the composite components [21,22]
Components р, kg/m3 V ав, МРа Е, НРа Ев, НРа г], MPa-s
Ероху 1210 0.33 20 2.8 1.114 250
Fiberglass 2520 0.26 4500 86.8 6.230 1000
Рассматриваются две структуры однородного армирования:
1) продольно-окружное 2Б-армирование (рис. 2, а) с интенсивностями армирования Ш\ = 0.1 и Ш2 = 0.3 в осевом и окружном направлениях;
2) пространственное 4Б-армирование, при котором два первых семейства волокон также уложены в продольном и окружном направлениях с плотностями армирования ш\ = 0.05 и ш2 = 0.3, а третье и четвертое семейства стекловолокон уложены наклонно с плотностями Шз = 004 = 0.025.
В случае 4Б-структуры направления армирования определяются следующими значениями углов сферической системы координат (см. рис. 3): 9\ = = д2 = 7г/2, вз = 7г/4, 04 = Зтг/4, = (рз = = 0 и (р2 = 7г/2 (на рис. 2, Ь
угол в = 7г/4 и оси ж1, Х'2 следует поменять местами). В обеих структурах армирования расход волокон одинаков.
При указанных условиях армирования, закрепления и нагружения цилиндрическая КМ-оболочка деформируется осесимметрично (решение начально-краевой задачи не зависит от окружной переменной Жд), а также симметрично относительно центрального сечения х.\ = Ь/2, например ж) = —ж), где ж = Ж1 — Ь/2 и |ж| ^ Ь/2.
На рис. 4 изображены зависимости от времени наибольшего абсолютного значения прогиба КМ-конструкции
и'ти) = тах |го(£,Ж1)|.
Номера кривых на рис. 4 соответствуют номеру структуры армирования. Кривые 1 и 2 получены при указанных выше входных данных. Кривая 1" рис. 4, а рассчитана при тех же условиях, что и кривая 1, но в рамках теории упругопластического деформирования материалов композиции (при
Г]
(к)
оо, 0 ^ к ^ И] см. (5)). Поведение кривой 1" показывает, что согласно упругопластическому расчету, КМ-конструкция колеблется неограниченно долго. Кривые же 1 и 2 свидетельствуют о том, что к моменту времени I = 400 мс осцилляции наибольшего по модулю значения прогиба практически прекращаются, а величина максимального остаточного прогиба имеет значение, равное примерно 1.3 мм. Сопоставление кривых 1 и 1" на рис. 4, а показывает, что упругопластические расчеты (кривая 1") не позволяют определить величину наибольшего остаточного прогиба цилиндрической КМ-обо-лочки даже весьма приближенно.
250 Time, t, ms
500
30
20
J 10
0 50 100 150
1"', 2"'
100 150 200 250 300 350 400 450 500 Time, i, ms
b
Рис. 4. Осцилляции максимального абсолютного значения прогиба: (а) — уточненные упругопластический и вязкоупругопластический расчеты КМ-оболочки с 2В-армировани-ем; (Ь) — уточненный вязкоупругопластический расчет оболочки с 4В-армированием
[Figure 4. Oscillations of the maximum absolute value of the deflection: (a) — refined ela.stopla.stic. and viscoela.stopla.stic calculations for the CM-sliell with 2D-reinforcement; (b) — refined viscoela.stopla.stic calculation for the shell with 4D-reinforcement]
На рис. 5 приведены зависимости ио(х) (х = х\ — Ь/2), рассчитанные при I = 500 мс. (Согласно поведению кривых 1 и 2 на рис. 4, осцилляции прогиба КМ-конструкции при этом практически уже прекратились.) Кривые 1 и 2 на рис. 5 определены при тех же условиях, что и кривые 1 и 2 на рис. 4. Как указывалось выше, решения исследуемой задачи обладают симметрией относительно центральных сечений конструкций х = 0 {х\ = Ь/2), поэтому на рис. 5 изображены только правые части зависимостей и){%) (|ж| ^ Ь/2|. Сравнение максимальных значений ординат точек на кривых 1 и 2 рис. 5, которые достигаются в зоне краевого эффекта (х ~ 1.84 м), показывает, что при t = 500 мс значения «>т на кривой 2 примерно на 2 % больше, чем на кривой 1. Сопоставление же этих кривых при 0 ^ ж ^ 1.3 м свидетельствует о том, что при использовании 4Б-структуры армирования (кривая 2) остаточный прогиб в центральной части даже относительно тонкой КМ-конструкции на 10-20% меньше аналогичной величины в цилиндрической оболочке с традиционным 2В-армированием.
Все кривые на рис. 4 и 5 получены по уточненной теории изгиба оболочек при М = 7 (см. (1) и (2)). Расчеты, выполненные по теории Амбарцумяна (М = 0), приводят к зависимостям *шт(£) и при I = 500 мс, которые визуально практически не отличаются от кривых, приведенных на рис. 4 и 5, поэтому не изображены.
На рис. 6 представлены зависимости от времени наибольших значений
(А')
интенсивности деформаций материалов композиции :
= тахг£°,
г
где |л| ^ /г, 0 -п /. для цилиндрической оболочки с 4В-структурой армирования в окрестности начального момента времени (рис. 6, а) и в окрестности I = 500 мс (рис. 6, Ь). Номера кривых соответствуют номеру к-то компонента композиции: 0, 0' — связующее (к = 0); 2, 2' — арматура второго (окружного) семейства (к = 2). Сплошные кривые 0 и 2 рассчитаны по
х = xi — L/2, ш
Рис. 5. Эпюры остаточных прогибов цилиндрических КМ-обо-
лочек, рассчитанные по уточненной теории [Figure 5. Diagrams of residual deflections of cylindrical CM-shells calculated by the refined theory]
уточненной теории (М = 7), а штриховые кривые 0' и 2'— по теории Ам-барцумяна (М = 0). На рис. 6, а кривые 2 и 2' визуально практически не различаются. Кривые 0, 0' и 2, 2' на рис. 6 получены при прежних входных данных.
Сравнение ординат глобальных максимумов на кривых 0 и О1 (при t = = 0.5 мс) на рис. 6, а показывает, что расчет по традиционной теории Ам-
барцумяна завышает наибольшее значение ет' в связующем оболочки с пространственной структурой армирования на 29 % по сравнению с расчетом, выполненным по уточненной теории. Сопоставление же кривых 0, 2 и 0', 2' на рис. 6, Ь демонстрирует, что теория Амбарцумяна примерно на 35 % завышает величину интенсивности остаточных деформаций в компонентах композиции по сравнению с расчетом, выполненным по уточненной теории.
На рис. 7 изображены зависимости рассчитанные по уточненной
теории (М = 7) в окрестности начального момента времени (рис. 7, а) и в ок-
2.0
1.6
ё? 1.2
0.4
0 2 4 8 10 12 14 16 18 20
Time, t, ms
(f /
0
/
/
\ 2
480 482 484 486 488 490 492 494 496 498 500 Time, t, ms
Рис. 6. Осцилляции наибольших значений интенсивности деформаций компонентов композиции оболочки с 4В-армированвем, рассчитанные по разным теориям: (а) — в окрестности начального момента времени; (Ь) —
в окрестности t = 500 мс [Figure 6. Oscillations of the highest values of the deformation intensities of components of shell compositions with 4D-reinforcement calculated by the different theories: (a) — close to the initial moment of time; (b) — close to
the time t = 500 ms]
рестности I = 500 мс (рис. 7, Ъ) для КМ-оболочки с 2Б-структурой армирования (сплошные кривые 0, 2) и с пространственной 4В-етруктурой (штриховые кривые 0', 2'). Как и на рис. 6, номера кривых соответствуют номеру к-то компонента композиции. Сравнение кривых 0 и О1 на рис. 7, а показывает, что для относительно тонкой цилиндрической оболочки замена 2В-струк-туры армирования на 4В-структуру приводит к увеличению максимального значения е^ примерно на 9% (см. кривые 0 и 0' при I ~ 0.5 мс). Однако сопоставление этих же кривых на рис. 7, Ь свидетельствует о том, что к моменту времени I = 500 мс конструкция с 4В-структурой армирования (кривая 0') практически перестала осциллировать, а оболочка с 2В-структу-рой (кривая 0) все еще колеблется. Сравнение ординаты точки кривой 0' при £ = 500 мс со средним значением ординат точек кривой 0 в окрестности этого момента времени демонстрирует, что замена 2С-структуры на 4Б-структу-ру армирования позволяет уменьшить величину интенсивности остаточных
8? 1.2 -
0.4 -
0
т
\
9 10
480 482 484 486 488 490 492 494 496 498 500 Time, t, ms
Рис. 7. Осцилляции наибольших значений интенсивности деформаций компонентов композиций оболочек с разными структурами армирования, рассчитанные по уточненной теории: (а) — в окрестности начального момента времени; (Ь) —в окрестности t. = 500 мс
[Figure 7. Oscillations of the highest values of the deformation intensities of components of shell compositions with different structures of reinforcement calculated according to the refined theory: (a) — close to the initial moment of time; (b) — close to the time t = 500 ms]
деформаций связующего примерно на 30%. Из поведения кривых 2 и 2' на рис. 7, b видно, что при такой замене структуры армирования максимальное значение интенсивности остаточных деформаций второго (окружного) семейства возрастает на 12%. Таким образом, замена «плоско»-перекрестной 2С-структуры армирования (рис. 2, а) на пространственную 4Б-структуру (рис. 2, Ъ) при динамическом нагружении тонкой длинной цилиндрической оболочки и фиксированном расходе волокон может приводить как к увеличению, так и уменьшению интенсивности остаточных деформаций компонентов композиции.
В связи с этим целесообразно проанализировать неупругое поведение таких цилиндрических КМ-оболочек при квазистатическом нагружении. Так, на рис. 4, b кривые 1"' и 2'" характеризуют зависимости wm(t), рассчитанные при разных структурах армирования и следующих параметрах нагрузки; Ртах = 3 МПа, im3X = 150 мс и imin = 300 мс (см. (15)). При принятом на рис. 4, b масштабе по оси ординат кривые 1'" и 2"' визуально не различаются, однако кривая 2'" все же лежит ниже кривой 1"'. В отличие от кривой 2, кривые 1"' и 2'" не осциллируют. Это свидетельствует о том, что указанное нагружение действительно можно рассматривать как квазистатическое.
На рис. 5 кривые 1"' и 2'" характеризуют зависимости w(x), рассчитанные при квазистатическом нагружении в момент времени t = 500 мс. Ординаты точек этих кривых уменьшены в 10 раз (см. правые участки кривых 1'" и 2'" на рис. 4, Ь). Из поведения кривых 1"' и 2"' на рис. 5 видно, что при квазистатическом нагружении цилиндрических КМ-оболочек внутренним давлением такие конструкции после снятия нагрузки не приобретают гофрированную остаточную форму, как это наблюдается при динамическом нагружении (см. кривые 1 и 2). На отрезке 0 < х < 1.6 м кривые 1"' и 2"' имеют практически горизонтальные участки. Это свидетельствует о том, что при квазистатическом нагружении в длинных цилиндрических КМ-оболоч-ках реализуется в основном безмоментное напряженное состояние и лишь в зонах краевых эффектов (при 1.6 < х ^ 2 м) возникает ярко выраженное моментное состояние. Сравнение правых участков всех кривых на рис. 5 показывает, что протяженность зон краевых эффектов как при динамическом (кривые 1 и 2), так и при квазистатическом (кривые 1'" и 2"') нагружениях примерно одна и та же.
Кривая 2"' на рис. 5 в основном лежит несколько ниже кривой 1"', т. е. при квазистатическом нагружении замена 2Б-структуры армирования (кривая 1'") на 4Б-структуру (кривая 2'") также приводит к уменьшению величины остаточного прогиба за пределами зон краевых эффектов. Однако это уменьшение мало (около 1.5%). Сравнение же интенсивности деформаций в компонентах композиции при этом показывает, что при t = 500 мс значения е^ = 3.14%, tiif = 1.11%) в случае 2Б-структуры армирования и £ш' = 4.08%, гт^ = 1.08% в случае использования 40-структуры. Следовательно, при квазистатическом нагружении тонкой длинной цилиндрической оболочки замена традиционной «плоско»-перекрестной структуры армирования на пространственную структуру приводит к увеличению интенсивности остаточных деформаций связующей матрицы на 29.9%) и к уменьшению аналогичной величины волокон второго семейства всего на 2.7%). Очевидно,
что при квазистатическом иагружеиии таких конструкций указанная замена структур армирования неэффективна.
Кривые V" и 2'" на рис. 4, Ъ и 5 получены по уточненной теории изгиба (М = 7). Расчеты, выполненные по теории Амбарцумяна, приводят практически к тем же результатам. Это относится и к определению деформированного состояния компонентов композиции оболочек. Поэтому неупругое поведение длинных тонких цилиндрических КМ-оболочек при их квазистатическом нагружении внутренним давлением вполне адекватно можно определять по традиционной неклассической теории Амбарцумяна.
Для расчета остаточных состояний (в частности, остаточных прогибов) динамически нагружаемых и неупруго деформируемых цилиндрических оболочек (в том числе и композитных, подкрепленных) ранее предлагалось использовать модель жесткопластического тела, что позволяет в ряде случаев получить достаточно простые и вполне обозримые аналитические решения [27]. Поэтому здесь целесообразно обсудить результаты расчетов по упру-гопластической модели, но приближенной к модели жесткопластического тела. С этой целью модули Юнга материалов композиции (см. таблицу) были увеличены: для связующей матрицы (к = 0) в 15 • 106 раз, а для арматуры (к ^ 1) в 106 раз. Для обеспечения устойчивости схемы «крест» шаг по времени А при этом был уменьшен в 10 раз. Динамический расчет проводился при исходных параметрах нагрузки: ршах = 8 МПа, £тах = 0.1 мс и ¿тт = 2 мс (см. (15)). Оказалось, что при таких входных данных величины
и>т имеют порядок Ю-9 м, а значения е^т (0 ^ к < И) — порядок 10_6 % как в случае 2В-структуры армирования, так и в случае 4В-структуры, т. е. пластичность в материалах композиции вообще не возникает Кроме того, в момент времени £ = £тах = 0.1 мс, когда внутреннее давление достигает своего максимального значения рШах> эпюра прогибов 1г(.е) имеет вид, аналогичный кривым 1"' и 2"' на рис. 5, т. е. качественно отличается от поведения остаточного прогиба, характеризуемого кривыми 1 и 2 на рис. 5. Следовательно, расчеты динамического неупругого поведения стеклопластиковых цилиндрических оболочек, выполненные по модели жесткопластического тела, ни количественно, ни качественно не позволяют определить и даже грубо оценить остаточные формы таких конструкций.
Приведенные выше результаты относятся к относительно тонким длинным цилиндрическим КМ-оболочкам (2h.fR = 1/50, Ь = AR). Дополнительные расчеты показывают, что для относительно толстых конструкций (2/г/ шт(й, Ь) = 1/10), а особенно для коротких цилиндрических КМ-оболочек, положительный эффект от замены «плоско»-перекрестной структуры армирования на пространственную 4В-структуру значительно возрастает, причем он начинает проявляться уже и при квазистатическом нагружении таких конструкций. В частности, величина наибольшего по модулю остаточного прогиба при такой замене может уменьшаться на несколько десятков процентов, а интенсивность остаточных деформаций в связующей матрице может уменьшаться даже в разы. Кроме того, с увеличением относительной толщины КМ-оболочки увеличивается различие между решениями, полученными по уточненной теории изгиба и традиционной неклассической теории Амбарцумяна. Разность наибольших по модулю значений остаточных прогибов при этом может превосходить 5%, а разность максимальных значений
интенсивности остаточных деформаций компонентов композиции может достигать даже сотен процентов как при «плоско»-перекрестном, так и при пространственном армировании конструкции.
Заключение. Разработана уточненная модель вязкоупругопластпческого деформирования гибких цилиндрических оболочек с пространственными структурами армирования, которая позволяет более точно, чем в рамках теории Амбарцумяна, рассчитывать остаточные перемещения и остаточные напряженно-деформированные состояния в компонентах композиции при действии интенсивных кратковременных нагрузок. Показано, что даже для относительно тонкой и длинной цилиндрической КМ-оболочки динамические расчеты, выполненные по теории Амбарцумяна, могут на несколько десятков процентов завышать максимальное значение и остаточную величину интенсивности деформаций компонентов композиции. Однако при квазистатическом нагружении таких неупруго деформируемых КМ-конструкций вполне оправданно можно использовать теорию Амбарцумяна.
Следовательно, для получения адекватных результатов динамических расчетов вязкоупругопластпческого поведения гибких КМ-конструкций целесообразно использовать теории изгиба более высокой точности, чем теория Амбарцумяна [11].
Проведенные расчеты продемонстрировали, что даже для тонких длинных оболочек из стеклопластика, нагруженных изнутри давлением взрывного типа, замена традиционной перекрестной структуры армирования по эквидистантным поверхностям (см. рис. 1, Ь и 2, а) на пространственную структуру (см. рис. 2, Ь) позволяет существенно уменьшить интенсивность остаточных деформаций связующей матрицы, а также снизить на 10—20% величину остаточных прогибов в центральной части КМ-оболочки. Наибольший остаточный прогиб, который возникает в зоне краевого эффекта, при такой замене структур армирования изменяется пренебрежимо мало. Однако осцилляции пространственно-армированной конструкции затухают быстрее, чем оболочки с «плоско»-перекрестной структурой армирования. При квазистатическом нагружении тонких длинных цилиндрических оболочек внутренним давлением такая замена структур армирования является неэффективной. В случаях же относительно толстых, и особенно коротких, цилиндрических оболочек положительный эффект от замены 2Б-структуры армирования (рис. 2, а) на 4Б-структуру (рис. 2, Ь) значительно возрастает, причем заметно проявляется и в случаях квазистатического нагружения таких конструкций.
Показано, что после неупругого динамического деформирования тонкие длинные стеклопластиковые цилиндрические КМ-оболочки приобретают гофрированную остаточную форму, причем образующиеся складки ориентированы в окружном направлении. При квазистатическом нагружении и неупругом деформировании таких конструкций складчатая остаточная форма не образуется (за пределами краевых эффектов цилиндрические КМ-оболочки при этом деформируются как безмоментные).
Продемонстрировано, что использование жесткойластической модели не позволяет адекватно рассчитывать форму остаточного прогиба цилиндрических стеклопластиковых оболочек при их динамическом нагружении. Такие расчеты на 6 порядков занижают максимальные значения прогибов и интен-сивностей деформаций компонентов композиции по сравнению с расчетами,
проведенными с использованием вязкоупругопластической модели. Поэтому модель жесткопластического тела нецелесообразно использовать для проведения динамических расчетов неупругого поведения тонкостенных стекло-пластиковых конструкций.
Конкурирующие интересы. У меня нет конфликта интересов в авторстве и публикации этой статьи.
Авторская ответственность. Я несу полную ответственность за предоставление окончательной версии рукописи в печать. Окончательная версия рукописи мною одобрена.
Финансирование. Работа выполнена в рамках государственного задания (№ госрегистрации 121030900260-6).
Библиографический список
1. Mouritz А. P., Gellert Е., Burchill P., ChaUis К. Review of advanced composite structures for naval ships and submarines// Compos. Struct., 2001. vol.63, no. 1. pp. 21-42. https:// doi.org/10.1016/S0263-8223(00)00175-6.
2. Reddy J. N. Mechanics of Laminated Composite Plates and Shells: Theory and Analysis. Boca Raton: CRC Press, 2003. xxiii+831 pp. https://doi.org/10.1201/bl2409.
3. Soutis C. Fibre reinforced composites in aircraft construction// Progress in Aerospace Sciences, 2005. vol.41, no. 2. pp. 143-151. https://doi.Org/10.1016/j.paerosci.2005.02. 004.
4. Gibson R. F. Principles of Composite Material Mechanics. Boca Raton: CRC Press, 2016. 700 pp. https://doi.org/10.1201/bl9626.
5. Vasiliev V. V., Morozov E. Advanced Mechanics of Composite Materials and Structural Elements. Elsever: Amsterdam, 2013. xii+412 pp. https://doi.org/10.1016/ C2011-0-07135-1.
6. Gill S. K., Gupta M., Satsangi P. S. Prediction of cutting forces in machining of unidirectional glass fiber reinforced plastics composite// Front. Mech. Eng., 2013. vol.8, no.2. pp. 187-200. https://doi.org/10.1007/sll465-013-0262-x.
7. Соломонов Ю. С., Георгиевский В. П., Недбай А. Я., Андрюшин В. А. Прик.яадные задачи механики композитных- цилиндрических, оболочек. М.: Физматлит, 2014. 408 с.
8. 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. https: //doi. org/ 10.1016/j.ijnonlinmec.2011.03.011.
9. Тарнопольский Ю. M., Жигун И. Г., Поляков В. А. Пространственно-армированные композиционные материалы: Справочник. М.: Машиностроение, 1987. 224 с.
10. Григоренко Я. М. Изотропные и анизотропные слоистые оболочки вращения переменной жесткости. Киев: Наукова думка, 1973. 228 с,
11. Амбардумян С. А. Общая теория анизотропных оболочек. М.: Наука, 1974. 446 с.
12. Reissner Е. On transverse vibrations of thin, shallow elastic shells// Quart. Appl. Math., 1955. vol.13, no. 2. pp. 169-176. https://doi.org/10.1090/qam/69715.
13. Богданович A. E. Нелинейные задачи динамики цилиндрических композитных оболочек. Рига: Зинатне, 1987. 295 с.
14. Абросимов Н. А., Баженов В. Г. Нелинейные задачи динамики композитных конструкций. Н. Новгород: Изд-во ННГУ, 2002. 400 с.
15. Янковский А. П. Моделирование вязкоупругопластического деформирования гибких армированных пластин с учетом слабого сопротивления поперечному сдвигу // Вычислительная механика сплошных сред, 2019. Т. 12, №1. С. 80-97. https://doi.org/ 10.7242/1999-6691/2019.12.1.8.
16. Янковский А. П. Моделирование вязкоупругопластического деформирования гибких пологих оболочек с пространственными структурами армирования // Вестн. Сам. гос.
техн. ун-та. Сер. Физ.-мат. науки, 2020. Т. 24, №3. С. 606-627. https://doi.org/ 10.14498/vsgtul709.
17. Андреев А. Н., Немировский Ю. В. Многослойные анизотропные оболочки и пластины. Изгиб, устойчивость и колебания. Новосибирск: Наука, 2001. 287 с.
18. Whitney J. М., Sun С. A higher order theory for extensional motion of laminated composites// J. Sound Vib., 1973. vol.30, no. 1. pp. 85-97. https://doi.org/10.1016/ s0022-460x(73)80052-5.
19. Lo К. H., Christensen R. M., Wu E. M. A high-order theory of plate deformation. Part 2: Laminated plates// J. Appl. Mech., 1977. vol. 44, no. 4. pp. 669-676. https://doi.org/ 10.1115/1.3424155.
20. Янковский А. П. Уточненная модель вязкоупругопластического деформирования армированных цилиндрических оболочек // Вестник Пермского национального исследовательского политехнического университета. Механика,, 2020. №1. С. 138-149. https://doi.org/10.15593/perm.mech/2020.1.11.
21. Композиционные материалы. Справочник / ред. Д. М. Карпинос. Киев: Наукова думка, 1985. 592 с.
22. Handbook of Composites/ ed. G. Lubin. New York: Van Nostrand Reinhold Company Inc., 1982. xi+786 pp. https://doi.org/10.1007/978-1-4615-7139-1.
23. Mohamed M. H., Bogdanovich A. E., Bickinson L. C., Singletary J. N., Lienhart R. R. A new generation of 3B woven fabric preforms and composites // Sampe J., 2001. vol. 37, no. 3. pp. 3-17.
24. Schuster J., Heider В., 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. https://doi.org/10.1007/sll029-009-9072-y.
25. Houlston R., BesRochers C. G. Nonlinear structural response of ship panels subjected to air blast loading// Comput. Struct., 1987. vol.26, no. 1-2. pp. 1—15. https://doi.org/ 10.1016/0045-7949(87)90232-X.
26. Хажинский Г. M. Модели деформирования и разрушения металлов. М.: Научный мир, 2011. 231 с.
27. Комаров К. Л., Немировский Ю. В. Динамик а жесткопластических элементов конструкций. Новосибирск: Наука, 1984. 234 с.
Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki
[J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2021, vol. 25, no. 2, pp. 343-364 ISSN: 2310-7081 (online), 1991-8615 (print) d https://doi.org/10.14498/vsgtul838
MSC: 74K20
A refined model of viscoelastic-plastic deformation of flexible spatially-reinforced cylindrical shells
© A. P. Yankovskii
Khristianovich Institute of Theoretical and Applied Mechanics,
Siberian Branch of the Russian Academy of Sciences,
4/1, Instit.utskaya. st., Novosibirsk, 630090, Russian Federation.
Abstract
A model of viscoelastic-plastic deformation of flexible circular cylindrical shells with spatial reinforcement structures is developed. The instant plastic behavior of the materials of the composition is determined by flow theory with isotropic hardening. The viscoelastic deformation of the components of the composition is described by the equations of the Maxwell-Boltzmann model. The geometric nonlinearity of the problem is taken into account in the Karman approximation. The used relations make it possible to calculate with varying degrees of accuracy the residual displacements of the points of the construction and the residual deformed state of the components of the composition. In this ease, a possible weak resistance of the reinforced shell to transverse shear is simulated. In the first approximation, the equations used, the initial and boundary conditions, are reduced to the formulas of the nonclassical Ambardzumyan theory.
The numerical solution of the formulated initial boundary-value problem is constructed according to the explicit "cross" scheme. Elastoplastic and viscoelastic-plastic dynamic deformation of thin fiberglass shells under the influence of internal pressure of an explosive type is investigated. Two reinforcement structures are considered:
1) orthogonal reinforcement in the longitudinal and circumferential directions;
2) spatial reinforcement in four directions.
It is shown that even for relatively thin composite shells the Ambardzumyan theory is unacceptable to obtain adequate results of calculations of their viscoelastic-plastic dynamic deformation. It has been demonstrated that a calculation according to the theory of elastoplastic deformation of reinforced structures does not allow even an approximate estimate of the residual states of composite shells after their dynamic loading. It is shown that even for a relatively thin and long cylindrical shell, the replacement of the traditional "flat'-cross-reinforcement structure with a spatial structure
Research Article
©® The content is published under the terms of the Creative Commons Attribution 4.0 International License (http: //creative commons . org/licenses/by/4.0/) Please cite this article in press as:
Yankovskii A. P. A refined model of viscoelastic-plastic deformation of flexible spatially-reinforced cylindrical shells, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2021, vol. 25, no. 2, pp. 343-364. https://doi.org/10.14498/vsgtul838 (In Russian).
Author's Details:
Andrei P. Yankovskii & https://orcid.org/0000-0002-2602-8357
Dr. Phys. & Math. Sci.; Leading Research Scientist; Lab. of Past Processes Physics;
e-mail: lab4nemir@rambler, ru
362
© Samara State Technical University
can significantly reduce the residual strain of the binder material. In cases of relatively thick and especially short shells, the positive effect of such a replacement of the reinforcement structures is manifested to a much greater extent.
Keywords: cylindrical shell, spatial reinforcement, reinforcement along equidistant surfaces, viscoelastic-plastic deformation, explosive-type loads, refined bending theories, Ambardzumyan theory, geometric nonlinearity, explicit numerical "cross" scheme type scheme.
Received: 8th December, 2020 / Revised: 12th March, 2021 J Accepted: 11th May, 2021 / First online: 5th June, 2021
Competing interests. I hereby declare that I have no conflict 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 carried out within the framework of a state assignment; state registration number — 121030900260-6.
References
1. Mouritz A. P., Gellert E., Burchill P., Challis K. Review of advanced composite structures for naval ships and submarines. Compos. Struct., 2001, vol.53, no. 1, pp. 21-42. https:// doi.org/10.1016/S0263-8223(00)00175-6.
2. Reddy J. N. Mechanics of Laminated Composite Plates and Shells: Theory and Analysis. Boca Raton, CRC Press, 2003, xxiii+831 pp. https://doi.org/10.1201/bl2409.
3. Soutis C. Fibre reinforced composites in aircraft construction, Progress in Aerospace Sciences, 2005, vol.41, no. 2, pp. 143-151. https://doi.Org/10.1016/j.paerosci.2005.02. 004.
4. Gibson R. F. Principles of Composite Material Mechanics. Boca Raton, CRC Press, 2016, 700 pp. https://doi.org/10.1201/bl9626.
5. Vasiliev V. V., Morozov E. Advanced Mechanics of Composite Materials and Structural Elements. Elsever, Amsterdam, 2013, xii+412 pp. https://doi.org/10.1016/ C2011-0-07135-1.
6. Gill S. K., Gupta M., Satsa.ngi P. S. Prediction of cutting forces in machining of unidirectional glass fiber reinforced plastics composite, Front. Mech. Eng., 2013, vol.8, no. 2, pp. 187-200. https://doi.org/10.1007/sll465-013-0262-x.
7. Solomonov Yu. S., Georgievskii V. P., Nedbai A. Ya,, Andryusliin V. A. Prikladnye zadachi mekhaniki kompozitnykh tsilindricheskikh obolochek [Applied Problems of Mechanics of Composite Cylindrical Shells], Moscow, Fizmatlit, 2014, 408 pp. (In Russian)
8. Kazanci Z. Bynamic response of composite sandwich plates subjected to time-dependent pressure pulses, Int. J. Nonlin. Mech., 2011, vol.46, no. 5, pp. 807—817. https://doi.org/ 10.1016/j.ijnonlinmec.2011.03.011.
9. 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)
10. Grigorenko Ya, M. Izotropnye i anizotropnye sloistye obolochki vrashcheniia peremennoi zhestkosti [Isotropic and Anisotropic Layered Shells of Revolution of Variable Stiffness], Kiev, Naukova dumka, 1973, 228 pp. (In Russian)
11. Ambartsumyan S. A. Obshchaia teoriia anizotropnykh obolochek [General Theory of Anisotropic Shells], Moscow, Nauka, 2002, 446 pp. (In Russian)
12. Reissner E. On transverse vibrations of thin, shallow elastic shells, Quart. Appl. Math., 1055, vol.13, no. 2, pp. 169-176. https://doi.org/10.1090/qam/69715.
13. Bogdanovich A. E. Nelineinye zadachi dinamiki tsilindricheskikh kompozitnykh obolochek [Non-Linear Dynamic Problems for Composite Cylindrical Shells], Riga, Zinatne, 1987, 295 pp. (In Russian)
14. Abrosimov N. A., Bazhenov V. G. Nelineinye zadachi dinamiki kompozitnykh konstrukt-sii [Nonlinear Problems of Dynamics of Composite Structures], Nizhni Novgorod, Nizhni Novgorod State Univ., 2002, 400 pp. (In Russian)
15. Yankovskii A. P. Modelling the viscoelastic-plastic deformation of flexible reinforced plates with account of weak resistance to transverse shear, Computational Continuum Mechanics, 2019, vol.12, no. 1, pp. 80-97 (In Russian), https://doi.org/10.7242/1999-6691/2019. 12.1.8.
16. Yankovskii A. P. Modeling of viscoelastoplastic deformation of flexible shallow shells with spatial-reinforcements structures, Vestn. Samar. Cos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.j, 2020, vol.24, no. 3, pp. 506-527 (In Russian), https://doi.org/10.14498/vsgtul709.
17. Andreev A. N., Nemirovskii Yu. V. Mnogosloinye anizotropnye obolochki i plastiny. Izgib, ustoichivost' i kolebaniia [Multilayer Anisotropic Shells and Plates. Bending, Stability, and Vibrations], Novosibirsk, Nauka, 2001, 287 pp. (In Russian)
18. Whitney J. M., Sun C. A higher order theory for extensional motion of laminated composites, J. Sound Vib., 1973, vol.30, no. 1, pp. 85-97. https://doi.org/10.1016/ s0022-460x(73)80052-5.
19. Lo K. H., Christensen R. M., Wu E. M. A high-order theory of plate deformation. Part 2: Laminated plates, J. Appl. Mech., 1977, vol.44, no.4, pp. 669-676. https://doi.org/ 10.1115/1.3424155.
20. Yankovskii A. P. The refined model of viscoelastic-plastic deformation of reinforced cylindrical shells, PNRPU Mechanics Bulletin, 2020, no. 1, pp. 138-149 (In Russian), https:// doi.org/10.15593/perm.mech/2020.1.11.
21. Kompozitsionnye materialy [Composite Materials], ed. D. M. Karpinos. Kiev, Naukova Dumka, 1985, 592 pp. (In Russian)
22. Handbook of Composites, ed. G. Lubin. New York, Van Nostrand Reinhold Company Inc., 1982, xi+786 pp. https://doi.org/10.1007/978-1-4615-7139-1.
23. Mohamed M. H., Bogdanovich A. E., Dickinson L. C., Singletary J. N., Lienhart R. R. A new generation of 3D woven fabric preforms and composites, Sampe J., 2001, vol.37, no. 3, pp. 3-17.
24. 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. https://doi.org/10.1007/sll029-009-9072-y.
25. Houlston R., DesRochers C. G. Nonlinear structural response of ship panels subjected to air blast loading, Comput. Struct, 1987, vol.26, no. 1-2, pp. 1-15. https://doi.org/10. 1016/0045-7949(87)90232-X.
26. Khazhinskii G. M. Modeli deformirovaniia i razrusheniia metallov [Deformation and Long-Term Strength of Metals], Moscow, Nauchnyi Mir, 2011, 231 pp. (In Russian)
27. Koma.rov K. L., Nemirovskii Yu. V. Dinamika zhestkoplasticheskikh elementov konstruktsii [Dynamics of Rigid-Plastic Structural Elements], Novosibirsk, Nauka, 1984, 234 pp. (In Russian)