Метод сингулярных интегральных уравнений в задачах дифракции упругих волн на цилиндрических включениях а. М. Назаренко, доц., Б. Е. Панченко, канд физ мат наук, А. М. Ложкин, студ icon

Метод сингулярных интегральных уравнений в задачах дифракции упругих волн на цилиндрических включениях а. М. Назаренко, доц., Б. Е. Панченко, канд физ мат наук, А. М. Ложкин, студ




Скачати 69.39 Kb.
НазваМетод сингулярных интегральных уравнений в задачах дифракции упругих волн на цилиндрических включениях а. М. Назаренко, доц., Б. Е. Панченко, канд физ мат наук, А. М. Ложкин, студ
Дата17.07.2012
Розмір69.39 Kb.
ТипЗадача


УДК 539.3


МЕТОД СИНГУЛЯРНЫХ ИНТЕГРАЛЬНЫХ УРАВНЕНИЙ В ЗАДАЧАХ ДИФРАКЦИИ УПРУГИХ ВОЛН НА ЦИЛИНДРИЧЕСКИХ ВКЛЮЧЕНИЯХ


А.М. Назаренко, доц., Б.Е. Панченко, канд. физ.-мат. наук,
А.М. Ложкин, студ.


Сумский государственный университет


ВВЕДЕНИЕ

Проблема концентрации динамических напряжений вблизи цилиндрической поверхности разрыва в упругой среде при прохождении плоских волн была и остается объектом интенсивных исследований.

Первоначальные исследования задач дифракции упругих волн на цилиндрических включениях основывались на методах разложения в ряд по собственным функциям [1]. В последнее время значительный прогресс достигнут благодаря эффективности метода интегральных
уравнений [2, 3]. Во-первых, этот метод может быть использован применительно к неоднородностям произвольной формы. Дополнительные преимущества метода интегральных уравнений заключаются в сокращении числа пространственных переменных и возможности применения эффективных численных методов решения интегральных уравнений.

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


^ ПОСТАНОВКА ЗАДАЧИ

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

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

(1)

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

(2)

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

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


U=U0+U1, V=V0+V1. (3)


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

. (4)


^ РЕШЕНИЕ ЗАДАЧИ

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

(5)

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

(6)

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

(7)

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

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

(8)

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

, (9)

здесь f1(s) и f2(s) – неизвестные функции; Gmn – компоненты матрицы Грина (m, n=1, 2), удовлетворяющие соотношениям



(10)





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

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

.

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

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

на 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 Рисунок 4


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


SUMMARY


Singular integral equation method is applied here to solve the problem of elastic waves diffraction on a fixed cylindrical insert. Stress distribution on the insert boundary is presented.


^ СПИСОК ЛИТЕРАТУРЫ


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

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

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

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


Поступила в редакцию 17 мая 2004 г.




«Вісник СумДУ», №8(67)’ 2004

Схожі:

