Numerical analysis of dynamic strength of composite cylindrical shells under multiple-pulse exposures

Abstract


The target of this research is studying the fiberglass cylindrical shellswith exposed ends resulted from the cross-winding of tapes made of an unidirectional composite material. The aim of the study was to develop a numerical technique to model a progressive fracture of laminated composite cylindrical shells under multiple-pulse loading with an internal pressure of a various intensity. Kinematic model of deforming the laminate package is based on the applied geometrically nonlinear theory of shells. The formulation of geometric dependencies is based on the relations of the simplest quadratic variant of the nonlinear elasticity theory. The physical relations of the elementary layer are formulated based on the generalized Hooke's law for the orthotropic material based on the hypotheses of the applied shells theory. The process of a progressive shell failure is described within the degradation model of stiffness characteristics in elementary layers of a multilayer package which is based on Hoffman's criteria for composite materials and on the criterion of maximum stresses for the fibers. The process of damage accumulation in the shell material due to a multiple application of impulse load is takeninto account by means of the computational scheme in which the calculation of the current stress strain state is carried out with stiffness characteristics obtained in the model of their degradation under previous loading. Energetically consistent system of motion equations of the applied shell theory is deduced from the stationarity condition related to the functional of the shell total energy. A numerical method for solving the formulated initial-boundary value problem is based on an explicit variational difference scheme. The accuracy of the considered techniqueis proved by comparing the obtained results with the known experimental data. The results of how the number of loadings affects the value of marginal circumferential deformationsare presented. It is established that the level of maximalring deformations is approximately ten times less than their limit values compared to a single loading.

Full Text

