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

книги / Переходные процессы в электрических цепях, содержащих машины переменного тока

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

По первому закону

Кирхгофа

имеем для узла

схемы

рис. 4-1:

 

 

 

 

 

 

1а +

1ы =

1*а-

(4-П)

 

+

^

=

V

(4-12)

Из векторной диаграммы (рис. 4-3) имеем:

 

6

= coe / +

ô - - J - ,

(4-13)

где 6— текущее значение угла

между продольной осью

d машины и магнитной осью фазы а ее статора;

со0 — электрическая

синхронная

угловая

частота

(иначе — синхронная скорость);

в

относитель­

ных единицах

в>0 =1;

 

 

 

 

 

ю0/ — текущее

значение угла

между

вектором

напря­

жения на шинах бесконечной мощности и осью

фазы а

статора генератора;

 

 

 

 

ô — угол между вектором

э. д. с. Е

генератора, совпа­

дающим с поперечной осью машины

и вектором

напряжения U на шинах бесконечной мощности.

Дифференцируя равенство

(4-13), получим:

 

 

É1

(On

rfS

 

d *

 

 

(4-14)

dt

=

1 + dt

 

 

 

dt

 

 

 

 

Поскольку угол à, являющийся искомой неизвестной

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

d *

О

и

d 8

 

(4-15)

dt

dt =

© — <v

Этим и объясняется, почему при записи уравнений

(4-1),

(4-2), (4-7) — (4-10)

мы

считали

угловую скорость

ю

по­

стоянной.

 

 

 

 

 

 

Если по окончании этого расчета окажется, что угол

ô

действительно очень мало меняется на протяжении 0,1 сек после возникновения аварии, то можно будет ограничиться только первым приближением и не решать снова систему исходных уравнений. Если же окажется, что угол ô будет

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

Мы получили

10

операторных

уравнений

(4-1)

- (4-4),

(4-7) — (4-12) с

10 неишесгнами и ы, U^,

4 d,

l d,

/„, 1Ы,

/ сq, / Hd и 1аг

Олень

легко

исключить

из этих уравнений

неизвестные I ad,

l u ,

Ula, Ulu/ и

Y ,.

 

 

 

 

Решим

оставшиеся пять уравнений относительно I d и / в:

 

 

 

 

а

д

д

 

 

 

 

Е

 

 

0

 

2 012/? +

1

0

 

 

Ud

 

1,25

 

 

Р

0,382/? +

0,0894

и я — 0.0С695— 1,25р

 

1

 

0,382

 

u d

— 0,617

 

 

0

0,999/7 +

1,084

 

0,617/? + 0,995

 

0

 

0,999

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

- 0,382

 

 

 

 

 

 

 

 

0.382/? +

0.С894

 

 

(4-16)

 

 

 

 

 

- 0,999

 

 

 

 

 

 

 

 

0,999/7 +

1.С84

 

 

 

 

 

 

 

1 116/»+

2,06

E 2 012/7 +

1

 

 

 

 

 

—0.CC695

и а

Р

 

 

I

= - * -

A

 

0

 

Vq

1

 

 

 

4

A

0,617p +

0,99ï

и а

0

 

 

 

 

 

 

0,617

V,

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

0

 

 

 

 

 

 

0,382p +

0,0894

 

 

0,382

 

 

(4-17)

 

 

 

0,382

0,382/» + 0.C894

 

 

 

0 ,9 9 9 /» + 1.C84

 

- 0 ,9 9 9

 

 

 

 

 

 

0,999

0 ,9 9 9 /» + 1.C84

 

 

Здесь

Д

представляет собой

следующий определитель

пятого

порядка:

 

 

 

 

 

 

 

 

 

 

1 116/7 +

2,06

 

 

О

 

 

2 012/»+1

Д -

 

— 0,00695

 

 

1,25

 

 

р

 

 

 

0

 

0,00095-1,25/»

1

 

 

 

0,617/» +

0,995

- 0 ,6 1 7

 

0

 

 

 

 

0,617

0,617/» +

0,995

О

 

 

 

 

 

О

 

 

 

О

 

 

 

 

 

 

0,382 р +

0,0894

 

— 0,382

 

 

(4-18)

 

 

 

0,382

0,382/» +

0,0894

 

 

 

0,999/»+

1,084

 

— 0,999

 

 

 

 

 

 

0,999

0,999/»+1,084

 

 

Раскрывая определители Д, Д/(г и

Д/д,

получим:

 

А = 2 358/»5 + 5 972/»4 +

8 942/»3 + 6 911/»2 +

5 691/» +

7,802;

 

 

 

 

 

 

 

 

 

 

 

(4-19)

=

Е (1,483/»* + 3.4С6/»3 +

5,01/»2 +

3,511/» + 3 ЗС4) —

— Uй (2 012/» + 1) (С,9154/»3 +

2,5876/»2 +

2,7036/» -

0,2716) —

— Uq (2 012/» + 1) (0,914/»2 + 2,9546/» + 2,6389);

(4-20)

д л =

Е (0.1864/?*-f- 0,1928/? +

0,2907) + t/ d(980,9/7*+3l61,4/7*+

+ 2954,6/7 + 4,016) — и ц (981,6/7* +

 

2876,8 /73 +

3191,9 /7я —

 

 

 

— 88,2/7 — 0,5348).

 

 

 

(4-21)

Введем обозначения:

 

 

 

 

 

 

 

 

 

 

 

Z ip) = 2 358/7* +

5 972/7* + 8 942/т* + б 911/7* +

 

 

 

 

+

5 691/7 +

7,802;

 

 

 

(4-22)

УliP) = (2 012/7 + 1) (0,9154/т* +

2,5875/7* +

2,7636/7 —

 

 

 

 

— 0,2716);

 

 

 

 

