SELF-SIMILARITY OF THE NEAR CRACK-TIP FIELDS UNDER CREEP REGIME TAKING INTO ACCOUNT DAMAGE AUGMENTATION
- Authors: Bykova Y.S1, Stepanova L.V1
- Affiliations:
- Samara National Research University, Samara, Russian Federation
- Issue: No 3 (2025)
- Pages: 54-69
- Section: ARTICLES
- URL: https://ered.pstu.ru/index.php/mechanics/article/view/4653
- DOI: https://doi.org/10.15593/perm.mech/2025.3.05
- Cite item
Abstract
The paper analyzes a series of computer experiments aimed at identifying the self-similar behavior of stresses and continuity (damage) near the crack tip in a steady-state creep mode in a damaged medium. Finite element computations of crack-tip fields under creep regime were carried out using the interdisciplinary, universal finite element platform SIMULIA Abaqus FEA using the UMAT utility, which integrates the process of damage augmentation into the FEM computational scenario. The paper implements computational simulations of uniaxial stretching of a plate weakened by a central horizontal crack in the steady-state creep mode, which includes damage growth that evolves over time according to the mathematical model of damage growth by Kachanov – Rabotnov (KR) according to a power law for various values of exponents of the kinetic equation and the power constitutive equations. The study and analysis of the FE for the crack-tip stress and continuity (integrity) fields for a number of material constants clearly reveal a self-similar behavior of stress and damage fields near the front tip. The structure of the solution is revealed and the values of the exponents in the self-similar variable and the self-similar representation of the solution are found, which can be interpreted as an intermediate self-similar solution of the second type according to the classification of G.I. Barenblatt. It is elucidated that the revealed property of self-similarity of the solution can be interpreted as the intermediate asymptotes of the far field of damage and stress. The obtained figures also clearly show the asymptotes of the near-field stress, characterized by the absence of a singularity in the immediate vicinity of the crack.
Full Text
Ползучесть – важное и сложное явление, часто наблюдаемое при анализе разрушения инженерных материалов и компонент, работающих при высоких температурах, включая материалы, используемые в аэрокосмических и авиационных двигателях, конструкциях энергетической промышленности [1-7]. В инженерных и научных исследованиях понимание и прогнозирование явлений ползучести имеют решающее значение для обеспечения безопасности, надежности и срока службы конструкций и материалов [1-7]. Сложность механики ползучести обусловлена разнообразием ее физических основ. Микроструктура и организация материалов, диффузия атомов, дефекты кристаллической решетки, движение дислокаций и многие другие факторы играют важную роль в процессе ползучести. Эти сложные физические процессы приводят к нелинейности и временной зависимости деформации ползучести, что затрудняет создание точных механических моделей и соответствующих им расчетных схем. При описании деформации ползучести принято выделять первую, вторую и третью стадии ползучести. Наряду с этим для учета тесного взаимодействия различных механизмов и их взаимосвязей в последнее время уделяется большое внимание интеракции ползучести и усталости и их влиянию на процессы разрушения. На накопленную деформацию и развитие поврежденности при ползучести в значительной мере оказывают влияние температура, напряжение, время и внешние факторы - окружающая среда. Численное моделирование аугментации повреждений при ползучести в твердых телах создает значительные трудности из-за присущей этим явлениям сложности и зависимости от времени [1-7]. Повреждения при ползучести, характеризующиеся постепенной деформацией материалов под действием постоянного напряжения, особенно важны в инженерных приложениях, где важна долговременная целостность конструкции [1-7]. За последние десятилетия численные методы получили существенное развитие и нашли успешное применение в инженерных задачах [1-7]. Стремительное развитие компьютерных технологий и программного обеспечения значительно расширило применение методов численного моделирования в исследованиях разрушения при ползучести. Начиная с метода конечных элементов (МКЭ), постоянные инновации в различных методах и технологиях численного моделирования позволили исследователям более точно моделировать и прогнозировать сложные процессы разрушения при ползучести [8-18]. Эти достижения значительно снизили затраты на прогнозирование разрушения при ползучести и продолжают совершенствоваться, повышая точность и эффективность моделирования и стимулируя постоянное развитие исследований механизмов и процессов разрушения при ползучести. Компьютерное моделирование, основанное на методе конечных элементов, остается одним из базисных подходов для симуляции роста повреждений [8-18]. Наряду с хорошо разработанными и апробированными подходами предлагаются и развиваются принципиально новые приемы вычислений [8,9]. Например, в [8] предложена усовершенствованная перидинамическая модель для моделирования вязкоупругой деформации и повреждений при ползучести. Ключевые элементы предлагаемой формулировки включают в себя перидинамические конституциональные уравнения зависящей от времени определяющей модели и неявную дискретизацию перидинамических уравнений. Перидинамическое моделирование деформации ползучести и повреждений, в котором векторы плотности сил в уравнениях равновесия перидинамики (ПД) получены путем рассмотрения модели ползучести Лю и Мураками с параметром повреждения, выполнено в [9]. Градиент деформации определяется с помощью дифференциального оператора перидинамики. Граничные условия сцепления и перемещения непосредственно задаются с помощью новой стратегии при решении уравнений равновесия ПД. ПД-форма элементов сцепления позволяет создавать условия сцепления в реальной области “пограничного слоя” без каких-либо нефизических перегибов вблизи границ. Этот подход был проверен для случаев одноосного и двухосного нагружений с учетом деформации ползучести, вызванной постоянными напряжениями при высоких температурах. Авторы показывают [9], что предсказания деформации при ползучести находятся в отличном согласии с экспериментальными данными и аналитическими решениями. Франкфорт и Мариго [10] усовершенствовали энергетическую теорию разрушения и предложили вариационную модель линейного упругого разрушения, которая разлагает общую энергию конструкции на сумму энергии деформации при деформировании и энергии разрушения при распространении трещины. В ней утверждается, что ранее существовавшие трещины и дефекты в материале будут расширяться под воздействием внешних нагрузок, что приведет к снижению общей энергии системы. На основе этой модели было разработано дополнительное скалярное фазовое поле, вводимое для описания процесса эволюции повреждений в материале через его вариации, что позволило сформировать основу метода фазового поля разрушения (МФП). Реализация дискретного расчета МФП обычно основана на МКЭ, но процесс генерации сетки в традиционном МКЭ является громоздким и часто использует полигональные сетки для дискретизации геометрических структур [11]. Это неизбежно приводит к геометрическим ошибкам, которые снижают точность расчетов напряжений и деформаций при обработке сложных траекторий трещин и детальных геометрических характеристик. Изогеометрический анализ напрямую использует геометрические модели автоматизированного проектирования для анализа и вычислений и позволяет избежать внесения геометрических ошибок и обеспечивая более точное геометрическое представление. Кроме того, МФП может легко удовлетворить требования к непрерывности более высокого порядка, что делает его весьма подходящим для моделирования тонких вариаций геометрии и сложных структур. Исследование [12] обеспечивает формирование и расчетно-экспериментальное обоснование подхода для прогнозирования живучести и долговечности элементов конструкций, базирующегося на принципах допустимой поврежденности и имитационного моделирования. Авторы [12] дополнительно отмечают, что подходы и методы, используемые в исследовании, тесно взаимосвязаны и в значительной степени дополняют друг друга. Несомненной инновационной новизной результатов является комплексный подход к оценке целостности конструкции, который сочетает в себе формулировку определяющих уравнений с функциями повреждения, численное моделирование на основе МКЭ и метода фазового поля разрушения, экспериментальные исследования усталости и роста трещин с прямыми измерениями полей деформаций с использованием метода DIC. Авторы [12] еще раз показывают, что численное моделирование имеет важное значение для практического применения, поскольку устраняет недостатки натурных испытаний конструкций, которые не обладают необходимой статистической надежностью и ограничены для прямых измерений процессов накопления и роста повреждений в процессе эксплуатации. Одним из быстро развивающихся направлений исследований является учет особенностей совместного влияния циклического нагружения и накопления деформаций ползучести [13,14]. В [13] предложена модель накопления усталостных повреждений на третьей стадии ползучести, учитывающая и описывающая двухстороннее, взаимное влияние усталостных повреждений и повреждений, аккумулируемых при ползучести в элементах конструкций, работающих в режимах периодического (циклического) нагружения при повышенных температурах. Экспериментально-численная модель, на основе выбранного эволюционного закона накопления повреждений Лангеберга [13] при ползучести, скомбинированной с усталостными повреждениями, применена для прогнозирования срока службы однонаправленных ламинатов с произвольными углами ориентации [13]. Следовательно, модель, развитая в [13], учитывает 1) воздействие температуры на прирост поврежденности при ползучести и 2) аккрецию повреждений при ползучести в различных условиях действия циклических нагрузок при повышенных температурах (т.е. при одновременном протекании процессов усталости и ползучести). Обзор современных вычислительных подходов может быть найден в [15], где детально описаны применяемые в самое последнее время численные методы механики разрушения и поврежденности. В этой работе представлено введение в популярные вычислительные модели механики повреждений и механики трещин, а также основные достижения и современное состояние применения вычислительных технологий в инженерии. Однако, несмотря на новейшие вычислительные подходы, обеспечивающие точное численное описание процессов деформирования и разрушения, занявшие доминирующее положение в исследованиях явлениях разрушения на различных масштабных уровнях, классический метод конечных элементов позволяет пролить свет на структуру аналитических решений задач механики разрушения и поврежденности. В настоящей работе для моделирования и анализа процессов накопления поврежденности применяется пользовательская процедура UMAT конечно-элементного комплекса SIMULIA Abaqus [16-20], встраивающая кинетическое уравнение накопления повреждений и степенной закон ползучести, соединяющий скорость деформации ползучести с истинными напряжениями в деформируемом теле, в вычислительный сценарий МКЭ, что дает возможность численно рассмотреть и оценить два взаимных, обоюдных процесса: 1) развитие поврежденности и эффекты его влияния на напряженно-деформированное состояние и 2) воздействие эволюции напряженно-деформированного состояния в условиях ползучести на рост повреждений в материале. В целом, идея внедрения пользовательских процедур в вычислительную реализацию МКЭ, не нова и активно применяется для расчета поля поврежденности у концентраторов напряжений [16-20]. В отличие от подходов, осуществленных в [16-20] и предыдущих работах авторов [21-25], ниже, на основании КЭ подхода с применением пользовательской утилиты UMAT, выявлен автомодельный характер процесса накопления повреждений. Опираясь на полученное МКЭ-решение, нам удалось обнаружить диапазоны изменения расстояния от вершины дефекта и временные интервалы, на которых справедливо автомодельное представление решения задачи о накоплении поврежденности у вершины трещины в условиях ползучести, подчиняющемся степенному эволюционному уравнению для среды со степенными конституциональными соотношениями.
About the authors
Y. S Bykova
Samara National Research University, Samara, Russian Federation
Author for correspondence.
Email: bykova.yus@ssau.ru
L. V Stepanova
Samara National Research University, Samara, Russian Federation
Email: stepanova.lv@ssau.ru
References
- Sandsrom, R. Basic Modelling and Theory of Creep of Metallic Materials/ R. Sandsrom. – Cham: Springer, 2024. – 317 p. https://doi.org/10.1007/978-3-031-49507-6
- Creep damage constitutive model of rock based on the mechanisms of crack-initiated damage and extended damage/ T. Li, C. Chen, F. Peng, C. Ma, M. Li, Y. Wang// Underground Space. – 2024. – Vol. 18. – P. 295-313. https://doi.org/10.1016/j.undsp.2023.12.008
- Chen, W. A thermodynamically consistent creep constitutive model considering damage mechanisms/ W. Chen, X. He, C. Sun, L. Meng// International Journal of Mechanical Sciences. – 2024. – Vol. 276. – 109373. https://doi.org/10.1016/j.ijmecsci.2024.109373
- Creep damage laws for bonded joints under pure mode II loading/ R.F.M. Couto, M.F.S.F. de Moura, A.G. Magalhães, R.D.F. Moreira// International Journal of Adhesion and Adhesives. – 2025. – Vol. 136. – 103879. https://doi.org/10.1016/j.ijadhadh.2024.103879
- Федоренков, Д. И. Реализация модели поврежденности Lemaitre с кинематическим упрочнением в конечно-элементном комплексе ANSYS / Д. И. Федоренков, Д. А. Косов // Вестник Пермского национального исследовательского политехнического университета. Механика. – 2022. – № 2. – С. 147-157. – doi: 10.15593/perm.mech/2022.2.12.
- Пантелеев, И. А. О выпуклости потенциала модели нелинейной упругой среды с тензорным параметром поврежденности / И. А. Пантелеев, В. Ляховский // Вестник Пермского национального исследовательского политехнического университета. Механика. – 2022. – № 1. – С. 89-101. – doi: 10.15593/perm.mech/2022.1.08.
- Хлыбов, А. А. О возможности спектрально-акустического контроля поврежденности поликристаллических материалов на базе имитационно-механической модели / А. А. Хлыбов, А. Л. Углов, Д. А. Рябов // Физическая мезомеханика. – 2023. – Т. 26, № 2. – С. 106-114. – doi: 10.55652/1683-805X_2023_26_2_106
- Wang, L. Peridynamic modelling of time-dependent behaviour and creep damage in hyper-viscoelastic solids with pre-cracks/ L. Wang, Z.-Y. Yin// Computer Methods in Applied Mechanics and Engineering. – 2025. – Vol. 433. – Part A. – 117512. https://doi.org/10.1016/j.cma.2024.117512
- Behera, D. Peridynamic simulation of creep deformation and damage/ D. Behera, P. Roy, E. Madenci// Continuum Mech. Thermodynamics. – 2024. – Vol. 36. – P. 1283–1304. https://doi.org/10.1007/s00161-024-01295-3
- Francfort, G.A. Revisiting brittle fracture as an energy minimization problem/ G.A. Francfort, J.-J. Marigo// Journal of Mechanics and Physics of Solids. – 1998. – Vol.46. – No. 8. – P. 1319-1342.
- Creep crack propagation using phase-field model within a multi-patch isogeometric framework/ Z. Si, Hirshikesh, T. Yu, W. Fang, S. Natarajan// Computer Methods in Applied Mechanics and Engineering. – 2025. – Vol. 435. – 117598. https://doi.org/10.1016/j.cma.2024.117598
- Khamidullin, R. Comprehensive study of the structural integrity of a damaged turbine disk using FEM, DIC and phase field methods/ R. Khamidullin, V. Shlyannikov, D. Kosov, A. Zakharov// International Journal of Fatigue. – 2025. – Vol. 192. – 108720. https://doi.org/10.1016/j.ijfatigue.2024.108720
- Fatigue-creep damage model for carbon fibre reinforced composites under high temperature cyclic loading/ Y.-E. Guo, D.-G. Shang, L.-X. Zuo, L.-F. Qu, C.-L. Chen// Composites Science and Technology. – 2024. Vol. 258. – 110909. https://doi.org/10.1016/j.compscitech.2024.110909
- Mechanisms of deformation, damage and life behavior of inconel 617 alloy during creep-fatigue interaction at 700 C/ Q. Wang, J. Yu, B. Li, Y. Li, K. Wang, X. Chen// International Journal of Fatigue. – 2025. – Vol. 190. – 108635. https://doi.org/10.1016/j.ijfatigue.2024.108635.
- Recent trends in computational damage models: An overview/ R. Piska, K. Sivadas, K. Boyina, A. Vuppuluri, A. Chaurasia, C. Parimi, T. Rabczuk// Theoretical and Applied Fracture Mechanics. – 2024. – Vol. 132. – 104494. https://doi.org/10.1016/j.tafmec.2024.104494
- Arruda, M.R.T. Jacobian vs. disturbance method for UMATs in ABAQUS: An application to isotropic damage mechanics/ M.R.T. Arruda, J. Shen// International Journal of Non-Linear Mechanics. – 2024. – Vol. 167. – 104928. https://doi.org/10.1016/j.ijnonlinmec.2024.104928.
- Lucarini, S. UMAT4COMSOL: An Abaqus user material (UMAT) subroutine wrapper for COMSOL/ S. Lucarini, E. Martínez-Pañeda// Advances in Engineering Software. – 20024. – Vol. 190. – 103610. https://doi.org/10.1016/j.advengsoft.2024.103610.
- OXFORD-UMAT: An efficient and versatile crystal plasticity framework/ E. Demir, A. Martinez-Pechero, C. Hardie, E. Tarleton// International Journal of Solids and Structures. – 2025. – Vol. 307. – 113110. https://doi.org/10.1016/j.ijsolstr.2024.113110.
- An insight into the creep failure mechanism of sliver defect in the second-generation nickel-based single crystal superalloy CMSX-4/ X. Yu, W. Xuan, C. Zhang, X. Zhang, X. Wang, Y. Zhao, B. Wang, H. Li, J. Bao, Z. Ren// Engineering Fracture Mechanics. – 2024. – Vol. 307. – 110307. doi: 10.1016/j.engfracmech.2024.110307
- Constitutive modelling and damage prediction of AlSi10Mg alloy manufactured by SLM technology with emphasis on ratcheting in LCF regime/ P. Das, R. Halama, A.V. Natarajan, N. Khutia, P.P. Dey, L. Kunčická, J. Hajnyš, R. Kocich, F. Sari// International Journal of Fatigue. – 2024. – Vol. 181. – 108115. https://doi.org/10.1016/j.ijfatigue.2023.108115
- Shlyannikov, V. Creep damage and stress intensity factor assessment for plane multi-axial and three-dimensional problems/ V. Shlyannikov, A. Tumanov// International Journal of Solids and Structures. – 2018. – Vol. 150. – P. 166-183. https://doi.org/10.1016/j.ijsolstr.2018.06.009
- Исследование асимптотики поля напряжений в окрестности вершины трещины в условиях ползучести с учетом поврежденности / Д. В. Чаплий, О. Н. Белова, Л. В. Степанова, Ю. С. Быкова // Вестник Пермского национального исследовательского политехнического университета. Механика. – 2024. – № 3. – С. 17-38. – doi: 10.15593/perm.mech/2024.3.02.
- Выявление асимптотического поведения поврежденности у фронта трещины при ползучести на основании конечно-элементных расчетов / Л. В. Степанова, О. Н. Белова, Д. В. Чаплий, Ю. С. Быкова // Прикладная механика и техническая физика. – doi: 10.15372/PMTF202415561
- Чаплий, Д. В. Параметрическое исследование полей, ассоциированных с вершиной трещины, в условиях ползучести с учетом процессов накопления поврежденности с использованием UMAT / Д. В. Чаплий, Л. В. Степанова, О. Н. Белова // Вестник Самарского государственного технического университета. Серия: Физико-математические науки. – 2023. – Т. 27, № 3. – С. 509-529. – doi: 10.14498/vsgtu2005.
- Быкова, Ю. С. Влияние процесса накопления повреждений на асимптотическое поведение полей напряжений в условиях ползучести образца с центральной трещиной / Ю. С. Быкова, Л. В. Степанова // Вестник Самарского университета. Естественнонаучная серия. – 2023. – Т. 29, № 4. – С. 7-25. – doi: 10.18287/2541-7525-2023-29-4-7-25
- ABAQUS/Standard User’s Manual, Version 2024, Dassault Systemes Simulia Corр, United States. – 2024
- On the prediction of creep behaviour of alloy 617 using Kachanov-Rabotnov model coupled with multi-objective genetic algorithm optimization/ J. Choi, L. Bortolan Neto, R.N. Wright, J.J. Kruzic, O. Muránsky// International Journal of Pressure Vessels and Piping. – 2022. – Vol. 199. – 104721. https://doi.org/10.1016/j.ijpvp.2022.104721
- Haque, M.S. Comparative analysis of the sin-hyperbolic and Kachanov–Rabotnov creep-damage model/ M. S. Haque, C. M. Stewart// International Journal of Pressure Vessels and Piping. – 2019. – Vol. 171. – P. 1-9. https://doi.org/10.1016/j.ijpvp.2019.02.001
- Viscoelastic–viscoplastic model with ductile damage accounting for tension–compression asymmetry and hydrostatic pressure effect for polyamide 66/ S. Satouri, G. Chatzigeorgiou, F. Meraghni, G. Robert// European Journal of Mechanics - A/Solids. – 2025. – Vol. 110. – 105491. https://doi.org/10.1016/j.euromechsol.2024.105491
- Качанов, Л.М. О времени разрушения в условиях ползучести // Изв. АН СССР. Отд. технических наук. – 1958. – № 8. – С. 26–31
- Работнов, Ю.Н. О механизме длительного разрушения // Вопросы прочности материалов и конструкций: сб. статей. Изд-во АН СССР. Москва. – 1959. – С. 5–7
- Meng Q., Wang Z. Creep damage models and their applications for crack growth analysis in pipes: A review// Engineering Fracture Mechanics. –2019. – Vol. 205. – P. 547-576
- Hutchinson, J.W. Singular behaviour at the end of a tensile crack in a hardening material/ J.W. Hutchinson// J Mech Phys Solids. – 1968. – V.16. – P. 13-3
- Rice, J. R. Plain strain deformation near a crack tip in a power-law hardening material/ J.R. Rice, G.F. Rosengren// J Mech Phys Solids. – 1968. – V. 16. – P. 1-12
- Hutchinson, J.W. Plastic stress and strain fields at a crack tip/ J.W. Hutchinson// J Mech Phys Solids. – 1968. – V. 16. – P. 337-347
- Barenblatt, G. I. Scaling, Self-similarity, and Intermediate Asymptotics. Cambridge University Press. 1996
- Barenblatt G.I. Flow, deformation and fracture lectures on fluid mechanics and mechanics of deformable solids for mathematicians and physicists. Cambridge: Cambridge University Press, 2014. 273 p
- Баренблатт Г.И. Подобие, автомодельность, промежуточная асимптотика. Теория и приложения к геофизической гидродинамике. М.: Ленанд, 2021. 256 с
- Riedel H. Fracture at High Temperature. Berlin: Springer. – 1987. – 418 p
- Астафьев, В.И., Радаев Ю.Н., Степанова Л.В. Нелинейная механика разрушения. – Самара: Издательство “Самарский университет”. – 2001. – 631 c
- Stepanova, L.V. Nonlinear eigenvalue problems arising from nonlinear fracture mechanics boundary value problems / L.V. Stepanova, E.M. Yakovleva // Procedia Structural Integrity. – 2022. – Vol. 37. – P. 908–919
- Степанова, Л.В. Смешанное деформирование пластины с трещиной в условиях плоского напряженного состояния/ Л.В. Степанова, Е.М. Яковлева// Вестник Пермского национального исследовательского политехнического университета. Механика. – 2014. – № 3. – С. 129-162. doi: 10.15593/perm.mech/2014.3.08
- Kong, W. A higher-order asymptotic solution for 3D sharp V-notch front tip fields in creeping solids/ W. Kong, Y. Dai, Y. Liu// Engineering Fracture Mechanics. – 2024. – Vol. 306. – 110256. doi: 10.1016/j.engfracmech.2024.11025