FEMM 5

Часть V · Электрические машины → § 37 · Турбогенераторы: от параметрической геометрии до характеристик, динамики и дефектов

Занятие 37.15

U-образные характеристики, реактивная мощность и угол нагрузки

Цель и задачи

Цель. Связать U-образную характеристику с реактивной мощностью.

Задачи

  • Выбрать неизвестные для условий U=\mathrm{const} и P=\mathrm{const}.
  • Проверить минимум тока около единичного коэффициента мощности.
Что нужно вспомнить перед началом

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

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

Расчётная лабораторная работа · Восстановление электрических характеристик синхронного генератора

Инженерный вопрос. Как минимум тока и знак реактивной мощности зависят от возбуждения?

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

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

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

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

  1. Обосновать постановку и способ измерения. Для характеристики запишите, что фиксировано: скорость, возбуждение, напряжение, нагрузка или мощность. Скоростная ЭДС восстанавливается из пространственного потокосцепления.

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

  3. Шаг 14: U-образные характеристики, реактивная мощность и угол нагрузки. Выполните действия практикума по его исходной таблице. Сведите возбуждение, ток статора, реактивную мощность и угол при согласованных внешних условиях.

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

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

Результат этого занятия. U-образная характеристика с отмеченным минимумом и проверкой знака реактивной мощности.

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

Оформление наблюдений. Возбуждение, скорость, фазные потокосцепления, ЭДС и согласованные токи/напряжения в объёме конкретного опыта; независимые условия каждой строки.

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

Проверка из практикума. Удержаны U и P; минимум тока связан с Q\approx 0; электрическая достижимость отделена от тепловой и системной допустимости.

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

Критерий выполнения. Удержаны U и P; минимум тока связан с Q\approx 0; электрическая достижимость отделена от тепловой и системной допустимости.

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

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

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

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

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

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

Зависимость I(If) строится при постоянных напряжении и активной мощности. Переход через Q=0 связывает форму U-образной характеристики с реактивным обменом; достижимые режимы сопоставляются с приближённой линейной моделью.

