Full Text

Введение Помимо эффективных упругопрочностных характеристик многие полимерные композитные материалы (ПКМ) обладают высоким удельным поглощением энергии при ударном разрушении. Это свойство важно в специальных демпфирующих элементах конструкций, которые в процессе разрушения поглощают кинетическую энергию. Удельное поглощение энергии является важным и достаточно общим параметром материала, оцениваемым на практике [1-3]. При проектировании зачастую пренебрегают зависимостью жесткостных и прочностных свойств полимерных композитов от скорости деформации, которая, однако, часто обнаруживается в экспериментальных исследованиях и рассматривается в моделях [4-15]. В результатах экспериментальных исследований по ударам в торец трубчатых образцов из слоистых композитов можно найти довольно противоречивые выводы для разных скоростей ударника. В работах [16, 17] отмечается обратная связь удельной поглощенной энергии и скорости удара, в то время как результаты [18, 19] говорят об ее увеличении со скоростью, а в [1, 20] изменения оказались незначительными. Результаты испытаний также чувствительны к ориентации армирующих волокон и проявляют масштабный эффект при варьировании размеров образца [21]. Таким образом, возникает потребность в модели, позволяющей учесть описанные эффекты при различных скоростях деформации, что позволит сократить объем экспериментальных исследований. 1. Определяющие соотношения для ортотропного слоистого композита с учетом скоростного упрочнения Определяющие соотношения для ортотропного композитного материала с повреждением волокна и связующего могут быть представлены в следующей форме [22-26]: (1) где параметры поврежденности и соответствуют формам повреждения, ассоциированным с волокном и матрицей. Индекс “0” означает использование констант неповрежденного материала. При этом (i = 1,2) соответствует полному разрушению, - исходному материалу. Будем использовать критерий максимальных напряжений, для которого в случае статического нагружения материал остается неповрежденным внутри заданной области и повреждается на ее границе: (2) где - предельные напряжения разрушения при сжатии в направлении волокна, получаемые из стандартных механических испытаний; - предельные напряжения при растяжении в направлении волокна; - предельные напряжения при сжатии в поперечном направлении; - предельные напряжения при растяжении в поперечном направлении; S - предельные напряжения сдвига, которые предполагаются равными для всех компонент сдвиговых напряжений. Отметим, что предельные значения для и приняты равными, что не будет ограничивать последующие рассуждения. При имеющихся экспериментальных данных для и могут быть приняты различные предельные напряжения. То же касается и предположения о пределах сдвиговых напряжений, одинаковых для всех компонент. При достижении одной из компонент напряжений предельных значений (2) начинается процесс накопления повреждений, который описывается диаграммами нагружения (рис. 1). Рис. 1. Общий вид диаграмм нагружения: а - вдоль волокна; б - в поперечном направлении; в - сдвиг Fig. 1. General form of loading diagrams, (a) along the fiber, (b) across the fiber, (c) shear В случае нагружения вдоль волокна после достижения предельных напряжений происходит резкое падение кривой, отражающее хрупкость разрушения, которое определяется разницей между деформациями инициации повреждений и полным разрушением (рис. 1, а). В случае нагружения в поперечном направлении и сдвига участок кривой 1-2 становится параллельным оси абсцисс (рис. 1, б), что соответствует способности связующего воспринимать нагрузку после достижения предельных напряжений, наблюдаемое экспериментально [24]. Диаграммы для компонент , и аналогичны представленным на рис. 1, б и рис. 1, в. Для перехода к случаю динамического нагружения, когда проявляется скоростное упрочнение материала, вводится зависимость предельных напряжений (2) от параметров {ψ1, ψ2, dψ1/dt, dψ2/dt}, т.е. от параметров поврежденности и скорости их изменения: (3) Предположим отсутствие зависимости от скорости нагружения предельных напряжений при одноосном растяжении в направлении волокон и поперечном растяжении , выявленной в [27, 28]. Каждую из функций , и S представим в виде произведения двух множителей, отвечающих за статическую и динамическую компоненту: (4) Точный вид функций , , , и S из (3) имеет сходство с соотношениями скоростного упрочнения Джонсона-Кука [29], однако представлены они не через деформации и скорость деформации, а через параметры поврежденности и скорость их изменения: (5) где параметры , , , , , , , , , Ci, Ni и (индекс i соответствует рассматриваемому виду нагружения: X, Y, S) подбираются таким образом, чтобы обеспечить удовлетворительное приближение к экспериментальным ударным и квазистатическим диаграммам деформирования. Для этого для различных скоростей деформаций записываются соотношения, которые рассматриваются как обыкновенные дифференциальные уравнения на функции поврежденности : (6) где верхний индекс “0” у упругих констант означает использование исходного значения неповрежденного материала. В левой части уравнений (6) записаны выражения для напряжений, которые можно определить из экспериментальных кривых, а справа - функции, которыми предлагается их описать. Путем серии численного решения перебором определяются значения констант из правой части, которые обеспечивают соответствие экспериментальным кривым. При численном решении динамической задачи из (6) после каждого шага по времени определены параметры и , которые используются в (1) для определения приращения деформаций. 2. Моделирование ударного нагружения трубчатого композитного образца Представленная в предыдущем разделе модель была реализована в комплексе конечно-элементного моделирования Abaqus для численного решения практических задач. В качестве примера тестовой задачи рассмотрим композитный трубчатый образец, ударяемый по свободному краю грузом массой 200 кг с начальной скоростью 5 м/c. Общая толщина стенки составляет 1,92 мм и состоит из восьми слоев равной толщины с ориентацией [90°/0°6/90°] (рис. 2), так что ось цилиндра направлена вдоль оси X, а слои материала укладываются по нормали к поверхности Z. Таким образом, слои с ориентацией 0° направлены вдоль оси образца, а слои с ориентацией 90° направлены по окружности. Образец имеет скос на переднем крае (фаску) под углом 30°. Слои образованы композитом на основе углеволокна и эпоксидного связующего IM6G/3501-6, параметры для определяющих соотношений (5) которого были подобраны в [30] и представлены в табл. 1-3. h = 1,92 мм Регулярная зона Зона фаски Рис. 2. Схема эксперимента по удару в торец трубчатого образца Fig. 2. Scheme of the impact experiment to the end of the tubular specimen Таблица 1 Параметры модели для нагружения сдвига Table 1 Parameters of the model for shear AS, MПа BS, MПа nS CS s-1 NS , MПа 75 76 1,5 0,06 2×10-5 1,5 7500 Таблица 2 Параметры модели для сжатия в поперечном направлении Table 2 Parameters of the model for transverse compression , MПа , MПа CY s-1 NY , MПа 155,2 135 0,2 10-8 2×10-6 5,9 13 500 Таблица 3 Параметры модели для сжатия в продольном направлении Table 3 Parameters of the model for longitudinal compression , MПа CX , s-1 NX , MПа 14295 13,5 6,2×10-6 2×10-5 0,11 171000 На рис. 3 показан общий вид объектов моделирования в Abaqus и трехмерные конечные элементы типа C3D8R, общее количество которых составляет 120 000. Одному слою по толщине (1,92 мм) соответствует один элемент. При этом два других размера элемента варьируются от 0,4 мм в зоне вблизи ударяемого края до 3 мм с обратной стороны. Между поверхностью ударника и элементами образца заданы условия механического контакта. Расчеты проводились с использованием модуля Abaqus Explicit с явной схемой интегрирования по временно́му шагу. Рис. 3. Общий вид разбиения объекта на конечные элементы Fig. 3. General view of the finite element mesh В процессе моделирования поврежденные элементы после потери жесткости могут испытывать экстремальные деформации, что влечет за собой уменьшение устойчивого шага интегрирования. Чтобы избежать этой проблемы, используется критерий удаления элемента: достижение 3 % деформаций по направлению вдоль волокна либо потеря 95 % жесткости из-за повреждения связующего, что может быть интерпретировано как почти полное разрушение. Кроме того, имеет место вырожденность соотношений (1) при нулевых значениях параметров поврежденности, которая также разрешается введением численного ограничения на их минимальное значение. На рис.4 показано распределение параметра после прохождения ударником 3,2 мм с момента удара, что соответствует размеру фаски вдоль оси цилиндра. ψ2 а ψ2 б Рис. 4. Параметр поврежденности связующего после прохождения ударником 3,2 мм после контакта с образцом: а - образец с фаской; б - образец без фаски Fig. 4. Damage parameter to the matrix after the impactor passing of 3.2 mm after the contact with the sample, (a) the sample with a chamfer, (b) the sample without a chamfer Как можно заметить, для образца с фаской повреждение распространяется существенно дальше перед фронтом разрушения, чем для образца с фаской. Это может объяснять снижение силы реакции при ударе у образца с фаской, так как у такого образца происходит большее повреждение материала перед фронтом разрушения, что снижает жесткость. Распределение не показано, так как повреждение волокна затрагивает лишь несколько элементов перед поверхностью ударника. На рис. 5 показаны диаграммы ударного нагружения для образца с концевой фаской, а для сравнения проведено моделирование образца без фаски. Как можно заметить, отсутствие фаски вызывает достижение силы реакции, более чем в два раза превышающее среднее значение после стабилизации процесса. В то же время разрушение образца без фаски происходит вблизи контакта с ударником, и сила реакции довольно быстро стабилизируется и осциллирует примерно около одинаковых значений с образцом с фаской. Рис. 5. Сила реакции в зависимости от перемещения ударника Fig. 5. Reaction intensity depending on the impactor displacement Удельное поглощение энергии на основе перемещения ударника на длину , при котором разрушается материал массой , рассчитывается по формуле , где - сила реакции; М - масса ударника; и - ускорение ударника и ускорение свободного падения соответственно. При использовании образца со скосом на свободном крае удельная энергия, рассчитанная для = 14 мм, равна 103 Дж/г, а для образца без скоса - 117 Дж/г. Заключение В работе было проведено моделирование ударного нагружения трубчатых образцов из слоистых композитов с учетом зависимости прочностных характеристик от скорости деформирования. Отмечено влияние скоса на ударяемом крае образца в момент удара, которое, однако, уменьшается после стабилизации процесса разрушения.

About the authors

A. N Fedorenko

Skolkovo institute of science and technology

B. N Fedulov

Skolkovo institute of science and technology; Lomonosov Moscow State University

E. V Lomakin

Lomonosov Moscow State University

References

Statistics

Views

Abstract - 286

PDF (Russian) - 149