Synchronous Machine pu Standard

Моделируйте динамику трехфазной синхронной машины с неявнополюсным ротором или линейным полюсом с использованием стандартных параметров в единицах измерения pu

  • Библиотека:
  • Simscape / Электрический / Специализированные Энергосистемы / Электрические Машины

  • Synchronous Machine pu Standard block

Описание

Блок Synchronous Machine pu Standard моделирует синхронную машину в генераторном или моторном режиме, используя стандартные параметры в единицах измерения pu. Режим работы диктует знак механической степени (положительный для режима генератора или отрицательный для режима мотора). Электрическая часть машины представлена моделью пространства состояний шестого порядка, и механическая часть та же, что и в блоке Simplified Synchronous Machine.

Модель учитывает динамику обмоток статора, поля и демпфера. Эквивалентная схема модели представлена в исходной системе координат ротора (qd система координат). Обмотки статора соединены в соединениях с внутренней нейтральной точкой. Все параметры ротора и электрические величины просматриваются из статора и идентифицируются грунтованными переменными. Нижние индексы:

  • d,q - d - и q-составляющая величина

  • R,s - Количество ротора и статора

  • l,m - Индуктивность утечек и намагничивания

  • f,k - количество обмотки возбуждения и демпфера

Электрическая модель машины показана на этих схемах.

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

Обычная теория синхронного машинного моделирования для анализа устойчивости принимает, что взаимные индуктивности между якорем, демпфером и полем на обмотках с прямой осью идентичны. Обычно обмотки демпфера находятся вблизи воздушного зазора, и в результате цепи демпфера, связывающие поток, почти равны якорю, связывающему поток. Эта гипотеза дает приемлемые результаты для широкой области значений исследований стабильности, особенно на стороне сети. Однако, когда дело доходит до полевых текущих исследований, существует значительная ошибка. Эквивалентная схемная динамическая модель синхронной машины может включать дополнительную индуктивность, представляющую различие между демпфером возбуждения и взаимными индуктивностями якоря возбуждения на оси D [1]. Эта индуктивность обычно называется Canay inductance. Индуктивность каната соответствует потоку утечек, ΦC, на следующем рисунке и интерпретируется как корректирующий элемент в эквивалентной модели, которая может иметь отрицательное значение [2].

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

Релевантными уравнениями являются:

Vd=idRsωψq+dψddtVq=iqRs+ωψd+dψqdtV0=i0R0+dψ0dtVfd=dψfddt+Rfdifd0=dψkddt+Rkdikd0=dψkq1dt+Rkq1ikq10=dψkq2dt+Rkq2ikq2[ψdψkdψfd]=[L+mdLlLmdLmdLmdLlkd+ Lf1d+ LmdLf1d+ LmdLmdLf1d+ LmdLlfd+ Lf1d+ Lmd][idikdifd][ψqψkq1ψkq2]=[Lmq+LlLmqLmqLmqLmq+Lkq1LmqLmqLmqLmq+Lkq2][iqikq1ikq2]

Допущения и ограничения

В дискретных системах, когда вы устанавливаете параметр Discrete solver model блока Synchronous Machine равным Trapezoidal non iterativeВозможно, вам придется подключить небольшую паразитную резистивную нагрузку к клеммам машины, чтобы избежать численных колебаний. Большие шаги расчета требуют больших нагрузок. Минимальная резистивная нагрузка пропорциональна шагу расчета. Как правило, помните, что с временным шагом 25 мкс на системе 60 Гц минимальная нагрузка составляет примерно 2,5% номинальной степени машины. Для примера синхронная машина 200 MVA в степень системе, дискретизированной с шагом расчета 50 мкс, требует приблизительно 5% сопротивления или 10 МВт. Если шаг расчета уменьшается до 20 мкс, должна быть достаточной резистивная нагрузка 4 МВт.

