FEMM 5

Часть VIII · Магнитная и вихретоковая дефектоскопия → § 58 · Вихретоковый контроль

Занятие 58.1

Импеданс вихретокового датчика: частота, проводимость и зазор

Цель и задачи

Цель. Выделить комплексный отклик вихретокового датчика.

Задачи

  • Вычислить глубину проникновения и компоненты разностного импеданса.
  • Получить Z из свойств цепи при одинаковом соглашении об амплитудах.
Что нужно вспомнить перед началом

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

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

Расчётная лабораторная работа · Согласование комплексного отклика и средней мощности

Инженерный вопрос. Как изменение зазора влияет на комплексный импеданс и отличие от дефекта?

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

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

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

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

  1. Обосновать постановку и способ измерения. Зафиксируйте частоту, фазовый нуль и амплитудную либо действующую нормировку; отдельно запишите действительную и мнимую части уравнения.

  2. Разбор измерительной модели и чисел. Проверьте комплексное произведение для мощности и физический смысл фазы на численном примере.

  3. Исследуйте влияние зазора на комплексный импеданс проверенного вихретокового датчика. Выполните действия практикума по его исходной таблице. Сохраните частоту, зазор, сопротивление и реактивную часть импеданса согласованных опытов.

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

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

Результат этого занятия. Комплексная диаграмма либо таблица с раздельным анализом компонент и мешающего зазора.

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

Оформление наблюдений. Частота, \operatorname{Re} и \operatorname{Im} входных и выходных величин, единицы, нормировка и средняя мощность; модуль без фазы недостаточен.

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

Проверка из практикума. Влияние зазора датчика отделено от влияния дефекта. Глубина проникновения приведена с условиями применимости, комплексное сопротивление — с единицами и частотой.

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

Критерий выполнения. Влияние зазора датчика отделено от влияния дефекта. Глубина проникновения приведена с условиями применимости, комплексное сопротивление — с единицами и частотой.

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

Основание сравнения. Опора — комплексные уравнения и определение средней мощности из занятия. Не смешивайте амплитудные и действующие значения и не заменяйте комплексную невязку одной разностью модулей.

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

Объясните, как фазовый отклик влияет на параметр цепи, реактивную составляющую или потери в исследуемом режиме.

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

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

Отклик катушки над проводящей пластиной характеризуется комплексным сопротивлением. Его зависимость от частоты и зазора рассматривается совместно с глубиной проникновения и балансом мощности на проверенном расчётном примере.

Соглашения об обозначениях и единицах
  • \(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\). При постоянных магнитах и взаимном возбуждении сначала учитывают ненулевое смещение потокосцепления.
  • Положительные направления, оси и специальные индексы определяются в конкретной задаче. Единицу берите из условия или заголовка таблицы; перед подстановкой в формулу СИ переводите миллиметры в метры, миллигенри в генри и градусы в радианы.
  • Численные результаты обычно показаны с четырьмя значащими цифрами; для сравнения близких значений в таблицах сходимости сохранено больше цифр. Полные значения остаются в исходных расчётных данных. Малые остатки вычисления координат на осях показаны как ноль; малые физические величины сохраняются в показательной записи.
Указатель формул этого занятия
  1. (58.1) · Индукционное возбуждение и измеряемая цепь
  2. (58.2) · Глубина проникновения как масштаб

Индукционное возбуждение и измеряемая цепь

Переменное поле катушки создаёт вихревые токи в мишени; их обратное поле меняет потокосцепление и напряжение катушки. Для линейной гармонической задачи Z = U/I. Сравнение выполняется при одинаковой частоте, геометрии обмотки и определении тока. В нелинейном материале дополнительно важен режим возбуждения.

В сохранённой учебной модели обмотка идеальна и не имеет собственной проводимости. Её активный входной отклик связан с потерями мишени; сопротивление провода, кабель, ёмкость и измерительная электроника не учтены. Поэтому полученное \operatorname{Re} Z нельзя без добавления этих элементов сравнивать с показанием реального LCR-метра.

Z=\frac{U}{I}=R+jX,\qquad \Delta Z=Z_{\mathrm{иссл}}-Z_{\mathrm{эт}},\qquad P=\frac12\Re(UI^*)
(58.1)
U и I — комплексные амплитуды; при использовании действующих значений множитель 1/2 не нужен.

Глубина проникновения как масштаб

Для однородного линейного проводящего полупространства характерная глубина проникновения \delta зависит от частоты, проницаемости и проводимости. На глубине \delta амплитуда простого плоского поля уменьшается в e раз. Конечная катушка и конечная стенка имеют неоднородное поле; \delta служит ориентиром выбора частоты и сетки, а не точной глубиной обнаруживаемого дефекта.

Если толщина стенки сопоставима с \delta, изменение толщины может существенно менять токовый путь. При высокой частоте исследуется преимущественно поверхность. В ферромагнетике проницаемость зависит от режима и может быть комплексной; подстановка постоянной \mu _{r} без обоснования даёт только условную оценку.

\delta=\sqrt{\frac{2}{\omega\mu\sigma}},\qquad \omega=2\pi f
(58.2)
\delta — м; f — Гц; \mu — Гн/м; \sigma — См/м. Формула относится к линейному однородному проводнику и гармоническому полю.

Проверенный датчик зазора — основа, а не модель трещины

