FEMM 5

Часть V · Электрические машины → § 32 · Коммутация и динамика вентильно-индукторной машины

Занятие 32.3

Динамика цепи и механики с квазистатическим восстановлением поля

Цель и задачи

Цель. Проверить связанную цепь и механику по балансу энергии.

Задачи

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

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

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

Расчётная лабораторная работа · Восстановление переходного режима по полевым данным

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

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

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

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

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

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

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

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

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

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

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

Результат этого занятия. Траектория с оценкой шага и невязки энергии в пределах карты.

Данные именно этого опыта. Сохраните ток, положение, скорость, момент и составляющие энергетического баланса.

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

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

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

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

Критерий выполнения. Входы и начальные условия полны; события и область карты соблюдены; энергетика и шаги проверены.

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

Основание сравнения. Опора — уравнения цепи, движения и аналитический предельный случай из занятия. Временная модель строится поверх квазистатических полевых данных; её допущения должны сохраняться на всей траектории.

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

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

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

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

Моделирование опытов Эрстеда и Фарадея

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

Лабораторные работы по методу: Опыт Эрстеда: момент стрелки, эквивалентный ток и движение; Измерительная рамка: от поля провода к ЭДС

Получаемые величины. момент стрелки M(\theta ); ошибка дискретизации токового слоя; эквивалентные \mu и намагниченность; угол и энергия движения при заданном Jмех

  1. Зафиксировать геометрию провода, стрелки и источник.
  2. Вычислить поле внешнего тока и силу каждого элемента эквивалентного токового слоя.
  3. Суммировать моменты и уточнить число элементов/толщину слоя при сохранении полных токов.
  4. Сравнить с FEMM при совместимых свойствах материала.
  5. Для динамики добавить момент инерции и оговорённое трение, проверить шаг времени и баланс энергии.

Страницы источника. PDF 1–19 полностью: 1–4 опыты и амперовская модель; 5–14 алгоритм сил и программные иллюстрации; 15–17 сравнение FEMM; 18–19 движение стрелки.

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

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

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

Лабораторные работы по методу: Трансформатор от потокосцеплений до схемы замещения и нагрузочной характеристики; Асинхронный двигатель от поля в зазоре до параметров и механической характеристики; Генератор с постоянными магнитами от конструкции ротора до ЭДС и нагрузочного напряжения; Определение коэффициентов поля явнополюсной машины; Синхронный генератор: от полевых параметров к внешней и угловой характеристикам; Зубцовые контуры: проводимости, коронки и взаимная связь при смещении; Контрольный опыт: от условия задачи к инженерному выводу; Опыт Эрстеда: момент стрелки, эквивалентный ток и движение; Измерительная рамка: от поля провода к ЭДС; Аналогии полей: проводимость, дуальная сеть и диффузия; От магнитной цепи к карте \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. (32.5) · Физическая и вычислительная постановка
  2. (32.6) · Связанная электрическая и механическая энергия
  3. (32.7) · Численный пример

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

Механическая динамика выводится из второго закона Ньютона для вращения: J \mathrm d\omega/\mathrm dt = T_{\text{эм}} - T_{\text{наг}} - T_{\text{пот}}. Положение удовлетворяет \mathrm d\theta/\mathrm dt = \omega. Момент читают из карты при текущих токах и угле. Одновременно электрическое состояние меняется по u - Ri. Совместная система позволяет исследовать пуск и коммутацию в рамках квазистатических карт, но не учитывает автоматически вихревые токи, магнитный гистерезис и демпферную память.

Для контроля энергии умножают механическое уравнение на \omega и интегрируют: изменение J \omega^{2}/2 должно соответствовать работе разности моментов. Электрическая энергия источника распределяется между потерями меди, изменением магнитной энергии и механической работой. Карта коэнергии полезна для согласованного момента; интерполяция момента независимо от \Psi может нарушить энергетический баланс. При изменении скорости частота периодического режима меняется, поэтому предположение медленного изменения поля и материальных характеристик также необходимо оценить.

J\frac{d\omega_{\text{мех}}}{dt}=T_{\text{эм}}-T_{\text{наг}}-T_{\text{пот}},\quad\frac{d\theta_{\text{мех}}}{dt}=\omega_{\text{мех}},\quad W_{\text{кин}}=\frac{J\omega_{\text{мех}}^2}{2}
(32.5)
J — момент инерции, \text{кг}\cdot \text{м}^2; \omega в \text{рад}/\text{с}; момент в \text{Н}\cdot \text{м}.

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

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

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

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

Связанная электрическая и механическая энергия

