FEMM 5

Часть VI · Электрические аппараты → § 38 · Устройство электромагнита и нелинейная магнитная цепь

Занятие 38.2

Нелинейная магнитная цепь: заданы ток или требуемый поток

Цель и задачи

Цель. Решить прямую и обратную задачи нелинейной магнитной цепи.

Задачи

  • Вывести требуемую МДС по заданному потоку и табличной H(B).
  • Решить нелинейную цепь при заданном токе без ложной постоянной \mu.
Что нужно вспомнить перед началом

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

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

Расчётная лабораторная работа · Построение карты отклика и локальных параметров

Инженерный вопрос. Как насыщение конкретной детали меняет потребную МДС и зависимость потокосцепления от тока?

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

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

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

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

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

  2. Ток для B=1{,}6 Тл в условной стали. Проверьте вычисление локального параметра по соседним точкам и его зависимость от шага.

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

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

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

Результат этого занятия. Характеристика \Psi (I) и сравнение секущей и дифференциальной индуктивностей с условной B(H) и объяснением расхождений цепи и поля.

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

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

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

Проверка из практикума. Кривая названа условной; значения H(B) не подменены постоянной \mu; различие цепной и полевой интерполяций учтено.

Независимый контроль. Сопоставьте прогноз с прямым расчётом в независимой точке и измените шаг производной или плотность карты отдельно от сетки.

Критерий выполнения. Кривая названа условной; значения H(B) не подменены постоянной \mu; различие цепной и полевой интерполяций учтено.

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

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

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

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

Границы вывода. Магнитная сила отдельного привода не подтверждает работу всего аппарата: контактные усилия, пружины, трение, динамика и нагрев должны быть заданы и проверены по назначению.

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

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

Сечение каждой детали определяет её рабочую точку

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

В учебной серии применена собственная монотонная B(H), без гистерезиса и частотных потерь. Например, B=1{,}45; 1,60; 1,75 Тл соответствуют H=600; 1400; 4500 А/м. Эти числа задают упражнение и не являются протоколом промышленной стали. Интерполируем H по B линейно между узлами; интерполяция FEMM может отличаться.

В нативной серии П-магнита используется своя условная B(H): например, при B=1{,}45 Тл задано H=600 А/м. В серии П6 при той же B задано H=1000 А/м. Это разные учебные материалы, не две оценки одной марки стали; сравнивать силы или NI между сериями как чистый эффект алгоритма нельзя.

Обратная задача: сколько нужно ампер-витков

Если поток задан, сначала находят индукции B_j=\Phi /S_j в каждой последовательной детали. Затем по характеристике материала получают H_j и суммируют падения магнитного потенциала. В линейных зазорах H=B/\mu _{0}. Такой расчёт не требует подбирать постоянную \mu всего сердечника.

NI=\sum_j H_j(\Phi/S_j)\ell_j+\sum_k\frac{\Phi\delta_k}{\mu_0S_k}
(38.4)
Формула относится к последовательному пути без разветвления потоков; для утечек нужны отдельные ветви.
  1. Для равных сечений, двух зазоров и принятой средней длины магнитного пути в ферромагнитном участке запишем уравнение одной неизвестной.

    f(B)=2\delta B/\mu_0+\ell_{\rm ж}H(B)-NI=0
    (38.5)
  2. При положительной производной H(B) левая часть монотонна, поэтому бисекция на охватывающем интервале имеет однозначный корень.

    f^{\prime}(B)=2\delta/\mu_0+\ell_{\rm ж}\,\mathrm dH/\mathrm dB>0
    (38.6)

Прямая задача и наклоны характеристики

При заданном I ищем B, затем \Psi=NBS. Секущая L ст=\Psi /I описывает отношение величин в рабочей точке. Для малого приращения тока нужна L дин=d\Psi /dI, зависящая от локального dH /dB. В нелинейном участке эти величины различны даже при неподвижном якоре.

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

L_{\rm дин}=\frac{N^2S}{2\delta/\mu_0+\ell_{\rm ж}\,\mathrm dH/\mathrm dB}
(38.7)
Получено дифференцированием цепного уравнения при фиксированной геометрии и общем сцепленном потоке.

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

Ток для B=1{,}6 Тл в условной стали

Дано. Два зазора по 1 мм, \ell_{\rm ж}=0{,}20\,\text{м}, N=300, B=1{,}6\,\text{Тл}, табличное H=1400\,\text{А/м}.

  1. Падение на железе равно 1400\cdot 0{,}2=280 А; в двух зазорах 2\cdot 0{,}001\cdot 1{,}6/\mu _{0}=2546{,}479089 А.

  2. Складываем МДС и делим на 300 витков. Наибольшая доля по-прежнему в воздухе, но железо уже нельзя исключать.

Ответ: NI=2826{,}479089\,\text{А}; I=9{,}421596965\,\text{А}. При 10 А цепная интерполяция даёт B\approx1{,}63031\,\text{Тл}.

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

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

Сравните линейную и нелинейную серии одного устройства

