Занятие 57.1

Намагничивание, насыщение и ориентация дефекта

Цель и задачи

Цель. Оценить влияние намагничивания и насыщения на сигнал дефекта.

Задачи

Рассчитать секущую и локальную дифференциальную проницаемости по заданным точкам.
Сравнить нормированный сигнал при нескольких уровнях намагничивания.

Что нужно вспомнить перед началом

Перед началом: 56.2 · Магнитное поле рассеяния над пазом ферромагнитного образца, 54.2 · Измерение B–H-характеристик и магнитных потерь

Если одно из понятий незнакомо, вернитесь к предыдущему занятию и воспроизведите его численный пример. После этого переходите к новой постановке.

Место занятия в инженерной методике

Расчётная лабораторная работа · Оценка чувствительности к дефекту и мешающим факторам

Инженерный вопрос. Можно ли сравнивать сигналы без учёта рабочего состояния материала?

Разностный сигнал и его зависимость от дефекта и зазора датчика позволяют оценить чувствительность схемы контроля и выбрать признак для проверенного диапазона геометрии.

Последовательность решения, отчёт и критерии результата

Результат получают из наблюдаемых величин и проверяют в условиях данного практикума.

Последовательность решения

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

Что включить в отчёт

Результат этого занятия. Сравнение сигналов с объяснением роли насыщения и ориентации паза.

Данные именно этого опыта. Сопоставьте возбуждение, рабочие точки и разностные профили линейной и нелинейной моделей.

Оформление наблюдений. Базовый и дефектный сигналы, их разность, координата сканирования, режим и один мешающий фактор с единицами.

Как обосновать результат

Проверка из практикума. Условная кривая не приписана марке стали; нормировка и контрольная точка заданы; линейное масштабирование и нелинейный режим различены.

Независимый контроль. Используйте согласованный бездефектный расчёт, отдельное изменение мешающего параметра и устойчивость выбранного признака.

Критерий выполнения. Условная кривая не приписана марке стали; нормировка и контрольная точка заданы; линейное масштабирование и нелинейный режим различены.

Точность и расхождения. Сравнивайте разностные сигналы на одинаковых координатах; сетку дефекта и шаг скана уточняйте раздельно. Численная различимость не задаёт вероятность обнаружения без шума и натурной валидации.

Основание сравнения. Опора — модель сигнала и геометрия дефекта из занятия. Плоский паз, кольцевое утонение и локальная трещина имеют разную размерность; амплитуда сигнала сама по себе не определяет размер дефекта.

Вывод для устройства

Обоснуйте пригодность признака для обнаружения или сравнения дефектов в проверенном диапазоне, сохраняя ограничения двумерной модели и датчика.

Границы вывода. Численная аномалия поля не равна доказанной обнаруживаемости дефекта. Геометрия, датчик, шум, свойства и натурные образцы необходимы для перехода к характеристикам средства контроля.

Теория и вывод формул

Режим намагничивания и принятая характеристика B(H) определяют разностный сигнал магнитного контроля. Исследование локального поля показывает, почему ток возбуждения или индукция вдали от образца не дают универсальной оценки его насыщения.

Соглашения об обозначениях и единицах

\(i(t)\), \(u(t)\) — мгновенные значения; \(\widehat I\), \(\widehat U\) — комплексные амплитуды; \(\underline I_{\text{д}}\), \(\underline U_{\text{д}}\) — комплексные действующие значения. Для синусоиды \(\widehat I=\sqrt{2}\,\underline I_{\text{д}}\); действующее значение тока равно \(|\underline I_{\text{д}}|\). Нормировка обязательно учитывается в формулах мощности.
\(\Phi=\int_S\mathbf B\cdot\mathrm d\mathbf S\) — поток через ориентированную поверхность; \(\Psi=\sum_k\oint_{C_k}\mathbf A\cdot\mathrm d\mathbf l\) — потокосцепление последовательных витков с согласованными замкнутыми контурами. Только при одинаковом сцепленном потоке каждого витка \(\Psi=N\Phi\). Обе величины измеряются в веберах: число витков безразмерно. Запись «Вб·виток» подчёркивает суммирование.
\(\theta_{\text{м}}\), \(\Omega\) (также \(\omega_{\text{мех}}\)) — механические угол и угловая скорость. Для основной гармоники \(\theta_{\text{э}}=p\theta_{\text{м}}\), \(\omega_{\text{э}}=p\Omega\), где \(p\) — число пар полюсов. Геометрический период зубцов задаётся отдельно. В производных по углу используйте радианы.
Статическая индуктивность \(L_{\text{ст}}=\Psi/i\) — наклон секущей из начала координат; это определение используется при \(\Psi(0)=0\) и \(i\ne0\). Динамическая индуктивность \(L_{\text{дин}}=\partial\Psi/\partial i\) — наклон касательной (дифференциальная, инкрементальная); положение и остальные токи фиксированы. Для линейной характеристики без смещения обе равны одной индуктивности \(L\). При постоянных магнитах и взаимном возбуждении сначала учитывают ненулевое смещение потокосцепления.
Положительные направления, оси и специальные индексы определяются в конкретной задаче. Единицу берите из условия или заголовка таблицы; перед подстановкой в формулу СИ переводите миллиметры в метры, миллигенри в генри и градусы в радианы.
Численные результаты обычно показаны с четырьмя значащими цифрами; для сравнения близких значений в таблицах сходимости сохранено больше цифр. Полные значения остаются в исходных расчётных данных. Малые остатки вычисления координат на осях показаны как ноль; малые физические величины сохраняются в показательной записи.

