Акустическое излучение в упругой среде от внутреннего дефекта с изломом
Анализ математической модели описания волнового процесса, порождаемого дефектом материала в некоторой неограниченной области, находящейся в состоянии пространственного сдвига. Разрешимость граничных интегральных уравнений и структура их решений.
Рубрика | Физика и энергетика |
Вид | статья |
Язык | русский |
Дата добавления | 12.05.2018 |
Размер файла | 573,1 K |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
Размещено на http://www.allbest.ru/
Донской казачий государственный институт пищевых технологий и бизнеса (филиал) «МГУТУ им. К.Г. Разумовского (ПКУ) » в г. Ростове-на-Дону
Акустическое излучение в упругой среде от внутреннего дефекта с изломом
Беркович В.Н.
ORCID: 0000-0003-0915-7170, Доктор физико-математических наук,
Аннотация
математический волновой дефект интегральный
Рассмотрена математическая модель описания волнового процесса, порождаемого дефектом материала в некоторой неограниченной области, находящейся в состоянии пространственного сдвига. Физический процесс возникновения колебаний изучается на стадии образования дефекта с изломом, появившимся при развитии внутреннего дефекта под воздействием нагрузок и формирующего дефект с изломом. Излучающим предполагается лишь новый дефект, появившийся в результате этого процесса. Проблема состоит в отыскании характеристик возникающей при этом акустической эмиссии (АЭ). Математическая постановка сформулированной проблемы приводит к смешанной краевой задаче математической физики. Последняя, в свою очередь, сводится к эквивалентной системе граничных интегральных уравнений (ГИУ). Установлена разрешимость ГИУ и структура их решений.
Предлагаемая к рассмотрению проблема связана с физико-математическим описанием волновых полей, порождаемых АЭ от дефектов в материалах.
Ключевые слова: неразрушающий контроль, акустическая эмиссия, дефект с изломом, граничное интегральное уравнение.
Abstract
Berkovich V.N.
ORCID: 0000-0003-0915-7170, PhD in Physics and Mathematics,
Don Cossack State Institute of Food Technology and Business branch of FSBEI of HE «K.G. Razumovsky Moscow State University of technologies and management (the First Cossack University)» in Rostov-on-Don
Acoustic radiation in flexible environment from internal defect with fracture
A mathematical model describing the wave process generated by a material defect in an unlimited region in the state of spatial displacement is considered. The physical process of the appearance of oscillations is studied at the stage of formation of a defect with a fracture that appeared when an internal defect developed under the influence of loads and formed a defect with a break. Only the new defect that appears as a result of this process is assumed to be radiating. The problem consists in finding the characteristics of acoustic emission (AE) arising in this process. The mathematical formulation of the formulated problem leads to a mixed boundary-value problem of mathematical physics. The latter, in turn, reduces to an equivalent system of boundary integral equations (BIE). The solvability of BIE and the structure of their solutions are established.
The problem proposed for consideration is related to the physical and mathematical description of the wave fields generated by AEs from defects in materials.
Keywords: non-destructive control, acoustic emission, defect with a fracture, boundary integral equation.
Задачи оценки и прогнозирования остаточного ресурса объектов технологического оборудования предприятий задачи приобретают особую актуальность, когда замена старого оборудования более современным требует серьезных материальных затрат, которые для целей производства в ряде случаев могут оказаться неоправданными.
Существующие способы анализа состояний предразрушения в процессе эксплуатации основаны на использовании различных методов неразрушающего контроля, например [1] и др. Одним из эффективных методов оценки стадий предразрушения конструкционных материалов является метод АЭ, детально исследованный и развитый в [2], [3]. Метод обладает уникальными возможностями, так как позволяет выявлять именно растущие, наиболее опасные дефекты. Для оценки параметров этих дефектов необходимо математическое моделирование волновых полей, зависящих от характера, формы и расположения дефектов в материале.
В работе [4] была построена математическая модель волнового поля при АЭ от двух независимо образовавшихся дефектов, выходящих на свободную поверхность тела. В предлагаемом исследовании моделируется фрагмент накопления повреждений, когда внутренний дефект, уже имеющийся в материале, порождает следующий, излучающий сигналы АЭ. В частности, этап ветвления трещины описан в [5]. Моделирование ситуации, изложенной выше, приводит к рассмотрению задачи динамической теории упругости об антиплоских колебаниях массивного тела, содержащего внутренний дефект J в форме ломаной , состоящей из двух звеньев длины и составляющих с направлением свободной границы Г тела углы соответственно . При этом массивное тело моделируется упругим полупространством Щ, граница Г которого предполагается свободной от нагрузок, а дефект J - разрезом в форме ломаной с изломом в точке A1. Ближайшая к свободной границе Г точка A0 удалена от нее на расстояние h. Расположение дефекта в полупространстве показано на Рис.1. Ненулевые когерентные источники гармонических колебаний смещений сдвига равной интенсивности задаются лишь на берегах разреза, ближайшего к границе в виде , на берегах смещения предполагаются равными нулю, волны смещений убегают на бесконечность, где их амплитуды затухают. При указанных условиях требуется восстановить волновое поле смещений во всей области Щ.
Рис.1 - Внутренний дефект J с изломом в упругом полупространстве Щ
Для решения проблемы рассматривается краевая смешанная задача для дифференциальных уравнений динамической теории упругости [5], [6] в упругом полупространстве Щ со свободной границей Г и разрезами J1,2 конечной длины, формирующими ломаную J. Представление результирующего колебания смещений в виде:
(1)
позволяет свести поставленную проблему к следующей задаче математической физики для амплитуд смещений полупространства Щ:
(2)
где D - плотность материала, м - модуль сдвига, v - внешняя нормаль к границе, - скачок функции на ломаном разрезе J. Во второй строке (2) указаны условия распространения волн на бесконечность, где их амплитуды исчезают Во второй строке (2) указаны условия распространения волн на бесконечность, где их амплитуды исчезают (условия излучения) [6]. Решения краевой задачи рассматривается в общем случае, как элементы пространства Соболева [7].
Метод исследования краевой задачи (2) основан на её сведении с помощью построения функций Грина для уравнения Гельмгольца (2). При этом функции Грина удовлетворяют первому граничному условию (2). Тогда имеет место представление комплексных амплитуд смещений полупространства в форме криволинейных интегралов по правым берегам разрезов :
(3)
В формуле (3) - безразмерные скачки напряжений на соответственно. Искомое поле смещений может быть восстановлено после определения скачков .
Помещая точку наблюдения на берега разрезов и переходя к локальным координатам вдоль берегов, приходим к следующей системе ГИУ относительно :
(4)
В выражениях (4) учтено, что вне разрезов, функция Макдональда [8], h - расстояние от точки A0 до свободной границы Г, -расстояние от точки A0 до точки пересечения направления звена дефекта J1 со свободной границей Г, величины связаны с расстояниями от точки излома A1 до свободной границы вдоль направлений звеньев J. Для удобства записи и описания дальнейших результатов введем вектор-функцию как продолжение вектор-функции на весь отрезок . Тогда систему (4) можно преобразовать к виду, аналогично рассмотренному в работе [4]:
(5)
При этом оказывается, что оператор левой части K представим в виде суммы однозначно обратимого оператора H и компактного оператора S, действующих в пространствах дробной гладкости.
Структура решения (5). детально исследована в [4] и имеет единственное представление:
(6)
При этом контур Г2 лежит выше -пространство вектор - функций , сходящихся к 0 с весом zу в некоторой полосе, содержащей -результат факторизации H(z) относительно .
Доказательство теоремы основано на результатах [9]. Последующая подстановка решения (6) в уравнение (5) приводит к отысканию неизвестных векторов из некоторой вспомогательной системы интегральных уравнений второго рода с компактным оператором в пространстве функций .
Компактность оператора системы позволяет представить его в виде суммы конечномерного и достаточно малого в соответствующей норме, что позволяет воспользоваться численными методами для решения системы уравнений с компактным оператором, а, следовательно, и для решения исходной краевой задачи.
Восстановление волнового поля в исследуемой области с дефектом основано на представлении (3) и решении (6) уравнений (5). Для восстановления волнового поля свободной поверхности необходимо. в формулах (3) устремить точку наблюдения (x,y) на границу Г полупространства.
Предлагаемый подход может быть применен при рассмотрении обратных задачи восстановления на основе анализа сигналов АЭ неизвестных физических параметров процесса развития внутренних дефектов внутри массивного упругого тела. Детальный анализ рассмотренной модели позволяет связать теоретические результаты с экспериментальными данными по обработки и анализу регистрируемых на свободной поверхности сигналов АЭ [10]. На основе построения необходимых амплитудно-частотных характеристик смещений точек свободной поверхности для восстановления требуемых параметров можно воспользоваться методами, разработанными в [11], [12] и др., для решения обратных задач динамической теории упругости.
Список литературы / References
1. Ермолов И.Н. Акустические методы контроля. Сер. « Неразрушающий контроль» / Ермолов И.Н., Алешин Н.П., Потапов А.И. // М. Высшая школа. Книга 2. - - 287 с.
2. Буйло С.И. Физико-механические, статистические и химические аспекты акустико-эмиссионной диагностики. - Ростов- на-Дону. Изд-во Южного федерального университета. - - 144 с.
3. Builo S.I. Physical, Mechanical and Statistical Aspects of Acoustic Emission Diagnostics. In: Physics and Mechanics of New Materials and Their Applications. New York: Nova Science Publishers. Chapter - 2013. -444 p.
4. Беркович В.Н. Колебания упругого полупространства с двумя излучающими дефектами конечной длины / Беркович В.Н., Бабкин А.В. // Международный научно-исследовательский журнал. -2017. - № 4. -С. 101-104. DOI: https://doi.org/10.23670/IRJ.2017.58.042
5. Рекач В.Г. Руководство к решению задач по теории упругости. - М. Высшая школа, 1977. - 275 с.
6. Исраилов М.Ш. Динамическая теория упругости и дифракция волн. - М. Изд. МГУ. - 1988. - 204 с.
7. Бесов О.В. Интегральные представления функций и теоремы вложения / Бесов О.В., Ильин В.П., Никольский С.М.// - М. Наука. 1975. - 478 с.
8. Лебедев Н.Н. Специальные функции и их приложения. - М.-Л.Наука. 1968. - 358 с.
9. Беркович В.Н. О точном решении одного класса интегральных уравнений смешанных задач упругости и математической физики.// Докл. АН СССР. -1982. -Т.267. -№2. - С.327- 330.
10. Буйло С.И. Диагностика предразрушающего состояния материалов по параметрам амплитудного распределения сигналов сопутствующего акустического излучения. - Дефектоскопия, 2012, № 11, С. 32-45.
11. Романов В.Г. Обратные задачи математической физики. - М. Наука. 1984. - 262с.
12. Ватульян А. О. Обратные задачи в механике деформируемого твердого тела. - М. Физматлит, - - 223 с.
Размещено на Allbest.ru
Подобные документы
Определяющие соотношения модели нелинейно упругой среды, вычисление компонент тензора напряжений. Определение автомодельного движения. Сведение модельных соотношений к системе дифференциальных уравнений. Краевая задача разгрузки нелинейно упругой среды.
курсовая работа [384,1 K], добавлен 30.01.2013Свойства и структура акустических волн. Дисперсионное соотношение для волн в неоднородной упругой среде с флуктуирующей плотностью: одномерный и трехмерный случаи. Корреляционные функции, метод релаксации для решения систем нелинейных уравнений.
контрольная работа [482,1 K], добавлен 02.01.2013Создание математической модели трехконтурной электрической схемы в среде табличного процессора Excel. Система уравнений для расчета контурных токов. Схема электрической цепи. Влияние изменения параметров схемы тяговой сети на токи тяговых подстанций.
контрольная работа [60,2 K], добавлен 14.12.2010Изучение понятия математической физики. Действительная и комплексная формы интеграла Фурье. Оригинал, изображение и операция над ними. Основные свойства преобразования Лапласа. Применение интегральных преобразований при интегрировании уравнений матфизики.
курсовая работа [281,3 K], добавлен 05.04.2014Решение уравнений, которые описывают совокупное волновое поле, создающее напряженно-деформированное состояние в окрестности кругового отверстия на безграничной тонкой упругой пластине. Основные методы применения цилиндрических функции Бесселя и Ханкеля.
курсовая работа [792,3 K], добавлен 25.11.2011Основные методы описания распространения электромагнитных волн в периодических средах с использованием волновых уравнений. Теории связанных волн, вывод уравнений. Выбор метода для описания генерации второй гармоники в периодически поляризованной среде.
дипломная работа [1,1 M], добавлен 17.03.2014Построение графика скорости центра масс фотона. Методы получения волнового уравнения Луи Де Бройля: выведение процесса описания движения центра масс фотона за рамки аксиомы. Основные математические модели, которые описывают главные характеристики фотона.
контрольная работа [628,3 K], добавлен 13.10.2010Сравнение процессов излучения и движения под действием гравитационного поля. Построение физической и математической модели окружающего нас мира. Различные положения частицы потока относительно центра потока. Увеличение длин волн линий в спектре источника.
статья [581,6 K], добавлен 15.06.2014Осциллографические методы измерения угла сдвига фаз. Измерение угла сдвига фаз методом линейной развертки. Измерение фазового сдвига путём преобразования во временной интервал. Цифровые фазометры с преобразованием фазового сдвига в постоянное напряжение.
контрольная работа [307,5 K], добавлен 20.09.2015Теоретическое описание разогрева жала паяльника с учетом потерь тепла на излучение. Средства среды MathCAD для моделирования исследуемого процесса. Решение задачи в данной среде. Составление графика зависимостей температуры, соответствующих параметрам.
контрольная работа [129,4 K], добавлен 17.12.2014