FEMM 5

Часть V · Электрические машины → § 29 · Идентификация асинхронной машины

Занятие 29.2

Проверка параметров на независимых частотах

Цель и задачи

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

Задачи

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

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

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

Расчётная лабораторная работа · Извлечение параметров ротора и восстановление рабочей характеристики

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Критерий выполнения. Наборы фиксированы заранее; комплексная ошибка показана; влияние сетки проверено; усложнение объяснено.

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

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

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

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

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

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

Проводимость рассеяния по коронкам зубцов

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

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

Получаемые величины. Φкоронок; Λкоронок и безразмерная λкоронок; отклонение от двух аналитических формул

  1. Выделяем поток зоны коронок из полного потока шлица при фиксированном пазовом токе.
  2. Вычитаем вклад условно равномерного участка шлица.
  3. Варьируем глубину шлица и отношение раскрытие/зазор.
  4. Сравниваем формулы и FEM на трёх сетках; объясняем область применимости формулы.

Страницы источника. PDF 1: определение; 2–3: формулы и численные примеры; 4: геометрическая чувствительность; 5: программа

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

Асинхронный двигатель: от схемы замещения к проверке поля и момента

Согласование схемы замещения и полевого расчёта позволяет проверить токи, магнитные оси, момент и баланс мощности асинхронного двигателя.

Лабораторные работы по методу: Асинхронный двигатель: механическая характеристика и вращающееся поле

Получаемые величины. I1, I2, фазовые углы и пазовые токи; B_{\delta}1, B зубцов, B ярма; T, P\delta, P мех; ошибки эквивалентного гладкого зазора

  1. Проверяем переход от Г- к Т-образной схеме и базовые единицы.
  2. Восстанавливаем мгновенные токи статора и ротора при выбранном скольжении.
  3. Согласуем магнитные оси отдельными расчётами; получаем нагрузочный режим.
  4. Сравниваем подробную геометрию с гладкой при \delta, k\delta \delta и k\delta k\mu \delta.
  5. Проверяем активную длину, заполнение, кольца, лобовые части и потери перед выводом о мощности на валу.

Страницы источника. PDF 1–10: схема; 11–18: токи и оси; 19–20: мощность; 20–26: упрощённые модели; 26–29: программа и входы; 30–34: Lua

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

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

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

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

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

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

Найденные R и L подставляют в Z_{\text{расч}} = R + j2 \pi fL на новой частоте. Сравнение с Z_{\text{опыт}} даёт невязку. Изменять параметры по той же проверочной точке и продолжать считать её независимой нельзя. Расхождение может быть связано с вытеснением тока, неверной поправкой на лобовые части, погрешностью сетки или недостаточной схемой замещения. Отдельный сеточный опыт на проверочной частоте помогает отличить численную погрешность от погрешности модели. Для проверки физической достоверности нужен внешний аналитический, экспериментальный или опубликованный результат. Собственная сохранённая таблица проверяет только повторяемость.

\varepsilon_Z(f)=\frac{|Z_{\text{расч}}(f)-Z_{\text{опыт}}(f)|}{\max(|Z_{\text{опыт}}(f)|,Z_{\text{масш}})}
(29.3)
Проверочная частота не используется для поиска параметров.

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

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

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

Переменная predicted содержит сопротивление, рассчитанное по ранее определённым параметрам. Переменная observed содержит независимый опорный результат на новой частоте. Модуль комплексной разности учитывает расхождение обеих составляющих сопротивления. Защитный масштаб ограничивает деление около нулевого результата. Параметры не изменяются по опорной точке, иначе проверка перестала бы быть независимой.

Проверка модели вне точек подгонки

Контрольные частоты выбираются до просмотра остатка и не участвуют в оценке параметров.

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

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

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

Таблица 29.4. Обозначения и единицы
ВеличинаФизический смыслЕдиница СИ
ZКомплексное электрическое сопротивлениеОм
fЧастота гармоникиГц
qПроверяемый результат; единица определяется выбором результатапо результату

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

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

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

Дано. По обучающим точкам получены R = 0{,}5 \,\text{Ом} и L = 0{,}02 \,\text{Гн}. Проверочная частота 2 \,\text{Гц}; Z_{\text{опыт}} = 0{,}51 + j0{,}25 \,\text{Ом}.

  1. Z_{\text{расч}} = 0{,}5 + j0{,}251327412 \,\text{Ом}.

    \varepsilon_Z=\frac{|(0{,}5+j0{,}251327)-(0{,}51+j0{,}25)|}{|0{,}51+j0{,}25|}=0{,}0177607
    (29.4)
  2. Разность -0{,}01 + j0{,}001327412 \,\text{Ом}; её модуль 0{,}0100877 \,\text{Ом}.

  3. Модуль Z_{\text{опыт}} = 0{,}567978873 \,\text{Ом}; относительная ошибка \approx 1{,}77607 \,\%.

Ответ: На независимой частоте ошибка около 1{,}776 \,\% при ненулевом Z_{\text{масш}} \le |Z_{\text{опыт}}|. Тогда знаменатель равен |Z_{\text{опыт}}|; больший защитный масштаб изменит приведённый процент.

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

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

Проверка схемы на отложенной частоте