Указатель формул этого занятия

(57.1) · Режим задаётся источником, состояние проверяется полем

Режим задаётся источником, состояние проверяется полем

В реальном магнитном контроле режим зависит от ярма, катушки, контакта и магнитной цепи. В учебной модели источник заменён граничным полем. Поэтому сохраняйте B внеш и дополнительно контролируйте B и H в заранее выбранных точках образца. Одно значение источника не доказывает одинаковую степень насыщения разных образцов.

Сравните профиль в линейном материале с серией условной нелинейной характеристики. Линейный ответ при фиксированной геометрии масштабируется с возбуждением; изменение нормированной формы или амплитуды в нелинейной серии связано с перераспределением проницаемости, а также должно быть отделено от численной ошибки.

\mu_{\mathrm{сек}}=\frac{B}{H},\qquad \mu_{\mathrm{диф}}=\frac{\mathrm dB}{\mathrm dH},\qquad \widehat S=\frac{S_y}{B_{\mathrm{внеш}}}

(57.1)

Обе проницаемости имеют единицу Гн/м; относительные значения получаются делением на \mu _{0}. Нормированный сигнал безразмерен.

Условная характеристика и насыщение

В опубликованных материалах задана авторская таблица B(H): (0;0), (0,5 Тл;100 А/м), (1;200), (1,4;400), (1,6;1000), (1,8;5000), (2;20000), (2,2;80000). Она нужна для изучения влияния формы характеристики и не объявлена измерением конкретной стали. Полная таблица сохраняется рядом с моделью.

При больших B рост H становится резким; дифференциальная проницаемость уменьшается. Поле около паза и в оставшейся перемычке может переходить в этот режим раньше выбранной контрольной точки вдали от паза. Проверяйте распределение, а не единственный цвет максимума. Сравните локальные B с верхней табличной точкой 2,2 Тл: выход за неё требует отдельного обоснования продолжения характеристики; контрольная точка вдали от паза не ограничивает всё поле. Насыщение не означает, что B перестаёт изменяться абсолютно.

Ориентация намагничивания

Паз наиболее явно нарушает путь потока, когда имеет существенную компоненту поперёк него. Если поток идёт вдоль протяжённого дефекта, контраст другого сечения может быть слабее и иной формы. Это качественный прогноз; величину нельзя получить простым косинусом без обоснованной модели.

Текущий плоский образец представляет паз, продолжающийся вдоль третьей координаты. Поворот вектора в плоскости меняет направление к поверхности и геометрию магнитной цепи одновременно. Поворот реальной локальной трещины в плоскости поверхности обычно требует 3D описания. В отчёте различите эти операции.

Гистерезис не появляется из одной B(H)-кривой

Однозначная нелинейная магнитостатическая характеристика не хранит магнитную историю, остаточную намагниченность и цикл размагничивания. Поэтому расчёт нескольких уровней источника не является петлёй гистерезиса. Не используйте такие модели для заявления о повторяемости реальной процедуры после отключения намагничивания.

Параметрическое исследование нужно завершить таблицей: источник, B/H в контрольной точке, размах сигнала, нормированный размах, тип материала, сетка. Вывод об увеличении либо уменьшении чувствительности формулируйте только для рассчитанного диапазона и выбранной компоненты.

Разобранный численный пример

Разбор измерительной модели и чисел

Дано. Локальный участок условной характеристики: B_{1} = 1{,}4 Тл при H_{1} = 400 А/м, B_{2} = 1{,}6 Тл при H_{2} = 1000 А/м.

Разностная оценка \mu диф = (B_{2} - B_{1})/(H_{2} - H_{1}) = 0{,}2/600 Гн/м.
Делим на \mu _{0} = 4\pi \cdot 10^{-7} Гн/м.
Получаем относительную дифференциальную проницаемость 265,258. Она не равна B_{2}/(\mu _{0}H_{2}) \approx 1273{,}24.

Ответ: Относительная дифференциальная проницаемость ≈ 265,258.

Проверка: Это секущая оценка производной на заданном интервале, не лабораторное измерение и не точная производная интерполяции FEMM.

Расчётный образец магнитной дефектоскопии

