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

книги из ГПНТБ / Болотин, Б. И. Инженерные методы расчетов устойчивости судовых автоматизированных электростанций

.pdf
Скачиваний:
19
Добавлен:
22.10.2023
Размер:
13.69 Mб
Скачать

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

 

«в =

/( //. и;

iy\ г'в)-

 

(11.66)

Линеаризуя это выражение, получим

 

 

Аив=

Ди +

д /, +

Д/у +

-^2 - At'B.

(11.67)

 

ди

d ld

d iy

dij,

 

Сравним выражение (11.67) с линеаризованным уравнением сило­ вой части регулятора напряжения, приведенным в [46],

Д«в = [7?iA — Д*'в— Д^у]. (11.68)

Сравнивая коэффициенты при одноименных членах уравнений

(11.67)

и (11.68),

получим

 

 

 

див

 

див

_див

див

(11.69)

ди

k

dld

div

d iB

 

Найдем значения коэффициентов в выражениях (11.68) и (11.69). Коэффициент, характеризующий отбор части тока, поступающего по каналу тока в канал напряжения и на

намагничивающую цепь, равен

k-

(11.70)

Рис. 11.11. Определение коэф­ фициента по данным харак­ теристики холостого хода

напряжению определяется

где реактивное сопротивление дросселя отбора и намагничивания ТФК;

хк — компаундирующее

сопротивление

системы

регулирования

напряжения.

Параметры xfl и хк

задаются

заво-

дами-изготовителями генераторов.

 

Можно принять,

что хк =

2 -з- 4 о. е.;

Хц = 5

20 о. е.,

тогда k =

0,3

0,7.

Коэффициент компаундирования по по следующей формуле

— *в тах (й + 1)

(11.71)

где мтах — максимальное напряжение генератора на холостом ходу

при отключенном корректоре; iBmax — ток возбуждения, соответст­ вующий «тах.

Значения нтах и t'Bтах могут быть непосредственно замерены при

отключенном корректоре. Ток iBmax можно также определить по ха­ рактеристике холостого хода. Зная значение напряжения генератора при отключенном корректоре umax, находим по характеристике хо­

лостого хода (рис. 11.11) соответствующий ему ток возбуждения tBmax.

60

Если учесть,

что и

шах

= 1,1 ч- 1,3

и

ном

и взять

значения

k в ин-

*

 

 

 

 

 

 

 

тервале от 0,3 до 0,7,

можно получить

 

 

 

 

 

 

R

('в т а х ( U - f - 1,7)_ =

1 |5 ^

2 о.

е.

(11.72)

 

 

(1,1 - 1 ,3 ) «ном

 

 

 

 

 

 

Коэффициент

R 1 характеризует

степень

компаундирования по

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

Для определения Д , необходимо иметь внешнюю характеристику генератора при cos ср = 0 и при отключенном корректоре напряже­ ния (рис. 11.12). Из этой характеристики при номинальном токе ста­

тора / ном = 1 о. е. находится нтах, а затем по формуле (11.73) — зна­ чение

Rl — (ив. х. х“Н Xd) №+

О ^ 2umax>

(11.73)

где ив. х. х — напряжение возбуждения

холостого хода,

выраженное

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

участку характеристики

холостого

хода

/

(рис. 11.11).

В зависимости от параметров

 

генератора

и настройки

регулятора

 

 

 

Д х ^ К б ч -З .

(11.74)

 

Коэффициент R t может быть определен из опыта короткого замыкания при отклю­ ченном корректоре (гу = 0), когда и = 0, а установившийся ток короткого замыка­

ния / к з равен

номинальному

току гене­

ратора

/ ном.

В

этом случае

выражение

(11.68)

примет вид

 

 

uB=

i - ( / ? 1/ , - i B).

(11.75)

 

 

К

 

0

1 I о.е

Рис. 11.12. Определение коэффициента Rx по внеш­ ней характеристике генера­ тора при отключенном кор­ ректоре

Учитывая, что в установившемся режиме гв = «в (в относительных единицах), a Id = / ном, получим

R __ ив О ~Ь к)

(11.76)

I ном

 

Определим далее значение тока управления корректора напряже­ ния. Ток гу есть функция напряжения генератора и (в случае парал­ лельной работы двух ГА) разности реактивных мощностей Qx и Q2 параллельно работающих генераторов.

Таким образом,

»у= П и; (Qi - Q 2)].

(11.77)

Линеаризуя это выражение, будем иметь

Aiy ^У-Ли-1------- ^ ----- A (Qi - Q 2).

(11.78)

ди

д (Qx— Q 2)

 

Подставив это выражение в формулы (II.67), получим для одного (например, первого) из двух параллельно работающих генераторов

д и в1 А,,

д и в1

А / I

5«В1- A/

-, d u B1

■ ^ Л и +

LJlIv

d ld i

 

^Bl

а *в1

 

ди

 

 

d i y i

ди

d u s l

d iy i

A(Qx

-Q>).

(11.79)

 

 

d ( Q i

d

i y i

Q2)

 

 

 

Определим частные

производные-^-

и ----- ^ ------, а также выра­

 

 

 

зи

д (Ql — Q2)

 

зим Д (Qx—Q2) через принятые независимые переменные, остальные частные производные этого уравнения определяются соотношениями

(11.69). Производная

является передаточной функцией коррек-

 

ди

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

Конкретный характер отставания тока управления определяется типом системы управления. Для системы управления с корректором напряжения и дросселем отбора эта зависимость может быть представ­

лена следующим образом:

 

 

diy

w K(P)

 

(II.80)

ди

(1 +

 

Тур) (1 -j- Ткр)

Коэффициент усиления

корректора

kK зависит от индивидуальной

настройки корректора, он может быть определен экспериментально путем снятия зависимости iy = f (и) непосредственно на объекте.

Постоянные времени дросселя отбора и усилителя корректора Ту и Тк могут также быть получены экспериментально (например, Ту можно определить путем осциллографирования iy при скачкообразном приложении сигнала к обмотке управления дросселя отбора, Т к — путем осциллографирования выходного напряжения корректора при скачкообразном изменении напряжения и). Частная производная

——— является общей передаточной функцией корректора напря-

d(Qi—Q2)

жения при наличии уравнительных связей по реактивным нагруз­ кам, осуществляемых с помощью блока параллельной работы (БПР), и представляет собой зависимость тока управления корректора iy от разности реактивных мощностей. Для системы управления с корре­ ктором напряжения эта зависимость

div

М б . п. р

(11.81)

д (Qi — Q2) = ^ к (Р) К . п. р =

(1 + ТуР) (1 + Т кр)

 

где k6_п. р — коэффициент передачи БПР.

Произведение kK /гб п р может быть определено экспериментально замером тока управления /у при искусственном разведении реактив­ ных нагрузок путем изменений уставки по напряжению одного из ге­ нераторов.

62

Подставив значения частных производных из выражений (II.69; 11.80; 11.81) в (11.79), получим

 

AuBl = - ^ - A u

+ - ^ A I d-

 

 

 

k

k

 

 

IAu + iB. выА П . рЛ (Qx-Q*)]