Однако, если вы задаете значение параметра Discrete solver model блока Synchronous Machine Trapezoidal iterative (alg. loop), можно использовать незначительную паразитную нагрузку (ниже 0,1% номинальной степени) при сохранении численной устойчивости. Эта итерационная модель создает алгебраический цикл и приводит к более медленной скорости симуляции.

Порты

Вход

расширить все

Механическая степень на валу машины, в пу. В режиме генератора этот вход может быть положительной константой или функцией или выходом блока первичного двигателя (см. блоки Hydraulic Turbine and Governor или Steam Turbine and Governor). В моторном режиме этот вход обычно является отрицательной константой или функцией.

Зависимости

Чтобы включить этот порт, установите параметр Mechanical input (Configuration вкладка) равным Mechanical power Pm.

Скорость машины, в рад/с.

Зависимости

Чтобы включить этот порт, установите параметр Mechanical input (Configuration вкладка) равным Speed w.

Напряжение возбуждения, в пу. Это напряжение может быть подано регулятором напряжения в режиме генератора (см. Excitation System блок). В моторном режиме это значение этого входа обычно является константой.

Выход

расширить все

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

Имя

Определение

Модули

МСФО

Ток статора is_a

pu

ibs

Ток статора is_b

pu

ics

Ток статора is_c

pu

IQ

Ток статора iq

pu

я бы

Идентификатор тока статора

pu

IFD

Ток поля ifd

pu

ikq1

Ток обмотки демпфера ikq1

pu

Ikq2

Ток обмотки демпфера ikq2

pu

ikd

Ток обмотки демпфера ikd

pu

phimq

Взаимный fimq потока

pu

phimd

Взаимный Фимд Потока

pu

vq

Напряжение статора vq

pu

vd

Напряжение статора vd

pu

lmq

Lmq насыщенная индуктивность

pu

lmd

Lmd насыщенная индуктивность

pu

dtheta

Отклонение угла ротора d_theta

рад

w

Скорость ротора wm

pu

Pe

Электрическая степень Pe

pu

собственный вес

Отклонение скорости ротора dw

pu

theta

Механический угол theta ротора

степень

Те

Электромагнитный крутящий момент Te

pu

дельта

Дельта угла нагрузки

степень

Pe0

Выходная активная степень Peo

pu

Qe0

Выходная реактивная степень Qeo

pu

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

Сохранение

расширить все

Специализированный электрический порт сопоставлен с электрическим выводом фазы А.

Специализированный электрический порт сопоставлен с электрическим выводом фазы B.

Специализированный электрический порт сопоставлен с электрическим выводом фазы С.

Механический вращательный порт сопоставлен с валом ротора машины.

Зависимости

Чтобы включить этот порт, установите параметр Mechanical input (Configuration вкладка) равным Mechanical rotational port.

Параметры

расширить все

Задайте стандартные параметры, также известные как рабочие параметры синхронной машины (такие как установившееся состояние, переходные и/или вычитающие реактивные напряжения и временные константы) в качестве входов для блока. Как правило, производители машин обеспечивают рабочие параметры.

В программном обеспечении анализа степеней уравнения синхронной машины обычно решаются с помощью метода преобразования прямой квадратуры и нуля. Метод преобразования данных вычисляет фундаментальные параметры из стандартных параметров [2], [3].

Строение

