Переходные процессы / Лабораторные работы / Лабораторная работа 14
.8.pdfМинистерство образования Российской Федерации
КУБАНСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНОЛОГИЧЕСКИЙ УНИВЕРСИТЕТ
Кафедра электроснабжения промышленных предприятий
ОПИСАНИЕ лабораторной работы № 14.8
по курсу «Переходные процессы в электрических системах»
ИССЛЕДОВАНИЕ УСТОЙЧИВОСТИ ПРОСТЕЙШЕЙ ЭЛЕКТРИЧЕСКОЙ СИСТЕМЫ
Краснодар 2003
УДК 621.311.1.001.45(076.5)
Печатается то решению Редакционно-издательского совета университета
ИССЛЕДОВАНИЕ УСТОЙЧИВОСТИ ПРОСТЕЙШЕЙ ЭЛЕКТРИЧЕСКОЙ СИСТЕМЫ
Лабораторная работа № 14.8 по курсу «Переходные процессы в электрических
системах»
Краснодар: изд. КубГТУ, 2003. – 18 с.
Исследуется динамическая устойчивость простейшей электрической системы. Рассматривается изменение запаса устойчивости в зависимости от различных факторов: наличия и отсутствия демпфирования, характера демпфирования, изменения времени срабатывания защиты, изменения загрузки турбины.
Составители: |
Б.А. Коробейников |
|
Л.И. Деревцова |
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 |
|