FEMM 5

Часть II · Теория четырёх физических модулей FEMM → § 8 · Магнитное поле, цепи и энергетические величины

Занятие 8.9

Знаки \(M\), взаимность и энергетическая проверка

Цель и задачи

Цель. Проверить знаки взаимной индуктивности через энергию.

Задачи

  • Вывести W_{\text{м}}=I^{\mathsf T}LI/2.
  • Проверить M^{2}\le L_{1} L_{2}.
Что нужно вспомнить перед началом

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

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

Расчётная лабораторная работа · Определение потокосцеплений, индуктивности и запаса энергии

Инженерный вопрос. Согласованы ли взаимность, направления обмоток и энергия?

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

Лабораторная работа:

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

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

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

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

  2. Допустимая матрица и встречные токи. Проверьте переход от тока и потокосцепления к параметру цепи и энергии в том же рабочем состоянии.

  3. Энергетический аудит двух обмоток. Выполните действия практикума по его исходной таблице. Запишите оба взаимных коэффициента и энергию при согласном и встречном токах.

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

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

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

Данные именно этого опыта. Запишите оба взаимных коэффициента и энергию при согласном и встречном токах.

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

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

Проверка из практикума. Взаимность в пределах 1\,\%; энергетическое сравнение с указанием абсолютного допуска для малой энергии.

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

Критерий выполнения. 6 \,\text{мГн} допустимо: M^{2}=36<40 \,\text{мГн}^2;7 \,\text{мГн} недопустимо:49>40. Обоснование выполнено через энергию.

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

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

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

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

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

Дополнительный расчётный метод

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

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

Лабораторные работы по методу: Зубцовые контуры: проводимости, коронки и взаимная связь при смещении

Получаемые величины. \Lambda =\Phi /F; \Lambda 12(x); влияние распределения тока; период зависимости связи

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

Страницы источника. PDF 1–2: определение; 3–4: программа и движение; 5–6: сравнение токовых моделей

Литература по методу

Энергия магнитного поля через векторный потенциал

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

Лабораторные работы по методу: Трансформатор от потокосцеплений до схемы замещения и нагрузочной характеристики; Асинхронный двигатель от поля в зазоре до параметров и механической характеристики; Генератор с постоянными магнитами от конструкции ротора до ЭДС и нагрузочного напряжения; Определение коэффициентов поля явнополюсной машины; Синхронный генератор: от полевых параметров к внешней и угловой характеристикам; Зубцовые контуры: проводимости, коронки и взаимная связь при смещении; Контрольный опыт: от условия задачи к инженерному выводу; Опыт Эрстеда: момент стрелки, эквивалентный ток и движение; Измерительная рамка: от поля провода к ЭДС; Аналогии полей: проводимость, дуальная сеть и диффузия; От магнитной цепи к карте \Psi (i,\theta ) и моменту SRM

Получаемые величины. AJ=\int A\cdot J dV; W и W'; матрица индуктивностей и перекрёстное насыщение; устойчивость параметров при сетке и перемещении

  1. Для линейной токовой системы сопоставить AJ, \Sigma i\Psi и 2W.
  2. Для нелинейной системы интегрировать H по B и B по H, сохраняя одну рабочую точку.
  3. Определить секущие и локальные параметры отдельными опытами.
  4. Проверить сетку в узких зазорах и перенос карты в уравнения цепи/механики.

Страницы источника. PDF 1–7: энергия и функционал; 8–9: нелинейная магнитная цепь; 10–18: постановка и параметры; 19–25: сетка и движение; 26–30: связанные модели.

Литература по методу

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

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

