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




НазваниеМетод сингулярных интегральных уравнений в задачах дифракции упругих волн на цилиндрических включениях
Дата11.10.2012
Размер69.2 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Программа Государственного экзамена по подготовке магистра по направлению «Физика акустических и гидродинамических волновых процессов» (510414)
Влияние границ и неоднородности среды на распространение волн и структуру волнового поля. Основные методы анализа волн (спектральный...
Метод сингулярных интегральных уравнений в задачах дифракции упругих волн на цилиндрических включениях iconСтатья 62 Некоторые интегральные тождества математической физики
Грина построены три группы интегральных соотношений, верных для функций, заданных на границе области. Тождества пригодны для исследования...
Метод сингулярных интегральных уравнений в задачах дифракции упругих волн на цилиндрических включениях iconВажным видом математической идеализации звуковых волн являются цилиндрические волны. Поверхность одинаковых фаз
Поверхность одинаковых фаз1 является в этом случае цилиндрической. Очевидно, что для описания данного типа волн удобно использовать...
Метод сингулярных интегральных уравнений в задачах дифракции упругих волн на цилиндрических включениях iconЛекция Содержание
Законы или правила Кирхгофа. Делители напряжений и токов. Возможные методы упрощения систем уравнений (метод узловых потенциалов...
Метод сингулярных интегральных уравнений в задачах дифракции упругих волн на цилиндрических включениях iconПрограмма вступительных испытаний в магистратуру по направлению 010100. 68 Математика Программа обсуждена на заседании кафедры ит
Системы линейных уравнений. Метод Гаусса решения систем линейных уравнений. Определитель матрицы. Свойства определителя. Метод Крамера...
Метод сингулярных интегральных уравнений в задачах дифракции упругих волн на цилиндрических включениях iconЗакономерности распространения волн в упругих средах
Задача дисциплины – получение теоретических знаний и практических навыков необходимых при изучении курса "Промышленная акустика"
Метод сингулярных интегральных уравнений в задачах дифракции упругих волн на цилиндрических включениях iconГоувпо «Самарский государственный архитектурно-строительный университет» Факультет информационных систем и технологий
Исследование, нлу, методы решения нелинейных уравнений, метод Ньютона, метод касательных, метод хорд, метод половинного деления,...
Разместите кнопку на своём сайте:
Библиотека


База данных защищена авторским правом ©lib.znate.ru 2014
обратиться к администрации
Библиотека
Главная страница