(4-23)

 

Уг(р) = (2 012/7 +

1) (0,914/?* +

 

2,9546р +

2,6389); (4-24)

 

Уzip) =

980,Эр3 + 3161,4р2 +

2954,6/7 +

4,016;

(4-25)

Уtip) = 981,6/7* +

2876,8/7* +

3191,9/7* — 88,2/7 — 0,5348;

 

 

 

 

 

 

 

 

 

 

 

 

 

(4-26)

У-iiP) = 1,483/7* + 3,406/7* +

5,01/7* +

3,51 Ip + 3,304;

 

 

 

 

 

 

 

 

 

 

 

 

 

(4-27)

 

Уsip) =

о, 1864/7* +

0 ,1928/7 +

0,2907.

(4-28)

Теперь токи

I d и I q могут быть написаны так:

 

 

I d =

E

y7(p)

U

, ™

 

 

*

Zip)

 

(4-29)

 

 

 

Zip)

 

 

 

m

 

 

 

 

 

I . =

E

Ytip) +

u

Ï & L

f,

Y4(p)

 

(4-30)

 

 

 

Z(p)

 

 

d

Zip)

 

 

щ Zip)

 

 

Подставив в уравнение

(4-3)

I d

из уравнения

(4-29),

G(p)

из уравнения

(4-5) и

 

xd (p)

из

уравнения (4-6), по­

лучим после преобразований потокосцепление с обмоткой статора по продольной оси:

(4-31)

где выражения Уъ(р), Yв(р)

и

Yt(p) определяются следую­

щим образом:

 

 

 

 

 

 

 

 

 

Yt(p) = Yх(р)

1

 

 

- =

0

116д + 2,06) (Г,9154/7» +

 

 

 

+

2,5875/7* +

2,7636р — 0,2716);

(4 32)

Y,{P) =

Yt{p)

112'c61P2~ ' - f -

=(1116/7 +

2,06) (0,914/7* +

 

 

 

 

 

+ 2,954/7 +

2,6389);

 

(4-33)

Y,(p) =

704/7b +

2 169/7* + 3 345/7* -f 2 984/7* +

1 997/7 +0,992.

 

 

 

 

 

 

 

 

 

 

 

 

(4-34)

При p = 0

в

соотношениях

(4-29) — (4-31)

мы должны

получить

 

iu•

и т|>Л при

стационарном (т. е. предшест-

вующем аварии,

нормальном) режиме:

 

 

.

_

j,

3,304

, ..

0,2716

..

2,6789

 

 

г40 —

0 ? 802

** 7,802

 

 

,0 7,802

 

 

 

 

 

 

 

 

 

 

 

 

 

-

0,2907

, „

4,016

 

,

0,5348

 

 

1а0 — Ел

7,802

-4- иаg

7,802

4-

7,802

 

 

*

 

*

 

 

 

 

.

 

с- 0,992

 

2,06-0,2716 .

2.06-2,6389

 

 

 

7,802

 

 

7,802

+ “<*

7.802

 

 

 

 

 

=

1,0247(1,0205).

 

 

В скобках приведены значения id0,

и i|)d0, полученные

из расчета нормального режима.

 

 

 

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

Поэтому

ud и uq получат

определенные изменения

Aud и

Auq. Так

как регуляторы

напряжения

учить ваться

не бу­

дут, то щ и Е будут величинами постоянными и

 

àuf = 0;

Д£ = 0;

LUf(p) = Q\

ЬЕ(р) = 0.

(4-35)

Изменениям

лапласовых изображений продольной и

поперечной составляющих напряжения на шинах бесконеч­ ной мощности U4 и Д Uf будут соответствовать опре-

деленн&е изменения лапласовых изображений продольной и поперечной составляющих токов Д I d и Д 1Ч и потокосцепления по продольной оси ДЧ^:

