Файл: Голембо, З. Б. Алгоритмизация и программирование электротехнических задач на электронных цифровых вычислительных машинах учеб. пособие.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 21.10.2024

Просмотров: 85

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

линии можно изменить направляемый по ней поток реактивной мощ­ ности, при этом изменяются потери активной мощности, связанные с перетоками по линии реактивной мощности.

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

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

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

денными в

гл. 6.

 

 

 

Для цепи, имеющей п узлов, может быть составлено п уравне­

ний, связывающих 2п переменных [Х^

и (5 г ) . Для

решения этой

системы в

каждом i-м узле дополнительно задают 5 г

или

Xt или

и

действительную или мнимую

часть мощности.

Уравне­

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

являются

нелинейными.

 

 

 

 

 

i задается

 

 

Для

рассматриваемого случая, в узле

и ве­

личина

Xi

берется равной |Хг +1|> а

начальное значение Sf выби­

рается таким, чтобы ReSj имело- заданное значение. Для

получения

очередных

значений

Х г

и Sj

сначала

определяется

St*:

 

St

= -

(Xt

+ 1) [ZYikXk

-

Xt

{Wik

+

27 i ; )

-

2K„.].

В результате

вычисляется

 

 

 

 

 

 

 

 

 

 

 

 

 

Re St

= ReS;;

\mSt

=

Im S] .

 

 

 

 

Затем определяется

значение

 

Xt*:

 

 

 

 

 

 

.

 

Р ;

Л + ( А т ) - ^

 

 

 

~ .

 

i x + i i

Y _ *

 

\ Л ; +1

/

;

 

у

 

/У*

I

п

\Л1

+

1

'"

 

 

YY-ь + ЪУ--

 

' A

i

=

( A ' + l>

y* i 1 I — 1-

 

 

 

A

J

^

 

 

 

 

 

 

 

\ x i

+

1 \

Следовательно,

 

= | Х г + 1|.

 

 

 

 

пересчет: один раз

Для

указанного

режима

требуется двойной

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

Полные мощности потерь AW в цепи

125


Величина небаланса H=HSi—AW служит для контроля вы­ числений и оценки точности решения.

Для узлов замещенной длинной линии электропередачи, в ко­

торых задано значение мощности в точках примыкания к

системам

и требуется определить напряжение, пользуясь методом

последо­

вательных

приближений, очередное приближение определяется

по формуле:

 

 

 

 

 

S K i f t X A - S 7 i ; - + S i / ( X £

+ l)

 

xt =

nYlk

+ i:Ytl

;(s = s,-),

где Sj—заданная величина мощности.

 

Однако

для указанного

режима при использовании, например,

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

Итерационный процесс, по Зейделю (при рассмотрении, напри­ мер, дальних линий передач в зависимости от параметров передачи

 

U„

и,

 

 

 

 

и,

о?

 

 

 

 

 

 

 

4=У"«

/ =

/ Г -Ук

/ 1

/

±Уг

fllift

 

 

 

 

 

- -

 

 

М

м

 

 

 

 

 

 

1 -

 

 

 

 

 

 

Рис. 7.12. Схема

длинной

линии

с отборами

мощности

 

и степени компенсации),

не всегда сходится и, следовательно, не

всегда можно

получить

искомое

решение.

 

 

 

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

входе и выходе четырехполюсника приводит к простым

выражени­

ям, широко используемым при расчете длинных линий

электропе­

редачи.

 

Для баланса потоков активной и реактивной мощностей приме­

нительно к k - щ узлу

(рис. 7.12)

справедливы уравнения

 

f L i + Qf-i Г, _ D _ д р .

Г к — Г/1-1

[12, ,

Кк

ГШ

1ЛГ k ,

126


где

Uh\ Uh_i—модули

 

напряжений в узлах

k и /г—1; Pk,

Pf c _i,

Qft>

Qh-i — активные

и реактивные составляющие

мощности

в

узлах k и /г—1; Poh

— активная мощность,

отбираемая в

узле;

APh

— активная мощность потерь в компенсирующих

устройствах;

Rh> Ха—активное

и

индуктивное сопротивление k

звена;

Yh

проводимость линии

на землю в узле k.

 

 

 

 

Для напряжения на выходе (применительно к узлу k) в общем

случае можно написать следующую функциональную

зависимость:

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

PkRk + QkXk

1 PkXk + QkRk

В зависимости от требований задаются те или иные граничные условия. Поэтому решение происходит по разным формулам. Рас­ смотрим основные случаи направления решения в зависимости от начальных условий:

1)заданы подведенная на вход мощность активной и реактивной составляющих W=P+jQ и модуль напряжения U на входе. Не­ обходимо найти напряжения в последующем узле и определить активные и реактивные составляющие мощности потерь;

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

тальных

узлах:

U3

Un,

лежащие

в определенных

допустимых

пределах

Umax

<J U <^ ит--

Надо найти величину

реактивной

мощности, обеспечивающей

указанные

выше условия.

Для узлов,

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

вида AQ2+BQ+C=0,

где А,

В, С—действительные

величины.

Параметры длинной

линии

электропередачи Rh

и Xk опреде­

ляются по алгоритму,

приведенному в настоящей главе.

При решении квадратного уравнения практический интерес представляет тот корень, который связан с устойчивой характе­ ристикой работы системы, т. е. для которого выполняется условие Re(Vf t + 1 )> 0. При таком условии логическое решение задачи сво­ дится к исследованию знака квадратичной формы.

в. Реализация алгоритма на ЭЦВМ

Принципиальная блок-схема программы расчета распределения реактивных мощностей длинных линий электропередач (рис. 7.13) состоит из отдельных арифметических и логических блоков. Неиз-

127


 

Ввод

исходной

 

информации

 

 

Образование команд,

зависящих

от

числа

п

 

 

 

(участков

линии)

 

 

Вычисление

 

yi,Zc

их

модулей

и

аргументов

Вычисление

 

и печать

 

элементов

схемы

 

 

замещение

четырехполюсника

 

 

 

 

 

 

X

 

 

 

 

 

 

Проверка

окончания

подсчета

параметров

 

 

 

Задание

/}0,\иг1

 

 

 

 

