FEMM 5

Часть V · Электрические машины → § 30 · Насыщенные генераторные режимы

Занятие 30.1

Насыщение и таблицы потокосцеплений по продольной и поперечной осям генератора

Цель и задачи

Цель. Получить нелинейные dq-карты с перекрёстным насыщением.

Задачи

  • Строить карты d/q и различать геометрическую явнополюсность и насыщение.
  • Извлечь \Psi_d и \Psi_q в нескольких независимых токовых направлениях.
Что нужно вспомнить перед началом

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

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

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

Инженерный вопрос. Зависит ли отклик одной оси от тока другой оси?

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

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

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

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

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

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

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

  3. Карты dq с перекрёстным насыщением. Выполните действия практикума по его исходной таблице. Сохраните токи d и q, оба потокосцепления и локальные производные контрольных точек.

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

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

Результат этого занятия. Карты с демонстрацией перекрёстного насыщения и диапазона применимости.

Данные именно этого опыта. Сохраните токи d и q, оба потокосцепления и локальные производные контрольных точек.

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

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

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

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

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

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

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

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

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

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

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

Коэффициенты поля возбуждения явнополюсного генератора

Коэффициенты поля связывают геометрию полюса с основной ЭДС, полным потоком и потребной МДС возбуждения.

Лабораторные работы по методу: Определение коэффициентов поля явнополюсной машины

Получаемые величины. k_f = B_{1}/B полюс; k_\Phi = Φполный/\Phi 1; \alpha \delta = B ср/B полюс; коэффициент насыщения k_z

  1. Рассчитываем ненасыщенную геометрию при заданном токе возбуждения.
  2. Интегрируем поток на полюсном делении и разложить B по гармоникам.
  3. Вводим зубцы, ярмо, проверенную B(H) и заполнение сталью.
  4. Сравниваем коэффициенты и требуемую МДС при одном целевом потоке.
\tau=\frac{\pi D_i}{2p},\qquad \alpha_{\mathrm{геом}}=\frac{b_p}{\tau}
\tau — длина полюсного деления на окружности расточки; для неперекрывающихся полюсов необходимо b_p\le \tau. Это необходимое, но не достаточное условие геометрической пригодности.
Численный пример применения метода

Для явнополюсной машины Di=6,50 м, p=24, ширина полюсного наконечника 306 мм. Определим полюсное деление и перекрытие.

  1. Число полюсов равно 2p=48.
  2. Полюсное деление \tau=πDi/(2p)=425{,}424 мм.
  3. Относительное перекрытие полюсным наконечником \alpha =306/425{,}424=0{,}71928. Межполюсное расстояние по дуге равно 119,424 мм.

Ответ. \tau \approx 425{,}4 мм; \alpha \approx 0{,}7193. Эти величины задают геометрию рабочего зазора перед расчётом формы поля.

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

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

Продольная и поперечная проводимости явнополюсной машины

Коэффициенты k_d и k_q объясняют различие продольного и поперечного магнитного отклика и помогают перейти к индуктивным параметрам явнополюсной машины.

Лабораторные работы по методу: Определение коэффициентов поля явнополюсной машины

Получаемые величины. k_d = B_{\mathrm{d\,}}{1}/Bsmooth; k_q = B_{\mathrm{q\,}}{1}/Bsmooth; отношение продольного и поперечного откликов

  1. Формируем гармоническую МДС по оси d и получаем B_{\mathrm{d\,}}{1}.
  2. Поворачиваем МДС на 90 электрических градусов при прежней амплитуде и получаем B_{\mathrm{q\,}}{1}.
  3. Выполняем опорный расчёт гладкого минимального зазора.
  4. Рассчитываем коэффициенты, затем сопоставляем их с картами \Psi _d(i_d, i_q), \Psi _q(i_d, i_q).

Страницы источника. PDF 1–2: определения; 3–4: картины d/q; 5: методика; 6: варианты

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

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

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

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

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

Продольную ось d выбирают вдоль оси поля ротора; q сдвигают на электрические 90^\circ. Для симметричного фазного набора преобразование координат выделяет составляющие, удобные для описания синхронной машины. Оно не меняет физические токи, а выражает их в другой системе. В линейном ориентире \Psi_d = L_d i_d + \Psi_{\text{возб}}, \Psi_q = L_q i_q. При насыщении \Psi_d и \Psi_q зависят от обоих токов и тока возбуждения, поэтому нужна карта и локальная матрица Якоби.

У неявнополюсного турбогенератора разница между d и q обычно меньше, чем у явнополюсного гидрогенератора; реальное отношение определяется геометрией и данными, а не названием машины. Карты строят при фиксированной температуре и выбранном токе возбуждения либо расширяют третьей осью. Из напряжений v_d = R_si_d - \omega_{\text{эл}} \Psi_q и v_q = R_si_q + \omega_{\text{эл}} \Psi_d видно, как скоростная ЭДС появляется через производную потокосцепления при вращении координат. Статическое FEMM-напряжение не заменяет эти члены.

\Psi_d=L_di_d+\Psi_{\text{возб}},\quad\Psi_q=L_qi_q,\quad T_{\text{эм}}=\frac32p(\Psi_di_q-\Psi_qi_d)
(30.1)
Амплитудно-инвариантная система dq с множителем 2/3; это момент синусоидальной модели основной гармоники без явной зависимости энергии от угла при фиксированных осевых токах. Токи положительны внутрь статора. В насыщенном случае используются обе согласованные карты потокосцеплений; зубцовый момент требует угловой проверки.

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

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

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

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

Таблицы генератора и перекрёстные производные

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

  1. Зафиксируйте преобразование фазных токов в осевые. Для каждой пары осевых токов восстановите фазные токи и выполните полевое решение.

  2. Получите оба потокосцепления и вычислите момент по общей билинейной формуле. Сравните его с полевым интегралом в выбранных контрольных состояниях.

  3. Частные производные таблиц по двум токам формируют локальную матрицу. Ненулевые смешанные элементы показывают перекрёстное насыщение и требуют совместного решения цепи.

Таблица 30.1. Обозначения и единицы
ВеличинаФизический смыслЕдиница СИ
\PsiПотокосцепление обмотки; тот же физический смысл, что у \PsiВб
LИндуктивность при явно указанной линеаризацииГн
TЭлектромагнитный моментН·м
pЧисло пар полюсов1
i_d,i_qОсевые токи в выбранной амплитудной нормировкеА

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

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

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

Дано. p = 2; L_d = 0{,}08 \,\text{Гн}; L_q = 0{,}05 \,\text{Гн}; \Psi_{\text{возб}} = 1 \,\text{Вб}; i_d = -2 \,\text{А}; i_q = 5 \,\text{А}.

  1. \Psi_d = 0{,}08 \cdot (-2) + 1 = 0{,}84 \,\text{Вб}; \Psi_q = 0{,}05 \cdot 5 = 0{,}25 \,\text{Вб}.

    T=\frac32\cdot2[0{,}84\cdot5-0{,}25(-2)]=14{,}1\ \text{Н}\cdot\text{м}
    (30.2)
  2. T_{\text{эм}} = 3[0{,}84 \cdot 5 - 0{,}25 \cdot (-2)] = 14{,}1 \,\text{Н}\cdot \text{м}.

  3. Та же формула раскрывается как 3[\Psi_{\text{возб}} i_q + (L_d - L_q)i_d i_q], показывая возбуждённую и реактивную части линейного ориентира.

Ответ: Момент линейного аналитического эталона равен +14{,}1 \,\text{Н}\cdot \text{м}. При положительной скорости и токах внутрь статора это двигательный знак момента. Пример проверяет алгебру общей карты; для генераторного режима в той же пассивной ориентации электромагнитный момент должен тормозить ротор.

Проверка: Знак момента меняется при изменении знака i_q; число пар полюсов относится к электрическому углу.

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

Карты dq с перекрёстным насыщением

Таблица 30.2. Исходные данные практикума
Исходная величинаЗначениеЕдиница
Сетка i_d-4; -2; 0; 2А
Сетка i_q0; 2; 5А
Ток возбуждениязаданныйА
  1. Выберите ранее рассмотренные модели турбогенератора и явнополюсного генератора; укажите оси d/q и число пар полюсов.
  2. Выведите обратное преобразование токов и задавайте физические фазные токи для каждой точки карты.
  3. Во встроенном Python считывайте три фазных \Psi и преобразуйте их в выбранную d/q-нормировку.
  4. В каждой модели постройте низкотоковый ориентир и насыщенную карту при одном токе возбуждения.
  5. Получите локальные производные, проверяя шаг тока. Сравните зависимость L_d и L_q без принудительного приравнивания.
  6. Добавьте второй ток возбуждения и отделите его влияние от изменения геометрии.

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

Как проверить. Прямое и обратное преобразования восстанавливают фазные данные; низкотоковый момент имеет независимый ориентир.

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

Зависит ли отклик одной оси от тока другой оси?

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

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

Итог работы. Карты с демонстрацией перекрёстного насыщения и диапазона применимости.

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

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

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

