Computational simulation of the damage accumulation processes in cracked solids by the user procedure UMAT of Simulia Abaqus

Abstract


The paper presents the experience of using the user subroutine UMAT for FEM package SIMULIA Abaqus/CAE for damage accumulation processes in the vicinity of the crack. A continuum damage mechanics model based on the constitutive relations of linear elastic isotropic materials with the incorporated damage tensor components is used to describe the material behavior. The material nonlinearity arising from the deformation process is modelled by introducing an anisotropic damage tensor of the second rank into the constitutive equation. The material model is described by means of user procedure UMAT of SIMULIA Abaqus. The finite element (FE) mechanical constitutive model is implemented in Abaqus/Standard via a UMAT routine. Numerical experiments for a large series of cracked specimens have been performed. Computed stress and damage tensor components were found. It is shown that they are not dependent on the FE mesh refinement. Distributions of the damage tensor components in the vicinity of the crack tip in cracked specimens of different configurations under mixed mode loading in a wide range of mixed mode loadings are found. The configurations of active damage accumulation process zone in the cracked specimens are obtained. It is shown that the damage accumulation process has a substantial influence on the stress-strain state in the vicinity of the crack tip and leads to a decrease of the stress concentration in cracked specimens.

Full Text

Введение В настоящее время разрушение (макроскопическое нарушение сплошности тела в результате воздействия на него внешнего окружения) рассматривается с учетом процессов скрытого разрушения. Исследование скрытого разрушения (зарождение и развитие микродефектов, рассеянных по объему тела) осуществляется с помощью методов и теорий механики поврежденности - динамично развивающегося раздела современной механики деформируемого твердого тела, теоретические основы которой заложены Л.М. Качановым и Ю.Н. Работновым [1-3]. В последнее время большой интерес вызывают вопросы компьютерного моделирования деформирования, повреждаемости и разрушения неупругих материалов и конструкций [4-6] и развития численных методов, учитывающих повреждаемость материала конструкции с течением времени, а также описания взаимного влияния эволюции поврежденности и напряженно-деформированного состояния [7-11]. В [7-11] предложена и используется одна из моделей анизотропной поврежденности, экспериментальное подтверждение которой обсуждается в [12-16]. В целом в последнее десятилетие и в самое последнее время вопросам компьютерного моделирования роста трещин и длительного разрушения с использованием различных мер поврежденности посвящается большое количество работ [17-20], появились фундаментальные монографии [21-23], а в 2014 г. вышло в свет второе издание основополагающей книги Ю.Н. Работнова «Ползучесть элементов конструкций» [24]. Континуальная механика поврежденности исходит из того, что изменение со временем механических свойств материалов феноменологически можно интерпретировать как результат накопления повреждений, различных дефектов, микропор, трещин. Когда повреждения достигают опасного уровня, происходит разрушение. Трещинообразование начинается на самых ранних этапах деформации и связано с ростом имеющихся и возникновением новых суб- и микродефектов. В материале всегда имеется большое число различных дефектов, приводящих к высоким местным напряжениям. Уменьшение прочности (свойства тела сопротивляться воздействиям со стороны внешнего окружения) деформируемых твердых тел часто может быть объяснено скрытым разрушением и микродефектной структурой тела. Таким образом, поскольку повреждения тела существенно влияют на характер его разрушения, становится очевидным, что и механика разрушения, и механика поврежденности призваны решить главную прикладную задачу об оценке запаса прочности твердого тела. Поэтому сейчас одной из актуальной задач современной механики деформируемого твердого тела и механики разрушения, в частности, является компьютерное имитационное моделирование процессов нелинейного деформирования и разрушения твердых тел под нагрузкой с учетом процессов накопления рассеянных повреждений. Под поврежденностью понимается относительное сокращение вследствие распределенных внутри образца микродефектов эффективной, несущей нагрузку площади поперечного сечения образца [22]. Деградацию материала можно трактовать как последовательное уменьшение эффективной площади, несущей приложенную нагрузку. Процесс накопления повреждений в твердых телах может быть описан с помощью введений скалярной или тензорной меры поврежденности [1-3, 22]. В простейшем варианте поврежденность можно описать некоторым скаляром, структурным параметром [1]. В начальном состоянии при отсутствии поврежденности , с течением времени функция убывает. Функцию можно интерпретировать как сплошность. Ю.Н. Работнов [2, 3] ввел функцию , равную нулю в начальном состоянии и единице в момент разрушения, которую можно принять за меру охрупчивания. Функцию естественно назвать поврежденностью, и можно считать, что . Величина может быть интерпретирована как относительная площадь поперечного сечения, занятая трещинами и вакансиями. Предполагается, что скорость изменения параметра поврежденности зависит от напряжения и от поврежденности . Простейшая гипотеза состоит в том, что есть степенная функция отношения . Полученное напряжение носит название истинного напряжения. Это отношение может быть истолковано как среднее напряжение на площади поперечного сечения, свободной от трещин. Истинное напряжение вводится в определяющие уравнения, построенные для первой и второй стадий ползучести, чтобы описать третий участок кривой ползучести. Впоследствии концепция истинного напряжения была распространена на упругие и пластические деформации [4-6], что позволило учитывать накопления повреждений в условиях упругого нагружения и пластического течения [4-6, 7-9]. В рамках связанной постановки задачи скалярный параметр поврежденности или компоненты тензора поврежденности входят в определяющие соотношения материала [22]. Исторически первой явилась математическая модель Качанова-Работнова [1-3] (1) В соотношениях (1) - компоненты тензора скоростей деформаций ползучести; - девиатор напряжений; - параметр сплошности; - интенсивность касательных напряжений; - константы материала. После пионерских работ Л.М. Качанова и Ю.Н. Работнова [1-3] континуальная механика поврежденности двигалась по пути 1) усложнения определяющих уравнений [4-6], 2) отражения анизотропии, вызванной накоплением повреждений [7-16, 22], 3) решения краевых задач в связках упругость-поврежденность, пластичность-поврежденность и ползучесть-поврежденность [25-31]. Существующие модели поврежденности для изотропных материалов были обобщены для анизотропных линейно-упругих материалов (бетон, горные породы), для функционально-градиентных материалов с переменными по координатам свойствами [32], для композиционных материалов [33-42]. Например, позднее определяющие уравнения модели Качанова-Работнова (1) были обобщены и представлены в форме [43, 44] (2) В целом закон роста поврежденности при сложном напряженном состоянии можно установить с помощью экспериментов, в которых находятся комбинации напряжений, приводящие к одному и тому же времени до разрушения. Геометрически такие результаты представляются в виде поверхностей в пространстве напряжений [45]. Математически изохронные поверхности напряжений описывают с помощью некоторого скалярного инварианта напряжений , где - максимальное главное напряжение; - интенсивность напряжений; - гидростатическое напряжение. В настоящее время изохронную поверхность чаще всего описывают с помощью следующей линейной интерполяции: Решение краевых задач механики деформируемого твердого тела для элементов конструкций, находящихся в реальных эксплуатационных условиях, требует использования многофункциональных комплексов, таких как SIMULIA Abaqus. Однако среди моделей материалов, представленных в конечно-элементном пакете SIMULIA Abaqus, нет возможности прямого учета процесса накопления повреждений. Возможность моделирования сложных сред в пакете SIMULIA Abaqus реализуется с помощью пользовательских процедур, таких как UMAT и VUMAT [17, 44, 46]. В настоящей работе представлен опыт использования процедуры UMAT для описания процессов накопления повреждений в образцах с концентраторами напряжений в материале с определяющими уравнениями, базирующимися на конституциональных уравнениях линейно-упругого изотропного материала и включающими компоненты тензора поврежденности второго ранга [7-11]. В работе найдены распределения поврежденности у вершины трещины и проанализированы зоны активного накопления повреждений в условиях смешанного нагружения на примере целого ряда задач для тел с различными концентраторами напряжений и разрезами в полном диапазоне смешанных форм нагружения. Выполнена серия вычислительных экспериментов для диска с вертикальной наклонной трещиной (с различными углами наклона надреза к хорде диска) в материале с учетом процессов накопления повреждений. 1. Определяющие уравнения материала с включенными компонентами тензора поврежденности второго ранга Для учета эффектов анизотропии был предложен целый ряд математических моделей [7-16, 21, 22]. Одна из них - модель для боросиликатного стекла, основанная на введении тензорной меры поврежденности [7-16]. Определяющие уравнения материала в общем случае имеют вид (3) где - компоненты тензора напряжений; - компоненты тензора деформаций; - компоненты тензора упругих модулей, в общем случае зависящих от температуры; - компоненты тензора упругих модулей, обусловленные накоплением повреждений; - температурные деформации. Компоненты тензора упругих модулей вычисляются как классические компоненты тензора упругих модулей линейно-упругого изотропного материала (4) где - постоянные Ламе. Компоненты тензора упругих модулей , обусловленные накоплением повреждений и учитывающие данный процесс, вычисляются с помощью соотношений [7-11] (5) где - компоненты тензора поврежденности; - константы материала, определяемые экспериментально. Предполагается, что Оригинальная модель анизотропной поврежденности (3)-(5) была предложена для боросиликатного стекла [10] и развита в [11]. Модель анизотропной поврежденности (3)-(5), как показано в [8], хорошо описывает растрескивание стекла при вдавливании штампа, что подтверждается экспериментальными данными, приведенными в [13, 14]. Подробное изложение экспериментальной процедуры исследования процессов растрескивания и локализации поврежденности приведено в [8]. Уравнения (3)-(5) изначально предложены для описания процессов деформирования вязкоупругого тела, переходящего в стеклообразное состояние при понижении температуры. В [8-10, 12-14] показано, что модель может быть обобщена на случай произвольного стекла путем выбора параметров модели. Главным преимуществом используемой модели является ее широкое экспериментальное подтверждение для различных классов материалов и условий разрушения [14-16]. Также одним из преимуществ данной модели является ее простота и возможность использования в расчетном комплексе SIMULIA Abaqus, в отличие от моделей, представленных уравнениями (1) и (2), для которых требуется написание достаточно сложной процедуры интегрирования кинетического уравнения по времени. В силу указанной причины находят широкое распространение модели поврежденности с алгебраическими уравнениями, описывающими поврежденность [38, 47]. В целях упрощения математической модели (3)-(5) можно предположить, что недиагональные компоненты тензора поврежденности равны нулю, а нормальные компоненты определяются выражениями (6) где - главные напряжения; - пороговое и критическое значения напряжений соответственно. При напряжениях ниже порогового значения накопления повреждений не происходит. Если главное напряжение достигает критического значения , то материал считается полностью поврежденным, Если поврежденность вызывается касательными напряжениями в плоскости, то компоненты тензора поврежденности вычисляются с помощью соотношений [7, 9, 10] (7) В развернутой форме определяющие уравнения материала (3)-(5) принимают вид (8) Определяющие соотношения (2)-(8), учитывающие анизотропный характер поврежденности, соответствуют подходу, предложенному для описания процессов накопления поврежденности и разрушения композиционных материалов в [35]. Подробные обзоры современной научной литературы, посвященной проблематике описания процессов накопления поврежденности и разрушения, могут быть найдены в [36-38]. Определение значений материальных констант и в уравнениях (3)-(8) представляет собой нетривиальную задачу, требующую выполнения целого ряда экспериментов. Тем не менее для простых случаев можно воспользоваться следующими соображениями. В случае двухосного растяжения квадратной пластины можно принять, что компоненты тензора напряжений и обращаются в нуль, когда в то время как для плоской задачи справедливо Для данного нагружения в силу симметрии должно выполняться : (9) (10) Из полученных равенств можно найти значения материальных констант и . В вычислениях, приведенных ниже, используется именно этот выбор констант. В [9] указывается и альтернативный выбор констант , , используемый в оригинальной работе [10], в которой два параметра модели находятся из условия обращения в нуль компоненты , когда для одноосного растяжения плоского образца в направлении оси . В [7] тензорная мера поврежденности применяется для описания разрушения боросиликатного стекла в условиях трехосного растяжения. Согласно уравнениям (6) компоненты тензора поврежденности связываются с главными напряжениями и с направлениями главных осей тензора напряжений. Авторы апеллируют к экспериментальным данным [8] и замечают, что в условиях температурного нагружения, рассматриваемого ими, нагружение является строго пропорциональным и в центральной части цилиндрического образца реализуется трехосное растяжение с главными осями, совпадающими с осями цилиндрической системы координат. Поэтому определяющие уравнения формулируются в цилиндрической системе координат: (11) где приняты обозначения Постоянные материала определяются из следующих условий. Компоненты тензора напряжений и обращаются в нуль, когда компоненты тензора поврежденности достигают своего критического значения В условиях трехосного растяжения . Тогда должны выполняться два условия Откуда легко найти значения материальных постоянных: Таким образом, выбор значений материальных постоянных представляет собой нетривиальную задачу и должен осуществляться в каждом случае отдельно, исходя из рассматриваемой задачи. 2. Пользовательская процедура UMAT МКЭ-пакета Simulia Abaqus В настоящем параграфе будут приведены конечно-элементные решения задач о трещинах в среде с поврежденностью в связанной постановке (упругость-поврежденность). Определяющие соотношения (3)-(7) описаны в многоцелевом расчетном комплексе Simulia Abaqus с помощью пользовательской процедуры UMAT. Пользовательские процедуры многофункционального расчетного комплекса Simulia Abaqus позволяют ввести в рассмотрение материалы, отсутствующие в стандартном наборе материальных моделей классического пакета Abaqus/CAE. Все конечно-элементные решения получены для определяющих уравнений (3)-(7). Интерес представляет оценка взаимного влияния эволюции напряженно-деформированного состояния в образце и процесса накопления повреждений, в особенности у вершины трещины или надреза. Первой задачей является задача о всестороннем и одноосном растяжении пластины с центральным круговым отверстием в среде с поврежденностью. Затем рассматриваются задачи определения напряженно-деформированного состояния и поврежденности у вершин трещин в условиях нормального отрыва и смешанного нагружения. 2.1. Распределение поврежденности в пластине с центральным круговым отверстием В качестве первого образца была выбрана пластина с центральным круговым отверстием под действием 1) всестороннего и 2) одноосного растяжения в материале с определяющими уравнениями (3)-(7). На рис. 1 и 2 показаны распределения поврежденности в пластине с круговым вырезом в условиях всестороннего и одноосного растяжения для определяющих уравнений с включенным тензором поврежденности. Из рис. 1 и 2 видно, что максимальные значения поврежденности достигаются на контуре кругового выреза в случае всесторонне растягиваемой пластины и в точках кругового отверстия, лежащих на оси ординат. Поскольку распределения поврежденности, полученные в ходе конечно-элементного эксперимента, соответствуют известным экспериментальным результатам [48, 49], рассматриваемая материальная модель была использована для решения более сложных задач, для анализа зон активного накопления повреждений у вершины трещины. Рис. 1. Распределение компоненты тензора поврежденности в пластине с центральным круговым отверстием под действием всестороннего растяжения Fig. 1. Distribution of damage tensor component in the plate under equi-biaxial tension of the specimen Рис. 2. Распределение компоненты тензора поврежденности в пластине с центральным круговым отверстием под действием одноосного растяжения по оси Fig. 2. Distribution of damage tensor component in the plate for uniaxial tension of the specimen 2.2. Растяжение пластины с центральной трещиной в среде с поврежденностью Особый интерес представляет определение конфигурации зон накопления повреждений у вершины трещины. Несмотря на достаточно богатый опыт работы с пользовательскими процедурами пакета SIMULIA Abaqus, накопленный исследователями в нашей стране и за рубежом [17-20, 32, 46, 50, 51], позволяющий моделировать широкий класс материалов с инкорпорированными скалярной и тензорной мерой поврежденности, отсутствуют результаты вычислений, показывающих процессы накопления повреждений вблизи кончика трещины, и их систематичный анализ. В силу указанных причин можно остановиться на простой модели пластины с центральной трещиной. Типичная схема разбиения области, окружающей трещины, показана на рис. 3. Рис. 3. Пример конечно-элементного разбиения пластины на конечные элементы: вершины трещины окружают сингулярные элементы Fig. 3. Details of finite element mesh pattern for the simulation the plate with a horizontal central crack and two crack tips with singular elements Модель анизотропной поврежденности с включенными компонентами тензора поврежденности в определяющие уравнения линейно-упругого материала была применена для оценки зон активного накопления повреждений у вершин трещины. На рис. 4-8 приведены результаты вычислений. На рис. 4-7 показаны распределения напряжений в пластине без учета процесса накопления повреждений (слева) и с учетом поврежденности материала (справа). Компоненты тензора поврежденности показаны на рис. 8. Из рис. 4-7 видно, что накопление повреждений приводит к снижению значений всех компонент тензора напряжений. Таким образом, влияние поврежденности в материале приводит к снижению особенности поля напряжений вблизи вершины трещины. На рис. 8, а и б приведены распределения компонент тензора поврежденности и в пластине, связанные с максимальным и минимальным главными напряжениями. В данном расчете компонента связана со средним главным напряжением а б Рис. 4. Распределение интенсивности напряжений в пластине без учета накоплений повреждений (а) и с учетом накопления повреждений (б) Fig. 4. Distribution of stress intensities in the plate with no consideration of damage accumulation (a), and with the consideration of damage accumulation (b) а б Рис. 5. Распределение компоненты напряжений в пластине без учета накоплений повреждений (а) и с учетом накопления повреждений (б) Fig. 5. Distribution of stress component in the plate with no consideration of damage accumulation (a), with the consideration of damage accumulation (b) а б Рис. 6. Распределение компоненты напряжений в пластине без учета накоплений повреждений (а) и с учетом накопления повреждений (б) Fig. 6. Distribution of stress component in the plate with no the consideration of damage accumulation (a), with the consideration of damage accumulation (b) а б Рис. 7. Распределение компоненты напряжений в пластине без учета накоплений повреждений (а) и с учетом накопления повреждений (б) Fig. 7. Distribution of stress component in the plate: in the plate with no the consideration of damage accumulation (a), with the consideration of damage accumulation (b) а б c Рис. 8. Распределение компонент тензора поврежденности: a - распределение компоненты ; б - распределение компоненты ; в - распределение компоненты Fig. 8. Distribution of damage tensor components: a - is the distribution of damage component ; b - is damage component distribution; c - is damage component distribution 2.3. Распределение поврежденности у вершин наклонной трещины в растягиваемой пластине Для моделирования смешанного нагружения в полном диапазоне смешанных форм деформирования в работе была выполнена серия конечно-элементных расчетов для наклонной трещины в растягиваемой по оси ординат пластине. На рис. 9-11 показаны распределения компонент тензора поврежденности для наклонной под углом трещины. На рис. 12-16 приведены распределения интенсивности напряжений и компонент тензора напряжений в пластине с наклонной трещиной без учета накоплений повреждений (слева) и с учетом накопления повреждений (справа). Сравнение результатов расчетов для неповрежденного и поврежденного материалов показало, что учет процессов накопления повреждений приводит к снижению значений компонент тензора напряжений и интенсивности касательных напряжений в образцах с концентраторами напряжений, что хорошо согласуется с аналитическими и численными решениями краевых задач о трещинах [25, 26]: учет процесса накопления повреждений приводит к ослаблению или полному устранению особенности поля напряжений у вершины трещины в изотропных линейно-упругих телах. Рис. 9. Распределение компоненты тензора поврежденности в растягиваемой пластине Fig. 9. Damage tensor component distribution in the cracked plate under tension Рис. 10. Распределение компоненты тензора поврежденности в растягиваемой пластине Fig. 10. Damage tensor component distribution in the cracked plate under tension Рис. 11. Распределение компоненты тензора поврежденности в растягиваемой пластине с наклонной трещиной Fig. 11. Damage tensor component distribution in the cracked plate under tension а б Рис. 12. Распределение интенсивности напряжений в пластине с наклонной трещиной без учета накоплений повреждений (а) и с учетом накопления повреждений (б) Fig. 12. Distribution of stress intensities in the plate with an inclined crack: with no consideration of damage accumulation (a), and with the consideration of damage accumulation (b) а б Рис. 13. Распределение компоненты напряжений в пластине без учета накоплений повреждений (а) и с учетом накопления повреждений (б) Fig. 13. Stress component distribution in the plate: with no consideration of damage accumulation (a), and with the consideration of damage accumulation (b) а б Рис. 14. Распределение компоненты напряжений в пластине без учета накоплений повреждений (а) и с учетом накопления повреждений (б) Fig. 14. Stress component distribution: with no consideration of damage accumulation (a), and with the consideration of damage accumulation (b) а б Рис. 15. Распределение компоненты напряжений в пластине без учета накоплений повреждений (а) и с учетом накопления повреждений (б) Fig. 15. Stress component distribution: with no consideration of damage accumulation (a), and with the consideration of damage accumulation (b) а б Рис. 16. Распределение компоненты напряжений в пластине без учета накоплений повреждений (а) и с учетом накопления повреждений (б) Fig. 16. Stress component distribution: with no consideration of damage accumulation (a), and with the consideration of damage accumulation (b) 2.4. Распределение поврежденности у вершины надреза в полудиске Одним из часто рассматриваемых образцов является полудиск с вертикальным и наклонным надрезом [52-55]. Геометрия образца и система приложенных нагрузок показана на рис. 17. Рис. 17. Схематичное изображение образцов с надрезом для исследования нормального отрыва (а) и смешанного нагружения (б) Fig. 17. General configuration of the semicircular bend (SCB) specimen with a cut for Mode I loading (a) and Mixed Mode loadings (b) На рис. 17 приняты следующие обозначения: - приложенная нагрузка; - длина выреза; - радиус диска; - угол наклона надреза. Меняя угол наклона надреза к вертикали, можно менять вид смешанного нагружения от нормального отрыва до нагружений, близких к поперечному сдвигу. В работе выполнена большая серия вычислительных экспериментов на трехточечный изгиб диска с надрезом в полном диапазоне значений угла . Результаты вычислений показаны на рис. 18-26. На примере данной конфигурации была исследована зависимость распределения напряжений и поврежденности от вида разбиения, что является характерной особенностью задач механики поврежденности [22]. Проведенные расчеты не показали такой зависимости, и полученные распределения компонент тензора поврежденности (рис. 18-21), и компонент тензора напряжений (рис. 22-26) не зависят от способа разбиения, от вида конечных элементов. В ходе вычислений были a б в Рис. 18. Распределение компонент тензора поврежденности в диске с надрезом (нормальный отрыв): а - распределение компоненты ; б - распределение компоненты ; в - распределение компоненты Fig. 18. Distribution of the anisotropic damage tensor components: a - is damage tensor component ; b - is damage tensor component ; c - is damage tensor component . проведены расчеты на различных сетках: менялось количество конечных элементов сетки и сравнивались значения искомых величин (компонент тензоров напряжений и повреждений) в совпадающих точках. Сравнение результатов показало, что измельчение сетки не ведет к изменению вычисленных значений для полей напряжений и поврежденности. Распределение компонент тензора поврежденности в образце с наклонным под углом надрезом показано на рис. 19-21. Рис. 19. Распределение компоненты тензора поврежденности Fig. 19. Distribution of the anisotropic damage tensor component Рис. 20. Распределение компоненты Fig. 20. Distribution of the anisotropic damage tensor component Рис. 21. Распределение компоненты Fig. 21. Distribution of the anisotropic damage tensor component а б Рис. 22. Распределение интенсивности напряжений в полудиске с наклонным надрезом без учета процесса накопления повреждений (а) и с учетом процесса накопления повреждений (б) Fig. 22. Distribution of the Mises equivalent stress in the SCB specimen with an inclined cut with no consideration of damage accumulation process (a) and with the consideration of damage accumulation processes (b) Рис. 23. Распределение компоненты тензора напряжений в полудиске с наклонным под углом надрезом Fig. 23. Distribution of the stress tensor component in the SCB specimen with an inclined cut with an angle of Рис. 24. Распределение компоненты тензора напряжений в полудиске с наклонным под углом надрезом Fig. 24. Distribution of the stress tensor component in the SCB specimen with an inclined angle of Рис. 25. Распределение компоненты тензора напряжений в полудиске с наклонным под углом надрезом Fig. 25. Distribution of the stress tensor component in the SCB specimen with an inclined angle of Рис. 26. Распределение компоненты тензора напряжений в полудиске с наклонным под углом надрезом Fig. 26. Distribution of the stress tensor component in the SCB specimen with an inclined angle of Выводы и обсуждение результатов В настоящей работе описана и реализована процедура расчета поля анизотропной поврежденности в образцах с трещинами посредством моделирования определяющих уравнений связанной постановки задачи в связке упругость-поврежденность с помощью модели анизотропной поврежденности. В конечно-элементном комплексе SIMULIA Abaqus с помощью пользовательской процедуры UMAT реализована математическая модель материала с включенными компонентами тензора поврежденности второго ранга. Учтен анизотропный характер накопления повреждений в образцах с различной геометрией и различными системами нагрузок: рассмотрены пластина с центральной горизонтальной и наклонной трещиной и полудиск с вертикальным и наклонными надрезами. Получены и проанализированы зоны активного накопления повреждений вблизи кончика трещины и надреза для целой серии экспериментальных образцов с дефектами в условиях смешанного нагружения в полном диапазоне смешанных форм деформирования от чистого нормального отрыва до нагружений, близких к поперечному сдвигу. На основании проведенных расчетов для различных образцов с трещинами и надрезами показано, что накопление поврежденности приводит к снижению уровня напряжений вблизи вершины трещины, что хорошо согласуется с имеющимися аналитическими и численными решениями. На основании проведенных вычислительных экспериментов можно заключить, что, по всей видимости, развитие практических приложений механики поврежденности будет идти по пути разработки комплексов вычислительных программ, позволяющих: 1) оценить напряженно-деформированное состояние и поле повреждений в элементах конструкции; 2) определить конфигурацию зон локализаций необратимых деформаций и зон активного накопления повреждений, что даст возможность оптимизировать массогабаритные характеристики элементов конструкций. Реализация имитационного (компьютерного) моделирования на основе экспериментальных и теоретических исследований для модельных и используемых в теплоэнергетике, авиа- и автомобилестроении элементов конструкций позволит выработать и сформулировать критерии разрушения, учитывающие изменение свойств материала в процессе работы элемента конструкции, дать современные варианты ремонтных технологий.

