Сравнительный анализ методов вычисления коэффициентов интенсивности напряжений для анизотропных материалов при смешанной моде разрушения
Журнал: Вычислительная механика сплошных сред @journal-icmm
Статья в выпуске: 3 т.19, 2026 года.
Бесплатный доступ
Представлены результаты сравнения четырех методов вычисления коэффициентов интенсивности напряжений (КИН) для анизотропных материалов в задачах линейной механики разрушения при смешанной моде: метод экстраполяции перемещений, метод экстраполяции напряжений, метод эквивалентного объемного интегрирования при вычислении Jk-интегралов (инвариантных интегралов J1 и J2) и метод вычисления Jk-интегралов на основе интегрирования по контуру (по определению). Детально рассмотрена специфика анизотропии материалов с кубической симметрией, актуальная для монокристаллических сплавов на никелевой основе, из которых изготавливаются рабочие и направляющие лопатки современных газотурбинных установок. Для этого класса материалов с использованием формализма Лехницкого получены аналитические выражения для определения КИН на основе Jk-интегралов при произвольной ориентации трещины и смешанной моде разрушения. В качестве тестовой задачи для сравнительного анализа возможностей методов решается задача об одноосном растяжении анизотропной плоскости с наклонной трещиной в предположении плоского напряженного состояния, допускающая аналитическое решение. Оцениваются точность и сходимость обсуждаемых методов вычисления КИН в зависимости от дискретизации области (числа конечных элементов), ориентации трещины по отношению к направлению нагружения и степени анизотропии материала. Задача решена для трех характерных углов наклона трещины: 15°, когда доминирует I мода (отрыва) в сравнении с II модой (продольного сдвига); 45°, когда КИН равны между собой; 75°, когда II мода (продольного сдвига) доминирует над I модой (отрыва). Исследована чувствительность погрешностей KI и KII к погрешностям вычисления J1 и J2, позволяющая установить ограничения на области применимости методов определения КИН на основе Jk-интегралов.
Короткий адрес: https://sciup.org/143186429
IDS: 143186429 | УДК: 539.42 | DOI: 10.7242/1999-6691/2026.19.3.21
Comparative analysis of computation methods for mixed-mode stress intensity factors of anisotropic materials
The article discusses the results of comparison of four numerical methods used to calculate stress intensity factors (SIFs) for anisotropic materials in problems of mixed-mode linear fracture mechanics: (i) displacement extrapolation method, (ii) stress extrapolation method, (iii) equivalent domain integral (EDI) method for computation of Jk-integrals (invariant integrals J1 and J2), and (iv) calculation of Jk-integrals using contour integration (by definition). The specific features of anisotropy in materials with cubic symmetry, which is inherent in single-crystal nickel-based alloys commonly used in the manufacture of blades of modern gas turbine units, are considered in detail. For this class of materials with arbitrary crack orientation and a mixed fracture mode, the analytical expressions for determining the SIF in terms of Jk-integrals are obtained using Lekhnitskii's formalism. The problem of uniaxial tension of an anisotropic plane with an inclined crack under the assumption of a plane stress state and admissibility of an analytical solution is considered as a test problem for a comparative analysis of the methods. The accuracy and convergence of the four considered methods for calculating the SIF are analyzed depending on the degree of domain discretization (number of finite elements), crack orientation with respect to the loading direction, and the degree of material anisotropy. The comparison is performed for three characteristic crack inclination angles: 15°, when the I mode (tension) dominates over the II mode (longitudinal shear); 45°, when the SIFs are equal to each other; 75°, when the II mode (longitudinal shear) dominates over the I mode (tension). The sensitivity of the errors in KI and KII to the errors in J1 and J2 is studied. This makes it possible to set limits on the scope of application of the methods for determining the SIF based on Jk-integrals.
