FEMM 5

Часть V · Электрические машины → § 23 · Асинхронные машины

Занятие 23.2

Частота токов ротора и серия неподвижных гармонических расчётов

Цель и задачи

Цель. Интерпретировать серию неподвижных частотных опытов асинхронной машины.

Задачи

  • Получить комплексное отношение потокосцепления к току.
  • Связать его вещественную и мнимую части с роторной постоянной времени.
Что нужно вспомнить перед началом

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

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

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

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

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

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

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

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

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

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

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

  3. Четыре ручных частотных опыта. Выполните действия практикума по его исходной таблице. Запишите частоту, токи и комплексные величины четырёх опытов.

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

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

Результат этого занятия. Частотная таблица с объяснением связи с принятой моделью скольжения.

Данные именно этого опыта. Запишите частоту, токи и комплексные величины четырёх опытов.

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

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

Проверка из практикума. Проверьте один токовый масштаб в линейной модели и три уровня сетки. Явно укажите, что полученные коэффициенты относятся к выбранной амплитуде и модели ротора.

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

Критерий выполнения. Указаны частота, ток, комплексное потокосцепление, нормировка сектора и остаточная ошибка независимой точки; серия не названа временным расчётом пуска.

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

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

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

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

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

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

Задание 4. Основная ЭДС и намагничивающий параметр асинхронного двигателя

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

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

Получаемые величины. Bp, Eф 1, Eline1; AJ/(3 I_{\mathrm{RMS}}^{2}); f2=s f1; после отдельной AC-серии R'2(f2), L_{\sigma}(f2), T(s)

Страницы источника. PDF1–11; протокол 1; геометрия 2; токи 3–5; поле/FFT7–9; AJ10

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

Магнитные поля в электрических машинах — введение (файл Гончарова)

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

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

Получаемые величины. \Phi =\Delta Az·depth, B=\Delta Az/\Delta n; L, M; Bp и ошибка обратного восстановления; kf, k_{\Phi}, kd, kq; E0(If), Xd, Xq, U, I,\varphi,\delta; Jротора(f2), T(s), отличие дефектного ротора

Страницы источника. PDF1–75; A_z/поле 1–16; L_{17}–20; FFT20–29; синхронные 30–60; AD60–74

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

Hrabovcová, Rafajdus, Makyš. Analysis of Electrical Machines (2020)

Расчёт сопротивлений и индуктивностей по геометрии, уравнения электрических машин в фазных и d/qкоординатах, энергетический метод силы и момента, полевые опыты трансформатора, АД, синхронных и реактивных преобразователей. Параметры поля связываются с напряжением, мощностью и динамикой при явных допущениях модели.

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

Получаемые величины. F(g), L(I, g); параметры трансформатора; R'2, L_{\sigma}, T(s); E_{0}, Xd, Xq; ΨPM, Tcog; \Psi (I,\theta ), W'(I,\theta ), T

Страницы источника. PDF21–38: параметры по геометрии; PDF39–120: общая теория, d/q и динамические уравнения; PDF130–136: электромагнит; PDF137–146: трансформатор; PDF147–163: асинхронная машина; PDF164–183: синхронные/PM; PDF184–189: SRM

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

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

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