(Р),

(Н*82)

Рис. 11.13. Структурная схема регулирования возбуждения при параллельной работе двух ГА

где коэффициенты R 1X и /?21 идентичны коэффициентам R х и R 2 (вто­ рые индексы указывают на принадлежность к первому генератору).

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

Определим приращение разности реактивных мощностей

A (QiQ2) = -glg i - e j ... АЕ 01

+

aF

'Q1”1

ocqx

 

 

д (Qi

Q2)

Д612-

^ ---- —д^еQ. 2 + otq2

(11.83)

63

Разность реактивных мощностей может быть выражена через не­ зависимые переменные следующим образом:

Qi— Q2

— cos а п ------

cos (б 1а— a i2)

 

 

 

 

 

zn

 

г 12

 

 

 

 

 

 

 

 

-Q2 cos а»

 

Eq\Eq2■c o s ( 6 12 + a 12) .

 

 

Тогда

 

 

 

 

Z 19

 

 

 

 

 

 

 

 

2ЕQ1

 

 

 

 

 

 

 

 

д ( Q i

Q 2)

cos a u —

cos

(6 l2 — a 12) +

 

 

 

Qi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ - ^ - c o s ( 6 12 + a 12);

 

 

 

(11.84)

 

 

 

zi2

 

 

 

 

 

 

 

 

а ( я р ~ 9а)

•=

Z22

cos a 22 ------7 — C0S (6X2

“ l 2) + 4 ^

cos

(612 + “ 12);

O E q 2

 

 

 

Z12

 

 

Z12

 

 

 

a(Q' r

Q2)

=

sin (6 1 2 - a 12) —

 

sin (fi12 + a 12).

 

ao12

 

 

zi2

 

 

 

zi2

 

 

 

Связь между отдельными составляющими уравнений регуляторов

напряжения

синхронных генераторов с

каналами амплитудно-фазо­

 

 

 

 

 

вого компаундирования (АФК) и кор­

 

 

 

 

 

ректорами напряжения при их па­

 

 

 

 

 

раллельной

работе

иллюстрируется

 

 

 

 

 

структурной схемой,

представленной

 

 

 

 

 

на рис.

11.13.

 

 

 

 

 

 

 

 

 

При непосредственном соединении

 

 

 

 

 

в параллель

 

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

 

 

 

 

 

генераторов (вместо действия блока

 

 

 

 

 

БПР с уравнительными соединениями)

 

 

 

 

 

в структурной схеме вместо ив 1 и ыв 2

 

 

 

 

 