About the authors

L V Stepanova

Samara University

References

  1. Качанов Л.М. О времени разрушения в условиях ползучести // Изв. АН СССР. ОТН. - 1958. - С. 26-31.
  2. Работнов Ю.Н. О механизме длительного разрушения // Вопросы прочности материалов и конструкций. - М.: Изд-во АН СССР, 1959. - С. 5-7.
  3. Работнов Ю.Н. Избранные труды. Проблемы механики деформируемого твердого тела. - М.: Наука, 1991. - 196 с.
  4. Кукуджанов В.Н. Компьютерное моделирование деформирования, повреждаемости и разрушения неупругих материалов и конструкций. - М.: Изд-во МФТИ, 2008. - 215 с.
  5. Kukudzhanov V.N. Numerical Continuum Mechanics. - Berlin: De Gruyer, 2012. - 425 p.
  6. Кукуджанов В.Н. Связанные модели упругопластичности и поврежденности и их интегрирование // Изв. РАН. МТТ. - 2006. - № 6. - C. 103-135.
  7. Fracture of a borosilicate glass under triaxial tension / V. Doquet, N. Ben Ali, A. Constantinescu, X. Boutillon // Mechanics of Materials. - 2013. - Vol. 57. - P. 15-29.
  8. Experimental and numerical study of crack healing in a nuclear glass / V. Doquet, N. Ben Ali, E. Chabert, F. Bouyer // Mechanics of Materials. - 2015. - Vol. 80. - P. 145-162.
  9. Modeling of thermal shock-induced damage in a borosilicate glass / M. Dube, V. Doquet, A. Constantinescu, D. George, Y. Remond, S. Ahzi // Mechanics of Materials. - 2010. - Vol. 42. - P. 863-872.
  10. Sun X., Khaleel M.A. Modeling of glass fracture damage using continuum damage mechanics-static spherical indentation // International Journal of Damage Mechanics - 2004. - Vol. 13. - P. 263-284.
  11. Sun X., Khaleel M.A., Davies R.W. Modeling of stone-impact resistance of monolithic glass ply using continuum damage mechanic // International Journal of Damage Mechanics. - 2005. - Vol. 14. - P. 165-178.
  12. Computational modeling of static indentation - induced damage in glass / J. Ismail, F. Zaiiri, M. Nait-Abdelaziz, Z. Azari // Computational Materials Science. - 2008. - Vol. 42. - P. 407-415.
  13. Analysis of damage in laminated automotive glazing subjected to simulated head impact / S. Zhao, L.R. Dharani, L. Chai, S. Barbat // Engineering Failure Analysis. - 2006. - Vol. 13. - P. 582-597.
  14. Damage interaction and angle effects on the erosion behavior of soda-lime-silica glass / Q. Jiang, J. Ismail, F. Zairi, M. Nait-Abdelaziz, Z. Qu, Z. Azari, X. Liu // International Journal of Damage Mechanics. - 2017.
  15. Experimental and numerical investigations on erosion damage in glass by impact of small-sized particles / J. Ismail, Zaïri Fahmi, M. Naït-Abdelaziz, S. Bouzid, Z. Azari // Wear. - 2011. - Vol. 271. - P. 817-826.
  16. Damage mechanisms in bioactive glass matrix composites under uniaxial compression / Q. Jiang, J. Ismail, F. Zairi, Z. Qu, X. Liu, F. Zairi // Journal of the Mechanical Behavior of Biomedical Materials. - 2018. - Vol. 79. - P. 264-272.
  17. Извеков О.Я., Крупеник А.М. Решение связанных задач континуального разрушения термо- и пороупругих сред в SIMULIA ABAQUS [Электронный ресурс]. - URL: http://www.tesis.com/software/abaqus/abaqus-exp.php.
  18. Извеков О.Я., Кондауров В.И. Модель пористой среды с упругим трещиноватым скелетом // Изв. РАН. Физика Земли. - 2009. - № 4. - С. 31-42.
  19. Кондауров В.И. Механика и термодинамика насыщенной пористой среды. - М.: Изд-во МФТИ, 2007. - 310 c.
  20. Извеков О.Я., Кондауров В.И. О рассеянном разрушении пористых материалов с хрупким скелетом // Изв. РАН. Механика твердого тела. - 2010. - № 3. - С. 164-187.
  21. Voyiadjis G.Z. Handbook of Damage Mechanics. Nano to Macro Scale for Materials and Structures. - Dordrecht: Springer, 2014. - 1000 p.
  22. Murakami S. Continuum Damage Mechanics. A Continuum Mechanics Approach to the Analysis of Damage and Fracture. - Dordrecht: Springer, 2012. - 423 p.
  23. Dormieux L., Kondo D. Micromechanics of Fracture and Damage (Mechanical Engineering and Solid Mechanics: Micromechanics Set). - Hebeken: Wiley-ISTE, 2016. - 312 p.
  24. Работнов Ю.Н. Ползучесть элементов конструкций. - М.: Наука, 2014. - 752 с.
  25. Степанова Л.В., Игонин С.А. Параметр поврежденности Ю.Н. Работнова и описание длительного разрушения: результаты, современное состояние, приложение к механике трещин // Прикладная механика и техническая физика. - 2015. - Т. 56, № 2(330). - С. 133-145.
  26. Степанова Л.В. О собственных значениях в задаче о трещине антиплоского сдвига в материале со степенными определяющими уравнениями// Прикладная механика и техническая физика. - 2008. - Т. 49, № 1(287). - С. 173-180.
  27. Степанова Л.В., Яковлева Е.М. Асимптотика собственных значений нелинейной задачи на собственные значения, следующей из проблемы определения напряженно-деформированного состояния у вершины трещины в условиях смешанного нагружения // Сиб. журн. вычислительной математики. - 2016. - Т. 19, № 2. - С. 207-222.
  28. Stepanova L., Yakovleva E. Stress-strain state near the crack tip under mixed-mode loading: Asymptotic approach and numerical solutions of nonlinear eigenvalue problems// AIP Conference Proceedings. - 2016. - Vol. 1785. - 030030.
  29. Степанова Л.В., Адылина Е.М. Напряженно-деформированное состояние в окрестности вершины трещины в условиях смешанного нагружения// Журнал прикладной механики и технической физики. - 2014. - Т. 55, № 5 (327). - С. 181-194.
  30. Степанова Л.В., Яковлева Е.М. Смешанное деформирование пластины с трещиной в условиях плоского напряженного состояния // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2014. - № 3. - С. 129-162.
  31. Stepanova L., Mironova E. Intermediate Asymptotic behavior of the stress and damage fields in the vicinity of the mixed-mode crack tip under creep regime // Ceur Workshop Proceedings. - 2017. - Vol. 1904. - P. 141-150.
  32. Моделирование локально-функциональных свойств материала, поврежденного полями дефектов / Е.В. Ломакин, С.А. Лурье, П.А. Белов, Л.Н. Рабинский // Докл. Акад. наук. - 2017. - Т. 472, № 3. - С. 282-285.
  33. Wen Z.X., Hou N.X., Yue Z.F. Creep damage and crack initiation behaviour of nickel-base single crystalline superalloys compact tension specimen with a void ahead of crack tip // Materials Science and Engineering A. - 2009. - Vol. 510-511. - P. 284-288.
  34. Development of a new nonlinear numerical material model for woven composite materials accounting for permanent deformation and damage / O. Cousigne, D. Moncayo, D. Coutellier, P. Camanho, H. Naceur, S. Hampel // Composite Structures. - 2013. - Vol. 106. - P. 601-614.
  35. Вильдеман В.Э., Соколкин Ю.В., Ташкинов А.А. Механика неупругого деформирования и разрушения композиционных материалов. - М.: Наука: Физматлит, 1997. - 288 с.
  36. Волегов П.С., Грибов Д.С., Трусов П.С. Поврежденность и разрушение: обзор экспериментальных работ // Физическая мезомеханика. - 2015. - Т. 18, № 3. - С. 11-24.
  37. Волегов П.С., Грибов Д.С., Трусов П.С. Поврежденность и разрушение: классические континуальные теории // Физическая мезомеханика. - 2015. - Т. 18, № 4. - С. 68-87.
  38. Волегов П.С., Грибов Д.С., Трусов П.С. Поврежденность и разрушение: модели, основанные на физических теориях пластичности // Физическая мезомеханика. - 2015. - Т. 18, № 6. - С. 12-23.
  39. Barbero E.J., Shahbazi M. Determination of material properties for ANSYS progressive damage analysis of laminated composites // Composite Structures. - 2017. - Vol. 176. - P. 768-779.
  40. Wan Y., Sun B., Gu B. Multi-scale structure modeling of damage behaviors of 3D orthogonal woven composite materials subject to quasi-static and high strain rate compressions // Mechanics of Materials. - 2016. - Vol. 94. - P. 1-25.
  41. Van Dongen B., Van Oostrum A., Zarouchas D. A blended continuum damage and fracture mechanics method for progressive damage analysis of composite structures using XFEM // Composite Structures. - 2018. - Vol. 184. - P. 512-522.
  42. Jin W., Arson C. Nonlocal enrichment of a micromechanical damage model with tensile softening: Advantages and limitations // Computers and Geotechnics. - 2018. - Vol. 94. - P. 196-206.
  43. Meng Q., Wang Z. Creep damage models and their applications for crack growth analysis in pipes: A review // Engineering Fracture Mechanics. - 2016.
  44. Fengxia Q. Abaqus implementation of creep failure in polymer matrix composites with transverse isotropy. - Ph. D. Thesis. 2005. - 90 p.
  45. Бойл Дж., Спенс Дж. Анализ напряжений в конструкциях при ползучести. - М.: Мир, 1986. - 360 с.
  46. Vershinin V.V. A correct form of Bai-Wierzbicki plasticity model and its extension for strain rate and temperature dependence // International Journal of Solids and Structures. - 2017. - Vol. 126-127. - P. 150-162.
  47. Gao Z., Zhang L., Yu W. A nonlocal continuum damage model for brittle fracture// Engineering Fracture Mechanics. - 2018. - Vol. 189. - P. 481-500.
  48. Hyakutake H., Yamamoto T. Damage near the notch root of notched FRP plates in static load-evaluation of damage by luminance-measuring system // Transactions on Engineering Sciences. - 1996. - Vol. 13. - P. 417-424.
  49. Strength degradation and stress analysis of composite plates with circular, square and rectangular notches using digital image correlation / A. Khechari, A. Tati, B. Guerira, A. Guettala, P.M. Mohite // Composite Structures. - 2018. - Vol. 185. - P. 699-715.
  50. Saeedi M.R., Morovvati M.R., Alizadeh-Vaghasloo Y. Experimental and numerical study of mode-I and mixed-mode fracture of ductile U-notched functionally graded materials// International Journal of Mechanical Sciences. - 2018. - P. 1-33.
  51. Ломакин Е.В., Федулов Б.Н. Растяжение полосы, ослабленной вырезами с круговым основанием, в условиях плоской деформации из материала с зависящими от вида напряженного состояния свойствами // Изв. РАН. Механика твердого тела. - 2013. - № 4. - С. 80-87.
  52. Fakhri M., Haghighat Kharrazi E., Aliha M.R.M. Mixed mode tensile - In plane shear fracture energy determination for hot mix asphalt mixtures under intermediate temperature conditions // Engineering Fracture Mechanics. - 2018. - Vol. 192. - P. 98-113.
  53. Malikova L., Vesely V., Seitl S. Crack propagation direction in a mixed mode geometry estimated via multi-parameter fracture criteria // International Journal of Fatigue. - 2016. - Vol. 89. - P. 99-107.
  54. Mirsayar M.M., Razmi A., Berto F. Tangential strain-based criteria for mixed-mode I/II fracture toughness of cement concrete // Fatigue Fract. Eng. Mater. Struct. - 2017. - P. 1-9.
  55. EMTSN criterion for evaluating mixed mode I/II crack propagation in rock materials / M.M. Mirsayar, A. Razmi, M.R.M. Aliha, F. Berto // Engineering Fracture Mechanics. - 2018. - Vol. 190. - No. 1. - P. 186-197.

Statistics

Views

Abstract - 625

PDF (Russian) - 284

Cited-By


PlumX


Copyright (c) 2018 Stepanova L.V.

Creative Commons License
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies