CC BY
179
Любутин П.С., Панин С.В., Титков В.В., Еремин А.В., Сундер Р. Развитие метода корреляции цифровых изображений для изучения процессов деформации и разрушения конструкционных материалов // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2019. - № 1. - С. 87-107. DOI: 10.15593/perm.mech/2019.1.08
Lyubutin P.S., Panin S.V., Titkov V.V., Eremin A.V., Sunder R. Development of digital image correlation method to study deformation and fracture processes of structural materials. PNRPUMechanics Bulletin, 2019, no. 1, pp. 87-107. DOI: 10.15593/perm.mech/2019.1.08
ВЕСТНИК ПНИПУ. МЕХАНИКА
№ 1, 2019 PNRPU MECHANICS BULLETIN
http://vestnik.pstu.ru/mechanics/about/inf/
Б01: 10.15593/регш.шеоЬ/2019.1.08 УДК 681.3
РАЗВИТИЕ МЕТОДА КОРРЕЛЯЦИИ ЦИФРОВЫХ ИЗОБРАЖЕНИЙ ДЛЯ ИЗУЧЕНИЯ ПРОЦЕССОВ ДЕФОРМАЦИИ И РАЗРУШЕНИЯ КОНСТРУКЦИОННЫХ МАТЕРИАЛОВ
П.С. Любутин1, С.В.Панин12, В.В. Титков1, А.В. Еремин1, Р. Сундер2
1Институт физики прочности и материаловедения СО РАН, Томск, Россия Национальный исследовательский Томский политехнический университет, Томск, Россия
О СТАТЬЕ АННОТАЦИЯ
Описан ряд задач, посвященных развитию метода корреляции цифровых изображений, как с позиции разработки и тестирования аппаратно-программных средств, так и решения проблем в области механики усталостного разрушения. Проведено тестирование разработанной компьютерной системы стереозрения на сериях стереопар, отражающих изменение положения объекта в пространстве, а также плоскостные и внеплоскостные деформации. Показано, что ошибка определения пространственных координат не превышает 0,75 единиц, а ошибка вычисления компонент тензора деформации в случае системы с одной камерой на два порядка больше, чем при использовании системы стереозрения. Предложен и протестирован алгоритм автоматического детектирования трещины на оптических изображениях и определения координат ее вершины. Показано, что при размере кадра 2000*1000 элементов он позволяет определить координаты вершины трещины со средней ошибкой порядка 56 пикселов, а средняя ошибка определения площади трещины не превышает 1,93 %. Предложен модифицированный инкрементальный алгоритм определения перемещений на серии стереопар, который позволяет оценивать перемещения большой величины при последовательной обработке изображений. Развит алгоритм измерения J-интеграла с использованием метода корреляции цифровых изображений. Показано, что отклонение рассчитываемых значений J-интеграла от модельных составляет в среднем 1,75 %. Проведена количественная характеризация процесса роста усталостной трещины в металлических сплавах с помощью методики, основанной на методе корреляции цифровых изображений (в терминах механики разрушения), включающей расчет скорости роста усталостной трещины da/dN; максимального значения деформации (£max) и эффективной асимметрии цикла (Reff). При испытании алюминиевого сплава Д16АТ показано, что в режиме Пэриса при единичной перегрузке минимальная скорость роста трещины достигается при ее удлинении на 30 % от размера пластической зоны, сформированной циклом перегрузки, и вызвана закрытием трещины, увеличивающим эффективную асимметрию цикла до Reff = 0,32.
© ПНИПУ
Получена: 03 декабря 2018 г. Принята: 11 марта 2019 г. Опубликована: 30 марта 2019 г.
Ключевые слова:
корреляция цифровых изображений, механика разрушения, усталость, рост трещины, оценка деформации, ^интеграл, стереозрение, инкрементальный метод.
© Любутин Павел Степанович - к.т.н., н.с., e-mail: [email protected], : 0000-0003-1732-0791 Панин Сергей Викторович - д.т.н., проф., зав. лаб., e-mail: [email protected], ¡D: 0000-0001 -7623-7360 Титков Владимир Викторович - к.т.н., м.н.с., e-mail: [email protected], ¡D: 0000-0002-5064-8142 Еремин Александр Вячеславович - к.т.н., м.н.с., e-mail: [email protected], ¡D: 0000-0002-1800-9818 Сундер Рамасуббу- к.т.н., доц., e-mail: [email protected], : 0000-0002-5339-0132
Pavel S. Lyubutin - CSc in Technical Sciences, Researcher о, e-mail: [email protected], : 0000-0003-1732-0791 Sergey V. Panin - Doctor in Technical Sciences, Professor, Head of Laboratory, e-mail: [email protected], j D: 0000-0001 -7623-7360
Vladimir V. Titkov - CSc in Technical Sciences, Junior Researcher, e-mail: [email protected], ¡D: 0000-0002-5064-8142
Alexander V. Eremin - CSc in Technical Sciences, Junior Researcher, e-mail: [email protected], : 0000-0002-1800-9818 Sunder Ramasubbu - CSc in Technical Sciences, Associate Professor, e-mail: [email protected], ¡D: 0000-0002-5339-0132
Эта статья доступна в соответствии с условиями лицензии Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)
DEVELOPMENT OF THE DIGITAL IMAGE CORRELATION METHOD TO STUDY DEFORMATION AND FRACTURE PROCESSES OF STRUCTURAL MATERIALS
P.S. Lyubutin1, S.V. Panin12, V.V. Titkov1, A.V. Eremin1, R. Sunder2
Institute of Strength Physics and Materials Science, Siberian Branch of the Russian Academy of Sciences, Tomsk, Russian Federation
2National Research Tomsk Polytechnic University, Tomsk, Russian Federation
ARTICLE INFO
ABSTRACT
Received: 03 December 2018 Accepted: 11 March 2019 Published: 30 March 2019
Keywords:
digital image correlation, fracture mechanics, fatigue, crack growth, strain estimation, J-integral, stereo vision, incremental method.
A number of problems related to the digital image correlation method from the standpoint of both the hardware and software development and testing, as well as solving problems in the field of fatigue fracture mechanics are described. Testing of the developed computer stereo vision system was carried out using a series of stereo pairs reflecting the change in the location of the object in space as well as plane and out-of-plane deformations. It is shown that the error in determining spatial coordinates does not exceed 0.75 units, while the error in computing the strain tensor components in case of a system with a single camera is two orders of magnitude larger than that at using the stereo machine vision system. An algorithm for automatic crack detection on optical images and calculating its tip coordinates was proposed and tested. It is shown that when the frame size is 2000x1000 pixels the coordinates of the crack tip might be determined with an average error of about 56 pixels, while the average error of the crack area determination does not exceed 1.93%. A modified incremental algorithm for calculating displacements on a series of stereo pairs is proposed, which allows one to estimate large magnitude displacements during serial processing of images.
An algorithm for measuring the J-integral using the digital image correlation method has been developed. It is shown that the deviation of the calculated J-integral values from the model ones is on average 1.75 %. A quantitative characterization of the fatigue crack growth process in metal alloys was carried out using the technique based on the digital image correlation (in terms of fracture mechanics) including the calculation of the fatigue crack growth rate da/dN, maximum strain (£max) and effective cycle asymmetry (Reff).
©PNRPU
Введение
Современный уровень развития технологий автоматической обработки, анализа и интерпретации визуальной информации обусловливает их широкое применение во многих отраслях: микробиологии, медицине, анализе структуры (в том числе наноматериалов), автоматизированных системах для выполнения широкого ряда производственных операций и др. [1]. В рамках экспериментальной механики также активно развивается и применяется ряд подходов, основанных на обработке и анализе изображений, получаемых в процессе нагружения объектов исследований.
К настоящему времени в научно-технической литературе получил широкое распространение термин «оптический поток» («optical flow»): кажущееся движение яркостной картины, наблюдаемое при перемещении объектов перед камерой либо камеры в неподвижной окружающей обстановке. Работы по данному направлению активно развиваются с начала 80-х гг. XX в. Среди основных методов определения оптического потока выделяют: дифференциальные алгоритмы [2-7], фазовые (частотные) алгоритмы [8-12]; блочные алгоритмы [1316]. Первые дифференциальные алгоритмы были предложены B.K. P. Horn и B.G. Schunck [4], а также B.D. Lucas и T. Kanade [5, 6].
Развитие алгоритмов в первую очередь было основано на усложнении модели движения, которая учитывает афинные преобразования, поворот [17-19]. Проблемы фильтрации векторных полей были затронуты в работах [20, 21]. Проблему ошибок определения поля оптического потока часто решают с использованием пространственной регуляризации, которая впоследствии была расширена до пространственно-временной [22-26]. Алгоритмы вычисления оптического потока активно реализуются с использованием параллельных вычислений [27-32]. В то же время известен ряд исследований, посвященных непосредственно алгоритмической оптимизации вычислений при определении оптического потока [33-35]. Также в работах [36-38] было рассмотрено применение многомасштабного иерархического подхода для обработки изображений (пирамиды изображений) при вычислении оптического потока.
Также на расчете оптического потока основан и получивший широкое распространение в экспериментальной механике для оценки деформации твердых тел метод корреляции цифровых изображений (DIC - Digital Image Correlation). В основе его действия лежит построение векторов перемещений и дальнейшей расчет компонент деформации с помощью процедуры численного дифференцирования. В настоящее время метод корреляции цифровых изображений является одним из
наиболее распространенных подходов к изучению процессов деформации и разрушения структурно-неоднородных материалов. Существенный вклад в его развитие внесли M.A. Sutton [13, 39-43], F. Hild, S. Roux [4447], B. Pan [48-52] и др.
В последние годы аппаратно-программные комплексы, принцип действия которых основан на корреляции цифровых изображений, получили широкое коммерческое распространение, прежде всего продукция Correlated Solutions Inc. (VIC 3D), LaVision GmbH (StrainMaster) и др. Указанные системы, как правило, оборудуются штатными оптическими системами, калибровка которых проводится перед проведением эксперимента в зависимости от условий видео (фото) наблюдения. Исследования с помощью указанных методов проводятся и в Центре экспериментальной механики ПНИПУ [53].
В то же время компьютерный код, поставляемый в составе таких систем, является закрытым, что не позволяет эффективно решать ряд прикладных задач. Поэтому рядом коллективов проводятся исследования по развитию собственных программно-алгоритмических решений, а также и аппаратных средств для изучения процессов деформации и разрушения, включая расчет стандартных параметров механики разрушения на основе метода корреляции цифровых изображений [46, 54-62].
В данной работе описан ряд задач подобного класса, решенных в последнее время в лаборатории механики полимерных композиционных материалов ИФПМ СО РАН.
1. Развитие аппаратных средств метода корреляции цифровых изображений
В методе корреляции цифровых изображений (DIC) алгоритмы определения движения на серии изображений в общем виде основаны на минимизации функционала:
E i w, S ) = E i w, S ) + oEs i w, S ),
(1)
где w - функция формы; - область изображения, по которой минимизируется Е; Ев - мера подобия блоков изображений; Ех - мера подобия векторов в оптическом потоке; а - регуляризационный коэффициент.
Меры близости Ев двух участков текущего изображения (ТИ) и эталонного изображения (ЭИ) могут определяться различными функциями, условием которых является наличие экстремума, обычно минимума. Наиболее простой мерой близости является сумма квадратов разности [63, 64]:
SSD = Xi F i X )- G i X, w))2
(2)
где Е (X), О (X, w) - значения яркости пикселов сравниваемых участков ТИ и ЭИ соответственно; X - координаты пикселов участка 5; w - функция формы.
Одной из наиболее сложных мер является нормированная сумма квадратов разности с нулевым средним [65, 66]:
ZNSSD = X
F i X )-F
vx x .s i F i X )-F )2 G i X, w)-G
(3)
y/x x.s (g i X, w)-G )2
где Е и О - среднеарифметические значения участков ТИ и ЭИ соответственно.
Существует два пути нахождения минимума функционала Е.
1. Блочный метод. Область 5 начального изображения сравнивается в скользящем окне со вторым изображением. Часто данный подход используется для предварительного поиска больших перемещений с пиксельной точностью. Однако он редко применяется при вычислениях с субпиксельной точностью и со сложной функцией формы, так как приводит к неоправданно большим вычислительным затратам.
2. Применение оптимизации. Используется какой-либо численный оптимизационный алгоритм, например алгоритм Ньютона-Рафсона или Гаусса-Ньютона. Оптимизационный подход, напротив, используется в субпиксельной области, поскольку обеспечивает устойчивую сходимость в диапазоне [-1.. .+1] пикселей.
Оптимизацию методом Ньютона-Рафсона [42, 67] можно записать в следующем виде:
pk+1 = pk -
VE i P ) VVE i P ) '
(4)
где к - номер итерации; р - параметрический вектор, содержащий коэффициенты функции формы.
Система уравнений, описывающая перемещения в области 5, основанная на аффинной функции формы, имеет следующий вид [13, 67]:
х(х,у) = р0+р2х + р4у, у(х,у) = р1 +р3х + р5у.
(5)
Был реализован лабораторный макет системы технического стереозрения (СТСЗ), который включает две цифровые фотокамеры Canon EOS 700D (разрешение КМОП (CMOS) матриц 5184x3456 пикселов при физическом размере сенсора 22,3x14,9 мм), закрепленные на штанге с возможностью регулировки расстояния и угла между камерами. На рис. 1 приведен лабораторный макет системы технического стереозрения с камерами на расстоянии 40 см друг от друга и углом поворота каждой камеры относительно центральной оптической оси системы равным 10°. Камеры были подключены к ЭВМ с целью их синхронизации и получения изображений в память ЭВМ в процессе работы системы. Система
позволяет записывать пару изображений (от левой и правой камеры) в текущий момент времени и может регистрировать серию изображений до 5 кадров в секунду при уставленных параметрах фокуса и выдержки.
Разработанная система технического стереозрения для задач оценки деформации включает следующие этапы функционирования (рис. 2): калибровка камер,
включающая поиск узловых точек калибровочного шаблона на изображениях; ректификация изображений стереопар; вычисление карт диспаратности; восстановление трехмерных координат точек стереопары; вычисление трехмерных перемещений точек; вычисление деформаций в терминах экспериментальной механики.
Рис. 1. Лабораторный макет системы технического стереозрения Fig. 1. The laboratory scale stereo machine vision system
Рис. 2. Алгоритм работы системы технического стереозрения Fig. 2. The operation algorithm of the stereo machine vision system
Сопоставление методик оценки деформации при 2D- и 3D-расчете. Сопоставление результатов расчета деформации в режиме моно- и стереозрения проведено на модели изгиба тонкой пластины, опертой по краям.
Такая модель выбрана в первую очередь по причине присутствия в ней внеплоскостной деформации. Была получена последовательность модельных пар изображений, отражающих внеплоскостные перемещения.
Оценивали среднюю абсолютную ошибку вычисления продольного и поперечного тензора деформации как для случая плоско-деформированной поверхности, так и для случая внеплоскостных деформаций. Среднюю абсолютную ошибку оценивали по формуле
1 N
N 518' -
(6)
где 8,. - заданное значение тензора деформации; Е -рассчитанное значение тензора деформации. Ошибки поперечной и продольной компоненты тензора обозначались Е^ и Е соответственно.
Эксперимент показал, что система с одной камерой не позволяет достоверно определять деформацию при наличии внеплоскостных перемещений (рис. 3). Ошибка вычисления компонент тензора деформации в случае системы с одной камерой на два порядка больше, чем при использовании системы стереозрения. Кроме того, в случае системы с одной камерой ошибка на порядок больше самой величины деформации.
б
Рис. 3. Зависимость ошибки определения компонент тензора деформации Eve (а) и Esy (б) от номера стереопары для системы технического зрения с одной (2D) и двумя камерами (3D)
Fig. 3. Error graph of determining the Eve (a) and Esy (b) strain tensor components vs the stereopair number for a vision system with one (2D) and two cameras (3D)
Пример восстановленной поверхности представлен на рис. 4.
Рис. 4. Поле глубины, отображенное в координатах левого изображения при изгибе поверхности опертой по краям
Fig. 4. The depth field displayed in the coordinates of the left image under bending of a surface supported over the edges
Пример поля компоненты szz тензора деформации в случае изгиба поверхности, опертой по краям, представлен на рис. 5.
Рис. 5. Поле компоненты Sxx тензора деформации, отображенное в координатах левого изображения при изгибе поверхности, опертой по краям
Fig. 5. The field of the Sxx strain tensor component displayed in the coordinates of the left image under bending of the surface supported over the edges
2. Инкрементальный алгоритм оценки перемещений
Определение пространственных координат заданных точек объекта по последовательности стереопар подразумевает определение положения этих точек на изображениях. Для этого необходимо знать положения этих точек в начальный момент времени. Определение перемещения точек между начальным и текущим изображением может быть затруднено по причине значительной деформации объекта, изменения условий съемки и т.д. При изменении положения объекта в пространстве на проекциях происходят значительные изменения. Поскольку между соседними изображениями в серии такие изменения значительно меньше, чем между начальным и текущим изображениями, то помехоустойчивость и точность определения перемещений значительно возрастает. По этой причине был предло-
а
жен инкрементальный алгоритм, расширенный для обработки последовательности стереопар.
Используется сквозная последовательная нумерация изображений стереопар (рис. 6).
Номер стереопары
Изображения стереопар
Левая Правая
камера камера
10 11
12 13
I{2ii I{2i+1}
Рис. 6. Нумерация изображений стереопар Fig. 6. The principle of a stereopair image numbering
Схема инкрементального алгоритма приведена на рис. 7, в ней используются обозначения, описанные ниже:
1. I0, I1, I2 и т.д. изображения, как на рис. 6.
2. VF01P, VF01S - поля векторов перемещений, посчитанные с пиксельной и субпиксельной точностью между изображениями I0 и I1.
3. Для обобщения п. 1, 2 и нумерации изображений и векторных полей в зависимости от текущего номера стереопары введены обозначения с фигурными скобками вида I{a}, VF{a}{b}P, VF{a}{b}S, где {a}, {b} -номера изображений пары в общем виде, например I{2(i-1)}, VF{2(i-1)}{2i}P.
4. VFP, VFS - процедуры определения полей перемещений с пиксельной и субпиксельной точностью соответственно.
10 II ±
0-я стереопара
10 И
VFP VFOIP^ VFS
VF0IS
VF!2i!{2i+l|S б
Рис. 7. Схема инкрементального алгоритма определения полей векторов перемещений для последовательностей стереопар
Fig. 7. The scheme of the incremental algorithm for constructing the displacement vectors fields for stereopairs in image series
Стереопара, соответствующая начальному моменту времени, обозначена как 0-я стереопара, изображения 10 и 11 (рис. 7, а). На первом этапе строится поле векторов перемещений с пиксельной точностью, по классическому алгоритму, основанному на быстром преобразовании Фурье и обозначенному на схеме как УТР. Затем поле векторов уточняется до субпиксельной точности (блок УТ^) с использованием дифференциального алгоритма. Поле векторов перемещений УТШ8 между кадрами стереопары является картой диспаратности, а также используется как начальное приближение для определения перемещений в паре 11, 13 с пиксельной точностью. Для первой стереопары входное векторное поле в левой части схемы УГ008 соответствует нулевому полю и не используется.
Блоки УТР используют предыдущее и текущее изображения для определения поля векторов перемещений. Перемещения определяются отдельно для левых изображений 1(2(1-1)}, 1(21} и правых 1{2(1-1)+1}, 1(21+1} (рис. 7, б). Также используется векторное поле, рассчитанное для предыдущей стереопары, перемещения которого используются для определения положения ядра корреляции.
Перемещения с субпиксельной точностью (блоки УГ8) определяются по изображению 10 (начальное левое изображение) и текущему изображению стереопары, левому 1(21}, либо правому 1(21+1}. Для начального приближения двух первых элементов р0 и р\ искомого вектора Р
х(х,у)=р0+р2х + р4у, у{х,у) = р1 +р,х + р5у
(7)
используется векторное поле с пиксельной точностью УР(2(1-1)}(21}Р либо УР(2(1-1)+1 }(21+1 }Р (см. рис. 7, б).
Вычитая поле УР0(21}Б из УР0(21+1}8, получаем векторное поле УР(21}(21+1}8, которое соответствует карте диспратности 1-й стереопары изображений 1(21} и 1(21+1}.
Была сгенерирована серия модельных изображений, отражающих изменение рельефа, при одновременном растяжении поверхности материала. Для изменения характера рельефа каждое изображение формировалось из пары изображений двух модельных серий для двухосного растяжения при заданном приращении деформации. Построенные поля векторов перемещений (рис. 8, а, б) показывают значительное преимущество инкрементального алгоритма над обычным подходом. Зависимость коэффициента корреляции модельного и построенного векторного поля от номера кадра серии изображений (рис. 8, в) показывает, что при изменении рельефа поверхности более чем на 50 % от исходного «классический» алгоритм не справляется с задачей сопоставления изображений. Более того, при изменении рельефа на 75 % и более сопоставить начальное и текущее изображения не представляется возможным. Инкрементальный подход позволяет сопоставлять изображения во всем диапазоне приращения деформаций.
а
б
9 13 17 21 Номер кадра
в
Рис. 8. Векторные поля, полученные по «классическому» алгоритму (а) и инкрементальному алгоритму (б). Зависимость коэффициента корреляции kr векторного поля от номера кадра серии изображений (в): 1 - классический алгоритм;
2 - инкрементальный алгоритм
Fig. 8. Vector fields calculated by the classical (a) and incremental algorithms (b). Dependence of the correlation coefficient kr of a vector field vs the frame number of an image series (c): 1 is the classical algorithm, 2 is the incremental algorithm
3. Развитие алгоритма выделения трещины
В рамках механики усталостного разрушения одной из ключевых является задача оценки параметров распространения трещины (длины, скорости роста, раскрытия берегов и т.д.), поэтому актуальной технической задачей является автоматизация процесса их измерения, что обусловливает необходимость слежения за ее вершиной. С позиции оценки сопротивления материалов зарождению и росту усталостных трещин [68] значительный интерес представляют задачи, когда длина трещины не превышает десятков микрон. В этом случае применение макрорежима видеофиксации не позволяет достичь требуемого пространственного разрешения [69].
Для решения указанной проблемы используются оптические микроскопы, оснащенные фото- или видеокамерой. Однако при больших увеличениях оптической системы в процессе роста (усталостной) трещины ее вершина будет постоянно выходить из поля зрения, поэтому необходимо создавать аппаратно-программные комплексы, позволяющие (в режиме реального времени) управлять механическим сканером, перемещающим микроскоп и видеокамеру относительно образца, на поверхности которого контролируется рост трещины.
Ключевыми работами, направленными на автоматическое детектирование трещин на оптических изображениях, можно считать [70, 71]. За прошедшее время сформировалось направление прикладных исследований, нацеленных на автоматическое выделение трещин на изображениях. Значительный практический интерес представляло детектирование и количественная харак-теризация трещин на поверхности дорожного полотна [71-73], в то время как другая заметная часть работ посвящена оценке состояния инженерных конструкций и выявлению трещин в них [74-76]. В данном разделе работы решалась задача разработки и тестирования алгоритма детектирования трещины на оптических изображениях поверхности нагруженного материала и определения координат ее вершины.
Алгоритм выделения трещины. Алгоритм детектирования трещины на серии изображений и определения координат ее вершины содержит следующие основные этапы: 1) анализ последовательности изображений; 2) выделение объекта максимальной площади, соответствующего области наибольших неравномерностей на картине оптического потока; 3) анализ индивидуальных изображений; 4) построение скелетона трещины; 5) поиск вершины трещины.
Анализ последовательности изображений. При анализе серии изображений (либо видеопоследовательности), отображающих распространение трещины, ключевую информацию содержат данные о перемещениях участков (фрагментов) на них. В этом случае характер текстуры на поверхности важен лишь для повышения достоверности (помехоустойчивости и точности) определения перемещений. Для дальнейшего анализа используются данные об оптическом потоке: Е(/(, 11_1), где I и /(1 - изображения в текущий и предыдущий
моменты времени. Таким образом, объектом анализа является оптический поток, визуализация которого в виде цветового кодирования показана на рис. 9, а [77]. При этом разрыв или резкое изменение перемещений в его пределах (на берегах раскрывающейся трещины) свидетельствует о наличии трещины.
Для определения оптического потока используется ранее описанный в подразд. 1 корреляционный алгоритм. Однако с точки зрения детектирования трещины алгоритм также должен обеспечивать устойчивость при поиске перемещений на ее границах (берегах). Данное требование, предъявляемое к подобного рода алгоритмам, обосновано необходимостью точнее выделять область интереса (здесь и далее под этим термином подразумевается фрагмент изображения, содержащий трещину). Для этого могут быть использованы алгоритмы, основанные на расчете взвешенной меры подобия блоков изображений, поскольку они позволяют в полной мере «сохранять» границы между объектами [78].
а
в
Рис. 9. Результаты анализа последовательности изображений на основе определения поля перемещений: изображение оптического потока (а); изображение дифференцированного оптического потока (б); алгоритмически определенное расположение трещины (в)
Fig. 9. Results of the image series analysis based on the constructing the displacement fields: optical flow pattern (a); differentiated optical flow pattern (b); algorithmically determined crack location (c)
Поскольку у «коротких» трещин (длиной до десятков микрон) величина раскрытия может составлять единицы пикселов, с целью дополнительного контрастирования и подчеркивания границ применяется процедура дифференцирования поля оптического потока dF (1, ) (рис. 9, б), что реализуется использованием
оператора Собеля для окна размером 3*3 [79].
Выделение объекта максимальной площади. Для определения области интереса полученное изображение дифференцированного оптического потока подвергается нормализации (приведению к диапазону [0, 255]) и пороговой бинаризации. В качестве порога бинаризации принимается значение яркости, равное единице, т.е. «отбрасываются» только самые малые перемещения. Далее к бинаризованному изображению применяется операция морфологического замыкания [79] с целью заполнения мелких промежутков внутри сплошных объектов. Процесс завершается поиском контуров, оставшихся на изображении объектов, и оценкой их формы. В качестве трещины принимается объект наибольшей площади. Данный метод применяется для первичного (грубого) детектирования фрагмента изображения, наиболее вероятно содержащего трещину (рис. 9, в).
Анализ индивидуальных изображений. При обработке одиночных изображений (без учета остальных изображений в серии/видеопотоке) в основе анализа данных лежит информация о текстуре поверхности. В процессе слежения за трещиной решаются две ключевых задачи: ее поиск на изображении и последующее определение местоположения ее вершины (детектирование и локация). Для этого проводится контурный анализ, основанный на разбиении изображения на отдельные объекты и фон, т.е. процедура сегментации [79]. Использование контурного анализа позволяет выделить объекты, максимально подобные трещине.
Априори считается, что трещина (как область, соответствующая отсутствующему материалу - несплошности) должна являться наиболее темным участком (объектом) на изображении. С целью снижения вычислительных затрат и повышения помехоустойчивости работы алгоритма детектирования трещин используется процедура бинаризации. Это позволяет заметно сократить количество выделяемых и анализируемых объектов. Для бинаризации применяется алгоритм, в котором для определения глобального порогового значения яркости анализируется соответствующая гистограмма.
Последующей операцией в обработке полученного черно-белого изображения является морфологическое замыкание [79]. Его целью является удаление мелких (темных) фрагментов на изображении (формо- и разме-розависимая фильтрация), соответствующих элементам текстуры на поверхности, не относящихся к трещине. Затем полученное черно-белое изображение подвергается контурному анализу с целью детектирования объекта (объектов), геометрические характеристики которого максимально соответствуют таковым для трещины. В качестве трещины принимается объект максимальной площади.
Построение скелетона трещины. Для поиска и определения координат вершины несплошности используется построение скелетона [80] объекта типа «трещина». Для этого оценивается кратчайшее расстояние от точки внутри трещины до ее границ. Применение данной операции основано на том, что скелетом несплошности является кривая, объединяющая точки, являющиеся равноудаленными от ее границ.
Поиск вершины трещины. Для определения координат вершины несплошности проводится оценка величины раскрытия вдоль оси скелетона. Решение данной задачи упрощается тем, что все смещения уже были рассчитаны в первой части работы алгоритма. Искомая точка принадлежит участку, характеризуемому наименьшими смещениями.
Тестирование алгоритма. Тестирование алгоритма проводили с использованием сгенерированной серии модельных изображений. На первом этапе функционирование предложенного метода исследовали в условиях отсутствия возмущающих воздействий (шум, размытие, имитирующее расфокусировку, и др.). Серию модельных изображений формировали в несколько этапов в
виде текстурированной поверхности, имитирующей напыленный спекл [13, 81], на которую программным способом накладывалась трещина (рис. 10, а): 1) создание модели поверхности материала (формирование изображения спекла); 2) задание параметров несплошности (длина, раскрытие, форма берегов, имитирующая ее ветвление); 3) изменение модели материала (пересчет изображения) с учетом изменения параметров трещины.
б
Рис. 10. Модельное (а) и экспериментально полученное (б) изображения, характеризующие распространение усталостной трещины
Fig. 10. The model (a) and experimentally obtained image (b) characterizing the fatigue crack propagation
Для определения устойчивости работы алгоритма на втором этапе тестирования было сгенерировано 8 пар модельных изображений, отражающих 4 варианта траектории роста трещины и 2 типа текстуры. Каждая траектория получена случайным разбиением линии раскрытия трещины. Длина последней составляла половину ширины изображения, а величина раскрытия между двумя соседними изображениями в последовательности, снятыми в моменты времени t и t + 1, составляла 30 и 60 пикселов. Использованные варианты текстуры различались размером элементов напыленного спекла, когда максимальный диаметр «пятен» равнялся 20 и 30 пикселов. Сгенерированные модельные изображения подвергали размытию и зашумлению (для имитации искажений, характерных для лабораторных экспериментов).
Зашумление изображений проводили путем наложения на каждый пиксель (значения яркости которого лежат в диапазоне [0; 255]) случайного значения, имеющего нормальное распределение (что подобно характеру шумов для реальных видеодатчиков) и различное стандартное отклонение s = [0;45]. Для размытия ис-
пользовали нормальное Гауссово размытие со среднеквадратичным отклонением а = [0; 4,5].
Дополнительно проведены эксперименты по детектированию трещины при обработке экспериментальных изображений [69], зарегистрированных при циклическом растяжении образца алюминиевого сплава Д16 толщиной 3 мм (рис. 10, б). Для оценки эффективности работы предложенного алгоритма были использованы два критерия:
1. Ошибка определения местоположения вершины трещины, равная расстоянию между ее координатами на модельном изображении (с их известным значением) и таковыми, определенными с помощью предложенного алгоритма:
(8)
где (С, Суг) и (С4, С^) - модельные и найденные
координаты соответственно. Заметим, что причиной данной ошибки в первую очередь может являться потеря информации в области вершины трещины, когда алгоритм в процессе бинаризации «отсекает» фрагмент, также относящийся к несплошности.
2. Ошибка определения площади трещины, выделенной с помощью предложенного алгоритма (см. выделение объекта максимальной площади), рассчитываемая из отношения найденной площади (5/) к ее величине, заданной в модели:
dES =
S
-l
100%.
(9)
О повышении точности функционирования алгоритма должно свидетельствовать достижение минимального значения каждого из этих двух показателей.
Результаты тестирования и их обсуждение. На рис. 11 представлены графики зависимостей ошибок определения параметров трещины от степени зашумле-ния (5) (рис. 11, а) и размытия (а) (рис. 11, б) исходных
изображений. Показано, что для незашумленных изображений (см. рис. 11; при стандартном отклонении шума (СО) и среднеквадратичном отклонении (СКО) фильтра размытия, равных нулю) средняя ошибка определения площади трещины составила 1,93 %, в то время как средняя ошибка определения координат ее вершины равняется 56 пикселей. При этом стандартное отклонение определения площади трещины и координат ее вершины составляет 0,21 % и 15 пикселей соответственно (согласно приведенным на графиках величинам разброса анализируемых параметров).
Зашумление изображений. Результаты анализа всех 8 пар изображений, представленных на графиках на рис. 11, а, показывают, что использование предложенного алгоритма в условиях значительного зашумления 5 = 45 обусловливает максимальную величину ошибки при детектировании (области) трещины, ёЕБ = 5,2 % (кри-
а
r
вая 1). При определении координат ее вершины (кривая 2) возрастание ошибки происходит практически линейно, и максимум данного параметра В = 90,2 (с учетом разброса) также наблюдается при максимальном СО шума 5 = 45. По мнению авторов, в рамках задачи слежения за трещиной при позиционировании оптической системы относительно поверхности образца с учетом поставленных в начале раздела технических требований подобные величины являются вполне приемлемыми.
5,5 5 4,5 4
£ 3,5
Т 1 А
ТтТт1
2,5 " 1
г
1,5
200
180
160
140
«
120 е л
е
100 ск и
п
80 С)"
60
40
20
0
20
40
60
Стандартное отклонение шума 1 2
аппрок. (1)
40 35 30 25
£ 20
К
§ 15 10 5 0 -5
аппрокс. (2)
т Т т
500 450 400 350 300 250 200 150 100 50 0
Среднеквадратическое отклонение фильтра размытия (а)
1
аппрокс. (1)
б
2
аппрокс. (2)
Рис. 11. Графики зависимости ошибок определения площади трещины и координат ее вершины от уровня шума (а) и размытия (б); 1 - ошибка определения площади трещины (dES); 2 - ошибка определения местоположения вершины трещины (D)
Fig. 11. Errors graphs on the crack area and crack tip coordinates determination vs the noise level (a) and blurring (b); The determination error of the crack tip location (D) - 2; the crack area (dES) - 1
рис. 11, б), которая возрастает в значительно большей степени и достигает ёЕБ = 36 % при размытии фильтром с о = 4,5 (кривая 1). При определении координат вершины трещины (кривая 2) возрастание ошибки происходит подобным образом; максимум данного параметра В = 280,2 также наблюдается при максимальном СКО фильтра размытия о = 4,5 .
В постановке задачи исследований было указано, что точность определения участка, содержащего вершину трещины, должна составлять порядка 20 % от горизонтального размера растра, т.е. при изображениях размером не менее 2000*1000 пикселей максимально допустимая величина ошибки определения координат составляет 400 пикселей. В случае всего исследованного диапазона зашумления (см. рис. 11, а) этот порог не оказывается превышенным, в то время как при размытии уже при о > 3 ошибка при работе алгоритма превышает допустимое значение (см. рис. 11, б). Это обусловлено тем, что большая часть трещины имеет расстояние между берегами большее, чем апертура фильтра размытия. Данное обстоятельство, в свою очередь, накладывает ограничения на минимальное значение раскрытия трещины, которое может быть определено с помощью предложенного алгоритма.
С другой стороны, введено ограничение на максимальную ошибку определения координат вершины трещины, не превышающую 15 % от их истинного положения. Результаты детектирования трещины в условиях зашумления полностью удовлетворяют данному условию, так как ошибка не превышает ёЕБ < 5,5 %. В то же время, при имитации расфокусировки, предельное значение этой ошибки достигается уже при о > 1. Это означает, что алгоритм может быть эффективно использован только для случаев небольших расфокусировок оптической системы.
Результат выделения трещины на экспериментальных изображениях представлен на рис. 12. Видно, что алгоритмически выделенный контур четко повторяет границы трещины, что свидетельствует в пользу возможности применения предложенного алгоритма для решения практических задач выделения трещин на реальных оптических изображениях.
Рис. 12. Выделенная трещина на экспериментальном изображении
Fig. 12. A crack detected on the experimental image
4. Использование DIC для изучения процессов деформирования и разрушения образцов с трещиной при статических и циклических нагрузках
Размытие изображений. Размытие изображений более критично влияет на величину ошибки (см.
Несмотря на более чем столетнюю историю исследований, многие аспекты усталостного разрушения ма-
2
0
а
териалов до сих пор остаются не до конца изученными. В частности, пристальное внимание уделяется процессам, происходящим впереди вершины трещины на микроуровне, а именно формированию пластических монотонной и циклической зон, их эволюции в процессе распространения трещины, особенно в случаях неоднородного нагружения.
Одним из наиболее распространенных методов изучения усталостного разрушения является фрактографиче-ский анализ изображений поверхностей изломов при помощи оптического и растрового электронного микроскопа. Однако этот метод имеет ряд существенных недостатков, прежде всего получение результатов только на уже разрушенном образце и, как следствие, невозможность мониторинга эволюции трещины и количественной оценки параметров ее роста во время циклического нагружения, а также необходимость дорогостоящего оборудования, в частности растрового электронного микроскопа, и высококвалифицированного оператора. При этом данный метод не позволяет напрямую количественно измерять многие параметры, используемые в механике разрушения.
Другим известным способом исследования процессов роста трещины является применение датчика раскрытия трещины (в англоязычной литературе - crack opening displacement, COD), который также не лишен ряда недостатков. Например, разрешающая способность COD-датчика не позволяет получить оценку изменений, происходящих в материале при неоднородном нагруже-нии, и для решения подобных задач данный подход не эффективен.
Перспективным способом преодоления описанных ограничений является применение метода корреляции цифровых изображений, который способен измерять механический отклик материала впереди вершины трещины при приложении как единичных циклов перегрузки, так и при случайном спектре нагружения. Однако для успешного применения данного метода необходима дальнейшая разработка аппаратно-программных средств с целью а) повышения точности и стабильности получаемых результатов (количественных оценок); б) разработки методик подготовки поверхности и регистрации изображений поверхности для дальнейшего расчета полей деформаций; в) разработки методов расчета параметров, используемых в механике разрушения, на основе построения полей деформаций для проведения последующей оценки состояния материала и прогнозирования остаточного ресурса.
Успешные исследования в области определения величины закрытия трещины при помощи корреляции цифровых изображений в процессе циклических испытаний (при наличии единичного цикла перегрузки, для образцов разной толщины, с целью исследования механизма распространения трещины) продемонстрированы в работах [54, 55, 82, 83].
Другой широко используемый подход основан на предположении о том, что на процесс роста трещины оказывает влияние зона пластичности и ее эволюция.
Концепция зоны пластичности служит для объяснения эффекта замедления трещины вследствие приложения единичного цикла перегрузки во время циклического испытания с постоянной амплитудой.
Аналитическая формулировка процесса роста трещины на микроуровне была разработана авторами [84]. Если при помощи традиционно используемых моделей закон роста трещины формулируется в зависимости от количества циклов нагружения da/dN, то описываемый подход основан на определении кинетики роста трещины в каждый конкретный момент времени в процессе нагружения, т.е. с использованием параметра da/dt. Основным предположением данной модели является то, что рост трещины зависит от взаимодействия монотонной и циклической пластических зон, размеры которых зависят от уровня закрытия трещины. Получение данных о распределении деформаций и напряжений вблизи вершины трещины при помощи оптических методов проводилось рядом авторов [56, 85-87].
Дальнейшим продолжением развития метода корреляции цифровых изображений является использование его не только для измерения величины закрытия или построения полей деформаций, но и для расчета параметров механики разрушения, таких как ./-интеграл и коэффициент интенсивности напряжений K.
/-интеграл, предложенный независимо Г.П. Черепановым и Дж. Райсом [88, 89], определяется как контурный интеграл, не зависящий от пути интегрирования, который характеризует интенсивность выделения энергии при распространении трещины. Обычно /-интеграл определяется по любому замкнутому контуру, включающему в себя вершину трещины. Преимуществом применения /-интеграла в расчетах является возможность его использования в случаях больших пластических деформаций материала (когда размер пластической зоны превышает 20 % от длины трещины).
В работах [58-60, 90] авторы использовали метод DIC для расчета полей деформаций на поверхности материала, а затем импортировали это поле в конечно -элементную модель для последующего расчета /интеграла.
Определение КИН с использованием DIC может быть выполнено с использованием трех подходов:
1) на основании полей смещений, полученных методом DIC, совмещенных с результатами решения степенных рядов Вестергарда (Westergaard) и Вильямса (Williams) [57, 61, 91-93];
2) на основании показаний программных датчиков деформаций в методе DIC, размещаемых на берегах раскрывающейся трещины [54, 62, 94];
3) на основе расчетов /-интеграла [46, 95, 96].
4.1. Методика расчета параметров механики усталостного разрушения при помощи метода корреляции цифровых изображений
Механический отклик материала. Одним из параметров, позволяющих качественно охарактеризовать поведение материала вблизи вершины трещины, является петля механического гистерезиса, где по оси Oy откладывается величина нагрузки, определенная по датчику нагрузки ИМ, а по оси Ox - величина раскрытия, определяемая по программному экстензометру, размещенному на поле смешений (ШС). Диаграмма нагружения и разгрузки принимает вид петли вследствие нелинейности отклика материала (рис. 13). При этом увеличение доли пластической деформации ведет к увеличению «ширины» петли, тогда как ее наклонный характер вызван влиянием открытия/закрытия трещины. Описываемая петля механического гистерезиса состоит из трех участков:
• начальный линейный участок, когда трещина закрыта;
• участок нелинейности, вызванный постепенным открытием трещины на полуцикле нагрузки и ее постепенным закрытием на полуцикле разгрузки;
• конечный линейный участок, когда трещина полностью открыта.
Таким образом, анализируя внешний вид петли, можно давать качественную характеристику происходящим в материале процессам. Численной оценкой может выступать максимальное значение деформации (Стек), достигаемое в каждой петле, что косвенно характеризует раскрытие трещины через податливость образца.
Оценка уровня открытия/закрытия трещины. Один из важных численных параметров, непосредственно влияющих на скорость распространения трещины, особенно на начальном этапе ее роста и в режиме Пэриса, является уровень открытия/закрытия трещины, который можно определить как эффективную асимметрию цикла:
P
D О mean
Reff - 2 -р—
-1.
P
R* - ■ ^
eff
P
(10)
где Pop - нагрузка, при которой трещина открывается; Pmax - максимальная нагрузка в цикле.
В идеальном случае, когда закрытие трещины отсутствует, петля механического гистерезиса должна иметь симметричную форму. То есть если провести вертикальною или горизонтальную прямую, которая будет делить петлю на две равные по площади части, то эта прямая и будет осью симметрии петли. При закрытии площадь нижней части петли станет меньше площади верхней. На этом предположении был разработан метод оценки уровня закрытия усталостной трещины на основе расчета площади петли механического гистерезиса.
Процедура расчета (см. рис. 13) уровня закрытия состоит:
1) из определения площади, очерченной петлей;
2) из определения такого уровня нагрузки Pmean, линия которой «разрезала» бы петлю на две равные по площади части - S1 и S2;
3) из расчета эффективной асимметрии цикла (Reff), которая находится из выражения
а,
Р
Р
АР = Р™*~ „ /1 / 1 / 1
/ 5, /
1
1 АР 1 ' efT
АР / S, / 1 1 1 1 1 1 1 1 1 1
Деформация (с)
Рис. 13. Схема расчета уровня открытия/закрытия трещины
Fig. 13. The scheme illustrating calculation principle of crack opening/closure level
Оценка параметров роста усталостной трещины после единичной перегрузки. Для проведения циклических испытаний использовался образец алюминиевого сплава Д16Т. Форма и основные размеры образца представлены на рис. 14.
Рис. 14. Форма образца. Толщина 3 мм
Fig. 14. Specimen's shape and dimension.
Thickness 3 mm
Данный образец имеет боковой надпил, в вершине которого происходит рост усталостной трещины. Перед началом испытаний трещина была предварительно выращена на длину 1,7 мм. Область фотографирования составляла 2,5 мм, что составляет в физическом масштабе 1,04 пикселей/мкм.
Схема нагружения представлена на рис. 15 и состоит из следующих этапов:
1) производится циклическое нагружение в базовом диапазоне нагрузок (0,25-2,5 кН) до тех пор, пока трещина не достигнет величины 2,5 мм (т.е. прирост должен составить 0,8 мм);
2) после проращивания трещины прикладывается цикл нагружения, в течение которого камера производит фотографирование поверхности образца с частотой 3 кадра/с, этот цикл - опорный;
3) затем прикладывается перегрузка с амплитудой на 50 % больше базовой, до 3,75 кН, во время которой также регистрируются изображения;
4) после приложения перегрузки происходит дальнейшее циклическое нагружение на базовом уровне,
однако в отдельных циклах запускается регистрация фотографий, чтобы проследить изменение характера роста усталостной трещины. Номера циклов, в которых производится фотографирование, после перегрузки: 1, 2, 3, 10, 25, 100, 500, 4000, 10000, 15000, 24000, 28000 и 37000 циклов.
Рис. 15. Схема нагружения для испытания с единичной перегрузкой. Параметры базового цикла Pmax=2,5 кН, Рпй1=0Д8 кН, асимметрия Л=0,1. Цикл перегрузки Povl=3,75 кН. Длина трещины a=14,5 мм. Коэффициенты интенсивности Kmax= 12,86 МПа^м, Kmin= 1,29 МПа^м, Kovl=19,3 МПа^м
Fig. 15. The loading sequence scheme for testing with a single overload. The parameters of the basic cycle Pmax = 2.5 kN, Pmin = 0.28 kN, asymmetry R = 0.1. The overload cycle Povl = 3.75 kN. The crack length is a = 14.5 mm. The stress intensity factors Kmax = 12.86 MPaVm, Kmn = 1.29 MPaVm, Kovl = 19.3 MPaVm
Циклы, в которых регистрировались фотографии образца, имели частоту нагружения 0,05 Гц с одновременным фотографированием с частотой 3 кадра/с. Таким образом, записывали 60 кадров в течение цикла «нагрузки-разгрузки». Между циклами измерения продолжалось «базовое» циклирование, но с большей частотой - 5 Гц.
Полученные серии испытания изображений далее обрабатывают при помощи метода корреляции цифровых изображений с целью вычисления полей деформаций и параметров механики разрушения.
Результаты исследований и их обсуждение. Параметром, характеризующим изменение реакции материала на прикладываемую нагрузку, является максимальное значение деформации в вершине трещины -£max, представленное на графике рис. 16. Видно, что Emax после перегрузки достаточно быстро снижается, пока трещина не достигнет длины в 20% rOvL, где стесненность материала вблизи вершины трещины максимальная - материал деформируется на минимальную величину при постоянных прикладываемых нагрузках. Затем следует достаточно быстрый рост значений emax до уровня выше начального (до перегрузки).
Результаты оценки уровня открытия/закрытия трещины представлены на рис. 17. Горизонтальной линией
отмечено значение асимметрии при базовом уровне нагрузок Я=0,1, т.е. значения меньше 0,1 означают, что трещина полностью открыта в течение всего цикла.
До перегрузки закрытие трещины минимально и составляет всего 0,115 от Ртах. Перегрузка вызывает притупление вершины трещины, вследствие чего трещина остается полностью открытой во время всего цикла нагрузки-разгрузки. Однако, как уже отмечалось выше, открытие и закрытие трещины процесс протяженный во времени, поэтому точки, лежащие вблизи Я = 0,1, могут оказывать влияние на скорость роста трещины, что и наблюдается в данном случае. Максимальный уровень закрытия достигается при -30% гоуь. При дальнейшем росте трещины закрытие снижается и остается на уровне немного выше базового (до перегрузки).
Рис. 16. Зависимость максимальной деформации Smax от расстояния, на которое трещина проросла в пластическую зону, вызванную перегрузкой
Fig. 16. Dependence of the maximum strain Smax vs the crack extension into the overloading induced plastic zone
Длина трещины/г0,,ь,
Рис. 17. Зависимость Reff от расстояния, на которое трещина проросла в пластическую зону, вызванную перегрузкой
Fig. 17. Dependence of the Reff vs the crack extension into the overloading induced plastic zone
4.2. Развитие алгоритма измерения J-интеграла
./-интеграл [88, 89] представляет собой независимый от пути (path-independent) линейный интеграл, который характеризует напряжения и деформации в вершине трещины, возникающие под действием упругих сил. Расчет /-интеграла происходит по заданному контуру интегрирования, который должен охватывать вершину трещины. Одним из преимуществ применения /-интеграла является возможность использования измеренных полей перемещений в областях, расположенных достаточно далеко от вершины трещины, в которых преобладают линейные упругие деформации и которые с достаточно высокой точностью могут быть определены при помощи DIC. Поэтому в расчетах могут быть использованы принципы Linear Elastic Fracture Mechanics (LEFM) и из уравнений теории упругости можно вычислить соответствующие напряжения и деформации, необходимые при оценке /-интеграла. В представленном исследовании описан алгоритм измерения /-интеграла с использованием метода корреляции цифровых изображений, и проведена оценка помехоустойчивости расчета в рамках применения метода корреляции цифровых изображений.
Процесс измерения /-интеграла и выделение пластической зоны на базе вычисления оптического потока при механических испытаниях материалов можно разделить на четыре этапа:
1) расчет полей перемещений, основанный на сопоставлении участков изображений, оригинальный алгоритм построения которых развит и описан авторами в работах [81, 97];
2) постобработка полей перемещений (сглаживание оптического потока);
3) расчет компонент деформации и напряжений;
4) вычисление /-интеграла и зоны пластичности в пределах построенного векторного поля.
Сглаживание оптического потока. Наличие шумов на изображении, искажение поверхности, неоднородность освещения и пр. обусловливают появление ошибок при оценке перемещений, что, в свою очередь, приводит к значительной зашумленности полей деформации вследствие применения процедуры численного дифференцирования. Частично данная проблема может быть решена посткорректировкой поля векторов перемещений (пространственная фильтрация и сглаживание) [98, 99].
Расчет компонент деформации, напряжений и J-интеграла. Дифференцируя векторное поле, можно вычислить компоненты деформации [100]. Выражения для продольной &xx, поперечной гуу и сдвиговой eXy компонент тензора дисторсии имеют вид
дПх дх
dUy
ду
1
,&ху 2
ди dUy} —х +—у
ду дх
(11)
где Ц., и - компоненты вектора перемещений участков изображения по осям х, у.
Вычисление производных проводили с привлечением процедуры численного дифференцирования. Безразностная формула для его реализации имеет вид [101]
у'( х )-К-1 )n
C'nyi d
t (t-l)...(t - n)
t - i
, (12)
п ! dt
где к - номер отсчета функции; п - величина, определяющая уровень приближения; ст - величина, рассчитываемая по следующей формуле:
с -
t (t-l)...(t - m +1)
m!
(13)
где т - натуральное число и / - произвольное действительное число.
При известных смещениях, деформациях и напряжениях ./-интеграл может быть рассчитан по известной формуле [88, 89]
(14)
J = [ | wdy — T — ds |,
Г1 dx )
где w - плотность энергии деформации,
£
w = jod £, (15)
о
Т - вектор нагрузок, T = ои; Г - контур произвольной формы, охватывающий вершину трещины; n - вектор нормали к Г; с, е, u - напряжения, деформации и поля перемещений соответственно.
В рамках линейно-упругой механики и в условиях плосконапряженного состояния тензор напряжений можно определить на основе закона Гука из тензора деформации:
ст- De-
ст
хх E
ст —
уу 1 -v2
X
[ ху _
1 v 0 v 1 0 0 0 1 -v
(16)
где О - матрица упругих модулей при плоском напряженном состоянии; Е - модуль упругости; V - коэффициент Пуассона. Таким образом, когда известны поля тензора напряжений и деформаций, /-интеграл можно оценить из уравнения (14) следующим образом [102]:
• = |1 + е * + У ху) ^ -
-[
ст-Пх + ххл ' Xxynx +стyyny _
ди дх äv .дх.
ds. (17)
Здесь Г - контур, охватывающий вершину трещины, обход по которому осуществляется против часов стрелки; п, П - компоненты вектора нормали к Г; и, V -компоненты вектора перемещения.
e -
Для верификации исследуемого метода измерения ./-интеграла использовался программный комплекс ANSYS для численного моделирования процесса статического растяжения образца с заданными геометрическими размерами (рис. 18). С помощью ANSYS были получены модельные поля перемещений (рис. 19), деформаций и напряжений, соответствующие величинам прикладываемой нагрузки, а также рассчитаны значения /-интеграла (таблица, рис. 20).
Рис. 18. Расчетная модель образца (расчет проведен с помощью пакета ANSYS)
Fig. 18. Calculated model of a specimen (the calculation was carried out using the ANSYS package)
0,40156 Mai
0,34864 0,29571 0,24279 0,18987 0,13694 0,084017 0,031093 -0,021831 0,074755 Min
0,19683 Mix
0,11667
0,036514
-1,143646
-1,12381
-1,21397
-1,28413
-1,36428
-1,44444
-0,5246 Min
б
Рис. 19. Результаты моделирования в ANSYS: directional deformation (X axis (а), Y axis (б))
Fig. 19. ANSYS simulation results : directional deformation (X axis (a), Y axis (b))
Оценка ошибки вычисления J-интеграла. Для получения количественной оценки «точности» и помехоустойчивости расчета ./-интеграла предложено использовать следующую расчетную величину:
Err =
Jмод Jpac4
J,,,
100%,
(18)
где /мод - величина /-интеграла, полученная с помощью пакета ANSYS; J - рассчитанное значение /-интеграла.
Из модельных перемещений по формуле (2) были рассчитаны значения /-интеграла (таблица, рис. 20). Видно, что с увеличением нагрузки возрастает отклонение измеряемой величины от модельных данных. Ошибка увеличивается от 0,06 до 3,95 %, при этом средняя ошибка для используемых величин нагрузки составляет 1,75 %, что говорит о применимости метода корреляции цифровых изображений для измерения /-интеграла.
Значения /-интеграла и ошибки вычисления /-интеграла
The values and calculation errors of the /-integral
Вели- Нагрузка, кН
чина 2 3 4 5 6 7 8 9 10 11 12 13
Jмод , Н/м 46 3 104 3 186 8 294 1 427 6 588 5 778 9 100 11 125 54 154 61 187 65 236 98
J , расч ' Н/м 46 3 104 1 185 9 291 9 423 0 580 1 764 6 978 3 122 24 149 94 181 33 227 61
Err, % 0,0 6 0,1 8 0,4 5 0,7 4 1,0 7 1,4 3 1,8 4 2,28 2,62 3,02 3,36 3,95
25000т-
Рис. 20. График изменения величины /-интеграла при возрастании прикладываемой нагрузки: 1 - расчет при помощи пакета ANSYS; 2 - по предлагаемой методике
а
Fig. 20. The J-integral variation with increasing applied load: 1 - calculation with the ANSYS; 2 - calculations based on the proposed method
Заключение
В рамках развития метода корреляции цифровых изображений разработана компьютерная система сте-реозрения и проведено тестирование на сериях стереопар, отражающих изменение положения объекта в пространстве, плоскостные и внеплоскостные деформации. Показано, что ошибка определения пространственных координат не превышает 0,75 ед. при угле между осями камер 20°. Предложенные алгоритмы позволяют определять плоскую и внеплоскостную деформацию с ошибкой, не превышающей 4-10-3. Ошибка вычисления компонент тензора деформации в случае системы с одной камерой на два порядка больше, чем при использовании системы стереозрения в случае внеплоскостных деформаций.
Предложен модифицированный инкрементальный алгоритм определения перемещений на серии стереопар, который позволяет оценивать большие перемещения при последовательной обработке изображений. Представлен алгоритм вычисления деформации на поверхности пространственного объекта с использованием восстановленных профилей поверхности объекта.
Предложен и протестирован алгоритм автоматического детектирования трещины на оптических изображениях и определения координат ее вершины. Результаты тестирования с использованием модельных изображений (включая размытие и зашумление) показали, что алгоритм позволяет определить координаты вершины трещины со средней ошибкой порядка 56 пикселов. Средняя ошибка определения площади трещины (включая выделение ее берегов) не превышает 1,93 %. Результаты тестирования алгоритма на экспериментальных изображениях показали возможности его эффективного применения для решения задач выделения трещин на реальных (оптических) изображениях.
Проведена количественная характеризация процесса роста усталостной трещины в металлических сплавах с помощью методики, основанной на методе корреляции цифровых изображений (в терминах механики раз-
Библиографический список
1. Методы компьютерной обработки изображений. - 2-е изд. / М.В. Гашников [и др.]; ред. В.А. Сойфер. - М.: Физмат-лит, 2003. - 784 с.
2. Barron J.L., Fleet D.J., Beauchemin S.S. Performance of optical flow techniques // International Journal of Computer Vision. - 1994. - Vol. 12. - No. 1. - P. 43-77.
3. Beauchemin S.S., Barron J.L. The computation of optical flow // ACM Computing Surveys. - 1995. - Vol. 27. - No. 3. -P. 433-466.
4. Horn B.K.P., Schunck B.G. Determining optical flow // Artificial Intelligence. - 1981. - Vol. 17. - No. 1-3. - P. 185203.
рушения), включающей расчет: а) скорости роста усталостной трещины da/dN; б) максимального значения деформации (emax) и в) эффективной асимметрии цикла (Ref). При испытании алюминиевого сплава Д16АТ показано, что в режиме Пэриса при единичной перегрузке минимальная скорость роста трещины достигается при длине трещины 30 % от размера пластической зоны, сформированной циклом перегрузки, и вызвана закрытием трещины, увеличивающим эффективную асимметрию цикла до Reff = 0,32.
Развит алгоритм измерения /-интеграла с использованием метода корреляции цифровых изображений, заключающийся в измерении векторов перемещений, вычислении деформаций и напряжений с последующим расчетом значения /-интеграла. Его тестирование показало, что отклонение рассчитываемых значений /-интеграла от модельных составляет в среднем 1,75 % при минимальной величине отклонения 0,06 % и максимальной 3,95 %.
Благодарности
Работа выполнена в рамках программы фундаментальных исследований государственных академий наук 2013-2020 гг. (направление 23.2), при финансовой поддержке РФФИ и администрации Томской области в рамках научного проекта № 18-47-700006, грантов Президента РФ № НШ-5875.2018.8, МК-6762.2018.9, СП-198.2018.4. Авторы благодарят А.А. Козулина за проведение расчетов /-интеграла с помощью пакета ANSYS.
Acknowledgements
The work was performed as part of the Fundamental Research Program of State Academies of Sciences during 20132020. (specialization 23.2), with the financial support of the Russian Foundation for Basic Research and the Administration of Tomsk Region within the Research Project Nr. 18-47700006, Grants of the President of the Russian Federation Nr. NSh-5875.2018.8, MK-6762.2018.9, SP-198.2018.4. The authors acknowledge A.A. Kozulin for performing J-integral calculations using the ANSYS package.
5. Lucas B.D., Kanade T. An iterative image registration technique with an application to stereo vision // Proceedings of Imaging Understanding Workshop. - 1981. - P. 121-130.
6. Lucas B.D. Generalized image matching by the method of differences. - Carnegie Mellon University, 1984. - 167 p.
7. Fennema C.L., Thompson W.B. Velocity determination in scenes containing several moving objects // Computer Graphics and Image Processing. - 1979. - Vol. 9. - No. 4. - P. 301-315.
8. Fleet D.J., Jepson A.D., Jenkin M.R.M. Phase-based disparity measurement // CVGIP: Image Understanding. - 1991. -Vol. 53. - No. 2. - P. 198-210.
9. Adelson E.H., Bergen J.R. Spatiotemporal energy models for the perception of motion // Journal of the Optical Society of America A. - 1985. - Vol. 2. - No. 2. - P. 284.
10. Heeger D.J. Optical flow using spatiotemporal filters // International Journal of Computer Vision. - 1988. - Vol. 1. -No. 4. - P. 279-302.
11. Watson A.B., Ahumada A.J. A look at motion in the frequency domain. - 1983.
12. Fleet D.J. The early processing of spatio-temporal visual information: thesis (M.Sc.). - University of Toronto, 1984. -136 p.
13. Schreier H., Orteu J.-J., Sutton M.A. Image correlation for shape, motion and deformation measurements. - Boston, MA: Springer US, 2009. - 321 p.
14. Обзор блочных методов оценки движения в цифровых видеосигналах / С.В. Гришин [и др.] // Программные системы и инструменты: тем. сб. - М.: Изд-во фак. ВМиК МГУ, 2008. - С. 50-62.
15. Zhu S., Ma K.K. A new diamond search algorithm for fast block-matching motion estimation // IEEE Transactions on Image Processing. - 2000. - Vol. 9. - №. 2. - P. 287-290.
16. Li W., Salari E. Successive elimination algorithm for motion estimation // IEEE Transactions on Image Processing. -1995. - Vol. 4. - №. 1. - P. 105-107.
17. Schreier H.W. Systematic errors in digital image correlation caused by intensity interpolation // Optical Engineering. -2000. - Vol. 39. - №. 11. - P. 2915-2921.
18. Wu Q.X. A correlation-relaxation-labeling framework for computing optical flow-template matching from a new perspective // IEEE Transactions on Pattern Analysis and Machine Intelligence. - 1995. - Vol. 17. - № 9. - P. 843-853.
19. Black M.J., Anandan P. The robust estimation of multiple motions: parametric and piecewise-smooth flow fields // Computer Vision and Image Understanding. - 1996. - Vol. 63. -№ 1. - P. 75-104.
20. Воскобойников Ю.Е., Гочаков А.В. Сравнительный анализ двух классов алгоритмов фильтрации изображений // Автометрия. - 2011. - Т. 47, № 1. - С. 17-29.
21. Воскобойников Ю.Е., Гочаков А.В. Построение алгоритмов вейвлет-фильтрации с двухпараметрическими пороговыми функциями // Автометрия. - 2012. - Т. 48, № 1. -С. 12-22.
22. Otte M., Nagel H.-H. Estimation of optical flow based on higher-order spatiotemporal derivatives in interlaced and noninterlaced image sequences // Artificial Intelligence. - 1995. -Vol. 78. - № 1-2. - P. 5-43.
23. Nagel H.-H. Displacement vectors derived from second-order intensity variations in image sequences // Computer Vision, Graphics, and Image Processing. - 1983. - Vol. 21. - №. 1. -P. 85-117.
24. Nagel H.-H., Enkelmann W. An Investigation of smoothness constraints for the estimation of displacement vector fields from image sequences // IEEE Transactions on Pattern Analysis and Machine Intelligence. - 1986. - Vol. PAMI-8. - №. 5. -P. 565-593.
25. Wang C.-M., Fan K.-C., Wang C.-T. Estimating Optical Flow by Integrating Multi-Frame Information // Journal of Information Science and Engineering. - 2008. - Vol. 24. - №. 6. -P. 1719-1731.
26. Irani M. Multi-frame optical flow estimation using subspace constraints // Proceedings of the Seventh IEEE International Conference on Computer Vision. - IEEE, 1999. - Vol. 1. -P. 626-633.
27. Gao W., Kemao Q. Parallel computing in experimental mechanics and optical measurement: A review // Optics and Lasers in Engineering. - 2012. - Vol. 50. - P. 608-617.
28. High accuracy digital image correlation powered by GPU-based parallel computing / L. Zhang [et al.] // Optics and Lasers in Engineering. - 2015. - Vol. 69. - P. 7-12.
29. Shao X., Dai X., He X. Noise robustness and parallel computation of the inverse compositional Gauss-Newton algorithm in digital image correlation // Optics and Lasers in Engineering. - 2015. - Vol. 71. - P. 9-19.
30. A multi-resolution approach for massively-parallel hardware-friendly optical flow estimation / F. Barranco [et al.] // Journal of Visual Communication and Image Representation. -2012. - Vol. 23. - No. 8. - P. 1272-1283.
31. Marzat J., Dumortier Y., Ducrot A. Real-time dense and accurate parallel optical flow using CUDA // 7th International Conference WSCG. - 2009. - P. 105-111.
32. SimpleFlow: A Non-iterative, Sublinear Optical Flow Algorithm / M. Tao [et al.] // Computer Graphics Forum. - 2012. -Vol. 31, 2pt1. - P. 345-353.
33. Rav-Acha A., Peleg S. Lucas-Kanade without iterative warping // 2006 International Conference on Image Processing. -IEEE, 2006. - P. 1097-1100.
34. Braspenning R.A., Haan G. de True-motion estimation using feature correspondences. Eds. S. Panchanathan, B. Vasu-dev. - 2004. - P. 396.
35. True-motion estimation with 3-D recursive search block matching / G. de Haan [et al.] // IEEE Transactions on Circuits and Systems for Video Technology. - 1993. - Vol. 3. - No. 5. -P. 368-379, 388.
36. Optimization of hierarchical 3DRS motion estimators for picture rate conversion / A. Heinrich [et al.] // IEEE Journal of Selected Topics in Signal Processing. - 2011. - Vol. 5. - No. 2. -P. 262-274.
37. Ahn T.G., Moon Y.H., Kim J.H. Fast full-search motion estimation based on multilevel successive elimination algorithm // IEEE Transactions on Circuits and Systems for Video Technology. - 2004. - Vol. 14. - No. 11. - P. 1265-1269.
38. Lopes F., Ghanbari M. Hierarchical motion estimation with spatial transforms // Proceedings 2000 International Conference on Image Processing (Cat. No.00CH37101). - IEEE, 2000. -Vol. 2. - P. 558-561.
39. Application of an optimized digital correlation method to planar deformation analysis / M. Sutton [et al.] // Image and Vision Computing. - 1986. - Vol. 4. - No. 3. - P. 143-150.
40. Applications of digital-image correlation techniques to experimental mechanics / T.C. Chu [et al.] // Experimental Mechanics. - 1985. - Vol. 25. - No. 3. - P. 232-244.
41. Determination of displacements using an improved digital correlation method / M. Sutton [et al.] // Image and Vision Computing. - 1983. - Vol. 1. - No. 3. - P. 133-139.
42. Digital image correlation using Newton-Raphson method of partial differential correction / H.A. Bruck [et al.] // Experimental Mechanics. - 1989. - Vol. 29. - No. 3. - P. 261-267.
43. Advances in two-dimensional and three-dimensional computer vision / M.A. Sutton [et al.] // Photomechanics. -Springer Berlin Heidelberg, 2000. - P. 323-372.
44. Hild F., Roux S. Digital image correlation: from displacement measurement to identification of elastic properties -a review // Strain. - 2006. - Vol. 42. - No. 2. - P. 69-80.
45. Besnard G., Hild F., Roux S. "Finite-element" displacement fields analysis from digital images: application to Portevin-Le Chatelier bands // Experimental Mechanics. - 2006. -Vol. 46. - No. 6. - P. 789-803.
46. Roux S., Hild F. Stress intensity factor measurements from digital image correlation: Post-processing and integrated approaches // International Journal of Fracture. - 2006. -Vol. 140. - №. 1-4. - P. 141-157.
47. Multiscale displacement field measurements of compressed mineral-wool samples by digital image correlation / F. Hild [et al.] // Applied Optics. - 2002. - Vol. 41. - №. 32. -P. 6815.
48. Two-dimensional digital image correlation for in-plane displacement and strain measurement: a review / B. Pan [et al.] // Measurement Science and Technology. - 2009. - Vol. 20. -№. 6. - P. 062001.
49. Performance of sub-pixel registration algorithms in digital image correlation / P. Bing [et al.] // Measurement Science and Technology. - 2006. - Vol. 17. - № 6. - P. 1615-1621.
50. Study on subset size selection in digital image correlation for speckle patterns / B. Pan [et al.] // Optics Express. -2008. - Vol. 16. - № 10. - P. 7037.
51. Digital image correlation using iterative least squares and pointwise least squares for displacement field and strain field measurements / B. Pan [et al.] // Optics and Lasers in Engineering. - 2009. - Vol. 47. - № 7-8. - P. 865-874.
52. Pan B. Recent progress in digital image correlation // Experimental Mechanics. - 2011. - Vol. 51. - № 7. - P. 1223-1235.
53. Механика материалов. Методы и средства экспериментальных исследований / В.Э. Вильдеман [и др.]. - Пермь: Изд-во Перм. нац. исслед. политехн. ун-та, 2011. - 164 с.
54. Investigation of fatigue crack closure using multiscale image correlation experiments / J. Carroll [et al.] // Engineering Fracture Mechanics. - 2009. - Vol. 76. - № 15. - P. 2384-2398.
55. Mathieu F., Hild F., Roux S. Identification of a crack propagation law by digital image correlation // International Journal ofFatigue. - 2012. - Vol. 36. - № 1. - P. 146-154.
56. High resolution digital image correlation measurements of strain accumulation in fatigue crack growth / J.D. Carroll [et al.] // International Journal of Fatigue. - 2013.
57. Zhang R., He L. Measurement of mixed-mode stress intensity factors using digital image correlation method // Optics and Lasers in Engineering. - 2012.
58. An approach to calculate the J-integral by digital image correlation displacement field measurement / T.H. Becker [et al.] // Fatigue and Fracture of Engineering Materials and Structures. -2012.
59. Analysis of standard fracture toughness test based on digital image correlation data / I. Jandejsek [et al.] // Engineering Fracture Mechanics. - 2017.
60. A J-integral approach using digital image correlation for evaluating stress intensity factors in fatigue cracks with closure effects / G.L.G. Gonzales [et al.] // Theoretical and Applied Fracture Mechanics. - 2017.
61. Yoneyama S., Morimoto Y., Takashi M. Automatic evaluation of mixed-mode stress intensity factors utilizing digital image correlation // Strain. - 2006.
62. Investigation of thermal effects on fatigue crack closure using multiscale digital image correlation experiments / M.C. Cas-person [et al.] // International Journal of Fatigue. - 2014.
63. Панин С.В., Сырямкин В.И., Любутин П.С. Оценка деформации твердых тел по изображениям поверхности // Автометрия. - 2005. - Т. 41, № 2. - С. 44-58.
64. Giachetti, A. Matching techniques to compute image motion // Image and Vision Computing. - 2000. - Vol. 18. -№ 3. - P. 247-260.
65. Tong W. An evaluation of digital image correlation criteria for strain mapping applications // Strain. - 2005. - Vol. 41. -№ 4. - P. 167-175.
66. Full-field strain measurement using a two-dimensional Savitzky-Golay digital differentiator in digital image correlation /
B. Pan [et al.] // Optical Engineering. - 2007. - Vol. 46, N 3. -P. 33601-33610.
67. Cofaru C., Philips W., Van Paepegem W. Improved New-ton-Raphson digital image correlation method for full-field displacement and strain calculation // Applied Optics. - 2010. -Vol. 49. - №. 33. - P. 6472.
68. Fundamental role of crystal structure curvature in plasticity and strength of solids / V.E. Panin [et al.] // Physical Me-somechanics. - 2015. - Vol. 18. - №. 2. - P. 89-99.
69. Synergy of crack closure, near-tip residual stress and crack-tip blunting in crack growth under periodic overloads -A fractographic study / R. Sunder [et al.] // International Journal of Fatigue. - 2016. - Vol. 93. - P. 18-29.
70. Song K.Y., Petrou M., Kittler J. Texture crack detection // Machine Vision and Applications. - 1995. - Vol. 8. - №. 1. -P. 63-75.
71. Kirschke K.R., Velinsky S.A. Histogram based approach for automated pavement crack sensing // Journal of Transportation Engineering. - 1992. - Vol. 118. - №. 5. - P. 700-710.
72. CrackTree: Automatic crack detection from pavement images / Q. Zou [et al.] // Pattern Recognition Letters. - 2012. -Vol. 33. - №. 3. - P. 227-238.
73. Sun L., Kamaliardakani M., Zhang Y. Weighted neighborhood pixels segmentation method for automated detection of cracks on pavement surface images // Journal of Computing in Civil Engineering. - 2016. - Vol. 30. - №. 2. - P. 04015021.
74. Yamaguchi T., Hashimoto S. Fast crack detection method for large-size concrete surface images using percolation-based image processing // Machine Vision and Applications. -2010. - Vol. 21. - №. 5. - P. 797-809.
75. Iyer S., Sinha S.K. A robust approach for automatic detection and segmentation of cracks in underground pipeline images // Image and Vision Computing. - 2005. - Vol. 23. - №. 10. -P. 921-933.
76. Lins R.G., Givigi S.N. Automatic crack detection and measurement based on image analysis // IEEE Transactions on Instrumentation and Measurement. - 2016. - Vol. 65. - №. 3. -P. 583-590.
77. A database and evaluation methodology for optical flow / S. Baker [et al.] // International Journal of Computer Vision. - 2011. - Vol. 92. - №. 1. - P. 1-31.
78. Панин С.В., Чемезов В.О., Любутин П.С. Комплексный алгоритм вычисления оптического потока с использованием взвешенной меры подобия // Вычислительные технологии. - 2017. - Т. 22, № 4. - С. 80-94.
79. Gonzalez R.C., Woods R.E. Digital image processing, 3rd Edition. - Prentice Hall, 2008. - 954 p.
80. Glantz R., Kropatsch W.G. Skeletons in the framework of graph pyramids // Visual Form 2001. Lecture Notes in Computer Science. Eds. C. Arcelli, L.P. Cordella, G.S. di Baja. - Springer Berlin Heidelberg, 2001. - P. 186-195.
81. Панин С.В., Титков В.В., Любутин П.С. Выбор параметров алгоритма трехмерного рекурсивного поиска при построении поля векторов перемещений с использованием иерархического подхода // Автометрия. - 2015. - Т. 51, № 2. -
C. 27-37.
82. Nowell D., De Matos P.F.P. Application of digital image correlation to the investigation of crack closure following overloads // Procedia Engineering. - 2010.
83. Yusof F., Lopez-Crespo P., Withers P.J. Effect of overload on crack closure in thick and thin specimens via digital image correlation // International Journal of Fatigue. - 2013. - Vol. 56. -P. 17-24.
84. Lu Z., Liu Y. Small time scale fatigue crack growth analysis // International Journal of Fatigue. - 2010.
85. Toribio J., Kharin V. Finite-deformation analysis of the crack-tip fields under cyclic loading // International Journal of Solids and Structures. - 2009.
86. Near-tip strain evolution under cyclic loading: In situ experimental observation and numerical modelling / J. Tong [et al.] // International Journal of Fatigue. - 2015.
87. Zhang W., Liu Y. Plastic zone size estimation under cyclic loadings using in situ optical microscopy fatigue testing // Fatigue and Fracture of Engineering Materials and Structures. - 2011.
88. Черепанов Г.П. О распространении трещин в сплошной среде // Прикладная математика и механика. - 1967. -Т. 31, № 3. - С. 476-488.
89. Rice J.R. A path independent integral and the approximate analysis of strain concentration by notches and cracks // Journal of Applied Mechanics. - 1968.
90. J-Integral calculation by finite element processing of measured full-field surface displacements / S.M. Barhli [et al.] // Experimental Mechanics. - 2017. - Vol. 57. - №. 6. - P. 9971009.
91. McNeill S.R., Peters W.H., Sutton M.A. Estimation of stress intensity factor by digital image correlation // Engineering Fracture Mechanics. - 1987. - Vol. 28. - № 1. - P. 101-112.
92. Lorenzino P., Beretta G., Navarro A. Application of digital image correlation (DIC) in resonance machines for measuring fatigue crack growth // Frattura ed Integrita Strutturale. - 2014.
93. Hos Y., Vormwald M., Freire J.L.F. Using digital image correlation to determine mixed mode stress intensity factors in fatigue cracks // Proceedings of COTEQ 2015. - 2015.
References
1. Metody komp'iuternoi obrabotki izobrazhenii 2-e izd. / M.V. Gashnikov [i dr.]; red. V.A. Soifer. - M.: FIZMATLIT 2003. -784 s. [Methods of computer image processing] 2nd ed. / M.V. Gashnikov et al.; ed. V.A. Soifer - M .: FIZMATLIT, 2003. - 784 p.
2. Barron, J.L. Performance of optical flow techniques / J.L. Barron, D.J. Fleet, S.S. Beauchemin // International Journal of Computer Vision. - 1994. - Vol. 12, N 1. - P. 43-77.
3. Beauchemin, S.S. The computation of optical flow / S.S. Beauchemin, J.L. Barron // ACM Computing Surveys. -1995. - Vol. 27, N 3. - P. 433-466.
4. Horn, B.K.P. Determining optical flow / B.K.P. Horn, B.G. Schunck // Artificial Intelligence. - 1981. - Vol. 17, N 1-3. -P. 185-203.
5. Lucas, B.D. An Iterative Image Registration Technique with an Application to Stereo Vision / B.D. Lucas, T. Kanade // Proceedings of Imaging Understanding Workshop. - 1981. -P. 121-130.
6. Lucas, B.D. Generalized image matching by the method of differences / Lucas Bruce David. - Carnegie Mellon University, 1984. - 167 p.
7. Fennema, C.L. Velocity determination in scenes containing several moving objects / C.L. Fennema, W.B. Thompson // Computer Graphics and Image Processing. - 1979. - Vol. 9, N 4. -P. 301-315.
8. Fleet, D.J. Phase-based disparity measurement / D.J. Fleet, A.D. Jepson, M.R.M. Jenkin // CVGIP: Image Understanding. -1991. - Vol. 53, N 2. - P. 198-210.
94. Matos P.F.P. de, Nowell D. Experimental and numerical investigation of thickness effects in plasticity-induced fatigue crack closure // International Journal of Fatigue. - 2009.
95. Estimation of mixed-mode stress intensity factors using digital image correlation and an interaction integral / J. Rethore [et al.] // International Journal of Fracture. - 2005.
96. Evaluating J-integral from displacement fields measured by digital image correlation / S. Yoneyama [et al.] // Strain. - 2014.
97. Панин С.В., Титков В.В., Любутин П.С. Инкремент-ный подход к определению перемещений фрагментов изображений при построении векторных полей // Автометрия. -2014. - Т. 50, № 2. - С. 39-49.
98. Панин С.В., Титков В.В., Любутин П.С. Сглаживание векторных полей с использованием поверхности Безье при оценке деформации методом корреляции цифровых изображений // Автометрия. - 2014. - Т. 50, № 1. - С. 74-81.
99. Панин С.В., Титков В.В., Любутин П.С. Исследование эффективности алгоритмов фильтрации векторных полей при оценке деформации материалов методом корреляции цифровых изображений // Автометрия. - 2013. - Т. 49, № 2. -С. 57-67.
100. Любутин П.С., Панин С.В. Измерение деформации на мезоуровне путем анализа оптических изображений поверхности нагруженных твердых тел // Прикладная механика и техническая физика. - 2006. - Т. 47, № 6. - С. 158-164.
101. Березин И.С., Жидков Н.П. Методы вычислений: учеб. пособие для вузов. - 2-е изд. - Т. 1. - М.: ГИФМЛ, 1962. - 464 с.
102. Vavrik D., Jandejsek I. Experimental evaluation of contour J integral and energy dissipated in the fracture process zone // Engineering Fracture Mechanics. - 2014. - Vol. 129. - P. 14-25.
9. Adelson, E.H. Spatiotemporal energy models for the perception of motion / E.H. Adelson, J.R. Bergen // Journal of the Optical Society of America A. - 1985. - Vol. 2, N 2. - P. 284.
10. Heeger, D.J. Optical flow using spatiotemporal filters / D.J. Heeger // International Journal of Computer Vision. - 1988. -Vol. 1, N 4. - P. 279-302.
11. Watson, A.B. A look at motion in the frequency domain. Vol. 1983 / A.B. Watson, A.J. Ahumada. - 1983.
12. Fleet, D.J. The early processing of spatio-temporal visual information : thesis (M.Sc.) / Fleet David J. - University of Toronto, 1984. - 136 p.
13. Schreier, H. Image Correlation for Shape, Motion and Deformation Measurements / H. Schreier, J.-J. Orteu, M.A. Sutton. -Boston, MA: Springer US, 2009. - 321 p.
14. Obzor blochnykh metodov otsenki dvizheniia v tsifrovykh video signalakh / S.V. Grishin (i dr.) [Overview of motion estimation block methods in digital video signals / S.V. Grishin [and others] Pro-grammnye sistemy i instrumenty. Tematicheskii sbornik. - M.: Izd-vo fakul'teta VMiKMGU, 2008. - pp. 50-62.
15. Zhu, S. A new diamond search algorithm for fast block-matching motion estimation / S. Zhu, K.K. Ma // IEEE Transactions on Image Processing. - 2000. - Vol. 9, N 2. - P. 287-290.
16. Li, W. Successive Elimination Algorithm for Motion Estimation / W. Li, E. Salari // IEEE Transactions on Image Processing. - 1995. - Vol. 4, N 1. - P. 105-107.
17. Schreier, H.W. Systematic errors in digital image correlation caused by intensity interpolation / H.W. Schreier // Optical Engineering. - 2000. - Vol. 39, N 11. - P. 2915-2921.
18. Wu, Q.X. A correlation-relaxation-labeling framework for computing optical flow-template matching from a new perspective / Q.X. Wu // IEEE Transactions on Pattern Analysis and Machine Intelligence. - 1995. - Vol. 17, N 9. - P. 843-853.
19. Black, M.J. The Robust Estimation of Multiple Motions: Parametric and Piecewise-Smooth Flow Fields / M.J. Black, P. Anandan // Computer Vision and Image Understanding. -1996. - Vol. 63, N 1. - P. 75-104.
20. Voskoboinikov, Y.E. Comparative analysis of two classes of image filtration algorithms / Y.E. Voskoboinikov, A.V. Go-chakov // Optoelectronics, Instrumentation and Data Processing. -
2011. - Vol. 47, N 1. - P. 12-22.
21. Voskoboinikov, Y.E. Construction of wavelet filtration algorithms with two-parameter threshold functions / Y.E. Vosko-boinikov, A. V. Gochakov // Optoelectronics, Instrumentation and Data Processing. - 2012. - Vol. 48, N 1. - P. 9-17.
22. Otte, M. Estimation of optical flow based on higherorder spatiotemporal derivatives in interlaced and non-interlaced image sequences / M. Otte, H.-H. Nagel // Artificial Intelligence. -1995. - Vol. 78, N 1-2. - P. 5-43.
23. Nagel, H.-H. Displacement vectors derived from second-order intensity variations in image sequences / H.-H. Nagel // Computer Vision, Graphics, and Image Processing. - 1983. -Vol. 21, N 1. - P. 85-117.
24. Nagel, H.-H. An Investigation of Smoothness Constraints for the Estimation of Displacement Vector Fields from Image Sequences / H.-H. Nagel, W. Enkelmann // IEEE Transactions on Pattern Analysis and Machine Intelligence. - 1986. -Vol. PAMI-8, N 5. - P. 565-593.
25. Wang, C.-M. Estimating Optical Flow by Integrating Multi-Frame Information / C.-M. Wang, K.-C. Fan, C.-T. Wang // Journal of Information Science and Engineering. - 2008. -Vol. 24, N 6. - P. 1719-1731.
26. Irani, M. Multi-frame optical flow estimation using subspace constraints / M. Irani // Proceedings of the Seventh IEEE International Conference on Computer Vision. - IEEE, 1999. -Vol. 1. - P. 626-633.
27. Gao, W. Parallel computing in experimental mechanics and optical measurement: A review / W. Gao, Q. Kemao // Optics and Lasers in Engineering. - 2012. - Vol. 50. - P. 608-617.
28. High accuracy digital image correlation powered by GPU-based parallel computing / L. Zhang [et al.] // Optics and Lasers in Engineering. - 2015. - Vol. 69. - P. 7-12.
29. Shao, X. Noise robustness and parallel computation of the inverse compositional Gauss-Newton algorithm in digital image correlation / X. Shao, X. Dai, X. He // Optics and Lasers in Engineering. - 2015. - Vol. 71. - P. 9-19.
30. A multi-resolution approach for massively-parallel hardware-friendly optical flow estimation / F. Barranco [et al.] // Journal of Visual Communication and Image Representation. -
2012. - Vol. 23, N 8. - P. 1272-1283.
31. Marzat, J. Real-time dense and accurate parallel optical flow using CUDA / J. Marzat, Y. Dumortier, A. Ducrot // 7th International Conference WSCG. - 2009. - P. 105-111.
32. SimpleFlow: A Non-iterative, Sublinear Optical Flow Algorithm / M. Tao [et al.] // Computer Graphics Forum. - 2012. -Vol. 31, 2pt1. - P. 345-353.
33. Rav-Acha, A. Lucas-Kanade without Iterative Warping / A. Rav-Acha, S. Peleg // 2006 International Conference on Image Processing. - IEEE, 2006. - P. 1097-1100.
34. Braspenning, R.A. True-motion estimation using feature correspondences / R.A. Braspenning, G. de Haan // eds. S. Pancha-nathan, B. Vasudev. - 2004. - P. 396.
35. True-motion estimation with 3-D recursive search block matching / G. de Haan [et al.] // IEEE Transactions on Circuits and Systems for Video Technology. - 1993. - Vol. 3, N 5. - P. 368-379, 388.
36. Optimization of Hierarchical 3DRS Motion Estimators for Picture Rate Conversion / A. Heinrich [et al.] // IEEE Journal of Selected Topics in Signal Processing. - 2011. - Vol. 5, N 2. -P. 262-274.
37. Ahn, T.G. Fast full-search motion estimation based on multilevel successive elimination algorithm / T.G. Ahn, Y.H. Moon, J.H. Kim // IEEE Transactions on Circuits and Systems for Video Technology. - 2004. - Vol. 14, N 11. - P. 1265-1269.
38. Lopes, F. Hierarchical motion estimation with spatial transforms / F. Lopes, M. Ghanbari // Proceedings 2000 International Conference on Image Processing (Cat. No.00CH37101). -IEEE, 2000. - Vol. 2. - P. 558-561.
39. Application of an optimized digital correlation method to planar deformation analysis / M. Sutton [et al.] // Image and Vision Computing. - 1986. - Vol. 4, N 3. - P. 143-150.
40. Applications of Digital-Image Correlation techniques to experimental mechanics / T.C. Chu [et al.] // Experimental Mechanics. - 1985. - Vol. 25, N 3. - P. 232-244.
41. Determination of displacements using an improved digital correlation method / M. Sutton [et al.] // Image and Vision Computing. - 1983. - Vol. 1, N 3. - P. 133-139.
42. Digital image correlation using Newton-Raphson method of partial differential correction / H.A. Bruck [et al.] // Experimental Mechanics. - 1989. - Vol. 29, N 3. - P. 261-267.
43. Advances in Two-Dimensional and Three-Dimensional Computer Vision / M.A. Sutton [et al.] // Photomechanics. -Springer Berlin Heidelberg, 2000. - P. 323-372.
44. Hild, F. Digital Image Correlation: from Displacement Measurement to Identification of Elastic Properties - a Review / F. Hild, S. Roux // Strain. - 2006. - Vol. 42, N 2. - P. 69-80.
45. Besnard, G. "Finite-element" displacement fields analysis from digital images: Application to Portevin-Le Châtelier bands / G. Besnard, F. Hild, S. Roux // Experimental Mechanics. -2006. - Vol. 46, N 6. - P. 789-803.
46. Roux, S. Stress intensity factor measurements from digital image correlation: Post-processing and integrated approaches / S. Roux, F. Hild // International Journal of Fracture. - 2006. -Vol. 140, N 1-4. - P. 141-157.
47. Multiscale displacement field measurements of compressed mineral-wool samples by digital image correlation / F. Hild [et al.] // Applied Optics. - 2002. - Vol. 41, N 32. - P. 6815.
48. Two-dimensional digital image correlation for in-plane displacement and strain measurement: a review / B. Pan [et al.] // Measurement Science and Technology. - 2009. - Vol. 20, N 6. -P. 062001.
49. Performance of sub-pixel registration algorithms in digital image correlation / P. Bing [et al.] // Measurement Science and Technology. - 2006. - Vol. 17, N 6. - P. 1615-1621.
50. Study on subset size selection in digital image correlation for speckle patterns / B. Pan [et al.] // Optics Express. -2008. - Vol. 16, N 10. - P. 7037.
51. Digital image correlation using iterative least squares and pointwise least squares for displacement field and strain field measurements / B. Pan [et al.] // Optics and Lasers in Engineering. - 2009. - Vol. 47, N 7-8. - P. 865-874.
52. Pan, B. Recent Progress in Digital Image Correlation / B. Pan // Experimental Mechanics. - 2011. - Vol. 51, N 7. -P. 1223-1235.
53. Mekhanika materialov. Metody i sredstva eksperimen-tal'nykh issledovanii [Mechanics of materials. Methods and means of experimental research] V.E. Vil'deman [et al.]. - Perm': Izd-vo PNIPU, 2011. - 164 p.].
54. Investigation of fatigue crack closure using multiscale image correlation experiments / J. Carroll [et al.] // Engineering Fracture Mechanics. - 2009. - Vol. 76, N 15. - P. 2384-2398.
55. Mathieu, F. Identification of a crack propagation law by digital image correlation / F. Mathieu, F. Hild, S. Roux // International Journal of Fatigue. - 2012. - Vol. 36, N 1. - P. 146-154.
56. High resolution digital image correlation measurements of strain accumulation in fatigue crack growth / J.D. Carroll [et al.] // International Journal of Fatigue. - 2013.
57. Zhang, R. Measurement of mixed-mode stress intensity factors using digital image correlation method / R. Zhang, L. He // Optics and Lasers in Engineering. - 2012.
58. An approach to calculate the J-integral by digital image correlation displacement field measurement / T.H. Becker [et al.] // Fatigue and Fracture of Engineering Materials and Structures. - 2012.
59. Analysis of standard fracture toughness test based on digital image correlation data / I. Jandejsek [et al.] // Engineering Fracture Mechanics. - 2017.
60. A J-integral approach using digital image correlation for evaluating stress intensity factors in fatigue cracks with closure effects / G.L.G. Gonzales [et al.] // Theoretical and Applied Fracture Mechanics. - 2017.
61. Yoneyama, S. Automatic evaluation of mixed-mode stress intensity factors utilizing digital image correlation / S. Yoneyama, Y. Morimoto, M. Takashi // Strain. - 2006.
62. Investigation of thermal effects on fatigue crack closure using multiscale digital image correlation experiments / M.C. Cas-person [et al.] // International Journal of Fatigue. - 2014.
63. Panin, S.V. Rigid body deformation estimation by surface images / S.V. Panin, V.I. Syryamkin, P.S. Lyubutin // Avtometriya. - 2005. - Vol. 41. - No. 2. - P. 44-58.
64. Giachetti, A. Matching techniques to compute image motion // Image and Vision Computing. - 2000. - Vol. 18, N 3. -P. 247-260.
65. Tong, W. An evaluation of digital image correlation criteria for strain mapping applications // Strain. - 2005. - Vol. 41. -No. 4. - P. 167-175.
66. Full-field strain measurement using a two-dimensional Savitzky-Golay digital differentiator in digital image correlation / B. Pan [et al.] // Optical Engineering. - 2007. - Vol. 46, N 3. -P. 33601-33610.
67. Cofaru C., Philips W., Van Paepegem W. Improved Newton-Raphson digital image correlation method for full-field displacement and strain calculation // Applied Optics. - 2010. -Vol. 49, N 33. - P. 6472.
68. Fundamental role of crystal structure curvature in plasticity and strength of solids / V.E. Panin [et al.] // Physical Me-somechanics. - 2015. - Vol. 18, N 2. - P. 89-99.
69. Synergy of crack closure, near-tip residual stress and crack-tip blunting in crack growth under periodic overloads -A fractographic study / R. Sunder [et al.] // International Journal of Fatigue. - 2016. - Vol. 93. - P. 18-29.
70. Song, K.Y. Texture crack detection / K.Y. Song, M. Petrou, J. Kittler // Machine Vision and Applications. - 1995. -Vol. 8, N 1. - P. 63-75.
71. Kirschke, K.R. Histogram-Based Approach for Automated Pavement-Crack Sensing / K.R. Kirschke, S.A. Velinsky // Journal of Transportation Engineering. - 1992. - Vol. 118, N 5. -P. 700-710.
72. CrackTree: Automatic crack detection from pavement images / Q. Zou [et al.] // Pattern Recognition Letters. - 2012. -Vol. 33, N 3. - P. 227-238.
73. Sun, L. Weighted Neighborhood Pixels Segmentation Method for Automated Detection of Cracks on Pavement Surface Images / L. Sun, M. Kamaliardakani, Y. Zhang // Journal of Computing in Civil Engineering. - 2016. - Vol. 30, N 2. -P. 04015021.
74. Yamaguchi, T. Fast crack detection method for large-size concrete surface images using percolation-based image processing / T. Yamaguchi, S. Hashimoto // Machine Vision and Applications. - 2010. - Vol. 21, N 5. - P. 797-809.
75. Iyer, S. A robust approach for automatic detection and segmentation of cracks in underground pipeline images / S. Iyer, S.K. Sinha // Image and Vision Computing. - 2005. - Vol. 23, N 10. - P. 921-933.
76. Lins, R.G. Automatic Crack Detection and Measurement Based on Image Analysis / R.G. Lins, S.N. Givigi // IEEE Transactions on Instrumentation and Measurement. - 2016. - Vol. 65, N 3. - P. 583-590.
77. A Database and Evaluation Methodology for Optical Flow / S. Baker [et al.] // International Journal of Computer Vision. - 2011. - Vol. 92, N 1. - P. 1-31.
78. Panin, S.V. Determination of optical flow by complex algorithm on weighted full search (in russian) / S.V. Panin, V.O. Chemezov, P.S. Lyubutin // Computational technologies. -2017. - V. 22, N. 4. - P. 80-94.
79. Gonzalez, R.C. Digital Image Processing (3rd Edition) / R.C. Gonzalez, R.E. Woods. - Prentice Hall, 2008. - 954 p.
80. Glantz, R. Skeletons in the Framework of Graph Pyramids / R. Glantz, W.G. Kropatsch // Visual Form 2001. Lecture Notes in Computer Science / eds. C. Arcelli, L.P. Cordella, G.S. di Baja. - Springer Berlin Heidelberg, 2001. - P. 186-195.
81. Panin, S.V. Selection of parameters of the three-dimensional recursive search algorithm in constructing displacement vector fields with the use of the hierarchical approach / S.V. Panin, V.V. Titkov, P.S. Lyubutin // Optoelectronics, Instrumentation and Data Processing. - 2015. - Vol. 51, N 2. - P. 124-133.
82. Nowell, D. Application of digital image correlation to the investigation of crack closure following overloads / D. Nowell, P.F.P. De Matos // Procedia Engineering. - 2010.
83. Yusof, F. Effect of overload on crack closure in thick and thin specimens via digital image correlation / F. Yusof, P. Lopez-Crespo, P.J. Withers // International Journal of Fatigue. -2013. - Vol. 56. - P. 17-24.
84. Lu, Z. Small time scale fatigue crack growth analysis / Z. Lu, Y. Liu // International Journal of Fatigue. - 2010.
85. Toribio, J. Finite-deformation analysis of the crack-tip fields under cyclic loading / J. Toribio, V. Kharin // International Journal of Solids and Structures. - 2009.
86. Near-tip strain evolution under cyclic loading: In situ experimental observation and numerical modelling / J. Tong [et al.] // International Journal of Fatigue. - 2015.
87. Zhang, W. Plastic zone size estimation under cyclic loadings using in situ optical microscopy fatigue testing / W. Zhang, Y. Liu // Fatigue and Fracture of Engineering Materials and Structures. - 2011.
88. Cherepanov, G.P. Crack propagation in continuous media / G.P. Cherepanov // Journal of Applied Mathematics and Mechanics. - 1967. - Vol. 31, N 3. - P. 503-512.
89. Rice, J.R. A Path Independent Integral and the Approximate Analysis of Strain Concentration by Notches and Cracks / J.R. Rice // Journal of Applied Mechanics. - 1968.
90. J-Integral Calculation by Finite Element Processing of Measured Full-Field Surface Displacements / S.M. Barhli [et al.] // Experimental Mechanics. - 2017. - Vol. 57, N 6. - P. 997-1009.
91. McNeill, S.R. Estimation of stress intensity factor by digital image correlation / S.R. McNeill, W.H. Peters, M.A. Sutton // Engineering Fracture Mechanics. - 1987. - Vol. 28, N 1. - P. 101-112.
92. Lorenzino, P. Application of digital image correlation (DIC) in resonance machines for measuring fatigue crack growth / P. Lorenzino, G. Beretta, A. Navarro // Frattura ed Integrita Strut-turale. - 2014.
93. Hos, Y. Using digital image correlation to determine mixed mode stress intensity factors in fatigue cracks / Y. Hos, M. Vormwald, J.L.F. Freire // Proceedings of COTEQ 2015. - 2015.
94. Matos, P.F.P. de. Experimental and numerical investigation of thickness effects in plasticity-induced fatigue crack closure / P.F.P. de Matos, D. Nowell // International Journal of Fatigue. - 2009.
95. Estimation of mixed-mode stress intensity factors using digital image correlation and an interaction integral / J. Rethore [et al.] // International Journal of Fracture. - 2005.
96. Evaluating J-integral from displacement fields measured by digital image correlation / S. Yoneyama [et al.] // Strain. - 2014.
97. Panin, S.V. Incremental approach to determination of image fragment displacements during vector field construction / S.V. Panin, V.V. Titkov, P.S. Lyubutin // Optoelectronics, Instrumentation and Data Processing. - 2014. - Vol. 50, N 2. - P. 139-147.
98. Panin, S.V. Smoothing of vector fields by using the Bezier surface for strain estimation by the method of digital image correlation / S.V. Panin, V.V. Titkov, P.S. Lyubutin // Optoelectronics, Instrumentation and Data Processing. - 2014. - Vol. 50, N 1. - P. 61-67.
99. Panin, S. V. Efficiency of vector field filtration algorithms in estimating material strain by the method of digital image correlation / S.V. Panin, V.V. Titkov, P.S. Lyubutin // Optoelec-
tronics, Instrumentation and Data Processing. - 2013. - Vol. 49, N 2. - P. 155-163.
100. Lyubutin, P.S. Mesoscale measurement of strains by analyzing optical images of the surface of loaded solids / P.S. Lyubutin, S.V. Panin // Journal of Applied Mechanics and Technical Physics. - 2006. - Vol. 47, N 6. - P. 905-910.
101. Berezin, I.S. Metody vychislenij: utchebnoe posobie dlya vuzov 2-e izd. T. 1 [Computing Methods 2-nd ed. Vol. 1]. Moscow, GIFML, 1962, 464 p.
102. Vavrik, D. Experimental evaluation of contour J integral and energy dissipated in the fracture process zone. Engineering Fracture Mechanics, 2014, vol. 129, pp. 14-25.