появляется

общее

напряжение

иъ,

Рис. 11.14. Структурная схема

а связь

через

БПР

исключается.

при

образования напряжения

 

Уравнения

для

напряжения

 

 

 

 

 

параллельной

работе

генераторов.

При параллельной работе уравнение для напряжения генераторов имеет вид [15]:

2 £ qi£ q2 COS (6 12— р ),

(11.85)

ziz2

 

где для явнополюсных генераторов при неучете активного сопротив­ ления статора г.

Xqli

 

~

»

z

= ____L _ _

полк

V ¥ T ^

гДе 2полн — модуль полного эквивалентного сопротивления нагрузки

С4

и генераторов; g, b— эквивалентные активная и

реактивная прово­

димости нагрузки и генераторов соответственно;

 

 

g = g i+ g 2 + g H;

g i _ c°spi

; gz-

cos р2

9 >

z2

 

zi

 

 

 

&==:&!+& 2 + ьн;

=

V

sin P2

 

ЛГН

 

 

 

P= Pi—P2 = a rc tg ^ arctg — = 0, (так как г! = га = 0).

 

 

 

Г\

r2

 

 

 

 

Линеаризуя выражение (11.85), получим

 

 

 

Аи -

ди

АЕ Q1-

ди

АЕ Q2 '

ди

Аб12.

( 11. 86)

дЕ

Qi

дЕ 2

дЬ,

 

 

 

Q

 

 

 

 

Частные производные выражения (11.86) могут быть определены из (11.85):

ди

 

 

 

 

dEQl

и 22

\

А

zlz2

ди .

