Full Text

Введение Методы, основанные на МРП, получили широкое распространение для оценки влияния накопленных повреждений на кривые зависимости напряжение-деформация. Последние исследования на эту тему включают работы Такеды и др. [1], Табией и Иванова [2], Сиды и др. [3]. Однако все эти исследования не в полной мере отвечают на вопросы, где в материале локализуются зоны повреждения и как они развиваются при нагрузке. Проблема развития повреждений (мониторинг распространения повреждения в материале) была отражена, в частности, в одной из последних работ Чжао и др. [4], где развитие повреждений в композите при растяжении исследуется с использованием конечно-элементного анализа. Применяются различные сочетания критериев прочности для волокна и матрицы и избирательный метод ухудшения прочности из работы Блэкеттера и др. [5]. Показано, что локализация первоначального повреждения, а также уровни нагрузки для начала разрушения и окончательного разрушения очень чувствительны к выбору критерия разрушения матрицы. Аналогичные вычисления были проведены Николетто и Ривой [6] для расчёта кривых зависимости деформации от напряжения и определения уровней напряжения, связанных с конкретными механизмами разрушения в ходе развития повреждений в тканях с саржевым переплетением при растяжении. У Такано и др. [7] предлагается метод трёхмерного многоуровневого моделирования развития повреждений в мезоструктурах композиционных материалов, армированных тканым и вязаным текстильным материалом, при растяжении. Авторы рассматривают растяжение до 1 % и принимают критерий напряжения Хоффмана [8] и модель разрушения Мураками-Оно [9, 10]. Тан и Уиткомб [11] исследуют чувствительность поведения материала при разрушении к структуре ткани. Проводится анализ разрушения с использованием критерия максимального напряжения и усовершенствованной модели разрушения Блэкеттера. Обнаружено, что текстильные материалы с большой волнистостью показывают сложное поведение, которое зависит от выбора модели разрушения, а композиты с низкой волнистостью подобной чувствительностью не обладают. У Зако и др. [12] моделирование развития повреждений методом конечных элементов выполняется для плоских композитов, армированных тканым материалом, при растяжении в направлении волокна. Используя критерий максимального напряжения и модель разрушения Мураками-Оно, авторы показывают, что зона повреждений начинается с поперечного растрескивания в зонах перекрещивания жгутов волокон, а затем прогрессирует в утке, пока не соединится с соседними зонами повреждений в областях перекрещивания. Модели, представленные в литературе, главным образом касаются композитов с взаимно перпендикулярным расположением волокон, которые нагружаются в направлении волокна. Это важно отметить, поскольку в таких моделях отсутствуют значительные сдвиговые деформации. В работах Ломова и др. [13, 14] авторы внедряли метод МРП для анализа развития повреждений в композитах с трёхосным плетением при растягивающей нагрузке, где хотя бы одно семейство волокон было ориентировано под углом 45° к направлению нагрузки. К основным недостаткам численных методик, основанных на МРП, относится сильная чувствительность кривых деформации-напряжения к размеру сетки и увеличению нагрузки [15]. Текущая работа направлена на определение пределов прочности трехмерно-армированного полимерного углепластика с эпоксидным связующим, отработку методики построения, наиболее приближенной к реальной геометрии волокон, и сравнение полученных результатов с данными, полученными путем томографического и акустического анализа, для дальнейшего получения идеализированного виртуального прототипа. Виртуальный прототип является цифровым макетом, может относительно быстро генерироваться и содержит достаточную информацию, чтобы позволить идентифицировать и решить большую часть потенциальных проблем не прибегая к трудоемким этапам проектирования. Для реализации одной из разновидностей метода МРП на языке APDL был разработан макрос. Макрос может использоваться как для однонаправленного, так и для тканого композиционного материала. Верификация метода выполнена ранее в работах [16, 17]. 1. Получение идеализированного виртуального прототипа В качестве объекта исследования выбран трехмерно-армированный композиционный материал из высокопрочного углеродного волокна и полимерной матрицы. Как уже отмечалось, не всегда имеется возможность определить свойства материала экспериментальным путем. В данном случае невозможно провести эксперимент на сдвиг в плоскости xz и yz из-за малой толщины материала, не позволяющей изготовить образцы необходимого размера. Сечение армирующего каркаса материала до пропитки связующим приведено на рис. 1. Армирующий каркас исследуемого материала до пропитки связующим представляет собой объемно извитую структуру с диагональным переплетением ортогональных волокон. Рис. 1. Фотография преформы Для определения эффективных свойств и пределов прочности композиционного материала выделен элемент периодичности, геометрия которого построена на основе томографии образцов (рис. 2). Геометрия элемента периодичности без матрицы строилась в программе TexGen, при этом вводилось допущение: площадь поперечного сечения идеализированного волокна равна средней поперечной площади соответствующих волокон, вычисленной по результатам томографии. Ввиду особенностей данного продукта в ходе его применения в ряде случаев геометрия генерируется с самопересечением нитей. С целью исключения самопересечений нитей решена контактная задача. Для этого в ANSYS WorkBench верс. 15 построена конечно-элементная модель, где на каждой грани волокна задается контактная поверхность. На противолежащих параллельных гранях элемента периодичности заданы уравнения связи в перемещениях (Coupled Sets в терминах ANSYS). а б X Z Y в г Рис. 2. Геометрия плетения образцов: а - по результатам томографии; б - идеализированный виртуальный прототип в плоскости сечения xz; в - конечно-элементная сетка на волокнах; г - на матрице Далее к координатам узлов начальной конечно-элементной сетки добавляются приращения, полученные из решения контактной задачи. Обновленные координаты являются начальной конфигурацией, исключающей самопересечение нитей, для решения основной задачи. На базе начальных положений узлов внешних граней нитей в дальнейшем, в NX, построены геометрические поверхности и объемы. Объёмная доля связующего в композите составляет 48,9 %, а объемная доля связующего в полученной геометрии составляет 49,47 %. Конечно-элементная сетка с матрицей и волокнами имеет 32,5 млн степеней свободы (рис. 2, в). Процедура численного определения эффективных упругих свойств основана на методике, апробированной и детально описанной в работе [18]. При этом свойства нити вычислялись с допущением, что она является однонаправленным композиционным материалом с объемным содержанием волокна 84% по эмпирической формуле Чамиса (Chamis) [19]. Общая математическая постановка задачи выглядит следующим образом: уравнение равновесия (1) обобщенный закон Гука (2) геометрические соотношения с учетом малых деформаций (3) где - тензоры напряжений, деформаций, перемещений и упругих постоянных соответственно, а - объемные силы. где - области, занятые армирующими нитями и связующим соответственно. Пределы прочности определялись с помощью вышеописанного метода механики рассеянных повреждений. Для моделирования разрушения на языке APDL был разработан макрос, который на каждом шаге нагружения анализирует напряжения в элементах. Если в каком-то элементе напряжения по какому-либо из направлений превышают предел прочности, то соответствующий материал считается разрушенным и в нем понижается модуль упругости. В рамках настоящей работы разрушение оценивалось по критерию максимальных напряжений: , (4) где - значение критерия максимальных напряжений, меняющееся от 0 до 1, 1 - соответствует разрушению; - растягивающие напряжения в соответствующем направлении; - сжимающие напряжения в соответствующем направлении; - предел прочности по напряжениям в соответствующих направлениях при растяжении или сжатии. В разрушенных элементах упругие свойства понижаются в 100 и 10 раз для волокна и матрицы соответственно. Свойства материала понижаются мгновенно. Примененная модель деградации представлена в табл. 1. Таблица 1 Модель деградации волокна Занижаемые свойства материала Направление, по которому был превышен предел прочности x y z xy yz xz Ex Ey Ez νxy νyz νxz Gxy Gyz Gxz На каждом шаге приращения граничных перемещений процесс деградации свойств в подвергнутых разрушению элементах производится единожды, новые свойства используются на следующем шаге нагружения. 2. Испытание на растяжение Расчетный график зависимости силы реакции от перемещений траверсы с помощью метода МРП при растяжении в направлении основы представлен на рис. 3. Рис. 3. Графики зависимости силы реакции от перемещений При численном анализе на ячейке периодичности были заданы следующие граничные условия: где - размеры представительного объема (см. рис. 2) по соответствующим направлениям. Начало разрушения наблюдается в матрице в направлении основы при перемещении 0,511 мм, что соответствует нагрузке в 19,24 кН. При перемещении 0,9 мм наблюдается критический рост повреждений в матрице, что приводит к изменению угла наклона кривой деформирования. Начало разрушения волокна происходит при перемещении 3 мм, что соответствует нагрузке в 86,7 кН. В дальнейшем из-за критического роста повреждений в волокне материал теряет свою несущую способность при перемещении 3,24 мм, что соответствует нагрузке в 89 кН. С момента образования первого повреждения до полной потери несущей способности материала нагрузка выросла в 4,63 раза. При проектировании конструкции предел прочности закладывается с коэффициентом 2, это означает, что даже при эксплуатационных нагрузках возможно критическое накопление повреждений, изменение свойств материала и его разрушение. Для подтверждения полученного предела прочности и упругих свойств материала проведена серия испытаний образцов на растяжение в направлении основы в соответствии со стандартом ASTM D3039. Испытания проводились с максимальным растягивающим усилием 200 кН и скоростью нагружения 2 мм/мин, при комнатной температуре. Испытывались стандартные образцы размером 250 × 25 мм со стеклопластиковыми концевыми накладками. Для регистрации значения локальной деформации образца его центральная область была препарирована тензорезисторами фольгового типа с базой измерений 5 × 3 мм, с двух сторон для компенсации изгиба. Для идентификации характера разрушения и определения морфологии накапливаемых повреждений в процессе испытаний образцов на растяжение проводилась запись событий акустической эмиссии. Для проведения акустической эмиссии применялась система регистрации AMSY-6 фирмы Vallen. Система позволяет производить регистрацию и подсчет событий акустической эмиссии, измерение амплитуд и частот акустического сигнала, сопутствующего локальному разрушению в структуре материала, означающих зарождение или развитие дефекта. Было использовано 4 датчика акустической эмиссии, расположенных парно вблизи накладок. Такое количество датчиков и их расположение позволяет производить локацию источников акустической эмиссии. Испытания были выполнены согласно методике проведения эксперимента для исследования феномена разрушения в углепластиках со сложной структурой, предложенной в [20]. Методика состоит из трех основных этапов: 1) проведение испытаний образцов до разрушения при квазистатической нагрузке с регистрацией сигналов акустической эмиссии; 2) определение характерных пороговых значений развития разрушения на основании данных акустической эмиссии; 3) исследование характера разрушения на этих этапах с помощью рентгеноскопии. По результатам анализа графиков зависимости кумулятивной энергии акустической эмиссии, полученным из испытаний образцов до разрушения, выбраны пороговые значения начала зарождения дефектов ε0, границы перехода от низкоэнергетических событий акустической эмиссии к высокоэнергетическим ε1 и обширного накопления повреждений ε2. Характерный график зависимости перемещений от нагрузки с наложенным на него графиком акустической эмиссии в логарифмической шкале приведен на рис. 4. Рис. 4. График зависимости нагрузки, энергии отдельных событий акустической эмиссии и кумулятивной кривой акустической эмиссии от перемещений траверсы для образца в направлении основы Чтобы определить характер разрушения, проведено томографическое исследование дефектов, возникших в образцах, испытанных на растяжение до пороговых значений ε1 и ε2. Для более отчетливого проявления дефектов образцы перед проведением томографии погружались в проникающую жидкость. На рис. 5 представлена томография образца перед началом испытаний в плоскости сечения XZ. Рис. 5. Томография образца перед началом испытаний в плоскости сечения XZ На рис. 6 представлена томография образца в плоскости сечения XZ, испытанного до порогового значения ε1, которое соответствовало нагрузке в 30 кН, перемещению 0,75 мм и значению акустической эмиссии 7,0e8 Eu. Также приведен рис. 6, б с разрушенными элементами при соответствующей нагрузке в двух сечениях плоскости XZ. Как видно, трещины зарождаются между волокнами утка, а также в местах переплетения волокон утка и основы. На рис. 7 представлена томография образца в плоскости сечения XZ, испытанного до порогового значения ε2, при нагрузке в 40 кН, перемещении 1,1 мм и соответствует значению акустической эмиссии 9.0e9 Eu. Также приведен рис. 7, б с разрушенными элементами при соответствующей нагрузке в двух сечениях плоскости XZ. В дальнейшем трещины продолжают свой рост в направлении Z. На рис. 8 представлена фотография образца в плоскости XY, испытанного до разрушения по первому этапу. а X Z Разршенные элементы б Рис. 6. Томография образца, испытанного до порогового значения ε1, в плоскости сечения XZ: а - эксперимент; б - расчет а X Z азрушенные элементы б Рис. 7. Томография образца, испытанного до порогового значения ε2, в плоскости сечения XZ: а - эксперимент; б - расчет Рис. 8. Фотография образца в плоскости XY после разрушения На снимках, полученных по результатам томографии, видны поперечные наклонные трещины, равномерно распределённые по всей длине образца. При этом в зоне накладок наблюдается незначительное увеличение плотности трещин, обусловленное геометрическим концентратором напряжений. У образцов, испытанных до уровня энергии ε1, наблюдаются зарождение характерных поперечных трещин на границе между нитями основы и утка, диагональное направление которых обусловлено извитостью нитей утка. На образцах, испытанных до уровня энергии ε2, трещины объединяются в магистральные поперечные. Плотность трещин в зависимости от нагрузки возрастает. Одной из особенностей испытанного материала является отсутствие продольных трещин, развивающихся в направлении прикладываемой нагрузки, что может являться следствием переплетения жгутов основы жгутами утка в направлении, перпендикуляром плоскости образца. Графики зависимости силы реакции от перемещений траверсы образцов и расчётный график были представлены на рис. 3. После проведения серии расчетов выявилось, что большой вклад в общую жёсткость при проведении испытаний вносят технологические накладки на образце, вследствие чего жёсткость накладок была учтена в расчете. Для этого в расчетную модель был добавлен упругий элемент с рассчитанной жесткостью траверсы. Средняя разрушающая нагрузка по испытаниям трех образцов составила 111 кН. Расчетная разрушающая нагрузка ниже экспериментальной на 18,9 % и составляет 89 кН. Перемещения при разрушающей нагрузке составили 3,24 мм, что на 2,56 % ниже, чем средние перемещения из эксперимента (3,325 мм). Сопоставление результатов акустической эмиссии с расчетом представлено на графике зависимости количества разрушенных элементов от перемещений (рис. 8). Черным маркером обозначено количество разрушенных элементов за один шаг нагружения, а красным сумма всех повреждений за все предыдущие шаги нагружения. Такая зависимость должна быть схожа по характеру с кумулятивной кривой, полученной акустической эмиссией. Из представленных на рис. 3 и 9 графиков видно, что после начала разрушения матрицы начинается бурный рост повреждений, вследствие чего кривая деформации меняет свой наклон, но в дальнейшем рост повреждений замедляется. В расчете, в отличие от эксперимента, рост повреждений после растрескивания матрицы не спадает, а, немного понизившись, продолжается вплоть до начала разрушения волокон. Рис. 9. График зависимости количества разрушенных элементов от перемещений 3. Испытание на сдвиг Также проведена серия испытаний образцов на сдвиг в плоскости XY в соответствии со стандартом ASTM D5379. Испытания проводились с использованием оптически измерительного комплекса Vic-3D. Фото образца с нанесенным покрытием для испытаний с Vic-3D и схематичным расположением тензодатчиков показано на рис. 10. Тензодатчики расположены под углом 90 градусов относительно друг друга и под углом 45 градусов относительно оси x. х y Рис. 10. Фото обраца с нанесенным покрытием x' y' Поле сдвиговых деформаций из расчета образца с полученными упругими эффективными свойствами, а также картина распределения сдвиговых деформаций, полученная с помощью Vic-3D, представлены на рис. 11. Как видно из рисунка, в центре образца образуется полоса равномерных сдвиговых деформаций. а б Рис. 11. Поле деформаций, полученное в расчете и с помощь Vic-3D: а - расчет; б - эксперимент После окончания испытаний на сдвиг из образца был вырезан участок с размерами 20 × 10 мм. Участок вырезан из зоны, в которой образовывались равномерные сдвиговые напряжения, после чего была выполнена его томография. Томография образца в плоскости сечения YZ представлена на рис. 12. Y Z Рис. 12. Томография в плоскости сечения YZ На томографии видны поперечные и наклонные трещины, проходящие по всей толщине образца. На рис. 13 показана картина распределения интенсивности напряжений по Мизесу в плоскости YZ, полученная из упругого расчета, а также картина с разрушенными элементами после потери несущей способности материала. Видно, что напряжения в матрице концентрируются между рядами волокон, направленных по основе. Вследствие этого данные концентраторы приводят к поперечному растрескиванию матрицы, что подтверждается проведённой томографией. Рис. 13. Напряжения по Мизесу в матрице в плоскости сечения YZ (а), картина разрушенных элементов в сечении YZ (б) Анализ полей перемещений, полученных с помощью Vic-3D и графика зависимости силы реакции от перемещений траверсы, показал, что оснастка, применяемая в эксперименте, имеет относительно малую жёсткость. В соответствии со стандартом ASTM D5379 была построена зависимость напряжений от деформаций. Напряжения вычислялись по формуле где P - сила реакции; S - площадь поперечного сечения. Деформация вычислялась как сумма деформаций двух тензодатчиков. Тензодатчики, приклеенные к поверхности образца, отклеились при определенном уровне деформации. Поэтому кривую деформации удалось построить не до полного разрушения образца. На рис. 14 показаны расчетные и экспериментальные кривые зависимости напряжений от деформаций. Усредненные по четырем образцам предельные напряжения показаны горизонтальной линией. При построении расчетной кривой брались усредненные сдвиговые напряжения по всему объему элемента периодичности и усредненные деформации , взятые с поверхности образца, перпендикулярной оси z в системе координат, , повернутой на 45 градусов вокруг оси z (см. рис. 10). Максимальные сдвиговые напряжения в эксперименте составляют 9,28 кгс/мм2, максимальные сдвиговые напряжения из расчёта составляют 9,18 кгс/мм2. Различие в напряжениях составляет 1 %. Данные расхождения следует считать удовлетворительными. Упругие участки кривых деформирования полностью совпадают, что подтверждает правильность рассчитанных упругих свойств. Выявлено значительное влияние оснастки на получаемые кривые зависимости силы реакции от перемещений. Данное влияние учтено при расчете. Рис. 14. Зависимость напряжений от деформаций Для эксперимента на сдвиг задавались следующие кинематические граничные условия, реализующие чистый сдвиг: В дальнейшем, для более точного моделирования процесса разрушения необходимо различать разрушение ортотропного и изотропного материала, для ортотропного применять более сложные критерии прочности, а для изотропного материала - иную модель деградации. Также необходимо ставить в зависимость коэффициенты деградации материала от объема разрушаемого элемента, так как вклад в общую жёсткость модели у элементов различен, и чем больше элемент, тем меньше должен быть коэффициент деградации материала для него. Выводы Проведены численные эксперименты с представительным объемом идеализированного виртуального прототипа трехмерно-армированного композиционного материала и определены недостающие эффективные упругие свойства и пределы прочности. Проведена серия испытаний образцов на растяжения по основе в соответствии со стандартом ASTM D3039, выполнена акустическая эмиссия и томография образцов на различных стадиях нагружения. Проведена серия испытаний образцов на сдвиг (в соответствии со стандартом ASTM D5379) с использованием оптического измерительного комплекса Vic-3D. Выполнена томография образца после испытания. Полученные данные сопоставлены с кривыми деформирования, томографией и акустической эмиссией. Расхождение между предельной расчетной нагрузкой и экспериментальной в случае растяжения по основе составило 18,9 %. При этом перемещения при разрушающей нагрузке составили 3,24 мм, что на 2,56 % ниже, чем средние перемещения из эксперимента (3,325 мм). В случае сдвига расхождение между предельными расчетными и экспериментальными сдвиговыми напряжениями составило 1 %.

About the authors

M V Tsepennikov

Perm National Research Polytechnic University

A A Strom

Perm National Research Polytechnic University

I A Povyshev

Perm National Research Polytechnic University

O Yu Smetannikov

Perm National Research Polytechnic University

References

Statistics

Views

Abstract - 196

PDF (Russian) - 95