The Refined Model of Viscoelastic-Plastic Deformation of Reinforced Cylindrical Shells

Abstract


The paper formulates the initial-boundary-value problem of the viscoelastic-plastic bending behavior of cylindrical circular shells cross-reinforced along equidistant surfaces. The instant elastoplastic deformation of the shell composition components is described by the governing equations of the theory of plastic flow with isotropic hardening. The viscoelastic deformation of these materials is described by the defining relations of the Maxwell - Boltzmann model of body. The geometric nonlinearity of the problem is taken into account in the Karman approximation. The used system of two-dimensional resolving equations and the corresponding initial and boundary conditions make it possible to determine displacements and stress-strain state (including residual one) in materials of the composition of flexible cylindrical shells with varying degrees of accuracy. In this case, the weak resistance of the considered composite structures to transverse shears is taken into account. In the first approximation, the equations are used, the initial and boundary conditions correspond to the relations of the widely used non-classical Reddy theory. A numerical solution of the initial-boundary-value problem posed is constructed using an explicit step-by-step "cross" scheme. The elastoplastic and viscoelastic-plastic dynamic deformation of a relatively thin long circular cylindrical shell is investigated. The structure is rationally reinforced in the circumferential direction and is loaded with an internal pressure of an explosive type. It has been demonstrated that under intense short-term loading even of a relatively thin cylindrical reinforced shell by internal pressure, the traditional Reddy theory does not guarantee that the maximum residual deflection and the intensity of residual deformations of the components of the composition are accurate to within 10% compared to calculations performed by the refined theory. The difference in the results of the corresponding calculations increases with an increase in the relative thickness of the composite shell. It was found that after plastic deformation of a long reinforced cylindrical shell in its residual state, not only appear zones of edge effects, but also a local zone of an intense deformation located in the vicinity of the central section of the shell. The length of the local central zone is comparable with the length of the zones of edge effects. It is shown that the amplitude of the transverse vibrations of the reinforced shell in the vicinity of the initial moment of time significantly (by an order of magnitude) exceeds the value of the maximum modulus of the residual deflection. Therefore, the calculations performed in the framework of the theory of elastoplastic deformation of composition materials do not allow a very approximate determination of the magnitude of the residual displacements and the magnitude of the residual deformed state of the components of the composition of the cylindrical shell.

Full Text

