О КОРРЕКТНОСТИ МЕТОДА МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ ПРОЦЕССА АКУСТИЧЕСКОГО ИЗЛУЧЕНИЯ ОТ ЛИНЕЙНОГО ДЕФЕКТА В УГЛОВОЙ ОБЛАСТИ | Опубликовать статью РИНЦ

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

Научная статья

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

ORCID 0000-0003-0915-7170,

Донской казачий государственный институт пищевых технологий и бизнеса (филиал) «МГУТУ им.К.Г. Разумовского

(ПКУ)», Ростов-на-Дону, Россия

* Корреспондирующий автор (vberkovich[at]mail.ru)

Аннотация

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

Ключевые слова: акустические свойства, излучение, граничное интегральное уравнение, неразрушающий контроль, угловая область,  пространства дробной гладкости.

ON CORRECTNESS OF MATHEMATICAL MODELING METHOD OF ACOUSTIC RADIATION PROCESS FROM LINEAR DEFECT IN ANGULAR AREA

Research article

Berkovich V.N. *

ORCID 0000-0003-0915-7170,

Don Cossack State Institute of Food Technologies and Business (branch), MSUTU named after K.G. Razumovsky

(PKU), Rostov-on-Don, Russia

* Corresponding author (vberkovich[at]mail.ru)

Abstract

This work is devoted to the analysis of vibrations in an elastic body in the form of a right dihedral angle excited by a linear defect generating acoustic radiation signals. It is assumed that during the growth process, the defect reaches the top of the angle, the boundaries of which are rigidly fixed. The author studied the problems of the correct application of the methods of the dynamic theory of elasticity for modeling this process. This situation arises when studying the problem of the reliability of the junction of elements of technological equipment operating in a dynamic mode. A similar situation arises in the geophysical analysis of wave processes in the angular blocks of the earth’s crust. The problem is reduced to the analysis of some boundary integral equation with respect to the jump in the amplitude of the stresses at the defect. The questions of solvability of the equation in spaces of fractional smoothness are studied.

Keywords: acoustic properties, radiation, boundary integral equation, non-destructive testing, angular region, spaces of fractional smoothness.

Введение

Целью настоящей  работы является  математическое моделирование одного из состояний предразрушения узла, представляющего стык угловых упругих элементов, а также изучение вопросов корректности применения традиционных методов динамической  теории упругости при математическом моделировании волновых процессов, возникающих в угловых элементах, диагностируемых методами неразрушающего контроля, как на предприятиях тяжелой индустрии, так и на предприятиях агропромышленного комплекса (предприятия хранения и переработки зерна, спиртодрожжевые и винодельческие предприятия, комбикормовые заводы и др.). В процессе длительной эксплуатации технологического оборудования в динамическом режиме появляется концентрация напряжений в угловой точке, и в её окрестности  возникает дефект, развивающийся в направлении угловой точки и генерирующий акустическое излучение (акустическая эмиссия). Аналогичная ситуация может обнаруживаться и при анализе волновых геофизических процессов при образовании излучающих дефектов в угловых стыкующихся блоках земной коры. В работе рассматривается предразрушающее состояние, когда  появление излучающего  дефекта имеет место только в одном из угловых элементов. Угловой элемент моделируется упругим телом в форме прямого двугранного угла, грани которого находятся в условиях жесткого сцепления с остальными угловыми элементами, а излучающий дефект – линейным разрезом 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°, а линейный дефект выходит не на вершину угла, также могут быть рассмотрены в рамках  предложенного подхода. Однако преобразования при этом оказываются достаточно громоздкими,  и поэтому указанные случаи в данной работе не рассматриваются.

Конфликт интересов Не указан.

Conflict of Interest None declared.