Соглашения об обозначениях и единицах
  • \(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. (37.122) · Линейная проверка
  2. (37.123) · Вывод U-кривой и системы нелинейного поиска
  3. (37.124) · Вывод U-кривой и системы нелинейного поиска
  4. (37.125) · Вывод U-кривой и системы нелинейного поиска
  5. (37.126) · Вывод U-кривой и системы нелинейного поиска
  6. (37.127) · Вывод U-кривой и системы нелинейного поиска
  7. (37.128) · Вывод U-кривой и системы нелинейного поиска

Почему появляется буква U

При фиксированных U и P активная составляющая тока постоянна, а реактивная меняется с возбуждением. Поэтому модуль тока минимален при Q=0, если такая точка достижима. Недовозбуждённые и перевозбуждённые ветви имеют противоположные знаки Q в принятом генераторном соглашении.

Готовая программа использует Q как параметр: при каждом заданном P,U задаются четыре Q, из Iout=(P−jQ)/(3U) вычисляется ток, а полевой поиск находит If и угол ротора. После этого точки строятся в координатах If,I. Таким образом условия U=\mathrm{const} и P=\mathrm{const} проверяются полем, хотя непосредственно перебирается Q. Четыре точки — редкая учебная сетка, её недостаточно для точного положения минимума.

Линейная проверка

Для линейной цилиндрической модели без R условие фиксированной активной мощности даёт E \sin \delta=PX/(3U). Реактивная составляющая определяется E \cos \delta -U. Условие |\sin \delta |\le 1 сразу показывает часть недостижимого диапазона возбуждения.

Формула минимума Imin=P/(3U) служит сильной проверкой нормировки фазного напряжения и действующего тока. В нелинейной модели с потерями и гармониками коэффициент мощности и минимальный ток надо определять по согласованной полной мощности; основной гармоники может быть недостаточно.

I=\frac{\sqrt{P^2+Q^2}}{3U},\qquad I_{min}=\frac{P}{3U}\ (Q=0),\qquad \sin\delta=\frac{PX_s}{3UE}
(37.122)
Симметричная синусоидальная линейная модель; U — действующее фазное В, P,Q — трёхфазные Вт и вар.
  1. При P=6900 Вт и U=230 В получаем Imin=6900/(3\cdot 230)=10 А.

  2. Если Q=6900 вар, ток возрастёт до 10√2=14{,}1421 А при той же P.

Допустимость не равна существованию решения

Сошедшаяся U-кривая показывает электрически согласованные точки учебной модели. Реальная область работы ограничивается нагревом ротора, статора, торцевых зон, устойчивостью и сетью. Без этих моделей нельзя называть всё семейство допустимой диаграммой мощности.

Сохраните U, P, Q, I, If, \delta, невязки и статусы ограничений. На графике не сглаживайте провалы от неудачных поисков. Для сравнения нескольких P используются одинаковые входы геометрии, сетки и сопротивления.

Вывод U-кривой и системы нелинейного поиска

Рассмотрим симметричный синусоидальный режим при заданных фазном U и трёхфазной активной мощности P. В линейной опоре пренебрегаем R и считаем внутреннюю ЭДС E управляемой возбуждением.

Каждая рассчитанная U-кривая требует одновременно выполнить два условия — U=\mathrm{const} и P=\mathrm{const}. Можно параметризовать возбуждением и искать ток либо параметризовать Q и заранее определить ток. В готовой программе используется второй путь: четыре заданных Q для каждого P, затем полевой поиск If и угла ротора.

Связь с программой. Шаг 14 перебирает четыре Q для двух P, определяет I и его фазу, затем target()/target_signed() находят If и угол ротора. Полученные точки строятся в координатах If,I. Для точного минимума нужна более частая сетка Q; минимум по четырём точкам не является непрерывным оптимумом.

  1. Из S=3UI* получаем квадрат модуля тока. При фиксированных P,U зависимость от Q входит только неотрицательным квадратом.

    I^2=\frac{P^2+Q^2}{9U^2}\ge\frac{P^2}{9U^2}
    (37.123)
  2. Равенство достигается при Q=0. Это даёт точный минимум в данной синусоидальной постановке, если существует соответствующий If.

    I_{min}=\frac{|P|}{3U},\qquad\cos\varphi=1
    (37.124)
  3. Из P=3UE \sin \delta /X выражаем \sin \delta. Для вещественного угла необходим нижний предел модуля E.

    \sin\delta=\frac{PX}{3UE},\qquad E\ge\frac{|P|X}{3U}
    (37.125)
  4. На выбранной устойчивой ветви \cos \delta \ge 0 подставляем корень в Q. Другая ветвь имеет противоположный знак \cos \delta и проверяется отдельно.

    Q(E)=\frac{3U}{X}\left(E\sqrt{1-\left(\frac{PX}{3UE}\right)^2}-U\right)
    (37.126)
  5. В реализованной параметризации Q известен, поэтому модуль и угол тока получаем из заданных мощностей и напряжения.

    \underline I_{out}=\frac{P-jQ}{3U_*},\quad I=\frac{\sqrt{P^2+Q^2}}{3U_*},\quad\beta=-\operatorname{atan2}(Q,P)
    (37.127)
  6. Затем по известному комплексному току решаем две невязки напряжения относительно U*=U*+j0. Искомые — If и угол ротора; U и P контролируются фактическими полевыми значениями.

    \mathbf r(I_f,\theta_m;P,Q)=\frac1{U_*}\binom{\Re U-U_*}{\Im U}=0,\qquad\mathbf J\Delta\mathbf x=-\mathbf r
    (37.128)
Таблица 37.44. Обозначения и единицы
ВеличинаФизический смыслЕдиница СИ
I_{min}Минимальный действующий фазный токА
PЗаданная трёхфазная активная мощностьВт
U_*Заданное действующее фазное напряжениеВ
Q(E)Реактивная мощность линейной опорывар
E,UВнутренняя ЭДС и фазное напряжение линейной опорыВ действующее
XСинхронное сопротивление линейной опорыОм
\delta,\theta_mЭлектрический угол нагрузки и механический угол роторарад
\betaЭлектрическая фаза задаваемого тока отдачирад в формуле; градус в интерфейсе

Границы применимости. Минимум основной гармоники не заменяет учёт гармонического RMS и всех потерь в полной машине.

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

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

Дано. Синусоидальный трёхфазный режим: P=6900 Вт, фазное U=230 В, Q=0.

  1. При Q=0 полная мощность S=P.

  2. Imin=P/(3U)=6900/690=10 А.

Ответ: Минимальный действующий фазный ток равен 10 А при единичном \cos \varphi.

Проверка: Это самостоятельный аналитический эталон. Сравнение с FEMM допустимо при совпадении всех входов, геометрии и определения величины.

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

Шаг 14: U-образные характеристики, реактивная мощность и угол нагрузки

Таблица 37.45. Исходные данные практикума
Исходная величинаЗначениеЕдиница
Учебная модельдвухполюсный цилиндрический ротор; 24 открытых паза статора, 12 открытых пазов ротора; отдельные немагнитные клинья
СредаFEMM 5 · встроенный MicroPython · магнитостатика позиций
Частота изменения положения50Гц
  1. Выполните run-14.py; проверьте два уровня P, четыре Q, вычисление комплексного Iout и полевой поиск If и угла ротора.
  2. Постройте I(If), Q(If) и \delta(If) для каждого P; невязки U и P сохраняйте рядом.
  3. Найдите точки смены знака Q и сравните ток с P/(3U) в пределах принятых потерь.
  4. Отметьте недостижимые точки, ограничения возбуждения и найденную ветвь.
  5. Сравните два уровня P; проверьте одинаковость R, частоты, геометрии и сетки.

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

Как проверить. Удержаны U и P; минимум тока связан с Q\approx 0; электрическая достижимость отделена от тепловой и системной допустимости.

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

Как минимум тока и знак реактивной мощности зависят от возбуждения?

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

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

Итог работы. U-образная характеристика с отмеченным минимумом и проверкой знака реактивной мощности.

Встроенный MicroPython: U-образные характеристики, реактивная мощность и угол нагрузки · Порядок запуска

Выполнить независимый шаг 14

Сохраните reproduce.py и выбранный run-NN.py в одну рабочую папку и прочитайте README.md. Укажите её абсолютный путь в TG_ROOT; пример на странице использует C:/FEMM5/turbogenerator. Запускайте через редактор Python FEMM 5. Каждый шаг восстанавливает необходимые входы без запуска предыдущих; внешний pyFEMM не требуется.

Python · встроенный MicroPython FEMM 5
"""U-образные характеристики, реактивная мощность и угол нагрузки.

Встроенный MicroPython FEMM 5; отдельный воспроизводимый шаг учебного турбогенератора.
Результат относится к адаптированной геометрии, а не к промышленной машине 340 МВт.
"""
# Укажите рабочий каталог; сохраните рядом reproduce.py и прочитайте README.md.
TG_ROOT = r"C:/FEMM5/turbogenerator"
STAGE = "14"
CASE = "standard"
# Программа заново строит нужные предпосылки и сохраняет модели, таблицы и квитанцию шага.
exec(open(TG_ROOT + "/reproduce.py").read(), globals())

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

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

Таблица 37.46. Результаты расчётного образца
ВеличинаРезультатЕдиница
Рассчитанные точки характеристики81
Заданное фазное напряжение80В
График сохранённых численных данных учебного образца; условия и формулы приведены в тексте. Это самостоятельная научная визуализация, а не снимок интерфейса FEMM.
Рисунок 37.11. График сохранённых численных данных учебного образца; условия и формулы приведены в тексте. Это самостоятельная научная визуализация, а не снимок интерфейса FEMM.

Нелинейная 2D магнитостатика активной части и отдельная редуцированная динамика. Проверка сетки относится к напряжению ХХ. Уточнение углового шага и баланс мощности приведены отдельно на этапе 15.

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

Постройте две U-кривые и оцените сдвиг If при Q=0. Для одной точки сравните повторный поиск из разных начальных приближений.

Критерий выполнения: Удержаны U и P; минимум тока связан с Q\approx 0; электрическая достижимость отделена от тепловой и системной допустимости.
Первоисточники и дополнительное чтение по занятию
  1. [65] — С. 1–12; ХХ, КЗ, регулировочные и U-образные характеристики; в статье рассмотрен турбогенератор 325 МВт.
  2. [33] — Магнитная плоская задача, цепи mi_/mo_, встроенный Python и постпроцессор.
Пояснения и типичные ошибки
Почему минимум I при фиксированных U и P соответствует Q=0?

Из I^{2}=(P^{2}+Q^{2})/(9U^{2}) минимум достигается при нулевой реактивной составляющей, если точка достижима.

Можно ли называть полученные кривые полной диаграммой допустимой мощности?

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

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

  • I(If) при фиксированном фазовом токе не является U-характеристикой.
  • Положительный статус поиска не подтверждает тепловую допустимость режима.