Соглашения об обозначениях и единицах
  • \(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. (8.79) · Вывод квадратичной энергии
  2. (8.80) · Откуда следует взаимность
  3. (8.81) · Неотрицательность и коэффициент связи
  4. (8.82) · Взаимность и положительность энергии
  5. (8.83) · Численный пример
  6. (8.84) · Численный пример

Начала обмоток определяют знак

Положительный ток каждой цепи имеет выбранную ориентацию. Если оба создают общий поток в одном направлении, взаимная индуктивность положительна; если в противоположных — отрицательна. При смене начала второй обмотки меняются знаки второй \Psi и второй координаты тока. Диагональные коэффициенты остаются положительными, а обе взаимные связи меняют знак. Отрицательный M не означает отрицательную энергию или ошибку: он выражает выбранное согласование ориентаций.

Вывод квадратичной энергии

В линейной обратимой системе без постоянных источников \Psi=LI. При медленном нарастании вектора токов от нуля по лучу I(s)=s\mathbf I электрическая работа равна интегралу I(s)^{\mathsf T}\mathrm d\Psi(s). Подстановка \Psi(s)=s\mathbf L\mathbf I даёт интеграл s от нуля до единицы и коэффициент 1/2. Общая энергия содержит перекрёстный член MI_1 I_{2}. Его знак может уменьшать энергию встречного возбуждения, но полная энергия остаётся неотрицательной.

W_\text{м}=\frac12\mathbf I^\mathsf T\mathbf L\mathbf I=\frac12L_{11}I_1^2+MI_1I_2+\frac12L_{22}I_2^2
(8.79)
Формула предполагает вещественную симметричную матрицу статической обратимой системы без постоянного смещения.

Откуда следует взаимность

Если \Psi_j=\partial W^{\prime}/\partial I_j, то равенство смешанных вторых производных коэнергии даёт \partial \Psi_{1}/\partial I_{2}=\partial \Psi_{2}/\partial I_{1}. В линейном случае это L_{12}=L_{21}. Материальный закон должен быть обратимым и гладким в рассматриваемом диапазоне. При численном извлечении допускается малая разность, которую проверяют на более мелкой сетке. Сильная разность обычно требует проверки источников, масштабов и ориентаций. Утверждение о взаимности не переносится автоматически на любое активное устройство либо произвольно определённый комплексный параметр.

L_{12}=\frac{\partial^2W^{\prime}_\text{м}}{\partial I_1\partial I_2}=\frac{\partial^2W^{\prime}_\text{м}}{\partial I_2\partial I_1}=L_{21}
(8.80)
Условие — существование скалярной гладкой коэнергии в выбранном обратимом состоянии.

Неотрицательность и коэффициент связи

Квадратичная форма должна быть неотрицательна при любых токах. Для двух обмоток это требует L_{11}\ge 0,L_{22}\ge 0 и неотрицательного определителя. Отсюда |M|\le \sqrt{L_{11}L_{22}}. При L_{11}>0 и L_{22}>0 нормированный коэффициент k=M/\sqrt{L_{11}L_{22}} находится между -1 и 1. Граница |k|=1 соответствует идеальному общему потоку без независимого рассеяния. Матрица тогда вырождена: определённое встречное возбуждение не создаёт поля. При нулевой собственной индуктивности соответствующая взаимная связь равна нулю, но коэффициент k не определён. В реальной системе рассеяние обычно даёт положительный остаток энергии.

\det\mathbf L=L_{11}L_{22}-M^2\ge0,\qquad k=\frac{M}{\sqrt{L_{11}L_{22}}},\quad |k|\le1
(8.81)
Определитель имеет единицу \text{Гн}^2, коэффициент k безразмерен. Проверка относится к энергетической статической матрице.

Как сравнить энергию FEMM и матрицы

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

Физический смысл встречного возбуждения

При двух равных магнитодвижущих силах противоположного направления основной поток уменьшается, но локальные поля рассеяния могут сохраняться. Это объясняет, почему реальный трансформатор при коротком замыкании не имеет строго нулевой энергии. Для первого опыта полезно подобрать токи около N_1I_1+N_2I_2=0, затем снять поля в окне и вне сердечника. Интерактивный элемент показывает общую энергию при смене знака M и токов; недопустимое |k|>1 должно отмечаться как нарушение энергетического условия.

Неравенство связи можно вывести завершением квадрата

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

Взаимность и положительность энергии

Рассматривается пассивная обратимая линейная магнитная система при неподвижной геометрии. Направления токов и потокосцеплений согласованы.

  1. Симметрия смешанных производных коэнергии даёт взаимность двух межобмоточных коэффициентов. Она проверяется независимо, а не навязывается усреднением до проверки.

  2. При L_{11}>0 запишите энергию как сумму полного квадрата и остатка. Неотрицательность при любых токах ограничивает взаимную индуктивность. В вырожденном случае L_{11}=0 пассивность требует M=0, и деление не применяется.

    2W=L_{11}\left(I_1+\frac{M}{L_{11}}I_2\right)^2+\left(L_{22}-\frac{M^2}{L_{11}}\right)I_2^2\ge0
    (8.82)
  3. Разворот выводов одной обмотки меняет знак её потокосцепления, тока и взаимных коэффициентов. Физическая энергия того же состояния при согласованной замене координат остаётся неизменной.

Таблица 8.22. Обозначения и единицы
ВеличинаФизический смыслЕдиница СИ
LИндуктивность при явно указанной линеаризацииГн
MВзаимная индуктивностьГн
IТок обмотки; в статике постоянныйА
WМагнитная энергияДж

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

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

Допустимая матрица и встречные токи

Дано. L_{11}=12 \,\text{мГн}; L_{22}=3 \,\text{мГн}; M=5 \,\text{мГн}; I_{1}=1 \,\text{А}; I_{2}=\pm 2 \,\text{А}.

  1. Проверка коэффициента связи и определителя.

    k=\frac5{\sqrt{12\cdot3}}=0{,}8333,\qquad\det\mathbf L=11\cdot10^{-6}\,\text{Гн}^2
    (8.83)
  2. Энергии для двух знаков второго тока.

    W_\text{м}(1,2)=0{,}022\,\text{Дж},\qquad W_\text{м}(1,-2)=0{,}002\,\text{Дж}
    (8.84)

Ответ: Матрица допустима; энергии 22 мДж и 2 мДж положительны.

Проверка: Встречное возбуждение уменьшает энергию, но не обнуляет её, поскольку k<1.

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

Энергетический аудит двух обмоток

Таблица 8.23. Исходные данные практикума
Исходная величинаЗначениеЕдиница
Тип задачиПлоскопараллельная (Planar)
Глубина модели20мм
Наружный квадрат магнитопровода60 \times 60мм
Окно магнитопровода40 \times 40мм
Относительная проницаемость линейного сердечника1000
Витки первой / второй обмотки100 / 50
Частота0Гц
Полуширина внешнего квадрата150мм
  1. Откройте точную сохранённую копию двухобмоточного магнитопровода 8.8 и создайте отдельную рабочую копию. Сохраните геометрию 60×60 мм, окно 40×40 мм, Depth=20 мм, \mu _{r}=1000, витки 100/50, направления цепей и внешнюю границу; заново строить тот же образец не требуется.
  2. Получите матрицу двумя независимыми возбуждениями и проверьте L_{12}\approx L_{21}, затем k и определитель.
  3. Рассчитайте энергию по матрице для токов (1;2) А и (1;-2) А; выполните два соответствующих решения и независимый интеграл полной энергии.
  4. Смените начала второй обмотки и подтвердите смену знака M при неизменных L_{11},L_{22}.

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

Как проверить. Взаимность в пределах 1\,\%; энергетическое сравнение с указанием абсолютного допуска для малой энергии.

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

Согласованы ли взаимность, направления обмоток и энергия?

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

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

Итог работы. Энергетический аудит с объяснением знака взаимной связи.

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

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

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

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

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

Таблица 8.24. Результаты расчётного образца
ВеличинаРезультатЕдиница
Частота0Гц
Активная длина20мм
Ток цепи Первичная обмотка1А
Потокосцепление цепи Первичная обмотка0,01393Вб
Ток цепи Вторичная обмотка0А
Потокосцепление цепи Вторичная обмотка0,00683Вб
Энергия поля0,006967Дж
Модуль B в контрольной точке0,6848Тл
Координаты контрольной точки (x, y)0; 25мм
Визуализация численных данных FEMM 5: Двухобмоточный магнитопровод: первичное возбуждение. Это экспорт результатов поля, а не снимок окна программы.
Рисунок 8.34. Визуализация численных данных FEMM 5: Двухобмоточный магнитопровод: первичное возбуждение. Это экспорт результатов поля, а не снимок окна программы.

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

Ориентиры в интерфейсе
Геометрия: две воздушные осесимметричные обмотки; вторичная разомкнута
Рисунок 8.35. Геометрия: две воздушные осесимметричные обмотки; вторичная разомкнута Открыть в полном размере
Распределение магнитной индукции: две воздушные осесимметричные обмотки; вторичная разомкнута
Рисунок 8.36. Распределение магнитной индукции: две воздушные осесимметричные обмотки; вторичная разомкнута Открыть в полном размере
Геометрия: двухобмоточный линейный сердечник, независимое первичное возбуждение
Рисунок 8.37. Геометрия: двухобмоточный линейный сердечник, независимое первичное возбуждение Открыть в полном размере
Распределение магнитной индукции: двухобмоточный линейный сердечник, независимое первичное возбуждение
Рисунок 8.38. Распределение магнитной индукции: двухобмоточный линейный сердечник, независимое первичное возбуждение Открыть в полном размере
Расчёт при возбуждении только вторичной обмотки двухобмоточного магнитопровода.
Рисунок 8.39. Расчёт при возбуждении только вторичной обмотки двухобмоточного магнитопровода. Открыть в полном размере
Ориентированный контур поперёк стержня и интеграл нормальной индукции.
Рисунок 8.40. Ориентированный контур поперёк стержня и интеграл нормальной индукции. Открыть в полном размере

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

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

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

Окна FEMM 5

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

Две обмотки: матрица индуктивностей, взаимность и энергия. Геометрия учебной модели; подписи скрыты для читаемости.
Две обмотки: матрица индуктивностей, взаимность и энергия. Геометрия учебной модели; подписи скрыты для читаемости. Открыть в полном размере
Двухобмоточный трансформатор: первое статическое базисное возбуждение для матрицы индуктивностей. Карта |B|, Т, и линии поля; рабочая область около сердечника и обмоток.
Двухобмоточный трансформатор: первое статическое базисное возбуждение для матрицы индуктивностей. Карта |B|, Т, и линии поля; рабочая область около сердечника и обмоток. Открыть в полном размере

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

Для L_{1}=10 \,\text{мГн},L_{2}=4 \,\text{мГн} проверьте варианты M=6 и 7 \,\text{мГн}.

Критерий выполнения: 6 \,\text{мГн} допустимо: M^{2}=36<40 \,\text{мГн}^2;7 \,\text{мГн} недопустимо:49>40. Обоснование выполнено через энергию.
Первоисточники и дополнительное чтение по занятию
  1. [33] — 46–48, 52
  2. [19] — 31–32 скачанного раздела Bianchi_FEMM_Lua_section.pdf; номер страницы оригинала на 35 больше
Пояснения и типичные ошибки
Является ли M<0 ошибкой?

Нет, это может быть правильная ориентация начал.

Почему |k|>1 подозрительно?

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

Почему \Psi взаимны?

Они являются производными одной коэнергии при принятых условиях обратимости.

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

  • Объявить отрицательный M физически невозможным.
  • Сравнивать энергии в разных Depth.
  • Проверять взаимность по преждевременно округлённым числам.