Метод сингулярных интегральных уравнений в задачах дифракции упругих волн на цилиндрических включениях а. М. Назаренко, доц., Б. Е. Панченко, канд физ мат наук, А. М. Ложкин, студ iconМетод сингулярных интегральных уравнений в задачах дифракции упругих волн на цилиндрических включениях а. М. Назаренко, доц., Б. Е. Панченко, канд физ мат наук, А. М. Ложкин, студ
Проблема концентрации динамических напряжений вблизи цилиндрической поверхности разрыва в упругой среде при прохождении плоских волн...
Метод сингулярных интегральных уравнений в задачах дифракции упругих волн на цилиндрических включениях а. М. Назаренко, доц., Б. Е. Панченко, канд физ мат наук, А. М. Ложкин, студ iconОсобенности решения обратной задачи в спектроэллипсометрических исследованиях в. Д. Карпуша, канд физ мат наук, доц.; У. С. Швец, асп
Целью настоящей работы является изучение особенностей применения метода спектральной эллипсометрии в исследовании оптических свойств...
Метод сингулярных интегральных уравнений в задачах дифракции упругих волн на цилиндрических включениях а. М. Назаренко, доц., Б. Е. Панченко, канд физ мат наук, А. М. Ложкин, студ iconЕлектрофізичні властивості та кристалічна структура плівок алюмінію а. Г. Басов, здобувач; А. О. Степаненко, аспірант; А. М. Чорноус, канд фіз мат наук, доц
А. Г. Басов, здобувач; А. О. Степаненко, аспірант; А. М. Чорноус, канд фіз мат наук, доц
Метод сингулярных интегральных уравнений в задачах дифракции упругих волн на цилиндрических включениях а. М. Назаренко, доц., Б. Е. Панченко, канд физ мат наук, А. М. Ложкин, студ iconЕлектрофізичні властивості та кристалічна структура плівок алюмінію а. Г. Басов, здобувач; А. О. Степаненко, аспірант; А. М. Чорноус, канд фіз мат наук, доц
А. Г. Басов, здобувач; А. О. Степаненко, аспірант; А. М. Чорноус, канд фіз мат наук, доц
Метод сингулярных интегральных уравнений в задачах дифракции упругих волн на цилиндрических включениях а. М. Назаренко, доц., Б. Е. Панченко, канд физ мат наук, А. М. Ложкин, студ iconПрограма вступного випробування з фізики для вступників на освітньо-кваліфікаційний рівень бакалавр
...
Метод сингулярных интегральных уравнений в задачах дифракции упругих волн на цилиндрических включениях а. М. Назаренко, доц., Б. Е. Панченко, канд физ мат наук, А. М. Ложкин, студ iconНаукові записки
Филипчук С. В., канд. іст наук, ст викл.; Цолін Д. В., канд філол наук, доц.; Шанюк В. І., канд філол наук, доц.; Шаправський С. А.,...
Метод сингулярных интегральных уравнений в задачах дифракции упругих волн на цилиндрических включениях а. М. Назаренко, доц., Б. Е. Панченко, канд физ мат наук, А. М. Ложкин, студ iconСумського державного педагогічного університету
Корнус А. О., канд геогр наук., доц. (відповідальний секретар); Нешатаєв Б. М., докт геогр наук, проф.; Злобін Ю. А., докт біол наук,...
Метод сингулярных интегральных уравнений в задачах дифракции упругих волн на цилиндрических включениях а. М. Назаренко, доц., Б. Е. Панченко, канд физ мат наук, А. М. Ложкин, студ iconГосударственный стандарт союза сср конструкции и изделия железобетонные радиационный метод определения толщины защитного слоя бетона, размеров и расположения
Л. Г. Родэ, канд техн наук; В. А. Клевцов, д-р техн наук; Ю. К. Матвеев; И. С. Лифанов; В. А. Воробьев, д-р техн наук; Н. В. Михайлова,...
Метод сингулярных интегральных уравнений в задачах дифракции упругих волн на цилиндрических включениях а. М. Назаренко, доц., Б. Е. Панченко, канд физ мат наук, А. М. Ложкин, студ iconО. Крамар, канд фіз мат наук; Ю. Скоренький, канд фіз мат наук; Ю. Довгоп’ятий

Метод сингулярных интегральных уравнений в задачах дифракции упругих волн на цилиндрических включениях а. М. Назаренко, доц., Б. Е. Панченко, канд физ мат наук, А. М. Ложкин, студ iconРоявление селективных процессов при стационарной неравновесной конденсации cr в. И. Перекрестов, канд физ мат наук, доц
В связи с этим в работе поставлена цель исследовать проявление селективных процессов при стационарной неравновесной конденсации Cr,...
Додайте кнопку на своєму сайті:
Документи


База даних захищена авторським правом ©zavantag.com 2000-2013
При копіюванні матеріалу обов'язкове зазначення активного посилання відкритою для індексації.
звернутися до адміністрації
Документи