Full Text

Введение Моделирование образования новых материальных поверхностей предполагает описание двух стадий процесса нагружения тела: докритической и посткритической. Докритическая стадия может быть описана в рамках известных уравнений механики деформируемого твердого тела как без учета нелинейных эффектов, так и в нелинейной (упругопластической) постановке. При этом основная проблема состоит в формулировке критериев перехода к посткритической стадии. Общая постановка задачи посткритического поведения при упругопластическом деформировании далека от завершения. Это связано, во-первых, с выбором условий прекращения взаимодействия между поверхностями, во-вторых, с необходимостью формулировки граничных условий на вновь образуемых поверхностях. Если разрушение локализуется на материальной поверхности меры ноль, то имеет место классическая модель математического разреза [1-3]. Если взаимодействие прекращается в результате разрушения слоя конечной толщины между поверхностями, то применяется концепция физического разреза [4, 5]. Наиболее изученными и обоснованными в настоящее время являются модели докритического деформирования. Так как в этом случае не происходит образования новых материальных поверхностей, то любой подход, согласующийся с экспериментом по определению предельной нагрузки, может быть рассмотрен в качестве модели предразрушения. Соответствующие модели применяются как для оценки трещиностойкости конструкций из однородных тел, биметаллов [6], так и из композиционных материалов, где адгезионные слои рассматриваются как область предполагаемого разрушения [7, 8]. Выделим модели, в которых явно присутствует параметр структуры среды [5, 8-12]. При выполнении того или иного критерия разрушения наступает фаза образования новых материальных поверхностей. Принципиальным моментом данного этапа является дискретность приращения поверхностей разрыва [13-15] и локализация процесса разрушения. При моделировании разрушения в пределах материальной поверхности меры ноль используются когезионные элементы [16-19] конечно-элементного пакета [20] или контактные элементы системы [21]. В работах [22-24] рассмотрен широко используемый билинейный закон распределения взаимодействия в когезионной зоне от раскрытия трещины с ниспадающим участком. Основным недостатком данного подхода является то, что когезионные элементы должны располагаться на траектории разрушения, которая зависит от новых материальных поверхностей. Кроме того, различные законы когезионного взаимодействия и материальные характеристики когезионных элементов существенно влияют на распределение напряжено-деформированного состояния [25] и требуют экспериментального подтверждения. Сложность описания разрушения материального объема при использовании падающего участка диаграммы растяжения [26, 28] связана с построением определяющих соотношений неустойчивого по Друкеру деформирования и их подтверждением в экспериментах. Основным подходом при конечно-элементном моделировании процесса разрушения является процедура kill element, подробно описанная в [21]. В этом случае по достижении в конечном элементе критериальной характеристики разрушения выделяемый материальный объем исключается из рассмотрения путем умножения локальной матрицы жесткости элемента на число, близкое к нулю. Отметим, что данная процедура является корректной при упругом деформировании, когда нагрузка и разгрузка определяются одинаковыми модулями в определяющих соотношениях. В случае упругопластического деформирования тела с удаляемым таким образом элементом необходимо учитывать возможную разгрузку. В данной работе предлагается рассматривать процесс разрушения δ-элемента как локальную разгрузку от напряжений, действующих на тело со стороны δ-элемента, при неизменной внешней нагрузке. Дискретизация задачи методом конечных элементов предполагает описание взаимодействия между частицами тела посредством узловых сил [29, 30]. Считаем, что при разрушении материального объема, взаимодействующего с несколькими конечными элементами, узловые силы взаимодействия будут разгружаться до нуля в рамках простого процесса при неизменной внешней нагрузке, достигнутой к моменту разрушения. Таким образом, неустойчивое по Друкеру деформирование исключается из рассмотрения. Отметим, что определение узловых сил взаимодействия с удаляемым материальным объемом во внутренних точках деформируемого тела является отдельной задачей. В статье предложена процедура повторного нагружения тела без разрушаемого объема, но с найденным по его границе полем перемещений для нахождения соответствующих узловых сил. 1. Постановка задачи На рис. 1 представлено тело, состоящее из трех частей с различными материальными свойствами, где область 3 ассоциируется с клеевой подложкой толщиной δ0, размер которой мал по сравнению с толщинами тел 1 и 2. Процесс нагружения предполагаем Рис. 1. Нагружение составного тела Fig. 1. The loading of the compound body квазистатическим и изотермическим. Условие равновесного протекания процесса [31] при малых деформациях и поворотах материальных волокон принимаем в виде (1) где - поле скоростей; - тензор деформации скорости; ; - скорость внешней нагрузки на контуре ; - скорость тензора напряжений; - параметр подобный времени; - внутренняя область составного тела. В области 3 средние по толщине δ0 скорости напряжений определяем в виде а средние скорости и скорости деформаций - через их граничные значения: (2) (3) (4) где , - векторы скоростей верхней и нижней границы области 3. Из выражений (3) приходим к представлению средней скорости сдвиговой деформации вдоль рассматриваемой области: (5) Принимаем, что векторы скоростей напряжений на сопряженных границах слоя 3 равны и противоположны векторам скоростей напряжений сопряженных границ тела. Кроме того, постулируется жесткое сцепление между границами области 3 и областями 1, 2: (6) Рассматривая скорости граничных векторов напряжений слоя 3 в качестве граничных условий для смежных с ним областей [5, 8], приходим к совместному решению вариационных уравнений равновесия для тела 1: (7) и тела 2: (8) Уравнения (7), (8) необходимо замкнуть конкретными определяющими соотношениями, связывающими скорости напряжений со скоростями деформаций. Поведение материала тел 1 и 2 при активном нагружении () определяем следующими физическими соотношениями: (9) (10) где ; - скорость девиатора тензора напряжений; - девиаторная составляющая тензора деформаций скоростей; ; - модуль объемного сжатия; - модуль сдвига при , при ; - предел текучести соответствующего материала; - интенсивность касательных напряжений. В состоянии разгрузки () определяющие соотношения запишем в виде , (11) . (12) В материале слоя определяющие соотношения считаем справедливыми для средних по толщине слоя скоростных характеристик НДС: (13) (14) (15) где - модуль объемного сжатия материала слоя; - модуль сдвига материала слоя при , при ; - предел текучести материала слоя. В результате подстановки в определяющие соотношения (13)-(15) выражений компонент средних деформаций скорости (2), (5) средние скорости напряжений определяются через граничные скорости и их производные. Таким образом, решение системы (7)-(12) сводится к определению поля скоростей в телах 1 и 2 (см. рис. 1). В этом случае в уравнении (7) будут присутствовать скорости границы NS тела 2 (см. рис. 1), а в уравнении (8) - скорости границы FC тела 1 (см. рис. 1). После определения поля скоростей для нахождения поля перемещений и деформаций используем эволюционные соотношения (16) (17) где - поле перемещений; - тензор деформации. С учетом определяющих соотношений (9)-(12) находим поле напряжений в телах 1 и 2: (18) а из (13)-(15), используя значения граничных скоростей слоя и связи (2), (5), определяем поле средних напряжений слоя: (19) Условия разрушения адгезионного слоя (АС) формулируем для δ-элементов слоя размером Это является следствием основного физического допущения - разрушение охватывает частицу материала с характерным размером [8, 11-15]. В качестве критерия когезионного разрушения АС будем использовать критерий Кулона, согласно которому разрушение δ-элемента происходит при достижении максимальным главным растягивающим усредненным напряжением критического значения: , где (20) а получены в результате решения (7)-(12), (19). 2. Когезионное дискретное разрушение Решение задачи докритического деформирования (7)-(15) строим на основе метода конечного элемента с квадратичным распределением поля скоростей на элементе. В этом случае поля скоростей напряжений и деформаций распределены по линейному закону в пределах каждого конечного элемента. Отметим, что размер грани конечного элемента может быть существенно меньше линейного размера и определяется вычислительной сходимостью дискретного решения задачи [8]. Интегрируя по площади конечного элемента компоненты соответствующих тензоров, определяем средние по элементу значения по которым, используя эволюционные соотношения (17), (18), находим напряжения и деформации . Переход в состояние пластичности соотносим со средними характеристиками НДС. Решение упругопластической задачи строим на основе метода «упругих решений» [32]. Секущий модуль сдвига каждого конечного элемента будем определять исходя из средних компонент считая секущий модуль в пределах конечного элемента постоянным. Узел элемента характеризуется узловой скоростью и скоростью узловой силы . Используя эволюционные соотношения (21) (22) в узлах можно определить узловые перемещения и узловые силы. При когезионном разрушении постулируем, что элемент слоя полностью разрушается в течение интервала времени . В этом случае на момент разрушения мысленно удалим из слоя элемент размером , где напряженное состояние, определяемое согласно (20), соответствует критериальной величине , а его влияние на тело компенсируем внешней узловой нагрузкой , со стороны элемента (рис. 2). Рис. 2. Эквивалентные нагрузки элемента слоя в момент предразрушения Fig. 2. The equivalent loadings of the element layer during the pre-fracture Полагаем, что в течение интервала δ-элемент будет полностью разрушен, если нагрузка, действующая со стороны элемента слоя на тело, станет равна нулю. В этом случае в течение интервала времени при неизменной внешней нагрузке необходимо разгрузить новые материальные поверхности. Будем считать процесс разгрузки простым, зависящим от одного параметра . Это соответствует заданию скоростей узловых внешних нагрузок: (23) Таким образом, рассматривая (23) в качестве внешней нагрузки, показанной на рис. 3, для конечно-элементной постановки задачи (7)-(15) можно смоделировать процесс формирования новых поверхностей на характерном размере. Процесс простой разгрузки на интервале может быть разбит на несколько интервалов, однако в рамках данной статьи ограничимся одним шагом. Рис. 3. Нагрузка в процессе разрушения элемента слоя Fig. 3. The loading during the fracture of the layer element Рис. 4. Схема повторного нагружения Fig. 4. The scheme of re-loading Основная проблема заключается в нахождении узловых сил , в момент предразрушения и, следовательно, в формировании граничных условий (23). Дискретное решение (7)-(15) для докритического деформирования по границам сопряжения слоя в качестве неизвестных содержит вектор узловых скоростей по которому из (16) может быть сформирован вектор перемещений , . Для определения узловых сил , предлагается использовать процедуру повторного нагружения, состоящую в замене действия разрушаемого элемента заданием процесса изменения со «временем» узловых перемещений элемента, при повторении закона внешнего воздействия на тело (рис. 4). Конечно-элементное решение повторного нагружения в узлах с заданным полем узловых скоростей , в качестве неизвестных будет содержать вектор скоростей узловых сил , по границам разрушаемого элемента. В зависимости от распределения основной внешней нагрузки в результате разрушения δ-элемента процесс разделения может продолжаться при неизменном значении основной внешней нагрузки (неустойчивое разрушение) либо прекратиться. При этом следует учитывать, что в результате δ-разгрузки, связанной с разрушением АС, в окружающих слоях может происходить догрузка с выходом в пластическую область. Данное обстоятельство не позволяет непосредственно использовать теорему о разгрузке А.А. Ильюшина [32]. 3. Результаты решения В качестве примера расчета рассмотрим композиционный материал в состоянии плоской деформации, состоящий из двух элементов с материальными характеристиками, близкими к сплаву Д16: Па; ; ; Па - предел упругости; Па - предел прочности, связанных эпоксидной смолой со следующими свойствами: Па; ; Па. Геометрические характеристики композита выбирались следующими: м; м; м; м; м. Скорость внешней нагрузки направлена по углом к оси . В качестве тестовой рассмотрим задачу разрушения элемента слоя в упругой постановке. Пусть критическое состояние в δ-элементе определяет внешняя нагрузка Вычислительная сходимость решения предполагает разбиение участка границы, сопряженного с δ-элементом, на четыре конечных элемента [8]. С учетом квадратичного закона распределения поля скоростей на конечном элементе девять узловых сил по верхней и нижней границе разрушаемого δ-элемента будут определять равновесие композита в момент предразрушения (см. рис. 2). Соответствующие значения узловых сил, определяемые методом повторного нагружения, представлены в таблице. Силы в таблице отнесены к минимальному значению модуля их проекций. Нумерация узлов выбрана слева направо. Значения проекций безразмерных узловых сил The values of the projections of the dimensionless nodal forces 8,9 3,2 1,3 2,3 1,1 2,1 1,0 2,0 1,0 -10,5 -9,8 -4,3 -8,3 3,9 -7,4 -3,5 -6,6 -3,3 7,5 -2,2 -1,5 -3,0 -1,5 -2,9 -1,4 -2,8 -1,4 -5,5 9,4 4,7 8,9 4,2 8,0 3,8 7,2 3,6 Без ограничения общности минимальный модуль проекций найденных узловых сил может быть рассмотрен в качестве параметра локальной разгрузки для определения граничных условий (23). График 2 на рис. 5, а показывает распределение безразмерных вертикальных перемещений по границе тела 1 и слоя 3 к концу докритического деформирования, график 3 - перемещений за счет процесса разрушения структурного элемента слоя, график 1 - суперпозицию полей перемещений. Перемещения отнесены к значению перемещения на начало процесса разрушения. Безразмерная координата в направлении оси абсцисс определяется как . График 1 на рис. 5, б определяет соответствующие вертикальные перемещения при решении задачи без первого структурного элемента слоя к концу интервала . Данное решение моделирует подход kill element, используемый в процедурах конечно-элементного моделирования [15]. График 2 на рис. 5, б повторяет кривую 1 на рис. 5, а. Как видно из рис. 5, б, совпадение графиков 1 и 2 показывает адекватность предлагаемой процедуры разгрузки новых поверхностей в процессе разрушения и подхода kill element в рамках упругого деформирования. Отметим, что по горизонтальным перемещениям тенденция аналогична. Максимальное главное растягивающее напряжение во втором элементе слоя при этом возрастает, что предполагает катастрофический характер разрушения в данной схеме. Проводя сравнение решения упругой задачи для заданной внешней нагрузки с первым структурным элементом слоя и без него, отметим, что для большинства конечных элементов интенсивность напряжений () уменьшается. На рис. 6 показана конечно-элементная область тел 1 и 2 (см. рис. 1), где интенсивность напряжений уменьшилась (элементы выделены заливкой) при решении задачи без первого структурного элемента слоя по сравнению с решением со структурным элементом. Данное обстоятельство дает право предположить, что при упругопластическом характере деформирования возможно образование зоны разгрузки пластических областей. Рис. 5. Распределения вертикальных перемещений: а - штрихпунктирная линия - докритическое деформирование, пунктирная линия - процесс разрушения структурного элемента, сплошная линия - суперпозиция решений; б - сплошная линия повторяет соответствующую линию (а), пунктирная линия описывает докритическое деформирование без первого структурного элемента Fig. 5. The distribution of vertical displacements: a - the broken line shows the subcritical deformation; the dashed line shows the fracture of the structural element; the solid line shows the superposition of the solutions; б - the solid line repeats the corresponding line (а); the dashed line shows the subcritical deformation without the first structural element l Рис. 6. Конфигурация конечно-элементной области, где уменьшается интенсивность напряжений при упругом решении задачи, построенная на основе подхода kill element. Элементы выделены цветом Fig. 6. The configuration of the finite element area where the stresses intensity decreases in the elastic solution of the problem based on the kill element approach. The elements are highlighted. Однако подход на основе процедуры kill element не позволяет корректно описать данный процесс. Рассмотрим решение данной задачи на основе предлагаемой процедуры. На рис. 7 выделена область пластических деформаций в состоянии предразрушения. Значение внешней нагрузки при этом обозначим через . Максимальное главное напряжение в адгезионном слое на характерном элементе равняется пределу прочности, а вне слоя предел прочности не достигается. Следовательно, разрушение в слое наступит быстрее, чем в сопряженных с ним материалах. Предполагаем, что адгезионная связь смола-сплав достаточно прочна и разрушение будет проходить по массиву адгезионного компонента. l Рис. 7. Зона пластичности в состоянии предразрушения. Область выделена цветом Fig. 7. The plasticity zone in the pre-fracture state. The area is highlighted Разрушение первого элемента приводит к перераспределению зоны пластичности и разгрузке ряда элементов композита. В этом случае рассмотрение разрушения как термомеханического процесса дает отличный от процедуры kill element результат. На рис. 8 темным показана зона пластической догрузки, а более светлым - элементы, где произошла упругая разгрузка из пластической области по завершении процесса разрушения первого δ-элемента. l - δ0 Рис. 8. Эволюция зоны пластичности и разгрузки в процессе разрушения первого элемента. Область разгрузки выделена серым цветом, а упругопластическая - черным Fig. 8. The evolution of the plasticity and unloading zones during the fracture of the first element. The unloading area is highlighted in gray; the elastoplastic area is in black Решая задачу нагружения композита критической нагрузкой без первого структурного элемента слоя, приходим к перераспределению пластической области (рис. 9). Сопоставляя рис.8 и рис. 9, видим, что разгрузка из пластической области не учитывается в подходе kill element, при этом область пластичности немного больше объединения областей разгрузки из пластической области и пластической догрузки, отображенных на рис. 8. l - δ0 Рис. 9. Зона пластичности в композите без первого структурного элемента адгезионного слоя. Область выделена цветом Fig. 9. Plasticity zone in the composite without the first structural element of the adhesive layer. The area is highlighted Разрушение первого элемента приводит к превышению предела прочности на втором элементе слоя, что означает его разрушение при фиксированной внешней нагрузке. Таким образом, при данном виде нагружения имеет место катастрофическое разрушение. Заключение В работе предложен вариант описания дискретного разрушения материального объема на основе гипотезы простой разгрузки. Для определения внутренних узловых сил, действующих на разрушаемый объём, использована процедура повторного нагружения с заданным полем узловых перемещений, найденным из решения задачи докритического деформирования. Результаты расчёта по предложенной модели в рамках упругого поведения материала не противоречат известной методике расчета kill element. В случае упругопластического материала напряженное состояние, получаемое при моделировании процесса разрушения методом локальной разгрузки, может существенно отличаться от состояния, определенного методом kill element. Предлагаемый подход позволяет учесть перераспределение пластических зон и возможность образования новых зон разрушения в результате локальной разгрузки. Предложенная методика может быть использована для моделирования процесса разрушения элементов в произвольном конечно-элементном континууме.

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 - 285

PDF (Russian) - 126