д / , = -

Уг(р) 1(J

Yzip) \ tr .

(4-36)

Zip)

 

Zip)

4

 

 

 

м ч =

д и а

Yiip) » it .

(+37)

Zip)

4'

 

Zip)

а

 

à w a =

Yd p )l u

+

Yeip) AU,.

(+38)

 

Zip)

 

Zip)

4

 

Приравняв затем нулю выражение Zip)

 

 

Zip) = 2 358р* +

5 972р* +

*- 942р* +

6 911/>2 +

 

+

5691/> +

7,802 = 0,

 

(4-39)

найдем корни этого характеристического уравнения:

р л = — 0ДЮ1375;

 

/?2 =

— 0,100656 + /1 ,0 2 5 ;

 

р г =

— 0,100656 —у 1 025;

(4-40)

P i =

— 1.1657

+ /0 ,9 4 7 3

;

/>5 =

— 1,1657

—/0,9473.

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

=

= 727 pûa в 2-315 "*■

(4-41)

Все остальные свободные поля обусловлены наличием обмотки статора. Так как в рассматриваемом случае цепь статора имеет две параллельные ветви (ветвь нагрузки и ветвь линия — шины бесконечной мощности) и на роторе

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

уравнения.

 

,

Токи и потокосцепления переходного режима

A I d,

и АТ^, определяемые

соотношениями (4-36) — (4-38), обу­

словлены внезапным

изменением напряжений

ud и ищпо

величине или, иначе говор», внезапным изменением напря­ жения и по величие от значения «0 до значения и0 + Аи и последующим непрерывным изменением его по фазе.

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

id =

iM +

Д iax

iq =

i„о +

Д V

(4-42)

Фл =

Ф.о +

Д Фл;

Ф, =

Ф,о + АФ<г

(4-43)

Изменения напряжений

Aud

и

A ы, представляют

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

нием

на шинах бесконечной мощности

(«0 +

Д ы) е‘ (90°“ 4)

и его

стационарным значением и0 е!(90°—ес) (рис.

4-44):

Aud =

 

СВ' — ОВ — СС sin ô — ОС sin ô0 =

(r/0-f A и) sin ô —

— «о sin ô0 = A «sin ô0 -f («0 +

A «) (sinô — sin ô0); (4-44)

Д uq =

 

CA' — OA =

C C cos ô —• OC cos ô0 =

(Uq -f- A u) cos ô —

—«o cos ô0 = A и cos ô0 -f (u0

A «) (cos ô — cos ô0). (4-45)

Переходя в (4-44) и (4-45)

к изображениям, получим:

Д а а =

----------- +

(«о + А и) L

sin ô — Sin ô0 = ---------- L +

 

 

P

 

 

 

 

P

+

(«o + A «) J

e ~ pt (sinô — sinô0) dt =

Д и sin S0

+ (uo + Д «) F

(4-46)

Д U . = —

-°*-*° +

(Ыо + Д и) L {cos6 — cosô0) = А“ С— - +

 

 

P

 

 

р

-f

(u0 +

Д « ) |

е~ pi(cos ô — cos ô0)cft=:

+

 

 

 

-f(w0 +

Д и ) ^ ^ ) .

(4-47)

С

учетом полученных

равенств перепишем

(4-36) —

- (4 - 3 8 ):

 

 

 

 

- i u c o s 8

, ^ ! - - ( «

. +

4 « ) - ^ f ! f,0> ):

(4-18)

д

Д U sin ôe ^

' p)-

+

(u0 + Au) Х М

Fi(p) —

 

 

 

pZip)

 

 

 

zip)

 

 

— Д и cos ôü Y~——

(u0 -f Д u)

Vj(P)

 

(4-49)

 

 

pZ(p)

 

 

 

Zip)

 

 

Д Y rf= A u s in ô e ^

+

 

(н» + А « )?“

Л (Р) +

 

 

 

pZ(p)

 

 

 

Zip)

 

 

-f

Д u cos ô0 ^ t~ L -f (Hg +

Д u)

Fa{p).

(4-50)

 

 

pZip)

 

 

 

Z(p)

 

 

Заметим, что в операторных выражениях

(4-48) — (4-50)

члены

 

 

 

 

 

 

 

 

 

— ^-Д ы sin ô 0; XML д и cos ô0;

(u0 -f Ди) будут рас-

pZ[p)

pZ{p)

 

 

 

Z{p)

 

 

 

врываться по единичной функции.

 

 

 

Вторые и четвертые слагаемые правых частей равенств

для Д I d,

Д 1 Ч и Д

раскроем с

помощью

теоремы

умно­

жения (иначе говоря, ишеграла

Дюамеля).

Полагая, на­

пример:

 

 

00

 

 

 

 

 

 

 

 

 

 

 

 

 

Д~~' {^(р)} =