Python · встроенный MicroPython FEMM 5
"""Насыщение и таблицы потокосцеплений по продольной и поперечной осям генератора.

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

# Линейные продольная и поперечная индуктивности, Гн; потокосцепление возбуждения, Вб.
# Токи по осям d и q заданы амплитудными значениями, А.
pole_pairs, ld_h, lq_h, excitation_wb = 2, 0.08, 0.05, 1.0
id_a, iq_a = -2.0, 5.0
# Ось d направлена вдоль возбуждения; ось q — поперёк. Потокосцепления выражены в Вб.
lambda_d = ld_h * id_a + excitation_wb
lambda_q = lq_h * iq_a
# Коэффициент 3/2 соответствует выбранному преобразованию трёх фаз и амплитудным токам.
torque_nm = 1.5 * pole_pairs * (lambda_d * iq_a - lambda_q * id_a)
# Порядок вывода: продольное потокосцепление, Вб; поперечное потокосцепление, Вб; момент, Н·м.
print(lambda_d, lambda_q, torque_nm)
# Скоростное напряжение добавляется отдельно:
# vd = resistance*id_a - omega_electrical*lambda_q
# vq = resistance*iq_a + omega_electrical*lambda_d

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

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

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

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

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

Таблица 30.3. Результаты расчётного образца
ВеличинаРезультатЕдиница
Частота0Гц
Активная длина80мм
Ток цепи Фаза А1А
Потокосцепление цепи Фаза А0,004146Вб
Ток цепи Фаза В1А
Потокосцепление цепи Фаза В0,03349Вб
Ток цепи Фаза С-2А
Потокосцепление цепи Фаза С-0,03512Вб
Ток цепи Обмотка возбуждения2А
Потокосцепление цепи Обмотка возбуждения0,1368Вб
Момент на исходной учебной сетке-0,04857Н·м
Энергия поля0,1908Дж
Модуль B в контрольной точке0,1114Тл
Координаты контрольной точки (x, y)56,42; 17,25мм
Момент на проверенной уточнённой сетке-0,04962Н·м
Визуализация численных данных FEMM 5: Синхронная машина с явнополюсным ротором. Это экспорт результатов поля, а не снимок окна программы.
Рисунок 30.1. Визуализация численных данных FEMM 5: Синхронная машина с явнополюсным ротором. Это экспорт результатов поля, а не снимок окна программы.

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

Ориентиры в интерфейсе
Геометрия: цилиндрический ротор с токовым возбуждением; неподвижный снимок, скоростная ЭДС отдельно
Рисунок 30.2. Геометрия: цилиндрический ротор с токовым возбуждением; неподвижный снимок, скоростная ЭДС отдельно Открыть в полном размере
Распределение магнитной индукции: цилиндрический ротор с токовым возбуждением; неподвижный снимок, скоростная ЭДС отдельно
Рисунок 30.3. Распределение магнитной индукции: цилиндрический ротор с токовым возбуждением; неподвижный снимок, скоростная ЭДС отдельно Открыть в полном размере
Геометрия: явнополюсный ротор с токовым возбуждением; отдельная геометрия
Рисунок 30.4. Геометрия: явнополюсный ротор с токовым возбуждением; отдельная геометрия Открыть в полном размере
Распределение магнитной индукции: явнополюсный ротор с токовым возбуждением; отдельная геометрия
Рисунок 30.5. Распределение магнитной индукции: явнополюсный ротор с токовым возбуждением; отдельная геометрия Открыть в полном размере
Синхронная машина с явнополюсным ротором: выбранная сетка воздушного зазора.
Рисунок 30.6. Синхронная машина с явнополюсным ротором: выбранная сетка воздушного зазора. Открыть в полном размере
Синхронная машина с явнополюсным ротором: момент ротора на выбранной сетке.
Рисунок 30.7. Синхронная машина с явнополюсным ротором: момент ротора на выбранной сетке. Открыть в полном размере
Встроенная консоль FEMM: Насыщение и таблицы потокосцеплений по продольной и поперечной осям генератора. Аналитическая проверка программы занятия. Показана консоль исходной редакции программы, зафиксированной при съёмке. Текущая программа с русскими пояснениями проверена отдельно во встроенном Python FEMM; этот снимок не является новым запуском текущей редакции.
Рисунок 30.8. Встроенная консоль FEMM: Насыщение и таблицы потокосцеплений по продольной и поперечной осям генератора. Аналитическая проверка программы занятия. Показана консоль исходной редакции программы, зафиксированной при съёмке. Текущая программа с русскими пояснениями проверена отдельно во встроенном Python FEMM; этот снимок не является новым запуском текущей редакции. Открыть в полном размере

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

Сопоставьте локальные d/q-матрицы турбо- и гидрогенератора в низкотоковой и насыщенной точках; объясните различие.

Критерий выполнения: Оси, нормировка и ток возбуждения указаны; получены две модели и два уровня возбуждения; шаг производной проверен.
Первоисточники и дополнительное чтение по занятию
  1. [33] — 10–15, 49–52, 105–110
  2. [81] — Раздел Packages · Локальная копия
  3. [19] — Выделенный Bianchi_FEMM_Lua_section.pdf, с. 68–71 и 84–103; исходный PDF, с. 103–106 и 119–138
Пояснения и типичные ошибки
Почему \Psi_d зависит от i_q при насыщении?

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

Можно ли считать L_d = L_q для любого турбогенератора?

Нет. Это возможное приближение, которое проверяется по конкретной геометрии и рабочей точке.

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

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