Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2023. Т. 27, № 1. С. 119—141_
ISSN: 2310-7081 (online), 1991-8615 (print) d https://doi.org/10.14498/vsgtu1958
EDN: YRWNPW
УДК 539.4
Моделирование неизотермического упругопластического поведения армированных пологих оболочек в рамках уточненной теории изгиба
А. П. Янковский
Институт теоретической и прикладной механики им. С. А. Христиановича СО РАН,
Россия, 630090, Новосибирск, ул. Институтская, 4/1.
Аннотация
Сформулирована динамическая задача неизотермического и неупругого деформирования гибких пологих многонаправленно армированных оболочек в рамках уточненной теории их изгиба. По толщине конструкций температура аппроксимируется полиномом 7-го порядка. Геометрическая нелинейность задачи моделируется в приближении Кармана. Решение поставленной связанной нелинейной двумерной задачи получено с использованием явной численной схемы. Исследован термоупруго-пластический отклик стеклопластиковых и металлокомпозитных цилиндрических удлиненных панелей с ортогональной структурой армирования, нагруженных фронтально воздушной взрывной волной. Показано, что в отличие от аналогичных по структуре и характерным размерам армированных пластин пологие оболочки при интенсивном кратковременном нагружении необходимо рассчитывать с учетом возникновения в них температурных полей. При этом следует использовать уточненную теорию изгиба искривленных панелей, а не ее простейший вариант — неклассическую теорию Амбарцумяна. Приращение температуры в отдельных точках пологих оболочек из стеклопластика может достигать 14-34 °С, а в аналогичных металлокомпозитных панелях — 50-150 °С. Цилиндрические пологие оболочки более интенсивно деформируются при их нагружении воздушной взрывной волной со стороны выпуклой лицевой поверхности.
Ключевые слова: гибкие пологие оболочки, многонаправленное армирование, динамическое деформирование, связанная термоупругопла-стичность, уточненная теория изгиба, явная численная схема.
Механика деформируемого твердого тела Научная статья
© Коллектив авторов, 2023 © СамГТУ, 2023 (составление, дизайн, макет) 3 ©® Контент публикуется на условиях лицензии Creative Commons Attribution 4.0 International (https://creativecommons.org/licenses/by/4.0/deed.ru) Образец для цитирования
Янковский А. П. Моделирование неизотермического упругопластического поведения армированных пологих оболочек в рамках уточненной теории изгиба // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2023. Т. 27, № 1. С. 119-141. EDN: YRWNPW. DOI: 10.14498/vsgtu1958. Сведения об авторе
Андрей Петрович Янковский & https://orcid.org/0000- 0002- 2602- 8357
доктор физико-математических наук; ведущий научный сотрудник; лаб. физики быстро-
протекающих процессов; e-mail: lab4nemir@rambler.ru
Получение: 26 сентября 2022 г. / Исправление: 28 декабря 2022 г. / Принятие: 13 февраля 2023 г. / Публикация онлайн: 27 марта 2023 г.
Введение. Изделия из композитных материалов (КМ) находят широкое применение в современной инженерной практике [1-15] и часто испытывают интенсивное термосиловое нагружение [2,3,5,8-10,12,15], при котором компоненты композиции деформируются неупруго [5,8,12,16-18]. Следовательно, актуальной является проблема математического моделирования термоупру-гопластического деформирования КМ-пластин и КМ-оболочек, находящаяся на сегодняшний день в стадии становления [5,8,12,19-23]. Так, в работе [23] разработана модель неизотермического упругопластического деформирования материала, многонаправленно армированного непрерывными волокнами, и выполнены расчеты динамического изгиба КМ-пластин из таких материалов.
Для учета ослабленного сопротивления композитных изгибаемых пластин и оболочек поперечным сдвигам и моделирования волновых процессов в них используют простейшие неклассические теории Рейсснера [2,5,10,12,24], Ам-барцумяна [1,23], Редди [6,9] или уточненные теории [3,5,10,22].
В работе [23] показано, что температуру по толщине изгибаемой КМ-пластины при ее динамическом нагружении следует аппроксимировать полиномами 6-7 порядка. Для описания же изгибного поведения армированных пластин в [23] использована теория Амбарцумяна [1]. При интенсивном кратковременном механическом нагружении конструкций главным источником выделения тепла является диссипация механической энергии, которая зависит от скорости деформаций и уровней напряжений [25]. Частота поперечных колебаний пологих оболочек значительно больше, чем пластин тех же характерных размеров, изготовленных из тех же материалов, поэтому и скорость деформаций в таких искривленных конструкциях много больше, чем в пластинах. Следовательно, выделение тепла в пологих оболочках может оказаться значительно большим, чем в аналогичных пластинах, а значит, и результаты работы [23] нельзя будет дословно переносить на армированные пологие оболочки. Кроме того, используя разные теории изгиба тонкостенных элементов конструкций (Рейсснера, Амбарцумяна, Редди или более точные), можно с разной степенью точности аппроксимировать скорости деформаций и напряжения по толщине КМ-пластин и КМ-оболочек. Так, в работе [22] было продемонстрировано, что для адекватного описания упругопластического (изотермического) динамического поведения армированных пологих оболочек необходимо использовать уточненную теорию их изгиба. Следовательно, применение уточненной теории изгиба может внести существенные поправки в расчет диссипации механической энергии в пологих КМ-оболочках, а значит, и в расчет температурных полей и механического отклика таких конструкций.
Для интегрирования нелинейных задач динамики тонкостенных элементов конструкций используют как явные [5,22,23], так и неявные [8,26] численные схемы.
Согласно вышеизложенному, настоящее исследование посвящено моделированию термоупругопластического динамического поведения армированных
пологих оболочек в рамках уточненной теории изгиба [22]. Нелинейная связанная термомеханическая задача численно интегрируется с использованием явной схемы шагов по времени [5,22,23].
1. Постановка задачи и метод расчета. Рассмотрим деформирование гибкой пологой КМ-оболочки толщиной 2К, с которой свяжем ортогональную систему координат Хг так, что ее срединная поверхность является отсчетной (|жэ| ^ К), а координатные линии х\ и х2 совпадают с линиями главной кривизны этой поверхности хз = 0 (рис. 1, где искривленность панели не изображена в силу ее малости). Оболочка многонаправленно армирована (возможно, и пространственно) К семействами волокон с плотностями армирования Шк (1 ^ к ^ К). В поперечном направлении конструкции Охз структура армирования однородна.
С каждым к-м семейством волокон свяжем ортогональную локальную си-
( к) (к) стему координат х\ , причем ось Ох\ сориентируем вдоль траектории армирования и ее направление зададим двумя углами сферической системы
координат Ok и ^к (рис. 2). Направляющие косинусы № осей Ох(') относи-
Лк)
г]
тельно осей Oxj (i, j = 1, 3, 1 ^ к ^ К) определяются формулами из [23]:
¿Ц = sin 9к cos ¡£к, 1\к) = sin 9к sin <рк, 1\к) = cos вк
¡(к) j(к)
'21 = — sin , ¿22 = COs ¡£к, l3l = — COs вк COs ¡£к, = — COs вк sin !£к,
(к )
(к )
1{ к) = 0 123 = 0,
(1)
® = sin ек, 1 < к < N.
Предполагается, что внешними касательными силами на лицевых поверхностях оболочки можно пренебречь, и в случае пространственной структуры армирования выполняются требования, изложенные в замечании в [23]. (В случаях плоско-перекрестного армирования, как на рис. 1, указанные требования заведомо выполняются.) При этом, согласно [22], усредненные деформации композиции е^ и перемещения точек гибкой пологой оболочки иг
Рис. 1. Элемент волокнистой пологой оболочки [Figure 1. Element of fibrous shallow shell]
Рис. 2. Локальная прямоугольная система координат, связанная с армирующим волокном fc-го семейства [Figure 2. Local cartesian coordinate system associated with fc-th family reinforcing fiber]
в рамках уточненной теории изгиба можно аппроксимировать так (геометрическая нелинейность задачи учитывается в приближении Кармана):
£4 (г, г) = 2 (дгЩ + дуЩ) - хздгд^Ю +
+ Л х^ ( ь х2 Л,. (ш) +д <ш)] + ,1д.д. (2)
+ Ъ Н2 и + 1 - т + 3Лд^3 +дзе" ^ + К + 2^" (2)
т=0
„2 М
ь2-
*«(<,г) = ^х), м = 1,2;
т=0
М хт+1 ь 2 х2
игг) =пг(г,X) -х3дгы + 2 £ (^^ - т+ь-)(<,х),
т=0
и3(г,г) = ю(<,X, ъ = 1,2, х <Ес,
|х3 | ^ ь, < ^ ¿0, X = (х1,х^ , Г = (х1,х2,хз),
(3)
где ю — прогиб; щ — перемещения точек отсчетной поверхности (хз = 0) в тангенциальных направлениях х^; К — главные радиусы кривизны срединной поверхности; <о — начальный момент времени <; — оператор частного дифференцирования по переменной (г = 1, 2); Ь^ —символ Кронекера; М — целое число, которое задает количество слагаемых, удерживаемых в частичных суммах по степеням хз; О — область, занимаемая искривленной панелью в плане. При М = 0 из соотношений (2) и (3) вытекают кинематические гипотезы теорий Амбарцумяна [1] и Редди [6,9]. В равенствах (2) и (3) подлежат
определению двумерные функции ю, щ, (г = 1, 2; 0 ^ т ^ М).
В данном исследовании изучается динамическое поведение пологой КМ-оболочки как гибкой тонкостенной термомеханической системы, следовательно, напряжение азз(<, г) с приемлемой для инженерных расчетов точностью можно линейно аппроксимировать по толщине конструкции [2]:
^3(г. г) = ^ х) -^С-х) х +
+ ^ х) х), X ^хк * о <0, (4)
где (X ^ з {<, х, ±Ь) — нормальные напряжения на верхней (+) и нижней (—) лицевых поверхностях панели, которые известны из силовых граничных условий.
К соотношениям (2)—(4) следует присоединить двумерные уравнения движения гибкой пологой оболочки (см. (10)—(12) в [22]) и определяющие уравнения, связывающие усредненные скорости деформаций е^, напряжений и температуры В в композиции конструкции. Последние уравнения целесообразно представить в матричном виде [23]:
& = Бё + р, (5)
где
& = (¿1,^2,^3,^4 ,&5,&б)Т = (¿11, ¿22, ¿33, ¿23, ¿31, ¿12)Т,
12)
£ = (£ 1, £2, ¿3, ¿4, ¿5, ¿6)Т = (£ 11, ¿22, ¿33, 2ё23, 2¿31, 2ё 12)Т;
(6)
В = (шоЪо + Н-1, р = { - Bg,
^ к=1 '
N N
{ = Шо ро + (рк + 2к Гк), Н = шо1 + Ек, к =1 к =1
N
g = ^ШкГк, Гк = Б-1 Як, Ек = Б-1Ск (1 О ^), к =1
N
Zк = 2к -СкZк, Рк = Ркв (0 О ^Ж), шо = 1 - ;
(7)
к =1
I — единичная бхб-матрица; Б, Ек, Ск, Zк, 2к, 2к — 6х6-матрицы; Б-1, Н-1 — матрицы, обратные 6х6-матрицам Бк, Н; р, Г, g, г к, Як, Р к, Рк — шестиком-понентные векторы-столбцы. Элементы матриц Ск = (с(к)'), Бк = (й^),
2к = (4^), 2к = (4?к)) и вектор-столбцов Як = (яР )), Рк = (Рг(к)) вычисляются по формулам [23]:
= 2$г]Ск + Ак, 2у(/ ) = Ск,
Ч] -
т/-( к ) 3 (к)
Рг(к) = ^ Е ^ + 3^к«к + С- [С* - Г(к)(Т(Мк - Т(вк)Ск)Ак] ,
т=1
(к)
Р(к) = С- № - тГЧтЫСк - т(@к)Ск)Ак], (г,3=1,3, г=4Т6), Ск
^ =Ак8^зк (г „7 = 176),
х Ук Ек п Ек Ек Ак = т^--ГГл-~—г, 2Ск = —-, оКк =
(1 + Ук)(1 - 2Ук)' 1+Ук' 1 - 2
А = 7к Ск К (к) = йКк С(к) = &Ск
Ак =(Ск + СкЬ(к)2, = ^ , ^ = ¿в ,
(к) 9г«(к) Л дт<(к) [ * .(к) .(к),.
ге ) = , ^^к = ^, Хк = ¡^Щ/Р^ЛЪ,
= 0, если Тк < т(к) или Тк = тв(к), Шк < 0, 7к = ^ 1, если Тк = тв(к), Шк > 0,
I
Шк = Ск*Т¿к + т(к)С-1(т(кЧк) - 4к)Ск)в,
36
Тк2 = 2 £ + £ ^ 0 О
6 6
$ = 4? = $ = Е № е ^ = Е $
6 ^ (9)
я[к) = о, 4к) = Е 9^ (р\0 - р\к)), ¿ = 276, 3=1,6, 1 < к < N;
I=1
(к) = (к) = ,(к) ,(к) (к) = (к) = ¡(к), (к) 9ц = 111 = 111 111 , 912 = 112 = 112 112 , ■ ■ ■ ,
п(к) = 2 (к) = 2,(к) ,(к)
916 = 2 Ч16 = 2112 111 , ■ ■ ■,
2 г,(к) = п(к) = 21 (к)](к)
2961 = 461 = 2121 111 , ■ ■ ■,
(к) _ (к) _ Лк)^к) Лк)^к) ,
966 = %6 = 111 122 + 1 12 121 , 1 < к < 1 ;
(10)
ак — коэффициент линейного температурного расширения к-го компонента композиции (к = 0 — связующая матрица, к ^ 1 — волокна к-го семейства);
(к)
Ек, Vк — модуль Юнга и коэффициент Пуассона того же материала; тз = = т(к\хк, В) —предел текучести того же материала при чистом сдвиге, зависящий от параметра Одквиста Хк и температуры В и при активном пластическом деформировании равный интенсивности касательных напряжений Тк; С к — касательный модуль сдвига при пластическом деформировании к-й фазы композиции; —напряжения в том же материале; Бк = (з^), ек = = (ё(к)) —шестикомпонентные векторы-столбцы с компонентами девиатора напряжений в^ и тензора скоростей деформаций е^ в к-м материале композиции, по структуре аналогичные векторам а и ё соответственно (см. (6)); (к)
•р\- — скорости пластических деформаций того же материала; 7к — параметр переключения, который при 7к = 0 определяет термоупругое деформирование, разгрузку или нейтральное нагружение, а при 7к = 1 — активное на-гружение при пластическом деформировании к-го компонента композиции;
р^ — элементы вектора-столбца рк (I = 1, 6, 0 < к < 1); точка — производная по времени ¿; индекс т —операция транспонирования. По повторяющимся индексам к и I в равенствах (8) суммирования нет. Не выписанные в (10)
элементы 6х6-матриц С к = (д(к^) и Qk = (0^), зависящие от направляющих косинусов (1), приведены в табл. 21.40 и 21.44 в [27]. Матрицы Ск и Qk определяют преобразования шестикомпонентных векторов-столбцов скоростей напряжений (Г к = ) и деформаций ек = (ё^) (I = 1,6) в к-м компоненте композиции при переходе от глобальной ортогональной системы
координат Xj к локальной системе х^ (см. рис. 2). Векторы Гк и ек, имеющие соответственно структуру, аналогичную (6), связаны определяющими соотношениями [23]:
Гк = Zfc¿к + Рк, о о <1, (11)
где матрица Zk и вектор-столбец р к вычисляются по формулам (7). В соотношениях (8) учтено, что упругие характеристики к-го материала композиции могут зависеть от температуры В (термочувствительность).
При выводе равенств (5) и (7) попутно получаются следующие структурные матричные соотношения:
¿0 = H-1е - H-1g, ek = Ek£0 + rk, 1 о < N. (12)
Первое выражение (12) определяет скорости деформаций связующей матрицы ¿о через скорости осредненных деформаций композиции ¿, а второе соотношение — скорости деформаций волокон k-го семейства ¿к через скорости деформаций связующего материала ¿0. По формулам же (11) при учете (12)
- -(к) 1
можно выразить и скорости напряжений т\- в k-м компоненте композиции
через осредненные скорости деформаций ( i ,j = 1,3, 0 ^ k ^ N).
В силу структуры векторов-столбцов & и ¿ (см. (6)) из третьего равенства системы (5) определим скорость линейной деформации в направлении хз:
¿33 = — Рз — Ьз1£11 — Ъз2^22 — 2Ьм£23 — 2Ъз^31 — 2Ъз^п) , (13)
где рз, b3i ( i = 1, 6) — элементы вектора-столбца p и матрицы B в соотношении (5); скорость &33 выражается из (4) путем дифференцирования по времени. Скорости деформаций ¿^ в правой части (13) вычисляются за счет дифференцирования по соотношений (2), т. е. определяются через двумерные функции w, W, щ, ¿^ (i = 1, 2, 0 ^ m ^ М).
Используя результаты работы [23], температуру конструкции В аппроксимируем по ее толщине полиномом порядка L:
L
B(t, r) - В0 = Е В(, x)xl3, x <eG, |хз| < h, t^ to, (14) 1=0
где В0 = const — температура естественного состояния искривленной панели; В[ —подлежащие определению двумерные функции (0 ^ I ^ L).
Для завершения формулировки связанной задачи неупругого термомеханического деформирования пологой КМ-оболочки к соотношениям (2)—(14) нужно присоединить уравнения теплового баланса. Так как метрика в искривленной панели с приемлемой для приложений точностью может быть отождествлена с метрикой в декартовой прямоугольной системе координат [1], в качестве двумерных уравнений, описывающих теплопроводность в рассматриваемой пологой КМ-оболочке, можно применять формулы (4.11), (4.13), (4.14) и (4.18) из работы [23]. Кроме того, следует использовать необходимые механические (см. (22)-(24) и (26) в [22]) и теплофизические (см. (3.20) в [23]) граничные условия, а также начальные условия, заданные при t = t0 (см. (25), (26) в [22] и (3.22) в [23]).
Для численного интегрирования поставленной нелинейной начально-краевой задачи используем явную схему шагов по времени, т. е. неизвестные функции будем определять в дискретные моменты времени t = tn (n = 0,1, 2,...). При этом считаем, что в моменты времени tm уже вычислены значения следующих функций [22,23]:
m(x) = w(tm, x), (x) = v!f\tm, x), CTij(r) = <Tij{tm, r),
¿3±)(х) =434^, х), им(х) = и(г)(1п, х), qг(г) = фп, г),
„,, 1
вв(х) = в8(гт, х), в в(х) = в3(1п-1, х), ^(х) = д^^п, х),
^(г) = ^^к)(г) = ^, Г), Х(к)(г) = Х(к)(^, Г), 1 = 1, 2, 1,э = 1, 3, т = п - 1, п, 0 ^ р ^М + 1, 0 - 2,
0 0 ^к^К, х еС, |жэ| < К,
(15)
где
и(р)(г, х) =/ иг(г, г)х^йх3, и(г)(г, х) =/ и (г, ф^, З-н З-н
1 = 1, 2, 0 ^р^М + 1, 0 - 2;
(16)
и — удельная внутренняя энергия композиции оболочки; о^7, е^ —напряжения и деформации в к-й фазе композиции (к = 0 — связующий материал, к ^ 1 — арматура к-го семейства); х(к) —параметр Одквиста в той же фа-
(±)
зе; дг — компоненты вектора теплового потока в композиции; ^ —заданные значения тепловых потоков через нижнюю (-) и верхнюю (+) лицевые поверхности оболочки. Искомые функции ш, иг и (г = 1, 2, 0 ^ т ^ М) в соотношениях (2) и (3) однозначно определяются через кинематические переменные и(р) (см. выражения (16)) за счет матричного равенства (30) в [22].
В механической составляющей рассматриваемой связанной неизотермической упругопластической задачи производные по времени £ аппроксимируем их конечно-разностными аналогами на трехточечном шаблоне {¿га-1, ^п, ¿п+1}. Это позволяет построить явную численную схему [22]. Так, после замены в двумерных уточненных уравнениях динамического равновесия гибкой пологой КМ-оболочки вторых производных по от кинематических переменных ш и и(р) их конечными разностями при учете соотношений (3), (4), (16) и обозначений, аналогичных (15), будем иметь [22]
^ ("Ш1 -2 ш + V) = (м>0 + £ м<0^ ш) -
]=1 ^ г=1 '
V" Р-1 ЛЖ (0) (+) - (-) - М гг +° 33 -а 33 ,
г=1
£("+'? -2и(') + V?) = ±е,(м™ -М-
3 = 1
- 3+ - (-1)г-3-)]^ш-ш£-1)Г« +я-lм<?,
х еС, г = 1, 2, 0 + 1, п = 1, 2, 3,...,
где
К К (D fh
р = PqWq + Е Pk Uk, Uo = 1 , (t, x) = aij(t, r)xl3dx3,
k=1 к=1 ''-h lMg-1\t, x) = I j^aii^, r)xl-1dx3 = | [(a3+ +^3-))(1 - (-1)) + (18)
+ ГТГ ^ -a3-))(l + (-1)i)], U = M, 0 ^l^M + 1;
Po, Pk — объемные плотности связующей матрицы и волокон к-го семейства; т = const > 0 — шаг по времени. Внешние объемные силы в равенствах (17) не учитываются. Осредненные напряжения aij в выражениях (18) связаны
с напряжениями в фазах композиции a(i,j = 1, 3, 0 ^ к ^ К) структурными соотношениями из [23] (по правилу простой смеси).
По формулам (18) с учетом предположений (15) в текущий момент времени tn можно определить силовые факторы М31 и внешние силы a^ в правых частях уравнений (17), поэтому, учитывая необходимые граничные условия [22], из равенств (17) по явной схеме можем вычислить искомые функции
n+1 n+1(l)
W и U i при t = tn+1.
Для численного интегрирования теплофизической составляющей исследуемой задачи также применим явную схему шагов по времени, но на двухточечном шаблоне {tn, tn+{}. При этом с учетом выражений (14), (16) и обозначений типа (15) двумерные уравнения теплового баланса пологой КМ-оболочки примут вид [23]:
P (^V) - Щ 3т)) = -^1 Q 1т) -d2% (т) -Q т + £ М, x eG, 0 ^m ^L - 2, n = 0,1, 2,...;
L
- E(-1)1 hl-1(l\3-] + ha3-))6i(i, x) =
=o
= a3-)(e3~) - в0)+ q3~)(t, x),
L
(19)
(20)
Y,hl-1(i\3V+ha+))ei(t, x)
=o
= а(+)(в3+) - в0) - q3+)(t, x), x eG, to;
Co ^^Я(г + т)в + C ^ +
=o 2 =o =o
+ T ££ £h(i + j + l + m)eie3ei
3 =o =o =o
= U 3m)(t, x), x eG, t^ to, 0 ^m ^L - 2, (21)
где
Н(з) = ---[1 - (-1)*+1 ], , х) = / Ф, (¿ = 1, 3),
8 + 1 З-н
/н
дздз(1, г)х^ = -тк3+) - (-1)т<?3-)] - тд3т-1)(;, х), ■н (22) /" н 1 ^^
Ш(m)(í, х) = аг]ег]х^(1х3, С1(х) = ■ V с\к)ркшк(х) (I = 0,1, 2), ]-Н Рк=0
А33 = Л33|©=©(4)х,±Н)' ^3 = 9оо ;
Х33 — коэффициент теплопроводности композиции в поперечном направлении Ох3, определяемый по структурным формулам из [28]; а(о) —коэффициент теплоотдачи на верхней (+) и нижней (—) лицевой поверхности оболоч-
(о)
ки; Но —температура окружающей среды со стороны той же поверхности; с(к) — коэффициент разложения удельной теплоемкости материала к-й фазы композиции Ск(в — 0°) по формуле (в случае термочувствительности):
Ск (в - в°) = с°к ) + с1к)(в - в°) + с2к)(в - в0)2, 0 <к<К. (23)
Равенства (20) — тепловые граничные условия на лицевых поверхностях искривленной КМ-панели, преобразованные с учетом представления температуры в виде (14). Соотношение (21) выражает двумерные функции
и Ы
(см. (16)) через коэффициенты разложения (14) с учетом равенства (23). Уравнения (20), (21) справедливы для любого момента времени ¿.
Используя соотношения (22) и предположения (15), в текущий момент времени можно вычислить правые части в равенствах (19). Учитывая при этом необходимые граничные (заданные на контуре Г, ограничивающем область С) и начальные тепловые условия [23], получаем из (19) по явной схеме
га+1( )
значения и (т) при Ь = Ьп+1. Далее, в следующий момент времени Ьп+1 из равенств (20), (21) (где правые части уже известны) при учете (22) определяем
п+1
коэффициенты разложения для температуры в I(х), 0 ^ I ^ Ь (см. (14)). В случае термочувствительности компонентов композиции (с^ ) = 0 и/или
с2к) =0 (0 ^ к ^ К) в равенстве (23)) система (20), (21) является нелинейной. Для ее линеаризации может быть использован метод переменных теплофизических параметров, аналогичный методу переменных параметров упругости. В остальных деталях разработанный численный метод (17), (19) при учете соотношений (2), (5)-(13), (18), (22) и структурных формул, полученных в [28] (см. выражения (20) и (22)), реализуется так же, как подробно описано в [22,23]. При этом в соотношениях (5), (7)-(9), (11) и (12) всем величинам следует приписать верхний индекс п (аналогично (17) и (19)), означающий номер текущего дискретного момента времени Ьп. (Величины же, входящие в равенства (10), от п не зависят.) Кроме того, в соотношениях (7) и (8) скорость температуры в необходимо выразить разностной формулой,
вытекающей из формулы трапеций [23]:
П 2 , п п-1/2ч п-1/2 п-1 ТП-1
0 — 0 = -(в - 0 ), 0 — 0 0 , п = 1,2,3,.... (24)
2
Согласно предположениям (15), при t = tn правые части в равенствах (24) известны, т. е. известна и скорость температуры 0 = 0 (tn, r).
В работе [23] показано, что в уравнениях (17), (19) шаг по времени т нужно задавать так, чтобы выполнялось необходимое условие устойчивости Куранта—Фридрихса—Леви [5].
2. Обсуждение результатов расчетов. Рассматривается термоупруго-пластический динамический изгиб относительно тонких пологих цилиндрических КМ-оболочек толщиной 2h = 2 см, имеющих в плане удлиненную прямоугольную форму (G: \xi\ ^ а, \х2\ ^ Ь, а = 3b; 1/Ri — 0, R2 — R = const; b = 50 см). Стрела подъема искривленных панелей f над продольными (в направлении оси Oxi) кромками равна 10 см. Радиус цилиндрической оболочки R при этом определяется формулой [22]: R = (b2 + f2)/(2/), 0 ^ f < 0.4b. По
п (т) п ^ -п
кромкам конструкции имеют жесткое закрепление: w = 0, щ =0, x £ 1, t ^ to (см. (16), (17)). До начального момента времени t = to = 0 панели покоятся (w = 0, Щ™ =0, x £ G, t = t0, i = 1, 2, 0 ^ m ^ M +1) в естественном состоянии при 0 = 00 = 20 °С (x £ G, |хз \ ^ h, t = t0). Начиная с момента времени to оболочки подвергаются воздействию избыточного давления p(t), порожденного воздушной взрывной волной [26]:
(+) (-) _ /,\ _ I pmax ' t/tmax, 0 ^ t ^ tmax,
033 033 — p( ) = |Pmax ■ exp[-a(i - tmax)], t > tmax, (25)
а = - ln(0.01)/(tmin - tmax) > 0, tmin » tmax, p(tmin) = 0.01pmax,
где смысл величин pmax, tmax, imin и а вполне очевиден и описан в [22,23]. В соответствии с экспериментальными данными [26] в расчетах задано max = = 0.1 мс и tmin = 2 мс. Из соотношений (25) получаем: при pmax > 0 конструкции нагружены со стороны нижней (вогнутой) лицевой поверхности
(при этом условно считаем — 0), а при pmax < 0 — со стороны верхней
(выпуклой) поверхности (при этом принимаем ^3) — 0).
Теплообмен оболочек с окружающей средой через их лицевые поверхности реализуется в условиях естественной конвекции (а1-±) = 30 Вт/(м2-K) [29],
q(±) — 0) при температуре воздуха, равной температуре естественного состояния конструкций: 0^±) = 00 (см. (20)). На торцевых поверхностях цилиндрических панелей задана температура, которая поддерживается равной 00.
Оболочки выполнены из эпоксисвязующего [17], армированного стеклянными волокнами [16] (стеклопластик), или магниевого сплава ВТ65 [30], усиленного стальной проволокой У8А [16] (металлокомпозиция). Упругопласти-ческое поведение компонентов композиции при постоянной температуре 0 и активном нагружении определяется билинейной диаграммой растяжения -сжатия:
Е(к)е, И < е{к) = о{к)/Е(к),
sign(e)aik) + Е{к)(е - sign(e)е{3к)), \е\ > e{sk), 0 < к < К,
t
где а, е — осевое напряжение и соответствующая ему деформация; Е(к) = = Е(к)(В), Е^ = Е^(В) — модули Юнга и линейного упрочнения материала
к-й фазы композиции; = а(к)(Э) —предел текучести этого же компонента. Физико-механические характеристики материалов фаз композиций представлены в таблице, где и — коэффициент Пуассона, Л — коэффициент теплопроводности, а — коэффициент линейного теплового расширения, с — удельная теплоемкость, а в скобках приведено значение температуры (В, С), при которой была определена данная характеристика. В расчетах использованы линейные аппроксимации по В для всех физико-механических характеристик по данным, указанным в таблице.
Для дискретизации задачи по пространственным переменным х\ и использована равномерная сетка с шагами Дх1 = Дх2 = 6/50, а шаг по времени т задавался равным 0.25 мкс. Необходимые условия устойчивости разработанной численной схемы при этом выполняются с запасом (см. соотношения (6.3) в [23]).
Пологие оболочки армированы двумя (К = 2) семействами непрерывных волокон, ориентированных по направлениям Ох1 и Ох2 (см. рис. 1), с плотностями (¿1 =0.1 и Ш2 = 0.3 соответственно. Углы армирования при этом имеют следующие значения (см. рис. 2): О1 = О2 = к/2, =0 и ^>2 = к/2.
В [23] были проведены расчеты для гибких КМ-пластин с рассматриваемой структурой армирования при тех же характерных размерах тонкостенных конструкций (Н, а, Ь и / = 0). При этом в соотношениях (2) и (3) принималось значение М = 0 (использовалась теория изгиба Амбарцумяна [1]),
Физико-механические характеристики материалов композиций [16,17,30] [Physico-mechanical characteristics of the components of composite [16,17,30]]
Характеристика материала [Material characteristic] Эпокси-связующее [Epoxy] Стекловолокно [Fiberglass] Магниевый сплав BT65(Mg) [VT65(Mg) magnesium alloy] Стальная проволока У8А [U8A steel wire]
p, kg/m3 1210.0 (20 °C) 1208.0 (40 C) 2520.0 (20 °C) 2519.6 (80 C) 1800.0 (20 °C) 1796.2 (100 C) 7800.0 (20 C) 7791.8 (100 C)
E, GPa 2.8 (20 C) 2.3 (40 C) 86.8 (20 C) 86.3 (80 C) 43.0 (20 C) 38.5 (100 C) 210.0 (20 C) 195.0 (100 C)
0.33 (20 C) 0.333 (40 C) 0.25 (20 C) 0.254 (80 C) 0.330 (20 C) 0.334 (100 C) 0.3 (20 C) 0.305 (100 C)
<js, MPa 20 (20 C) 18 (40 C) 4500 (20 C) 4400 (80 C) 267 (20 C) 219 (100 C) 3968 (20 C) 3971 (200 C)
Es, GPa 1.114 (20 C) 0.783 (40 c) 6.230 (20 C) 5.168 (80 C) 0.379 (20 C) 0.367 (100 C) 6.973 (20 C) 5.014 (200 C)
A, W/(m • K) 0.243 (20 C) 0.240 (40 C) 0.89 (20 C) 0.86 (80 C) 117.23 (20 C) 121.42 (100 C) 42.7 (20 C) 41.7 (100 C)
a • 106, K-1 68.1 (20 C) 70.3 (40 C) 2.5 (20 C) 2.6 (80 C) 20.9 (20 C) 22.6 (100 C) 12.3 (20 C) 13.2 (100 C)
c, kJ/(kg • K) 1.54 (20 C) 1.60 (40 c) 0.800 (20 C) 0.839 (80 C) 1.032 (20 C) 1.054 (100 C) 0.485 (20 C) 0.488 (100 C)
а в разложении температуры (14) варьировалось значение Ь ^ 2. Было показано, что в динамически изгибаемых КМ-пластинах температурные поля следует определять по формуле (14) при задании в них Ь = 6 или Ь = 7.
В работе же [22] было продемонстрировано, что в случае изотермического изгибного динамического деформирования пологих КМ-оболочек для адекватного описания их неупругого механического поведения в кинематических соотношениях (2) и (3) нужно принимать не М = 0, а М = 6 или М = 7 (т. е. следует применять уточненную теорию изгиба).
Как уже отмечалось ранее, в данной работе исследуется влияние использования уточненных соотношений (2) и (3) (при М = 7) на результаты расчета термомеханического отклика динамически изгибаемых искривленных КМ-панелей при упругопластическом деформировании материалов компонентов их композиций.
В связи с этим на рис. 3 изображены зависимости максимальных значений температуры Эт(¿; М) = шахв(£, г; М) в рассматриваемых пологих КМ-
г
оболочках. При этом на основании результатов, полученных в [23], в формуле (14) было принято Ь = 7. Кривые 1 и 2 рассчитаны соответственно при М = 0 (теория Амбарцумяна) и М = 7 (уточненная теория изгиба) для стек-лопластиковой панели, нагруженной снизу (рис. 3, а) при ртах = 8 МПа (см. (25)) и сверху (рис. 3, Ь) при ртах = —8 МПа, а также металлокомпозит-ной конструкции, нагруженной сверху при ртах = —30 МПа (рис. 3, с). Для случая нагружения металлокомпозитной оболочки снизу (ртах = 30 МПа) зависимости вт(£; М) здесь не приводятся, так как кривые 1 и 2 при этом на рассматриваемом интервале времени 0 ^ Ь ^ 24 мс практически совпадают (незначительное различие наблюдается лишь при £ > 22 мс).
Сравнение кривых 1 и 2 на рис. 3 показывает, что температурные поля, рассчитанные по простейшей неклассической теории изгиба искривленной КМ-панели — теории Амбарцумяна (см. кривые 1) —и по уточненной теории изгиба (см. кривые 2), с течением времени все более различаются. Следовательно, есть основания предполагать, что теория Амбарцумяна, широко используемая в расчетной практике, в случае динамического термоупругопла-стического деформирования тонкостенных КМ-конструкций типа оболочек может оказаться непригодной для проведения таких расчетов. Окончательное же решение этого вопроса возможно только после сравнения с результатами соответствующих экспериментов, которые автору неизвестны. (Здесь следует отметить, что дополнительные расчеты показали: в случаях квазистатического нагружения и динамического термоупругого деформирования пологих КМ-оболочек результаты, полученные по теории Амбарцумяна и по уточненной теории изгиба, практически совпадают.)
Сделанный выше вывод подтверждают и расчеты поперечных колебаний рассматриваемых панелей. Так, на рис. 4 представлены осцилляции прогиба центральной точки стеклопластиковой пологой КМ-оболочки (1о(Ь ) = = , 0, 0)), определенные при ртах = 8 МПа (рис. 4, а) и ртах = —8 МПа (рис. 4, Ь). Кривые 1 и 2 на рис. 4 получены при тех же значениях М, что и на рис. 3, а кривые 2' на рис. 4 изображены для сравнения: они рассчитаны по уточненной теории изгиба и соответствуют изотермическому упру-гопластическому случаю, когда изменение температуры в пологой КМ-обо-лочке не учитывается. Поведение кривых 1 и 2 на рис. 4 свидетельствует
О 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32
Time t, msec
а
55 F ,
Р 50 ®S 45 а 40 | 35 t 30
I 25 20
О 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
Time t, msec
Ь
о160 ° „140
ф 120
2100
| 80
I 60
¡3 40 20
0 2 4 6 8 10 12 14 16 18 20 22 24
Time i, msec
с
Рис. 3. Зависимость от времени максимального значения температуры в стеклопластико-вой пологой оболочке, нагруженной снизу (а) и сверху (b), а также в металлокомпозитной
панели, нагруженной сверху (с)
[Figure 3. Time dependence of the maximum temperature value in a fiberglass shallow shell loaded from below (a) and from above (b), as well as in a metal-composite panel loaded from
above (с)]
о том, что с увеличением времени зависимости wq(t), рассчитанные по теории Амбарцумяна (кривые 1) и по уточненной теории (кривые 2), существенно различаются. Аналогичное сопоставление кривых 2 и 2' демонстрирует, что с увеличением t зависимости Wo(t), определенные с учетом (кривые 2) и без учета (кривые 2') температурного отклика в искривленной стеклопластико-вой панели, все более отличаются друг от друга. Согласно этим результатам, в еще большей степени различаются напряженно-деформированные состояния в компонентах композиции таких тонкостенных конструкций, рассчитанные по теории Амбарцумяна (М = 0) и уточненной теории изгиба (М = 7), а также с учетом и без учета температурных полей, возникающих в них. Соответствующие зависимости здесь не представлены.
3
3
2'
а 1
о 1
о
3 -1
-7
О
5
10 15 20 25 30 35 40 45 50
Time t, msec
1
,У
-ю
о
5
10 15 20 25 30 35 40 45 50 Time t, msec
Ъ
Рис. 4. Осцилляции прогиба центральной точки стеклопластиковой пологой оболочки, нагруженной снизу (а) и сверху (b), рассчитанные по разным теориям изгиба
[Figure 4. Deflection oscillations of the central point of a fiberglass shallow shell loaded from below (a) and from above (b) calculated using different bending theories]
На рис. 5 изображены аналогичные кривые 1, 2 и 2', рассчитанные для ме-таллокомпозитной пологой оболочки при pmax = 30 МПа (рис. 5, а) и pmax = = -30 МПа (рис. 5, b). Как видно из рис. 5, b, при нагружении искривленной Mg-УБА-панели сверху уже на рассматриваемом малом интервале времени (0 ^ t ^ 24 мс) кривые 1 и 2, а также 2 и 2' значительно различаются между собой. Кривые 1 и 2 на рис. 5, а (нагружение снизу) визуально незначительно различаются между собой лишь при t ~ 23.5 мс. Объяснятся это тем, что, как отмечалось выше, при нагружении металлокомпозитной пологой оболочки снизу зависимости Qm(t; М), рассчитанные по теории Амбарцумяна (кривая 1) и по уточненной теории изгиба (кривая 2), заметно различаются только при t > 22 мс. Однако сравнение кривых 2 и 2' на рис. 5, а свидетельствует о том, что даже в этом случае нагружения расчеты Mg-УБА-панели, выполненные по уточненной теории изгиба с учетом (кривая 2) и без учета (кривая 2') теплового отклика в ней, приводят к значительно различающимся зависимостям wo(t).
Несмотря на то, что на рис. 5, а кривые 1 и 2 визуально почти не различаются, деформированные состояния компонентов композиции, рассчитанные по теории Амбарцумяна и по уточненной теории изгиба, в этом случае на рассматриваемом интервале времени различаются весьма заметно. Так, на рис. 6 изображены осцилляции максимальных значений интенсивности деформаций £*fc) (£<m(t; М) = maxe*fc) (t, r; M)) связующего материала (рис. 6, а; k = 0)
и арматуры второго семейства (рис. 6, b; к = 2), которая испытывает наи-
r
4
■6 с-'-*-'-1-*-'-1-'-1-'-1-'-1-'-1-'-1-'-1-'-1-'-1
О 2 4 6 8 10 12 14 16 18 20 22 24
Time t, msec
а
О 2 4 6 8 10 12 14 16 18 20 22 24
Time i, msec
Ь
Рис. 5. Осцилляции прогиба центральной точки металлокомпозитной пологой оболочки, нагруженной снизу (а) и сверху (b), рассчитанные по разным теориям изгиба
[Figure 5. Deflection oscillations of the central point of a metal-composite shallow shell loaded from below (a) and from above (b) calculated using different bending theories]
большее деформирование в металлокомпозитной панели, нагруженной снизу. Кривые 1 и 2 на рис. 6 получены при тех же условиях, что и аналогичные кривые на рис. 5, а. Кривые 1 и 2 на рис. 6 (особенно на рис. 6, b) показывают, что расчет деформированного состояния материалов фаз композиции даже в этом случае необходимо проводить по уточненной теории изгиба (кривые 2). Действительно, поведение кривой 1 на рис. 6, b (расчет по теории Амбарцумяна) свидетельствует о том, что наибольшее значение интенсивности деформаций в волокнах второго семейства einax = max £т (t) достигается
при t = 19.1 мс и равно 4.99%, а согласно поведению кривой 2 на этом ри-
(2)
сунке (уточненная теория изгиба), значение £m&x = 4.78 % достигается при t = 6.7 мс.
Согласно поведению кривой 2 на рис. 3, с, температура в металлокомпо-зитной панели, нагруженной сверху, в отдельных точках в некоторые моменты времени может достигать порядка 150-170 С. При нагружении же такой пологой КМ-оболочки снизу температура не превышает 70С. (Как уже отмечалось, соответствующие зависимости Qm(t; М) здесь не приводятся.) Целесообразно предположить, что высокий уровень нагрева динамически изгибаемой металлокомпозитной конструкции обусловлен тем, что в проведенных расчетах предполагалось: теплообмен КМ-панели с окружающей средой через лицевые поверхности реализуется при естественной конвекции воздуха (а(±) = 30 Вт/(м2 ■ K)). Однако при нагружении искривленной панели взрывной волной она осциллирует в поперечном направлении с высокой частотой (см. кривые на рис. 5), поэтому теплообмен через ее лицевые поверхности
^ 0 '-1-1-'-1-'-1-'-1-'-1-1-1-'-1-1-1-1-1-'-1-1-1-'-1
" О 2 4 6 8 10 12 14 16 18 20 22 24
Time t, msec
а
0
1 о -»-,-,-,-,-1-,-,-,-1-,-1-»-,-,-,-,-,-,-,-,-1-,-1
О 2 4 6 8 10 12 14 16 18 20 22 24
Time t, msec
Ь
Рис. 6. Осцилляции максимальных значений интенсивности деформаций связующего материала (а) и волокон второго семейства (b) в металлокомпозитной пологой оболочке, нагруженной снизу, рассчитанные по разным теориям изгиба [Figure 6. Oscillations of the maximum values of the intensity of strains of the binder material (a) and fibers of the second family (b) in a metal-composite shallow shell loaded from below, calculated according to different bending theories]
следовало бы моделировать не при естественной, а при вынужденной конвекции. Согласно этому допущению, были выполнены дополнительные расчеты при значении коэффициентов теплоотдачи а( ± = 300 Вт/(м2 ■ K) (см. равенства (20)), которое характеризует вынужденную конвекцию для газов [29]. Определенные таким образом значения Qm(t; М) отличаются от полученных ранее, т. е. при естественной конвекции (см. рис. 3), лишь в четвертой значащей цифре. А значит, увеличение коэффициентов теплоотдачи на лицевых поверхностях изучаемых искривленных КМ-панелей в 10 раз практически не оказывает влияния на результаты расчетов их неизотермического упругопластического поведения на рассматриваемых временных интервалах (0 ^ t ^ 50 мс). Следовательно, на таких, достаточно малых, интервалах времени теплофизические процессы в пологих КМ-оболочках можно считать адиабатическими, т. е. можно не учитывать кондуктивную составляющую процесса теплопереноса в них.
Заключение. Разработанная математическая модель неизотермического упругопластического поведения гибких пологих оболочек с произвольными структурами армирования непрерывными волокнами позволяет аппроксимировать тангенциальные перемещения и температуру в поперечном направлении с разной степенью точности.
Сравнительный анализ термоупругопластической динамики искривленных КМ-панелей показал, что в отличие от гибких армированных пластин [23]
как стеклопластиковые, так и металлокомпозитные пологие оболочки, нагруженные в поперечном направлении, необходимо рассчитывать, учитывая тепловой отклик в них. При этом целесообразно использовать не простейшую неклассическую теорию изгиба таких тонкостенных элементов конструкций (теорию Амбарцумяна [1,23]), а уточненную теорию изгиба [22].
Выполненные расчеты показали, что при нагружении поперечной взрывной волной относительно тонкие искривленные панели из стеклопластика в отдельных точках кратковременно (менее 1 мс) могут дополнительно нагреваться на 14-34 С (аналогичные по структуре и характерным размерам стеклопластиковые пластины нагреваются не более чем на 14 С [23]), а металло-композитные пологие оболочки в отдельных точках способны дополнительно нагреваться на 5G С при их нагружении со стороны вогнутой лицевой поверхности и на 15G С при нагружении со стороны выпуклой лицевой поверхности (аналогичные металлокомпозитные гибкие пластины нагреваются не выше, чем на 3G-32 С [23]). Более интенсивное тепловыделение и, как результат, более активное воздействие температурных полей на механический отклик всегда наблюдается при динамическом нагружении искривленных КМ-пане-лей со стороны их выпуклой лицевой поверхности.
При нагружении пологой КМ-оболочки взрывной волной с любой лицевой стороны наблюдается ее прощелкивание в сторону вогнутости. Максимальный модуль прогиба конструкции при таком прощелкивании, которое реализуется позже прекращения действия нагрузки, в случае нагружения пологой оболочки со стороны ее выпуклой лицевой поверхности значительно превосходит аналогичную величину при нагружении искривленной панели со стороны ее вогнутой лицевой поверхности. Это обстоятельство и объясняет более интенсивное тепловыделение в первом случае. Наибольшее значение температуры в отдельных точках пологой КМ-оболочки достигается именно в момент ее прощелкивания.
Конкурирующие интересы. Конкурирующих интересов не имею. Авторский вклад и ответственность. Я несу полную ответственность за предоставление окончательной версии рукописи в печать. Окончательная версия рукописи мною одобрена.
Финансирование. Работа выполнена в рамках государственного задания (№ госрегистрации 121030900260-6).
Библиографический список
1. Амбарцумян С. А. Общая теория анизотропных оболочек. М.: Наука, 1974. 446 с.
2. Богданович А. Е. Нелинейные задачи динамики цилиндрических композитных оболочек. Рига: Зинатне, 1987. 29б с.
3. Куликов Г. М. Термоупругость гибких многослойных анизотропных оболочек // Изв. РАН. МТТ, 1994. №2. С. 33-42. EDN: SLRUCB.
4. Bannister M. Challenges for composites into the next millennium — a reinforcement perspective// Compos. Part A: Appl. Sci. Manuf., 2001. vol.32, no. T. pp. 901-910. DOI: https://doi.org/10.1016/S1359-835X(01)00008-2.
б. Абросимов H. А., Баженов В. Г. Нелинейные задачи динамики композитных конструкций. Н. Новгород: ННГУ, 2002. 400 с.
6. Reddy J. N. Mechanics of Laminated Composite Plates and Shells. Theory and Analysis. Boca Raton: CRC Press, 2004. xxiii+831 pp. DOI: https://doi.org/10.1201/b12409.
7. Qatu M. S., Sullivan R. W., Wang W. Recent research advances on the dynamic analysis of composite shells: 2000-2009// Compos. Struct., 2010. vol.93, no. 1. pp. 14-31. DOI: https://doi.org/10.1016/j.compstruct.2010.05.014.
8. Kazanci Z. Dynamic response of composite sandwich plates subjected to time-dependent pressure pulses // Int. J. Non-Linear Mech., 2011. vol. 46, no. 5. pp. 807-817. DOI: https:// doi.org/10.1016/j.ijnonlinmec.2011.03.011.
9. Андреев А. Упругость и термоупругость слоистых композитных оболочек. Математическая модель и некоторые аспекты численного анализа. Saarbrücken: Palmarium Academic Publ., 2013. 93 с. EDN: QZAPNP.
10. Vasiliev V. V., Morozov E. Advanced Mechanics of Composite Materials and Structural Elements. Amsterdam: Elsever, 2013. xii+412 pp. EDN: UERHXD. DOI: https://doi.org/10. 1016/C2011-0-07135- 1.
11. 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. DOI: https://doi.org/10.1007/s11465-013-0262-x.
12. Соломонов Ю. С., Георгиевский В. П., Недбай А. Я., Андрюшин В. А. Прикладные задачи механики композитных цилиндрических оболочек. М.: Физматлит, 2014. 408 с. EDN: UGLCQJ.
13. Moriniere F. D., Alderliesten R. C., Benedictus R. Modelling of impact damage and dynamics in fibre-metal laminates — A review// Int. J. Impact Eng., 2014. vol. 67. pp. 27-38. DOI: https://doi.org/10.1016/j.ijimpeng.2014.01.004.
14. Gibson R. F. Principles of Composite Material Mechanics. Boca Raton: CRC Press, 2016. xxiii+700 pp. DOI: https://doi.org/10.1201/b19626.
15. Димитриенко Ю. И. Механика композитных конструкций при высоких температурах. М.: Физматлит, 2019. 448 с https://www.rfbr.ru/rffi/ru/books/o_2079143#1.
16. Композиционные материалы: Справочник/ ред. Д. М. Карпинос. Киев: Наук. думка, 1985. 592 с.
17. Handbook of Composities/ ed. G. Lubin. New York: Van Nostrand Reinhold, 1982.
18. Leu S.-Y., Hsu H.-C. Exact solutions for plastic responses of orthotropic strain-hardening rotating hollow cylinders// Int. J. Mech. Sci., 2010. vol.52, no. 12. pp. 1579-1587. DOI: https://doi.org/10.1016/j.ijmecsci.2010.07.006.
19. Vena P., Gastaldi D., Contro R. Determination of the effective elastic-plastic response of metal-ceramic composites // Int. J. Plasticity, 2008. vol. 24, no. 3. pp. 483-508. DOI: https://doi.org/10.1016Zj.ijplas.2007.07.001.
20. Brassart L., Stainier L., Doghri I., Delannay L. Homogenization of elasto-(visco) plastic composites based on an incremental variational principle // Int. J. Plasticity, 2012. vol. 36. pp. 86-112. DOI: https://doi.org/10.1016/j.ijplas.2012.03.010.
21. Ахундов В. М. Инкрементальная каркасная теория сред волокнистого строения при больших упругих и пластических деформациях// Мех. композ. матер., 2015. Т. 51, №3. С. 539-558.
22. Янковский А. П. Уточненная модель упругопластического изгибного деформирования гибких армированных пологих оболочек, построенная на основе явной схемы типа «крест»// Вычисл. мех. сплош. сред, 2017. Т. 10, №3. С. 276-292. EDN: ZHZVIN DOI: https://doi.org/10.7242/1999-6691/2017.10.3.22.
23. Янковский А. П. Моделирование термоупруговязкопластического деформирования гибких армированных пластин// ПММ, 2022. Т. 86, №1. С. 121-150. EDN: EKGCRN DOI: https://doi.org/10.31857/S003282352201009X.
24. Reissner E. On transverse vibrations of thin shallow elastic shells// Quart. Appl. Math., 1955. vol.13, no.2. pp. 169-176. DOI: https://doi.org/10.1090/qam/69715.
25. Грешнов В. М. Физико-математическая теория больших необратимых деформаций металлов. М.: Физматлит, 2018. 232 с.
26. Houlston R., DesRochers C. G. Nonlinear structural response of ship panels subjected to air blast loading// Comp. Struct., 1987. vol. 26, no. 1-2. pp. 1-15. DOI: https://doi.org/ 10.1016/0045-7949(87)90232-X.
27. Малмейстер А. К., Тамуж В. П., Тетерс Г. А. Сопротивление жестких полимерных материалов. Рига: Зинатне, 1972. 500 с.
28. Янковский А. П. Моделирование процессов теплопроводности в пространственно-армированных композитах с произвольной ориентацией волокон // Прикл. физ., 2011. №3. С. 32-38. ЕБЫ: ЫШУНг.
29. Луканин В. Н., Шатров М. Г., Камфер Г. М. [и др.] Теплотехника/ ред. В. Н. Луканин. М.: Высш. шк., 2003. 671 с. ЕБЫ: дИНУВН.
30. Безухов Н. И., Бажанов В. Л., Гольденблат И. И. [и др.] Расчеты на прочность, устойчивость и колебания в условиях высоких температур / ред. И. И. Гольденблат. М.: Машиностроение, 1965. 567 с.
Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki
[J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2023, vol. 27, no. 1, pp. 119-141
d https://doi.org/10.14498/vsgtu1958
ISSN: 2310-7081 (online), 1991-8615 (print)
MSC: 74K20
Modeling of non-isothermal elastic-plastic behavior of reinforced shallow shells in the framework of a refined bending theory
A. P. Yankovskii
Khristianovich Institute of Theoretical and Applied Mechanics,
Siberian Branch of the Russian Academy of Sciences,
4/1, Institutskaya st., Novosibirsk, 630090, Russian Federation.
Abstract
The dynamic problem of non-isothermal and inelastic deformation of flexible shallow multidirectionally reinforced shells is formulated in the frameworks of the refined theory of bending. The temperature is approximated by a 7th order polynomial over the thickness of constructions. The geometric nonlinearity of the problem is modeled by the Karman approximation. The solution of the formulated coupled nonlinear two-dimensional problem is obtained using an explicit numerical scheme. The thermo-elastic-plastic response of fiberglass and metal-composite cylindrical elongated panels with an orthogonal reinforcement structure, loaded frontally with an air blast wave, has been studied. It is shown that, unlike reinforced plates similar in structure and characteristic dimensions, shallow shells under intense short-term loading must be calculated taking into account the occurrence of temperature fields in them. In this case, the refined theory of bending of curved panels should be used instead of the simplified version (the non-classical theory of Ambartsumyan). The temperature increment at separate points of shallow fiberglass shells can reach 14-34 °C, and in similar metal-composite panels can reach 50-150 C. Cylindrical shallow shells are more intensively deformed when they are loaded by an air blast wave from the side of a convex front surface.
Keywords: flexible shallow shells, multidirectional reinforcement, dynamic deformation, coupled thermo-elastic-plasticity, refined bending theory, explicit numerical scheme.
Received: 26th September, 2022 / Revised: 28th December, 2022 / Accepted: 13th February, 2023 / First online: 27th March, 2023
Mechanics of Solids Research Article
© Authors, 2023
© Samara State Technical University, 2023 (Compilation, Design, and Layout) 9 ©® The content is published under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/) Please cite this article in press as:
Yankovskii A. P. Modeling of non-isothermal elastic-plastic behavior of reinforced shallow shells in the framework of a refined bending theory, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2023, vol. 27, no. 1, pp. 119-141. EDN: YRWNPW. DOI: 10.14498/vsgtu1958 (In Russian). Author's Details:
Andrei P. Yankovskii https://orcid.org/0000-0002-2602-8357
Dr. Phys. & Math. Sci.; Leading Research Scientist; Lab. of Fast Processes Physics;
e-mail: lab4nemir@rambler .ru
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. Ambartsumian S. A. Obshchaia teoriia anizotropnykh obolochek [The General Theory of Anisotropic Shells]. Moscow, Nauka, 1974, 446 pp. (In Russian)
2. Bogdanovich A. E. Nelineinye zadachi dinamiki tsilindricheskikh kompozitnykh obolochek [Nonlinear Problems of the Dynamics of Cylindrical Composite Shells]. Riga, Zinatne, 1987, 295 pp. (In Russian)
3. Kulikov G. M. Thermoelasticity of flexible multilayer anisotropic shells, Mech. Solids, 1994, vol. 29, no. 2, pp. 27-35. EDN: TVBRLF.
4. Bannister M. Challenges for composites into the next millennium — a reinforcement perspective, Compos. Part A: Appl. Sci. Manuf., 2001, vol. 32, no. 7, pp. 901-910. DOI: https:// doi.org/10.1016/S1359-835X(01)00008-2.
5. Abrosimov N. A., Bazhenov V. G. Nelineinye zadachi dinamiki kompozitnykh konstruktsii [Nonlinear Problems of Dynamics Composites Designs]. Nizhniy Novgorod, Nizhniy Novgorod State Univ., 2002, 400 pp. (In Russian)
6. Reddy J. N. Mechanics of Laminated Composite Plates and Shells. Theory and Analysis. Boca Raton, CRC Press, 2004, xxiii+831 pp. DOI: https://doi.org/10.1201/b12409.
7. Qatu M. S., Sullivan R. W., Wang W. Recent research advances on the dynamic analysis of composite shells: 2000-2009, Compos. Struct., 2010, vol. 93, no. 1, pp. 14-31. DOI: https:// doi.org/10.1016/j.compstruct.2010.05.014.
8. Kazanci Z. Dynamic response of composite sandwich plates subjected to time-dependent pressure pulses, Int. J. Non-Linear Mech., 2011, vol.46, no. 5, pp. 807-817. DOI: https:// doi.org/10.1016/j.ijnonlinmec.2011.03.011.
9. Andreev A. Uprugost' i termouprugost' sloistykh kompozitnykh obolochek. Matematicheskaia model' i nekotorye aspekty chislennogo analiza [Elasticity and Thermoelasticity of Layered Composite Shells. Mathematical Model and Some Aspects of Numerical Analysis]. Saarbrücken, Palmarium Academic Publ., 2013, 93 pp. (In Russian). EDN: QZAPNP.
10. Vasiliev V. V., Morozov E. Advanced Mechanics of Composite Materials and Structural Elements. Amsterdam, Elsever, 2013, xii+412 pp. EDN: UERHXD. DOI: https://doi.org/10. 1016/C2011-0-07135- 1.
11. 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. DOI: https://doi.org/10.1007/s11465-013-0262-x.
12. Solomonov Yu. S., Georgievskii V. P., Nedbai A. Ya., Andryushin V. A. Prikladnye zadachi mekhaniki kompozitnykh tsilindricheskikh obolochek [Applied Problems of Mechanics of Composite Cylindrical Shells]. Moscow, Fizmatlit, 2014, 408 pp. (In Russian)
13. Moriniere F. D., Alderliesten R. C., Benedictus R. Modelling of impact damage and dynamics in fibre-metal laminates — A review, Int. J. Impact Eng., 2014, vol. 67, pp. 27-38. DOI: https://doi.org/10.1016Zj.ijimpeng.2014.01.004.
14. Gibson R. F. Principles of Composite Material Mechanics. Boca Raton, CRC Press, 2016, xxiii+700 pp. DOI: https://doi. org/10.1201/b19626.
15. Dimitrienko Yu. I. Mekhanika kompozitnykh konstruktsii pri vysokikh temperaturakh [Mechanics of Composite Structures under High Temperatures]. Moscow, Fizmatlit, 2019, 448 pp. (In Russian). https://www.rfbr.ru/rffi/ru/books/o_2079143#1
16. Kompozitsionnye materialy: Spravochnik [Composite Materials: Handbook], ed. D. M. Kar-pinos. Kiev, Nauk. Dumka, 1985, 592 pp. (In Russian)
17. Handbook of Composities, ed. G. Lubin. New York, Van Nostrand Reinhold, 1982.
18. Leu S.-Y., Hsu H.-C. Exact solutions for plastic responses of orthotropic strain-hardening rotating hollow cylinders, Int. J. Mech. Sci., 2010, vol.52, no. 12, pp. 1579-1587. DOI: https://doi.org/10.1016/j.ijmecsci.2010.07.006.
19. Vena P., Gastaldi D., Contro R. Determination of the effective elastic-plastic response of metal-ceramic composites, Int. J. Plasticity, 2008, vol. 24, no. 3, pp. 483-508. DOI: https:// doi.org/10.1016/j.ijplas.2007.07.001.
20. Brassart L., Stainier L., Doghri I., Delannay L. Homogenization of elasto-(visco) plastic composites based on an incremental variational principle, Int. J. Plasticity, 2012, vol. 36, pp. 86-112. DOI: https://doi.org/10.1016Zj.ijplas.2012.03.010.
21. Akhundov V. M. Incremental carcass theory of fibrous media under large elastic and plastic deformations, Mech. Compos. Mater., 2015, vol.51, no. 3, pp. 383-396. EDN: MCFUEJ. DOI: https://doi.org/10.1007/s11029-015-9509-4.
22. Yankovskii A. P. A refined model of elastic-plastic bending deformation of flexible reinforced shallow shells based on explicit "cross" scheme, Comput. Cont. Mech., 2017, vol.10, no. 3, pp. 276-292 (In Russian). EDN: ZHZVIN DOI: https://doi.org/10.7242/1999-6691/2017. 10.3.22.
23. Yankovskii A. P. Modeling of thermoelastic-visco-plastic deformation of flexible reinforced plates, Mech. Solids, 2022, vol.57, no. 7, pp. 111-133. DOI: https://doi.org/10.3103/ S0025654422070184.
24. Reissner E. On transverse vibrations of thin shallow elastic shells, Quart. Appl. Math., 1955, vol.13, no. 2, pp. 169-176. DOI: https://doi.org/10.1090/qam/69715.
25. Greshnov V. M. Physico-Mathematical Theory of High Irreversible Strains in Metals. Boca Raton, CRC Press, 2019, xii+242 pp. DOI: https://doi.org/10.1201/9780429259791.
26. Houlston R., DesRochers C. G. Nonlinear structural response of ship panels subjected to air blast loading, Comp. Struct., 1987, vol.26, no. 1-2, pp. 1-15. DOI: https://doi.org/ 10.1016/0045-7949(87)90232-X.
27. Malmeister A. K., Tamuzh V. P., Teters G. A. Soprotivlenie zhestkikh polimernykh materi-alov [Resistance of Rigid Polymeric Materials]. Riga, Zinatne, 1972, 500 pp. (In Russian)
28. Yankovskii A. P. Simulation of processes of heat conductivity in spatially reinforced composites with an arbitrary orientation of fibers, Prikl. Fiz., 2011, no. 3, pp. 32-38 (In Russian). EDN: NURWHZ.
29. Lukanin V. N., Shatrov M. G., Kamfer G. M., et al. Teplotekhnika [Heat Engineering], ed. V. N. Lukanin. Moscow, Vyssh. Shk., 2003, 671 pp. (In Russian). EDN: QMHYSH.
30. Bezukhov N. I., Bazhanov V. L., Gol'denblat I. I., et al. Raschety na prochnost', ustoichivost' i kolebaniia v usloviiakh vysokikh temperatur [Calculations for Strength, Stability Vibrations under High Temperatures], ed. I. I. Gol'denblat. Moscow, Mashinostroenie, 1965, 567 pp. (In Russian)