Таблица 29.5. Исходные данные практикума
Исходная величинаЗначениеЕдиница
Обучающие частоты1; 3Гц
Проверочные частоты2; 4Гц
Порог параметровзаранее задан
  1. Загрузите найденные параметры предыдущего занятия, не меняя их.
  2. Рассчитайте предсказание на двух отложенных частотах и сохраните комплексные остатки.
  3. Повторите полевой опыт на уточнённой сетке для одной проверочной частоты.
  4. Если остаток устойчив, определите, какой механизм отсутствует в выбранной схеме; не добавляйте параметры без физического объяснения.
  5. Повторите подбор с исключением одной обучающей точки и сравните изменение параметров.
  6. В отчёте разделите ошибки подбора, предсказания, сетки и независимой проверки физической достоверности.

Ожидаемый результат. Модель проверена на частотах, не использованных при определении параметров; устойчивость параметров исследована повторным отбором исходных данных.

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

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

Предсказывает ли схема независимую частотную точку?

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

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

Итог работы. Ошибки действительной и мнимой частей и вывод о переносимости параметров.

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

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

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

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

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

import math

# Полное сопротивление на новой частоте 2 Гц вычисляем по ранее определённым R и L.
predicted = 0.5 + 1j * 2 * math.pi * 2 * 0.02
# Независимая опорная точка, Ом; её не используют для повторного определения R и L.
observed = 0.51 + 0.25j
# Модуль комплексной разности учитывает расхождение активной и реактивной частей.
# Защитный масштаб 0,01 Ом предотвращает неустойчивое деление около нуля.
relative = abs(predicted - observed) / max(abs(observed), 0.01)
print("Предсказание, Ом:", predicted)
print("Проверочная ошибка, %:", 100 * relative)

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

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

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

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

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

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

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

Таблица 29.6. Результаты расчётного образца
ВеличинаРезультатЕдиница
Частота1Гц
Активная длина100мм
Ток цепи Фаза А1 + j0А
Потокосцепление цепи Фаза А0{,}0421 - j0{,}03968Вб
Ток цепи Фаза С-0{,}5 - j0{,}866А
Потокосцепление цепи Фаза С-0{,}05517 - j0{,}01806Вб
Ток цепи Фаза В-0{,}5 + j0{,}866А
Потокосцепление цепи Фаза В0{,}01251 + j0{,}05637Вб
Средняя энергия поля0,03197Дж
Норма комплексных амплитуд B в контрольной точке0,4584Тл
Координаты контрольной точки (x, y)50; 0мм
Визуализация численных данных FEMM 5: Асинхронный двигатель: сектор модели. Это экспорт результатов поля, а не снимок окна программы.
Рисунок 29.5. Визуализация численных данных FEMM 5: Асинхронный двигатель: сектор модели. Это экспорт результатов поля, а не снимок окна программы.

Учебная модель по Meeker 2horse. Частота — частота тока ротора в данном неподвижном опыте. Геометрия представляет четверть машины; результаты приведены без дополнительного множителя 4, который выводят отдельно. Пуск рассматривается в отдельном расчётном опыте. Комплексные токи и потокосцепления заданы амплитудными значениями; они отличаются от действующих значений источника.

Ориентиры в интерфейсе
Геометрия: адаптированная модель Meeker: 36 пазов статора, 28 стержней, неподвижный AC-опыт
Рисунок 29.6. Геометрия: адаптированная модель Meeker: 36 пазов статора, 28 стержней, неподвижный AC-опыт Открыть в полном размере
Распределение магнитной индукции: адаптированная модель Meeker: 36 пазов статора, 28 стержней, неподвижный AC-опыт
Рисунок 29.7. Распределение магнитной индукции: адаптированная модель Meeker: 36 пазов статора, 28 стержней, неподвижный AC-опыт Открыть в полном размере
Встроенная консоль FEMM: Проверка параметров на независимых частотах. Аналитическая проверка программы занятия. Показана консоль исходной редакции программы, зафиксированной при съёмке. Текущая программа с русскими пояснениями проверена отдельно во встроенном Python FEMM; этот снимок не является новым запуском текущей редакции.
Рисунок 29.8. Встроенная консоль FEMM: Проверка параметров на независимых частотах. Аналитическая проверка программы занятия. Показана консоль исходной редакции программы, зафиксированной при съёмке. Текущая программа с русскими пояснениями проверена отдельно во встроенном Python FEMM; этот снимок не является новым запуском текущей редакции. Открыть в полном размере

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

Заранее отложите две частоты и сравните постоянную RL-схему со схемой, содержащей одну физически обоснованную частотную поправку.

Критерий выполнения: Наборы фиксированы заранее; комплексная ошибка показана; влияние сетки проверено; усложнение объяснено.
Первоисточники и дополнительное чтение по занятию
  1. [33] — 10–15, 39–52, 110
  2. [81] — Раздел Packages · Локальная копия
  3. [4] — Скачанный RapidIM_Report.pdf, с. 1–16
  4. [37]
Пояснения и типичные ошибки
Можно ли после неудачной проверки подогнать к ней модель?

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

Достаточен ли повтор собственного CSV?

Он проверяет воспроизводимость, но не физическую достоверность модели.

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

  • Малое число точек не позволяет надёжно идентифицировать большое число параметров.
  • Модуль импеданса без фазы может скрыть компенсирующие ошибки R и X.