Катушка среднего радиуса 20 мм, сечения 4×6 мм, 100 витков и амплитуды 1 А работает на 1000 Гц. Мишень — немагнитный диск радиуса 60 мм, толщины 8 мм, \sigma = 35 МСм/м. В опубликованных материалах посчитаны зазоры 2, 5 и 10 мм на трёх сетках. Эти исходные модели и их квитанции сохраняются без переименования и перезаписи.

В этом опыте изменяется зазор между катушкой и образцом, а дефект не вводится. Расчёт позволяет оценить влияние зазора (lift-off) как мешающего фактора и обосновать калибровку положения. Для исследования утонения далее используется отдельная пара моделей: трубы без дефекта и трубы с уменьшенной толщиной стенки. Обнаружение локальной трещины требует иной допустимой постановки.

Комплексная плоскость и мешающие параметры

Изменение зазора, проводимости и дефекта может давать разные направления \Delta Z в комплексной плоскости. Два канала \operatorname{Re} Z и \operatorname{Im} Z могут помочь разделению, если направления изменений достаточно независимы и их амплитуды превышают неопределённость. Это проверяется локальным анализом чувствительности и набором вариантов.

Частотный ряд добавляет наблюдения, но одновременно требует нового решения поля и нового контроля сетки по \delta. Частоты не должны использоваться как дополнительные независимые точки одной калибровки без учёта общего неизвестного материала, температуры и геометрии. Для нового ряда сохраняйте новые каталоги и квитанции.

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

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

Дано. Немагнитный проводник: \sigma = 35\cdot 10^{6} См/м, \mu = \mu _{0}, f = 1000 Гц.

  1. Вычислите \omega = 2\pi \cdot 1000 рад/с.

  2. Подставьте \mu _{0} = 4\pi \cdot 10^{-7} Гн/м и \sigma в формулу \delta.

  3. \delta \approx 0{,}00269021 м = 2,69021 мм. Толщина исходного диска 8 мм составляет около 2{,}97\delta.

Ответ: Глубина проникновения \delta \approx 2{,}69021 мм.

Проверка: Оценка полупространства не является результатом FEMM и не задаёт вероятность обнаружения дефекта.

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

Исследуйте влияние зазора на комплексный импеданс проверенного вихретокового датчика

Таблица 58.1. Исходные данные практикума
Исходная величинаЗначениеЕдиница
Частота1000Гц
Проводимость диска35МСм/м
Зазоры2; 5; 10мм
Ток: амплитуда1А
  1. Откройте модели работы eddy-current-sensor; проверьте тип Axisymmetric, \sigma в МСм/м и связь обмотки с цепью.
  2. На финальной сетке каждого зазора получите ток и напряжение цепи; вычислите R, X и кажущуюся L = X/\omega.
  3. Сопоставьте среднюю входную активную мощность 0,5 Re(UI*) с объёмными потерями мишени.
  4. Постройте точки Z на комплексной плоскости и отдельно графики R(h), X(h); не соединяйте три точки как доказанную точную калибровочную кривую.
  5. Сформулируйте отдельный новый опыт по частоте или \sigma; изменяйте один параметр за серию и заново контролируйте сетку.

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

Как проверить. Исходные квитанции сохранены; зазор не назван дефектом; число \delta приведено с условиями применимости, а Z — с единицами и частотой.

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

Как изменение зазора влияет на комплексный импеданс и отличие от дефекта?

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

С чем сравнить. Опора — комплексные уравнения и определение средней мощности из занятия. Не смешивайте амплитудные и действующие значения и не заменяйте комплексную невязку одной разностью модулей.

Итог работы. Комплексная диаграмма либо таблица с раздельным анализом компонент и мешающего зазора.

Программа с русскими пояснениями: Вихретоковый импеданс · Порядок запуска

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

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

Python · встроенный MicroPython FEMM 5
"""Масштаб проникновения для немагнитного полупространства."""
import math
# Единицы СИ: частота в Гц, проводимость в См/м, проницаемость в Гн/м.
frequency_Hz = 1000.0
conductivity_S_m = 35e6
mu_H_m = 4 * math.pi * 1e-7
# В миллиметры переводим только итоговую длину.
delta_mm = math.sqrt(2 / (2 * math.pi * frequency_Hz * mu_H_m * conductivity_S_m)) * 1000
print("Глубина проникновения, мм:", delta_mm)

Расчётные исследования по теме

Изучите физическую постановку, выполните программу во встроенном Python FEMM 5 и проверьте результаты по сохранённым моделям и независимым критериям.

  • Вихретоковый датчик зазора над проводящей пластиной

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

    Решений FEMM: 9 · Численных проверок: 32

Все темы учебного практикума

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

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

Критерий выполнения: Исходные квитанции сохранены; зазор не назван дефектом; число \delta приведено с условиями применимости, а Z — с единицами и частотой.
Первоисточники и дополнительное чтение по занятию
  1. [33] — Редакция от 1 октября 2026 г.; главы по магнитным задачам и команды mi_/mo_
  2. [10] — Катушка над пластиной с трещиной: изменение сопротивления и индуктивности. Служит примером внешнего эталона, а не квитанцией моделей курса.
Пояснения и типичные ошибки
Почему \operatorname{Re} Z положительна для пассивной проводящей мишени?

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

Почему импеданс не даёт глубину дефекта сам по себе?

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

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

  • Не смешивайте амплитуды с действующими значениями в балансе мощности.
  • Не называйте X/\omega статической индуктивностью нелинейного магнитопровода.