Моделирование динамики трехфазной круглороторной или явно полюсной синхронной машины с использованием фундаментальных параметров в блоках СИ
Simscape / Электрический / Специализированные Энергосистемы / Электрические Машины

Основной блок Synchronous Machine SI моделирует синхронную машину в режиме генератора или двигателя с использованием основных параметров в единицах СИ. Режим работы диктуется признаком механической мощности (положительной для режима генератора или отрицательной для режима двигателя). Электрическая часть машины представлена моделью пространства состояния шестого порядка, а механическая часть такая же, как в блоке упрощенной синхронной машины.
Модель учитывает динамику обмоток статора, поля и демпфера. Эквивалентная схема модели представлена в опорной рамке ротора (qd-рамке). Обмотки статора соединены трепом с внутренней нейтральной точкой. Все параметры ротора и электрические величины просматриваются со статора и идентифицируются по заполненным переменным. Подстрочные индексы:
d, q - количество по осям d- и q
R, s - Количество ротора и статора
l, m - Индуктивность утечки и намагничивания
f, k - Количество обмоток возбуждения и демпфера
Электрическая модель машины показана на этих схемах.

Традиционная теория синхронного моделирования машин для анализа устойчивости предполагает, что взаимные индуктивности между якорем, демпфером и полем на обмотках с прямой осью идентичны. Обычно обмотки демпфера находятся вблизи воздушного зазора, и в результате схемы демпфера, связывающие поток, почти равны якорю, связывающему поток. Эта гипотеза дает приемлемые результаты для широкого спектра исследований стабильности, особенно на стороне сети. Однако, когда речь идет о текущих полевых исследованиях, существует значительная ошибка. Эквивалентная схемная динамическая модель синхронной машины может включать в себя дополнительную индуктивность, представляющую разность между взаимными индуктивами «поле-демпфер» и «поле-якорь» на оси D [1]. Эта индуктивность обычно называется индуктивностью Каная. Индуктивность Каная соответствует потоку утечки, ФС, на следующем рисунке и интерпретируется как корректирующий элемент в эквивалентной модели, который может иметь отрицательное значение [2].

В стандарте IEEE 1110-2002 [3] представлены прямые и квадратичные оси динамической модели синхронной машины, как показано на диаграммах.


Соответствующими уравнениями являются:
Блок Synchronous Machine SI Fundamental позволяет задать основные параметры синхронной машины. Параметры обмотки поля и демпфера (сопротивления, индуктивности утечки и взаимные индуктивности) вводятся в SI (Λ, H). Параметры RL обмоток возбуждения и демпфера являются не фактическими значениями RL поля машины, а значениями RL, относящимися к статору.
В этой таблице отображаются базовые значения статора.
Vn23
| Напряжение базового статора = пиковое номинальное напряжение линии-нейтрали (В) |
Pn2Vn3
| Ток базового статора (A) |
Vn2Pn
| Импеданс основного статора (Λ) |
fn
| Базовая угловая частота (рад/с) |
Zsbaseü base
| Индуктивность основного статора (H) |
где:
Pn - трехфазная номинальная мощность (VA).
Vn - номинальное межлинейное напряжение (Vrms).
fn - номинальная частота (Гц).
ifn - номинальный ток поля, вырабатывающий номинальное напряжение статора при отсутствии нагрузки (A).
В блоке задаются следующие параметры статора:
| RS | Сопротивление статора на фазу (Λ) |
| Ll | Индуктивность утечки статора (H) |
| Lmd | Индуктивность намагничивания по прямой оси от статора (Н) |
| Lmq | Квадратурная индуктивность намагничивания со статора (H) |
На этом рисунке показана одна фаза обмотки статора, соединенной с обмоткой возбуждения.

