Full Text

Введение Наряду с экспериментальными методами построения многопараметрического разложения М. Уильямса у вершины трещины в настоящее время коэффициенты асимптотического разложения находят с помощью вычислительного эксперимента в программных комплексах, реализующих метод конечных элементов [1-14]. В данных работах предложены различные схемы расчета, нацеленные на извлечение параметров разрушения из конечно-элементного анализа. Следует заметить, что в большинстве работ определяются первые два коэффициента многопараметрического разложения М. Уильямса, а именно - коэффициенты интенсивности напряжений и Т-напряжение. Определение коэффициентов высших приближений вызывает математические сложности и имеются отдельные результаты [8-11] для ряда конфигураций тел с трещинами и надрезами. В настоящей статье с помощью конечно-элементного моделирования вычислены коэффициенты многопараметрического асимптотического разложения М. Уильямса поля напряжений у вершины трещины в пластине с центральной трещиной и у вершины надреза в круговом полудиске, предназначенном на испытание на трехточечный изгиб. В разложении М. Уильямса удержаны высшие приближения (пятнадцать слагаемых), что позволяет 1) дать более точную оценку компонент тензоров напряжений, деформаций и вектора перемещений; 2) расширить область доминирования асимптотического представления М. Уильямса. 1. Проблема определения коэффициентов многопараметрического асимптотического разложения М. Уильямса Целью настоящего исследования является численное определение коэффициентов многопараметрического асимптотического разложения М. Уильямса поля напряжений у вершины трещины: (1.1) где индекс принимает значение, равное 1 для нормального отрыва, и равное 2 - для поперечного сдвига; - амплитудные, масштабные множители, отражающие влияние геометрии образца и системы приложенных нагрузок; - угловые распределения компонент тензора напряжений [15, 16]. Универсальные угловые распределения хорошо известны из решения краевых задач и приведены в [15, 16]. Как отмечалось в первой части работы, все разнообразие задач о трещинах проявляется в масштабных, амплитудных множителях , зависящих от конфигурации образца и комбинации приложенных нагрузок. В настоящее время сложилось четкое представление о необходимости удержания высших приближений в асимптотическом разложении (1.1). В целом ряде работ [17-29] показано, что удержание высших приближений в асимптотическом разложении М. Уильямса позволяет расширить область действия асимптотического разложения. Поэтому в последнее время появляются работы, в которых авторы пытаются определить коэффициенты разложения в замкнутой форме, используя известные теоретические решения [15-22]. Так, для бесконечной пластины с центральной трещиной в [15] впервые найдены аналитические выражения для всех коэффициентов асимптотического разложения М. Уильямса: - для трещины нормального отрыва (1.2) - поперечного сдвига (1.3) Целью настоящей работы является численное определение амплитудных коэффициентов многопараметрического асимптотического разложения М. Уильямса для прямоугольной пластины с центральной трещиной. Формулы (1.2) будут использованы ниже для верификации результатов извлечения масштабных множителей из конечно-элементного анализа. Соотношения (1.2) будут служить начальным приближением для определения коэффициентов разложения (1.1) для пластины конечных размеров с помощью процедуры минимизации функционала. 2. Извлечение параметров механики разрушения из конечно-элементного эксперимента: восстановление асимптотического разложения М. Уильямса из конечно-элементных вычислений 2.1. Особенности конечно-элементного эксперимента С целью проверки результатов обработки экспериментальных данных, полученных с помощью метода фотоупругости, была проведена серия конечно-элементных вычислений в программном комплексе Simulia Abaqus. Вычислительный эксперимент был проведен для восстановления коэффициентов многопараметрического асимптотического разложения поля напряжений в окрестности вершины трещины в изотропном линейной упругом материале (1.1). Для проведения вычислений была использована пластина высотой 40 см и шириной 20 см с центральной трещиной длиной 1 см в условиях одноосного растяжения. Следует отметить, что такое соотношение геометрических параметров пластины было выбрано для последующего сравнения коэффициентов многопараметрического разложения М. Уильямса, получаемых из конечно-элементного решения, с коэффициентами разложения М. Уильямса, известными из имеющегося аналитического решения задачи (1.2) и (1.3). Модуль Юнга и коэффициент Пуассона материала пластины в серии вычислений полагались равными 2350 МПа и 0,35 соответственно. В рамках конечно-элементного эксперимента для моделирования окрестности вершины трещины и сингулярности поля напряжений были использованы метод контурного интеграла и сингулярные конечные элементы. В качестве базового конечного элемента использовался восьмиузловой биквадратный конечный элемент CPE8R. Общее количество конечных элементов варьировалось от 12 000 до 44 320 элементов. Для построения детальной сетки выбрано достаточно мелкое разбиение, состоящее из 72 секторов, раствор каждого сектора равен 5 градусам. Внутри круговой области имеются 75 концентрических окружностей. Типичное конечно-элементное разбиение вблизи трещины показано на рис. 1. Рис. 1. Типичная конечно-элементная сетка вблизи двух вершин трещин в пластине с центральной трещиной Fig. 1. Typical mesh containing singular elements near the crack tips Следует также отметить, что проведение вычислительного эксперимента для данного типа геометрии образца связано также с имеющимися многочисленными оценками коэффициентов интенсивности напряжений [30-33]. Наряду с автоматическим вычислением коэффициентов интенсивности напряжений, многофункциональный комплекс Simulia Abaqus дает возможность вычислить также Т-напряжения в окрестности вершины трещины. Однако многочисленные исследования вклада высших приближений в асимптотическое разложение М. Уильямса [15-29] показали необходимость удержания и следующих регулярных слагаемых в (1.1). Таким образом, основной целью данной части исследования является извлечение коэффициентов асимптотического разложения из конечно-элементных расчетов. Результаты конечно-элементных вычислений приведены на рис. 2-16. На рис. 2, 3 представлены полученные распределения интенсивности напряжений и нормальные компоненты тензора напряжений. На рис. 4-16 показаны сравнения аналитического решения задачи для растягиваемой бесконечной пластины с центральной трещиной и конечно-элементное решение, а также аналитические решения, построенные с помощью различного количества удерживаемых слагаемых в асимптотическом решении. Рис. 2. Распределение интенсивности напряжений в окрестности вершины трещины малой длины по сравнению с шириной пластины (конечно-элементное решение) Fig. 2. Distribution of the von Mises equivalent stress in the vicinity of the right crack tip Рис.3. Распределение компонент тензора напряжений и у вершины трещины Fig. 3. Distribution of the stress components and in the vicinity of the crack tip На рис. 4, 5 приведены распределения компоненты тензора напряжений в прямоугольной пластине с центральной трещиной длины, много меньшей по сравнению с шириной пластины, на расстоянии от вершины трещины, где красным цветом изображено аналитическое решение (1.1) с коэффициентами, вычисленными по формуле (1.2). Синими точками отображено конечно-элементное решение, полученное методом конечных элементов в программном комплексе Simulia Abaqus. На рис. 4, а изображено асимптотическое решение, в котором удерживается лишь одно слагаемое с коэффициентом интенсивности напряжений. На рис. 4, б красным цветом показано двучленное асимптотическое разложение. Видно, что на указанном расстоянии двучленного асимптотического разложения недостаточно для описания тензора напряжений. На рис. 5 красным цветом показано трехчленное асимптотическое разложение компоненты тензора напряжений . Видно, что на данном расстоянии трех слагаемых достаточно для описания поля напряжений. а б Рис. 4. Распределение компоненты тензора напряжений у вершины трещины на расстоянии : красным цветом показано одночленное асимптотическое разложение М. Уильямса (а), двучленное асимптотическое разложение М. Уильямса (б), синим цветом показано конечно-элементное решение Fig. 4. Distribution of the stress component at the distance from the crack tip : the red curve shows the one-term asymptotic solution (a), the two-term asymptotic solution (b), the blue points show the FEM solution Рис. 5. Распределение компоненты тензора напряжений у вершины трещины на расстоянии : красным цветом показано трехчленное асимптотическое разложение М. Уильямса, синим цветом - конечно-элементное решение Fig. 5. Distribution of stress component at the distance from the crack tip : the red curve shows the three-term asymptotic solution, the blue points show the FEM solution На рис. 6 представлено распределение компоненты на расстоянии от кончика трещины. На рис. 6, а показано двучленное разложение, на рис. 6, б - семичленное асимптотическое разложение (1.1). Из рис. 6 ясно, что на расстоянии от вершины трещины Т-напряжений недостаточно, трех- и пятичленные асимптотические разложения все еще отличаются от конечно-элементного решения и необходимо удерживать семь слагаемых в асимптотическом разложении. Из рис. 7-9 видно, что на расстоянии от кончика трещины следует удерживать семь слагаемых в многопараметрическом асимптотическом разложении М. Уильямса. На рис. 9 показано, что семичленное асимптотическое разложение полностью совпадает с конечно-элементным решением. Проведенные сравнения четко показывают, что наряду с коэффициентом интенсивности напряжений и Т-напряжениями в асимптотическом решении М. Уильямса следует удерживать высшие приближения и правильно описать напряженно-деформированное состояния в окрестности вершины трещины можно, имея коэффициенты при высших приближениях и удерживая несингулярные слагаемые. На рис. 10-13 показаны распределения компоненты тензора напряжений на расстоянии от вершины трещины, полученные в многоцелевом комплексе программ Simulia Abaqus, реализующем метод конечного элемента. Численное решение показано синими точками. На рис. 10-13 также приведено аналитическое асимптотическое решение задачи - многопараметрическое разложение М. Уильямса, в котором удерживается различное количество слагаемых: на рис. 10, а красная кривая построена с помощью одночленного асимптотического разложения. Очевидно, что на указанном расстоянии от вершины трещины однопараметрического разложения недостаточно и следует удерживать несингулярные слагаемые. На рис. 10, б красным цветом показано двучленное асимптотическое разложение М. Уильямса. Таким образом, в решении учитываются Т-напряжения. На рассматриваемом расстоянии от кончика трещины двухпараметрического разложения, учитывающего Т-напряжения, также недостаточно. а б Рис. 6. Распределение компоненты тензора напряжения на расстоянии от вершины трещины, построенное с помощью двучленного асимптотического разложения М. Уильямса (а) и семичленного асимптотического разложения (б), синим цветом изображено конечно-элементное решение Fig. 6. Distribution of stress component at the distance from the crack tip : the red curve shows the two-term asymptotic solution (a), the seven-term asymptotic solution (b), the blue points show the FEM solution а б Рис. 7. Распределение компоненты тензора напряжения на расстоянии от вершины трещины, построенное с помощью одночленного асимптотического разложения М. Уильямса (а) и двучленного асимптотического разложения (б), синим цветом изображено конечно-элементное решение Fig. 7. Distribution of stress component at the distance from the crack tip : the red curve shows the one-term asymptotic solution (a), the two-term asymptotic solution (b); the blue points show the FEM solution а б Рис. 8. Распределение компоненты тензора напряжения на расстоянии от вершины трещины, построенное с помощью трехчленного асимптотического разложения М. Уильямса (а) и пятичленного асимптотического разложения (б), синим цветом изображено конечно-элементное решение Fig. 8. Distribution of stress component at the distance from the crack tip : the red curve shows the three-term asymptotic solution (а), the five-term asymptotic solution (b); the blue points show the FEM solution Рис. 9. Распределение компоненты тензора напряжения на расстоянии от вершины трещины, построенное с помощью семичленного асимптотического разложения М. Уильямса, синим цветом изображено конечно-элементное решение Fig. 9. Distribution of stress component at the distance from the crack tip : the red curve shows the seven- term asymptotic solution, the blue points show the FEM solution На рис. 11 приведены трех- и пятичленные асимптотические разложения поля напряжений. Видно, что кривые приближаются к численному решению, но следует в асимптотическом представлении поля напряжений удерживать следующие слагаемые. Поэтому на рис. 12 построены семи- и девятичленные асимптотические разложения компоненты тензора напряжений у вершины трещины, а на рис. 13 красные кривые иллюстрируют распределения поля напряжений, полученные с помощью одиннадцатичленного и тринадцатичленного асимптотических разложений у вершины трещины. Только для последнего случая конечно-элементное решение и асимптотическое решение, содержащее 13 слагаемых, совпадают. Таким образом, можно заключить, что на указанном расстоянии для асимптотического описания поля напряжений необходимо знать 13 коэффициентов многопараметрического разложения М. Уильямса. Необходимость удержания высших приближений также иллюстрируют рис. 14-16, где приведены распределения касательного напряжения на том же расстоянии от вершины трещины. Можно было бы ожидать, что существенным будет Т-напряжение в окрестности вершины, фигурирующее в нормальном напряжении . Однако распределения касательного напряжения еще раз доказывают необходимость удержания высших приближений в асимптотическом разложении М. Уильямса. В силу этой причины, возвращаясь к экспериментальным данным, можно утверждать, что, если точка изохроматической полосы находится на расстояниях порядка полудлины трещины и более, следует удерживать высшие приближения, вплоть до 13-15 коэффициентов. Даже для расстояний, не превышающих 10 % от половины длины трещины, уже недостаточно одночленного и двучленного асимптотического разложений (см. рис. 4, 5). Поэтому при обработке всей совокупности экспериментальной информации нужно иметь в виду, что, ограничиваясь только одночленным или двучленным асимптотическими разложениями, мы не сможем правильно описать поле напряжений в окрестности вершины трещины и следует сразу прибегать к многопараметрическому асимптотическому разложению. Очевидно, что поле изохроматических полос, приведенное на рис. 2-4 первой части статьи, не может быть правильно описано с помощью одночленного или двучленного асимптотического разложения. С целью извлечения коэффициентов многопараметрического разложения из конечно-элементного эксперимента были собраны экспериментальные точки, принадлежащие концентрическим окружностям с разными радиусами: всего для вычислений параметров механики разрушения были использованы различные множества точек: от 219 до 1314. Для вычислений коэффициентов высших приближений асимптотического разложения М. Уильямса были использованы соотношения (1.1)-(1.3), которые в матричной форме могут быть представлены в виде (2.1) где - вектор-строка, состоящий из компонент тензора напряжений; - вектор-столбец, состоящий из амплитудных множителей - коэффициентов асимптотического разложения поля напряжений. Уравнения (2.1) представляют собой систему линейных алгебраических уравнений относительно масштабных множителей , которые подлежат определению (где - число взятых из вычислительного эксперимента точек). а б Рис. 10. Распределение компоненты тензора напряжения на расстоянии от вершины трещины, построенное с помощью одночленного асимптотического разложения М. Уильямса (а) и двучленного асимптотического разложения (б), синим цветом изображено конечно-элементное решение Fig. 10. Distribution of stress component at the distance from the crack tip : the red curve shows the one-term asymptotic solution (а), the two-term asymptotic solution (b) а б Рис. 11. Распределение компоненты тензора напряжения на расстоянии от вершины трещины, построенное с помощью трехчленного асимптотического разложения М. Уильямса (а) и пятичленного асимптотического разложения (б) Fig. 11. Distribution of stress component at the distance from the crack tip : the red curve shows the three-term asymptotic solution (а), the five-term asymptotic solution (b) а б Рис. 12. Распределение компоненты тензора напряжения на расстоянии от вершины трещины, построенное с помощью семичленного асимптотического разложения М. Уильямса (а) и девятичленного асимптотического разложения (б) Fig. 12. Distribution of stress component at the distance from the crack tip : the red curve shows the seven-term asymptotic solution (а) and, the nine-term asymptotic solution (b) а б Рис. 13. Распределение компоненты тензора напряжения на расстоянии от вершины трещины, построенное с помощью одиннадцатичленного асимптотического разложения М. Уильямса (а) и тринадцатичленного асимптотического разложения (б), синим цветом изображено конечно-элементное решение Fig. 13. Distribution of stress component at the distance from the crack tip : the red curve shows the eleven-term asymptotic solution (a), the thirteen-term asymptotic solution (b); the blue points show the FEM solution а б Рис. 14. Распределение компоненты тензора напряжения на расстоянии от вершины трещины, построенное с помощью одночленного асимптотического разложения М. Уильямса (а) и трехчленного асимптотического разложения (б), синим цветом изображено конечно-элементное решение Fig. 14. Distribution of stress component at the distance from the crack tip : the red curve shows the one-term asymptotic solution (a) and the three-term asymptotic solution (b) а б Рис. 15. Распределение компоненты тензора напряжения на расстоянии от вершины трещины, построенное с помощью пятичленного асимптотического разложения М. Уильямса (а) и семичленного асимптотического разложения (б), синим цветом изображено конечно-элементное решение Fig. 15. Distribution of stress component at the distance from the crack tip : the red curve shows the five-term asymptotic solution (a), the seven-term asymptotic solution (b); the blue points show the FEM solution а б Рис. 16. Распределение компоненты тензора напряжения на расстоянии от вершины трещины, построенное с помощью девятичленного асимптотического разложения М. Уильямса (а) и одиннадцатичленного асимптотического разложения (б), синим цветом изображено конечно-элементное решение Fig. 16. Distribution of the stress component at the distance from the crack tip : the red curve shows the nine-term asymptotic solution (a), the eleven-term asymptotic solution (b) В общем случае число известных из натурного или вычислительного эксперимента значений компонент тензора напряжений может быть больше числа неизвестных масштабных множителей . Очевидно, что, выбрав какое-то число экспериментальных точек, мы всегда можем выбрать следующее дополнительное множество точек, лежащих на следующей изохроматической полосе, или узлов из конечно-элементного анализа, в которых известны все компоненты тензора напряжений. Следовательно, при обработке всей совокупности экспериментальной информации или анализе конечно-элементного расчета возникает задача решения переопределенной системы линейных алгебраических уравнений (3.1). Следует отметить, что система линейных алгебраических уравнений (2.1) имеет замкнутое решение: . (2.2) Задачу отыскания амплитудных множителей можно свести к задаче минимизации функционала - целевой функции . (2.3) В ходе вычислительного эксперимента целевая функция (2.3) была использована для нахождения значений амплитудных масштабных множителей . В качестве вектора были использованы значения, вычисленные для пластины с центральной трещиной. Как было указано выше, число компонент вектора менялось от 219 до 1314, что отвечает разному количеству исследованных круговых контуров, охватывающих вершину трещины. Всего были проанализированы 60 концентрических окружностей с 73 точками на каждой. Результаты вычислений сведены в табл. 1, 2. В табл. 1 приведены данные обработки результатов конечно-элементного расчета для пластины с размерами, совпадающими с размерами пластины в натурном эксперименте (в эксперименте на фотоупругость, обсуждаемом в первой части статьи). Таблица 1 Коэффициенты асимптотического разложения М. Уильямса для пластины, идентичной натурному эксперименту с помощью метода фотоупругости [34-45], определенные с помощью метода конечных элементов Table 1 Coefficients of William’s expansion obtained by the FEM analysis В табл. 2 приведены данные расчета для пластины с малой по сравнению с шириной пластины трещиной. Для этой пластины значения масштабных множителей сравнивались с имеющимся точным аналитическим решением для бесконечной пластины. Вычисленные с помощью процедуры минимизации функционала значения коэффициентов асимптотического разложения М. Уильямса представлены в табл. 3. Во втором столбце табл. 2 приведена погрешность, с которой вычислены определенные значения (отличие от значений аналитического решения для бесконечной пластины). Таблица 2 Значения масштабных множителей - коэффициентов многопараметрического асимптотического разложения М. Уильямса для пластины с трещиной малой длины Table 2 Coefficients of multi-parameter William’s series expansion for a plate with a central crack of a small length Амплитудный коэффициент, размерность коэффициентов Погрешность, % 0,01 0,09 0,13 0,22 0,31 0,35 0,44 0,54 0,68 Таблица 3 Коэффициенты асимптотического разложения М. Уильямса у вершины трещины в полудиске с вертикальным надрезом для испытаний на трехточечный изгиб Table 3 Coefficients of the William’s series expansion in SCB specimen Конечно-элементный анализ Коэффициенты с четными номерами, как показал вычислительный эксперимент, обращаются в нуль. Таким образом, можно заключить, что подход, основанный на минимизации функционала (2.3), предоставляет удобное средство для вычисления всех коэффициентов асимптотического разложения М. Уильямса поля напряжений у вершины трещины. Как показывает табл. 2, данный подход позволяет построить пятнадцатичленное асимптотическое разложение, коэффициенты в котором определяются с погрешностью, не превышающей 0,68 %, для пятнадцатичленного асимптотического разложения, и с точностью, не превышающей 0,54 %, для тринадцатичленного асимптотического разложения. Поэтому подход для вычисления коэффициентов многопараметрического асимптотического разложения М. Уильямса, основанный на минимизации функционала (2.3), следует считать эффективным. Проведя данные вычисления на примере пластины с малой центральной трещиной, можно обратиться к компьютерному моделированию пластины с разрезом, изображенной на рис. 2-4 первой части статьи. Результаты вычислений показаны в табл. 1. 2.2 Вычисление Т-напряжений Многофункциональный комплекс Simulia Abaqus позволяет провести вычисления Т-напряжений для образцов с различным отношением полудлины трещины к ширине пластины и дать аппроксимационную оценку Т-напряжений в пластине с центральной трещиной. Для получения оценки Т-напряжений была выполнена большая серия расчетов для прямоугольных пластин с одинаковыми размерами сторон, но с различным отношением полудлины трещины к ширине пластины. Отношение полудлины трещины к половине ширины пластины менялось от 0,01 до 0,9. Для каждого значения отношения вычислялось Т-напряжение. Результаты вычислений и приближенная кривая показаны на рис. 17. Синим цветом показаны точки, полученные методом конечно-элемента в пакете Simulia Abaqus. Красным цветом показана аппроксимирующая кривая. Рис. 17. Зависимость Т-напряжений от отношения длины трещины к ширине образца Fig. 17. Variation of T-stress as function of ratio Для вычисления Т-напряжений использовалось представление (2.4) где - тарировочная функция, которая в случае пластины с центральной трещиной зависит только от отношения полудлины трещины к половине ширины пластины L, (2.5) C помощью выражений (2.4) и (2.5) можно найти Т-напряжение в пластине шириной с центральной трещиной длины , не прибегая к конечно-элементному анализу. 3. Численное определение коэффициентов асимптотического разложения М. Уильямса для полудиска с надрезом из фотоупругого эксперимента и конечно-элементного анализа Наряду с пластиной с центральной трещиной в настоящей работе рассмотренная методика была применена к сравнительно новому образцу с надрезом - полудиску с вертикальным надрезом, который используется для испытаний на трехточечный изгиб. Целью настоящего экспериментально-численного подхода является определение коэффициентов высших приближений многопараметрического разложения М. Уильямса для данной геометрии. Размеры образца следующие: радиус полудиска 5 см, длина надреза 1,92 см, толщина образца 4,9 мм. Рис. 18. Распределение интенсивности касательных напряжений в полудиске с вертикальным надрезом Fig. 18. Distribution of the von Mises equivalent stress Для проверки полученных экспериментальных результатов был выполнен расчет напряженно-деформированного состояния в полудиске с надрезом с теми же геометрическими размерами в конечно-элементном пакете Simulia Abaqus. Конечно-элементный пакет позволяет непосредственно (сразу после расчетов, не прибегая к написанию дополнительных процедур) получить значения коэффициента интенсивности напряжений и Т-напряжения. Полученные в Simulia Abaqus значения коэффициента интенсивности напряжений и Т-напряжения были сравнены с экспериментальными значениями. Различия не превышают 0,01 %. Коэффициенты высших приближений были вычислены с помощью соотношений (2.2) и с помощью минимизации функционала (2.3). Результаты вычислений сведены в табл. 3. Результаты минимизации функционала (2.3) и значения коэффициента интенсивности напряжений и Т-напряжения полностью совпадают, что свидетельствует о достоверности расчетов для коэффициентов высших приближений. Сравнение результатов, полученных с помощью методов цифровой упругости и конечных элементов, показывает их хорошее соответствие. Выводы В работе с помощью конечно-элементного анализа, проведенного в многофункциональном комплексе Simulia Abaqus, вычислены коэффициенты многопараметрического асимптотического разложения М. Уильямса поля напряжений у вершины трещины для двух различных конфигураций: для пластины с центральной трещиной в поле одноосного растяжения и полудиска с надрезом в условиях трехточечного изгиба. Из конечно-элементного расчета с помощью описанного алгоритма извлечены коэффициенты высших приближений многопараметрического разложения М. Уильямса. Алгоритм основан на рассмотрении точек, лежащих на концентрических окружностях с центром в вершине трещины. Методика позволяет в целом рассматривать любое количество концентрических окружностей и получать переопределенную систему линейных уравнений. Показано, что для определения коэффициентов разложения может быть введена целевая функция - квадратичная относительно определяемых коэффициентов, процедура ее минимизации также приводит к коэффициентам высших приближений разложения М. Уильямса. Вычислены первые пятнадцать коэффициентов асимптотического разложения поля напряжений у вершины трещины как для пластины с центральной трещиной, так и для полукругового диска с надрезом. Дано сравнение вычисленных коэффициентов разложения М. Уильямса с теоретическим решением, имеющимся для бесконечной пластины с центральной трещиной. Извлеченные из конечно-элементного решения коэффициенты совпадают с точным решением, что свидетельствует о достоверности алгоритма расчета. Получена аппроксимационная формула для Т-напряжений в пластине с трещиной, позволяющая вычислить Т-напряжений, не прибегая к экспериментальному или конечно-элементному анализу. Проведенный анализ полей напряжений у вершины трещины и острого надреза показал эффективность экспериментального метода цифровой фотоупругости, подкрепленного конечно-элементным моделированием.

About the authors

L. V Stepanova

Don State Technical University

References

Statistics

Views

Abstract - 395

PDF (Russian) - 354