Механическое состояние задаётся углом и скоростью; инерция относится ко всем вращающимся частям, включённым в модель.

  1. По электрическому состоянию и углу восстановите ток, поле и момент. Разность электромагнитного, нагрузочного и тормозящего моментов задаёт ускорение.

  2. Умножение уравнения механики на скорость даёт скорость изменения кинетической энергии.

    \frac{\mathrm d}{\mathrm dt}\frac{J\Omega^2}{2}=(T_{\text{эм}}-T_{\text{наг}}-T_{\text{пот}})\Omega
    (32.6)
  3. За интервал интегрирования сравните электрическую входную энергию с изменением магнитной и кинетической энергий, механической работой нагрузки и потерями. Затем уточните временной шаг.

Таблица 32.7. Обозначения и единицы
ВеличинаФизический смыслЕдиница СИ
\OmegaМеханическая угловая скоростьрад/с
\thetaМеханический угол при дифференцировании энергиирад
TЭлектромагнитный моментН·м
WМагнитная энергияДж
\Delta tШаг интегрирования по временис
JМомент инерции вращающихся частейкг·м²

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

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

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

Дано. J = 0{,}001 \,\text{кг}\cdot \text{м}^2; T_{\text{эм}} = 1 \,\text{Н}\cdot \text{м}; T_{\text{наг}} = 0{,}2 \,\text{Н}\cdot \text{м}; T_{\text{пот}} = 0; \omega(0) = 100 \,\text{рад}/\text{с}; \Delta t = 0{,}1 мс.

  1. Ускорение \alpha = (1 - 0{,}2)/0{,}001 = 800 \,\text{рад}/\text{с}^2; новая скорость 100{,}08 \,\text{рад}/\text{с}.

    \Delta W_{\text{кин}}=\frac{0{,}001}{2}(100{,}08^2-100^2)=0{,}0080032\ \text{Дж}
    (32.7)
  2. При постоянном ускорении \Delta\theta = \omega_{0} \Delta t + \alpha\Delta t^2/2 = 0{,}010004 \,\text{рад}.

  3. \Delta W_{\text{кин}} = 0{,}001(100{,}08^{2} - 100^{2})/2 = 0{,}0080032 \,\text{Дж}; работа 0{,}8 \cdot 0{,}010004 даёт то же число.

Ответ: За 0,1 мс скорость увеличилась на 0{,}08 \,\text{рад}/\text{с}, угол — на 0{,}010004 \,\text{рад}, кинетическая энергия — на 8,0032 мДж.

Проверка: Работа моментов совпадает с приростом энергии; этот аналитический тест независим от карты вентильно-индукторной машины.

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

Цепь, механика и контроль энергии

Таблица 32.8. Исходные данные практикума
Исходная величинаЗначениеЕдиница
Момент инерции0,001кг·м²
Начальная скорость100рад/с
Шаг0,1мс
  1. Проверьте механический интегратор по постоянным моментам указанного эталона.
  2. Подключите фазную \Psi-карту и карту коэнергии, задайте нагрузочную характеристику.
  3. Интегрируйте электрические и механические переменные согласованным шагом, разрешая угловые события.
  4. В выбранных состояниях восстановите поле отдельным статическим решением FEMM и сравните с интерполяцией.
  5. Рассчитайте энергию источника, меди, магнитного поля и механической работы по одной границе.
  6. Повторите с уменьшенным шагом времени и карты; укажите диапазон, в котором квазистатическое приближение принято.

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

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

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

Объясняется ли изменение скорости источником, потерями и полевой энергией?

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

С чем сравнить. Опора — уравнения цепи, движения и аналитический предельный случай из занятия. Временная модель строится поверх квазистатических полевых данных; её допущения должны сохраняться на всей траектории.

Итог работы. Траектория с оценкой шага и невязки энергии в пределах карты.

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

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

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