Самостоятельная пара плоских моделей MFL: сохранение, повторное открытие и 16 решений FEMM 5 с русскими именами свойств. Линейная базовая серия проверена на высотах 1, 2 и 4 мм; нелинейные варианты сохранены для дальнейшей проверки.

Таблица 57.1. Сохранённый разностный сигнал и проверка сетки
Величина	Результат	Единица
Размах \Delta By при h = 1 мм	243,8	мкТл
Размах \Delta By при h = 2 мм	97,76	мкТл
Размах \Delta By при h = 4 мм	34,02	мкТл
Сеточное изменение размаха при h = 1 мм	4,314	%
Сеточное изменение размаха при h = 2 мм	1,288	%
Сеточное изменение размаха при h = 4 мм	1,503	%

Самостоятельная учебная модель MFL. Сходимость размаха проверяется на всех трёх высотах линейной базовой серии. Нелинейные варианты служат отдельным исследованием на финальной сетке; для них индивидуальная сходимость не заявляется. Расчётная проверка не является физической аттестацией дефектоскопа.

Практикум в FEMM 5

Сравните намагничивание и разностный сигнал в линейном и условном нелинейном образцах

Таблица 57.2. Исходные данные практикума
Исходная величина	Значение	Единица
Условная характеристика	8 точек B(H)	Тл; А/м
Нелинейные режимы B внеш	0,01; 0,1; 0,2; 0,25	Тл

Откройте сохранённые нелинейные FEM; проверьте таблицу B(H), нулевую частоту и отсутствие собственной проводимости.
По result.json составьте таблицу B внеш, B_{\mathrm{sample\,T}} и Sy(h = 1 мм) для каждой пары целого образца и паза.
Нормируйте Sy на B внеш. Объясните отличие линейного масштабирования от нелинейной серии; укажите разрешимость сигнала относительно сеточной разности.
В отдельной новой серии проверьте промежуточный режим, не меняя характеристику и геометрию. Исходные решения и квитанцию не перезаписывайте.
Сформулируйте, какие измерения требуются для переноса условной кривой на реальный образец и что останется неизвестным о магнитной истории.

Ожидаемый результат. Таблица связывает режим источника с состоянием материала; выводы о насыщении сопровождаются распределением B и областью применимости.

Как проверить. Условная кривая не приписана марке стали; нормировка и контрольная точка заданы; линейное масштабирование и нелинейный режим различены.

Обработка и анализ результатов

Можно ли сравнивать сигналы без учёта рабочего состояния материала?

Сопоставьте возбуждение, рабочие точки и разностные профили линейной и нелинейной моделей.
Вычтите согласованный бездефектный результат, сравните разностные сигналы и влияние одного мешающего параметра.

С чем сравнить. Опора — модель сигнала и геометрия дефекта из занятия. Плоский паз, кольцевое утонение и локальная трещина имеют разную размерность; амплитуда сигнала сама по себе не определяет размер дефекта.

Итог работы. Сравнение сигналов с объяснением роли насыщения и ориентации паза.

Программа с русскими пояснениями: Намагничивание образца · Порядок запуска

Обработка во встроенном Python FEMM 5

Сохраните программу и выполните встроенным Python FEMM 5. Этот фрагмент обрабатывает явно заданные данные; полевые модели рассчитываются отдельной программой со страницы работы.

"""Проницаемость на участке условной учебной характеристики."""
import math
# Отдельно задаём приращения, чтобы не спутать B/H с производной.
delta_b_T = 0.2
delta_h_A_m = 600.0
# Относительная величина получается делением на магнитную постоянную.
mu_r_diff = (delta_b_T / delta_h_A_m) / (4 * math.pi * 1e-7)
print("Относительная дифференциальная проницаемость:", mu_r_diff)

Самостоятельная работа

Сравните намагничивание и разностный сигнал в линейном и условном нелинейном образцах. Сохраните файлы и таблицу; отделите модельный вывод от экспериментального.

Критерий выполнения: Условная кривая не приписана марке стали; нормировка и контрольная точка заданы; линейное масштабирование и нелинейный режим различены.

Первоисточники и дополнительное чтение по занятию

[33] — Редакция от 1 октября 2026 г.; главы по магнитным задачам и команды mi_/mo_
[60] — 2010; влияние геометрии и ориентации паза; переход от 2D к 3D. Учебная модель курса разработана независимо и не воспроизводит установку статьи.

Пояснения и типичные ошибки

Почему увеличение возбуждения не гарантирует пропорциональный рост сигнала?

Нелинейная проницаемость меняет распределение поля в образце и около паза; нормированный размах нужно вычислить и проверить.

Почему три магнитостатических состояния не являются петлёй гистерезиса?

Однозначная B(H)-кривая не хранит историю и не задаёт разные ветви намагничивания и размагничивания.

На что обратить внимание

Нельзя объявлять условную B(H)-таблицу протоколом реального материала.
Нельзя считать секущую проницаемость равной дифференциальной в нелинейной области.