L~x I J

 

 

(sinô — sinô0)rf^J =

 

 

 

 

о

 

 

 

 

 

 

 

— Л (0 =

si.i Ô — sin ô0;

 

(4-51)

 

 

 

 

оо

 

 

 

 

 

 

 

 

 

L~l {^з(P)} — £- '

| J 6~p(cos ô — cos ô0) dt j =

 

 

 

 

— U t)

— cos ô — cos ô0;

 

 

 

(4-52)

 

 

' Y'(P) )

_

 

Y dPk)'k* .

 

(4-53)

 

 

.Z(P) f

f î( )

L

 

Z(Pk)

 

 

 

 

 

 

k= 1

 

 

 

 

 

 

 

 

' Y*(P) 1 _ $ /a

_ VI

Уг(Рк) A* .

 

(4-54)

 

 

. Zip) )

fi{

Ъ

 

Z’ipù

 

получим:

 

 

 

k= 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

Д U)

Л 0»)} =

(«o +

 

*

 

 

 

 

 

L 1{(«• +

A U) j 7 , ( т )Л (/~ T) <fr

 

 

* 5

 

,

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

(«o -f Au) f Y] Pk

e *

rsin ô (/ — t) — sin ô0] dx^z

 

 

 

J J * .

PkZ (Pk)

 

 

 

 

 

 

 

 

 

 

0 k—1

 

 

 

 

 

 

 

 

 

 

K

+ Д U) [sin 8 ( / - 0

- s t a

8, ] j j -

[ j

^

f

]

dx

:

 

 

 

 

 

 

 

k= 1

p j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(«0 -f Д u) [sin ô (/) — sin ô0] ^

Yi(Pk) e *

 

(4-55)

 

 

 

 

 

 

k= 1 PkZ (Pk)

 

 

 

где угол 6(1 — г) за

время 0 < т < / изменяется очень мало,

если, как будет видно дальше, само /-^(0,1

сек.

Поэтому

квадратную скобку

|sin Ô ( / — т) — sin ô0]

считаем

постоян­

ной

выносим за знак интеграла и, кроме того,

в

силу

тех

же причин

заменяем ô(t — т)

через

б(/).

 

 

 

 

 

Точно так же имеем:

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I - 1|

(и, +

Д и) Ш

- F3{p) j =

(«0 +

Ди) J * A ( r ) M t - x ) d x ^

 

 

 

 

 

 

 

О

 

 

 

 

 

 

 

 

 

 

 

 

5

 

P I

 

 

 

~

(«о + А и) [cos Ô [t) — cos ôe] V

 

e *

. /4.56J

 

 

 

 

 

 

 

PkZ (Pft)

'

Аналогично поступим в выражениях для Д /, и ДЧ^,. С учетом сказанного перейдем в равенсгоах (4-48) — (4-50) к оригиналам:

 

 

 

 

 

 

 

 

 

 

= 1

 

 

'о +

 

 

 

 

 

 

5

 

V

 

 

Д и) [sin ô (t) — sin ô0] V

—l^Ü S ------- Д и cos ô0 x

 

(«<

 

 

 

 

 

 

L l

pt Z (pt)

 

 

 

 

 

 

 

 

 

k= i

 

 

 

X ГM L

+

\ 1

:(P*)

e *

— («о +

Д «) [cos ô (/) — cos ô0] x

PtO

 

LZ(0)

^

L

PkZ’(Pt)

 

 

 

 

 

 

 

 

 

к = 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

pJ

 

 

 

 

 

 

 

X

S

Y'ÀPk) e * .

 

H-57)

 

 

 

 

 

Т Т 7 Г Г •

 

 

 

 

 

 

 

k= 1

Pt Z '(рк)

 

 

-

‘« +

à

l , -

v +

4 0 Sta«. [ - S f - + X

/ * ' ] +

 

 

 

 

 

 

 

 

 

к= 1

 

 

+

(«о +

A «) [sin ô (t) -

sin ô0] J ]

Ks'Pi) ***

-

 

 

 

 

 

 

 

 

 

* = i

PkZ (Pt)

 

- A U COSÔ0 [

m

+

j

^

 

/ * '] - < « . +

4 U) x

 

 

 

 

 

 

k= 1

 

 

 

 

 

 

 

 

 

 

 

 

5

 

V

(4-58)

 

 

 

X [cos ô (*) — cos ô0' V

Yj(p-~ e— ;

 

 

 

 

 

 

 

 

*=1 p*Z ipè

v

 

♦rf =

^ 0 + Д ♦ |- Ч > л +

Awsinôo

+

 

 

5

 

 

 

 

 

 

 

 

 

 

+

S

 

 

 

]

+ ("o +

A U) [sin ô

(/) - sin Ô01 x

 

к—1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

Ü x Pt £ (Pt)

+