лучевая модель Траектории упругих волн в неоднородной среде
лучевая модель Траектории упругих волн в неоднородной среде
Аннотация
В настоящее время в научной литературе имеется множество теоретических и экспериментальных работ, посвященных исследованию закономерностей распространения упругих волн в различных неоднородных средах. Рассматриваются как слоисто неоднородные (или дискретно неоднородные) среды, так и непрерывно неоднородные. Слоисто неоднородные среды представляют собой многослойную структуру, в которой сами слои являются однородными, а их физические свойства резко меняются на границе между слоями. Непрерывно неоднородные среды – это среды, физические характеристики которых являются непрерывными функциями координат. Примерами непрерывно неоднородных сред являются: атмосфера Земли с непостоянной по высоте температурой; вода в океане с различной соленостью по глубине; твердотельные изделия с разной твердостью по толщине и др. Для всех подобных сред наблюдается одинаковая особенность, заключающаяся в зависимости скорости распространения упругих волн от координат. Наличие такой особенности приводит к искривлению их траектории распространения. Поэтому при исследовании закономерностей распространения упругих волн в непрерывно неоднородных средах, необходимо учитывать зависимость скорости этих волн от координат.
1. Введение
В настоящее время одной из актуальных задач, связанных с распространением упругих волн в непрерывно неоднородных средах, является проблема неразрушающего ультразвукового контроля металлоизделий с упрочненным слоем, получаемым при их термической обработке
, . В таком изделии между упрочненным и неупрочненным слоями имеется переходная область, в которой и наблюдается непрерывное изменение твердости материала изделия, а значит и скорости ультразвука, в направлении от упрочненного слоя к неупрочненному (твердость уменьшается, а скорость ультразвука увеличивается) . Наличие переходной области не позволяет определить толщину упрочненного слоя стандартными методами ультразвукового контроля, т.к. отсутствует контрастная граница, от которой наблюдалось бы отражение ультразвука. Один из выходов в этом случае заключается в использовании особенности распространения упругих волн в переходной области. Эта особенность заключается в искривлении траектории ультразвука при направлении его распространения отличном от направления градиента скорости. Поэтому необходимо теоретически и экспериментально исследовать особенность распространения упругих (в частности, ультразвуковых) волн в среде, в которой их скорость зависит от координаты. А также необходимо исследовать влияние различных физических параметров данной среды на траекторию упругих волн.В известной научной литературе практически отсутствуют сведения о математической модели, описывающей траекторию упругих волн в непрерывно неоднородных средах. Поэтому целью настоящей работы является разработка этой модели и ее теоретическое исследование.
2. Разработка модели
![Структурно-геометрическая схема рассматриваемой среды](/media/images/2024-05-06/dfff9339-b512-4997-9f3c-2c8cca8d9265.png)
Рисунок 1 - Структурно-геометрическая схема рассматриваемой среды
Для поиска уравнения, описывающего подобную траекторию, отвлечемся пока от слоистой структуры и рассмотрим по подробнее распространение упругих волн в переходной области. При этом будем пока полагать, что она является полубесконечной в направлении оси Oz и начало системы координат совпадает с ее границей, лежащей в плоскости xOy.
Для простоты рассуждений будем считать, что градиент скорости является постоянной величиной. Пусть при этом скорость будет зависеть только от координаты z. Тогда зависимость скорости от координаты можно записать в виде
![Траектория акустического луча, прошедшего насквозь переходную область](/media/images/2024-05-06/0cbc0daf-6ecc-48fe-8e48-d8a9cd6e61a4.png)
Рисунок 2 - Траектория акустического луча, прошедшего насквозь переходную область
где α – угол между направлением скорости упругой волны и осью Ox в произвольной точке траектории, а v – модуль скорости упругой волны в этой же точке. Также при выводе уравнения траектории воспользуемся свойством производной от функции в некоторой произвольной точке, а именно
Таким образом, из системы уравнений (1)-(3) после некоторых математических преобразований получаем выражение
из которого после разделения переменных, интегрирования и математических преобразований получаем искомое уравнение траектории z(x) упругой волны в переходной области
![Схематичный график уравнения (4): траектория упругой волны в переходной области](/media/images/2024-05-06/1743802f-b467-4def-ae35-0b0479ba5053.png)
Рисунок 3 - Схематичный график уравнения (4): траектория упругой волны в переходной области
Следует заметить, что формулы (5) и (6) справедливы только для переходного слоя. В упрочненном и неупрочненном слое зависимости z(x) являются линейными, и их можно представить в виде
где – координата по оси Oz границы между переходной областью и неупрочненным слоем; xD – координата по оси Ox точки пересечения акустическим лучом границы между переходной областью и неупрочненным слоем, которую можно определить, приняв в уравнении (6) z=D.
В результате из выше изложенного следует, что траектория акустического луча, изображенная на рис. 2, будет описываться кусочной функцией в виде
![Траектория акустического луча, развернувшегося в переходной области](/media/images/2024-05-06/731104a5-8410-40b7-8625-aa9071aa0fe6.png)
Рисунок 4 - Траектория акустического луча, развернувшегося в переходной области
Здесь x2 – координата выхода акустического луча из переходной области обратно в упрочненный слой; b – координата по оси Ox выхода акустического луча из рассматриваемой трехслойной структуры.
Исследуем функцию (8) на экстремум и найдем координаты вершины графика этой функции:
С другой стороны, исходя из симметрии траектории, координату вершины можно найти как точку пересечения продолжения прямолинейных участков траектории из формулы (8):
Тогда, подставив последнее в формулу (9), получаем выражение для вычисления координаты выхода акустического луча из трехслойной структуры:
Уравнения (9) и (10), как и ожидалось, показывают, что координаты вершины траектории зависят от толщины упрочненного слоя, градиента скорости и угла α0, под которым акустический луч входит в рассматриваемую среду. Но так как толщина переходной области ограничена, то вершина траектории акустического луча может быть, либо действительной при zв≤D, либо мнимой при zв>D. Действительная вершина располагается в переходной области и лежит на траектории луча, а мнимая вершина располагается в неупрочненном слое и лежит на продолжении криволинейного участка траектории луча в неупрочненный слой.
Другими словами расчеты по формулам (9) и (10) дают не только информацию о координатах вершины траектории акустического луча, но и позволяют понять каким будет характер этой траектории – возвратный или сквозной. Очевидно, что при zв>D траектория будет сквозной, а при zв≤D – возвратной.
3. Результаты модельных расчетов траектории ультразвукового луча и их анализ
Для проведения модельных расчетов и построения траектории ультразвукового луча возьмем границы упрочненного и переходного слоев и значения скорости ультразвука характерные для стали: d=0,04 м, D=0,14 м, v0=5920 м/с, v=5980 м/с. Здесь v0 – скорость ультразвука в упрочненном слое и на его границе с переходной областью; v – скорость ультразвука в неупрочненном слое и на его границе с переходной областью. Тогда для градиента скорости ультразвука в переходной области можно записать выражение . Все расчеты будем проводить на основании уравнений, полученных в предыдущем разделе.
![Траектории ультразвука при значениях угла α0 = 1º–10º с шагом 1º](/media/images/2024-05-06/017b7cc6-fa3b-4907-aa06-21ca962de482.jpg)
Рисунок 5 - Траектории ультразвука при значениях угла α0 = 1º–10º с шагом 1º
Кроме этого анализ также показывает, что углы α0 и β будут равны друг другу, если толщина упрочненного слоя будет соответствовать двум третьим толщины всей двухслойной структуры, т.е.
![Траектории ультразвука с одной точкой выхода](/media/images/2024-05-06/cd971a91-6b8e-4b0d-bc3e-3730c7013e0c.jpg)
Рисунок 6 - Траектории ультразвука с одной точкой выхода
Примечание: горизонтальным линиям соответствуют границы упрочненного и переходного слоев
Далее проведем расчет зависимости траектории ультразвукового луча от толщины упрочненного слоя и толщины переходной области. Для этих вычислений в качестве исходных данных будем использовать те же значения скоростей. Значение угла α0 возьмем 8º, что соответствует траектории из рис. 5, у которой вершина почти касается границы z=D.
На рис. 8 представлены результаты расчета траектории ультразвука в зависимости от различных значений толщины переходного слоя D. Здесь также точками указана граница упрочненного слоя, а пунктирными линиями указаны задаваемые положения границы переходного слоя. Из рисунка видно, что значение D очень сильно влияет на траекторию ультразвука и, как следствие, сильно влияет на положение точки обратного выхода ультразвука из рассматриваемой среды. При этом чем меньше толщина переходного слоя, тем меньше расстояние между точками входа и выхода траектории ультразвука в рассматриваемой среде. Из результатов расчета следует, что изменение толщины переходного слоя на 3 см приводит к смещению точки выхода траектории на величину порядка 80–90 cм.
![Траектории ультразвука при d = 2D/3 и значениях угла α0=1º–10º с шагом 1º](/media/images/2024-05-06/f7653479-dda8-446c-87c6-9b1cff93790d.jpg)
Рисунок 7 - Траектории ультразвука при d = 2D/3 и значениях угла α0=1º–10º с шагом 1º
![Траектории ультразвука при значениях толщины переходного слоя: D = 0,05–0,14 м, шаг равен 0,03 м](/media/images/2024-05-06/55019949-dc4f-4d15-aa13-802533ff58c9.jpg)
Рисунок 8 - Траектории ультразвука при значениях толщины переходного слоя: D = 0,05–0,14 м, шаг равен 0,03 м
![Траектории ультразвука при значениях толщины упрочненного слоя: d = 0,02–0,08 м, шаг равен 0,02 м](/media/images/2024-05-06/b8c51380-eb09-4b8f-9b03-3acd23247ad1.jpg)
Рисунок 9 - Траектории ультразвука при значениях толщины упрочненного слоя: d = 0,02–0,08 м, шаг равен 0,02 м
4. Заключение
Из всех возможных траекторий ультразвука в рассматриваемой структуре наибольший интерес представляет такая, у которой вершина касается границы переходной области с неупрочненным слоем, т.е. координата вершины которой удовлетворяет условию zв=D. При проведении анализа этой траектории выяснилось, что независимо от толщины упрочненного слоя и переходной области вершина траектории всегда будет касаться границы z=D, если ввести луч в рассматриваемую структуру под углом, удовлетворяющим уравнению:
где . Т.е. этот угол зависит только от скоростей ультразвука на границах переходного слоя z=d и z=D.