ЭКСПЕРИМЕНТАЛЬНОЕ ОПРЕДЕЛЕНИЕ ПОЛЕЙ У ВЕРШИНЫ ТРЕЩИНЫ: МЕТОД ГОЛОГРАФИЧЕСКОЙ ИНТЕРФЕРОМЕТРИИ И МЕТОД КОРРЕЛЯЦИИ ЦИФРОВЫХ ИЗОБРАЖЕНИЙ
- Авторы: Анисимов Г.С1, Степанова Л.В1
- Учреждения:
- Самарский национальный исследовательский университет имени академика С.П. Королева, Самара, Российская Федерация
- Выпуск: № 3 (2024)
- Страницы: 39-56
- Раздел: Статьи
- URL: https://ered.pstu.ru/index.php/mechanics/article/view/4292
- DOI: https://doi.org/10.15593/perm.mech/2024.3.03
- Цитировать
Аннотация
Рассматриваются вопросы реконструкции поля напряжений у вершин системы горизонтальных и наклонных трещин в линейно упругой среде с помощью результатов экспериментальных исследований, проведенных интерференционно-оптическими методами: методом голографической интерферометрии и методом корреляции цифровых изображений. Эксперимент нацелен на построение многокомпонентного асимптотического разложения М. Уильямса с удержанием регулярных (неособых) слагаемых ряда для пластин, ослабленных двумя взаимодействующими трещинами. Для восстановления коэффициентов ряда М. Уильямса были использованы картины абсолютной разности хода (изодромы), позволяющие с помощью закона Фавра определить главные напряжения в окрестности вершины трещины. С помощью метода корреляции цифровых изображений для ряда конфигураций образцов с трещинами определены поля деформаций. Экспериментальная информация, полученная двумя интерференционно-оптическими методами, использована для вычисления коэффициентов многопараметрического асимптотического разложения М. Уильямса. Предложена новая вариация переопределенного метода, ориентированная на линеаризованный закон Фавра, и позволяющая отыскать посредством итерационной процедуры коэффициенты асимптотического ряда М. Уильямса (обобщенные коэффициенты интенсивности напряжений), опираясь на результаты поляризационно-оптических измерений. Для верификации результатов обработки всех совокупности экспериментальных данных дополнительно проведен вычислительный эксперимент с помощью метода конечных элементов, что позволило вычислить обобщенные коэффициенты интенсивности напряжений, опираясь на поля напряжений, найденные методом конечных элементов. Предложена модификация переопределенного метода, основанная на применении исключительно полей напряжений, ассоциированных с вершиной трещины, определяемых из конечно-элементного решения. Показано, что коэффициенты ряда М. Уильямса, определенные с помощью натурного и вычислительного эксперимента, хорошо согласуются между собой.
Полный текст
Оценка структурной целостности, прочности, ресурса и живучести элементов конструкций инженерных сооружений является важной частью философии возможного использования конструкций с дефектами, размеры которых не превышают своих критических значений [1-4]. Коэффициенты интенсивности напряжений, используемые для характеристики интенсивности сингулярного упругого поля и полей напряжений вблизи вершины трещины, поддаются аналитическому вычислению для достаточно простых геометрий при определенных системах нагружения [5,6]. В силу указанной причины для оценки коэффициентов интенсивности напряжений и Т-напряжений применяют экспериментальные подходы и вычислительные технологии [7-18]. Несмотря на то, что определение коэффициентов интенсивности напряжений, характеризующего поля, ассоциированные с вершиной трещины в линейно упругом теле, представляет собой самостоятельную задачу [7-10], в настоящее время, выстроилось уверенное понимание необходимости сохранения регулярных слагаемых в многокомпонентном разложении М. Уильямса с коэффициентами, получившими название обобщенных коэффициентов интенсивности напряжений [11-18]. Нахождение обобщенных коэффициентов интенсивности напряжений с помощью теоретических подходов возможно и, как правило, опирается на решения теории функций комплексного переменного [5,6]. Идентификация коэффициентов интенсивности напряжений, Т-напряжений и коэффициентов неособых (регулярных) слагаемых ряда М. Уильямса (обобщенных коэффициентов интенсивности напряжений) по экспериментальным данным обеспечивает надежную основу для описания стационарного состояния и распространения трещин [4-11]. Хотя коэффициенты интенсивности напряжений являются наиболее важными параметрами механики разрушения, многие исследования показали [12-18], что члены разложения Уильямса более высокого порядка также могут вносить значительный вклад в поля напряжений и деформаций вблизи вершины трещины. В последние десятилетия развитие оптических методов измерения напряжений, деформаций или смещений позволило решить проблемы реконструкции многокомпонентного ряда Уильямса, описывающего поля у вершины трещины, и прояснить некоторые нерешенные или, по крайней мере, неясные вопросы [14,19-23]. Следует подчеркнуть, что для экспериментального определения полей, ассоциированных с вершиной надреза или трещины, традиционно используются интерференционно-оптические методы, которые позволили получить фундаментальные результаты в механике хрупкого разрушения [21-26]. Метод фотоупругости, трансформировавшийся в последние десятилетия в метод цифровой фотоупругости, долгое время применялся и применяется в настоящее время [27,28] для экспериментального вычисления коэффициентов интенсивности напряжений и Т-напряжений, а также коэффициентов при регулярных слагаемых ряда М, Уильямса [21-26]. Однако, в самое последнее время фокус внимания исследователей смещается в сторону реконструкции интерференционных картин с помощью численных методов и подходов искусственного интеллекта. Например, в [27] отмечается, что численный анализ сложных задач неизменно требует использования исходных данных из экспериментальных методов, таких как фотоупругость, для проверки граничных условий, а также численного моделирования. Авторы подчеркивают, что в нашу эпоху, когда инженерное сообщество придерживается культуры открытого исходного кода для решения численных задач, важным представляется возможность построения фотоупругих контуров по результатам конечно-элементного анализа. В [27] предлагается методика построения картин изохроматических полос, отточенная на эталонной задаче о диске при диаметральном сжатии. Результаты восстановленной интерференционной картины сопоставлены с лабораторными экспериментальными результатами и далее развитая методология применена к пластине с краевой трещиной. В целом, школа К.Рамеша в области фотоупругого анализа достигла убедительных успехов, пройдя путь от теоретических работ до внедрения своих достижений в практической области. В апреле 2023 года Индийский технологический институт Мадраса с помощью отраслевого партнера внедрил в коммерческую эксплуатацию четыре самых современных пакета программного обеспечения в области фотоупругого анализа и моделирования. Программное обеспечение может найти и находит применение в самых разных областях: от сельского хозяйства до движения биологических микроорганизмов, выявление зон концентрации напряжений или выявления дефектов в 3D-электронике. Новые сферы применения фотоупругости включают такие области, как биомедицинский и традиционный анализ напряжений, включающий сложные нагрузки и граничные условия, а также аэрокосмическую, авиационную, автомобильную и обрабатывающую отрасли промышленности. К. Рамеш отмечает [21,22,26-28], что анализ фотоупругих напряжений претерпел кардинальные изменения с появлением цифровых компьютеров и технологий обработки изображений. До сих пор не существовало комплексной программной платформы для внедрения подобных разработок в исследовательских лабораториях и отраслях промышленности, использующих этот метод. Авторы разработали первое комплексное программное обеспечение для экспериментального анализа фотоупругих напряжений. Врачи, агрономы и биологи в настоящее время все чаще обращаются к использованию фотоупругости для решения своих предметных задач. Они могут достоверно обрабатывать записанные изображения с помощью программного обеспечения для получения важных выводов из своих исследований. В [29] предлагается эволюционное развитие метода фотоупругости для измерения напряжений, основанного на применении кольцевых фотоупругих датчиков, а также регистрации картины напряжений с помощью цифровой камеры и ее распознавания с помощью искусственных нейронных сетей. Представлен анализ современного применения метода фотоупругости для решения различных задач в рамках теории прочности. Рассмотрен принцип работы фотоупругих датчиков, основанный на эффекте фотоупругости. Представлены оптические схемы в кольцевом фотоупругом датчике для различных значений компонент напряжения. В результате лабораторного эксперимента был собран набор данных из 1500 изохроматических изображений. В качестве алгоритма машинного обучения был применен подвид нейронной сети, а именно сверточная нейронная сеть. Были использованы различные комбинации моделей и оптимизаторов. Предложено применение скважинных датчиков для непрерывного мониторинга изменений напряженного состояния массива горных пород и интеграция этих данных в цифровую модель месторождения. В [30] авторы отмечают, что фотоупругость широко используется для анализа изменений свойств объектов, в том числе и объектов с микроразмерами, под действием сложных систем сил. При прогнозирующем обслуживании, основанном на машинном обучении, прогнозирование отказов оборудования осуществляется с использованием данных датчиков и системных сигналов. Однако о прогнозируемом обслуживании с использованием фотоупругости и машинного обучения пока не сообщается из-за трудностей с получением фотоупругих изображений с течением времени. Прогнозирование будущего фотоупругого изображения с использованием фотоупругих изображений может облегчить выявление неисправностей оборудования в ближайшем будущем, чтобы можно было своевременно принять превентивные меры для механических или структурных компонентов. Такое прогнозирование представляет собой пространственно-временную проблему, которая требует значительных повторяющихся данных для временных соотношений и обработки пространственной информации на основе изображений. В [30] впервые сгенерированы наборы данных синтетических фотоупругих изображений с временной обработкой. Было исследовано несколько комбинаций моделей глубокого обучения, включая автокодирование, рекуррентную нейронную сеть и сверточную нейронную сеть. Наконец, в [30] предлагается новая многоступенчатая модель сверточного автокодирования для прогнозирования сбоев оборудования. Эффективность предложенной модели с точки зрения потерь и точности была проверена с помощью численных тестов. Результаты показали, что предложенная модель является эффективной основой для прогнозирования будущего поведения микро-/наноразмерных материалов на основе изображений. Кроме того, предложенная модель может быть распространена на другие приложения, которым требуются пространственно-временные данные на основе интерференционных изображений. Как отмечается в [31], количественная оценка поля напряжений, создаваемого в изделии при его нагружении, важна для инженерных областей, поскольку позволяет охарактеризовать механическое поведение и отказы элементов конструкций, вызванные действием сложных систем термомеханических воздействий. Для решения этой задачи долгое время использовалась цифровая фотоупругость [21-28] в силу ее визуальной способностью представлять информацию о напряжении с помощью изображений изохроматических полос. К сожалению, демодуляция таких полос остается сложным процессом, который в некоторых случаях зависит от нескольких факторов, например, попиксельного сравнения, динамических условий загрузки приложений, исправления несоответствий, зависимости пользователей, процессов развертывания полос и т.д. [31]. Учитывая эти недостатки и используя преимущества мощных результатов, полученных в ходе машинного глубокого обучения, таких как процесс разворачивания изохроматических полос, авторы в [31] предлагают глубокую сверточную нейронную сеть для восстановления поля напряжений, представленного в виде интерференционной картины изохроматических полос, полученных с помощью цифровых исследований фотоупругости. Развитая модель [31] основана на необученной сверточной нейронной сети для точной демодуляции карт напряжений путем ввода только одного изображения фотоупругости. В [31] продемонстрировано, что предложенный метод точно восстанавливает поле напряжений сложных распределений полос на смоделированных изображениях. Работа [32] посвящена процедуре метода фотоупругости и получению качественных изображений. При получении экспериментальных изображений фотоупругости спектральное взаимодействие между источником света и используемым датчиком влияет на визуальную информацию об изохроматической картине полос на полученных изображениях. Такое взаимодействие может привести к появлению интерференционных полос с общим высоким качеством, но также может привести к получению изображений с неразличимыми полосами и плохой реконструкции поля напряжений. В [32] представлена стратегия оценки такого взаимодействия, основанная на измерении значения четырех дескрипторов, созданных вручную: контраста, дескриптора изображения, который одновременно учитывает размытость и шум, дескриптора на основе Фурье для измерения качества изображения и энтропии изображения. Полезность предложенной стратегии была подтверждена путем измерения выбранных дескрипторов на вычислительных изображениях фотоупругости и порядков полос, полученных при оценке поля напряжений, из 240 спектральных конфигураций: 24 источников света и 10 датчиков. Было обнаружено, что высокие значения выбранных дескрипторов могут быть связаны со спектральными конфигурациями, которые приводят к лучшей реконструкции поля напряжений. В целом, результаты показывают, что выбранные дескрипторы могут быть полезны для идентификации плохих и хороших спектральных взаимодействий, что могло бы помочь разработать более совершенные протоколы для получения изображений фотоупругости. Возвращаясь к теме машинного обучения необходимо отметить работу [33], в которой снова подчеркивается, что оценка поля напряжений на основе фотоупругости имеет жизненно важное значение в инженерных областях. Для достижения цели эффективной демодуляции распределения напряжений и преодоления ограничений традиционных методов важно разработать метод глубокого обучения, который упростит и ускорит процесс получения и обработки изображений. В [33] Предлагается подход, нацеленный на повышение точности прогнозирования. Авторы используют Resnet в качестве основы и, применив архитектуру U-Net и добавив модуль физических ограничений, разрабатывают модель, обеспечивающую восстановление поля напряжений с более высоким структурным сходством. Показано, что для различных условий нагружения и геометрий рассматриваемых экспериментальных образцов модель работает надежно, несмотря на сложную геометрию и большой диапазон напряжений. Результаты доказывают универсальность и эффективность нашей модели и дают возможность мгновенного определения напряжения. Таким образом, в современной экспериментальной механике техника цифровой фотоупругости с появлением новых возможностей обработки изображений остается рабочим инструментом исследователя. Наряду с цифровой фотоупругостью в современной механике разрушения для экспериментального исследования полей, ассоциированных с вершиной трещины или выреза, наиболее часто прибегают к методу корреляции цифровых изображений (DIC) [14,33], который позволяет по экспериментальным полям деформаций восстановить полное разложение М. Уильямса. Например, в [14] было проведено численное и экспериментальное исследование характеристик полей у вершин криволинейных трещин. Это исследование сочетает в себе новый подход дифференциальной геометрии в сочетании с результатами 3D-DIC и КЭА для определения параметров механики трещин по полю смещения на изогнутом элементе с трещиной. Работа направлена на проверку и сравнение некоторых аспектов, которые не могут быть решены с помощью 3D-DIC. Значения коэффициентов интенсивности напряжений для различных длин трещин определяются с использованием данных DIC и КЭА и сравниваются. Кроме того, проводится всестороннее сравнение экспериментальных и численных полей с использованием метода декомпозиции изображений, основанного на полиномах Чебышева. В [33] метод цифровой корреляции изображений используется для получения коэффициентов членов высшего порядка в разложении Уильямса в образцах компактного растяжения, изготовленных из полиметилметакрилата (ПММА). Поле смещения определяется корреляцией между эталонным изображением (т.е. до деформации) и деформированным изображением. Часть поля перемещений, возникающих в результате движения и вращения твердого тела, исключается из поля перемещений. Для большого числа точек вблизи вершины трещины собирается переопределенный набор одновременных линейных уравнений, и с использованием фундаментальных концепций метода наименьших квадратов рассчитываются коэффициенты разложения Уильямса для условий чистого режима I. Затем экспериментальные результаты сравниваются с численными результатами, рассчитанными методом конечных элементов (МКЭ). Показано, что существует очень хорошее соответствие между результатами DIC и КЭА, подтверждающее эффективность метода DIC в получении коэффициентов членов более высокого порядка разложения в ряд Уильямса по полю смещений вокруг вершины трещины. Несмотря на широкое использование экспериментальных техник для реконструкции асимптотического поля напряжений у вершины трещины, многие проблемы остаются открытыми. Например, определение коэффициента интенсивности напряжений является тривиальной задачей, если экспериментатор может быть уверен, что все данные о поле смещений находятся в пределах зоны с преобладанием сингулярности, окружающей вершину трещины. К сожалению, это редко случается из-за двух противоречивых ограничений. Во-первых, практические ограничения на размер образца уменьшают размер зоны, в которой доминирует сингулярность, до неудобно маленькой области, которая может иметь неизвестные размер и форму. Во-вторых, экспериментатору не следует брать данные слишком близко к вершине трещины из-за наличия нелинейной зоны. Это ограничение является наиболее строгим, когда данные берутся с поверхности образца, поскольку экспериментатор должен находиться на расстоянии не менее половины толщины пластины от вершины трещины, чтобы гарантировать, что смещения поверхности соответствуют смещениям по всей толщине образца. Описанная проблема приводит к необходимости построения многокомпонентного разложения и устойчивой процедуры извлечения обобщенных коэффициентов интенсивности напряжений из экспериментальных техник и вычислительных расчетов. Второй проблемой является задача привлечения нескольких экспериментальных методов для валидации результатов обработки экспериментальных данных. В этом отношении является возможность применения двух методов, обеспечивающих одновременное нахождение и полей напряжений, и полей смещений. Третьем направлением поиска, на наш взгляд, может стать использование метода голографической интерференции, который в силу более сложной оптической схемы, более сложной экспериментальной процедурой и более сложной обработкой данных, менее часто применяется в механике разрушения, несмотря на то, что закон Фавра, на котором базируется голографическая фотоупругость обеспечивает возможность получение систем уравнений для разделения главных напряжений, без привлечения дополнительных методов разделения напряжений. Ответам на сформулированные вопросы и посвящена настоящая работа.Об авторах
Г. С Анисимов
Самарский национальный исследовательский университет имени академика С.П. Королева, Самара, Российская Федерация
Л. В Степанова
Самарский национальный исследовательский университет имени академика С.П. Королева, Самара, Российская Федерация
Список литературы
- Матвиенко Ю.Г. Основы физики и механики разрушения. – М.: Физматлит, 2022. – 144 с
- Матвиенко Ю.Г. Двухпараметрическая механика разрушения. – М.: Физматлит, 2021. – 208 с
- Ritchie R.O., Liu D. Introduction to fracture mechanics. – Amsterdam: Elsevier, 2021. – 168 p
- Saxena A. Basic Fracture Mechanics and its Applications. – Boca Raton: CRC Press, 2022. – 342 p
- Hello G. Derivation of complete crack-tip stress expansions from Westergaard-Sanford solutions// International Journal of Solids and Structures. – 2018. – Vol. 144-145. – P. 265-275
- Hello G., Tahar, M.B., Roeland, J.-M. Analytical determination of coefficients in crack-tip stress expansions for a finite crack in an infinite plane medium// International Journal of Solids and Structures. – 2012. – Vol. 49. – P. 556-566.
- Dugnani R., Ma L. Stress intensity factor by combined perturbation method and universal weight functions// Theoretical and Applied Fracture Mechanics. – 2023. – Vol. 126. – 103994.
- Zhang X., Zhao T., Liu Y., Chen Q., Wang Z., Wang Z. A data-driven model for predicting the mixed-mode stress intensity factors of a crack in composite// Engineering Fracture Mechanics. – 2023. – Vol. 288. – 109385
- Ru M., Liu C., Wei Y. On enrichment strategies and methods to extract stress intensity factors using extended finite element method for biomaterials// Engineering Fracture Mechanics. – 2023. – Vol.281. – 109060
- Rao Q-h., Zhao C-c., Yi W., Sun D., Liu Z. Stress intensity factor calculation of the cracks interacted by the oval-holes in anisotropic elastic solids under remote and non-uniform surface stresses// Theoretical and Applied Fracture Mechanics. – 2022. – Vol. 121. – 103475.
- Profant T., Sladek J., Sladek V., Kotoul M. Assessment of amplitude factors of asymptotic expansion at the crack tip in flexoelectric solid under mode I and II loadings// International Journal of Solids and Structures. – 2023. – Vol. 269. – 112194
- Nejati M., Gholi S., Ayatollahi M.R. Crack tip asymptotic fields in anisotropic planes: Importance of higher order terms// Applied Mathematical Modelling. – 2021. – Vol. 91. – P. 837-862
- Melching D., Breitbarth E. Advanced crack tip characterization using conjugate work integrals// International Journal of Fatigue. – 2023. – Vol. 169. – 107501
- Camacho-Reyes A., Antunes F.V., Vasco-Olmo J.M., Diaz F.A., Felipe-Sese L.A. On the experimental characterization of the crack tip displacement fields on nonplanar elements: Numerical and experimental analysis// Theoretical and Applied Fracture Mechanics. – 2023. – Vol. 124. – 103803
- Stepanova L.V., Belova O.N. Coefficients of the Williams power expansion of the near crack tip stress field in continuum linear elastic fracture mechanics at the nanoscale // Theoretical and Applied Fracture Mechanics. – 2022. – V. 119. – 103298.
- Stepanova L.V., Belova O.N. Stress intensity factors, T-stresses and higher order coefficients of the Williams series expansion and their evaluation through molecular dynamics simulations// Mechanics of Advanced Materials and Structures. – 2022. – Vol. 30. – Issue 19. – P. 3862-3884
- Степанова Л.В., Семенов Д.А., Анисимов Г.С. Применение метода голографической интерферометрии для реконструкции ряда М. Уильямса у вершины трещины //Вестник Самарского университета. Естественнонаучная серия. 2023. – Т. 29. – No 1. – С. 15–46.
- Степанова Л. В., Белова О.Н. Идентификация коэффициентов интенсивности напряжений, T-напряжений и коэффициентов регулярных слагаемых высокого порядка в разложении Уильямса с помощью молекулярно-динамического моделирования// Вестник Пермского национального исследовательского политехнического университета. Механика. – 2023. – № 2. – С. 47-77.
- Alavi S.K., Ayatollahi M.R., Daneshfar M., Bahrami B. Experimental stress determination of blunt notches under combinations of modes I and II loading// Engineering Structures. – 2023. – Vol. 278. – 115517
- Shuai J., Zhao J., Lei L. Simple crack tip and stress intensity factor determination method for mode I crack using digital image correlation// nanoscale // Theoretical and Applied Fracture Mechanics. – 2022. – Vol. 122. – 103621
- Ramesh K. Digital Photoelasticity. Advanced Techniques and Applications. – Sham: Springer, 2013. – 430 p
- Ramesh K., Sasikumar S. Digital photoelasticity: Recent developments and diverse applications// Optics and Lasers in Engineering. – 2020. – 106186
- Разумовский И..А. Интерференционно-оптические методы механики деформируемого твердого тела. – М.: Издательство МГТУ им. Н.Э. Баумана, 2007. – 240 с
- Vivekanandan A., Ramesh K. Study of crack interaction effects under thermal loading by digital photoelasticity and finite elements// Experimental Mechanics. – 2020. – Vol. 60(3). – P. 295-316
- Vivekanandan A., Ramesh K. Study of interaction effects of asymmetric cracks under biaxial loading using digital photoelasticity// Theoretical and Applied Fracture Mechanics. – 2019. – V. 99. – P. 104-117
- Ramesh K. Developments in Photoelasticity-A renaissance. – IOP Publishing, 2021. – 225 p
- Anand C., Natarajan S., Ramesh K. Simulating Isochromatic Fringes from Finite Element Results of FEniCS// Experimental Techniques. 2023.
- Sasikumar S., Ramesh K. Framework to select refining parameters in Total fringe order photoelastcity (TFP)// Optics and Lasers in Engineering. – 2023. – Vol. 160. – 107277.
- Konurin A., Neverov S., Neverov A., Orlov D., Zharov I., Konurina M. Application of artificial neural networks for state analysis based on the photoelastic method// Geohazard Mechanics. – 2023.
- Lee H., An H., Lee D.-W. Time-staged photoelastic image prediction using multi-stage convolutional autoencoders// Engineering Applications of Artificial Intelligence. – 2022. – Vol. 116. – 105265.
- Brinez-de Leon J.C., Rico-Garcia M., Restrepo-Martinez A. PhotoelastNet : a deep convolutional neural network for evaluating the stress field by using a single color photoelasticity image// Applied Optics. – Vol. – Issue 7. – P. D50-62
- Fandino-Toro H., Aristizabal-Lopez Y., Restrepo-Martinez A., Brinez-de Leon J. Fringe pattern analysis to evaluate light sources and sensors in digital photoelasticity// Applied Optics. – 2023. – Vol. 62. – Issue 10. – P. D2560-2568.
- Zhao W., Zhang G., Li J. Accuracy improvement of demodulating the stress field with StressUnet in photoelasticity// Applied Optics. – 2022. – Vol. 61. – P. 8678-8687
- Ayatollahi M.R., Moazzami M. Digital image correlation methos for calculating coefficients of Williams expansion in compact tension specimen// Optics and Lasers in Engineering. – 2017. – Vol. 90. – P. 26-33
- Favre H. Sur une novelle methode optic de determination des tensions interieres// Revue d’Optic. – 1929. – Vol. 8(193). – P. 241-289
- Nocedal J., Wright., S.J. Numerical Optimization. – Cham: Springer, 2006. – 686 p
- Bahrami B., Talebi H., Ayatollahi M.R., Khosravani M.R. Artificial neural network in prediction of mixed-mode I/II fracture load// International Journal of Mechanical Sciences. – 2023. – Vol. 248. – 10821