Python · встроенный MicroPython FEMM 5
"""Динамика цепи и механики с квазистатическим восстановлением поля.

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

# Момент инерции, кг·м²; электромагнитный момент и момент нагрузки, Н·м.
inertia_kgm2, torque_nm, load_nm = 0.001, 1.0, 0.2
# Начальная угловая скорость, рад/с; шаг времени, с. Моменты в пределах шага постоянны.
omega_start, step_s = 100.0, 0.0001
# Угловое ускорение определяется разностью моментов, делённой на момент инерции, рад/с².
acceleration = (torque_nm - load_nm) / inertia_kgm2
omega_end = omega_start + acceleration * step_s
# Изменение угла учитывает начальную скорость и постоянное ускорение, рад.
angle_change = omega_start * step_s + 0.5 * acceleration * step_s**2
# Независимо вычисляем приращение кинетической энергии и работу разности моментов, Дж.
energy_change = 0.5 * inertia_kgm2 * (omega_end**2 - omega_start**2)
mechanical_work = (torque_nm - load_nm) * angle_change
# Порядок вывода: конечная скорость; изменение угла; изменение энергии; механическая работа.
print(omega_end, angle_change, energy_change, mechanical_work)
print("Остаток энергии, Дж:", energy_change - mechanical_work)

Расчётное исследование и разбор результатов

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

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

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

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

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

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

Таблица 32.9. Результаты расчётного образца
ВеличинаРезультатЕдиница
Частота0Гц
Активная длина80мм
Ток цепи Фаза А3А
Потокосцепление цепи Фаза А0,06916Вб
Ток цепи Фаза С0А
Потокосцепление цепи Фаза С317,1 × 10⁻⁶Вб
Ток цепи Фаза В0А
Потокосцепление цепи Фаза В-0,001422Вб
Момент на исходной учебной сетке0,4394Н·м
Энергия поля0,1037Дж
Модуль B в контрольной точке0,03358Тл
Координаты контрольной точки (x, y)59,27; 5,186мм
Момент на проверенной уточнённой сетке0,439Н·м
Визуализация численных данных FEMM 5: Вентильно-индукторная машина: смещённый ротор. Это экспорт результатов поля, а не снимок окна программы.
Рисунок 32.15. Визуализация численных данных FEMM 5: Вентильно-индукторная машина: смещённый ротор. Это экспорт результатов поля, а не снимок окна программы.

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

Ориентиры в интерфейсе
Геометрия: вентильно-индукторная машина 6/4, первая фаза, согласованное положение
Рисунок 32.16. Геометрия: вентильно-индукторная машина 6/4, первая фаза, согласованное положение Открыть в полном размере
Распределение магнитной индукции: вентильно-индукторная машина 6/4, первая фаза, согласованное положение
Рисунок 32.17. Распределение магнитной индукции: вентильно-индукторная машина 6/4, первая фаза, согласованное положение Открыть в полном размере
Вентильно-индукторная машина: отдельное смещённое положение ротора; угол и знак проверяются по исходному файлу, ток фазы сохранён.
Рисунок 32.18. Вентильно-индукторная машина: отдельное смещённое положение ротора; угол и знак проверяются по исходному файлу, ток фазы сохранён. Открыть в полном размере
Вентильно-индукторная машина: смещённый ротор: выбранная сетка воздушного зазора.
Рисунок 32.19. Вентильно-индукторная машина: смещённый ротор: выбранная сетка воздушного зазора. Открыть в полном размере
Вентильно-индукторная машина: смещённый ротор: момент ротора на выбранной сетке.
Рисунок 32.20. Вентильно-индукторная машина: смещённый ротор: момент ротора на выбранной сетке. Открыть в полном размере
Встроенная консоль FEMM: Динамика цепи и механики с квазистатическим восстановлением поля. Аналитическая проверка программы занятия. Показана консоль исходной редакции программы, зафиксированной при съёмке. Текущая программа с русскими пояснениями проверена отдельно во встроенном Python FEMM; этот снимок не является новым запуском текущей редакции.
Рисунок 32.21. Встроенная консоль FEMM: Динамика цепи и механики с квазистатическим восстановлением поля. Аналитическая проверка программы занятия. Показана консоль исходной редакции программы, зафиксированной при съёмке. Текущая программа с русскими пояснениями проверена отдельно во встроенном Python FEMM; этот снимок не является новым запуском текущей редакции. Открыть в полном размере

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

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

Критерий выполнения: Входы и начальные условия полны; события и область карты соблюдены; энергетика и шаги проверены.
Первоисточники и дополнительное чтение по занятию
  1. [33] — 10–15, 49–52, 90–94, 105–111
  2. [81] — Раздел Packages · Локальная копия
Пояснения и типичные ошибки
Почему поля восстанавливают выборочно?

Карта ускоряет динамику, а отдельные решения проверяют интерполяцию; полный временной полевой решатель этим не создаётся.

Что происходит при выходе тока за карту?

Расчёт должен остановиться или карта должна быть обоснованно расширена; молчаливая экстраполяция недопустима.

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

  • Угловой период карты не отменяет изменение абсолютной механической позиции.
  • Динамика по статическим картам не воспроизводит автоматически потери от вихретоковой памяти.