Дослідження траєкторії тріщини змішаного режиму руйнування за допомогою неявної схеми

Abstract

Досліджується траєкторія похилої тріщини в плоскому елементі конструкції з використанням неявної схеми, оптимізованої для досягнення точності в аналізі змішаного режиму руйнування. Розроблена схема забезпечує точне виконання умови рівності нулю коефіцієнта інтенсивності зсувних напружень (КІН) на кожному кроці алгоритму, що дозволяє точно відслідковувати напрямок розвитку тріщини. Розроблено адаптивний алгоритм для генерації регулярної сітки в межах кругової зони, що оточує вершину тріщини, змодельованої кусково-лінійною кривою. Такий підхід до створення сітки дозволяє покращити точність обчислення КІН, що підтверджується шляхом порівняння з відомими з літератури розв’язками для руйнування нормальним відривом. Крім того, проаналізовано залежність зсувного КІН від кута повороту другого сегмента тріщини в умовах критичного стану. Результат аналізу дозволив спростити алгоритм неявної схеми. Для подальшої перевірки було проведено паралельне дослідження з використанням нерегулярної сітки, згенерованої вбудованими функціями MATLAB, що продемонструвало порівнянну точність у визначенні КІН для різних типів сітки. Проілюстровано розраховану траєкторію тріщини та розподіл поля напружень фон Мізеса навколо її вершини, що надає змогу скласти повне уявлення про поведінку напружень та вплив уточнення сітки на прогнозування траєкторії.This paper presents a study of crack trajectory determination for an inclined crack in a planar structural element using an implicit investigation scheme optimized for the accuracy of mixed mode fracture analysis. The developed scheme ensures that the condition that the shear stress intensity factor (SIF) is zero at each algorithmic step holds strictly, enabling accurate tracking of crack propagation direction. For this purpose, an adaptive algorithm was developed to generate a regular mesh in the circular zone around the crack tip, modeled as a piecewise linear curve. This meshing approach improves the accuracy of SIF computation, which is validated through comparison with established mode I solutions in the literature. In addition, the dependence of the shear SIF on the rotation angle of the second crack segment under critical state conditions is analyzed to simplify the implementation of the implicit algorithm. For further validation, a parallel study using an irregular mesh generated by MATLAB’s built-in meshing functions was conducted, which demonstrated comparable accuracy of SIF determination in different types of meshes. Finally, this paper shows the calculated crack trajectory and the von Mises stress field distribution around the crack tip for the first four propagation steps, which provides a comprehensive understanding of the stress behavior and the effect of mesh refinement on the trajectory prediction.