Набор предопределенных электрических и механических параметров для различных номиналов синхронной машины степени (кВА), напряжения преобразования фазы в фазу (V), частоты (Гц) и номинальной скорости (об/мин).

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

  • 01: 50Hz 400V 8.1kVA 1500RPM

  • 02: 50Hz 400V 16kVA 1500RPM

  • 03: 50Hz 400V 31.3kVA 1500RPM

  • 04: 50Hz 400V 42.5kVA 1500RPM

  • 05: 50Hz 400V 60kVA 1500RPM

  • 06: 50Hz 400V 85kVA 1500RPM

  • 07: 50Hz 400V 250kVA 1500RPM

  • 08: 50Hz 400V 325kVA 1500RPM

  • 09: 50Hz 400V 670kVA 1500RPM

  • 10: 50Hz 400V 910kVA 1500RPM

  • 11: 50Hz 400V 1320kVA 1500RPM

  • 12: 50Hz 400V 2000kVA 1500RPM

  • 13: 60Hz 460V 10.2kVA 1800RPM

  • 14: 60Hz 460V 20kVA 1800RPM

  • 15: 60Hz 460V 37.5kVA 1800RPM

  • 16: 60Hz 460V 52.5kVA 1800RPM

  • 17: 60Hz 460V 72.5kVA 1800RPM

  • 18: 60Hz 460V 100kVA 1800RPM

  • 19: 60Hz 460V 300kVA 1800RPM

  • 20: 60Hz 460V 406kVA 1800RPM

  • 21: 60Hz 460V 800kVA 1800RPM

  • 22: 60Hz 460V 1075kVA 1800RPM

  • 23: 60Hz 460V 1588kVA 1800RPM

  • 24: 60Hz 460V 2500kVA 1800RPM

Выберите No (по умолчанию) если вы не хотите использовать предустановленную модель или хотите изменить некоторые параметры предустановленной модели, как описано ниже.

При выборе предустановленной модели электрические и механические параметры на вкладке Parameters затемняются. Чтобы начать с предустановленной модели и затем изменить параметры машины:

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

  2. Измените параметр Preset model на No. Это действие не изменяет параметры машины, но разрывает связь с предустановленной моделью.

  3. Измените параметры машины так, как вы хотите.

Представляет ли механическая степень, приложенная к валу, или скорость ротора как Simulink® вход блока или для представления вала машины с Simscape™ вращательным механическим портом.

Выберите Mechanical power Pm чтобы задать механический вход степени, в pu, и показать порт Pm. Скорость машины определяется постоянной инерцией H и различием между механическим крутящим моментом Tm, который является результатом приложенной механической Pm степени, и внутренним электромагнитным крутящим моментом Te. Когда скорость положительная, положительный механический сигнал степени указывает на режим генератора, а отрицательный сигнал указывает на режим мотора.

Выберите Speed w чтобы задать вход в pu и открыть w порт. Скорость машины налагается, и механическая часть модели (постоянное H инерции) игнорируется. Использование скорости в качестве механического входа позволяет вам смоделировать механическую связь между двумя машинами.

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

Выход скорости машины 1 (двигатель) соединяется со входом скорости машины 2 (генератор). На этом рисунке крутящий момент трения проигнорирован в машине 2. Поэтому выходной Te электромагнитного крутящего момента соответствует Tm механического крутящего момента, приложенному к валу машины 1. Соответствующая механическая входная степень машины 1 вычисляется как Pm = Tm * w. Коэффициент Kw учитывает модули скорости обеих машин (pu) и отношение коробки передач w2/w1. Коэффициент KT учитывает модули крутящего момента обеих машин (pu) и номинальные значения машины. Кроме того, поскольку J2 инерции проигнорирована в машине 2, J2 относится к скорости машины 1 и должно быть добавлено к инерции машины 1, J1.

Выберите Mechanical rotational port S, чтобы открыть механический вращательный порт Simscape, который позволяет вам соединить вал машины с другим валом машины или с другими блоками Simscape, которые имеют механические вращательные порты.

Рисунок показывает, как соединить блок Ideal Torque Source из библиотеки Simscape с портом вала машины, чтобы представлять машину в режиме двигателя или режиме генератора, когда скорость ротора положительная.

Задайте тип ротора следующим образом Salient-pole или Round (цилиндрический). Настройка влияет на количество схем ротора в q-оси (обмотки демпфера).

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

Когда этот флажок снят, выходы измерения используют определение сигнала для идентификации меток шины. Метки содержат неалфавитно-цифровые символы, которые несовместимы с некоторыми приложениями Simulink.

Параметры

