О корректности метода математического моделирования процесса акустического излучения от линейного дефекта в угловой области

Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

О корректности метода математического моделирования процесса акустического излучения от линейного дефекта в угловой области

Беркович В.Н.

Работа посвящена анализу колебаний в упругом теле в форме прямого двугранного угла, возбуждаемых линейным дефектом, генерирующим сигналы акустической излучения. Предполагается, что в процессе роста дефект выходит на вершину угла, границы которого жестко закреплены. Изучаются вопросы корректности применения методов динамической теории упругости для моделирования указанного процесса. Указанная ситуация возникает при изучении проблемы надежности стыка элементов технологического оборудования, работающего в динамическом режиме. Аналогичная ситуация возникает при геофизическом анализе волновых процессов в угловых блоках земной коры. Проблема сведена к анализу некоторого граничного интегрального уравнения относительно скачка амплитуды напряжений на дефекте. Изучены вопросы разрешимости уравнения в пространствах дробной гладкости.

Целью настоящей работы является математическое моделирование одного из состояний предразрушения узла, представляющего стык угловых упругих элементов, а также изучение вопросов корректности применения традиционных методов динамической теории упругости при математическом моделировании волновых процессов, возникающих в угловых элементах, диагностируемых методами неразрушающего контроля, как на предприятиях тяжелой индустрии, так и на предприятиях агропромышленного комплекса (предприятия хранения и переработки зерна, спиртодрожжевые и винодельческие предприятия, комбикормовые заводы и др.). В процессе длительной эксплуатации технологического оборудования в динамическом режиме появляется концентрация напряжений в угловой точке, и в её окрестности возникает дефект, развивающийся в направлении угловой точки и генерирующий акустическое излучение (акустическая эмиссия). Аналогичная ситуация может обнаруживаться и при анализе волновых геофизических процессов при образовании излучающих дефектов в угловых стыкующихся блоках земной коры. В работе рассматривается предразрушающее состояние, когда появление излучающего дефекта имеет место только в одном из угловых элементов. Угловой элемент моделируется упругим телом в форме прямого двугранного угла, грани которого находятся в условиях жесткого сцепления с остальными угловыми элементами, а излучающий дефект - линейным разрезом I конечной длины l, на берегах которого расположены источники колебаний (Рис. 1). Разрез I выходит на вершину углового элемента под углом в к его границе.

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

Учитывая сказанное выше, будем рассматривать случай, когда границы упругой угловой области Щ являются неподвижными, а на разрезе I находятся источники гармонических колебаний пространственного сдвига, моделирующие сигналы АИ, задаваемые с помощью функции , где щ - круговая частота колебаний. При этом предполагается, что источники колебаний на берегах разреза когерентны и имеют одинаковую интенсивность.

Рис. 1. Линейный излучающий разрез I в угловой области, смещения границ разреза перпендикулярны плоскости чертежа

Известно, что в этом случае колебания смещений U упругой среды во времени ф могут быть описаны с помощью волнового уравнения (a-скорость распространения волн сдвига):

(1)

Считая колебания установившимися и отыскивая смещения U в виде для изучения акустических свойств среды, приходим к следующей краевой задаче для уравнения Гельмгольца [1] относительно комплексной амплитуды смещений в области прямого двугранного угла Щ с выходящим на его ребро линейным разрезом I конечной длины l:

(2)

где м - модуль сдвига, D - плотность материала угла, Г - граница угловой области Щ, оси координат x, y направлены вдоль прямолинейных границ Г. Будем предполагать, что при функция и выполняются условия излучения Зоммерфельда [2].

Ставится задача о восстановлении волнового поля во всей рассматриваемой области Щ c помощью амплитуды смещений берегов разреза .

Замечание. Отметим, что, несмотря на локальный характер процесса образования дефекта, генерирующего сигналы акустического излучения (АИ), методика анализа акустических свойств и моделирования процесса антиплоскими колебаниями сдвига, тем не менее, приводит к удовлетворительным качественным результатам, как отмечается в работе [3] при моделировании явления акустической эмиссии, порождаемой линейными излучающими дефектами в предразрушающих состояниях упругих материалов.

Описание метода и полученные результаты

1. Метод исследования краевой задачи (2) основан на ее сведении к некоторому граничному интегральному уравнению (ГИУ) используя функцию Грина, с помощью которой любое регулярное решение уравнения Гельмгольца в конечной области S представимо в виде:

где C - граница области S, v - нормаль к границе, - функция Грина для прямого угла без разреза, удовлетворяющая следующей краевой задаче ( - дельта - функция Дирака):

(4)

(5)

Кроме того, при выполняется условие

(6)

и условия излучения Зоммерфельда.

Для нахождения функции Грина использовался метод интегрального преобразования Фурье [4]:

(7)

и свойства д - функции Дирака

(8)

Выбирая функцию Грина, а также д - функцию Дирака в форме интегрального преобразования Фурье (7), (8) и, осуществляя преобразования образующихся при этом контурных интегралов, приходим к следующему выражению вспомогательной функции G₀ для полупространства, удовлетворяющей только одному граничному условию

(9)

При этом полученные контурные интегралы были преобразованы с использованием интегрального представления функции Ханкеля [2]:

Метод отражений позволяет с помощью вспомогательной функции G₀ (9) для полупространства построить функцию Грина G для прямого угла, которая будет удовлетворять граничным условиям (4) - (6) и имеет вид:

(10)

2. Для получения граничного интегрального уравнения (ГИУ)