Соглашения об обозначениях и единицах
  • \(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. (23.4) · Серия в неподвижной геометрии
  2. (23.5) · Вывод формы частотного отклика
  3. (23.6) · Почему комплексная индуктивность зависит от частоты

Серия в неподвижной геометрии

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

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

\mathcal L(\omega_2)=\frac{\widehat\Psi}{\widehat I}=L_{\text{расс}}+\frac{L_{\text{м}}}{1+j\omega_2\tau_{\text{р}}}
(23.4)
Комплексная передаточная величина в генри не является мгновенной вещественной энергетической индуктивностью.

Вывод формы частотного отклика

В принятой схеме ветвь намагничивания связана с короткозамкнутой роторной ветвью. Её реакция имеет первый порядок по параметру \omega_2\tau_{\text{р}}. Умножив дробь на комплексно-сопряжённый знаменатель, получаем вещественную и мнимую части. Мнимая часть отрицательна при соглашении \exp(j\omega t); она отражает фазовый сдвиг наведённого роторного отклика, а не отрицательную запасаемую энергию.

При малой частоте роторная реакция мала и отношение \Psi/I стремится к L_{\text{расс}}+L_{\text{м}}. При большой частоте в данной простой схеме остаётся преимущественно рассеяние. Максимум модуля мнимой части возникает при \omega_2\tau_{\text{р}} = 1. Это даёт первый ориентир постоянной времени, который следует проверить дополнительной частотной точкой, не использованной для подбора.

К странице приложена одна фактически рассчитанная точка на 1 \,\text{Гц} исходной четверти Meeker-модели. Табличная частотная серия, температурные изменения и работа на скорости являются учебными заданиями и отдельными аналитическими проверками. Их нельзя считать уже выполненными по этой единственной карте поля. Полное фазное потокосцепление и момент требуют выведенного секторного восстановления.

\operatorname{Re}\mathcal L=L_{\text{расс}}+\frac{L_{\text{м}}}{1+x^2},\qquad\operatorname{Im}\mathcal L=-\frac{L_{\text{м}}x}{1+x^2},\qquad x=\omega_2\tau_{\text{р}}
(23.5)
Частота в герцах переводится в \omega_{2} = 2\pi f_2, иначе размерность постоянной времени нарушается.

Почему комплексная индуктивность зависит от частоты

Эквивалентная роторная ветвь описывается постоянной времени, а поле решается в неподвижном сечении на частоте относительно ротора.

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

  2. Умножьте дробь на комплексно-сопряжённый знаменатель. Действительная часть убывает, мнимая отрицательна при выбранном временном соглашении.

    \frac1{1+jx}=\frac{1-jx}{1+x^2},\qquad x=\omega_2\tau_{\text{р}}
    (23.6)
  3. Оцените низкочастотный и высокочастотный пределы и положение максимума модуля мнимой части. По серии частот идентифицируйте параметры, затем проверьте их на дополнительных частотах.

Таблица 23.4. Обозначения и единицы
ВеличинаФизический смыслЕдиница СИ
LИндуктивность при явно указанной линеаризацииГн
\widehat IКомплексная амплитуда токаА
\PsiПотокосцепление обмоткиВб
\omega=2\pi fУгловая электрическая частотарад/с
\tau_{\text{р}}Электромагнитная постоянная времени роторас

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

Гармонический параметр и квазистатическая производная

Отношение комплексных амплитуд \widehat\Psi/\widehat I зависит от частоты, потерь и принятой гармонической постановки. Оно не становится автоматически динамической индуктивностью \partial\Psi/\partial I магнитостатической карты. В линейном без потерь квазистатическом пределе оба способа согласуются с единственной L; при вихревых токах или эффективной нелинейной гармонической модели сравнивать их можно только после согласования режима.

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

Четыре точки простой схемы

Дано. Учебные параметры L_{\text{м}} = 0{,}30 \,\text{Гн}; L_{\text{расс}} = 0{,}015 \,\text{Гн}; \tau_{\text{р}} = 0{,}15 \,\text{с}.

  1. Для частот 0,5; 1; 2; 3 \,\text{Гц} вычисляем x = 2\pi f_2\cdot 0{,}15.

  2. Комплексные значения \Psi/I соответственно: 0{,}260486-j0{,}115683; 0{,}173876-j0{,}149737; 0{,}080890-j0{,}124199; 0{,}048354-j0{,}094307 \,\text{Гн}.

  3. Частота максимума модуля мнимой части равна 1/(2\pi \tau_{\text{р}}) = 1{,}06103 \,\text{Гц}; там \operatorname{Im}(\Psi/I) = -0{,}150 \,\text{Гн}.

Ответ: Четыре точки показывают изменение фазового отклика; значение около 1 \,\text{Гц} наиболее чувствительно к постоянной времени.

Проверка: При f_{2} \to 0 получаем 0{,}315 \,\text{Гн}, при f_{2} \to \infty — 0{,}015 \,\text{Гн}. Это аналитическая проверка, а не опубликованный набор результатов FEMM.

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

Четыре ручных частотных опыта

Таблица 23.5. Исходные данные практикума
Исходная величинаЗначениеЕдиница
Число пар полюсов21
Пазы статора / стержни ротора36 / 281
Диаметр ротора80мм
Радиальный зазор0,375мм
Активная длина100мм
Удельная электрическая проводимость алюминия при исходной температуре34,45МСм/м
  1. Откройте учебную модель Асинхронный-двигатель-09-Готовая-модель.fem из официального набора Induction Motor Example. Сохраните новую копию и проверьте распознавание всех материалов, цепей, знаков витков и граничных свойств в FEMM 5; перенос старого файла не считается выполненным расчётом.
  2. Проверьте плоскую постановку, миллиметры и глубину 100 \,\text{мм}. Исходная модель описывает четверть машины: восстановление результата зависит от величины, а не от единого множителя для всех полей.
  3. Назначьте комплексные амплитуды сохранённой модели: A = 1 \,\text{А}, B = -0{,}5 + j0{,}866025 А, C = -0{,}5 - j0{,}866025 А (углы 0^\circ, +120^\circ, -120^\circ). В линейном учебном варианте оставьте проводимость статора нулевой, а алюминиевые стержни проводящими.
  4. Сравните частотную постановку с графическим отображением материалов: проверьте, что именно стержни создают наведённые токи и принят путь их замыкания. Сохраните журнал параметров каждого ручного расчёта.
  5. Сохраните четыре копии с частотами 0,5; 1; 2; 3 \,\text{Гц}. Во всех используйте одну и ту же амплитуду тока и материалы.
  6. После каждого расчёта перепишите комплексное потокосцепление и ток одной фазы из результатов цепей. Выполните комплексное деление и сохраните вещественную/мнимую части.
  7. Выведите восстановление полного фазного потокосцепления для своей четверти машины. Не умножайте локальную B или напряжение той же цепи автоматически на четыре.
  8. Сравните форму зависимости со схемой. Подберите ориентировочные L_{\text{м}}, L_{\text{расс}} и \tau_{\text{р}} и проверьте пятую независимую точку 1{,}5 \,\text{Гц}.

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

Как проверить. Проверьте один токовый масштаб в линейной модели и три уровня сетки. Явно укажите, что полученные коэффициенты относятся к выбранной амплитуде и модели ротора.

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

Как частота ротора меняет получаемый отклик?

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

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

Итог работы. Частотная таблица с объяснением связи с принятой моделью скольжения.

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

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

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

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

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

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

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

Таблица 23.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: Асинхронный двигатель: сектор модели. Это экспорт результатов поля, а не снимок окна программы.
Рисунок 23.5. Визуализация численных данных FEMM 5: Асинхронный двигатель: сектор модели. Это экспорт результатов поля, а не снимок окна программы.

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

Ориентиры в интерфейсе
Геометрия: адаптированная модель Meeker: 36 пазов статора, 28 стержней, неподвижный AC-опыт
Рисунок 23.6. Геометрия: адаптированная модель Meeker: 36 пазов статора, 28 стержней, неподвижный AC-опыт Открыть в полном размере
Распределение магнитной индукции: адаптированная модель Meeker: 36 пазов статора, 28 стержней, неподвижный AC-опыт
Рисунок 23.7. Распределение магнитной индукции: адаптированная модель Meeker: 36 пазов статора, 28 стержней, неподвижный AC-опыт Открыть в полном размере

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

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

  • Асинхронная машина: частота токов ротора, импеданс и потери

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

    Решений FEMM: 4 · Численных проверок: 9

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

Окна FEMM 5

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

Асинхронная машина: расчётный сектор 90° со статорными пазами, воздушным зазором и стержнями ротора. Подписи скрыты для читаемости.
Асинхронная машина: расчётный сектор 90° со статорными пазами, воздушным зазором и стержнями ротора. Подписи скрыты для читаемости. Открыть в полном размере
Асинхронная машина: поле расчётного сектора при частоте скольжения 2 Гц. Норма комплексной амплитуды магнитной индукции B, Т; контуры A скрыты.
Асинхронная машина: поле расчётного сектора при частоте скольжения 2 Гц. Норма комплексной амплитуды магнитной индукции B, Т; контуры A скрыты. Открыть в полном размере

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

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

Критерий выполнения: Указаны частота, ток, комплексное потокосцепление, нормировка сектора и остаточная ошибка независимой точки; серия не названа временным расчётом пуска.
Первоисточники и дополнительное чтение по занятию
  1. [33] — 24–39; 45–52
  2. [39] — Разделы Model Formulation, Identification и Determination of Operating Point; исходная модель FEMM 4, перенос проверяется отдельно
  3. Скачанный оригинал Induction Motor Example — 6–9: расчёты при разных частотах и параметры
  4. Carbonieri, Bianchi. Rapid Procedure for Squirrel Cage Induction Machine Finite Element Analysis — 1–16: дополнительный методический источник по геометрии и режимам; при ручном выполнении расчёта внешний код MATLAB не требуется
  5. [37]
Пояснения и типичные ошибки
Отрицательная мнимая часть \Psi/I означает отрицательную магнитную энергию?

Нет. Это комплексная передаточная величина с фазовым запаздыванием в системе с потерями.

Почему в частотной серии нельзя менять ток одновременно с частотой?

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

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

  • Вычисление \tau_{\text{р}} в секундах из частоты в герцах без 2\pi.
  • Использование AC-отношения \Psi/I как динамической индуктивности магнитостатической карты.
  • Молчаливый повторный секторный множитель.