На схеме:
Ns и Nf - эквивалентное число синусоидально распределенных витков обмотки статора и обмотки поля соответственно.
Rs и L1 - сопротивление статора и индуктивность утечки.
Rf ′ и Llfd ′ являются сопротивлением поля и индуктивностью утечки.
Когда три обмотки статора запитаны напряжением трехфазной положительной последовательности и обмотка возбуждения разомкнута, индуктивность намагничивания статора равна Lmd. Однако, когда под напряжением находится только одна фаза и обмотка возбуждения разомкнута, индуктивность намагничивания составляет 2/3 Lmd, как показано на рисунке.
При отсутствии нагрузки, когда обмотка поля вращается с номинальной скоростью и несет номинальный постоянный ток поля ifn, амплитуда напряжения переменного тока (пиковое значение), индуцированного на одной фазе статора, равна Vsbase.
Максимальная взаимная индуктивность между одной обмоткой статора и обмоткой возбуждения получается при совмещении двух обмоток и определяется:
startbase
из чего выводим коэффициент трансформации:
startbaseVsbase.
Коэффициент трансформации может быть также выражен как:
23IfbaseIsbase,
где Ifbase - ток базового поля. В таблице показано, как вычислять базовые значения полей.
Lmd_pu
| Ток базового поля (A) |
PnIfbase
| Напряжение основного поля (В) |
VfbaseIfbase
| Импеданс основного поля (Λ) |
Zfbaseü base
| Индуктивность основного поля (H) |
Фактические параметры поля:
Rf ′ - сопротивление поля (Λ)
Llfd ′ - индуктивность утечки поля (Н)
Можно задать параметры поля для сопротивления поля и индуктивности утечки, относящейся к статору (Rf, Llfd).
Если номинальный ток поля ifn известен, коэффициент преобразования Ns/Nf вычисляется с использованием того же уравнения, что и для коэффициента преобразования статор/поле.
По Краузе [4] напряжение поля и ток, относящийся к статору (Vf, If), выводятся из фактических значений Vf ′, If ′ следующим образом:
NsNfVf '
23NfNsIf'
Когда известны фактическое сопротивление поля Rf ′ и индуктивности утечки Llfd ′ (видимые из ротора), значения, относящиеся к статору:
NsNf) 2
NsNf) 2
Когда номинальный ток поля неизвестен и если известны значения pu сопротивления поля и индуктивностей утечки (Rf_pu, Llfd_pu), соответствующие значения SI, относящиеся к статору, вычисляются следующим образом:
Zsbase
Lsbase
Те же преобразования используются для параметров RL обмоток демпфера.
При указании номинального тока поля сигнал на входном порте Vf соответствует фактическому напряжению поля, как в реальном времени. Ток поля, возвращаемый выходным сигналом измерения, также соответствует фактическому току поля If.
Номинальное полевое напряжение, создающее номинальное напряжение статора при отсутствии нагрузки, определяется:
ifn
Если номинальный ток поля не указан, сигнал на входном порте Vf соответствует фактическому напряжению поля статора. В этом случае номинальное напряжение поля, относящееся к статору, создающему номинальное напряжение статора при отсутствии нагрузки, составляет:
startbaseVsbase.
Ток поля, возвращаемый выходным сигналом измерения, представляет собой ток поля, относящийся к статору. Номинальный ток поля статора:
IsbaseIfbaseifn.
В дискретных системах при установке для параметра модели дискретного решателя блока синхронной машины значения Trapezoidal non iterative, возможно, вам придется подключить небольшую паразитную резистивную нагрузку на клеммах машины, чтобы избежать цифровых колебаний. Большие времена выборки требуют больших нагрузок. Минимальная резистивная нагрузка пропорциональна времени выборки. Как правило, следует помнить, что с шагом времени 25 мкс в системе с частотой 60 Гц минимальная нагрузка составляет приблизительно 2,5% от номинальной мощности машины. Например, синхронная машина 200 МВА в энергосистеме, дискретизированная с временем выборки 50 мкс, требует приблизительно 5% резистивной нагрузки или 10 МВт. Если время выборки уменьшено до 20 мкс, достаточно резистивной нагрузки 4 МВт.
Однако если для параметра модели дискретного решателя блока синхронной машины задано значение Trapezoidal iterative (alg. loop), можно использовать ничтожно малую паразитную нагрузку (ниже 0,1% номинальной мощности) при сохранении численной устойчивости. Эта итеративная модель создает алгебраический цикл и приводит к более низкой скорости моделирования.
power_SM_Fundamental В примере используется основной блок SI синхронной машины и основной блок pu синхронной машины для моделирования синхронного генератора 555 МВА, 24 кВ, 60 Гц, 3600 об/мин. Здесь показано, как задать параметры SI и pu, а также как вычислить параметры обмотки ротора поля и демпфера, которые относятся к статору. Кроме обмотки поля, круглый ротор этой машины имеет три обмотки демпфера: один демпфер по прямой оси и два демпфера по квадратурной оси.
Три схемы имитируют одну и ту же синхронную машину:
Схема 1: Основные параметры задаются в СИ. Номинальный ток поля указан (ifn = 1300 A).
Схема 2: Основные параметры задаются в SI. Номинальный ток поля не указан (ifn = 0).
Схема 3: Основные параметры указаны в pu.
Параметры машины взяты из примера в Kundur [5].
Ознакомьтесь с разделом обратного вызова Model Properties/PreLoad Fcn модели, чтобы увидеть спецификации машины и расчет статора и полевых оснований, параметров ротора RL, относящихся к статору, коэффициента преобразования и номинального напряжения и тока поля.
Машины первоначально работают в установившемся состоянии практически без нагрузки (нагрузка = 0,1% номинальной мощности) при постоянном напряжении поля и механической мощности. Двухфазный шестицикльный отказ применяется при t = 0,1 с. Блок Scope показывает сравнение напряжения AB линии к линии, тока статора фазы A и тока поля трех машин.
Для моделирования дискретной модели в блоке powergui задайте для параметра Тип моделирования значение Discrete. Модель дискретизируют со временем выборки Ts = 50 мкс. Чтобы получить стабильную модель с такой небольшой нагрузкой (0,1% номинальной мощности), на вкладке «Дополнительно» каждого блока синхронной машины установите значение «Дискретная модель решателя». Trapezoidal iterative (alg. loop).
power_syncmachine В примере используется блок Синхронная машина (Synchronous Machine) в режиме двигателя. Моделируемая система состоит из синхронного двигателя мощностью 150 л.с. (112 кВА), 762 В промышленного класса, подключенного к сети с уровнем короткого замыкания 10 МВА. Машина инициализируется для выходной электрической мощности − 50 кВт (отрицательное значение для моторного режима), соответствующей механической мощности − 48,9 кВт. Соответствующие значения механической мощности и напряжения поля определяются блоком Pm Step и блоком Vf Constant. Блок Pm Step применяет внезапное увеличение механической мощности с − 48,9 кВт до − 60 кВт в момент времени t = 0,1 с.
Запустите моделирование.
После увеличения нагрузки с 48,9 кВт до 60 кВт при t = 0,1 с скорость машины колеблется перед стабилизацией до 1800 об/мин. Угол нагрузки (угол между напряжением клеммы и внутренним напряжением) увеличивается от − 21 градуса до − 53 градуса.
[1] Канай, И.М. «Причины расхождений при расчете количества ротора и точных эквивалентных диаграмм синхронной машины». Транзакции IEEE ® по устройствам и системам питания. PAS-88, № 7 (1969): 1114-1120.
[2] Moeini, A., et al. «Синхронная модель Machine Stability, Обновление Метода Перевода Данных о Станд. IEEE 1110-2002». Сеансы панели стандартов IEEE. 2018.
[3] Руководство IEEE по методам моделирования синхронных генераторов и их применению при анализе стабильности энергосистемы. IEEE Std 1110-2002 (Редакция IEEE Std 1110-1991 [2003]): 1-72.
[4] Краузе, P.C. Анализ электрических машин. Раздел 12.5. Нью-Йорк: Макгроу-Хилл, 1986.
[5] Kundur, P. Стабильность и управление энергосистемой. Нью-Йорк, Макгроу-Хилл, 1994.
Система возбуждения | Гидравлическая турбина и регулятор | powergui | Упрощенная синхронная машина | Паровая турбина и регулятор | Синхронная машина pu Основные | Синхронная машина pu Standard