и(Гполн

 

1

/

Eq2 I EQl .

dEQ2.

^О^полн

\

z2

 

 

z2

 

J 12

>

J 12

(11.87)

 

ди

1

 

2EqiEQ2

 

P6i2

и 72

 

Z1Z2

 

 

и(гполн

 

где

и0 — значение напряжения

в соответствии с выражением (11.85)

при

исходных фиксированных значениях переменных Eq\\ E q2 ‘, б12.

Схема образования напряжения дана на рис. 11.14.

Таким образом, система дифференциальных уравнений, описы­ вающая параллельную работу двух разнотипных генераторных аг­ регатов, может быть представлена в следующем виде:

Уравнения моментов

 

7> Ф 1 =

 

фх +

^ p i

(Р ) T i------тг^~ X

 

 

dfc

 

3(Ч>1+ 1)

 

0[с

 

х

Ф1 '

 

dfc

ф2 -----~ ASJ2— D n (Р) рА612-

 

д (фх-г 1)

д (ф2+ 1)

 

 

 

d6 i2

 

 

 

 

 

 

АЕ QX-

дМх АЕ'Q2>г

 

 

 

 

 

дЕ<3х

 

дЕ<32

 

 

Т’дгРфг —

М ^ ф 2 + ^р2 (Р)Ф2— ^ - Х

 

 

dfc

 

 

d (ф2

 

дМ,

dfc

 

X

фх-

 

dfc

 

Д^Х2+ Т)22 (р) рА812—

(II.8 8 а)

 

д (Фх+1)

3(ф2+ 1)

 

6 ,,

 

 

 

 

дМг АЕ Q 1 -

■АЕ Q2-

 

 

 

 

дЕQI

 

дЕ<32

 

 

65

Уравнения электромагнитных процессов в роторах ГА

(pTdoi~\~ 1) AtBi = AuBi-{-pTd01\ I <ц (xdi Xdi),

{pTd02~\~ 1) Al'gj = AuB2-i' pTd02^Id2 (-^d2

Xdi) t

А/dl'-

ЗЛп . д E Q1

■AE Q2 ' i

^di

лл12i

 

dE Q1

dEq2

a6 i2

 

А/ d2 '

91d2 _AE

a / d 2 ■A£ <3 1 -

a/<ja

• A6'12-

Q2

аб,

 

dE Q2

Qi

 

 

Уравнения регуляторов напряжения

А« в1 = - ^ Аи + -4 ^ А / д -

 

 

 

«1

«1

 

 

 

— — [Ам -Г^б. п. рА (Qi Q2 )] ^к1 (Р)>

 

h

 

 

 

Амв2 = - ^ Ам+ - 7 ^ Л /^ -

 

 

 

«2

" 2

 

 

— [Аи— ^б. п. рА (QiQ2)l ^ к2 (Р)-

«2

 

 

 

 

 

Уравнения для напряжения

 

Эи

Д£ Q1"

ди

А£ Q2 "

д«

■А6 12-

Аи = -<Э£

2

36,.

Q1

 

q

 

 

 

Уравнение связи

6и = у -(ф 1—Фа)-

(11.886)

(П.88в)

(И .8 8 г)

(1 1 .8 8 д)

§ 8. Выбор метода исследования статической устойчивости (устойчивости «в малом»)

Несмотря на сделанные ранее допущения, что обобщенная функ­ циональная схема, предназначенная для изучения переходных про­ цессов в САЭС, содержит только два ГА с соответствующими систе­ мами регулирования, задача исследования такой электростанции ос­ тается достаточно сложной. Это обусловливается высоким порядком дифференциальных уравнений, описывающих переходные процессы

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

вквадратурах.

Анализ устойчивости такой системы даже в наиболее простых слу­ чаях при значительных допущениях труден.

В настоящее время еще не разработаны аналитические, инженер­ ные методы расчета устойчивости САЭС с параллельно работающими ГА. Исследование таких систем производится, как правило, на ЭВМ

6 6

1 2 1 , 2 2 , 23], хотя ее использование для решения подобных задач имеет ряд отрицательных моментов:

сложность определения вида и места включения корректирую­ щих связей. Такие задачи на ЭВМ могут быть решены только методом проб, что не может дать оптимального решения;

невозможность проведения быстрого прикидочного решения

впроизводственных (заводских) условиях;

невозможность контролировать правильность решения и пр.

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

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

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

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

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

Впервые исследование работы синхронной машины с помощью структурных схем было проведено М. М. Ботвинником [11 ] и В. М. Ма­ тюхиным [33]. Структурные схемы использовались ими для нагляд­ ного представления взаимных связей между величинами, входящими

67

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

Целесообразные варианты структурной схемы одиночно работаю­ щего генератора предложил О. М. Костюк [28]. Анализ этих схем дал возможность автору рассмотреть общие вопросы построения ра­ циональных систем регулирования возбуждения. В то же время опи­ санная в работе структурная схема генератора, работающего в па­ раллель с мощной сетью, не наглядна, проведение анализа по ней весьма сложно.

Более удобная структурная схема параллельной работы синхрон­

ного генератора

с сетью бесконечной мощности была предложена

И. Д. Урусовым

[52]. Использование в качестве выходной коорди­

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

Применение полученного критерия позволяет непосредственно по структуре установить, каким параметром машины или системы ре­ гулирования обусловлена неустойчивость.

Структурные схемы параллельно работающих генераторов (при­ менительно к береговым электростанциям) наиболее полно разрабо­ таны Г. В. Михневичем [38]. В этой работе на основании исследова­ ния структур дан всесторонний анализ динамических свойств системы автоматического регулирования возбуждения параллельно работаю­ щих синхронных машин и разработана методика синтеза ее рацио­ нальной структуры.

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

Наиболее полно отмеченные особенности проявляются в струк­ турных схемах, предназначенных для исследования динамики парал­ лельной работы ГА на ЭВМ [21, 22, 23]. При построении подобных схем использованы известные уравнения Парка-Горева в осях dq.

В качестве независимых переменных в уравнениях взяты напря­ жения, ток и скорость. Выбор этих переменных обусловил появление

ряда нелинейностей (типа х2; У х; ху\ х!у;

У х 2 + г/2; и др.), что при­

вело при моделировании к необходимости

использовать значительное

количество блоков нелинейностей. Кроме того, применение уравне­ ний в осях потребовало преобразования осей одного генератора к осям

68

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

Если в качестве независимых переменных взять не ток, а поток, структурная схема получается еще более сложной.

Проведение структурного анализа по этим схемам практически невозможно.

Как показано в литературе [18] для решения большинства задач, связанных с анализом динамики параллельной работы и синтезом оптимального регулирования, в качестве независимых переменных целесообразно выбрать электрическую координату — э. д. с. по про­ дольной оси Ed и механическую координату — угол б12 между рото­ рами генераторных агрегатов, а все остальные переменные предста­ вить как функции независимых переменных. Действительно, выбор в качестве независимой переменной напряжения неудобен для оценки устойчивости параллельной работы, так как при колебаниях синхрон­ ных генераторов напряжение на шинах (так же как и мощность на­ грузки) практически остается постоянным.

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

Структурная схема параллельной работы генератора с сетью, в ко­ торой за независимую переменную принят угол б12, была использо­ вана в работе [48] для определения причин возникновения и поисков путей устранения обменных колебаний мощности при параллельной работе ГА.

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

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

При исследовании непосредственно параллельной работы ГА в су­ довых условиях более рационально в качестве независимой перемен­ ной брать не угол между роторами, а скорости самих ГА (отклонение

скоростей) <рх и ф2. В данном случае угол б12 = (фх— ф2). Это по­

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

Структурная схема параллельной работы двух однотипных ГА, составленная на основании уравнений, где за независимые переменные приняты Edl\ Ейг и фх; ф2 была предложена в работах [7, 10]. Однако по структурной схеме, предложенной в этих работах, нельзя было

Соседние файлы в папке книги из ГПНТБ