воспользуемся интегральным представлением (3). Выберем в (3) в качестве области S часть угловой области Щ, ограниченной контуром , где - окружность радиуса R замыкающая прямой угол с центром в его вершине. Перейдем в соотношении (3) к пределу при с учетом условий убывания (6), условий Зоммерфельда и асимптотических представлениями функций Ханкеля [4]:

В результате указанных преобразований в соотношении (3) остается лишь интеграл по разрезу I, и выражение смещений в точках области Щ принимает вид:

(11)

Устремим точку на границу разреза I в соотношении (11). Тогда, переходя к локальным координатам по формулам:

приходим к следующему ГИУ относительно скачка нормальной производной (безразмерных амплитуд напряжений) на разрезе I по заданным на нём амплитудам смещений

(12)

где функция G определяется равенством (10).

Перейдем в представлении (10) функции G от функций Ханкеля к функциям Макдональда [4]:

Применение формул сложения [4] для функций Макдональда преобразует ГИУ (12) к следующему виду:

(13)

3. Результат исследования ГИУ (13) может быть сформулирован в форме следующей теоремы.

Теорема. Оператор K левой части ГИУ (12) однозначно обратим, как оператор, действующий в пространствах Соболева-Слободецкого [6] по правилу:

Следуя известному в теории дифракции методу Фока, будем временно считать и исследуем ГИУ (13), рассмотренное выше. Тогда оператор K левой части ГИУ оказывается положительно определенным в пространстве обобщенных решений уравнения (13), которое вводится с помощью нормы

со скалярным произведением в виде (* означает комплексное сопряжение):

Далее используется классическая теорема Рисса о единственности представления линейного непрерывного функционала в гильбертовом пространстве [5], которое приводит к следующему условию разрешимости ГИУ (13):

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

(14)

Левая и правая части (14) представляют норму в пространстве дробной гладкости Соболева-Слободецкого , которому должна принадлежать функция правой части (12). Из теоремы Рисса и оценок (14) следует , где пространство является сопряженным [6] к пространству . Доказательство обратного вложения приводит к результату теоремы.

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

Описанный результат согласуется с известными представлениями о граничных свойствах функций из пространств Соболева, в которых отыскиваются решения краевых задач динамической теории упругости. Для построения решения (13) можно использовать результаты работы [8] и восстановить волновое поле смещений в Щ с помощью представления (11), которое будет выражено через амплитуды излучения с границ дефекта I.

Это позволяет утверждать однозначную разрешимость исходной краевой задачи (2) в пространстве Соболева .

Результаты, полученные для рассмотренного примера, позволяют утверждать корректность применения метода математического моделирования волнового процесса АИ на основе использования классических моделей динамической теории упругости. Восстановленное волновое поле в Щ при решении прямой задачи может служить основой для постановки и решения обратных задач [3], [9], [10] восстановления параметров излучения дефекта I по характеристикам сигналов АИ, регистрируемых акустико-эмиссионной либо геофизической аппаратурой.

Предложенный выше подход позволяет рассматривать и случаи других граничных условий на полубесконечных границах угловой области. Случаи, когда углы раствора угловых областей не равны 90°, а линейный дефект выходит не на вершину угла, также могут быть рассмотрены в рамках предложенного подхода. Однако преобразования при этом оказываются достаточно громоздкими, и поэтому указанные случаи в данной работе не рассматриваются.

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

колебание акустический излучение упругость

1. Рекач В.Г. Руководство к решению задач по теории упругости / В.Г. Рекач. - М.: Высшая школа. -1977. -275 с.

2. Купрадзе В.Д. Граничные задачи теории колебаний и интегральные уравнения/ В.Д. Купрадзе. - М.: ГИТТЛ. - 1958. -280 с.

3. Berkovich V.N. Amplitude Reconstruction of the Defect Radiation by Acoustic Emission Signals on the Unloaded Boundary of a Massive Body / V.N. Berkovich, S.I. Builo. // Russian Journal of Nondestuctive Testing. - 2019. - №4. - pp. 262-267.

4. Лебедев Н.Н. Специальные функции и их приложения / Н.Н. Лебедев. - М.-Л.: Наука. - 1968. -358 с.

5. Люстерник Л.А. Элементы функционального анализа / Л.А. Люстерник, В.И. Соболев. - М.: Наука. 1965. -519 с.

6. Бесов О.В. Интегральные представления функций и теоремы вложения / О.В. Бесов, В.П. Ильин, С.М. Никольский. - М.: Наука. - 1975. -478 с.

7. Привалов И.И. Введение в теорию функций комплексного переменного. М.: Наука. - 1967. -376 с.

8. Беркович В.Н. О точном решении одного класса интегральных уравнений смешанных задач упругости и математической физики/ В.Н. Беркович // Докл.АН СССР. -1982. - Т.267. - №2. - С. 327-330.

9. Ватульян А.О. Прямые и обратные задачи для однородных и неоднородных упругих и электроупругих тел/ А.О. Ватульян, А.Н. Соловьев. - Ростов-на-Дону: Изд-во ЮФУ. - 2008. - 175 с.

10. Романов В.Г. Некоторые обратные задачи для уравнений гиперболического типа/ В.Г. Романов. - Новосибирск: Наука. -1972. -315 с.

11. Рябой В.З. Кинематические и динамические характеристики глубинных волн, связанных с границами в земной коре и верхней мантии/ Рябой В.З. // Изв. АН СССР. - Сер.геофиз. - 1966 - №3. - С. 175-182

Размещено на Allbest.ru