Full Text

Введение Моделирование процесса разрушения композиционных материалов, представляющих собой совокупность двух тел, соединенных посредством адгезива конечной толщины, определяется типом нарушения связи между компонентами. Рассматриваются модели адгезионного разрушения в случае, когда адгезив целиком отделяется от материала, или когезионного разрушения по массиву адгезива или связанных им тел. Для адгезионных моделей [1-4], как правило, пренебрегают толщиной адгезива, а его механические свойства сводятся к силам взаимодействия склеенных материалов, которые могут иметь разные механические свойства [5-7]. При когезионном разрушении [8-13] рассматривается процесс зарождения трещины в материале с конкретными механическими свойствами. Модель расслоения композита с одинаковыми физическими свойствами, но с разными прочностными характеристиками рассмотрена в работе [14]. В этом плане представляется перспективной разработка таких моделей, которые бы учитывали как механические свойства материалов композита, так и тип разрушения в зависимости от напряженно-деформированного состояния адгезива и соединенных им тел. В настоящее время наиболее распространенными являются подходы на основе метода конечных элементов [15, 16]. Недостатком данных подходов при прямом моделировании геометрии связующего компонента является наличие точек сингулярности различных видов [17, 18]. Принимая ту или иную форму окончания связующего слоя, в решении задачи приходим к существенно различным распределениям напряжений, в том числе и сингулярным. Отметим, что в случае малости толщины связующего геометрия его окончания в зоне обрыва соединений является неопределенной и связь соответствующей области с заданной геометрией является не вполне корректной. Таким образом, разработка моделей, в которых напряженно-деформированное состояние (НДС) соответствующего композита не зависит от геометрии окончания адгезионного слоя (АС), является достаточно актуальной. Считаем, что в предлагаемой модели адгезионный слой и соединяемые им материалы могут иметь разные механические свойства. Толщина адгезива мала по сравнению с соответствующими характеристиками соединяемых тел. Для описания НДС адгезионного слоя применим концепцию слоя взаимодействия, используемую в модели трещиноподобного дефекта [19, 20]. Средние напряжения в слое выражаются через компоненты тензора напряжений на границах АС. Соответствующие граничные напряжения естественным образом формируют граничные условия для сопрягаемых со слоем материалов. Граничные напряжения в данном случае будем ассоциировать с адгезионными силами сцепления и достижение ими критического значения трактовать как адгезионное разрушение. Использование средних по толщине напряжений позволяет отказаться от конкретизации геометрии торца АС. Таким образом, задача определения предельного состояния композиционного материала под воздействием внешней нагрузки сводится к нахождению предела прочности в соответствующих материалах и по границам слоя. По достижении той или иной критериальной характеристики разрушение трактуется либо как адгезионное, либо как когезионное. 1. Постановка задачи На рис. 1 представлено тело, состоящее из трех областей, в общем случае с различными материальными свойствами, где область 3 ассоциируется с клеевой подложкой толщиной , размер которой мал по сравнению с толщинами тел 1 и 2. Процесс нагружения предполагаем квазистатическим и изотермическим. Условие равновесия запишем в вариационной форме: (1) где - внешняя нагрузка на контуре ; - тензор напряжений; - тензор деформаций; - поле перемещений. Так как система из трех тел находится в равновесии, запишем уравнения равновесия для каждого из них: (2) (3) (4) где , - векторы напряжений, действующие по верхней и нижней границе АС; - граничные напряжения АС; , - векторы перемещения верхней и нижней границы АС; - внешние контуры тел 1 и 2 без учета границ с АС. Торцы слоя считаем свободными от напряжений. При этом постулируется жесткое сцепление между границами АС c областями 1, 2. Кроме того, принимаем, что векторы напряжений на сопряженных границах АС равны и противоположны векторам напряжений сопряженных границ тела. 2 1 3 0 B A C D F X2 X1 M N S T R Рис. 1. Нагружение составного тела Работу внутренних напряжений в АС выразим посредством средних характеристик НДС, полагая (5) где , - соответственно тензоры средних напряжений и деформаций в слое, Отметим, что в силу симметрии касательных напряжений средние касательные напряжения также симметричны: . Средние деформации и перемещения определяем через их граничные значения следующим образом [19, 20]: (6) (7) (8) Из выражений (7) получим представление средней сдвиговой деформации вдоль слоя: (9) Из (4) и (5) с учетом (6) и (9) получаем вариационное условие равновесия АС: (10) При отсутствии торцевых нагрузок на АС, интегрируя по частям, имеем (11) (12) где . Из (10), с учетом (11) и (12), приравнивая слагаемые при одинаковых вариациях, приходим к представлениям граничных напряжений АС: , (13) , (14) Подставив (13) в (2), а (14) в (3) и выполнив преобразования, обратные (11), (12), приходим к вариационным уравнениям равновесия для тела 1: (15) и тела 2: (16) Отметим, что соотношения (15), (16) получены без ограничения на свойства материалов и форму тела. Уравнения (15), (16) необходимо замкнуть конкретными определяющими соотношениями. Поведение материала при нагружении определяем законом Гука для тела 1: (17) (18) для тела 2: (19) (20) где - девиатор тензора напряжений; - девиатор тензора деформаций; - модуль объемного сжатия; - модуль сдвига ; В материале слоя определяющие соотношения считаем справедливыми для средних по толщине слоя характеристик НДС. (21) (22) В результате подстановки в определяющие соотношения (21) и (22) выражений компонент средних деформаций (6), (9) средние напряжения определяются через граничные перемещения и их производные. В этом случае в уравнении (15) будут присутствовать перемещения границы NS тела 2 (см. рис. 1), а в уравнении (16) - перемещения границы FC тела 1 (см. рис. 1). Таким образом, совместное решение системы (15)-(22) сводится к определению поля перемещений в телах 1 и 2. После определения полей перемещений в телах 1 и 2, в том числе и по границам со слоем, из (21), (22) находим средние напряжения в слое. Зная распределение средних напряжений вдоль слоя, из уравнений (13), (14) находим граничные напряжения по его границе. 2. Метод дискретного решения Еще раз отметим, что прямое моделирование композитных материалов с явным заданием геометрии торцевой поверхности АС сопряжено с зависимостью НДС в окрестности торцевой поверхности от условий сопряжения слоя и смежных материалов. В частности, если торцевая поверхность является плоскостью, то модель будет сингулярной. Предлагаемый вариант рассмотрения модели композита за счет введения средних по толщине характеристик НДС универсален в плане независимости НДС от геометрии окончания АС. Для численного решения системы (15)-(22) можно использовать различные дискретные методы, обеспечивающие численную сходимость. В частности, будем использовать МКЭ для дискретизации задачи. Существенно, что построение решения формально не накладывает ограничений на размер конечного элемента, который может быть меньше параметра для обеспечения численной сходимости решения. Однако условия разрушения АС формулируются для δ-элементов слоя размером и на их верхней и нижней границах с использованием процедуры осреднения характеристик НДС, полученных в результате решения системы (15)-(22) и (13), (14): (23) Это является следствием основного физического допущения - разрушение охватывает частицу материала с характерным размером [21, 22]. В качестве критерия разрушения АС будем использовать критерий Кулона, согласно которому разрушение δ-элемента происходит при достижении максимальным главным растягивающим усредненным напряжением критического значения: . Для адгезионного разрушения будем сравнивать граничные напряжения , с адгезионной прочностью на отрыв, а граничные напряжения , - с адгезионной прочностью на сдвиг. Решение задачи (15)-(22) дает распределение поля перемещений в узловых точках тела 1, в том числе и по границе с АС. После нахождения соответствующего решения возникает задача определения НДС в АС. При решении (15)-(22) используем метод конечного элемента с квадратичным распределением поля перемещений. 3. Результаты численного расчета В качестве примера использования предлагаемой модели рассмотрим композиционный материал (см. рис. 1) в состоянии плоской деформации, состоящий из адгезионного слоя, близкого по характеристикам к эпоксидной смоле с механическими характеристиками: Па; ; Па - предел прочности; материал 1 и 2 соответствует характеристикам алюминиевого сплава Д16Т: Па; Па. Адгезионную прочность по границе смола-сплав принимаем равной Па на отрыв и сдвиг. Геометрические характеристики образца выбираем следующими: м; м; м; м; м. На границе задаем распределенную нагрузку Па, направленную под углом к оси . Граница жестко закреплена от перемещений. Остальная поверхность свободна. В качестве исследуемых параметров задачи будем рассматривать граничные и средние напряжения в адгезионном слое. На рис. 2 показана вычислительная сходимость решения для отношения максимального главного растягивающего усредненного напряжения на δ-элементе АС к соответствующему значению при грани конечного элемента равной параметру в зависимости от размера в концевой области слоя относительно параметра . Как видно из графика (см. рис. 2), результат расчета имеет вычислительную сходимость. Дальнейшие вычисления проводились при характерном размере конечного элемента . На рис. 3 показаны распределения напряжений (23) по структурным δ-элементам слоя. Все напряжения отнесены к граничному напряжению на первом элементе . Графики рис. 3, a определяют граничные напряжения, где кривая 1 задает напряжение , кривая 2 - , кривая 3 - , кривая 4 - . На рис. 3, б приведены средние напряжения, где кривая 1 задает напряжение , кривая 2 - , кривая 3 - . Рис. 2. Сходимость решения 1 2 3 4 3 2 1 а б Рис. 3. Распределение напряжений в слое В результате решения задачи на первом δ-элементе АС были получены следующие характеристики: 6 Па, 6,3 Па, -1,7 Па, -2,1 Па, 6,9 Па. Следовательно, разрушение соответствующего композита произойдет при нагрузке Па посредством адгезионного отслоения по нижней границе слоя в силу достижения предела адгезионной прочности на отрыв. Когезионное разрушение в соответствующих композитах будет наблюдаться в случае, когда адгезионная прочность на отрыв или сдвиг будет приближаться к пределу прочности АС. Заключение Использование предлагаемого подхода позволяет естественным образом учитывать характер разрушения при механических воздействиях на композиты с адгезионным слоем. В этом случае задача сводится к решению вариационного уравнения в рамках тех или иных определяющих соотношений.

About the authors

V V Glagolev

Tula State University

A A Markin

Tula State University

A A Fursaev

Tula State University

References

Statistics

Views

Abstract - 210

PDF (Russian) - 104