Восстановление формы дефекта ферромагнитной пластины путем решения обратной задачи магнитостатики и серии прямых задач
- Авторы: Никитин А.В.1, Михайлов Л.В.1, Михайлов А.В.1, Гобов Ю.Л.1, Костин В.Н.1, Смородинский Я.Г.1
-
Учреждения:
- Институт физики металлов им. М.Н. Михеева УрО РАН
- Выпуск: № 9 (2024)
- Страницы: 67-72
- Раздел: По материалам XXXV Уральской конференции «Физические методы неразрушающего контроля (Янусовские чтения)»
- URL: https://journals.rcsi.science/0130-3082/article/view/265161
- DOI: https://doi.org/10.31857/S0130308224090086
- ID: 265161
Цитировать
Полный текст
Аннотация
Представлена методика верификации решения обратной геометрической задачи магнитостатики в пластине из магнитомягкого ферромагнетика. Методика состоит из решения ряда прямых задач, в которых в качестве первого приближения используют форму дефекта, полученного в результате решения обратной геометрической задачи магнитостатики, а далее, увеличивая либо уменьшая глубину дефекта, не меняя при этом форму граничной поверхности, — сравнения топографий компонент магнитного поля, полученных в ходе измерений над поверхностью пластины и вычисленных (в результате решения прямой задачи) в тех же точках компонент магнитного поля рассеяния от реконструированного трехмерного дефекта. В результате применения методики также могут быть уточнены геометрические параметры исследуемого дефекта. Получение начальных условий для решения обратной задачи и решение прямых задач магнитостатики осуществляется с помощью метода конечных элементов в программе ELMER. Методика работает при одностороннем доступе к любой поверхности пластины (бездефектной поверхности, либо поверхности с дефектом).
Полный текст
На начальном этапе развития неразрушающего контроля методом регистрации магнитного потока рассеяния (MFL-метод) предполагалось, что при помощи различных подходов, путем сравнения топографии измеренного магнитного поля над поверхностью объекта контроля, содержащего дефект, с топографией поля эталонного дефекта удастся определить тип и геометрические параметры формы обнаруженного дефекта [1—6]. Однако большое разнообразие дефектов, а также нелинейный отклик ферромагнетика на внешнее намагничивающее поле, не позволяют устойчивым образом решать такие задачи. Области применения полученных аналитических формул для вычисления компонент магнитного поля эталонных дефектов в зависимости от конкретных параметров задачи ограничены и требуют большого числа поправочных коэффициентов [7], а численные расчеты прямых задач магнитостатики позволяют получать топографии магнитных полей только от конкретной геометрии границы металла. По этим причинам для решения практических задач по определению геометрических параметров формы исследуемого дефекта становится важной разработка методик решения обратной геометрической задачи магнитостатики для трехмерного случая с учетом нелинейного характера зависимости магнитной индукции в ферромагнетике от внешнего намагничивающего поля. Решение обратной геометрической задачи магнитостатики не требует иной априорной информации кроме информации о характере нелинейной зависимости величины магнитной индукции ферромагнитного материала от величины внешнего намагничивающего поля. Впервые такая обратная геометрическая задача магнитостатики для двумерного линейного случая была решена в работе [8]. В работе [9] были получены расчетные формулы для определения граничных поверхностей дефекта потери сплошности металла для трехмерного линейного случая и при условии, что ферромагнетик намагничен до насыщения, а именно, вектор намагниченности во всем ферромагнетике постоянен по величине и направлению.
В работах [10—12] была разработана методика решения обратной геометрической задачи магнитостатики в трехмерном случае с учетом нелинейного отклика ферромагнетика на внешнее намагничивающее поле для класса дефектов, чьи граничные поверхности можно описать с заданной точностью отрезком двойного ряда Фурье (рис. 1). Из этого класса дефектов исключены трещины и трещиноподобные дефекты (определяем такие дефекты согласно [13]). В дальнейшем будем использовать результаты именно этой методики решения обратной геометрической задачи магнитостатики.
Рис. 1. Пример функции f(x, y), представимой в виде отрезка двойного ряда Фурье в исследуемой области.
Целью данной работы является реализация алгоритма методики по проверке правильности решения обратной геометрической задачи магнитостатики и, возможно, последующего уточнения геометричских параметров формы дефекта на численном эксперименте.
МЕТОДИКА ИССЛЕДОВАНИЙ
Представленная в работах [10—12] методика вне зависимости от взаимного расположения дефекта и магнитных преобразователей (рис. 2) позволяет решить обратную геометрическую задачу магнитостатики и найти приближенную форму дефекта без использования каких-либо поправочных коэффициентов.
Рис. 2. Сечение исследуемой ферромагнитной пластины в поперечном направлении. Намагничивающее поле H0 направлено вдоль оси ОХ.
Для определенности выберем конфигурацию 1 на рис. 2 и на ее примере покажем, как работает предлагаемая методика. В результате решения обратной геометрической задачи получим приближенную форму дефекта, как на рис. 3—4.
Рис. 3. Поперечное сечение пластины с реконструкцией границы металла в области дефекта, полученное в результате решения обратной геометрической задачи.
Рис. 4. Граница куполообразного дефекта с размерами 25×25×3 мм, полученная в результате решения обратной геометрической задачи магнитостатики для пластины из стали 20 толщиной 10 мм.
Для получения начальных условий, т.е. «измеренных» значений компонент магнитного поля над бездефектной поверхностью, использовалось решение прямой задачи в программе ELMER.
Заметим, что размеры дефекта на поверхности пластины определяются по топографии измеренных компонент магнитного поля рассеяния. Поэтому реконструируем дефект с учетом полученной в результате измерений магнитного поля информации о краях дефекта на поверхности пластины (рис. 5). Затем сдвинем поверхность дефекта, целиком углубляя и уменьшая глубину дефекта.
Рис. 5. Поперечное сечение пластины. Реконструкция формы дефекта выполнена с учетом его реальных размеров на поверхности пластины.
Будем считать, что форма реконструированного дефекта наилучшим образом соответствует форме реального дефекта, если достигается минимум положительного среднего отклонения (1) икс-компонент напряженности магнитных полей рассеяния упомянутых ранее дефектов в точках на уровне измерений, на линии, параллельной оси ОХ и проходящей через центр дефекта (рис. 6). В нашем случае центр дефекта также легко определить по топографии компонент «измеренного» магнитного поля:
(1)
где N — общее число точек на указанной в тексте линии; i — номер точки на линии; Hx0,i, Hxr,i — икс-компоненты напряженности магнитного поля в точках на линии исследуемого и реконструированного дефектов соответственно.
Рис. 6. График значений икс-компонент напряженности магнитного поля на поверхности измерений, на линии, проходящей через центр дефекта по оси ОХ. 1 — значения Hx0 исследуемого дефекта; 2 — значения Hxr реконструированного дефекта.
РЕЗУЛЬТАТЫ ИССЛЕДОВАНИЙ И ИХ ОБСУЖДЕНИЕ
На рис. 7 приведен график изменения среднего отклонения D в зависимости от смещения поверхности дефекта на величину S; S = 0 соответствует форме дефекта, полученной в результате решения обратной геометрической задачи.
Рис. 7. Зависимость среднего отклонения D от смещения по вертикали поверхности дефекта. Положительный знак S соответствует увеличению глубины дефекта.
В результате решения обратной геометрической задачи магнитостатики был реконструирован куполообразный дефект со следующими параметрами формы 25×25×2,5 мм. Глубина реконструированного дефекта отличается от глубины реального дефекта на 16,7 %. При этом величина среднего отклонения D составила 150 А/м. Уменьшение глубины дефекта на 0,25 мм приводит к росту D. Минимальное положительное значение D достигается при увеличении глубины дефекта на 0,5 мм, что соответствует геометрическим параметрам реального дефекта.
ВЫВОДЫ
С помощью предложенной методики верифицируется решение обратной геометрической задачи магнитостатики и возможно уточнение значений геометрических параметров формы исследуемого дефекта.
Л.В. Михайлов благодарит Институт физики металлов имени М.Н. Михеева за поддержку его работы по государственному заданию Минобрнауки России по теме «Диагностика», которая выполнялась в рамках молодежного проекта ИФМ УрО РАН № М8-23.
Работа выполнена в рамках государственного задания МИНОБРНАУКИ России (тема «Диагностика» №122021000030–1).
Об авторах
А. В. Никитин
Институт физики металлов им. М.Н. Михеева УрО РАН
Автор, ответственный за переписку.
Email: an@imp.uran.ru
Россия, 620108 Екатеринбург, ул. Софьи Ковалевской, 18
Л. В. Михайлов
Институт физики металлов им. М.Н. Михеева УрО РАН
Email: mikhaylov_lv@imp.uran.ru
Россия, 620108 Екатеринбург, ул. Софьи Ковалевской, 18
А. В. Михайлов
Институт физики металлов им. М.Н. Михеева УрО РАН
Email: mikhaylov_lv@imp.uran.ru
Россия, 620108 Екатеринбург, ул. Софьи Ковалевской, 18
Ю. Л. Гобов
Институт физики металлов им. М.Н. Михеева УрО РАН
Email: mikhaylov_lv@imp.uran.ru
Россия, 620108 Екатеринбург, ул. Софьи Ковалевской, 18
В. Н. Костин
Институт физики металлов им. М.Н. Михеева УрО РАН
Email: mikhaylov_lv@imp.uran.ru
Россия, 620108 Екатеринбург, ул. Софьи Ковалевской, 18
Я. Г. Смородинский
Институт физики металлов им. М.Н. Михеева УрО РАН
Email: mikhaylov_lv@imp.uran.ru
Россия, 620108 Екатеринбург, ул. Софьи Ковалевской, 18
Список литературы
- Янус Р.И. Некоторые вопросы теории магнитной дефектоскопии // Труды ИФМ УФАН СССР. 1948. Вып. 7. С. 23—37.
- Вонсовский С.В. Простейшие расчеты для задач магнитной дефектоскопии // Журнал технической физики. 1938. Т. 8. № 16. С. 1453—1467.
- Аркадьев В.К. О развитии теоретических основ дефектоскопии // Известия АН СССР. Отделение технических наук. 1937. № 2. С. 233—239.
- Зацепин Н. Н., Щербинин В.Е. К расчету магнитостатического поля поверхностных дефектов. I. Топография полей моделей дефектов // Дефектоскопия. 1966. № 5. С. 50—58.
- Ферстер Ф. Неразрушающий контроль методом магнитных полей рассеяния. Теоретические и экспериментальные основы выявления поверхностных дефектов конечной и бесконечной глубины // Дефектоскопия. 1982. № 11. С. 3—25.
- Загидулин Р.В., Мужицкий В.Ф., Курозаев В.И. Расчет магнитостатического поля внутреннего дефекта и дефекта внутренней поверхности в ферромагнитной пластине. II Магнитное поле дефекта в воздухе. Нелинейная среда // Дефектоскопия. 1997. № 1. С. 55—62.
- Халилеев П.А., Патраманский Б.В., Лоскутов В.Е., Корзунин Г.С., Щербинин В.Е.. Выявляемость дефектов в трубопроводах из различных марок стали в зависимости от их конфигурации // Дефектоскопия. 2000. № 8. С. 22—33.
- Кротов Л.Н. Реконструкция границы раздела сред по пространственному распределению магнитного поля рассеяния. II. Постановка и метод решения обратной геометрической задачи магнитостатики // Дефектоскопия. 2004. № 6. С. 36—44.
- Дякин В.В., Кудряшова О.В., Раевский В.Я. Обратная задача магнитостатики в полях насыщения // Дефектоскопия. 2019. № 10. С. 35—44.
- Гобов Ю.Л., Никитин А.В., Попов С.Э. Решение обратной геометрической задачи магнитостатики для дефектов коррозии с учетом нелинейных свойств ферромагнетика // Дефектоскопия. 2018. № 12. С. 31—37.
- Никитин А.В., Гобов Ю.Л., Михайлов А.В., Михайлов Л.В. Методика решения обратной геометрической задачи магнитостатики для поверхностных дефектов магнитомягкого ферромагнетика // Дефектоскопия. 2022. № 58. С. 24—34.
- Никитин А.В., Михайлов А.В., Михайлов Л.В., Гобов Ю.Л., Костин В.Н., Смородинский Я.Г. Область применимости методики построения линий магнитной индукции для дефектометрии протяженных объектов // Дефектоскопия. 2023. № 59. С. 51—59.
- Specifications and requirements for in-line inspection of pipelines. 2016. https://pipelineoperators.org/
Дополнительные файлы










