АВТОМОДЕЛЬНЫЙ ХАРАКТЕР ПОЛЕЙ ВБЛИЗИ ВЕРШИНЫ РАЗРЕЗА В ПРОЦЕССЕ ПОЛЗУЧЕСТИ В УСЛОВИЯХ АУГМЕНТАЦИИ ПОВРЕЖДЕННОСТИ

Аннотация


Статья посвящена обработке и анализу конечно-элементных (КЭ) расчетов выполненного цикла вычислительных экспериментов, на основании которых выявлен автомодельный характер распределения функции сплошности (поврежденности) и составляющих напряжения в непосредственной близости вершины разреза на третьей стадии ползучести в поврежденной среде. КЭ расчеты полей вблизи кончика дефекта проведены с использованием междисциплинарной, универсальной конечно-элементной платформы SIMULIA Abaqus с привлечением утилиты UMAT, интегрирующей процесс аугментации повреждений в вычислительный сценарий МКЭ. В работе реализовано компьютерное моделирование одноосного растяжения пластины, ослабленной центральным горизонтальным разрезом в режиме ползучести, в рамках которого в вычислительные алгоритмы включен прирост повреждений, прогрессирующий с времени согласно классической механической модели нарастания повреждений Качанова-Работнова (КР) по степенному закону для различных значений показателей степеней кинетического уравнения и степенного определяющего уравнения с концепцией истинного напряжения в связанной постановке. Исследование и анализ полученных КЭ полей напряжений и сплошности в окрестности кончика трещины для целого ряда постоянных материала явственно обнаруживает автомодельный характер распределения полей напряжений и поврежденности вблизи кончика дефекта степенного типа. Выявлена структура решения и найдены значения показателей степеней в автомодельной переменной и автомодельном представлении решения, которое можно интерпретировать как промежуточное автомодельное решение второго типа согласно классификации Г.И. Баренблатта. Показано, что полученное свойство автомодельности решения можно трактовать как асимптотику дальнего поля повреждений и напряжений. Также на извлеченных из МКЭ-расчетов зависимостях напряжения от расстояния от кончика разреза, воспроизведенных в двойных логарифмических координатах, наглядно проявляется асимптотическое поведение, отвечающее ближнему полю напряжений, характеризующееся полным отсутствием сингулярности в непосредственной близости вершины разреза.


Полный текст

