FEMM 5

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

Занятие 37.14

Регулировочная характеристика и потребное возбуждение

Цель и задачи

Цель. Определить потребное возбуждение для поддержания напряжения.

Задачи

  • Поставить регулировочный опыт с фиксированным U и \cos \varphi.
  • Различать регулировочную кривую и алгоритм автоматического регулятора возбуждения.
Что нужно вспомнить перед началом

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Проверка из практикума. U и \cos \varphi удержаны в объявленных допусках; недостижимые точки выделены; статическая кривая не названа динамикой AVR.

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

Критерий выполнения. U и \cos \varphi удержаны в объявленных допусках; недостижимые точки выделены; статическая кривая не названа динамикой AVR.

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

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

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

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

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

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

Регулировочная характеристика If(I) определяет возбуждение, необходимое для поддержания напряжения при заданном характере нагрузки. Продолжение по режиму сопровождается проверкой ограничений и невязок в каждой расчётной точке.

Соглашения об обозначениях и единицах
  • \(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.116) · Продолжение и ограничения
  2. (37.117) · Вывод потребного возбуждения и корректного фазового угла
  3. (37.118) · Вывод потребного возбуждения и корректного фазового угла
  4. (37.119) · Вывод потребного возбуждения и корректного фазового угла
  5. (37.120) · Вывод потребного возбуждения и корректного фазового угла
  6. (37.121) · Вывод потребного возбуждения и корректного фазового угла

Условия регулировочной кривой

Регулировочная характеристика показывает потребный ток возбуждения для поддержания заданного напряжения при изменении тока нагрузки. Частота и коэффициент мощности задаются явно; при вещественном целевом U фаза тока \beta I фиксирована, а механический угол ротора \alpha r ищется вместе с If. Это набор установившихся точек, а не модель быстродействия регулятора.

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

Продолжение и ограничения

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

Физические пределы тока возбуждения и статора потребовали бы заявленных изготовителем тепловых данных. Числовая граница поиска If=0,1…150 А в программе является учебным ограничением алгоритма, а допустимая невязка U — его критерием. Если ограничение достигнуто раньше требуемого режима, точка не считается успешно регулируемой.

\Re U(I_f,\alpha_r;I,\beta_I)=U_*,\qquad \Im U=0,\qquad \Delta U=100\frac{U_0-U_N}{U_N}
(37.116)
Первые два условия определяют If и механический угол \alpha r при заданной фазе тока \beta I; \Delta U — отдельный показатель изменения напряжения без регулирования при неизменном возбуждении.
  1. Если при неизменном If напряжение ХХ 250 В, а под нагрузкой 230 В, \Delta U=100\cdot 20/230.

  2. Получаем 8,69565 %. Это не наклон регулировочной характеристики If(I).

Переход к регулятору как отдельная задача

Чтобы моделировать AVR, нужны измеритель напряжения, фильтры, усиление, ограничения, возбудитель и динамика цепи возбуждения. Кривая If(I) задаёт статическую потребность и может служить опорой, но не даёт постоянной времени и устойчивости регулятора.

Проверяйте чувствительность результата к допуску обратного поиска и полевой сетке отдельно. Слишком грубое поле способно дать малую невязку U относительно самого себя, сохраняя заметную ошибку фактической кривой.

Вывод потребного возбуждения и корректного фазового угла

Для независимого аналитического примера берём U вещественным, ток I(\cos \varphi -j \sin \varphi ) и схему E=U+(R+jX)I. Знак \varphi различает отстающий и опережающий ток.

В статье 2019 в одном напечатанном скалярном модуле встречается знак минус между квадратами. Здесь модуль выводится непосредственно из комплексного вектора и использует сумму квадратов; угол определяется atan2, сохраняющей квадрант.