Общая трехфазная видимая степень (VA), линейное напряжение RMS (pu), частота (Гц) и ток возбуждения (pu).

Этот параметр идентичен параметру Nominal power, voltage, frequency, field current [ Pn(VA) Vn(Vrms) fn(Hz) ifn(A) ] блока Synchronous Machine SI Fundamental, кроме того, что вы не задаете номинальный ток поля. Это значение не требуется, потому что блоку Synchronous Machine pu Standard не нужен коэффициент преобразования. Поскольку величины ротора просматриваются из статора, они преобразуются в pu с помощью базовых величин статора, выведенных из предыдущих трех номинальных параметров.

d ось синхронное реактивное Xd, переходное реактивное Xd 'и вычитающее реактивное Xd "; q -ось синхронного реактивного Xq и вычитающего реактивного Xq "; и реактивное Xl утечек (все в pu).

Зависимости

Чтобы включить этот параметр, установите Rotor type равным Salient-pole.

d ось синхронное реактивное Xd, переходное реактивное Xd 'и вычитающее реактивное Xd "; q - синхронное реактивное Xq, переходное реактивное Xq 'и вычитающее реактивное Xq "; и реактивное Xl утечек (все в pu).

Зависимости

Чтобы включить этот параметр, установите Rotor type равным Round.

Временная константа для оси d.

Временная константа для оси q.

