Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Математическое моделирование газотурбинных мини-электростанций и мини-энергосистем

..pdf
Скачиваний:
7
Добавлен:
12.11.2023
Размер:
12.91 Mб
Скачать

ние относительных единиц дает некоторые преимущества

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

вотносительных единицах.

Все величины, характеризующие режим работы энерго­ системы в системе относительных единиц, выражаются в до­ лях некоторых величин (той же размерности), принятых за единицы величин, называемых базисными (базовыми). В ка­ честве базисных величин выбираются номинальные напря­ жение и мощность [85]: C/e = U„0№Sç =S„0M.

Таким образом, при выбранных базисных условиях от­ носительные значения ЭДС, напряжения, тока, мощности и сопротивления следующие:

Е* =E/U6; U* =U/U6; / ’ = / / / 6; 5* =S/S6; z* =z/z6.

Для номинального режима в относительных единицах при линейных токах и напряжениях имеем:

^ном — Аюм — ^ком — Ь ^ном ~ Фном ’ Ô HOM — Фном ’

**

где P , Q - активная и реактивная мощности в относи­

тельных единицах.

В принятой системе относительных единиц базисное сопротивление определяется выражением z6=U6/S6, базис­

ный ток

/б = 5б l4bUb. Таким образом, базовый ток / 6 = / н,

так как

5Н=>/3-С/Н-/Н, и

56 = >/з-t/ 6 -76. При этом для ак­

тивного,

реактивного и

полного сопротивлений r * = r/z 6;

x*=x/z6; z* = z /z 6.

Частоту вращения ю, с-1, относят к синхронной номи­ нальной Шб=о»ном, которая при/ = 50 Гц соответствует 314 с-1.

В качестве базисного момента (электромагнитного и ме­ ханического) рассматривается такой момент, который при синхронной номинальной частоте вращения соответствует активной мощности, равной S6: Мб= S6/со„ом .

В относительных единицах справедливо равенство

^ном —^ном —Софтом- В качестве базисного времени принимается время, в те­

чение которого ротор при синхронной номинальной частоте вращения поворачивается на угол 1 рад: / 6 = 1 /шном.

При/ = 50 Гц t\ - t/t6=314/, аналогично пересчитыва­ ются в относительные единицы все постоянные времени: t -

в секундах, /j - в радианах.

Относительные единицы для обмотки возбуждения син­ хронных машин записываются с условием их приведения к статорным цепям [163]. При этом принимается:

базовый ток обмотки возбуждения

I/S=I/0Xad'y

базовое напряжение обмотки возбуждения

базовое сопротивление обмотки возбуждения

где Iff, - ток возбуждения холостого хода синхронной ма­ шины.

Вращающиеся электрические машины. Под вращаю­ щимися электрическими машинами здесь понимаются:

• синхронные машины (генераторы и двигатели);

• асинхронные машины (двигатели).

Вопросам математического моделирования электриче­ ских машин переменного тока посвящено много работ [25, 30, 71, 77, 99 и др.]. К настоящему времени предложено и ап­ робировано большое число математических моделей, описы­ вающих поведение электрических машин в различных режи­ мах с той или иной степенью полноты и достоверности. Эти модели могут быть получены различным путем: теория це­ пей, теория поля, идентификация. Алгоритмы могут быть ориентированы на использование аналоговых [147] и цифро­ вых [102, 164] вычислительных машин, отличаться по спосо­ бу использования и назначению. По этой причине выбор ма­ тематического описания нуждается в обосновании.

Известно [85], что выбор математической модели, преж­ де всего, определяется той задачей, для решения которой вы­ полняется моделирование. Исходя из задачи создания мате­ матической модели мини-ЭЭС, необходимо выбрать такую математическую модель электрической машины, которая по­ зволит имитировать все основные режимы, существенно влияющие на поведение мини-ЭЭС. Из всей совокупности нормальных и аварийных режимов электрических машин можно выделить группу характерных режимов, переходные процессы в которых являются определяющими как для самой машины, так и для узла нагрузки, к которому она подключе­ на. Перечислим эти режимы: внезапное короткое замыкание; внезапный наброс нагрузки на двигатель; пуск и самозапуск асинхронных двигателей; выбег после отключения двигателя от сети; внезапное отключение нагрузки; синхронизация ге­ нераторов. Кроме того, необходимо исследовать влияние ре­ гулирования возбуждения синхронных двигателей и син­ хронных генераторов.

Начнем рассмотрение с математического описания син­ хронных машин, поскольку, как справедливо отмечается в работе [62], для электрических систем главными являются

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

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

Исходные допущения. При выводе дифференциальных уравнений применяются следующие общеизвестные допуще­ ния [85]:

а) статор имеет трехфазную симметричную обмотку; б) распределение магнитного поля каждой обмотки

вдоль окружности воздушного зазора синусоидально; в) отсутствуют потери в стали;

г) наличие пазовой неравномерности воздушного зазора не влияет на магнитную проводимость.