Связь с программой. target()/target_signed() в шаге 13 удерживают комплексное напряжение. Программа использует фундаментальную компоненту поля; линейное выражение выше служит независимым расчётом знаков и тенденций.

  1. Раскрываем произведение (R+jX)I(\cos \varphi -j \sin \varphi ) и собираем вещественную и мнимую составляющие.

    E_x=U+RI\cos\varphi+XI\sin\varphi,\qquad E_y=XI\cos\varphi-RI\sin\varphi
    (37.117)
  2. Умножаем комплексную ЭДС на сопряжённую: мнимые перекрёстные члены сокращаются, поэтому под корнем именно плюс.

    |E|^2=(E_x+jE_y)(E_x-jE_y)=E_x^2+E_y^2
    (37.118)
  3. Положительный модуль и угол восстанавливаются из двух ортогональных проекций. \arctan(Ey/Ex) без квадранта способен выбрать неверную ветвь.

    E=\sqrt{E_x^2+E_y^2},\qquad\delta=\operatorname{atan2}(E_y,E_x)
    (37.119)
  4. При линейной связи E=kE If потребное возбуждение получается делением модуля. При насыщении используется обратная полевая задача.

    I_f=\frac1{k_E}\sqrt{(U+RI\cos\varphi+XI\sin\varphi)^2+(XI\cos\varphi-RI\sin\varphi)^2}
    (37.120)
  5. В нелинейной реализации целевое U вещественно, фаза тока \beta I задана коэффициентом мощности; две невязки определяют If и механический угол \alpha r. Продолжение передаёт именно эту пару от предыдущей нагрузки.

    \mathbf r(I_f,\alpha_r;I,\beta_I)=0,\qquad\mathbf x_0(I_{k+1})=\mathbf x_*(I_k),\quad\mathbf x=(I_f,\alpha_r)^T
    (37.121)
Таблица 37.41. Обозначения и единицы
ВеличинаФизический смыслЕдиница СИ
E_x,E_yОртогональные проекции ЭДСВ
k_EЛинейный коэффициент возбужденияВ/А
\deltaФазовый угол, определённый atan2рад
I_fПотребное возбуждениеА
U,IФазное напряжение и ток отдачиВ действующее; А действующее
R,XПараметры линейной опорыОм
\varphiУгол напряжения относительно токарад
\beta_IЗаданная фаза тока в неподвижной системе программыградус в интерфейсе
\alpha_rИскомый механический угол ротораградус в интерфейсе

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

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

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

Дано. При одном If напряжение ХХ U_{0}=250 В, при нагрузке UN=230 В; используется фазное действующее напряжение.

  1. Разность равна 20 В.

  2. \Delta U=100\cdot (250-230)/230=8{,}695652 %.

Ответ: Изменение напряжения без регулирования равно 8,69565 % относительно UN.

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

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

Шаг 13: Регулировочная характеристика и потребное возбуждение

Таблица 37.42. Исходные данные практикума
Исходная величинаЗначениеЕдиница
Учебная модельдвухполюсный цилиндрический ротор; 24 открытых паза статора, 12 открытых пазов ротора; отдельные немагнитные клинья
СредаFEMM 5 · встроенный MicroPython · магнитостатика позиций
Частота изменения положения50Гц
  1. Выполните run-13.py; проверьте целевое напряжение, заданный токовый ряд и ветвь \cos \varphi.
  2. Постройте If(I), угол и невязку напряжения; ограничения покажите на той же таблице.
  3. Для каждой точки проверьте фактический \cos \varphi из U·Iout*, а не только введённое значение.
  4. Сравните два характера нагрузки в отдельных каталогах при одинаковых U, f и геометрии.
  5. Для одной точки повторите поиск без использования предыдущего приближения и сравните результат.

Ожидаемый результат. Регулировочная характеристика с картой ограничений и сохранёнными невязками.

Как проверить. U и \cos \varphi удержаны в объявленных допусках; недостижимые точки выделены; статическая кривая не названа динамикой AVR.

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

Как нагрузка меняет необходимый ток возбуждения?

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

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

Итог работы. Регулировочная характеристика с проверкой поддерживаемого условия и допустимости найденных режимов.

Встроенный MicroPython: Регулировочная характеристика и потребное возбуждение · Порядок запуска

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

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

Python · встроенный MicroPython FEMM 5
"""Регулировочная характеристика и потребное возбуждение.

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

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

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

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

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

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

Сравните регулировочные кривые для отстающего и опережающего токов; укажите источник ограничения If и объясните изменение реакции якоря.

Критерий выполнения: U и \cos \varphi удержаны в объявленных допусках; недостижимые точки выделены; статическая кривая не названа динамикой AVR.
Первоисточники и дополнительное чтение по занятию
  1. [65] — С. 1–12; ХХ, КЗ, регулировочные и U-образные характеристики; в статье рассмотрен турбогенератор 325 МВт.
  2. [107] — С. 36–41; Методическая связь: согласование поля и выходного режима; реализован собственный алгоритм учебной модели.
  3. [33] — Магнитная плоская задача, цепи mi_/mo_, встроенный Python и постпроцессор.
Пояснения и типичные ошибки
Почему маленькая невязка U ещё не доказывает точность характеристики?

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

Можно ли по If(I) определить постоянную времени возбудителя?

Нет. Установившаяся кривая не содержит временного уравнения и параметров возбудителя.

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

  • Регулировочная характеристика и относительное изменение напряжения — разные определения.
  • Продолжение не доказывает отсутствия другой ветви решения.