Введение Цилиндрические оболочки из композиционных материалов (КМ) широко используются в инженерных приложениях в качестве тонкостенных силовых элементов [1-5]. Поэтому актуальной является проблема математического описания механического поведения таких оболочек, особенно при воздействии высокоинтенсивных динамических нагрузок [5]. В рамках классической теории деформирования слоистые и анизотропные оболочки изучались в [6-8], но при этом игнорируется экспериментально установленный факт слабого сопротивления КМ-оболочек поперечным сдвигам, которые традиционно учитываются в рамках теории Рейсснера - Миндлина [5, 9-14] или теории Редди [7, 15-18], реже используются теории более высоких порядков точности [19-23]. При этом в подавляющем большинстве публикаций рассматривается линейно-упругое или линейно-вязкоупругое деформирование тонкостенных конструкций из КМ [5-7, 9-11, 13-17, 19-22]. Однако известно, что компоненты композиции могут деформироваться и пластически [4, 8, 18, 23-26]. Проблема моделирования неупругого поведения КМ-конструкций при высокоинтенсивных нагрузках находится сейчас на стадии становления. Так, в [8] исследовалось предельное состояние армированных цилиндрических оболочек на базе модели жесткопластического деформирования компонентов композиции; в [12] изучалось динамическое упругопластическое поведение слоистых тонкостенных конструкций с изотропными слоями; в [23] моделировалось упругопластическое деформирование волокнистых цилиндрических оболочек. Ранее в [15] было продемонстрировано, что при линейно-упругом деформировании компонентов композиции оболочек не требуется уточнять результаты расчетов, выполненные по теории Редди. Но в [23] показано, что для адекватного описания динамического неупругого поведения армированных цилиндрических оболочек необходимо использовать более точные, чем теория Редди, модели деформирования таких тонкостенных элементов КМ-конструкций. Механическое поведение материалов композиции в [12, 23] описывалось определяющими соотношениями теории Прандтля - Рейсса - Хилла (ТПРХ) [27], в которых не учитываются демпфирующие свойства реальных материалов при осциллирующем деформировании [28]. Следовательно, модели, использованные в [12, 23], не позволяют рассчитывать остаточные перемещения и остаточное напряженно-деформированное состояние (НДС) в компонентах композиции при интенсивных динамических воздействиях на КМ-оболочки. Демпфирующие свойства материалов композиции волокнистых пластин учтены в [18], где была разработана структурная модель вязкоупругопластического поведения перекрестно армированного по эквидистантным поверхностям тонкостенного КМ-элемента. Для численного решения динамических задач механики тонкостенных элементов КМ-конструкций используют явные схемы, например типа «крест» [12, 18, 23], или неявные методы, как правило, Ньюмарка [29-32]. В связи с вышеизложенным настоящее исследование посвящено математическому моделированию вязкоупругопластического поведения гибких цилиндрических волокнистых оболочек в рамках уточненной теории их изгиба при использовании явной численной схемы типа «крест». 1. Постановка задачи Рассмотрим тонкую замкнутую круговую цилиндрическую оболочку радиусом R, длиной L и толщиной (рис. 1), с которой связана цилиндрическая система координат , где - осевая координата ( ), - угловая окружная координата ( ), - радиальная координата ( ). Оболочка усилена N семействами волокон с плотностями и углами ( ) армирования, отсчитываемыми от направления . Волокна уложены по поверхностям, эквидистантным срединной поверхности конструкции ( ), и структура армирования однородна в радиальном направлении. Как и в [23], ограничимся исследованием практически важного частного случая, когда на лицевых поверхностях оболочки можно пренебречь действием касательных внешних распределенных сил. Согласно этому предположению для учета слабого сопротивления рассматриваемой КМ-оболочки поперечным сдвигам аппроксимируем усредненные деформации композиции и перемещения точек конструкции следующими соотношениями [23] (геометрическая нелинейность задачи учитывается в приближении Кармана): (1) (2) (3) (4) где (5) Рис. 1. Цилиндрические КМ-оболочки: а - жестко закрепленная только на левой кромке; б - жестко закрепленная на обеих кромках Fig. 1. Cylindrical shells clamped only along the left edge (a) and clamped on both edges (b) - коэффициенты частичных сумм разложений деформаций поперечных сдвигов по степеням переменной z; K - целое число, равное количеству слагаемых, удерживаемых в этих разложениях; z - введенная для удобства новая радиальная координата, отсчитываемая от срединной поверхности; w - прогиб точек этой поверхности ( ); , - продольное и окружное перемещения точек срединной поверхности; - начальный момент времени t; r - вектор места точки конструкции; x - вектор места точки срединной поверхности; - область, занимаемая этой поверхностью в координатах ; суммирование по индексу s проводится от 0 до K, как это указано в (1). Равенство (4) - традиционная для тонкостенных элементов конструкций кинематическая гипотеза, согласно которой изменяемостью поперечного перемещения в радиальном направлении можно пренебречь [5-18, 21-23, 29]. Таким образом, в равенствах (1)-(4) при учете выражений (5) неизвестны функции w, и ( ), зависящие от времени t и только двух пространственных переменных и . При из (1)-(5) вытекают кинематические соотношения теории Редди [15-17, 23]. Поскольку исследуется динамическое поведение КМ-оболочки как гибкой тонкостенной механической системы, то напряжение с приемлемой для практических приложений точностью можно аппроксимировать линейно по поперечной координате z [11]: (6) где - нормальные напряжения на внешней (+) и внутренней (-) лицевых поверхностях, известные из соответствующих силовых граничных условий (так, на рис. 1, а изображен случай, когда и ). Двумерные уравнения движения тонкой цилиндрической оболочки, соответствующие аппроксимации перемещений (3) и (4), при учете (5) и (6) имеют вид [23]: (7) где (8) - усредненные напряжения в композиции; , - объемная плотность материала связующего и арматуры k-го семейства; - компоненты объемной нагрузки, действующей на КМ-оболочку; точка - производная по времени t. Согласно четвертому равенству (8), в рамках гипотез, используемых для тонкостенных элементов конструкций, некоторым из введенных силовых факторов можно дать механическую трактовку: , - приближенные выражения для мембранных усилий и изгибающих и крутящих механических моментов; ( ) - приближенно представляют собой поперечные силы; остальные силовые факторы , и , определенные в (8), - моменты высших порядков от напряжений, вычисленные с весами , и соответственно. Материал k-го компонента композиции предполагается однородным и изотропным, а его определяющие уравнения вязкоупругопластического деформирования при учете (6) можно записать в матричной форме [18]: (9) Здесь и далее: (10) (11) , , , - 5×5-матрицы, , , , - пятикомпонентные вектор-столбцы, ненулевые элементы которых определяются так: (12) , - компоненты тензора напряжений и деформаций; , - модуль Юнга и коэффициент Пуассона; - коэффициент линейной вязкости при сдвиге; - коэффициент объемной вязкости; - значение интенсивности касательных напряжений , при котором впервые достигается пластическое состояние; - максимальное значение , достигнутое за всю предыдущую историю деформирования элемента материала; - символ Кронекера; индекс T - операция транспонирования; индекс k - номер компонента композиции ( - связующее, - волокна k-го семейства). Соотношения (11) задают соответствия между пятью компонентами ( ) некоторого вектор-столбца и компонентами симметричного тензора второго ранга ( , ). Соотношения (9) при учете (10) и (12) получены в предположении, что малые деформации могут быть разложены на вязкоупругие и пластические (несжимаемые) составляющие. При этом мгновенное пластическое деформирование определяется уравнениями теории течения с изотропным упрочнением [27], а вязкоупругое деформирование - уравнениями модели тела Максвелла - Больцмана [18, 28]: (13) где (14) Если вязкие свойства в (13) не учитываются ( и ), то эти соотношения при учете (14) редуцируются в закон Гука, а равенства (9) - в определяющие уравнения ТПРХ [23, 27], так как согласно (10) и (12) имеем и . Используя соотношения (9) и повторяя рассуждения из [18], для текущего момента времени t получим следующее определяющее уравнение, записанное в матричной форме и описывающее вязкоупругопластическое деформирование волокнистой композиции рассматриваемой цилиндрической оболочки: (15) где (16) , - пятикомпонентные вектор-столбцы скоростей усредненных напряжений и деформаций в композиции, по структуре аналогичные (11); B, , - 5×5-матрицы; , - матрицы, обратные 5×5-матрицам и H; p, f, g, , - пятикомпонентные векторы-столбцы. Элементы , матриц , и компоненты вектор-столбцов определяются по формулам (55) и (56) из [18] и зависят от механического состояния k-го материала композиции и направлений армирования, т.е. от углов (см. рис. 1, б). Для однозначного интегрирования исследуемой начально-краевой задачи необходимо задать начальные и граничные условия. Если на торцевой поверхности КМ-оболочки ( и/или ) заданы силовые граничные условия, то для силовых факторов, введенных в (8), получим [23]: (17) где (18) - заданные на торцевой поверхности внешние распределенные нагрузки, действующие в направлениях ( ); при и при . Если на торцевой поверхности оболочки заданы кинематические граничные условия, то для обобщенных кинематических переменных, введенных в (8), при учете (4) имеем равенства [23]: (19) где (20) - заданный на кромке прогиб; - заданные на торцевой поверхности оболочки тангенциальные перемещения в направлениях При следует использовать начальные условия (см. (4) и (8)) [23]: (21) где (22) , ( ) - известные в начальный момент времени перемещения и скорости точек КМ-оболочки. Кинематические переменные w, и ( ), использованные в (1)-(4), и обобщенные кинематические переменные ( , ), введенные в (8), связаны матричными равенствами [23] (23) где (24) - матрицы, обратные -матрицам , - -компонентные вектор-столбцы, элементы которых определяются по формулам (31) и (33) в [23] независимо от решения исследуемой задачи. Если в текущий момент времени t из каких-либо соображений известны значения функций w и ( ), то из уравнений (23) при учете (24) можно вычислить значения функций , ( , ), характеризующих усредненные деформации композиции (1), (2) и тангенциальные перемещения (3) точек цилиндрической КМ-оболочки. 2. Метод расчета Как отмечалось выше, для численного интегрирования рассматриваемой задачи используем метод шагов по времени [12, 18, 23, 29-32], т.е. будем определять неизвестные функции в дискретные моменты времени ( ), где - шаг по времени. При этом, согласно [18, 23], предполагаем, что в моменты времени уже известны значения следующих функций: (25) Используя выражения (8) при учете (25), в данный момент времени можно вычислить все силовые факторы , и , входящие в уравнения (7), и граничные условия (17). Как и в [18, 23], производные по времени t аппроксимируем центральными конечными разностями на трехточечном шаблоне { }. Это позволяет построить явную численную схему типа «крест». Замена вторых производных по времени в левых частях равенств (7) при учете обозначений, аналогичных (25), приводит к уравнениям (26) На основании (8) и (25) получаем, что правые части в (26) известны. Добавив к этим уравнениям необходимые граничные условия (17) и/или (19) при учете (18) и (20), можем определить по явной схеме значения функций и ( , ) в следующий момент времени . Далее по формулам (23) при учете (24) вычисляем функции и ( , ), используя которые по формулам (1) и (2) получаем усредненные деформации композиции в каждой точке конструкции. Последующее решение рассматриваемой задачи при учете соотношений (15) и (16) строится совершенно так же, как и в [18, 23]. Согласно структуре левых частей уравнений (26), для начала расчета по разработанной явной схеме требуется знать не только значения функций и , известные из начальных условий (21) при учете (22), но и значения функций , (см. (26) при ). Значения этих функций определим по формуле Тейлора при учете начальных условий (21) и уравнений движения (7), рассмотренных в момент времени [23]: В случае замкнутой круговой цилиндрической КМ-оболочки область Ω (см. (4)) является прямоугольной, поэтому после аппроксимации в равенствах (17) и (26) производных по пространственным переменным их конечно-разностными аналогами от уже известных (см. (8) и (25)) в данный момент времени сеточных функций окончательно получим явную численную схему «крест» [12, 18, 23]. Необходимые условия устойчивости схемы «крест» вытекают из критерия Куранта - Фридрихса - Леви и для тонкостенных элементов конструкций приведены в [12]. В случае волокнистой оболочки эти условия выполняются с запасом, если выполняются для каждого компонента композиции в отдельности. 3. Обсуждение результатов расчетов В качестве конкретного примера рассмотрим динамическое вязкоупругопластическое и упругопластическое поведение тонкой длинной цилиндрической оболочки радиусом , длиной и толщиной ( ). Конструкция жестко закреплена на обеих кромках: и , , , и (см. (19), (20) и рис. 1, б), и нагружается равномерным внутренним давлением по закону (см. (6)-(8), (25) и (26)) [29]: (27) где (28) - момент времени, в который достигает наибольшего значения ( ); - момент времени, при превышении которого можно пренебречь по сравнению с (согласно (28), ). Примем , и . Эти значения условно соответствуют давлению, вызванному воздушной взрывной волной [29]. Внешние массовые нагрузки не учитываем: , (см. (8)). В начальный момент времени оболочка покоится: и , (см. (21) и (22)). Конструкция изготовлена из эпоксисвязующего и усилена стеклянными волокнами марки S-994. Мгновенное упругопластическое деформирование материалов композиции при активном нагружении определяется билинейной диаграммой растяжения - сжатия где и - осевое напряжение и деформация; - модуль линейного упрочнения материала k-го компонента композиции; , - условный предел текучести и соответствующая ему деформация того же материала. Физико-механические характеристики компонентов композиции приведены в таблице, в которой принято , т.е. объемная вязкость материалов не учитывается (см. (12) и (13)). Физико-механические характеристики компонентов композиции [24, 25] Physical and mechanical characteristics of components of the composition [24, 25] Материал ν Эпоксисвязующее 1210 0,33 20 2,8 1,114 150 Стекловолокно S-994 2520 0,25 4500 86,8 6,230 600 Конструкция армирована в окружном направлении ( ; см. рис. 1, б) с постоянной плотностью . Для длинной тонкой цилиндрической оболочки такая структура армирования является рациональной [8, 23, 33]. При указанных условиях закрепления, нагружения и армирования рассматриваемой оболочки она деформируется не только осесимметрично (решения соответствующей начально-краевой задачи не зависит от окружной координаты ), но и симметрично относительно центрального сечения (т.е., например, , , и т.п.). На рис. 2 изображены поперечные колебания точек центрального сечения оболочки ( , ), полученные по ТПРХ (рис. 2, а, который полностью совпадает с рис. 3, б в [23]) и по теории вязкоупругопластического деформирования материалов композиции конструкции (рис. 2, б). Рис. 2. Осцилляции прогибов точек центрального сечения тонкой длинной цилиндрической КМ-оболочки, рассчитанные по ТПРХ (а) и по вязкоупругопластической теории (б) деформирования компонентов композиции Fig. 2. Oscillations of the deflections of the central section points of a thin long cylindrical CM-shell calculated by the Prandtl-Reuss-Hill theory (a) and by the viscoelastic theory (b) of the deformation of the composition components а б Рис. 3. Осцилляции максимального по модулю значения прогиба вязкоупругопластической КМ-оболочки в окрестности начального момента времени (а) и в окрестности момента времени (б), рассчитанные по разным теориям изгибного деформирования Fig. 3. Oscillations of the maximum modulus of deflection of the viscoelastic-plastic CM-shell in the vicinity of the initial instant of time (a) and in the vicinity of the instant of time (b), calculated according to different theories of bending deformation Расчеты проводились по уточненной теории изгиба КМ-оболочек при (см. (1)). Согласно поведению зависимости , приведенной на рис. 2, а, при расчете рассматриваемой конструкции по ТПРХ она неограниченно долго осциллирует, что не позволяет определить величину остаточного прогиба точек центрального сечения оболочки. Поведение же кривой на рис. 2, б свидетельствует о том, что при расчете по вязкоупругопластической модели колебания прогиба точек центрального сечения к моменту времени практически полностью прекращаются и остаточный прогиб имеет значение порядка 1 мм. На рис. 3 изображены зависимости от времени наибольшего по модулю значения прогиба КМ-оболочки ( ), рассчитанные в окрестности начального момента времени (рис. 3, а) и в окрестности момента времени (рис. 3, б). Кривая 1 получена по теории Редди ( ), а кривая 2 - по уточненной теории при (расчеты проводились по вязкоупругопластической модели деформирования компонентов композиции). На рис. 3, а кривые 1 и 2 визуально практически неразличимы. Поведение этих кривых на рис. 3, а показывает, что с течением времени поперечные осцилляции конструкции затухают и к моменту времени почти полностью прекращаются (см. рис. 3, б). Из рис. 3, б видно, что при кривые 1 и 2 заметно различаются, причем остаточный прогиб, определенный по теории Редди (кривая 1), примерно на 10,4 % больше аналогичной величины, рассчитанной по уточненной теории (кривая 2). В связи с этим представляет интерес рассмотрение эпюр остаточных прогибов, рассчитанных по разным теориям деформирования КМ-оболочки. На рис. 4 изображены зависимости w от x ( ), определенные при , когда поперечные колебания конструкции почти полностью уже прекратились. Кривые на рис. 4 рассчитаны при тех же условиях, что и на рис. 3. Поскольку решение исследуемой задачи обладает симметрией относительно центрального сечения оболочки ( ), то на рис. 4 приведены только правые участки соответствующих кривых (при ). Сравнение максимальных значений ординат точек на кривых рис. 4, которые достигаются при (т.е. в зоне краевого эффекта), показывает, что при значение , полученное по теории Редди (кривая 1), на 11,7 % больше аналогичной величины, определенной по уточненной теории (кривая 2). Рис. 4. Эпюры прогибов армированной вязкоупругопластической оболочки в момент времени , рассчитанные по разным теориям изгибного деформирования Fig. 4. Diagrams of deflections of a reinforced viscoelastic-plastic shell at a time instant calculated according to different theories of bending deformation На рис. 5 изображены рассчитанные при тех же условиях осцилляции максимальных значений интенсивности деформаций армирующих волокон ( , , ). Рис. 5. Осцилляции максимальных значений интенсивности деформаций арматуры вязкоупругопластической КМ-оболочки, рассчитанные по разным теориям изгибного деформирования Fig. 5. Oscillations of the maximum values of the strain intensity of the reinforcement of a viscoelastic-plastic CM-shell, calculated according to different theories of bending deformation Сравнение кривых на рис. 5 при свидетельствует о том, что расчет по теории Редди (кривая 1) примерно на 13 % завышает остаточную величину по сравнению с расчетом по уточненной теории (кривая 2). При кривые 1 и 2 монотонно возрастают, стремясь к соответствующим горизонтальным асимптотам и практически уже не осциллируя. Такая особенность поведения этих кривых объясняется ползучестью материала волокон, описываемой вязкими свойствами модели тела Максвелла - Больцмана (см. (13) и (14)). Эта ползучесть вызвана тем, что после пластического деформирования компонентов композиции в них возникают различные остаточные НДС, которые с течением времени стремятся выровняться. Заключение Разработанная уточненная математическая модель вязкоупругопластического поведения цилиндрических оболочек, армированных по эквидистантным поверхностям, позволяет более точно, чем в рамках широко используемой неклассической теории Редди, рассчитывать их затухающие колебания при кратковременных интенсивных механических воздействиях и определять остаточные перемещения и остаточное НДС в компонентах композиции после их пластического динамического деформирования. Выполненные расчеты показали, что даже для весьма тонких (с относительной толщиной порядка 1/100) и длинных цилиндрических КМ-оболочек традиционная неклассическая теория Редди, в первом приближении учитывающая искривление поперечной нормали конструкции, не гарантирует получения величины наибольшего остаточного прогиба и интенсивности остаточных деформаций компонентов композиции с точностью до 10 %. Как известно, для относительно толстых и/или коротких цилиндрических армированных оболочек (с относительной толщиной порядка 1/10) такое различие решений, полученных по уточненной теории и теории Редди, становится еще больше [23] и для остаточных деформаций может составлять несколько десятков процентов. Таким образом, для проведения адекватных динамических расчетов вязкоупругопластически деформируемых цилиндрических КМ-оболочек следует рекомендовать применение теорий повышенного порядка точности, в частности, предложенных в [23] и использованных настоящем исследовании. Продемонстрировано, что в тонкой длинной КМ-оболочке остаточный прогиб может сложным образом зависеть от продольной координаты. Так, помимо краевых эффектов в окрестности центрального сечения конструкции может возникнуть локальный эффект, протяженность которого сопоставима с протяженностью зон краевых эффектов (см. рис. 4). Такая сложная форма остаточного прогиба является следствием сложных неупругих волновых процессов, протекающих в динамически нагруженной тонкой КМ-конструкции. Высокопрочные армирующие волокна, имеющие большие пределы текучести и большие значения деформаций, им соответствующие, при деформировании запасают в большом количестве упругую энергию, поэтому амплитуды осцилляций прогиба точек цилиндрической КМ-оболочки в окрестности начального момента времени оказываются существенно (даже на порядок) больше максимальных значений величин остаточных прогибов (см. рис. 2, б и рис. 3). Аналогичное поведение наблюдается и для зависимостей от времени максимальных значений интенсивности деформаций компонентов композиции.