Временные константы все в s для следующего:

  • d разомкнутой цепи с переходной осью (Tdo ')

  • d -ось, вычитающая разомкнутую схему (Tdo ")

  • q -ось, вычитающая разомкнутую схему (Tqo ")

Зависимости

Чтобы включить этот параметр, установите Rotor type равным Salient-pole, d axis к Open-circuit, и q axis к Open-circuit.

Временные константы все в s для следующего:

  • d разомкнутой цепи с переходной осью (Tdo ')

  • d -ось, вычитающая разомкнутую схему (Tdo ")

  • q разомкнутый контур с переходной осью (Tqo ')

  • q -ось, вычитающая разомкнутую схему (Tqo ")

Зависимости

Чтобы включить этот параметр, установите Rotor type равным Round, d axis к Open-circuit, и q axis к Open-circuit.

Временные константы все в s для следующего:

  • d разомкнутой цепи с переходной осью (Tdo ')

  • d -ось, вычитающая разомкнутую схему (Tdo ")

  • q -ось, вычитающая короткую схему (Tq ")

Зависимости

Чтобы включить этот параметр, установите Rotor type равным Salient-pole, d axis к Open-circuit, и q axis к Short-circuit.

Временные константы все в s для следующего:

  • d разомкнутой цепи с переходной осью (Tdo ')

  • d -ось, вычитающая разомкнутую схему (Tdo ")

  • q короткая короткая схема с переходной осью (Tq ')

  • q -ось, вычитающая короткую схему (Tq ")

Зависимости

Чтобы включить этот параметр, установите Rotor type равным Round, d axis к Open-circuit, и q axis к Short-circuit.

Временные константы все в s для следующего:

  • d -ось переходной короткой схемы (Td ')

  • d -ось вычитает короткую схему (Td ")

  • q -ось, вычитающая разомкнутую схему (Tqo ")

Зависимости

Чтобы включить этот параметр, установите Rotor type равным Salient-pole, d axis к Short-circuit, и q axis к Open-circuit.

Временные константы все в s для следующего:

  • d -ось переходной короткой схемы (Td ')

  • d -ось вычитает короткую схему (Td ")

  • q разомкнутый контур с переходной осью (Tqo ')

  • q -ось, вычитающая разомкнутую схему (Tqo ")

Зависимости

Чтобы включить этот параметр, установите Rotor type равным Round, d axis к Short-circuit, и q axis к Open-circuit.

Временные константы все в s для следующего:

  • d -ось переходной короткой схемы (Td ')

  • d -ось вычитает короткую схему (Td ")

  • q -ось, вычитающая короткую схему (Tq ")

Зависимости

Чтобы включить этот параметр, установите Rotor type равным Salient-pole, d axis к Short-circuit, и q axis к Short-circuit.

Временные константы все в s для следующего:

  • d -ось переходной короткой схемы (Td ')

  • d -ось вычитает короткую схему (Td ")

  • q короткая короткая схема с переходной осью (Tq ')

  • q -ось, вычитающая короткую схему (Tq ")

Зависимости

Чтобы включить этот параметр, установите Rotor type равным Round, d axis к Short-circuit, и q axis к Short-circuit.

Сопротивление статора Rs, в pu.

Постоянные H инерции, которые являются отношением энергии, сохраненной в роторе с номинальной скоростью по номинальной степени машины, коэффициентом трения F (крутящий момент pu/скорость pu) и количеством p пар полюсов. Крутящий момент трения Tf пропорциональен скорости вращения ротора ω (Tf = F., где все величины выражены в pu).

Зависимости

Чтобы включить этот параметр, установите Mechanical input равным Mechanical power Pm или Mechanical rotational port.

Количество пар полюсов синхронной машины.

Зависимости

Чтобы включить этот параметр, установите Mechanical input равным Speed w.

Начальное отклонение скорости Δω (процент номинальной скорости), электрический угол Θe ротора (градусы), линейные переменные тока ia, ib, ic (pu) и углы фазы pha, phb, phc (градусы) и начальное напряжение возбуждения Vf (pu). Можно вычислить эти значения автоматически с помощью инструмента Load Flow или инструмента Machine Initialization блока powergui.

Моделируется ли магнитное насыщение ротора и железа статора.

Установите этот флажок, чтобы предоставить матрицу параметров для симуляции насыщения.

Снимите этот флажок, чтобы не моделировать насыщение в симуляции. В этом случае отношение между ifd и Vt линейное (без насыщения).

Параметры кривой насыщения без нагрузки. Магнитное насыщение статора и роторного железа моделируется кусочно-линейной зависимостью, задающей точки на кривой насыщения без нагрузки. Первая строка этой матрицы содержит значения токов поля (в pu). Вторая строка содержит значения соответствующих контактных напряжений (в pu). Первая точка (первый столбец матрицы) должна отличаться от [0,0]. Эта точка соответствует точке, где начинается эффект насыщения. Для насыщения номинальный ток возбуждения и номинальное линейное напряжение RMS являются базовыми значениями для тока возбуждения и контактного напряжения, соответственно.

Щелкните Plot, чтобы просмотреть кривую насыщения без нагрузки.

Зависимости

Чтобы включить этот параметр, выберите Simulate saturation.

Расширенный

Чтобы включить вкладку Advanced, в блоке powergui задайте Simulation type Discrete и очистить Automatically handle discrete solver.

Шаг расчета, используемого блоком. Чтобы наследовать шаг расчета, заданную в блоке powergui, установите этот параметр равным −1.

Метод интегрирования, используемый блоком, когда параметр Solver type блока powergui равен Discrete.

Значение Discrete solver model автоматически устанавливается на Trapezoidal robust при выборе параметра Automatically handle Discrete solver and Advanced tab solver settings of blocks блока powergui.

The Trapezoidal non iterative и Trapezoidal iterative (alg. loop) методы больше не рекомендуются для дискретизации блока Synchronous Machine pu Standard. Trapezoidal non iterative требует, чтобы вы добавили незначительные нагрузки шунта на терминалах машины, чтобы гарантировать стабильность симуляции, и Trapezoidal iterative (alg. loop) может не сходиться и привести к остановке симуляции, когда количество машин увеличивается в модели.

The Trapezoidal robust и Backward Euler robust методы позволяют устранить необходимость использования паразитарных нагрузок и моделировать машину без нагрузок. Для устранения топологических ошибок машин, подключенных к индуктивной схеме (для примера, выключатель, соединенный последовательно с машиной); машина моделирует незначительную внутреннюю нагрузку в 0,01% от номинальной степени.

The Trapezoidal robust метод немного точнее, чем Backward Euler robust метод, особенно когда модель моделируется при больших шагах расчета. The Trapezoidal robust способ может привести к легким демпфированным численным колебаниям на напряжении машины при условиях отсутствия нагрузки, в то время как Backward Euler robust способ предотвращает колебания и поддерживает хорошую точность.

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

Поток нагрузки

Параметры потока нагрузки используются, чтобы задать параметры блоков для использования с инструментом Load Flow блока powergui. Эти параметры потока нагрузки используются только для инициализации модели. Они не имеют никакого влияния ни на модель блока, ни на эффективность симуляции.

Тип генератора машины. Выбор:

  • swing - реализует генератор, управляющий величиной и углом фазы своего терминального напряжения. Ссылки величины напряжения и угол заданы параметрами Swing bus or PV bus voltage и Swing bus voltage angle блока Load Flow Bus, соединенного с клеммами машины.

  • PV - реализует генератор, управляющий его выходной активной P степени и V величины напряжения. P задается параметром Active power generation P (W) блока. V задается параметром Swing bus or PV bus voltage блока Load Flow Bus, соединенного с клеммами машины. Можно управлять минимальной и максимальной реактивной степенью, генерируемой блоком, используя параметры Minimum reactive power Qmin (var) и Maximum reactive power Qmax (var).

  • PQ - Реализует генератор, управляющий его выходной активной P степени и реактивной Q степени. P и Q заданы параметрами Active power generation P (W) и Reactive power generation Q (var) блока, соответственно.

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

Зависимости

Чтобы включить этот параметр, установите Generator type равным PV или PQ.

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

Зависимости

Чтобы включить этот параметр, установите Generator type равным PQ.

Минимальная реактивная степень, которое может быть сгенерировано машиной при сохранении напряжения на клемме на ссылку значении. Это опорное напряжение задается параметром Swing bus or PV bus voltage блока Load Flow Bus, соединенного с клеммами машины. Значение по умолчанию -inf, что означает отсутствие нижнего предела на выходной реактивной степени.

Зависимости

Чтобы включить этот параметр, установите Generator type равным PV.

Максимальная реактивная степень, которое может быть сгенерировано машиной при сохранении напряжения на ссылку значении. Это опорное напряжение задается параметром Swing bus or PV bus voltage блока Load Flow Bus, соединенного с клеммами машины. Значение по умолчанию inf, что означает отсутствие верхнего предела на выходной реактивной степени.

Зависимости

Чтобы включить этот параметр, установите Generator type равным PV.

Примеры

The power_thermal пример использует блок Synchronous Machine pu Standard для моделирования субсинхронного резонанса (SSR) в паровой турбине и регуляторе в сети с последовательной компенсацией.

Ссылки

[1] Канай, И.М. «Причины расхождений в вычислении количества роторов и точных эквивалентных схем синхронной машины». IEEE®Транзакции на энергетических аппаратах и системах. PAS-88, № 7 (1969): 1114-1120.

[2] Moeini, A., et al. «Синхронная модель устойчивости машины, обновление до IEEE Std 1110-2002 Data Translation Technique». Сеансы панели стандартов IEEE. 2018.

[3] Руководство IEEE по методикам моделирования синхронных генераторов и применениям в анализе устойчивости системы степеней. IEEE Std 1110-2002 (Редакция IEEE Std 1110-1991 [2003]): 1-72.

[4] Krause, P.C. Анализ электрической техники. Раздел 12.5. Нью-Йорк: McGraw-Hill, 1986.

[5] Кундур, П. Устойчивость системы Степени и Управление. Нью-Йорк, McGraw-Hill, 1994.

Расширенные возможности

Генерация кода C/C + +
Сгенерируйте код C и C++ с помощью Coder™ Simulink ®

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