Введение Композитные материалы благодаря высокой удельной прочности находят все более широкое применение при создании различных конструкций, испытывающих интенсивные динамические нагрузки. В частности, к таким конструкциям относятся взрывозащитные камеры, для которых требуется с высокой степенью надёжности исключить их разрушение в условиях многократного воздействия импульсной нагрузки. Экспериментально установлено [1], что динамическая прочность и деформируемость стеклопластиковых цилиндрических оболочек существенно зависят от уровня импульсной нагрузки и числа нагружений. Кроме того, при повторном нагружении наблюдалось снижение величины предельной деформации оболочки [2]. Данные эффекты, связанные с неоднократным нагружением, изучались экспериментально при нагрузках, близких к предельным, что приводило к значительным повреждениям после первого нагружения и полному разрушению всего слоистого пакета при последующих нагружениях. Для обоснованного применения стеклопластиков при проектировании взрывозащитных контейнеров необходимо провести исследование их усталостной прочности при многократных импульсных воздействиях, которые существенно меньше их предельных значений. Некоторые аспекты экспериментального изучения малоцикловой усталостной прочности стеклопластиков рассматривались в [3]. Обзор критериев прогрессирующего разрушения многослойных композитов содержится в работе [4]. Различные варианты моделей прогрессирующего разрушения слоистых композитов рассматривались в работах [5-15]. Предлагаемая работа посвящена численному моделированию прогрессирующего разрушения стеклопластиковых цилиндрических оболочек при многократном взрывном нагружении. 1. Постановка и метод решения задачи Рассматривается стеклопластиковая цилиндрическая оболочка, изготовленная перекрестной намоткой лент однонаправленного композитного материала. Отнесем оболочку к ортогональной криволинейной системе координат , совпадающей с линиями главных кривизн и внешней нормалью к внутренней поверхности оболочки. Коэффициенты Ламе и главные кривизны: где R - радиус внутренней поверхности оболочки. Кинематическая модель деформирования многослойного пакета базируется на прикладной геометрически нелинейной теории оболочек [16], в рамках которой распределение перемещений по толщине оболочки запишется в виде (1) где - перемещения точек внутренней поверхности в направлениях осей ; - углы поворота нормали к внутренней поверхности. Компоненты нелинейного тензора деформаций прикладной теории цилиндрических оболочек можно представить в виде [17] (2) где (3) Символ здесь и далее обозначает, что из записанного соотношения можно получить еще одно одновременной заменой 1 на 2, 2 на 1. Физические соотношения элементарного слоя с учетом гипотез прикладной теории оболочек запишутся в виде (4) где - ступенчато изменяющиеся функции переменной , которые вычисляются через модули упругости и коэффициенты Пуассона элементарного слоя композитного материала. В процессе нестационарного деформирования допускается локальное разрушение элементарных слоев в пакете многослойного композита. В результате разрушения отдельных элементарных слоев композита происходит перераспределение напряжений между слоями, а многослойный пакет продолжает оказывать сопротивление дальнейшему деформированию. Процесс прогрессирующего разрушения оболочки описывается в рамках модели деградации жесткостных характеристик элементарных слоев многослойного пакета, которая базируется на критерии Хоффмана [18] для композитного материала: и на критерии максимальных напряжений для волокон , , вычисленных в осях ортотропии элементарного слоя. Здесь - пределы прочности однонаправленного композитного материала и волокон при растяжении и сжатии. В процессе деформирования и трещинообразования связующего предполагается жесткое сцепление между соседними слоями. В зависимости от знака объемной деформации и моды разрушения композитного материала и волокон модули упругости монослоя с трещинами принимают одно из возможных значений, приведенных в таблице. Схема редуцирования модулей упругости слоя Reduction scheme of layer elasticity moduli Характер разрушения Условия на значение объемной деформации 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 Здесь , , - значения модулей упругости и сдвига в результате деградации материала оболочки, вычисляемые по формулам [19] , ; , , - значения модулей упругости в начальном неповрежденном состоянии; , - коэффициенты Пуассона; - текущие значения деформаций ; , - максимальные значения деформаций при растяжении за предысторию деформирования. Процесс накопления поврежденности в материале оболочки, обусловленный многократным приложением импульсной нагрузки, учитывается посредством вычислительной схемы, в которой расчет текущего напряженно-деформированного состояния осуществляется с жесткостными характеристиками, полученными в рамках модели их деградации при предыдущем нагружении. Будем считать, что исчерпание несущей способности оболочки происходит в момент обнуления в каком-либо поперечном сечении всех жесткостных характеристик. Энергетически согласованная система уравнений движения прикладной теории оболочек выводится из условия стационарности функционала полной энергии оболочки. Вариационное уравнение движения цилиндрической оболочки со свободными торцами, нагруженной импульсом внутреннего давления, запишется в виде (5) где где - нагрузка по направлению координатной оси ; S - область интегрирования на внутренней поверхности оболочки; - эффективные жесткостные характеристики слоя, которые вычисляются через модули упругости и коэффициенты Пуассона [16]; - координаты слоев, отсчитываемые от внутренней поверхности оболочки. Минимизация функционала полной энергии оболочки (5) позволяет получить систему уравнений движения (6) и естественные граничные условия (7) Дополняя соотношения (6), (7) необходимым числом начальных условий (8) получим полную систему уравнений для анализа нелинейных волновых процессов деформации и послойного разрушения цилиндрических оболочек, полученных перекрестной намоткой однонаправленного композитного материала. Численный метод решения сформулированной задачи основывается на явной вариационно-разностной схеме [20]. 2. Результаты исследования Рассматривалась оболочка радиусом внутренней поверхности R = 0,15 м, толщиной = 0,00819 м и длиной L = 4R, выполненная из стеклопластика со следующими жесткостными и прочностными характеристиками элементарного слоя и волокон[21]: = 55,0 ГПа, = 9,04 ГПа, = 0,276, = 0,3, = = 5,08 ГПа, = 5,2 ГПа, = 2010 кг/м3, = 1,5 ГПа, = 0,397 ГПа, = = 0,0267 ГПа, = = 0,0719 ГПа, = = 0,0357 ГПа, = 0,04 ГПа, = 2,8 ГПа, = 0,7425 ГПа. В расчетах импульс внутреннего давления задавался с помощью эмпирической зависимости где m - масса заряда; = 4,77×106 Дж/кг - теплотворная способность взрывчатого вещества (ВВ); - расстояние от центра заряда до точки внутренней поверхности оболочки. Достоверность рассматриваемого подхода подтверждалась сравнением численных расчетов с экспериментальными данными [2] по значениям максимальных кольцевых деформаций в центральном сечении неразрушившейся оболочки и периоду радиальных колебаний T, полученным для оболочки со структурой пакета при нагружении однократным зарядом ВВ массой m = 0,0206 кг. В эксперименте , T×106 = 180 с, в расчете , T×106 = 190 с.Наблюдается достаточно хорошее совпадение результатов как по амплитуде, так и периоду колебаний. Далее рассматривалась задача динамической прочности стеклопластиковой цилиндрической оболочки толщиной = 0,0102 м со структурой армирования при многократных импульсных воздействиях, которые существенно меньше их предельных значений. Типичные зависимости радиальных колебаний в центральном сечении оболочки показаны на рис. 1 для пяти последовательных нагружений зарядом ВВ массой m = 0,105 кг. Видно, что с ростом числа нагружений наблюдается увеличение амплитуды и периода колебаний, обусловленное снижением жесткостных характеристик многослойного пакета вследствие разрушения некоторых элементарных слоев. Рис. 1. Осциллограммы окружных деформаций в центральном сечении при многократном нагружении Fig. 1. Waveform circumferential deformations in the central section with repeated loading Иллюстрация процесса разрушения оболочки при многократном нагружении зарядом ВВ массой m = 0,105 кг представлена на рис. 2-4. На рис. 2 для двух моментов времени показаны деформированные конфигурации оболочки с нанесенными зонами разрушения внешнего слоя. t = 1,0×10-4 с t = 7,0×10-4 с t = 1,0×10-4 с t = 7,0×10-4 с t = 1,0×10-4 с t = 7,0×10-4 с а б в Рис. 2. Деформированные конфигурации оболочки (прогибы увеличены в пять раз): а - первое нагружение; б - четвертое нагружение; в - пятое нагружение Fig. 2. The deformed configuration of the shell (deflections magnified five times): a - first loading; b - loading the fourth; с - fifth loading На рис. 3 показана соответствующая кинетика развития зон разрушения по толщине пакета слоев в центральном сечении оболочки для пяти последовательных нагружений. Здесь светлый фон соответствует зонам, в которых разрушения отсутствуют; серый фон - зоны разрушения матрицы; в зонах черного фона происходит разрушение волокон. Процесс разрушения стеклопластиковых оболочек при многократных импульсных воздействиях аналогичен ранее полученным результатам при однократном нагружении [22], что свидетельствует о слабом влиянии скорости деформации на характер возникающих в композите повреждений. д г в б а t = 1,0×10-4 с t = 4,0×10-4 с t = 7,0×10-4 с Рис. 3. Кинетика зон разрушения по толщине оболочки: а - первое нагружение; б - второе нагружение; в - третье нагружение; г - четвертое нагружение; д - пятое нагружение Fig. 3. Kinetics of fracture zones across the shell thickness: a - first loading; b - second loading; с - loading the third; d - fourth loading; е - the fifth loading На рис. 4 показаны экспериментальные [3] и расчетные зависимости предельных значений окружных деформаций от логарифма числа нагружений N. Анализ полученных результатов показал, что для данного типа стеклопластиковых оболочек кольцевая деформация линейно зависит от логарифма числа нагружений до разрушения. Разброс расчетных данных для при lgN = 0 объясняется тем, что результаты получены для различных масс ВВ, при которых разрушение происходило при однократном нагружении. Кроме того, существует граничный уровень предельной нагрузки, при котором разрушение не наблюдается, а число нагружений может быть достаточно большим. Отмеченный эффект может быть связан с тем, что рассматриваемая модель не учитывает накопление повреждений в материале при нагрузках меньше критических. При этом уровень максимальных кольцевых деформаций ниже Рис. 4. Максимальные окружные деформации в зависимости от числа нагружений: - разрушение; - нет разрушения; 1 - эксперимент; 2 - расчет Fig. 4. The maximum circumferential strain depending on the number of loadings: - destruction; - no destruction; 1 - experiment; 2 - calculation их предельных значений по сравнению с однократным нагружением примерно в десятьраз. Полученные результаты удовлетворительно коррелируют с экспериментальными данными [3]. Заключение Предложенная расчетная модель позволяет с удовлетворительной точностью описывать процесс послойного прогрессирующего разрушения стеклопластиковых цилиндрических оболочек при многократных импульсных воздействиях. Значения предельных окружных деформаций при многократном импульсном нагружении примерно на порядок меньше по сравнению с соответствующими предельными значениями при однократном нагружении, что удовлетворительно согласуется с результатами экспериментальных испытаний.