Уравнения синхронной машины.

1. Уравнения в преобразованных координатах d. а. В на­ стоящее время для моделирования синхронных машин обще­ принятыми являются полные или упрощенные уравнения Парка-Горева [25], эти уравнения составляются для коор­ динат d,q, вращающихся синхронно с ротором машины. Такое представление помогает освободиться от переменных коэффициентов, являющихся функциями от углового поло­ жения ротора. Для синхронного генератора, имеющего на ро­ торе два эквивалентных демпферных контура, полную сис­ тему уравнений Парка-Горева удобно записать в следую­ щем виде:

Ud =-'Vq<ù-d '¥ d/d t- idr;

Uq ='¥d(ù-d'¥q/d t-iqr; Uf= d4ff /dt + if rf \

0 = d'VD/dt + iDrD;

0 = dxVQ/dt + iQrQ\

(2.2)

d(ù/dt = j ( M T-M);

dy/dt = (ù; d5/dt = (o-(û0; M ^ diq- ^ qid,

где ©о - угловая частота вращения поля статора; со - угловая частота вращения ротора; ô - внутренний угол машины (угол нагрузки); у - угол поворота оси d по отношению к непод­ вижной оси А (см. рис. 2.3); Мт- момент турбины; М - элек­

тромагнитный момент генератора;

'F, - потокосцепления

статора по продольной и поперечной осям; 'Vu,

- потокос­

цепления демпферных контуров по продольной и поперечной осям; id, iq - токи статора по продольной и поперечной осям; /о, ÎQ - токи демпферных контуров по продольной и попереч­ ной осям; г - активное сопротивление обмотки статора; rD, rQ- активные сопротивления демпферных контуров по про­ дольной и поперечной осям; 'Vf, if, rf - потокосцепление, ток и сопротивление обмотки возбуждения; J - момент инерции электрической машины; Ц/, Uq- напряжения на обмотке ста­ тора по продольной и поперечной осям.

Уравнения записаны для следующего взаимного распо­ ложения координатных осей (рис. 2.3). Это обстоятельство специально подчеркивается, так как в литературе существует 128 вариантов записи уравнений (2 .2 ), отвечающей различ­ ным комбинациям направлений d, q, со, у [25].

d

а

P

Рис. 2.3. Соотношение между координатами (А, В, C),d,q и а, р

Уравнения (2.2) полностью характеризуют переходный процесс в том случае, если напряжение Ù на шинах не зави­ сит от режима системы (шины неизменного напряжения). Ес­ ли это напряжение зависит от режима системы, то необходи­ мо или составить дополнительные уравнения, выявляющие эту зависимость [25], или рассчитать это напряжение в от­ дельных программных модулях, имитирующих электриче­ скую нагрузку генератора. Система уравнений (2.2) дает связь между ЭДС, напряжениями, токами, потокосцеплениями и параметрами рассматриваемой машины. Отметим, что

всилу принятой записи уравнений (2 .2 ), матрицы прямого

иобратного преобразования координат d, q н а, Ь, с имеют следующий вид:

cosy (

^ 2 . Г = — smу

1

1

1

2

2

2

J

 

 

 

cosy smy

(

2п

. ( 2я

(2.4)

г = cos у -----

sml у — —

V

3

 

 

cos(y

4п

sin(y - ■— )

 

~ )

 

Уравнения связи электрических контуров при отсутст­ вии насыщения записываются в векторной форме:

¥ = Х1, (2.5) где X - матрица относительных индуктивных сопротивле­ ний; i и векторы тока и потокосцеплений, или:

 

| Ч

*ad

*ad

0

0

^

| V

 

 

 

 

 

 

 

Xad

XD

*ad

0

0

 

*D

 

 

 

 

 

 

=

Xad

Xad

Xf

0

0

 

X h

(2.6)

%

0

0

0

**

 

 

l4

 

 

 

0

0

* *

XQ J

KlQ j

 

 

 

 

 

 

 

 

Поскольку матрица X постоянна, то, найдя один раз об­ ратную ей матрицу X"1, умножив вектор переменного столб­ ца свободных членов Y на постоянную матрицу X”1, можно вычислить вектор i [85]:

i = X "'Y

(2.7)

Уравнения (2.2) называются «полными», так как они учитывают все основные компоненты электромагнитного процесса. Упрощенные уравнения получаются из полных за счет ряда допущений (пренебрежение демпферными конту­ рами, апериодическими процессами в обмотке статора и др.). Различные варианты уравнений подробно рассмотрены в ра­ ботах [25,96,175, 189,192-194,195].

Обычно для синхронного двигателя используются урав­ нения с теми же знаками переменных величин, что и для ге-

нераторов [см. уравнения (2 .2 )]. При этом механический мо­ мент на валу Мс (вместо Мт) принимается отрицательным. Ес­ тественно, что и угол нагрузки ô в этом случае будет отрица­ тельным (ротор отстает от магнитного поля статора).

Напряжение возбуждения является внешней переменной модели. Закон его изменения должен описываться одним из блоков программного комплекса, который реализует автома­ тический регулятор возбуждения (АРВ). Характер изменения механического момента (Мт) на валу должен описываться специальными блоками вычислительного комплекса.

Получим математическое описание синхронной машины в обобщенной векторно-матричной форме (2 .1 ).

Используя выражения (2.2) и (2.6), получаем путем под­ становки (2 .6 ) в (2 .2 ) при условии непрерывности всех токов и постоянства всех индуктивных сопротивлений:

f xd

0

Xad

Xad

0

'

 

0

 

0

0

 

 

 

 

 

 

х <ч

 

Xod

0

Xf

Xad

0

 

Xad

0

Xad

XD

0

 

I

0

 

0

0

 

,

 

 

x ô

/

 

(0

 

0

 

 

Г

 

 

 

со

xd

г

 

 

 

 

—1

 

 

 

-0» Xad

 

 

 

 

 

 

 

 

0

 

0

 

rf

 

 

 

 

 

 

 

 

0

 

0

 

0

 

 

0

 

( )

 

0

 

 

\

 

 

 

 

 

 

 

 

г

V I

( U

é

 

 

 

 

 

V ,

 

 

 

X

h

-

 

f

 

 

 

- V

 

 

 

I D

 

0

 

V /

P h

P h

P^D

A ,

0

1 Xad

0

JrD

0

\

=

и v

0

0

0

к !QJ

l о y

Поскольку нас интересуют статорные токи (Id и Iq),

выделим из системы (2 .8 ) производные статорных токов:

 

rpld

 

 

 

(ÙX„

о

сох.aq

О

яу кр1я)

~(ùX,

r

-(ùxad

0

 

 

 

 

 

 

 

 

 

f

т\

 

(2.9)

 

ad

 

\

P1/

rU,

 

*ad

 

 

0

0

 

ptD -

\ U„

 

 

V

e

J

 

 

 

 

 

 

\J Q J

Избавимся от производных роторных токов в правой части системы (2.9), для этого из системы (2.6) выразим про­ изводные потокосцеплений роторных контуров и приравняем их к правым частям системы (2 .2 ):

(

Xad

A \

f

Aad

 

V n 7 N

' p V f '

0

 

0

PJf

 

 

V

/

 

 

PlD

P ^ D

Xad

0

+ Ы XD

0

[ P V Q ) 4

0

 

V

0

XQ ) \ plQ)

 

 

' U f - R f b '

(2.10)

- J Q R Q

Следовательно, выражая производные роторных токов через производные статорных токов, получим

V

/

( 4

Xad

0 "ч-1

P^D =

Xad

XD

0

A

>

1 °

0

 

 

 

Kad

О

pJd

 

U , - Rг,I/A

 

о

 

b Ro

 

 

 

\ pI«j

~ J Q R Q

Подставляем производные роторных токов из (2.11) в (2.9):

О '

(хА

О

т„

О

О сох,

Ad

 

 

 

ад

О

х.Ч)\ Р^Ч )

~(ùxd г

- m ad

-o>xad

K*QJ

 

 

 

 

 

 

- 1

 

 

0 '

xf

Xad

0

4

+ ' Xad

Xad

 

 

 

 

 

x4)

Xad

XD 0

 

K J l 0

0

 

 

 

 

 

 

0

xQd

 

 

l

0

 

 

 

 

 

 

 

 

o '

- 1

xad

Xad

О )

xf

Xad

 

 

 

 

 

 

 

О

 

гад J

Xad

XD

0

 

0

 

 

 

 

 

 

(x

O

Xad

 

 

<

Xad

U

 

l Xaq j

0

%- V f '

~J DrD

' p l d '

^ 4 J

. (2. 12)

\ 0 0 xe, < “ V < ? >

Преобразуем систему (2.12) так, чтобы в левой части выделить вектор производных статорных токов p l d и р 1 ч:

(

7 \

 

 

 

 

 

 

f

 

 

0]

 

 

0 " x r

Xad

PU

 

f Xad Xad

XD

 

 

 

Xn

 

 

0

JC

Xad

 

 

о

X,

J

H »

 

0

Л

;

 

Я

 

aq J

l 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ч-1

----/"■

O

 

<u

 

Т \

 

-1

 

 

 

 

Xad

0 '

сох„

 

О

шх.ад

Xad

0

 

K~(ÙXd R

~(ÙXnrl

-шх,

x„„

0

 

 

ad

Kad

 

aq )

 

 

 

 

 

 

 

 

 

\ lQ)

 

 

Xf

Xad

0

■wU f - r fif \ \

 

Xad Xad

0

 

 

~lûrD

 

Xad

XD

0

,U4J

 

0 X.aq J

 

X,Q)

 

 

 

0

0

~IerQ J)

Соседние файлы в папке книги