CC BY
37
Янковский А.П. Моделирование динамики армированных пологих оболочек из нелинейно-упругих материалов // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2017. - № 2. -С. 226-245. DOI: 10.15593/perm.mech/2017.2.13
Yankovskii A.P. Modeling of the dynamics of reinforced shallow shells made of nonlinear elastic materials. PNRPU Mechanics Bulletin, 2017, no. 2, pp. 226-245. DOI: 10.15593/perm.mech/2017.2.13
ВЕСТНИК ПНИПУ. МЕХАНИКА
№ 2,2017 PNRPU MECHANICS BULLETIN
http://vestnik.pstu.ru/mechanics/about/ini7
001 10.15593/регш.шесЬ/2017.2.13 УДК 539.4
МОДЕЛИРОВАНИЕ ДИНАМИКИ АРМИРОВАННЫХ ПОЛОГИХ ОБОЛОЧЕК ИЗ НЕЛИНЕЙНО-УПРУГИХ МАТЕРИАЛОВ
А.П. Янковский
Институт теоретической и прикладной механики им. С.А. Христиановича СО РАН, Новосибирск, Россия
О СТАТЬЕ АННОТАЦИЯ
Сформулирована задача динамического поведения гибких армированных пологих оболочек, изготовленных из нелинейно-упругих материалов фаз композиции. Армирование осуществляется перекрестно по поверхностям, эквидистантным срединным поверхностям конструкций. Структура армирования по толщине оболочек является квазиоднородной. Геометрическая нелинейность задачи рассматривается в приближении Кармана, а ослабленное сопротивление композитных конструкций поперечному сдвигу описывается соотношениями неклассической теории Ред-ди, в рамках которой учитывается искривление поперечных нормалей к срединным поверхностям пологих оболочек. Численное интегрирование поставленной начально-краевой задачи осуществляется на основе метода шагов по времени с привлечением явной схемы типа «крест». Проведены конкретные расчеты динамического поведения относительно тонких и толстых пологих сферических оболочек и пластин, имеющих в плане кольцевую форму, при наличии жесткой внутренней шайбы и нагружении избыточным давлением, вызванным взрывной воздушной волной. Тонкостенные конструкции имеют жесткое закрепление по внешней кромке и осе-симметричное армирование по логарифмическим спиралям в плане. Изучается влияние углов армирования на податливость и напряженно-деформированное состояние в материалах фаз композиции гибких пластин и пологих оболочек. Выявлено, что на множестве рассматриваемых структур армирования рациональным является армирование по радиальным (меридиональным) направлениям, так как такая структура обеспечивает минимальную податливость и наименьшее напряженное состояние в материале связующей матрицы композиции. Показано, что в силу геометрической и физической нелинейности исследуемой задачи динамический отклик композитных пологих оболочек существенно зависит от того, к какой лицевой поверхности (выпуклой или вогнутой) прикладывается избыточное давление взрывного типа.
© ПНИПУ
Получена: 26 февраля 2017 г. Принята: 06 мая 2017 г. Опубликована: 30 июня 2017 г.
Ключевые слова:
пологие оболочки, структуры армирования, динамическое деформирование, нелинейная упругость, геометрическая нелинейность, теория Редди, взрывные нагрузки, схема « крест».
© Янковский Андрей Петрович - доктор физико-математических наук, ведущий научный сотруд ник, e-mail: [email protected], [email protected]
Andrey P. Yankovskii - Doctor of Physical and Mathematical Sciences, Leading Research Scientist e-mail: [email protected], [email protected]
MODELING THE DYNAMICS OF REINFORCED SHALLOW SHELLS MADE OF NONLINEAR ELASTIC MATERIALS
A.P. Yankovskii
Khristianovich Institute of Theoretical and Applied Mechanics SB RAS, Novosibirsk, Russian Federation
ABSTRACT
The problem of dynamic behavior of flexible reinforced shallow shells made of nonlinear elastic materials of the composition phases is formulated. The geometric nonlinearity of the problem considered in the Karman approximation, and the weakened resistance to shear is taken into account in the framework of non-classical Reddy theory. The numerical integration of the formulated initial-boundary value problem is carried out on the basis of the method of steps in time with the involvement of an explicit "cross" scheme. Specific calculations of the dynamic behavior are carried out for relatively thin and thick shallow spherical shells and plates having an annular shape in plan and with the rigid inner insertion which are under the pressure caused by the blast air wave. Thin-walled structures are clamped on the outer edge and axially reinforced by logarithmic spirals in the plan. The influence of the angles of reinforcement on the flexibility and stressstrain state in the materials of the phase of the composition of elastic plates and shallow shells is studied. It is obtained that on the set of considered structures of reinforcement, the reinforcement in the radial (meridional) directions is rational, as this structure provides the minimal flexibility and the minimal stressed state in the material of the binder matrix of compositions. It is shown that because of the geometric and physical nonlinearity of the studied problem, the dynamic response of composite shallow shells essentially depends on which front surface (convex or concave) the overpressure of the explosive type is applied.
© PNRPU
Введение
Тонкостенные элементы конструкций типа пластин и оболочек из композиционных материалов находят все более широкое применение в инженерной практике [1-3]. Активное использование оболочечных конструкций объясняется их высокой эффективностью, так как они обладают легкостью в сочетании с высокой прочностью [4, 5]. Особенно ярко эта эффективность проявляется при внедрении в инженерные изделия композиционных оболочек волокнистой структуры [5-8], поэтому актуальной является проблема адекватного расчета механического поведения композитных панелей и оболочек как при квазистатическом, так и динамическом их нагружении [4, 5, 9, 10].
На сегодняшний день подавляющее число публикаций по проблемам механического поведения анизотропных и композитных конструкций базируется на предположении о линейно-упругом поведении материалов компонентов композиции [4-11]. Однако современные композиты могут быть изготовлены из материалов фаз, проявляющих неупругие или нелинейно-упругие свойства [12-15]. Анализ справочных данных [14, 15] показывает, что при расчетах в рамках линейно-упругого поведения материалов всех компонентов композиции уровень напряжений, возникающих при этом в высокопрочных армирующих волокнах, не превосходит, как правило, 10-20 % от предела их прочности, т.е. уникальные прочностные свойства арматуры используются весьма неэффективно. При неупругом или нелинейно-упругом деформировании материала связующей матрицы композита несущую способность волокон можно использовать практически полностью [16].
Известно, что при армировании композитных пластин и оболочек высокомодульными волокнами необходимо при изучении их механического поведения учитывать ослабленное сопротивление поперечным сдвигам [4, 5, 9-11]. Как правило, этот учет проводится в рамках теории Рейсснера [9, 10, 17, 18] или теории Редди [4, 5, 11, 16, 19, 20] (в русскоязычной литературе их называют первым и вторым вариантами теории Тимошенко [19]).
ARTICLE INFO
Received: 26 February 2017 Accepted: 06 May 2017 Published: 30 June 2017
Keywords:
shallow shell, structure of reinforcement, dynamic deformation, nonlinear elasticity, geometric nonlinearity, Reddy theory, explosive loads, "cross" scheme.
В монографии [9] на основе гипотез Рейсснера исследовалось динамическое поведение гибких тонкостенных элементов композитных конструкций, изготовленных из линейно-упругих материалов компонентов композиции. В работе [21] проводилось сравнение динамического поведения гибких композитных балок-стенок из линейно-упругих материалов, рассчитанного по классической теории и двум вариантам теории Тимошенко, причем показано, что результаты расчетов по теориям Тимошенко могут существенно различаться. В [16] на базе гипотез теории Редди в геометрически линейной постановке построена модель нелинейно-упругого деформирования армированных пластин, нагруженных статически. Однако динамическое поведение тонкостенных композитных конструкций из таких материалов до настоящего времени не изучалось.
Для численного интегрирования динамических задач механики деформируемого твердого тела используют как явные, так и неявные схемы шагов по времени. Из явных численных схем наиболее часто используются схемы типа «крест» [9, 21-23], так как они достаточно просто поддаются алгоритмизации [9], а из неявных схем наибольшее распространение получили методы Ньюмарка [24, 25].
В связи с вышеизложенным настоящее исследование посвящено моделированию в рамках теории Редди динамического поведения гибких армированных пологих оболочек, изготовленных из нелинейно-упругих материалов компонентов композиции. Моделирование производится на базе метода шагов по времени с привлечением явной численной схемы типа «крест».
1. Постановка задачи
Рассматривается пологая оболочка толщиной 2h, с которой связана криволинейная ортогональная система координат xi так, что отсчетная поверхность Ox1 x2 (x3 = 0)
совпадает со срединной поверхностью оболочки (< h); координатные линии x1 и x2 совпадают с линиями главной кривизны отсчетной поверхности; ось x3 направлена
по толщине оболочки. Структура армирования конструкции в направлении x3 квазиод-
нородна. На рис. 1 изображен малый элемент пологой оболочки, причем искривленность отсчетной поверхности в силу ее малости не изображена.
Ослабленное сопротивление тонкостенного армированного элемента поперечному сдвигу будем учитывать в рамках теории Редди [4, 5, 11, 16, 19, 20], так как с точки зрения математического моделирования она является более точной [21], чем теория Рейсснера [9, 10, 17, 18]. Согласно теории Редди сдвиговые деформации в поперечном направлении вг3 аппроксимируем по переменной x3 так [4]:
xl [>0000000000[
Рис. 1. Элемент армированной композитной пологой оболочки или пластины Fig. l. Element of the reinforced composite shallow shell or a plate
(, г ) =
х + И
.(+)
(*, x)-
х - И
( x) +
1 2 2 И — х
3- в0
■2 ¡3
(*, x),
2И "" 4 ' 7 2И "" 4 ' 7 И2 ^ '' (1)
x е О, <И, * > *о, x = { Х2}, г = { Х2, Х3}, г = 1,2,
где в(3г), в03 - функции только двух пространственных переменных х1, х2 и времени *, подлежащие определению и имеющие следующий смысл: в(3) - деформации поперечных сдвигов на верхней (х3 = И) и нижней (х3 = — И ) лицевых поверхностях оболочки, при в(3) = 0 функции в03 определяют деформации поперечных сдвигов на срединной поверхности х3 = 0 ; *0 - начальный момент времени; О - область, занимаемая пологой оболочкой в плане.
Согласно кинематической гипотезе теории Редди изменяемостью перемещений точек оболочки в направлении х3 пренебрегаем, т.е. принимаем
и (*, г) = w(*, x), x е О, < И, * > *0,
(2)
где — - прогиб точек отсчетной поверхности х3 = 0 .
На основании дифференциальных соотношений, связывающих деформации с перемещениями, с учетом (1) и (2) в приближении Кармана [4] для пологой оболочки получим
и, (*, г) = и, (*, х^) — х3^н +
2 х
(
И2
2 Л
И2 — ^ 3
в 03 + х1 (4 + И |в(3+) — ^ [ Х^—ИИ—), , = 1,2; (3)
И V 2
И V 2
г) = 5,Щ — х352— +
2 х,
И2
И2 — ^
2
5,4 + ^ ( ^ + и 15,.в(3+) — Х3 ( ^ — И кв(3—) +— + г г3 И V 2 | , ,3 И V 2 I ' ,3 Я
+ 2 (5г— )
В12 (*, Г)= 2(51Щ2 +52Щ1 ) — Х35152— + И2
С х 2 ^
И2 — х 3
у
(23 +52В°°3) + 2 х (4)
х^ + И|(51в2+) + 52в(+)) — 2и— И|(51в2—) +52в(—)) +)5—5— xеО, |х3|<И, г = 1,2,
где иг - перемещения точек оболочки в тангенциальных направлениях хг; иг - то же для точек срединной поверхности (х3 = 0 ); Я - главные радиусы кривизны отсчетной поверхности; 5,. - оператор дифференцирования по переменной х (¡ = 1, 2).
Таким образом, соотношения (1), (3) и (4) содержат следующие неизвестные функции: —, и,, в(3), в03 (/ = 1, 2), зависящие от двух пространственных переменных х1, х2 и времени *.
Согласно принципу Даламбера, уравнения динамического равновесия пологих оболочек в приближении Кармана с учетом равенств (2) и (3) имеют вид [4]
22 £5; (( — ^35,—) + Я + X, (*, x) + о(+) — о(—) = 2ИрЩ + ^ (в(+) — в)),
}=1
£
1=1
5,
Е +£ Е 5
—
Я
(5)
+ Х3 (*, x) + а3+) — а33) = 2Ир—, % = 1, 2;
Усткоувкп А.Р. /РЫЯРиМесНап¡ев БиПеШ 2 (2017) 226-245
£ам,-ъз + т(, х)н(о(з)+о(-))=3н3р{8в0 -ер]+2н3р(в£)+(-)),
1=1
х е О, I > I = 1,2,
где, согласно [14]
К К
р = ро®о +Ерк®к, ®о =1 -Еюк, (7)
к=1 к=1
р0, рк - объемная плотность материала связующего и арматуры к-го семейства; шк -плотность армирования волокнами к-го семейства; К - количество семейств арматуры; а(з) = oi3 (, х, ± Н) - заданные напряжения (поверхностные нагрузки) на верхней (+) и нижней (-) лицевых поверхностях оболочки; Х1 - приведенные внешние нагрузки, действующие в направлении х1 и порожденные объемными силами (I = 1, 3 ); т1 - распределенные внешние моменты от объемных сил; /¡3 - поперечные силы; - мембранные усилия; М1 - изгибающие и крутящие моменты; точка означает частное дифференцирование по времени
Материалы всех фаз композиции оболочки предполагаются изотропными и однородными, а их нелинейно-упругое поведение описывается тензорно-квазилинейными определяющими соотношениями [16, 26, 27]
о«-5^ = gik) (в(к), 80к)))-5^) (I, 1 = 1, 2, 3),
3 (8)
о0к) = g0k) (в(к), е0к))к), о(к) = 3) (в(к), е0к))в(к), 0 < к < К,
где о(к), в(к) - компоненты тензоров напряжений и деформаций в материале к-го компонента композиции (индекс к = 0 соответствует связующему, а индекс к = 1, 2, ..., К - арматуре к-го семейства); о0к), 80к) - среднее нормальное напряжение и средняя линейная деформация в к-м компоненте композиции; о(к), в(к) - интенсивности напряжений и деформаций в тех же материалах (см. (4) в [27]); g(к) (в(к), 80к)) - функция, известная из эксперимента и являющаяся при умножении ее на 3/2 коэффициентом пропорциональности между о(к) и в(к); g(0k) (в(к), в0к)) - утроенный секущий модуль объемного расширения
материала к-го компонента композиции, известный из эксперимента; 5, - символ Кроне-
кера. Зависимости функций g(k) и g0k) в соотношениях (8) от аргумента в0к) позволяют учесть эффект разносопротивляемости материала к-й фазы композиции [26, 28].
Равенства (8) запишем в матричной форме [19]:
о к = Ак е к, к = 0,1,2,..., К, (9)
где
о го(к) о(к) о(к) о(к) о(к) о(к)) в
°к 1°11 > 22 > °33 > 23 > °31 > °12 ( ' к
/в( к) в( к) в(к) в( к) в(к) в( к) Г (10)
\Ь11 ' 22 > Ь33 ' 23 > Ь31 ' Ь12 ) '
Лк = (аг(,к}) - симметричные 6 х 6-матрицы, ненулевые компоненты которых согласно (8)
с учетом (10) определяются по формулам (24) из [27]; звездочка в верхнем индексе означает операцию транспонирования.
Как и в [27], предполагаем, что определяющие соотношения (8) (а значит, и (9) с учетом (10)) удовлетворяют достаточным условиям сходимости метода последовательных приближений [26, с. 199], аналогичного методу переменных параметров упругости [29]. Поэтому далее в фиксированный момент времени I равенства (8) и (9) на текущей итерации считаем линеаризованными. В этом случае матричные соотношения (9) формально полностью совпадают с равенствами (22) из [27] (где ядра ползучести нужно принять тождественно равными нулю). Повторяя соответствующие рассуждения из [27], на данной итерации получим линеаризованные определяющие соотношения для армированной среды, записанные в матричной форме:
а = Лв, (11)
где а, в - шестикомпонентные векторы-столбцы средних напряжений и деформаций в композиции, имеющие структуру, аналогичную (10); Л = ( ау ) - известная 6 х 6-матрица,
которую на текущей итерации можно трактовать как матрицу эффективных жесткостей армированной среды. Элементы ау матрицы Л в (11) определяются из матричных соотношений (27) в [27] и зависят от деформированного состояния композиции в у, известного из решения на предыдущей
итерации, и от параметров армирования: плотностей шк и направлений армирования, определяемых углами сферической системой координат 9к и фк . На рис. 2 согласно [27] изображена взаимная ориентация глобальной х и локальной х(к)
(г = 1, 3 ), связанной с волокном к-го семейства, систем координат в случае пространственного армирования.
Пусть из решения начально-краевой задачи для композитного тела (в частности, пологой оболочки), механическое поведение которого описывается определяющим соотношением (11), в текущий момент времени I на данной итерации известны осреднен-ные деформации в, тогда из матричных равенств (31) и (32) в [27] можно вычислить и деформации вк в к-й фазе композиции (9) и (10). После чего согласно методу последовательных приближений можно уточнить значения коэффициентов g¡(k) и g0k) в равенствах (8) (или, что то же самое, значения элементов матриц Лк в равенствах (9))
и по приведенному выше алгоритму (подробнее см.: [27]) можно определить следующее приближение решения и т. д., пока итерационный процесс не сойдется с требуемой точностью.
Рис. 2. Локальная система координат, связанная с волокном k-го семейства Fig. 2. The local system of coordinates connected to a fiber of k-th family
В случае структуры армирования, изображенной на рис. 1 (9к =л:/2 (см. рис. 2))
и рассматриваемой в настоящем исследовании, матрица А в соотношении (11) имеет такую структуру, что ее компоненты в 4-й и 5-й строках и в 4-м и 5-м столбцах, кроме диагональных, тождественно равны нулю, т.е. согласно обозначениям, аналогичным (10), имеют место равенства
С23 = С44В23, С13 = С55„13. (12)
На основании силовой гипотезы теории Редди с33 (г, г) « 0 [4, 5, 11, 16, 19-21] систему
шести линейных алгебраических соотношений (11) с учетом (12) можно преобразовать к традиционному виду, исключив из рассмотрения приближение для деформации в33:
о = Б е, с 3 = Ъг в^ г =1,2 (13)
где
° = {с11, С22, С12} , ё={в11, В22, В12} , В33 СС31В11 + С32В22 )/С33 ; (14)
Ъ] = Сг] - аг3а3 у/с33 , Ъг3 = аг6, Ъ3 ] = С6 ] , Ъ33 = ^ Ъ1 = ^ Ъ2 = а44, К ] = 1,2 (15)
а. - элементы матрицы А в равенстве (11); Ъ. - элементы 3 х 3-матрицы Б; с., в. - компоненты средних напряжений и деформаций в композиции (элементы векторов-столбцов ст и в в равенстве (11)). Согласно (15) величины Ъ. и Ъ в соотношениях (13) в текущий
момент времени г предполагаются известными из решения рассматриваемой задачи на предыдущей итерации метода последовательных приближений.
Для удобства дальнейшего изложения введем в рассмотрение новые неизвестные функции
8
У,3 (, Х)= 54 -дгw,
отсюда
5
вг°3 (, х) = -(уй + г = 1,2. (16)
Уравнения движения (6) с учетом (16) принимают вид
X д]М. 3 +Щ (, X ) + к (с(з) + 0(3) ) = 3 к3р(уг 3 + ё(3) + В«), хеа, г >о, / = 1,2. (17)
]=1 3
Равенства (5) и (17) образуют систему уравнений движения рассматриваемой гибкой пологой оболочки, к которым следует добавить четыре силовых граничных условия на лицевых поверхностях (см. (13) с учетом выражений (1) и (15)):
г( ±)
(, X) = Ъ (г, X, ± к)в(3±) (г, X), г = 1,2, (18)
где левые части известны. (В рамках классической теории [4, 24, 28] и теории Рейсснера [9, 10, 17, 18] граничные условия (18) удовлетворить не удается.) Из равенств (18) получаем
в(±) (, X) = с(±) (г, X)/Ъг (г, X, ± к), г = 1, 2. (19)
Следовательно, в данный момент времени ^ текущие приближения деформаций поперечного сдвига на лицевых поверхностях пологой оболочки можно считать известными, так как приближения «эффективных жесткостей» Ъ согласно равенствам (15) предполагаются уже известными из решения, полученного на предыдущей итерации метода последовательных приближений. В дальнейшем функции в(±) считаем известными из (19)
в каждый рассматриваемый момент времени.
Используя соотношения (13) с учетом выражений (1), (4) и (14)-(16), на текущей итерации можем определить все внутренние силовые факторы в пологой оболочке:
п 2
F =| аdx3 = ^
-к г=1
Bгди -I ^ -4^ ]д>+4^д,у¿3+2c(+)aгв(í) -2c(-)aгв(з) +
— и 1 — / ч
+2 в' (ди)
+ ^ (д2^1 + д1и2 ) - 2 ^ A3 - "4 C3 д1д2и + "4 C3 (д2У13 + д1У23 ) +
+2C3+) (д2в(;) + д1в2з)) - 2C3-) (д2в(-) + д^) + д^д2и,
к2 M =| ах3йХ3 = ^
Aг-ди -I Eг --4Hl ]д> + -4Hгдгуг3 + 2H(+)дгв(3-) -в(-) + (20)
а и 1 \2
+^+2 (д
+ ^ (д2и1 + д1и2 ) - 2 I E3 - "4 H3 }д1д2и + 4 ^ (д2У13 + д1У23 ) +
+2^+) (д 2в(;) + д^«) - 2H3-) (д 2в(-) + д^«) + -3 дид2 и,
V3 = } агзdХз = 8dl (уг3 + ди) + а(+) - а
- аг(-)в(-), г = 1,2, x е G, I >
где согласно (14)
F-{^11, ^22, , M - (Мп, М22, М12},
(21)
на основании
Aг., Bг., Cг., C(+), Eг., ^, H(+) - трехкомпонентные векторы-столбцы, элементы которых совпадают с компонентами г-х столбцов 3 х 3 -матриц Л = (а. ), В = (Ъуг), С = (ер ),
С(±) = (с(±), Е = (.), н = (кг), И(±) = (к(±) (г, у = 1,3), определяемых соотношений (1), (4), (13) и (20) так:
к к к I 2 ] к
Л Вх3аХ3, В ВаХ3, С вЦ к -аХ3, Е Вх32аХ3,
- к - к - к к V 3 У - к
к „2 (
с(±) - Г В—I^ ± н\ аХ3, И - Г ВХ2 ) 2к V 2 I 3 ) к2
- к
-к
к2
2
к2 - ^
V 3 У
к 2 . Х^ | Х-}
аХ3, н(±)-ГВ±н]аХ, (22)
к к2 - Х32
аг - Гъг аХ3, а(±) - Гъг Х3±каХ3, г = 1,2.
г ^ г г,2 3' г J г пи 3
- к
к
2к
В - 3 х 3 -матрица из соотношений (13) (см. (15)).
Для однозначного интегрирования рассматриваемой динамической задачи необходимо использовать силовые [9, 10, 17-20]
= F
1 n3>
2 2 f 2 ^ Z n (( - F^w) =F Z n F13 +Z Fдj
1=1 1=1 V j=1 У
M„n,2 + М22П22 + 2М12П1П2 = Mnn, (m22 - Mu )nin2 + M12 (2 - П22) = Mnx, (23) n1 = cosa, n2 = sin a, x ёГ, t > t0, i = 1, 2 и кинематические (см. (2) и (3))
w(t, x) = w.(t, x), xёГ, t > to; (24)
U (t, r) = U„. (t, r), xёГ, ^ <h, t > to, i = 1,2 (25)
граничные условия, а также начальные условия в момент времени t = t0:
w ( x) = uo3 (x), w (^ x) = vo3 (x); (26)
U (to, r) = Uoi (r), U (to, r) = Voi (r), x e G, |х,| < h, i = 1,2. (27)
Здесь Г - контур, ограничивающий область G, занимаемую пологой оболочкой в плане; Fi - заданные на Г мембранные силы, действующие по направлениям xi; Fn3 - заданная
на контуре Г поперечная сила; Mnn, Mm - заданные на контуре Г изгибающий и крутящий моменты; w„ - заданный на Г прогиб; U„i - заданные на торцевой поверхности оболочки перемещения в тангенциальных направлениях xi ; a - угол, задающий направление внешней нормали к контуру Г; uo3, vo3, Uoi, Voi (i = 1, 2) - заданные в начальный момент времени to перемещения и скорости точек пологой оболочки.
Согласно выражениям (3) граничные (25) и начальные (27) условия нельзя удовлетворить в каждой точке по толщине оболочки при произвольных зависимостях функций Usi,
Uoi, Voi от поперечной координаты x3. В силу этого обстоятельства кинематические граничные условия (25) и начальные условия (27) подобно силовым граничным условиям (23) будем удовлетворять в интегральном смысле, т. е. проинтегрируем равенства (25) и (27) по толщине оболочки с весами 1 и x3, тогда с учетом соотношений (3) и (16) получим
h2 h2 2hü ( x) + у в« ( x)-—е(3-) ( x ) = u.i (t, x),
2
3h3 [y,3 (t, x) + в(3+)(t, x) + в(3-)(t, x)] = ü, (t, x), xёГ, t > to, i = 1,2;
h2 h2 -2hü (to, x)+ у 4 ) (to, x)-y в(3-) (to, x) = uor (x),
2h3 [yi3 ( x) + s(3+) ( x) + s(3-) (to, x)] = ñ0i (x)
2hü (to, x)+ y[s(3+) (to, x)-B(3-) (to, x)] = Voi (x), 2
-h3 [yi3 (to, x) + é(3+) ( x) + s(3) (to, x)] = voi (x), x e G, i = 1,2,
(28)
где
h h h
ич (t, x) = J U,t (t, r)dx3, йщ (t, x) = J U,t (t, r)X3dx3, u01 (x) = J U0i (r)dx3,
-h -h -h (30)
h h h v '
U0i
(x) = J u0i (r)x3dx3. ^ (x) = J V0i (r)dxз, ^ (x) = J V0i (r)x3dx3. i = 1, 2.
Таким образом, для однозначного интегрирования исследуемой динамической задачи в каждой точке области G в момент времени t0 должны быть заданы начальные условия
(26) и (29) с учетом (30), а в каждой точке контура Г - силовые граничные условия (23) или кинематические граничные условия (24) и (28) с учетом (30). Возможно задать и пять смешанных из (23), (24) и (28) граничных условий, например в случае свободного опира-ния кромки пологой оболочки.
После подстановки внутренних силовых факторов (20) с учетом (21) в уравнения движения (5), (17) и в силовые граничные условия (23) получим систему разрешающих уравнений и соответствующие ей силовые граничные условия, записанные в кинематических переменных w, ui и yi3 (см. (16)). Эти соотношения громоздки, поэтому не будем их
приводить, помня о том, что далее предполагается использовать явную схему численного интегрирования рассматриваемой динамической задачи, для реализации которой (схемы) не требуется использовать уравнения движения и силовые граничные условия, выраженные явно через зависимые переменные w, ui и уi3 (i = 1, 2).
2. Метод расчета
Для численного интегрирования сформулированной в п. 1 начально-краевой задачи применим метод шагов по времени [9, 21-25, 27], т.е. будем разыскивать решение в дискретные моменты времени tn+1 = tn + т (n = 0, 1, 2...), где т = const > 0 - шаг по времени.
Предполагается, что в моменты времени tm известны значения функций
m m m m
W (x )= W (tm , x) U (x )= U (tm , x) Уг3 (x) = Уг 3 (tm , x), ^ (x) = j (tm , x), (31)
aa (r) = akl (tm, r), i = 1,2, j = 1,3, k, l = 1,6, m = n-1, n, x e G, < h.
Здесь akl - компоненты 6 x 6 -матрицы A в равенстве (11)), тогда по формулам (20) с учетом (15), (19), (21) и (22) в момент времени tn можем определить значения всех внутренних силовых факторов Fj, Mij и Fi3, входящих в уравнения движения (5), (17) и силовые
граничные условия (23).
Далее для аппроксимации встречающихся производных по времени будем использовать центральные конечные разности [9, 21-23], что позволяет получить явную численную схему для интегрирования рассматриваемой динамической задачи. Согласно этому конечно-разностные аналоги уравнений (5) и (17) с учетом обозначений, аналогичных (31), примут вид
2hp (n+1 „ n n-1 ^ Л
1 w - 2 w+ w I = ^
Т2 V J 7=1
F13 + Z Fj д Jw
j=1
+R—1 fu
+ X3 + ^3) — ^3),
2кр
т2
2к3р
п+1 п п-1Л 2 / т п п\ п п п п
»•-2<+»■ I = £г, I Р,,-з г, н | + Л,-1 -+X, + о(+ >-о»
У '=' 1 ^ (32)
Сп+1 п п-А 2 п п п { п п Л
3т2 V У ,=1
У, -2У, + У, 1 = ЕгЛ --3 + тг + к
о(+)+о(-)
3 ^ 3
хеС, , = 1, 2, и = 1, 2, 3...,
У
где
у, (Г, х)-у,3 + в(+) + в(-), пгв(, х)и, + -к ( - в(-)), , = 1, 2. (33)
Так как правые части в уравнениях (32) известны, из них с учетом (19), (31) и (33)
п+1 п+1 п+1
можем определить неизвестные функции н, п, и у,3 в следующий момент времени 1п+1. Строго говоря, при I = 1п к системе (32) следует присоединить необходимые граничные
п+1 п+1
условия (23), (24) и (28) с учетом обозначений, аналогичных (31). При известных н , п, и
п+1
у,3 по формулам (1), (4) и (14) с учетом соотношений (16) и (19) можем вычислить сред-
п+1
ние деформации композиции в,, во всех точках пологой оболочки при I = tn+1. Используя
п+1
известные функции в,, (,', , = 1, 2, 3), по формулам (31) и (32) из [27] в момент времени tn+1 можно определить текущие приближения для деформаций фаз композиции, а затем по формулам (24) и (27) из [27] можно вычислить и компоненты 6 х 6 -матрицы А в равенстве (11). При этом в каждой точке пологой оболочки требуется реализовать итерационный процесс, аналогичный методу переменных параметров упругости. Для начала такого процесса можно взять известные при t = tn матрицы, используемые в соотношениях (31), (32) и (27) из [27]. Это вполне оправданно при достаточно малом шаге по времени т .
После окончания такого итерационного процесса известными при т = п и т = п +1 становятся все функции, указанные в (31), поэтому по предложенной схеме можно построить решение исследуемой задачи в следующий момент времени tn+2 и т.д.
Структура левых частей уравнений (32) с учетом выражений (33) свидетельствует о том, что для начала расчетов по описанной выше явной численной схеме необходимо
0 0 0
знать значения как функций н, п, и у,3, известных из начальных условий (26) и (29)
1 1 1
с учетом (19) и (30), так и функций н , п, и у,3 (см. (32) и (33) при п = 1 и t = Значения этих функций получим, используя формулу Тейлора и учитывая начальные условия (26), (29) и уравнения движения (5) и (17) при t = t0 [22, 23]:
1 0 0 т2 0
п, (х) = п, (х)+тп, (х)+-2п, (х)+0 (т3) * 0
2
у,3 (х) = у, 3 (х) + ту, 3 (х) + - у, 3 (х) + О (т3) * 0, (34)
1 0 0 т2 0 , .
н (х) = н (х) + т Н (х) +—Н (х) + О (т3 )* 0, , = 1,2, х е С,
где приближенные равенства выполняются с точность порядка т3, а нулевые правые части соответствуют естественному состоянию, когда в начальный момент времени t0 пологая оболочка покоится и внешние нагрузки отсутствуют.
Если область О, занимаемая пологой оболочкой в плане, является прямоугольной, то, аппроксимируя в уравнениях (32) производные дг (•) по пространственным переменным хг
их конечно-разностными аналогами, получим в конечном счете явную численную схему типа «крест» [9, 22, 23]. Если же область О является неканонической, то дискретизацию уравнений (32) по пространственным переменным хг (г = 1, 2) можно провести на основе
вариационно-разностного подхода, использованного в [9], где показано, что для тонкостенных элементов конструкций необходимые условия устойчивости схемы типа «крест» определяются неравенствами, вытекающими из условия Куранта-Фридрихса-Леви [22]:
Ае /т> ае, 2И/т> а3, (35)
где
ае = VЕее / Р, а3 = V Е33 / Р, (36)
р определено первым соотношением (7); Ае - максимальный размер ячейки дискретизации области О в направлении вектора е, параллельного отсчетной поверхности пологой оболочки; Еее, Е33 - эффективные модули упругости первого рода в направлениях е и х3,
которые определяются с использованием компонент матрицы эффективных податливо-стей композиции, обратной матрице А (см. равенство (11) при линейно-упругом поведении материалов фаз композиции оболочки) [19].
Так как первое неравенство (35) должно выполняться при любом направлении е, для неоднородных, а тем более и криволинейных структур армирования пологих оболочек выбор шага по времени т из условий (35) представляет самостоятельную трудоемкую задачу. Однако нетрудно показать, что если для каждого компонента композиции выполняются неравенства, аналогичные (35), то автоматически выполняются и условия (35) для армированной пологой оболочки. В силу этого далее используем именно такой необходимый критерий устойчивости предложенной явной численной схемы.
Используя известные формулы перехода от декартовой прямоугольной системы координат к цилиндрической [30], полученные выше уравнения и граничные условия можно записать в полярной системе координат. Это удобно, если пологая оболочка является, например, сферической (что часто встречается в инженерной практике [31]), а область О представляет собой круг, кольцо или сектор круга либо кольца.
3. Обсуждение результатов расчетов
В качестве конкретных примеров рассмотрим динамическое поведение пологих сферических оболочек (Е1 = Е = Я ), для которых область О является кольцом, ограниченным контурами радиусов г0 = 10 см и г1 = 60 см. Оболочки имеют постоянную толщину 2Н = 1 см или 2Н = 5 см и разные величины подъема стрелы / полюсных точек над их внешними опорными контурами (х1 = г1, где х1 - полярный радиус). Из геометрических соображений получаем
R =
£±fi 2f
0 < f <-r. 5 1
Объемными нагрузками пренебрегаем, т.е. X, - 0 и т, - 0, I = 1, 3 , , = 1, 2 (см. (5), (6), (17) и (32)). Касательные нагрузки на лицевых поверхностях оболочки отсутствуют
(огз - 0), поэтому из (19) следует в(3) - 0, I = 1, 2. Конструкции нагружаются фронтальной нагрузкой, порожденной воздушной взрывной волной (см. (5)) [24]:
(+) „(-) _
°33 °33
- р (t )=
I Pmaxt / tmax '
0 < t < t
/maxeXP [-a( - tmax )] , t > tm
где
a = -ln(0,01)/(t. -1 )> 0, t. »t ,
V ? / ' \ mm max / 5 mm max'
(38)
(39)
tmax - момент времени, в который нагрузка р (t) достигает максимального по модулю значения pmax; tmin - момент времени, в который нагрузка р (t) по модулю становится пренебрежимо малой по сравнению с |pmax| (так, формула (39) соответствует случаю Р (tmin ) = 0,01 pmax). Согласно экспериментальным данным, приведенным в [24], в расчетах примем tmax = 0,1 мс и tmin = 200 мс.
Предполагается, что пологие оболочки усилены осесимметрично по логарифмическим спиралям двумя (K = 2) семействами волокон с одинаковыми плотностями армирования (рис. 3), т.е.
®1 (Х1 ) = ®2 (Х1 ) ф1 (х1 ) = -ф2 (х1 ) = ф = const, (40)
r0 < X1 < ri.
Так как волокна имеют постоянные поперечные сечения [14], интенсивности армирования &к в осесимметричном случае
п _ Т, ~ при укладке волокон по логарифмическим
Рис. 3. Кольцевая область, занимаемая пла- r J r т
стиной или пологой сферической оболочкой спиралям при любых углах Ф определяются
в плане, с осесимметричными траекториями так [16]-
армирования по логарифмическим спиралям ^ (х ) = r ^ / х
Fig. 3. The ring area occupied with a plate or a shal- 0 0k 1
low spherical shell in the plan, with axial-symmetric ®0 k = ®k (Г), (41)
trajectories of reinforcement on logarithmic spirals ю = ю k = 12
где ш0к - плотность армирования волокнами к-го семейства, заданная на внутренней кромке оболочки х1 = г0. Из равенств (41) вытекает, что при любых углах армирования ф (см. (40)) общий расход арматуры в пологих оболочках фиксированных размеров при заданных значе-
ниях ш0к будет один и тот же. В расчетах принимается ш01 = ш02 = 0,35 . В этом случае на внутренней кромке (х1 = г0) суммарная плотность армирования равна 0,7 (ш01 +ю02 = 0,7), что на практике соответствует предельно допустимому значению этой величины. Согласно (41) с удалением от внутренней кромки (х1 > г0) суммарная интенсивность армирования
(ш1 (х1) + ю2 (х1) ) уменьшается и становится меньше предельно допустимого значения 0,7.
Внешние кромки пологих оболочек предполагаем жестко закрепленными, т.е. в соотношениях (24) и (28) согласно (30) нужно принять w„= 0, ищ = 0, й„. = 0 (х1 = г1)
и учесть, что в(3) = 0, . = 1, 2. В отверстия же оболочек вставлены абсолютно жесткие
плоские шайбы, которые жестко прикреплены к пологим оболочкам по их внутренним кромкам, т.е. имеют место условия сопряжения
и 1х= ,= 0, У131х=г = 0, Wo (г) = w| г > г0, (42)
|х1 = г0 |х1 = г0 х ' |х1 =г0
где w0 - осевое перемещение шайбы как жесткого целого в направлении х3.
При рассматриваемых типах нагружения (см. (38) и (39)), закрепления (см. (42)) и структурах армирования (см. (40) и (41)) в пологих оболочках реализуется случай осе-симметричного изгиба, т. е. решения соответствующих начально-краевых задач не зависят от полярного угла х2. Кроме того, поперечный сдвиг и тангенциальное перемещение
в окружном направлении отсутствуют, т.е. у23 = 0 и и2 = 0 .
Чтобы в рассматриваемом случае замкнуть систему разрешающих уравнений (см. п. 1), необходимо использовать дополнительное уравнение движения жесткой шайбы, которое с учетом (38) имеет вид
М0^ = Р (г) + 2^0 ^ (г, х )|х =Г0 , г > гo, (43)
где М0 - масса шайбы; _Р13 - поперечная сила, определенная по формулам (20) на внутренней кромке оболочки ( х1 = г0).
В начальный момент времени г = г0 = 0 конструкции находятся в состоянии покоя (см. (26) и (29) при ит = 0, У0г. = 0, и0у = 0, У0у = 0, I = 1, 3 , у = 1, 2 ), а внешняя нагрузка
отсутствует (см. (38)), поэтому справедливы соотношения (34). При замене в уравнении (43) второй производной по времени от w0 ее конечно-разностным аналогом исследуемая задача с учетом условий (42) может быть проинтегрирована численно по явной схеме «крест» (см. п. 2).
Предполагаем, что нелинейно-упругая зависимость напряжения а от деформации в в каждом компоненте композиции оболочки при одноосном растяжении-сжатии характеризуется идеализированной диаграммой с линейным упрочнением [29]:
а = <
екв, Ы<в<к) =ак)/ ек,
(44)
^п(в)а^) + Ек) (в-^п(в)в^)), |в|>в(к), 0 < к< К,
где Ек, Ек(к) - модули Юнга и линейного упрочнения материала к-й фазы композиции; а( к), в^к) - напряжение и соответствующая ему деформация, при превышении которых
материал ведет себя нелинейно. Физико-механические характеристики материалов фаз композиции пологих оболочек приведены в таблице, где V - коэффициент Пуассона. Характеристики материала связующей матрицы условно соответствуют эпоксидной смоле, отвержденной ароматическим аммиаком [15], а характеристики волокон - арматуре из Ь-стекла [14]. В последнем столбце таблицы указаны значения скорости звука а в материалах фаз композиции, вычисленные по формулам, аналогичным (36).
Физико-механические характеристики материалов фаз композиции [14, 15] Physical and mechanical characteristics of the materials of the composition phase [14, 15]
Материал р, кг/м3 V cs, МПа E, ГПа Es, ГПа a, м/с
Связующая матрица 1210 0,4 20 2,8 1,114 1521,2
Армирующие волокна 4300 0,25 1500 51,0 10,851 3443,9
Для шайб условно примем р = 2710 кг / м3, что соответствует алюминиевым сплавам, а их толщины зададим вдвое большими толщин 2к соответствующих пологих оболочек. По этим данным можно рассчитать массы шайб М 0 (см. (43)).
При проведении расчетов вдоль полярного радиуса х1 (см. рис. 3) вводилась регулярная сетка с шагом Ах1 = ( - г0) /100 = 5 мм, а шаг по времени т = 1 мкс. На основании этих данных для первого неравенства (35) получаем отношение Ах1 / т = 5 км /с. Кроме
того, в случае пологих оболочек рассматриваемых толщин ( 2к = 1 и 5 см) для второго неравенства (35) имеем 2к/т = 10 и 50 км/с. Эти отношения превышают значения а, указанные в таблице для материалов фаз композиции, поэтому необходимые условия устойчивости (35) используемой схемы «крест» выполняются со значительным запасом.
На рис. 4 и 5 изображены зависимости н0 (t), т.е. осцилляции жестких шайб, рассчитанные для пластин и пологих оболочек разной толщины, при разных структурах армирования и при разных уровнях динамического нагружения. Кривые 1, 2 и 3 на этих рисунках определены при углах армирования ф = 0, л/4 и л/2 соответственно (см. (40) и рис. 3, где изображен случай ф = л / 4 - спиральное ортогональное армирование). На рис. 4 приведены результаты расчетов для относительно тонких ( 2к = 1 см, 2к / ( - г0) = 1 / 50 ), а на рис. 5 - для относительно толстых ( 2к = 5 см, 2к / (г1 - г0) = 1/10) композитных конструкций. Кривые на рис. 4, а и рис. 5, а характеризуют динамическое поведение шайб в случаях армированных пластин (/ = 0 и Я ^ю), а на рис. 4, б и 4, в, а также на рис. 5, б и 5, в -в случаях пологих оболочек со стрелой подъема / = 0,2 м (см. (37)).
Кривые на рис. 4, а и рис. 5, а определены при значениях ртах = 0,3 МПа и Ртах = 1,6 МПа соответственно, т.е. согласно (38) конструкции нагружаются избыточным давлением р () со стороны нижних лицевых поверхностей (х3 =-к). Так как эти случаи соответствуют пластинам, изменение знака нагрузки р ^) при сохранении ее интенсивности приводит к тому, что кривые на этих рисунках зеркально отражаются относительно горизонтальной оси н0 = 0, поэтому такие кривые не изображены. Сравнение кривых с одинаковыми номерами на рис. 4, а и рис. 5, а показывает, что в случае относи-
тельно тонких гибких пластин (см. рис. 4, а) шайбы совершают более сложные осцилляции, чем в случае относительно толстых конструкций (рис. 5, а).
vv0, см 10
w0, см 1,01 1
0,5 0,0 -0,5 -1,0
-1,5 -2,0
10 15 20 25 30 35 U мс
10 15 20 25 30 35 f, мс
Рис. 4. Осцилляции шайбы в случае относительно тонких композитных конструкций при трех структурах армирования: гибкой пластины (а); пологой оболочки, нагруженной снизу (б); пологой оболочки, нагруженной сверху (в) Fig. 4. Oscillations of the puck in the case of relatively thin composite structures with three reinforcement structures: a flexible plate (a); a shallow shell loaded from the bottom (b); a shallow shell loaded from the top (c)
Кривые на рис. 4, б и рис. 5, б рассчитаны при pmax = 0,5 МПа и pmax = 3,4 МПа соответственно, т.е. при нагружении пологих оболочек избыточным давлением со стороны нижних (вогнутых) лицевых поверхностей (х3 = -И ). Кривые же на рис. 4, в и рис. 5, в определены при pmax = -0,5 МПа и pmax = -3,4 МПа, т.е. при нагружении пологих оболочек избыточным давлением со стороны верхних (выпуклых) лицевых поверхностей ( х3 = И).
Сопоставление кривых с одинаковыми номерами на рис. 4, б и рис. 4, в, а также на рис. 5, б и рис. 5, в демонстрирует, что в отличие от гибких пластин (см. рис. 4, а и рис. 5, а) динамическое поведение гибких армированных пологих оболочек, изготовленных из нелинейно-упругих материалов фаз композиции, существенно зависит от того, к какой лицевой поверхности (выпуклой или вогнутой) прикладывается избыточное давление p (t) .
Сравнение же кривых с разными номерами на каждом рис. 4 и 5 свидетельствует о том, что структуры армирования пластин и пологих оболочек, а точнее изменения углов армирования ф при фиксированном распределении плотностей армирования шк (см. (40)
и (41)) существенно влияют на динамическое поведение этих тонкостенных композитных конструкций.
Поведение кривых на рис. 4, а и рис. 5, а показывает, что в случае гибких пластин max |w0 (t|| достигается при первом локальном экстремуме зависимости w0 (t| . Однако, например, поведение пунктирной кривой 2 на рис. 4, б свидетельствует о том, что в случае пологих оболочек max |w0 (t | может достигаться и не при первом локальном экстре-
муме зависимости w0 (t) , т.е. динамическое поведение гибких композитных пологих оболочек является более сложным, чем гибких армированных пластин, изготовленных из нелинейно-упругих материалов компонентов композиции.
ч>п, см
3
). » . ( .2 |1
4 i ; J, s Î i {\ Î с i ^. .
(iiÉMililiM
10 15
(
25 30 35 t, mc
35 t, mc
Рис. 5. Осцилляции шайбы в случае относительно толстых композитных конструкций при трех структурах армирования: гибкой пластины (а); пологой оболочки, нагруженной снизу (б); пологой оболочки, нагруженной сверху (в) Fig. 5. Oscillations of the puck in the case of relatively thick composite structures with three reinforcement structures: a flexible plate (a); a shallow shell loaded from the bottom (b); a shallow shell loaded from the top (c)
Во всех приведенных выше примерах максимальные значения интенсивности деформаций в^) в компонентах композиций тонкостенных конструкций имели порядок 5 %
и менее, поэтому все представленные решения, полученные на базе определяющих соотношений (8), являются вполне корректными с механической точки зрения. При этом наибольшие по модулю деформации возникают в связующем материале при окружном армировании (ф = л/2) пластин и пологих оболочек (см. штриховые кривые 3 на рис. 4 и рис. 5). Наилучшими же с точки зрения минимизации податливости и напряженно-деформированного состояния в компонентах композиции являются структуры с радиальным (меридиональным) армированием: ф = 0 (см. сплошные кривые 1 на рис. 4 и рис. 5).
Выполненные расчеты показали, что при нелинейно-упругом поведении материалов фаз композиции для получения значений в^) (к = 0, 1, 2) с точностью до трех значащих
цифр на каждом шаге интегрирования по времени требуется делать до пяти итераций по методу переменных параметров упругости (см. п. 1 и 2).
Заключение
Проведенное исследование динамического поведения гибких композитных пологих оболочек и пластин, изготовленных из нелинейно-упругих материалов фаз композиции, продемонстрировало, что изменение направлений армирования существенно влияет как на
податливость таких тонкостенных конструкций, так и на напряженно-деформированное состояние в материалах фаз их композиций, особенно в материале связующей матрицы. Следовательно, на множестве структур армирования целесообразно осуществлять поиск оптимальных или рациональных проектов армирования. В частности, для гибких пластин и пологих сферических оболочек, имеющих в плане кольцевую форму, при наличии внутренней жесткой шайбы наилучшим с точки зрения минимума податливости тонкостенной конструкции является армирование в радиальном (меридиональном) направлении. В силу геометрической и физической нелинейности рассмотренной задачи динамический отклик армированных пологих оболочек существенно зависит от того, к каким лицевым поверхностям (выпуклым или вогнутым) прикладывается избыточное давление взрывного типа, так как при этом значительно различаются не только амплитуды, но и частоты колебаний.
Работа выполнена при финансовой поддержке РФФИ (грант № 17-01-00156-а). Библиографический список
1. Review of advanced composite structures for naval ships and submarines / A.P. Mouritz, E. Gellert, P. Burchill, K. Challis // Compos. Struct. - 2001. - Vol. 53. - No. 1. - P. 21-42.
2. Bannister M. Challenger for composites into the next millennium - a reinforcement perspective // Composites. - 2001. - Part A 32. - P. 901-910.
3. Gibson R.F. Principles of composite material mechanics / 3rd ed. - Boca Raton: CRC Press, Taylor & Francis Group, 2012. - 686 р.
4. Амбарцумян С.А. Общая теория анизотропных оболочек. - М.: Наука, 1974. - 446 с.
5. Андреев А. Упругость и термоупругость слоистых композитных оболочек. Математическая модель и некоторые аспекты численного анализа. - Saarbrucken (Deutschland): Palmarium Academic Publishing, 2013. - 93 c.
6. Cao Z., Cardew-Hall M. Interference-fit riveting technique in fiber composite laminates // Aerospace Sci. and Technol. - 2006. - Vol. 10. - P. 327-330.
7. Muc A., Ulatowska A. Desing of plates with curved fiber format // Compos. Struct. - 2010. -Vol. 92. - No. 7. - P. 1728-1733.
8. Muc A., Muc-Wierzgon M. An evolution strategy in structural optimization problems for plates and shells // Compos. Struct. - 2012. - Vol. 94. - No. 4. - P. 1461-1470.
9. Абросимов Н.А., Баженов В.Г. Нелинейные задачи динамики композитных конструкций. -Н. Новгород: Изд-во Нижегород. гос. ун-та, 2002. - 400 с.
10. Баженов В.А., Кривенко О.П., Соловей Н.А. Нелинейное деформирование и устойчивость упругих оболочек неоднородной структуры: Модели, методы, алгоритмы, малоизученные и новые задачи. - М.: ЛИБРОКОМ, 2012. - 336 с.
11. Белькаид К., Тати А., Бумараф Р. Простой конечный элемент с пятью степенями свободы в узле, основанный на теории сдвигового деформирования третьего порядка // Механика композитных материалов. - 2016. - Т. 52, № 2. - С. 367-384.
12. Macko W., Kowalewski Z.L. Mechanical properties of A359/SiCp metal matrix composites at wide range of strain rates // Appl. Mech. Mater. - 2011. - Vol. 82. - P. 166-171.
13. Gill S.K., Gupta M., Satsangi P. Prediction of cutting forces in machining of unidirectional glass-fiber-reinforced plastic composites // Frontiers of Mechanical Eng. - 2013. - Vol. 8. - No. 2. - P. 187-200.
14. Композиционные материалы: справ. / под ред. Д.М. Карпиноса. - Киев: Наук. думка, 1985. - 592 с.
15. Справочник по композитным материалам: в 2 кн. Кн. 1 / под ред. Дж. Любина; пер. с англ. А.Б. Геллера, М.М. Гельмонта; под ред. Б.Э. Геллера. - М.: Машиностроение, 1988. - 448 с.
16. Янковский А.П. Упругопластическое деформирование изгибаемых армированных пластин при ослабленном сопротивлении поперечному сдвигу // Прикладная математика и механика. -2013. - Т. 77. - Вып. 6. - С. 853-876.
17. Reissner E. The effect of transverse shear deformations on the bending of elastic plate // J. Appl. Mech. - 1945. - Vol. 12. - No. 2. - P. 69-77.
18. Шкутин Л.И. Нелинейные деформации и катастрофы тонких тел. - Новосибирск: Изд-во СО РАН, 2014. - 139 с.
19. Малмейстер А. К., Тамуж В. П., Тетерс Г. А. Сопротивление жестких полимерных материалов. - Рига: Зинатне, 1972. - 500 с.
20. Reddy J.N. Energy and Variational Methods in Applied Mechanics. - N.Y.: John Wiley, 1984. -604 р.
21. Янковский А.П. Практическая устойчивость схемы «крест» при численном интегрировании уравнений динамики для гибких тонкостенных элементов конструкций, подчиняющихся гипотезам теории Тимошенко // Мат. методи та фiз.-мех. поля. - 2015. - Т. 58, № 1. - С. 65-83.
22. Рихтмайер Р., Мортон К. Разностные методы решения краевых задач. - М: Мир, 1972. -418 с.
23. Самарский А.А. Теория разностных схем. - М.: Наука, 1989. - 616 с.
24. Houlston R., DesRochers C.G. Nonlinear structural response of ship panels subjected to air blast loading // Computers & Structures. - 1987. - Vol. 26. - No. 1/2. - P. 1-15.
25. Zeinkiewicz O.C., Taylor R.L. The finite element method. - Oxford: Butterworth-Heinemann, 2000. - 707 p.
26. Ильюшин А.А. Труды. Т. 3. Теория термовязкоупругости / сост. Е.А. Ильюшина, В.Г. Тунгускова. - М.: Физматлит, 2007. - 288 с.
27. Янковский А.П. Моделирование механического поведения композитов с пространственной структурой армирования из нелинейно-наследственных материалов // Конструкции из композиционных материалов. - 2012. - № 2. - С. 12-25.
28. Трещёв А.А. Изотропные пластины и оболочки, выполненные из материалов, чувствительных к виду напряженного состояния. - М.; Тула: РААСН; Изд-во ТулГУ, 2013. - 249 с.
29. Малинин Н.Н. Прикладная теория пластичности и ползучести. - М.: Машиностроение, 1968. - 400 с.
30. Бидерман В.Л. Механика тонкостенных конструкций. Статика. - М.: Машиностроение, 1977. - 488 с.
31. Зубчанинов В.Г. Основы теории упругости и пластичности. - М.: Высш. шк., 1990. - 368 с. 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.
2. Bannister, M. Challenger for composites into the next millennium - a reinforcement perspective. Composites, 2001, part A 32, pp. 901-910.
3. Gibson R.F. Principles of composite material mechanics / 3rd ed. Boca Raton, CRC Press, Taylor & Francis Group, 2012, 686 p.
4. Ambarcumian S.A. Obshchaya teoriya anizotropnykh obolochek [The general theory of anisotropic shells]. Moskva, Nauka, 1974, 446 p.
5. Andreev A.N. Uprugost' i termouprugost' sloistykh kompozitnykh obolochek. Matematicheskaia model i nekotorye aspekty chislennogo analiza [Elasticity and thermo-elasticity layered composite shells. Mathematical model and some aspects of the numerical analysis.]. Saarbrucken (Deutschland), Palmarium Academic Publishing, 2013, 93 p.
6. Cao, Z., Cardew-Hall, M. Interference-fit riveting technique in fiber composite laminates. Aerospace Sci. and Technol., 2006, vol. 10, pp. 327-330.
7. Muc, A., Ulatowska, A. Desing of plates with curved fiber format. Compos. Struct., 2010, vol. 92, no. 7, pp. 1728-1733.
8. Muc, A., Muc-Wierzgon, M. An evolution strategy in structural optimization problems for plates and shells, Compos. Struct., 2012, vol. 94, no. 4, pp. 1461-1470.
9. Abrosimov N.A., Bazhenov V.G. Nelinejnye zadachi dinamiki kompositnykh konstrukcij [Nonlinear problems of dynamics composites designs]. Nizhniy Novgorod, Publ Nizhniy Novgorod State University, 2002, 400 p.
10. Bazhenov V.A., Krivenko O.P., Solovei N.A. Nelineinoe deformirovanie i ustoichivost' uprugikh obolochek neodnorodnoi struktury: Modeli, metody, algoritmy, maloizuchennye i novye zadachi [Nonlinear deformation and stability of elastic shells of non-uniform structure: Models, methods, algorithms, the insufficiently studied and new problems]. Moscow, Knizhnyi dom "LIBROKOM", 2012, 336 p.
11. Belkaid, K., Tati, A., Boumaraf, R. A simple finite element with five degrees of freedom based on Red-dy's third-order shear deformation theory, Mechanics of Composite Materials, 2016. vol. 52, no. 2, pp. 367-384.
12. Macko, W., Kowalewski, Z.L. Mechanical properties of A359/SiCp metal matrix composites at wide range of strain rates, Appl. Mech. Mater., 2011, vol. 82, pp. 166-171.
13. Gill, S.K., Gupta, M., Satsangi, P. Prediction of cutting forces in machining of unidirectional glass-fiber-reinforced plastic composites, Frontiers of Mechanical Eng., 2013, vol. 8, no. 2, pp. 187-200.
14. Karpinos D.M. Kompozitsionnye materialy. Spravochnik [Composite materials. Reference Book]. Naukova dumka, 1985, 592 p.
15. Lubin G. Handbook of composites. New York: Van Nostrand Reinhold Company Inc., 1982, 442 p.
16. Yankovskii, A.P. Elastoplastic deformation of flexible reinforced plates with a reduced shear strength. J. Appl. Math. Mech., 2013, vol. 77, no. 6, pp. 613-628.
17. Reissner, E. The effect of transverse shear deformations on the bending of elastic plate, J. Appl. Mech., 1945, vol. 12, no. 2, pp. 69-77.
18. Shkutin L.I. Nelineinye deformacii i katastrofy tonkikh tel [Nonlinear deformations and accidents of thin bodies]. Novosibirsk, SBRASPubl, 2014, 139 p.
19. Malmeister A.K., Tamuzh V.P., Teters G.A. Soprotivlenie zhyestkikh polimernykh materialov [Resistance rigid polymeric materials]. Riga, Zinatne, 1972, 500 p.
20. Reddy J.N. Energy and Variational Methods in Applied Mechanics. New York, John Wiley, 1984, 608 p.
21. Yankovskii, A.P. Practical stability of the "cross" scheme in the numerical integration of dynamic equations for flexible thin-walled structural elements under the hypotheses of the Timoshenko theory. Journal of Mathematical Sciences, 2016, vol. 207, no. 1, pp 65-83.
22. Richtmyer R.D., Morton K.W. Difference methods for initial-value problems. New York, John Wiley & Sons, 1967, 420 p.
23. Samarakii A.A. Teoriya raznostnykh skhem [The theory of finite difference schemes]. Moscow, Nauka, 1989, 616 p.
24. 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.
25. Zeinkiewicz O.C., Taylor R.L. The finite element method. Oxford, Butterworth-Heinemann, 2000, 707 p.
26. Il'yushina E.A., Tunguskova V.G. Trudy. Tom 3. Teoriya termovyazkouprugosti [Works. Vol. 3. The theory thermo-visco-elasticity]. Moskva, Fizmatlit, 2007, 288 p.
27. Yankovskii, A.P. Modelirovanie mekhanicheskogo povedenia kompozitov s prostranstvennoi strukturoi armirovania iz nelineino-nasledstvennykh materialov [Modelling of mechanical behaviour of the is space-reinforced composites from is nonlinear-hereditary materials]. Konstrukcii iz kompozicionnykh materialov - Constructions from composite materials, 2012, no. 2, pp. 12-25.
28. Treshchev A.A. Izotropnye plastiny i obolochki, vypolnennye iz materialov, chuvstvitel'nykh k vidu napriazhennogo sostoyania [Isotropic plates and shells made from materials, sensitive to a kind of an stresses condition]. Moscow - Tula, RAASNPubl.; Publ. TulSU, 2013, 249 p.
29. Malinin N.N. Prikladnaya teoriya plastichnosti i polzuchesti [The applied theory of plasticity and creep]. Moscow, Mashinostroenie, 1968, 400 p.
30. Biderman V.L. Mekhanika tonkostennykh konstruktsiy. Statika [Mechanics of thin-walled designs. Statics]. Moscow, Mashinostroenie, 1977, 488 p.
31. Zubchaninov V.G. Osnovy teorii uprugosti i plastichnosti [Bases of the theory of elasticity and plasticity]. Moscow, Higher school, 358 p.