About the authors

N A Abrosimov

Research Institute of Mechanics of National Research Nizhny Novgorod State University n.a. N.I. Lobachevsky

A V Elesin

Research Institute of Mechanics of National Research Nizhny Novgorod State University n.a. N.I. Lobachevsky

References

  1. Деформация и разрушение цилиндрических оболочек из стеклопластика при внутреннем импульсном нагружении / В.И. Цыпкин [и др.] // Механика композитных материалов. - 1981. - №2. - С. 249-255.
  2. Федоренко А.Г., Сырунин М.А., Иванов А.Г. Влияние структуры армирования ориентированных стеклопластиков на прочность круговых цилиндрических оболочек при взрывном нагружении изнутри // Механика композитных материалов. - 1991. - № 4. - С. 631-640.
  3. Сырунин М.А., Федоренко А.Г., Иванов А.Г. Динамическая прочность цилиндрических оболочек из стеклопластика при многократном взрывном нагружении // Физика горения и взрыва. - 1997. - Т. 33. № 6. - С. 102-107.
  4. Mark R. Garnich and Venkata M.K. Akula Review of degradation models for progressive failure analysis of fiber-reinforced polymer composites // Appl. Mech. Rev. - 2009. - Vol. 62. - No. 1. - P. 1-33. doi: 10.1115/1.3013822
  5. Veera Skytta, Olli Saarela, Markus Wallin Progressive failure of composite laminates; analysis vs experiments // Fracture of Nano and Engineering Materials and Structures. Proceedings of the 16th European Conference of Fracture, Alexandroupolis, Greece, July 3-7, 2006. - Part B. - P. 341-342. doi: 10.1007/1-4020-4972-2_168
  6. Elisa Pietropaoli Progressive failure analysis of composite structures using a constitutive material model (USERMAT) developed and implemented in ANSYS © // Applied Composite Materials. - 2012. - Vol. 19. - No. 3. - P. 657-668. doi: 10.1007/s10443-011-9220-0
  7. Computational modeling of complex failure mechanisms in laminates / Meer F.P. van der, Sluys L.J., Hallett S.R., Wisnom M.R. // J. Composite Materials. - 2012. - Vol. 46. - No. 5. - P. 603-623. doi: 10.1177/0021998311410473
  8. Gunel M., Kayran A. Non-linear progressive failure analysis of open-hole composite laminates under combined loading // Journal of Sandwich Structures & Materials. - 2013. - Vol. 15. - No. 3. - P. 309-339. doi: 10.1177/1099636213483651
  9. Ellul B., Camilleri D., Betts J.C. A progressive failure analysis applied to fiber-reinforced composite plates subject to out-of-plane bending // Mechanics of Composite Materials. - 2013. - Vol. 49. - No. 6. - P. 605-620. doi: 10.1007/s11029-013-9377-8
  10. Chen X., Li Z., Wang H. Progressive failure analysis of an open-hole composite laminate by using the s-version finite-element method // Mechanics of Composite Materials. - 2014. - Vol. 50. - No. 3. - P. 279-294. doi: 10.1007/s11029-014-9414-2
  11. Shokrieh M.M., Karamnejad A. A investigation of strain rate effects on the dynamic response of a glass/epoxy composite plate under blast loading by using the finite-difference method // Mechanics of Composite Materials. - 2014. - Vol. 50. - No. 3. - P. 295-310. doi: 10.1007/s11029-014-9415-1
  12. Liu P.F., Xing L.J., Zheng J.Y. Failure analysis of carbon fiber/epoxy composite cylindrical laminates using explicit finite element method // Composites Part B-Engineering. - 2014. - Vol. 56. - P. 54-61. doi: 10.1016/j.compositesb.2013.08.017
  13. Jayashree Sengupta, Arghya Ghosh, Dipankar Chakravorty Progressive failure analysis of laminated composite cylindrical shell roofs // Journal of Failure Analysis and Prevention. - 2015. - Vol. 15. - No. 3. - P. 390-400. doi: 10.1007/s11668-015-9951-6
  14. Finite element analysis of progressive failure and strain localization of carbon fiber/epoxy composite laminates by ABAQUS / P.F. Liu, Y.H. Yang, Z.P. Gu, J.Y. Zheng // Applied Composite Materials. - 2015. - Vol. 22. - No. 6. - P. 711-731. doi: 10.1007/s10443-014-9432-1
  15. Shroff S., Kassapoglou C. Progressive failure modelling of impacted composite panels under compression // Journal of Reinforced Plastics and Composites - 2015. - Vol. 34. - No. 19. - P. 1603-1614. doi: 10.1177/0731684415592485
  16. Васильев В.В. Механика конструкций из композиционных материалов. - М.: Машиностроение, 1988. - 272 с.
  17. Шаповалов Л.А. Об учете поперечного обжатия в уравнениях нелинейной динамики оболочек // Изв. РАН МТТ. - 1997. - № 3. - С. 156-168.
  18. Браутман Л., Крок Р. Композиционные материалы. Т. 2. Механика композиционных материалов. - М.: Мир, 1978. - 568 с.
  19. Композиционные материалы: справочник / В.В. Васильев [и др.]; под общ. ред. В.В. Васильева, Ю.М. Тарнопольского. - М.: Машиностроение, 1990. - 512 с.
  20. Абросимов Н.А., Баженов В.Г. Нелинейные задачи динамики композитных конструкций. - Н. Новгород: Изд-во Нижегород. гос. ун-та, 2002. - 400 с.
  21. Предельная деформируемость и прочность базальтопластиковых оболочек при внутреннем взрывном нагружении / В.Н. Русак [и др.] // ПМТФ. - 2002. - Т. 43, № 1. - С. 186-195.
  22. Абросимов Н.А., Елесин А.В., Пирогов С.А. Численный анализ неосесимметричного деформирования и прогрессирующего разрушения слоистых композитных цилиндрических оболочек при импульсном нагружении // Проблемы прочности и пластичности. - 2015. - Т. 77, № 1. - С. 23-32.

Statistics

Views

Abstract - 272

PDF (Russian) - 74

Cited-By


PlumX


Copyright (c) 2016 Abrosimov N.A., Elesin A.V.

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