Таблица 38.4. Исходные данные практикума
Исходная величинаЗначениеЕдиница
Семействоu-core
Каждый зазор1мм
Нелинейные токи1; 3; 6; 10А
Средняя длина цепной оценки200мм
  1. Из geometry-contract.json перенесите всю таблицу B(H), не заменяя её тремя приведёнными точками. Проверьте монотонность и единицы.
  2. Решите цепное уравнение бисекцией при 1, 3, 6, 10 А. Сохраните таблицу B, \Psi, долей МДС в железе и воздухе.
  3. Встроенным reproduce.py получите соответствующие нелинейные поля. Сопоставьте нормальный поток сечения и \Psi цепи с цепной оценкой.
  4. Уточните сетку в самом насыщенном случае. Для выбранного тока возьмите независимые соседние точки и оцените d\Psi /dI; половинный шаг проверяется отдельно.

Ожидаемый результат. График \Psi (I), таблица секущей и динамической L, объяснение начала насыщения конкретной детали.

Как проверить. Кривая названа условной; значения H(B) не подменены постоянной \mu; различие цепной и полевой интерполяций учтено.

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

Как насыщение конкретной детали меняет потребную МДС и зависимость потокосцепления от тока?

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

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

Итог работы. Характеристика \Psi (I) и сравнение секущей и дифференциальной индуктивностей с условной B(H) и объяснением расхождений цепи и поля.

Встроенный MicroPython: Нелинейная магнитная цепь: заданы ток или требуемый поток · Порядок запуска

Найти ток по требуемой индукции и условной B(H)

Выполните короткий контроль во встроенном Python FEMM 5. Он проверяет формулу и единицы. Полевые модели строит отдельный reproduce.py из материалов занятия; сравнивайте только одинаковые параметры и определения. Печатный вывод этой программы не является результатом решения поля.

Python · встроенный MicroPython FEMM 5
"""Найти ток по требуемой индукции и условной B(H).

Среда: встроенный MicroPython FEMM 5.
Самостоятельный аналитический контроль; полевой решатель не вызывается.
"""
import math
mu0 = 4 * math.pi * 1e-7
# Точка взята из условной учебной таблицы: B в Тл, H в А/м.
b, h = 1.6, 1400.0
n, gap, iron_length = 300, 1e-3, 0.20
mmf_air = 2 * gap * b / mu0
# В обратной задаче падения МДС складываются до деления на число витков.
mmf_iron = iron_length * h
print("МДС зазоров, А:", mmf_air)
print("МДС железа, А:", mmf_iron)
print("Требуемый ток, А:", (mmf_air + mmf_iron) / n)

Место занятия в расчёте магнитной системы

Нелинейная МДС стали и решение магнитной цепи. Найдите рабочую точку по H(B) и оцените, какую долю полной МДС расходует сталь при выбранном потоке или токе.

Опорная теория: 8.1 · Материальные уравнения B(H); 38.1 · Геометрия участков магнитной цепи.

Применение: 44.4 · Разбить насыщенный полюс на участки; 48.1 · Получить нелинейную карту реле.

Маршрут магнитных систем в основном оглавлении · Каталог рассчитанных конструкций

Расчётный образец FEMM 5

Нелинейный режим собственной П-образной модели: N=300, I=6 А, два зазора по 1 мм. Использована условная B(H), не протокол марки стали; простая аналитическая цепь занятия имеет отдельные исходные допущения. Откройте модель, воспроизведите указанный режим и сопоставьте свои измерения с таблицей результатов.

Таблица 38.5. Результаты расчётного образца
ВеличинаРезультатЕдиница
Сила притяжения192,5Н
Потокосцепление0,09429Вб
Коэнергия0,2912Дж
Конечные элементы636901
Самостоятельные конструкции, размеры в миллиметрах. Схема построена по паспорту геометрии; это не снимок интерфейса FEMM.
Рисунок 38.2. Самостоятельные конструкции, размеры в миллиметрах. Схема построена по протоколу геометрии; это не снимок интерфейса FEMM.

Используются линейная сталь и отдельная условная B(H). Серии и проверка силы приведены в исследовании электромагнитов.

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

Сравните один и тот же магнит при заданном I и при требовании заданного \Phi; покажите, какая часть МДС уходит в железо и какая в воздух.

Критерий выполнения: Кривая названа условной; значения H(B) не подменены постоянной \mu; различие цепной и полевой интерполяций учтено.
Первоисточники и дополнительное чтение по занятию
  1. [100] — С. 30–43 (печатная нумерация). Полная предоставленная копия изучена отдельно; ссылка ведёт на фрагмент издателя. Изложение и численные примеры курса созданы самостоятельно.
  2. [33] — Магнитная плоская и осесимметричная задачи, цепи, интегралы силы, встроенный Python.
Пояснения и типичные ошибки
Почему увеличение тока не гарантирует пропорциональный рост потока?

При росте H в насыщаемых деталях увеличивается магнитное сопротивление; зависимость определяется B(H), зазорами и сечениями.

Можно ли получить поле заданного потока, просто назначив каждому участку одинаковую B?

Нет: в цепной оценке сначала учитывают разные площади, а поле должно удовлетворять уравнениям, источникам и граничным условиям. Заданная B всей детали не является обычным источником FEMM.

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

  • Условная B(H) не подтверждает свойства марки стали.
  • Секущая \mu =B/H не заменяет дифференциальную проницаемость в расчёте приращений.