Метод сингулярных интегральных уравнений в задачах дифракции упругих волн на цилиндрических включениях

Постановка задачи

Рассмотрим в неограниченной изотропной среде бесконечный вдоль оси OZ неподвижный цилиндр, поперечное сечение которого ограничено замкнутым контуром L типа Ляпунова. Пусть также внешнее поле перемещений действует перпендикулярно оси OZ.

При таких предположениях мы находимся в условиях плоской деформации. В качестве внешнего воздействия будем рассматривать набегающую на цилиндр из бесконечности монохроматическую волну расширения-сжатия (Р-случай)

(1)

или волну сдвига (SV-случай)

(2)

здесь ? - амплитуда падающей волны; с₁ и с₂ - скорости продольной и поперечной волн; ? - частота колебаний; t - время; ? и ? - постоянные Лямэ; ? - плотность среды; i - мнимая единица ().

При взаимодействии приходящей волны с цилиндром возникают отраженные волны двух типов (продольные и поперечные), причем другие типы волн не образуются. Пусть - смещения отраженного поля. Тогда общее поле амплитуд перемещений равно

U=U₀+U₁, V=V₀+V₁. (3)

Будем предполагать, что поперечное сечение цилиндра описывается гладкой замкнутой кривой L, в точках которой будем удовлетворять граничные условия. В случае неподвижного цилиндра они имеют вид

. (4)

Решение задачи

В случае установившихся волновых движений упругого тела (зависимость от времени выражается множителем ) амплитудные значения отраженных волн перемещений удовлетворяют соотношениям

(5)

Амплитудные значения напряжений связаны с амплитудами перемещений U и V формулами ():

(6)

Пусть L - некоторая кривая в поперечном сечении цилиндра. Обозначим через S₁ и S₂ амплитуды тангенциальной и нормальной компонент вектора напряжений на L. Тогда в произвольной точке кривой эти напряжения выражаются через компоненты тензора амплитуд напряжений следующим образом:

 (7)

где - угол положительной касательной к L в точке с осью Ох.

На границе тела представляют интерес распределения компонент тензора амплитуд напряжений , которые будем находить по формулам

(8)

Будем строить интегральные представления амплитуд перемещений U₁ и V₁, чтобы они автоматически удовлетворяли уравнениям движения (5) и условиям излучения на бесконечности, т. е. чтобы они представляли собой расходящиеся волны. Следуя [3], представим U₁ и V₁ в виде потенциалов типа простого слоя:

, (9)

здесь f₁(s) и f₂(s) - неизвестные функции; G_mn - компоненты матрицы Грина (m, n=1, 2), удовлетворяющие соотношениям

(10)

где - функция Ханкеля 1-го рода j-го порядка.

Анализ формул (10) показывает, что функции и непрерывны в нуле, а функции и обладают логарифмической особенностью

По этой причине подстановка представлений (9) в граничные условия (4) сводит краевую задачу к системе двух сингулярных уравнений с логарифмическими ядрами, численная реализация которых затруднительна.

С целью получения сингулярных интегральных уравнений с ядром типа Коши [4] представления (9) дифференцировались по дуговой координате s₀. Имеем

на L,(11)

,

Вычисление необходимых для (11) производных дает

(12)

Можно показать, что ядро , определенное в (10), является непрерывным, а ядро - сингулярно (, где - непрерывно).

Подставляя (12) в граничные условия (11), приходим к следующей системе сингулярных интегральных уравнений

(13)

, , ,

в P-случае,

в SV-случае.

Необходимые дополнительные условия для разрешимости сингулярных интегральных уравнений 1-го рода (13) вытекают, например, из равенства нулю смещений на L в некоторой фиксированной точке или из равенства нулю средних смещений на L. В последнем случае имеем (l - длина контура L)

(14)

в Р-случае,

в SV-случае.

Численная реализация алгоритма

Для численной реализации алгоритма в настоящей работе использован метод, теоретически обоснованный в работе [4] и основанный на приближении плотностей интегральных уравнений тригонометрическими многочленами и в последующем точном вычислении интегралов с непрерывными и сингулярными ядрами.

Проведем параметризацию контура L по формулам

Интерполяционный многочлен для неизвестных плотностей интегральных уравнений (13) имеет вид

(15)

где в случае четного числа узлов N=2n и в случае нечетного N=2n+1.

Подстановка (15) в интегралы с сингулярными ядрами дает

(16)

если - непрерывное ядро и

(17)

в случае ядра Гильберта, причем квадратурные формулы (16), (17) имеют место как при четном, так и нечетном числе узлов разбиения контура L.

Отметим, что формула (17) аналогична правилу приближенного вычисления регулярных интегралов (16). По этой причине при численной реализации сингулярных интегралов ядро Гильберта выделять из сингулярного ядра необязательно. В работе как для регулярных, так и для сингулярных интегралов использовалась квадратурная формула (16).

В качестве примера рассматривалось пространство, содержащее цилиндрическое неподвижное включение эллиптического поперечного сечения

(18)

На контуре включения проводилось вычисление напряжений

где компоненты амплитуд напряжений , находились по формулам (8), Р - максимальное значение напряжения в падающей волне, равное в случае излучения Р-волны (1) и - в случае излучения SV-волны (2).

Отметим, что в случае неподвижного включения напряжение всегда меньше и связано с последним соотношением

На рис. 1, 2 и 3, 4 показано изменение напряжений , на контуре эллиптического неподвижного включения для случаев излучения из бесконечности продольной волны и волны сдвига соответственно.

На рис. 1, 2 кривые 1, 2, 3, 4 отвечают случаям =0,4; 0,7; 1,0; 1,3 соответственно при b/a=2, =0,3. Видно, что в довольно большой окрестности точки соскальзывания () значение практически не влияет на величину напряжения . Напряжение вблизи точки принимает максимальное значение. Преобладающим здесь является напряжение .Кривые 1, 2 и 3 отвечают значениям параметра =0,2; 0,3 и 0,4 соответственно (b/a=2, ). Здесь в окрестности теневой () и лобовой () точек преобладающим является напряжение . В освещенной зоне с увеличением параметра , значения напряжений также увеличиваются. Напряжение принимает свое максимальное значение вблизи точки , а напряжение - в лобовой точке ().

Список литературы

1. Гузь А.Н., Кубенко В.Д., Черевко М.А. Дифракция упругих волн. - Киев: Наукова думка, 1978. - 307с.

2. Фильштинский Л.А. Дифракция упругих волн на трещинах, отверстиях, включениях в изотропной среде //Изв. АН СССР. Мех. твердого тела. - 1991. - №4. - С.19-127.

3. Назаренко А.М., Панченко Б.Е. Дифракция волн сдвига на цилиндрических неоднородностях произвольного поперечного сечения // Динам. и прочность машин. Респ. межвед. научно-техн. сб., 1991.- Вып. 52. - С. 38-45.

4. Панасюк В.В., Саврук М.П., Назарчук З.Т. Метод сингулярных интегральных уравнений в двумерных задачах дифракции. - Киев: Наукова думка, 1984. - 344 с.