Вычисление

U2,

 

 

 

Вычисление

 

Qt,

по заданным

 

Q0,U1,P1

 

по

заданным\Ui\,\U2\fi

 

 

 

Печать

 

 

 

 

 

 

 

Вычисление

UK+u по заданным

UK,PK,QK

 

 

 

Щуп задания

 

UK+1

 

 

 

 

 

Контроль

 

 

 

 

Вычисление

по

 

min\U\u\UK+,\£max\L

 

 

заданным^,

1<иЛ

Печать к-го

участка'

 

ВычислениеAQuP=PK-P'/o\AQ\

 

при&о>оилиР=Рк-д%Щ!№

 

 

 

 

 

Проверка

окончания

\-*-\ Вычисление

 

QK+1,

Q'M,

PK+i

Вычисление

и

печать

 

 

Щуп отбора

Рк+1

 

Qn+U

Q'nM,PnU,AP

 

 

 

 

 

 

 

 

 

 

 

.СТОП

 

 

 

 

Отбор

 

Рк+1

 

 

Рис. 7.13 . Блок-схема программы расчета распределе­ ния реактивных мощностей дальних линий электропе­ редач

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

128


включают в состав общей программы. В логических блоках машина;

избирает путь дальнейшего

расчета.

xkr

Программа

построена

таким

образом, что все параметры

rh, yh,

IIh , Ph,

Qft, входящие в

расчетные формулы, хранятся

в

рабочих

ячейках памяти

машины.

 

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

сылка величин

Ui, Pi, xi,

п,

у2

в рабочие ячейки. Затем вычисля­

ется величина

Qi' = Ui2yi.

 

 

 

 

В связи с тем, что граничные условия первого участка могут быть различными (на входе задаются или Q0 и Ux, или Ui и £/,)>. и машина должна выбирать один из двух возможных случаев расчета.. Выбор осуществляется при помощи логической шкалы напряжений_

Случай 1. Пусть задан модуль напряжения U2 во втором

узле;

В этом случае по заданным Ui и (У2 определяется величина

реак­

тивной мощности на входе первого участка. Для этого уравнениеразрешается относительно Q (берется то решение, для которого» ReU,,+1> 0). Далее вычисляется угол <р между векторами напряже­ ний в первом и втором узлах, после чего выводятся на печать номеррассчитываемого участка, модуль напряжения U2, ф, Qi, Qi' и Pi_ Затем машина переходит к расчету следующего участка. В случае,,

если

для обоих

решений

R&Uh+1<i 0, машина останавливается,,

такой

режим

неосуществим.

 

 

Случай 2.

Пусть задана

Q0 — реактивная

мощность на

входе

линии. Вначале

вычисляем

Q i = Q 0 + Q i . Затем

определяем

напря­

жение во втором узле и угол ср. Дальнейший расчет протекает в. зависимости от того, выполняется ли неравенство Um-in<z U2 <.ilmsx_ или не выполняется. В том случае, когда неравенство удовлетво­ ряется, на печать выводится номер участков линии U2, q>i, Qi', Р * и машина переходит к расчету следующего участка. При невыпол­ нении неравенства расчет протекает более сложно.

Когда U2<.Um\n

(или U2>Umax),

на

печать выводится

номер

участка

линий U2,

ср, Qi,, затем

вторично производится

расчет-

первого

участка. При этом перерасчете

определяется величина Qb

на входе первого участка при 02т-ш-

Вычисляется новое зна­

чение угла ср' и величина отбора реактивной мощности A Q = Qi—Qi._

Мощность активных потерь в устройствах отбора первого узла, вычисляется по формулам A P = | A Q | Р, если A Q > 0 и A P = | A Q | если AQ < 0 .

Затем выводятся на печать величина отбора реактивной мощ­

ности AQ и новое значение

активной

мощности

с учетом

потерь.

( P i = P i — А Р ) , номер участка

линии и

минимум

11, ср, Qi,

Qi', Pi ..

После этого первый участок рассчитывается в третий раз; при этом ^2 = i 7 m i n ; P = P i . Определяется новое значение реактивной мощ­

ности на его входе. На

печать вновь выводится номер участка л и ­

нии U2, ф 1 , Qi, Qi', Pi'

и машина переходит к расчету следующего*

участка.

 

Случай 3. Рассмотрим расчет /е-го участка линии электропе­ редачи. В зависимости от того, задано ли напряжение в (&+1)-м уз-

129>