Ползучесть – важное и сложное явление, часто наблюдаемое при анализе разрушения инженерных материалов и компонент, работающих при высоких температурах, включая материалы, используемые в аэрокосмических и авиационных двигателях, конструкциях энергетической промышленности [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, выявлен автомодельный характер процесса накопления повреждений. Опираясь на полученное МКЭ-решение, нам удалось обнаружить диапазоны изменения расстояния от вершины дефекта и временные интервалы, на которых справедливо автомодельное представление решения задачи о накоплении поврежденности у вершины трещины в условиях ползучести, подчиняющемся степенному эволюционному уравнению для среды со степенными конституциональными соотношениями.

Об авторах

Ю. С Быкова

Самарский национальный исследовательский университет имени академика С.П. Королева, Самара, Российская Федерация

Автор, ответственный за переписку.
Email: bykova.yus@ssau.ru

Л. В Степанова

Самарский национальный исследовательский университет имени академика С.П. Королева, Самара, Российская Федерация

Email: stepanova.lv@ssau.ru

Список литературы

  1. 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
  2. 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
  3. 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
  4. 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
  5. Федоренков, Д. И. Реализация модели поврежденности Lemaitre с кинематическим упрочнением в конечно-элементном комплексе ANSYS / Д. И. Федоренков, Д. А. Косов // Вестник Пермского национального исследовательского политехнического университета. Механика. – 2022. – № 2. – С. 147-157. – doi: 10.15593/perm.mech/2022.2.12.
  6. Пантелеев, И. А. О выпуклости потенциала модели нелинейной упругой среды с тензорным параметром поврежденности / И. А. Пантелеев, В. Ляховский // Вестник Пермского национального исследовательского политехнического университета. Механика. – 2022. – № 1. – С. 89-101. – doi: 10.15593/perm.mech/2022.1.08.
  7. Хлыбов, А. А. О возможности спектрально-акустического контроля поврежденности поликристаллических материалов на базе имитационно-механической модели / А. А. Хлыбов, А. Л. Углов, Д. А. Рябов // Физическая мезомеханика. – 2023. – Т. 26, № 2. – С. 106-114. – doi: 10.55652/1683-805X_2023_26_2_106
  8. 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
  9. 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
  10. 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.
  11. 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
  12. 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
  13. 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
  14. 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.
  15. 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
  16. 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.
  17. 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.
  18. 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.
  19. 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
  20. 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
  21. 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
  22. Исследование асимптотики поля напряжений в окрестности вершины трещины в условиях ползучести с учетом поврежденности / Д. В. Чаплий, О. Н. Белова, Л. В. Степанова, Ю. С. Быкова // Вестник Пермского национального исследовательского политехнического университета. Механика. – 2024. – № 3. – С. 17-38. – doi: 10.15593/perm.mech/2024.3.02.
  23. Выявление асимптотического поведения поврежденности у фронта трещины при ползучести на основании конечно-элементных расчетов / Л. В. Степанова, О. Н. Белова, Д. В. Чаплий, Ю. С. Быкова // Прикладная механика и техническая физика. – doi: 10.15372/PMTF202415561
  24. Чаплий, Д. В. Параметрическое исследование полей, ассоциированных с вершиной трещины, в условиях ползучести с учетом процессов накопления поврежденности с использованием UMAT / Д. В. Чаплий, Л. В. Степанова, О. Н. Белова // Вестник Самарского государственного технического университета. Серия: Физико-математические науки. – 2023. – Т. 27, № 3. – С. 509-529. – doi: 10.14498/vsgtu2005.
  25. Быкова, Ю. С. Влияние процесса накопления повреждений на асимптотическое поведение полей напряжений в условиях ползучести образца с центральной трещиной / Ю. С. Быкова, Л. В. Степанова // Вестник Самарского университета. Естественнонаучная серия. – 2023. – Т. 29, № 4. – С. 7-25. – doi: 10.18287/2541-7525-2023-29-4-7-25
  26. ABAQUS/Standard User’s Manual, Version 2024, Dassault Systemes Simulia Corр, United States. – 2024
  27. 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
  28. 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
  29. 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
  30. Качанов, Л.М. О времени разрушения в условиях ползучести // Изв. АН СССР. Отд. технических наук. – 1958. – № 8. – С. 26–31
  31. Работнов, Ю.Н. О механизме длительного разрушения // Вопросы прочности материалов и конструкций: сб. статей. Изд-во АН СССР. Москва. – 1959. – С. 5–7
  32. 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
  33. 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
  34. 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
  35. 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
  36. Barenblatt, G. I. Scaling, Self-similarity, and Intermediate Asymptotics. Cambridge University Press. 1996
  37. 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
  38. Баренблатт Г.И. Подобие, автомодельность, промежуточная асимптотика. Теория и приложения к геофизической гидродинамике. М.: Ленанд, 2021. 256 с
  39. Riedel H. Fracture at High Temperature. Berlin: Springer. – 1987. – 418 p
  40. Астафьев, В.И., Радаев Ю.Н., Степанова Л.В. Нелинейная механика разрушения. – Самара: Издательство “Самарский университет”. – 2001. – 631 c
  41. 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
  42. Степанова, Л.В. Смешанное деформирование пластины с трещиной в условиях плоского напряженного состояния/ Л.В. Степанова, Е.М. Яковлева// Вестник Пермского национального исследовательского политехнического университета. Механика. – 2014. – № 3. – С. 129-162. doi: 10.15593/perm.mech/2014.3.08
  43. 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

Статистика

Просмотры

Аннотация - 403

PDF (Russian) - 50

Cited-By


PlumX


© Быкова Ю.С., Степанова Л.В., 2025

Creative Commons License
Эта статья доступна по лицензии Creative Commons Attribution-NonCommercial 4.0 International License.

Данный сайт использует cookie-файлы

Продолжая использовать наш сайт, вы даете согласие на обработку файлов cookie, которые обеспечивают правильную работу сайта.

О куки-файлах