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

Statistics

Views

Abstract - 222

PDF (Russian) - 57