About the authors

A P Yankovskii

Khristianovich Institute of Theoretical and Applied Mechanics SB RAS

References

  1. Bannister M. Challenger for composites into the next millennium - a reinforcement perspective // Composites. - 2001. - Part A 32. - Pp. 901-910.
  2. 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.
  3. Gibson R.F. Principles of composite material mechanics. - 3rd ed. - Boca Raton: CRC Press, Taylor & Francis Group, 2012. - 686 p.
  4. 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.
  5. Прикладные задачи механики цилиндрических оболочек / Ю.С. Соломонов, В.П. Георгиевский, А.Я. Недбай, В.А. Андрюшин. - М.: Физматлит, 2014. - 408 с.
  6. Григоренко Я.М. Изотропные и анизотропные слоистые оболочки вращения переменной жесткости. - Киев: Наукова думка, 1973. - 228 с.
  7. Амбарцумян С.А. Общая теория анизотропных оболочек. - М.: Наука, 1974. - 446 с.
  8. Васильев В.В., Елпатьевский А.Н. Прочность цилиндрических оболочек из армированных материалов. - М.: Машиностроение, 1972. - 168 с.
  9. Reissner E. On transverse vibrations of thin shallow elastic shells // Quarterly of Applied Mathematics. - 1955. - Vol. 13, no. 2. - P. 169-176.
  10. Mindlin R.D., Bleich H.H. Response of an elastic cylindrical shells to a transverse step shock wave // Trans. ASME. Ser. E. J. Appl. Mech. - 1953. - Vol. 20, no. 2. - P. 189-195.
  11. Богданович А.Е. Нелинейные задачи динамики цилиндрических композитных оболочек. - Рига: Зинатне, 1987. - 295 с.
  12. Абросимов Н.А., Баженов В.Г. Нелинейные задачи динамики композитных конструкций. - Н. Новгород: Изд-во ННГУ, 2002. - 400 с.
  13. Баженов В.А., Кривенко О.П., Соловей Н.А. Нелинейное деформирование и устойчивость упругих оболочек неоднородной структуры: Модели, методы, алгоритмы, малоизученные и новые задачи. - М.: ЛИБРОКОМ, 2012. - 336 с.
  14. Шкутин Л.И. Нелинейные деформации и катастрофы тонких тел. - Новосибирск: Изд-во СО РАН, 2014. - 139 с.
  15. Андреев А.Н., Немировский Ю.В. Многослойные анизотропные оболочки и пластины. Изгиб, устойчивость и колебания. - Новосибирск: Наука, 2001. - 287 с.
  16. Reddy J.N. Mechanics of laminated composite plates and shells: Theory and analysis. - 2nd ed. - Boca Raton: CRC Press, 2004. - 831 p.
  17. Андреев А. Упругость и термоупругость слоистых композитных оболочек. Математическая модель и некоторые аспекты численного анализа. - Saarbrucken (Deutschland): Palmarium Academic Publishing, 2013. - 93 c.
  18. Янковский А.П. Моделирование вязкоупругопластического деформирования гибких армированных пластин с учетом слабого сопротивления поперечному сдвигу // Вычислительная механика сплошных сред. - 2019. - Т. 12, № 1. - С. 80-97. doi: 10.7242/1999-6691/2019.12.1.8.
  19. Whitney J., Sun C. A higher order theory for extensional motion of laminated composites // J. of Sound and Vibration. - 1973. - Vol. 30, no. 1. - P. 85-97.
  20. Lo K.H., Christensen R.M., Wu E.M. A higher-order theory of plate deformation. Part 2: Laminated plates // Trans. ASME, J. Appl. Mech. - 1977. - Vol. 44. - P. 669-676.
  21. Куликов Г.М. Термоупругость гибких многослойных анизотропных оболочек // Изв. РАН. МТТ. - 1994. - № 2. - С. 33-42.
  22. Пикуль В.В. Механика оболочек. - Владивосток: Дальнаука, 2009. - 536 с.
  23. Янковский А.П. Моделирование динамического поведения армированных цилиндрических оболочек при упругопластическом деформировании материалов компонентов композиции // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2018. - № 2. - С. 133-146.
  24. Композиционные материалы: справочник / под ред. Д.М. Карпиноса. - Киев: Наук. думка, 1985. - 592 с.
  25. Справочник по композитным материалам: в 2 кн. Кн. 1 / под ред. Дж. Любина; пер. с англ. А.Б. Геллера, М.М. Гельмонта; под ред. Б.Э. Геллера. - М.: Машиностроение, 1988. - 448 с.
  26. Maćko 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.
  27. Зубчанинов В.Г. Механика процессов пластических сред. - М.: Физматлит, 2010. - 352 с.
  28. Фрейденталь А., Гейрингер Х. Математические теории неупругой сплошной среды. - М.: Физматгиз, 1962. - 432 с.
  29. 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.
  30. Zeinkiewicz O.C., Taylor R.L. The finite element method. - Oxford: Butterworth-Heinemann, 2000. - 707 p.
  31. Librescu L., Oh S.-Y., Hohe J. Linear and non-linear dynamic response of sandwich panels to blast loading // Composites. - 2004. - Part B 35. - P. 673-683.
  32. Kazanci Z. Dynamic response of composite sandwich plates subjected to time-dependent pressure pulses // International Journal of Non-Linear Mechanics. - 2011. - Vol. 46. - P. 807-817.
  33. Muc A., Muc-Wierzgoń M. An evolution strategy in structural optimization problems for plates and shells // Compos. Struct. - 2012. - Vol. 94, no. 4. - P. 1461-1470.

Statistics

Views

Abstract - 350

PDF (Russian) - 217

Cited-By


PlumX


Copyright (c) 2020 Yankovskii A.P.

Creative Commons License
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies