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

Переходные процессы / Лабораторные работы / Лабораторная работа 14

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

Министерство образования Российской Федерации

КУБАНСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНОЛОГИЧЕСКИЙ УНИВЕРСИТЕТ

Кафедра электроснабжения промышленных предприятий

ОПИСАНИЕ лабораторной работы № 14.8

по курсу «Переходные процессы в электрических системах»

ИССЛЕДОВАНИЕ УСТОЙЧИВОСТИ ПРОСТЕЙШЕЙ ЭЛЕКТРИЧЕСКОЙ СИСТЕМЫ

Краснодар 2003

УДК 621.311.1.001.45(076.5)

Печатается то решению Редакционно-издательского совета университета

ИССЛЕДОВАНИЕ УСТОЙЧИВОСТИ ПРОСТЕЙШЕЙ ЭЛЕКТРИЧЕСКОЙ СИСТЕМЫ

Лабораторная работа № 14.8 по курсу «Переходные процессы в электрических

системах»

Краснодар: изд. КубГТУ, 2003. – 18 с.

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

Составители:

Б.А. Коробейников

 

Л.И. Деревцова

Pd.эк

3

1. ОБЩИЕ СВЕДЕНИЯ

1.1. Д и ф ф е р е н ц и а л ь н о е у р а в н е н и е п р о с т е й ш е й с и с т е м ы

Одной из основных задач, решаемых на аналоговой вычислительной машине (АВМ), является интегрирование обыкновенных дифференциальных уравнений.

Динамика простейшей электрической системы описывается неоднородным дифференциальным уравнением 2-го порядка с постоянными коэффициентами [1]

T

 

d2δ

+P

 

dδ

P = ω P

,

(1.1)

dt

dt

j

 

d.эк

 

0

0 тδ

 

 

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

δ- угол поворота ротора относительно вращающегося магнитного потока статора, рад;

- демпферный коэффициент, относительные единицы (о/е); ω0 - угловая промышленная частота;

P - электрическая внутренняя мощность генератора, о/е; Pтδ - механическая мощность первичного двигателя, о/е;

Уравнение (1.1) справедливо при ddtδ << ω0 , т.е. при скоростях ротора ге-

нератора, близких к синхронной.

Решение уравнения заключается в определении зависимости δ(t) и

скорости ddtδ при заданных параметрах Tj, Pd.эк, P и заданном воздейст-

вии Pтδ.

Электрическая мощность генератора определяется соотношением

P = P +

E U

sin δ,

(1.2)

 

11

X

 

 

 

 

 

где P11 - мощность, характеризующая составляющую вращающего момен-

та генератора, обусловленную потерями в электропередаче; E - ЭДС генератора;

U - напряжение системы;

X - взаимное сопротивление между генератором и системой.

4

Слагающая Pd.эк ddtδ в уравнении (1.1) учитывает влияние демпфи-

рования (т. е. появления мощностей, зависящих от скорости движения ротора).

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

1.2. Д и н а м и ч е с к а я у с т о й ч и в о с т ь э л е к т р и ч е с к о й с и с т е м ы б е з у ч е т а д е м п ф и р о в а н и я

Система описывается уравнением

пр о с т е й ш е й

вл и я н и я

T

 

d2δ

=(P

P) ω .

(1.3)

dt2

j

 

тδ

0

 

На рис. 1.1 изображены две системы, причем система I соответствует нормальному, а система II – послеаварийному режиму.

На рис. 1.2а показан ход рабочей точки при переходе системы из положения I в положение II. Вследствие изменения параметров системы и соответствующего изменения внутренней электрической мощности генератора происходит переход рабочей точки из положения "a" в "b". Воз-

никает небаланс P0 электрической и механической мощности, который

приводит к появлению ускорения α =

P и относительному перемеще-

 

Tj

нию ротора ω. При значении угла δ = δ0II (точка «c») электрическая и

механическая мощности уравновешиваются, однако наличие кинетической энергии, запасенной ротором в процессе ускорения и пропорциональной площадке «abca», приводит к дальнейшему увеличению угла δ. При движении от точки «c» к точке «d» ротор испытывает торможение под действием электрической мощности, которая теперь больше ускоряющей механической мощности. При движении от «c» к «d» кинетическая энергия ускорения переходит в потенциальную энергию. В точке «d» вся кинетическая энергия израсходовалась, скорость движения ротора равна синхронной, ротор начинает тормозиться за счет накопленной потенциальной энергии и действия на ротор избыточного отрицательного момента. Под действием избыточного момента ротор снова подходит к точке «с», обладая кинетической энергией, пропорциональной площадке «cdec». В точке «b» относительная скорость движения ротора и кинетическая энергия ротора равны нулю и процесс начинается сначала.

 

 

 

 

 

 

5

I

x

 

 

II

x

 

 

 

U = const

 

U = const

 

 

 

 

 

E

x

 

 

E

 

 

 

 

 

 

 

 

 

а

 

 

 

б

 

 

а) PI =

E U

sin δ

б) PII =

E U

sin δ

 

 

0,5 x

 

 

x

 

Рис. 1.1 – Схемы эквивалентные системы

 

P

 

 

 

 

 

 

 

 

PII

PI

 

 

 

d

 

 

 

 

 

 

P

a

c

 

d'

P

2

e

т

P

 

 

 

 

 

 

 

 

 

0

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

δ

0

δI

δII

δ

max

0

 

0

0

 

 

 

dδ

 

 

a)

 

 

 

 

 

 

 

dt

 

 

 

 

0

 

 

 

δ

 

 

 

 

б)

0

 

 

 

 

 

0

 

 

 

δ

0

 

 

 

 

в)

 

 

t

 

 

 

 

 

 

 

Рис. 1.2

 

P

 

 

 

 

 

 

 

 

 

PI

 

 

e

 

 

Pт

a

с

d

 

b

 

 

 

 

 

 

 

 

PII

 

 

 

 

 

 

 

δ

δI

δII

δ

max

= δ

кр

0

0

 

 

dδ

a)

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

δ

 

 

 

 

б)

 

 

 

 

 

 

δ

 

 

 

 

 

Выпадение из

 

в)

 

синхронизма

 

 

 

 

 

t

 

 

 

 

 

Рис. 1.3

 

 

 

 

6

В зависимости от соотношения площадок ускорения и торможения различают три возможных режима системы:

а) если при колебаниях ротор не достигает критического угла, соответствующего точке «d» на рис. 1.2а), переход от режима I к режиму II устойчив. Вся энергия, полученная при ускорении ротора («авса») уравновешивается энергией торможения («сdec») до перехода к точке «d». Площадка возможного торможения «cdd’ec» больше площадки ускорения;

б) в критическом режиме (рис. 1.3) система находится на границе устойчивости (точки «d» и «d’» совпадают). В этом случае равновероятно как сохранение устойчивости, так и ее нарушение при переходе oт режима I к режиму II (см. рис. 1.1);

в) неустойчивый режим (рис. 1.4) наблюдается при условии, когда площадка ускорения больше площадки возможного торможения. При этом в точке «d» не вся кинетическая энергия ротора израсходована и угол продолжает возрастать. Далее рабочая точка попадает в область ускоряющего момента, что приводит к прогрессирующему увеличению угла δ и выпадению генератора из синхронизма.

1.3. Д и н а м и ч е с к а я у с т о й ч и в о с т ь э л е к т р и ч е с к о й с и с т е м ы п р и у ч е т е д е м п ф и р о в а н и я

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

ной ω = ωс −ω) и появлением дополнительной асинхронной мощности Рас, вызванной скольжением ротора. Демпфирование при малых изменениях скорости приближенно учитывается выражением

Рас = Рdэк s ;

(1.4)

Рdэк = KR Pd ,

(1.5)

где Рас - составляющая электрической мощности генератора, учитываю-

щая демпфирование; s – скольжение;

KR , Pd - составляющие демпферного коэффициента.

7

P

 

 

 

P

 

 

 

f

 

 

P

 

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Pтδ

 

 

 

 

I

a1

b

d c

 

 

 

 

 

 

 

 

e

Pтδ

 

 

k

e1

f1

a

c

 

 

 

 

d

 

 

 

 

 

b

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

δ

0

 

 

 

δ

δ0I δ0II

 

δd

 

δ0II

δmax1 δmax2

 

 

δ0I

 

dδ

 

 

 

dδ

 

 

 

 

dt

 

 

 

dt

 

 

 

 

0

δ

0

δ

 

 

0

δ

 

 

 

 

 

t

Рис. 1.4

Рис. 1.5

8

Обычно Рdэк имеет отрицательное значение ( KR < Pd ), что соответствует случаю положительного демпфирования. Положительное Рdэк

( KR > Pd ), соответствующее отрицательному демпфированию, имеет ме-

сто при неправильной настройке регуляторов, либо ори больших сопротивлениях в цепи статора. Электрическая мощность с учетом демпфирования

P = P

+

E U

sin δ+P

 

s .

(1.6)

 

 

11

 

X

dэк

 

 

 

 

 

 

 

 

Наличие положительного демпфирования

( Pdэк) изменяет

харак-

теристику P(δ), которая при отрицательном скольжении пройдет выше характеристики, не учитывающей демпфирование. При этом площадка ускорения уменьшается от значения «авса» до «авfа» (рис. 1.5). Соответственно уменьшается и площадка торможения. Это приводит к уменьшению размаха колебаний, а, следовательно, к увеличению запаса устойчивости.

При отрицательном демпфировании ( Pdэк > 0 ) кривая элек-

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

1.4. П р и в е д е н и е у р а в н е н и я к в и д у , м о д е л и р о в а н и я

ди ф ф е р е н ц и а л ь н о г о

уд о б н о м у д л я

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

Преобразование уравнения (1.1)

d2δ

 

P

P

 

Pd

эк

 

dδ

 

 

 

= ω

Т

 

 

 

.

(1.7)

 

 

 

 

dt2

0

 

TJ

 

TJ

dt

 

 

 

 

 

 

 

Система уравнений, эквивалентная уравнению 2-го порядка (1.1):

ds

 

P

 

P

 

Pd

эк

 

 

 

=

Т

 

 

s ;

(1.8)

dt

 

 

 

 

 

 

T

 

T

 

 

 

 

 

J

 

 

J

 

9

 

P

 

 

 

 

f1

e

 

 

 

 

 

 

 

 

 

d1

a1

 

c

 

 

 

Pтδ

 

 

 

 

c1

b

f

e1

 

 

a

 

 

 

 

 

 

 

d

 

 

 

 

 

 

 

δ

0

δ

 

 

δ

 

δ

δ

0

 

max2

max1

 

max3

 

 

 

 

 

dδ

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

δ

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 1.6

 

 

 

10

ds

= −ω s .

(1.9)

 

dt

0

 

 

 

Уравнением (1.7) удобно пользоваться для составления структурной схемы, а системой уравнений (1.6 – 1.9) – для выбора коэффициентов передачи операционных блоков.

1.5. С о с т а в л е н и е с т р у к т у р н о й с х е м ы

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

Момент начала интегрирования определяется моментом подачи соответствующих напряжений на входы интегрирующих блоков, что делается с помощью контактов реле управления. При размыкании контактов интегрирование прекращается. Повторным замыканием можно продолжить интегрирование далее с той же точки, ибо емкости обратной связи «запомнили» на интегрированное напряжение. Для возвращения схемы в первоначальное состояние необходимо разрядить емкости, закоротив их.

1.6. В ы б о р м а с ш т а б о в з а в и с и м ы х п е р е м е н н ы х и в р е м е н и

Масштабы переменных величин выбирается, исходя из максимальных значений переменных и шкалы напряжений в АВМ, отобразивших эти переменные (обычно ± 100 В). В соответствии с этим используются формулы:

m y =

100

;

(1.10)

 

 

ymax