Digital processing of the results of optoelectronic measurements. The photoelasticity method and its application for determination of coefficients of the multiparameter asymptotic Williams expansion of the stress field

Full Text

Введение. О методе цифровой фотоупругости. В современной фундаментальной механике деформируемого твердого тела, в механике разрушения в частности, и ее промышленных приложениях первоочередной задачей, возникающей в теоретическом анализе, при расчете и конструировании инженерных сооружений, является необходимость учета нелинейностей различного вида и сингулярных особенностей, вызванных частичным разрушением у концентраторов напряжений в виде трещин, острых углов, вырезов и включений. Эти задачи не решаются инженерными методами с помощью обычных нормативных методик и требуют разработки более сложных подходов к их решению с использованием математического аппарата механики разрушения, сложных программных комплексов. Высокоточные экспериментальные методики часто оказываются единственным средством верификации различных теорий механики деформируемого твердого тела и механики разрушения как одной из ее составных частей. Весьма эффективными следует признать экспериментальные поляризационно-оптические методы определения напряжений. Поэтому поляризационно-оптические методы исследования полей напряжений и перемещений в механике деформируемого твердого тела используются как в исследованиях отечественных ученых, так и за рубежом [1-3]. Экспериментальной оценке полей напряжений вблизи острых вырезов и трещин в нашей стране и за рубежом посвящается значительное количество публикаций, среди которых можно отметить исследования [4-8]. В целом можно выделить три основных направления развития интерференционно-оптических методов в современной механике деформируемого твердого тела. 1. Автоматизация обработки данных фотоупругих измерений и компьютерные программы и коды, автоматизирующие получение и обработку экспериментальной информации, полученной с помощью интерференционно-оптических методов (метода фотоупругих измерений) [9-11]. 2. Развитие методов обработки цифровых изображений, получаемых из интерференционно-оптических методов (метода фотоупругости) [12-17]. 3. Приложения метода фотоупругости к механике разрушения и многопараметрическому описанию полей у вершины трещины; решению контактных задач, приложениям фотоупругости в смежных областях (задачи механики разрушения в биомеханике, медицине) [18-30]. В настоящей работе метод фотоупругости используется для вычисления коэффициентов многопараметрического асимптотического поля напряжений у вершины трещины. В последние годы в механике хрупкого разрушения сложилось четкое и ясное понимание необходимости удержания высших приближений в полном асимптотическом разложении поля напряжений М. Уильямса у вершины трещины [31-35]. В этих работах показано, что чем больше расстояние от кончика трещины, тем больше слагаемых необходимо удерживать в асимптотическом представлении М. Уильямса поля напряжений. Целью работы [33] являлось описание полей напряжений и перемещений в окрестности вершины трещины и на ее удалении в ряде эксперименталь718 Цифровая обработка результатов оптоэлектронных измерений ных образцов для определения параметров механики разрушения. В статье предложены образцы с новой геометрией, позволяющей рассмотреть широкий диапазон видов нагружения и оценить размеры областей нелинейного деформирования, возникающих у вершины трещины. Предложены надрезанные образцы призматической формы для испытаний на четырехточечный изгиб. Наряду с проведенными экспериментальными исследованиями выполнено численное моделирование этих же образцов. В асимптотических разложениях поля напряжений удерживались сингулярное слагаемое, T -напряжение и высшие приближения (несингулярные слагаемые). В статье выполнено сравнение экспериментальных результатов, численного анализа и предложенных асимптотических разложений в зависимости от расстояния до кончика трещины и от количества удерживаемых слагаемых в асимптотических разложениях. В [34] представлен анализ полей напряжений и перемещений у вершины трещины в образцах с краевой трещиной, предназначенных для испытаний на расклинивание. В англоязычной литературе такие образцы носят название WST-образцов. Выполнен многопараметрический анализ поля напряжений, для чего в полном асимптотическом разложении М. Уильямса удерживались несингулярные слагаемые. Изучено влияние краевых условий (условий закрепления образца) на коэффициент интенсивности напряжений и на T -напряжения. В целом проведенный анализ позволяет учесть до 12 слагаемых в асимптотическом разложении. Коэффициенты высших приближений асимптотических разложений определялись из экспериментальных данных. С помощью определенных коэффициентов асимптотических разложений было реконструировано поле напряжений и далее выполнялось сравнение с вычислениями, базирующимися на методе конечного элемента. Проведенное сравнение ясно показывает, что для длинных трещин и для больших расстояний от вершины трещины аккуратное описание поля напряжений требует удержания нескольких слагаемых в полном асимптотическом разложении поля напряжений у вершины трещины. В [35] получены аналитические выражения для амплитудных коэффициентов полного асимптотического разложения поля напряжений у вершин двух коллинеарных трещин конечной длины в бесконечной пластине, находящейся в трех условиях: чистое растяжение; чистый поперечный сдвиг; смешанное нагружение. Показано, что наряду с главным членом асимптотического разложения и T -напряжением высшие приближения в полном асимптотическом разложении являются важными параметрами механики разрушения, учет которых необходим для аккуратной оценки полей напряжений в окрестности вершины дефекта и, следовательно, для аккуратной оценки несущей способности и долговечности образца. Аналитические выражения коэффициентов асимптотического разложения поля напряжений для любого наперед заданного числа удерживаемых слагаемых доступны лишь для пластины с одной трещиной конечной длины, поэтому важно расширить класс образцов, для которых имеются аналитические зависимости амплитудных множителей от приложенной нагрузки и геометрических параметров образца. Полученное в работе асимптотическое решение для бесконечной пластины с двумя трещинами равной длины может быть использовано в следующих случаях: 719 С т е п а н о в а Л. В., Д о л г и х В. С. 1) для проверки численной обработки экспериментальных данных, полученных с помощью поляризационно-оптических методов механики деформируемого твердого тела для образцов с трещинами; 2) для проверки численных алгоритмов, активно разрабатываемых в настоящее время для вычисления амплитудных коэффициентов многопараметрического поля напряжений. Особенно важное значение высших приближений в полном асимптотическом разложении М. Уильямса приобретает при обработке экспериментальной информации вне зависимости от метода исследования полей у вершины трещины. Это остается справедливым как для классического метода фотоупругости, так и для новой техники, вошедшей в употребление - метода корреляции цифровых изображений. Поэтому в настоящей работе с помощью метода фотоупругости наряду с коэффициентами интенсивности напряжений и T -напряжениями определены коэффициенты высших приближений асимптотических разложений поля напряжений у вершины трещины в пластине с двумя коллинеарными трещинами одинаковой и различной длины. 1. Экспериментальная процедура и экспериментальное определение полей напряжений у вершины трещины методом фотоупругости. В настоящей работе была проведена серия экспериментов с помощью метода фотоупругости на проекционно-поляризационной установке ППУ-7 на различных образцах с трещинами: пластины с двумя коллинеарными горизонтальными и наклонными трещинами одинаковой и различной длины. Образцы были изготовлены из эпоксидной смолы ЭД20. Образцы с трещинами различной конфигурации испытывались при разной амплитуде нагружения. На рис. 1 изображены картины изохроматических полос в пластинах с двумя коллинеарными трещинами одинаковой длины, находящихся под действием растягивающей нагрузки в 120, 125, 145 и 175 кг. На рис. 2 показаны картины изохроматических полос в пластинах с двумя коллинеарными трещинами различной длины. Интерференционная картина изохроматических полос с наклонными трещинами показана на рис. 3. Образец с наклонной трещиной является удобной моделью для экспериментального изучения смешанного нагружения, когда присутствуют обе моды нагружения. Следует отметить, что корректная интерпретация интерференционных картин, получаемых при просвечивании в поляроскопе, вызывает значительные сложности. В настоящей работы полученные картины изохроматичсеких полос будут использованы для определения коэффициентов полного асимптотического разложения М. Уильямса поля напряжений у вершин двух трещин в линейно-упругой пластине (рис. 2). Для определения цены изохроматической полосы (оптической постоянной материала) fs использовался диск из того же материала (эпоксидная смола ЭД20), нагруженный поочередно силами 60, 70, 120, 140, 180 и 210 кг. Первый шаг любого экспериментального исследования с помощью метода фотоупругости предполагает определение оптической константы материала. Процесс определения цены изохроматической полосы (оптической константы материала) называют тарировкой. Для каждого шага нагружения фиксируется порядок полосы и приложенная нагрузка. Проводится измерение диаметра 720 Цифровая обработка результатов оптоэлектронных измерений Рис. 1. Картина изохроматических полос в пластинах с двумя коллинеарными трещинами одинаковой длины [Figure 1. Isochromatic fringes for plates with two collinear cracks of equal lengths] Рис. 2. Картина изохроматических полос в пластинах с двумя коллинеарными трещинами различной длины [Figure 2. Isochromatic fringes for plates with two collinear cracks of different lengths] и толщины диска, а также поиск координат точки, лежащей на выбранной изохроматической полосе. Качественная связь между порядком полос и разностью главных напряжений устанавливается по следующей формуле [1, 2], представляющей собой математическое выражение закона Вертгейма [1]: σ1 - σ2 = N fs /h, (1) где N - порядок полосы, h - толщина образца, а σ1 , σ2 - главные напряжения. Отсюда следует, что fs выражает значение разности главных напряжений при заданном номере изохроматической полосы и толщине образца. Это значение и называют ценой полосы материала. Определение значения fs для данного материала осуществляется из экспериментов на образцах, для которых известно теоретическое решение и напряжения могут быть вычислены. Для определения fs были проведены эксперименты на сжатие кругового дис721 С т е п а н о в а Л. В., Д о л г и х В. С. Рис. 3. Картины изохроматических полос в пластинах с двумя наклонными трещинами и пластине с двумя боковыми надрезами [Figures 3. Isochromatic fringes for plates with two inclined cracks of equal lengths and two edge cracks] ка парой сосредоточенных сил, направленных строго в вертикальной плоскости при условии, что начало координат совпадает с центром диска радиуса R. Для выбранной полосы порядка N , которая проходит через точку (x1 , 0) диаметра диска, совпадающего с горизонтальном направлением оси абсцисс заданной системы координат, и для нагружения силами P , при просвечивании диска монохроматическим светом будет иметь место формула fs = |σ22 - σ11 |h/N, так как в точках диаметра, совпадающего с осью абсцисс, напряжения σ11 и σ22 - главные. Вычисленное значение оптической константы для рассматриваемого материала - fσ = 18.23 кг/см. 2. Многопараметрическое разложение М. Уильямса поля напряжений у вершины трещины. Асимптотическое разложение поля напряжений в окрестности вершины трещины в изотропном линейно упругом материала имеет вид 2 ∞ m,ij am (θ)rk/2-1 , k fk σij (r, θ) = m=1 k=-∞ где σij - компоненты тензора напряжений; r, θ - полярные координаты с полюсом в вершине трещины; индекс m отвечает моде нагружения: нормальному отрыву соответствует значение m = 1, а поперечному сдвигу - m = 2; am k - амплитудные (масштабные) множители, зависящие от системы приложенных сил и геометрии испытываемого образца; fkm,ij (θ) - угловые распределения компонент тензора напряжений, определяемые в ходе решения краевых задач [36]: fk1,11 (θ) = k (2 + k/2 + (-1)k ) cos(k/2 - 1)θ - (k/2 - 1) cos(k/2 - 3)θ /2, fk1,22 (θ) = k (2 - k/2 - (-1)k ) cos(k/2 - 1)θ + (k/2 - 1) cos(k/2 - 3)θ /2, fk1,12 (θ) = k -(k/2 + (-1)k ) sin(k/2 - 1)θ + (k/2 - 1) sin(k/2 - 3)θ /2, 722 Цифровая обработка результатов оптоэлектронных измерений fk2,11 (θ) = -k (2 + k/2 - (-1)k ) sin(k/2 - 1)θ - (k/2 - 1) sin(k/2 - 3)θ /2, fk2,22 (θ) = -k (2 - k/2 + (-1)k ) sin(k/2 - 1)θ + (k/2 - 1) sin(k/2 - 3)θ /2, fk2,12 (θ) = k -(k/2 - (-1)k ) cos(k/2 - 1)θ + (k/2 - 1) cos(k/2 - 3)θ /2. Поле перемещений в окрестности вершины трещины может быть представлено в форме [37, 38]: 2 ∞ ui (r, θ) = m=1 k=-∞ 1 m m,i a u ˜ (θ)rk/2 , µ k k где u ˜k1,1 (θ) = (κ + k/2 + (-1)k ) cos(θ/2) - (k/2) cos(k/2 - 2)θ, u ˜k1,2 (θ) = (κ - k/2 - (-1)k ) sin(θ/2) + (k/2) sin(k/2 - 2)θ, u ˜k2,1 (θ) = -(κ + k/2 - (-1)k ) sin(θ/2) + (k/2) sin(k/2 - 2)θ, u ˜k2,2 (θ) = (κ - k/2 + (-1)k ) cos(θ/2) + (k/2) cos(k/2 - 2)θ, κ - константа плоской задачи теории упругости (κ = 3 - 4ν для случая плоского деформированного состояния, κ = (3 - ν)/(1 + ν) для плоского напряженного состояния); µ = E/(2(1 + ν)) - модуль сдвига; E и ν - модуль Юнга и коэффициент Пуассона соответственно. Для каждого образца с трещиной и системы нагрузок существует своя система коэффициентов (масштабных, амплитудных) множителей akm . В большинстве инженерных приложений в асимптотическом разложении М. Уильямса удерживается только первый или первые два слагаемых, коэффициенты при которых носят название коэффициентов интенсивности напряжений и T -напряжения: KI = √ (1) 2πa11 f1,22 (0), KII = √ (1) 2πa21 f2,12 (0), (2) T = a12 f1,11 (0). Однако в последние годы в механике хрупкого разрушения сложилось ясное понимание необходимости удержания высших приближения в асимптотическом разложении М. Уильямса [31-41]. В этих работах отмечается потребность развития методов и подходов, позволяющих получать аналитические зависимости коэффициентов от геометрии образца и системы приложенных нагрузок, и, следовательно, расширить таблицы, содержащие найденные оценки для коэффициентов интенсивности напряжений для различных образцов, выражениями для коэффициентов высших приближений. Численные значения коэффициентов многопараметрических разложений могут быть найдены с помощью компьютерного имитационного моделирования в пакетах, реализующих метод конечного элемента. Однако значение наличия коэффициентов высших приближений, определенных экспериментально, сложно переоценить. Поэтому в настоящей работе с помощью основного закона фотоупругости и серии экспериментов определены коэффициенты высших приближений полного асимптотического разложения М. Уильямса поля напряжений у вершин двух трещин в пластине с двумя коллинеарными трещинами. 723 С т е п а н о в а Л. В., Д о л г и х В. С. Создан и отлажен комплекс программ для анализа и численной обработки всей совокупности экспериментальной информации. В первой программе данные фотоупругих измерений собраны для выбранного количества точек в исследуемом образце. Алгоритм программы базируется на том факте, что каждый пиксель фотографии имеет значение интенсивности в диапазоне от 0 до 255. Поэтому пиксель, соответствующий значению 0, на фотографии изображается черным цветом, а пиксель со значением 255 изображается белым цветом. Следовательно, пиксели с более низкими значениями интенсивности отвечают более темному цвету. Созданная программа используется для определения значений интенсивности точек, лежащих внутри изохроматических полос. Эта процедура выполняется несколько раз для различных линий в радиальных направлениях. В результате определяется наиболее темная точка каждой изохромы. Вычисленные координаты этих наиболее темных точек используются как входные данные для работы второй программы. Программы объединены в единый комплекс, позволяющий определить поля напряжений в исследуемом образце. В алгоритме реализована возможность учета высших членов в асимптотическом решении задачи (в полном решении М. Уильямса) и применения высших приближений в оптико-механическом законе. Для определения коэффициентов высших приближений был использован переопределенный метод, суть которого описана ниже. Если в асимптотическом разложении М. Уильямса удерживать K слагаемых, отвечающих нормальному отрыву, и M слагаемых, отвечающих поперечному сдвигу, то оптико-механический закон (1) примет вид N fσ h 2 1 2 = (σ11 - σ22 )2 + σ12 = 4 K = - (k) a1k rk/2-1 f1,11 (θ) + M (k) a2k rk/2-1 f2,11 (θ)- k=1 k=1 K M 2 1 k/2-1 (k) 2 k/2-1 (k) ak r f1,22 (θ) - ak r f2,22 (θ) + k=1 k=1 K M 2 (k) (k) + a1k rk/2-1 f1,12 (θ) + a2k rk/2-1 f2,12 (θ) . k=1 k=1 (2) Таким образом, получена система K + M нелинейных алгебраических уравнений относительно масштабных (амплитудных) множителей a11 , a12 , . . . , a1K и a21 , a22 , . . . , a2M . Решение сформулированной системы нелинейных алгебраических уравнений представляет самостоятельный интерес. Для решения системы уравнений (2) был использован переопределенный метод [39], широко применяемый в настоящее время [40-42]. В соответствии с переопределенным методом в рассмотрение вводится функция σ11 - σ22 2 Nm f σ 2 gm = . + (σ12 )2m - 2 h m Очевидно, что если для масштабных множителей a1k , a12 , . . . , a1K и a21 , a22 , . . . , a2M выбраны некоторые начальные приближения, то, в общем случае, введенная 724 Цифровая обработка результатов оптоэлектронных измерений функция gm не будет обращаться в нуль, поскольку априори амплитудные множители неизвестны. Для уточнения значений масштабных множителей a1k , a12 , . . . , a1K и a21 , a22 , . . . , a2M в рамках переопределенного метода вводится итерационная процедура, основанная на разложении функции gm в ряд Тейлора. Разложение функции gm в ряд Тейлора имеет вид (gm )i+1 = (gm )i + ∂gm ∂gm ∂gm (△a11 )i + (△a12 )i + · · · + 1 (△a1K )i + 1 1 ∂a1 ∂a2 ∂aK ∂gm ∂gm ∂gm + (△a21 )i + (△a22 )i + · · · + 2 (△a2M )i , (3) ∂a21 ∂a22 ∂aM где i - номер итерации; разности между амплитудными множителями (△a11 )i , (△a12 )i , . . . , (△a1K )i , (△a21 )i , (△a22 )i , . . . , (△a2L )i определяются выражениями (△a1k )i = a1k(i+1) - a1k(i) , (△a2m )i = a2m(i+1) - a2m(i) , k = 1, 2, . . . K; m = 1, 2, . . . , M. Уточненные значения масштабных множителей a11 , . . . , a1K и a21 , . . . , a2M определяются таким образом, чтобы (gm )i+1 = 0, и в силу этого уравнение (3) приводит к системе линейных алгебраических уравнений: - (gm )i = ∂gm ∂gm ∂gm (△a11 )i + (△a12 )i + · · · + 1 (△a1K )i + 1 1 ∂a1 ∂a2 ∂aK ∂g ∂gm ∂gm m (△a21 )i + (△a22 )i + · · · + 2 (△a2M )i . (4) + 2 2 ∂a1 ∂a2 ∂aM Разрешая полученную систему линейных алгебраических уравнений (4), можно получить новые уточненные значения масштабных множителей a1k , a12 , . . . , a1K и a21 , a22 , . . . , a2M . Следует отметить, что в ходе вычислений легко найти все частные производные, фигурирующие в системе (4): 1 ∂gm ∂σ11 ∂σ22 = (σ11 - σ22 )m - 1 2 ∂aN ∂a1N ∂a1N m + 2 σ12 ∂σ12 ∂a1N m . Частные производные компонент тензора напряжений по амплитудным множителям вычисляются по формулам     ∂σ11      k + k cos( k - 1)θ - k - 1 cos k - 3 θ      2 + (-1) 1  ∂a          2 2 2 2  ∂σ k   k k-2    22 k k k k k 2 = r 2 - (-1) - - 1 θ + - 1 cos - 3 θ cos ∂a1k     2 2 2 2   2          k k k k ∂σ12     k     sin -1 θ+ - 1 sin -3 θ  - (-1) +   1 2 2 2 2 ∂ak 725 С т е п а н о в а Л. В., Д о л г и х В. С. и   ∂σ11        ∂a2k      ∂σ   k k-2 22 =- r 2 2  ∂ak   2     ∂σ12        ∂a2k  k k   2 - (-1)k + sin -1 θ-    2 2  k k 2 + (-1)k - -1 θ+ sin  2 2    k k   - (-1)k - -1 θ+ cos 2 2 k - 1 sin 2 k - 1 sin 2 k - 1 cos 2  k  - 3 θ    2  k -3 θ .  2    k  -3 θ  2 Итерационная процедура завершается при выполнении принятого критерия сходимости расчета. На практике применялись два критерия. Первый критерий основан на оценке значений разности между амплитудными множителями на i-той и i+1-й итерациями. Итерационный процесс завершается, если выполнено неравенство K M (△a1k )2i k=1 (△a2m )2i < 10-6 , + (5) m=1 Проведенные вычисления показали, что если в качестве нулевого приближения в рассматриваемом алгоритме используются значения амплитудных множителей, вычисленные с помощью аналитического решения для бесконечной пластины, ослабленной двумя коллинеарными трещинами, то итерационный процесс быстро сходится и достаточно всего 8-10 итераций для выполнения условия (5). Выбранные начальные приближения для масштабных множителей am k позволяют найти коэффициенты высших слагаемых в асимптотическом разложении М. Уильямса. Второй критерий [39] основан на дополнительной оценке порядка изохроматической полосы N l , где l = 1, 2, . . . , L; L - общее количество анализируемых изохроматических полос. Новые уточненные значения масштабных множителей используются для вычисления порядка изохроматической полосы на l каждой итерации Ntheory и сравниваются с экспериментально определенным l . Итерационный процесс завершаетпорядком изохроматической полосы Nexp ся, если выполнен критерий 1 Q L l=1 l l Ntheory - Nexp 10-6 , где Q - общее количество экспериментальных точек. 3. Результаты вычислений. Принимая в качестве начального приближения для масштабных множителей - коэффициентов полного асимптотического разложения М. Уильямса am k точное аналитическое решение задачи о деформировании бесконечной пластины с двумя коллинеарными трещинами [35], можно получить экспериментально определенные значения амплитудных множителей. Результаты вычислений сведены в таблицу. Принятые обозначения и размеры приведены на рис. 4. На рис. 5 синим цветом показаны выбранные экспериментальные точки на изохроматических полосах с номерами N = 2 и N = 3. Красным цветом показаны построенные с помощью вычисленных коэффициентов линии 726 Цифровая обработка результатов оптоэлектронных измерений k z=a z=c z=b z=d 1 2 3 5 7 9 11 13 15 17 19 21 23 25 27 29 31 33 35 37 39 41 43 45 47 0.5139 -0.250 0.2512 -0.0646 0.0330 -0.0208 0.0145 -0.0109 0.0086 -0.0069 0.0058 -0.0049 0.00426 -0.0037 0.0033 -0.0029 0.0026 -0.0024 0.0022 -0.0020 0.0018 -0.0017 0.0016 -0.0015 0.0014 0.52398 -0.2500 0.27936 -0.04384 0.04700 -0.00820 0.00234 -0.00279 0.001468 -0.00126 0.01027 -0.00068 0.00769 -0.00041 0.00603 -0.00027 0.00489 -0.00018 0.00407 -0.00013 0.000345 -0.00009 0.00298 -0.00007 0.00260 0.52397 -0.250 0.2512 -0.0646 0.0330 -0.0208 0.0145 -0.0109 0.0086 -0.0069 0.0058 -0.0049 0.00426 -0.0037 0.0033 -0.0029 0.0026 -0.0024 0.0022 -0.0020 0.0018 -0.0017 0.0016 -0.0015 0.0014 0.56275 -0.2500 0.31993 -0.02563 0.06089 -0.00065 0.00299 -0.00117 0.001853 0.00122 0.01286 0.00106 0.00958 0.00901 0.00749 0.00766 0.00060 0.00065 0.00504 0.00057 0.00027 0.00050 0.00368 0.00044 0.00321 ∞ Коэффициенты a ˜1k = a1k /σ22 -(k+2)/2 (см ) полного асимптотического разложения М. Уильямса для вершины трещины z ∞ [Coefficients a ˜1k = a1k /σ22 -(k+2)/2 (cm ) of the asymptotic Williams’ series expansion in the vicinity of the crack tip z Рис. 4. Схематичное изображение образца и приложенной нагрузки [Figure 4. Schematic view of cracked specimen and the load applied] Рис. 5. Экспериментальные точки и теоретически построенные картины изохроматических полос (онлайн в цвете) [Figure 5. Experimental points and theoretically reconstracted isochromatic fringe pattern (color online)] 727 С т е п а н о в а Л. В., Д о л г и х В. С. равных значений интенсивности напряжений. Сравнение экспериментальных и теоретических изохроматических полос показывает, что полосы проходят ровно через экспериментальные точки. На рис. 6 показаны теоретически реконструированные изохроматические полосы, построенные у вершины трещины нормального отрыва с помощью двучленного асимптотического разложения интенсивности напряжений (слева), с помощью пятичленного разложения интенсивности напряжений (в центре) и девятичленного асимптотического разложения интенсивности напряжений (справа). Видно, что построенные картины полностью воспроизводят экспериментально полученные интерференционные картины (рис. 1). Из рис. 1 видно, что многопараметрическое асимптотическое разложение М. Уильямса поля напряжений у вершины трещины приводит к более аккуратному описанию поля напряжений. Рис. 6. Теоретически построенные картины изохроматических полос [Figure 6. Theoretically reconstracted isochromatic fringe pattern] Выводы. В работе выполнена серия экспериментов на смешанное нагружение оптически чувствительных образцов с трещинами, надрезами и угловыми вырезами с помощью интерференционно-оптического метода механики деформируемого твердого тела - метода фотоупругости. В результате проведенных экспериментов получены картины интерференционных полос: изохром и изоклин для образцов с двумя горизонтальными и наклонными коллинеарными трещинами одинаковой и различных длин. Выбранные типы образцов на смешанное нагружение позволяют эффективно исследовать смешанное деформирование образцов и моделировать смешанное нагружение образца с трещиной или угловым вырезом для всего диапазона значений параметра смешанности нагружения Me = σ22 (r, θ) 2 arctg π σ12 (r, θ) θ=0 , характеризующего вид нагружения: параметр смешанности нагружения изменяется от 0 (что отвечает поперечному сдвигу) до 1 (что соответствует трещине нормального отрыва). Дана интерпретация изохроматических полос для каждого типа образца. Проведена цифровая обработка результатов интерференционно-оптического метода (метода фотоупругости). В работе впервые экспериментально определены коэффициенты высших приближений полного асимптотического разложения М. Уильямса поля на728 Цифровая обработка результатов оптоэлектронных измерений пряжений в окрестности вершины трещины в пластине с двумя коллинеарными трещинами. С помощью основного закона фотоупругости получена система уравнений относительно коэффициентов разложения М. Уильямса. Данная система линеаризована и решена относительно коэффициентов асимптотического разложения, в котором удерживаются высшие приближения. Конкурирующие интересы. Заявляем, что в отношении авторства и публикации этой статьи конфликта интересов не имеем. Авторский вклад и ответственность. Все авторы принимали участие в разработке концепции статьи и в написании рукописи. Авторы несут полную ответственность за предоставление окончательной рукописи в печать. Окончательная версия рукописи была одобрена всеми авторами.

About the authors

Larisa V Stepanova

Samara National Research University

Email: StepanovaLV2015@yandex.ru
http://orcid.org/0000-0002-6693-3132 Dr. Phys. & Math. Sci.; Professor; Dept. of Dept. Mathematical Modelling in Mechanics 34, Moskovskoye shosse, Samara, 443086, Russian Federation

Vadim S Dolgich

Samara National Research University

Email: Vadek13@rambler.ru
Postgraduate Student; Dept. of Dept. Mathematical Modelling in Mechanics 34, Moskovskoye shosse, Samara, 443086, Russian Federation

References

Supplementary files