FEMM 5

Часть III · Встроенный Python и воспроизводимая автоматизация → § 17 · Общие нелинейные карты, индуктивности и энергетические производные

Занятие 17.1

Повтор катушки и электромагнита в реальной геометрии

Цель и задачи

Цель. Использовать карты катушки и электромагнита для сравнения режимов.

Задачи

  • Автоматизировать базовую магнитостатическую серию и отделять геометрию сердечника от цепных поправок.
  • Сохранить карты \Psi(I,\delta) и B_\delta(I,\delta) для геометрически согласованной серии.
Что нужно вспомнить перед началом

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

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

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

Инженерный вопрос. Как ток и зазор меняют полезную характеристику устройства?

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

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

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

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

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

  2. Вычисленный аналитический эталон. Проверьте вычисление локального параметра по соседним точкам и его зависимость от шага.

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

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

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

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

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

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

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

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

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

Критерий выполнения. Есть карта, независимый ориентир и сеточная проверка; критерий 10 \,\% обозначен учебным и не превращён в отраслевую норму.

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

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

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

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

Границы вывода. Параметр относится к указанной рабочей точке и определению индуктивности. Торцевые поля, сопротивление и нагрев не возникают автоматически из магнитостатической энергии.

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

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

Соглашения об обозначениях и единицах
  • \(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. (17.1) · Физическая и вычислительная постановка
  2. (17.2) · Численный пример

Физическая и вычислительная постановка

Для первой оценки электромагнита магнитную цепь разделяют на сердечник и воздушный зазор. МДС определяется циркуляцией \oint\mathbf H\cdot\mathrm d\mathbf l. Для выбранного магнитостатического контура закон полного тока даёт NI как сумму магнитных напряжений — произведений H на длины соответствующих участков. Если сердечник имеет большую проницаемость и одинаковая площадь потока допустима, основная магнитодвижущая сила приходится на зазор: B_\delta \approx \mu_0NI/\delta. Затем \Phi = B_\delta S и \Psi = N\Phi, откуда L \approx \mu_0 N^2 S/\delta. Этот вывод явно показывает зависимость от площади, числа витков и зазора.

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

NI\approx H_\delta\delta=\frac{B_\delta\delta}{\mu_0},\quad L\approx\frac{\mu_0N^2S}{\delta}
(17.1)
Вывод для одного эквивалентного зазора, линейного сердечника и пренебрежимо малого рассеяния.

Порядок проверки результата

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

Пояснение программы

Магнитная постоянная задана в единицах СИ, площадь — в квадратных метрах, зазор — в метрах. В каждом режиме изменяется только зазор. Из индуктивности и того же тока вычисляется потокосцепление. Аналитическое выражение предполагает однородное поле и пренебрежение рассеянием и магнитным сопротивлением железа. При переходе к полевой модели необходимо учесть число физических зазоров и пространственное распределение поля.

Проверка масштабных законов в реальной геометрии

Для предварительного эталона сопротивление железа мало, а поле сосредоточено в одном эффективном зазоре. Реальная модель может содержать несколько зазоров.

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

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

  3. Рост расхождения с простой схемой объясняйте после анализа рассеяния потока, выпучивания поля и насыщения. Повторите контроль силы независимым методом.

Таблица 17.1. Обозначения и единицы
ВеличинаФизический смыслЕдиница СИ
NЧисло витков с учётом принятого направления1
IТок обмотки; в статике постоянныйА
gДлина одного воздушного зазорам
LИндуктивность при явно указанной линеаризацииГн
FКомпонента электромагнитной силыН

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

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

Вычисленный аналитический эталон

Дано. N = 200; S = 4 \cdot 10^{-4} \,\text{м}^2; \delta = 1 \,\text{мм}; I = 1 \,\text{А}; \mu_0 = 4\pi \cdot 10^{-7} \,\text{Гн}/\text{м}.

  1. L = 4\pi \cdot 10^{-7} \cdot 200^{2} \cdot 4 \cdot 10^{-4} / 10^{-3} = 0{,}020106 \,\text{Гн}.

    L=\frac{4\pi\cdot10^{-7}\cdot200^2\cdot4\cdot10^{-4}}{10^{-3}}=0{,}0201062\ \text{Гн}
    (17.2)
  2. B_\delta = 4\pi \cdot 10^{-7} \cdot 200 / 10^{-3} = 0{,}251327 \,\text{Тл}.

  3. \Phi = 0{,}000100531 \,\text{Вб}; \Psi = 200\Phi = 0{,}020106 \,\text{Вб}.

Ответ: Аналитический ориентир: L = 20{,}106 \,\text{мГн} и B_\delta = 0{,}25133 \,\text{Тл} при 1 \,\text{А}.

Проверка: Удвоение зазора в этой упрощённой модели уменьшает L и B вдвое; FEMM покажет отклонение из-за рассеяния.

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

Карты катушки и электромагнита по току и зазору

Таблица 17.2. Исходные данные практикума
Исходная величинаЗначениеЕдиница
Витки200
Площадь полюса400мм²
Зазоры0,5; 1; 2мм
Токи0; 0,5; 1; 2А
  1. Откройте модель электромагнита в графическом интерфейсе; проверьте, что сумма физических зазоров соответствует \delta в аналитическом ориентире.
  2. Сохраните рабочую копию, назначьте одному якорю отдельную группу и перемещайте только её.
  3. Для каждого зазора выполните проход по токам; извлекайте \Psi из цепи и B в заранее определённой средней точке зазора.
  4. Выполните низкотоковый линейный опыт, затем верните нелинейную B\text{–}H характеристику сердечника.
  5. Оцените длину меди и сопротивление при заданной температуре отдельным расчётом; не изменяйте NI этим сопротивлением при источнике тока.
  6. Для крайних зазоров повторите сеточную проверку и зафиксируйте геометрическое состояние.

Ожидаемый результат. Карты \Psi(I,\delta) и B_\delta(I,\delta), сопоставленные с предельной магнитной цепью; насыщение исследуется отдельно.

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

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

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

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

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

Итог работы. Карта результатов с проверенными предельными точками и объяснением насыщения.

Программа с русскими пояснениями: Численная проверка выведенного эталона · Порядок запуска

Численная проверка выведенного эталона

Код выполняется во встроенном Python FEMM 5. Он проверяет приведённый аналитический вывод; получение полевых данных выполняется отдельным практикумом с указанной моделью.

Python · встроенный MicroPython FEMM 5
"""Повтор катушки и электромагнита в реальной геометрии.

Среда: встроенный MicroPython FEMM 5.
Аналитический учебный пример; полевой решатель не вызывается.
Параметры, единицы величин и порядок вывода пояснены по стадиям.
"""

import math

# Магнитная постоянная выражена в Гн/м. Все геометрические величины ниже — в СИ.
mu0 = 4 * math.pi * 1e-7
# Число витков безразмерно; площадь в м²; ток в А.
turns, area_m2, current_a = 200, 4e-4, 1.0
# Меняем только длину эквивалентного воздушного зазора, м.
# Формулы предполагают однородное поле, пренебрежение рассеянием и сопротивлением железа.
for gap_m in [0.0005, 0.001, 0.002]:
    # Индуктивность рассчитываем из магнитного сопротивления воздушного зазора.
    inductance_h = mu0 * turns**2 * area_m2 / gap_m
    flux_density_t = mu0 * turns * current_a / gap_m
    linkage_wb = inductance_h * current_a
    # Порядок вывода: зазор, м; индуктивность, Гн; индукция, Тл; потокосцепление, Вб.
    print(gap_m, inductance_h, flux_density_t, linkage_wb)

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

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

Опорная теория: 8.7 · Измеряемые потокосцепление и энергия; 8.1 · Физическая постановка и свойства среды.

Применение: 38.1 · Связать программу с устройством аппарата; 40.1 · Повторить конструкцию П-магнита.

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

Пошаговая настройка расчётного примера

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

Все расчётные примеры курса

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

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

Таблица 17.3. Результаты расчётного образца
ВеличинаРезультатЕдиница
Частота0Гц
Активная длина20мм
Ток цепи Катушка5А
Потокосцепление цепи Катушка0,0153Вб
Сила по оси y-29,65Н
Энергия поля0,03824Дж
Модуль B в контрольной точке0,2922Тл
Координаты контрольной точки (x, y)20; 20,5мм
Визуализация численных данных FEMM 5: Электромагнит с отдельным якорем. Это экспорт результатов поля, а не снимок окна программы.
Рисунок 17.1. Визуализация численных данных FEMM 5: Электромагнит с отдельным якорем. Это экспорт результатов поля, а не снимок окна программы.

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

Ориентиры в интерфейсе
Геометрия: электромагнит с двумя зазорами и отдельным якорем
Рисунок 17.2. Геометрия: электромагнит с двумя зазорами и отдельным якорем Открыть в полном размере
Распределение магнитной индукции: электромагнит с двумя зазорами и отдельным якорем При открытии сохранённого решения FEMM сообщает Multiply Defined Regions. Кадр показывает интерфейс и картину поля; количественные результаты требуют проверки задания областей.
Рисунок 17.3. Распределение магнитной индукции: электромагнит с двумя зазорами и отдельным якорем При открытии сохранённого решения FEMM сообщает Multiply Defined Regions. Кадр показывает интерфейс и картину поля; количественные результаты требуют проверки задания областей. Открыть в полном размере
Встроенная консоль FEMM: Повтор катушки и электромагнита в реальной геометрии. Аналитическая проверка программы занятия. Показана консоль исходной редакции программы, зафиксированной при съёмке. Текущая программа с русскими пояснениями проверена отдельно во встроенном Python FEMM; этот снимок не является новым запуском текущей редакции.
Рисунок 17.4. Встроенная консоль FEMM: Повтор катушки и электромагнита в реальной геометрии. Аналитическая проверка программы занятия. Показана консоль исходной редакции программы, зафиксированной при съёмке. Текущая программа с русскими пояснениями проверена отдельно во встроенном Python FEMM; этот снимок не является новым запуском текущей редакции. Открыть в полном размере

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

Определите диапазон тока, где аналитический линейный ориентир и поле отличаются не более заранее выбранных 10 \,\%. Объясните физические причины отклонения.

Критерий выполнения: Есть карта, независимый ориентир и сеточная проверка; критерий 10 \,\% обозначен учебным и не превращён в отраслевую норму.
Первоисточники и дополнительное чтение по занятию
  1. [33] — 10–15, 39–52, 99–110
  2. [81] — Раздел Packages · Локальная копия
  3. [19] — Выделенный Bianchi_FEMM_Lua_section.pdf, с. 24–32; исходный PDF, с. 59–67
Пояснения и типичные ошибки
Почему зависимость L \propto 1/\delta перестаёт быть точной?

Изменяются выпучивание и рассеяние потока, а при насыщении существенна магнитодвижущая сила сердечника.

Можно ли сравнивать один зазор модели с суммой двух реальных зазоров?

Нет. В аналитической магнитной цепи \delta должна обозначать суммарную длину воздушных участков по пути потока.

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

  • Площадь в \text{мм}^2 необходимо перевести в \text{м}^2 перед подстановкой \mu_0.
  • Не использовать половину симметричной геометрии без корректного восстановления витков и энергии.