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

Учебное пособие 1220

.pdf
Скачиваний:
8
Добавлен:
30.04.2022
Размер:
884.01 Кб
Скачать

2. Сравнить режимы активной и реактивной мощности генераторов при изменении угла.

Параметры генераторов взять из задачи 1.1, а значения ЭДС Eq и напряжения Uг из задачи 2.1. Принять значения Eq и Uг независимыми от режима работы генераторов.

Рекомендации. Внутренняя реактивная мощность генератора находится из

уравнения

 

=

г

 

 

 

 

 

 

 

=

2

 

δ

 

cos г .

(2.26)

 

 

 

 

 

 

 

 

 

Реактивная мощность, выдаваемая генератором:

 

 

 

 

Г

 

г

 

Г -

г2 .

(2.27)

Изменяя угол

г

от 0° до 360°, в формулах (2.26) и (2.27) можно построить уг-

ловые

характеристики реактивной мощности синхронных машин. Сопоставляя их с

 

 

 

 

 

 

 

 

 

 

 

угловой характеристикой , можно сравнить режимыРг и Qггенератора [1].

 

Типовая задача 2.5

 

 

Характеристическое уравнение простейшей электрической системы при

0 3 + 1 2 + 2 + 3

= 0.

(2.31)

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

 

 

Значения коэффициентов данного уравнения приведены в табл. 2.2. Требуется проверить устойчивость рассматриваемой системы с помощью критериев Михайлова, Гурвица и Рауса.

Таблица 2.2 Значения коэффициентов характеристического уравнения

№ варианта

0

1

2

3

№ варианта

0

1

2

3

1

0,2

1,55

4,0

0,52

18

0,3

1,52

5,3

0,80

2

0,3

1,45

3,8

0,53

19

0,4

1,42

5,1

0,82

3

0,4

1,35

3,6

0,55

20

0,5

1,32

4,7

0,84

4

0,5

1,25

3,4

0,58

21

0,6

1,22

4,3

0,86

5

0,6

1,15

3,3

0,59

22

0,7

1,12

4,1

0,88

6

0,7

1,05

3,2

0,6

23

0,8

0,91

3,9

0,9

7

0,8

1,0

3,1

0,62

24

0,9

0,81

3,7

0,91

8

0,9

0,95

3,0

0,64

25

1,0

0,71

3,4

0,93

9

1,0

0,85

2,8

0,66

26

1,1

0,63

3,1

0,94

10

1,1

0,75

2,6

0,68

27

1,2

0,53

2,9

0,95

11

1,2

0,65

2,4

0,70

28

1,3

0,43

2,8

0,96

12

1,3

0,55

2,2

0,71

29

1,4

0,21

2,6

0,97

13

1,4

0,45

1,78

0,72

30

1,5

0,15

2,3

0,98

14

1,5

0,35

1,61

0,74

31

0,8

1,22

3,8

0,66

15

1,6

0,25

1,56

0,75

32

0,9

1,12

3,6

0,68

16

0,1

1,72

6,0

0,76

33

1

0,91

3,4

0,70

17

0,2

1,62

5,7

0,78

34

1,1

0,81

3,3

0,71

21

Рекомендации. Для оценки статической устойчивости системы необходимо в характеристическое уравнение подставить значение р = jω и получить

где

D()= 0(jω) + 1

(jω)−1

+ + = (ω) + (ω),

(2.32)

комплексный полином

 

 

 

 

(ω)==

− −2ω2 ++ −4ω4

+

 

 

(ω) −1

ω − −3ω3

−5ω5

+ .

 

Если задаться серией значений ω в пределах от 0 до +∞, то каждому из них будет соответствовать некоторое значение полинома (2.32), которое на комплексной плоскости определяет точку. Геометрическое место этих точек образует кривую, которая называется годографом характеристического многочлена или кривой Михайлова, по которой можно судить об устойчивости системы [1].

Критерий устойчивости Михайлова можно записать в следующей форму-

лировке: для того чтобы система была устойчивой, необходимо и достаточно, чтобы вектор кривой Михайлова D ( ) при изменении от 0 до + монотонно

поворачивался против часовой

стрелки на угол

n

 

, где n – степень характе-

 

 

 

ристического уравнения, или, что то же самое,

чтобы кривая Михайлова обхо-

 

/2

 

дила только против часовой стрелки последовательно n квадратов координат-ной плоскости [2],[4]. При этом модуль кривой0Михайлова при всех значениях

должен быть отличным от нуля, т.е. D ( ) .

Кривая Михайлова для устойчивых систем всегда имеет плавную спиралевидную форму, проходит последовательно (без пропусков) все n квадратов, причем конец ее уходит в бесконечность в том квадрате координатной плоскости, номер которого равен степени характеристического уравнения.

На рис.2.6 показаны типичные кривые Михайлова для устойчивых систем, описываемых характеристическими уравнениями, начиная от первого (n=1) до

пятого (n=5) порядка. Для удобства сравнения коэффициенты

во всех случа-

ях приняты одинаковыми.

 

 

Признаками неустойчивости системы является нарушение числа и после-

угол n

.

 

 

 

довательности пройденных кривой Михайлова квадрантов координатной плос-

 

/2

 

) оказывается меньше, чем

кости, вследствие чего угол поворота вектора D (

22

Рис. 2.6. Кривые Михайлова для устойчивых систем На рис. 2.7 показаны кривые Михайлова для неустойчивых систем. Кривая

1 начинается на отрицательной вещественной полуоси, в кривой 2 нарушена последовательность прохождения квадрантов, а кривая 3 находится вся в одном квадранте при степени уравнения n=5.

нице устойчивости. При = 0, что соответствует границе апериодической устойчивости, кривая Михайлова выходит из начала координат, свидетельствуя о

По виду кривой Михайлова можно судить о нахождении системы на гра-

наличии нулевого корня (рис.2.8, а). При нахождении на границе колебатель-

пары

= 0

=

(рис. 2.8, б).

ной устойчивости кривая Михайлова при угловой скорости незатухающих ко-

лебаний

 

,+1 ± 0

 

 

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

мнимых корней

Рис. 2.7. Кривые Михайлова для неустойчивых систем

23

Рис. 2.8. Кривые Михайлова на границе устойчивости: а- апериодической; б - колебательной

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

Для оценки устойчивости системы по критерию Рауса составляется специальная таблица Рауса (табл. 2.3) [2], [4].

 

 

 

 

 

 

 

Таблица Рауса

 

 

 

 

Таблица 2.3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Номер i-

 

 

 

Номер к-го столбца

 

 

 

 

Коэффициент

той

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ

 

строки

 

1

 

2

 

 

3

 

4

 

 

….

 

 

−2

 

 

 

 

 

 

 

 

 

 

 

1

 

11

0

21

2

 

31

4

41

6

 

 

 

-

 

2

 

 

 

 

 

-

 

3

 

12

1

22

3

 

32

5

42

7

 

 

λ

0

1

 

13

 

23а4

 

 

33

 

43

 

 

 

 

1

 

=

λ2

 

=λ

 

 

=λ 6

 

=λ 8

 

 

 

 

 

 

 

4

1 3

 

1 5

 

 

1 7

 

1 9

 

 

 

2

1

13

 

=14 3

23

=24 5

2

=34 7

43

=44 9

 

 

 

 

 

2с23

− λ2с23

− λ2с43

− λ2с53

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Элементами первой строки этой таблицы служат коэффициенты характе-

ристического уравнения (2.31) с четными индексами, начиная с

 

элемента-

ми второй строки – с нечетными индексами. Элементы каждой 0последующей, а

строки находятся по формуле

+1,−2 λ−2 +1,−1 ,

 

 

 

 

 

где k - номер столбца;

 

 

 

 

(2.33)

 

λ−2= 1, −2/ 1,−1

-

 

 

24

 

 

 

 

 

 

 

 

 

i- номер строки, i ≥ 3;

 

 

 

 

 

 

 

 

 

коэффициент.

Критерий устойчивости Рауса формулируется следующим образом: для то-

го, чтобы система была устойчива, необходимо и достаточно, чтобы все эле-

11 > 0 , 12 > 0, 13

> 0,…, 1,+1

> 0 .

(2.34)

менты первого столбца были положительны:

 

 

 

Алгебраический критерий Гурвица широко используется во многих исследованиях устойчивости и базируется на анализе коэффициентов характеристического уравнения и полученных на их основе соотношений в виде неравенств [2], [4].

Для оценки устойчивости из коэффициентов характеристического уравнения составляется квадратная матрица Гурвица n-го порядка:

Правило составления матрицы следующее: по главной диагонали сверху и направо вниз записываются последовательно1 коэффициенты уравнения (2.31) в порядке их нумерации, начиная с до . Столбцы вверх от главной диагонали дополняют коэффициентами характеристического уравнения с последовательно возрастающими индексами, а столбцы вниз – с убывающими индексами. Все недостающие коэффициенты, т.е. коэффициенты с индексами больше n или меньше нуля заменяют нулями.

Отчеркивая в матрице диагональные миноры, получаем определители Гурвица:

Как видим, последний определитель

 

включает в себя всю матрицу Гур-

вица целиком. Если его раскрыть по

элементам последнего столбца, содержа-

 

 

 

 

 

щего только коэффициент

называемый обычно свободным членом характе-

ристического уравнения, то можно,

записать

 

=

 

.

Номер определителя Гурвица

соответствует номеру коэффициента по диа-

 

 

 

 

−1

 

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

т.е. , , …, и , , …, .

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

25

 

−1

> 0

 

.

= 0

 

устойчивости в нуль обратиться, прежде всего, определитель

Если при

этом

 

, то граница устойчивости определяется условием

 

 

. Это –

граница апериодической устойчивости, так как один действительный корень

находится на мнимой оси плоскости корней. Если

, то в нуль обращается

через

 

1,2

 

1

 

 

 

 

 

 

 

границе пары чисто мнимых

определитель

 

 

 

 

 

 

 

 

 

 

 

 

 

, что соответствует наличию на

> 0

 

 

корней

 

= ±−1. Это – граница колебательной устойчивости. При переходе

 

эту границу начинается самораскачивание системы с частотой

. Если и

дальше продолжать изменять коэффициенты,

то могут стать

отрицательными и

 

1

другие определители Гурвица, а

 

 

вновь может стать положительным. По-

вости системы: должны

и

−1

значит и

, еще не гарантирует устойчи-

этому положительность

 

, а−1

 

 

 

 

 

 

быть положительными также и остальные определите-

ли Гурвица.

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

В качестве примера рассмотрим систему, описанную линейным дифференциальным уравнением третьего порядка с постоянными коэффициентами (2.31).

Определитель Гурвица для данного уравнения имеет вид

1

=

Согласно;

критерию;

;

 

=

1 20 3

> 0

;

.

0 > 0;

 

1 > 0 2

> 0 3 > 0 2

 

 

∆3= 3∆2> 0

 

 

 

 

 

Гурвица эта система будет устойчивой, если

 

 

 

Следовательно, кроме положительности всех коэффициентов характери-

раскачиванию 1 20 3

> 0

 

 

 

 

 

 

 

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

соотношения

системы.

 

. Нарушение этого соотношения приведет к само-

 

 

 

 

 

 

 

 

 

 

 

 

 

Следует отметить, что в случае положительности всех коэффициентов ха-

рактеристического уравнения, необходимые и достаточные условия устойчивости сводятся к тому, чтобы среди определителей Гурвица были положительными все определители с четными или нечетными индексами. Такой критерий носит название критерия Льенара – Шипара; он требует раскрытия меньшего числа определителей, чем обычный критерий Гурвица, а поэтому особенно удобен при исследовании систем по уравнениям высокого порядка.

Решить задачу 2.5

26

Характеристическое уравнение простейшей электрической системы при отсутствии регулирования возбуждения0,12 3 + 0,1имеет2 вид

+1,72p+0,5=0.

Требуется проверить устойчивость системы с помощью:

1)критерия Михайлова;

2)критерия Гурвица;

3)критерия Рауса.

Типовая задача 2.6

Для электрической системы, принципиальная схема которой приведена на рис. 1.1, требуется исследовать зависимости коэффициента запаса статической устойчивости от коэффициента мощности нагрузки, длины линии электропередачи, числа проводов в фазе. Исходные данные принять в соответствии с заданным вариантом по табл. 1.1 и 1.2. Расчеты выполнить при условии, что на генераторах станции установлены АРВ пропорционального

действия. Принять Ро = Pн, Qo = Qн.

При исследовании влияния на коэффициент запаса числа проводов в фазе принять для линий номинальным напряжением 110 кВ марку провода АС

150/24, радиус провода

= 7,9 мм, среднегеометрическое расстояние между

фазами

 

 

 

линий номинальным напряжением 220 кВ - АС 300/39,

Dср = 5 м; для пр

 

 

 

 

 

 

 

= 12 мм, Dcp = 8 м; для линий номинальным напряжением 330 кВ - АС

400/51,

 

 

 

 

 

 

 

 

 

 

пр

 

= 13,75 мм, Dср = 11 м. Для линий всех напряжений шаг расщепления

в одной пр

 

 

 

= 400 мм.

 

 

 

 

 

 

 

фазе

 

 

 

 

 

 

 

 

 

Рекомендации.

 

 

 

 

 

 

 

Реактивная мощность нагрузки определяется по формуле

 

 

 

 

 

 

 

н

н

н

 

(2.35)

 

Эквивалентный радиус

расщепленной фазы

 

 

 

 

=

 

 

 

 

 

 

 

 

прэ =

пр ср

 

 

 

 

 

 

 

 

 

 

−1 ,

 

(2.36)

 

где N - число проводов в расщепленной фазе;

 

 

 

 

 

радиус провода;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

пр --среднегеометрическое расстояние между проводами одной фазы.

 

 

 

 

 

 

0 = 0,1441 прэср + 0,0157

,

(2.37)

 

Удельное индуктивное сопротивление линии электропередачи

 

где Dcp - среднегеометрическое расстояние между фазами.

27

Типовая задача 2.7

Для электрической системы, принципиальная схема которой приведена на рис. 1.1, требуется исследовать влияние шунтирующего реактора, включаемого в начале линии, на статическую устойчивость системы. Исходные данные принять в соответствии с заданным вариантомхр по табл. 1.1 и 1.2. Сопротивление шунтирующего реактора принять = 500 Ом. Рассмотреть следующие случаи:

без АРВ на генераторах станции и отсутствии шунтирующего реактора;

при неизменной ЭДС генераторов станции без АРВ и включении шунтирующего реактора;

при увеличении ЭДС генераторов станции для поддержания напряжения на шинах генераторов и наличии шунтирующего реактора.

Построить угловые характеристики мощности для каждого случая.

Рекомендации.

 

+ + л + АТ +

р

 

 

 

Сопротивление системы при включении шунтирующего реактора в начале линии

 

=

 

( + )( л+ АТ),

(2.38)

 

 

 

 

где xр - сопротивление шунтирующего реактора.

ЭДС генераторов станции без АРВ при отсутствии шунтирующего реактора в начале линии определяется по формуле (2.2). Напряжение на шинах генераторов при отсутствии шунтирующего реактора находится по формуле (2.13).

Уравнение мощности при неизменной ЭДС и включении шунтирующего

 

 

 

 

реактора

 

 

= sinδ.

 

 

г

 

 

 

(2.39)

 

+ ( л+ АТ)

 

 

Уравнение мощности при увеличении ЭДС

sinδ.

(2.40)

 

=

 

 

Типовая задача 2.8

Для электрической системы, схема замещения которой представлена на рис. 2.4, требуется определить предел мощности и предел устойчивости. Необходимые исходные данные взять из решения задачи 2.3. Принять номинальную мощность системы 7800 МВт, постоянную инерции 10 с.

Рекомендации

11sin 11

+ 12.

(2.41)

P= 2

Предельная мощность станции определяется по формуле

 

28

Предельный угол по условию устойчивости

 

δ12 = (2 12 1−1+

с

),

(2.42)

где Tj - постоянная инерции генераторов станции

; с

 

 

Tjc - постоянная инерции генераторов системы.

 

Типовая задача 2.9

Для электрической системы, принципиальная схема которой представлена на рис. 1.1 с установленными на электростанции явнополюсными генераторами без АРВ, с АРВ ПД и АРВ СД, требуется построить векторную диаграмму и угловые характеристики мощности. Параметры элементов системы взять из решения задачи 1.1. Принять Xq = Xd/1,5.

Рекомендации.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ЭДС явнополюсного генератора и напряжение на шинах станции находят-

ся по формулам:

 

 

= ( + 0 ∑)2

+ ( 0 ∑)2;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= ( +

0 ∑)2 + ( 0 ∑)2 ;

(2.43)

Углы

 

 

 

г= ( +

0

)2

+ ( 0 )2,

 

= Xq + Xc .

 

Eq, , Uг, U

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

0 ∑

 

 

 

 

 

 

 

между векторами

 

 

 

 

 

 

 

 

 

 

 

 

;

 

 

 

 

 

 

 

 

δ =

+ 0

 

 

(2.44)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

δ

 

= + 0

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

δ

 

=

 

0

 

.

 

 

 

ЭДС холостого хода

 

 

+ 0

 

 

 

 

 

 

 

 

 

 

=

−′

,

 

(2.45)

 

 

=

 

 

− ′

 

 

 

 

где

) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos( -

 

 

 

 

 

 

2

∑ ∑

 

 

 

 

 

 

 

=

 

 

Характеристики мощности

 

sin

+ 2

 

 

 

 

 

 

 

 

 

 

 

 

 

sin2 .

(2.46)

 

 

 

 

 

 

 

 

 

 

29

 

 

 

 

 

 

 

 

 

или ориентировочно

 

=

 

 

 

 

 

 

 

 

sinδ.

 

 

(2.47)

=

 

 

2

∑ ′ ∑

 

г

 

sinδ

- 2

sin2δ.

(2.48)

 

 

 

2

 

∑ с

 

 

 

 

=

г

sinδ -

 

2*

с

sin2δ.

(2.49)

 

 

 

 

Типовая задача 2.10

Для системы, схема которой представлена на рис. 1.1, требуется исследовать статическую устойчивость без учета и с учетом демпфирования, найти частоту и период(собственных колебаний, а также построить зависимость изменения угла ) от времени при отклонении ротора на 1 градус от положения установившегося режима при 0°, 60°, 90°, 100°.

Параметры элементов схемы замещения принять из решения задачи 1.1. Демпфирующую мощность для вариантов 1–5 принять 50, для 6–10 – 60, для

11–15 – 70, для 16–20 – 80, для 21 – 100.

Под демпфированием понимают процесс подавления электромеханических колебаний в системе, либо уменьшение их амплитуды до допустимых пределов с помощью устройств или приспособления, поглощающих энергию колебанийдемпферов [2]. Демпфирование может быть обусловлено как естественными факторами – потерями на трение, так и действием регуляторов, и появлением дополнительного асинхронного момента и соответственно мощности из-за протекания токов в демпферной обмотке и демпферных контурах, вызванных скольжением ротора синхронной машины относительно поля статора.

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

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

ностях может быть записано в виде

М Д .

 

22 = 0

(2.50)

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

ности синхронной машины по скольжению= , т/.е.

Д ас .

30