Решение систем нелинейных уравнений установившегося режима методом Зейделя

Суть метода - найденное на текущей (k+1)-й итерации приближение напряжения узла номер (i-1) Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru сразу же используется для нахождения (к+1)-го приближения напряжения следующего i-го узла Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru .

Нелинейное уравнение установившегося режима в форме баланса токов для i-го узла имеет вид:

Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru (1)

Левая часть уравнения – линейна, правая – не линейна. Т.е. уравнение нели-нейно.

Преобразуем это уравнение к виду, удобному для решения итераци-онными методами, т.е. решим его относительно Ui :

Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru (2)

Запишем это уравнение в итерационной форме, т.е. с указанием номеров при-ближений:

Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru (3)

Уравнение (3) соответствует вычислительной схеме метода простой итера-ции.

Если при вычислении Ui(к+1) в правую часть выражения (3) подставлять найденные на текущей итерации приближения Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru , то получим выражение, соответствующее вычислительной схеме метода Зей-деля:

Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru (4)

Для сети, состоящей из n узлов, нужно записать систему n таких уравне-ний, которая в итерационной форме, соответствующей методу Зейделя, имеет вид:

Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru (5)

Алгоритм решения СНАУ установившегося режима методом

Зейделя

1. Задание условий и параметров расчета:

точность расчета E, предельное количество итераций nпр, счетчик

итераций к=0;

2. Задание начальных приближений неизвестных:

принимаем напряжение U(0)=Uном.

3. Выполнение итерации расчета для всех узлов, кроме опорного, в соот-

ветствии с (5).

При расчете U1(к+1) в правую часть 1-го уравнения подставляем к-е

приближения всех неизвестных. При расчете U2(к+1) в правую часть

2-го уравнения подставляем только что найденное значение U1(к+1) и

к-е приближения остальных неизвестных и т.д. При расчёте послед-

него напряжения Un(к+1) в правую часть последнего уравнения подстав

ляем найденные ранее на этой итерации (к+1)-е приближения осталь-

ных напряжений;

4. Проверка завершения итерационного процесса в соответствии с усло-

вием:

Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru (6)

Если это условие не выполняется – возврат к п.3. и повторение расчета

при новых приближениях неизвестных.

В модели реальной электрической сети могут присутствовать специ-альные узлы, например, узлы с заданным модулем напряжения (узлы с фик-сацией модуля напряжения ФМ). В таких узлах заданными параметрами яв-ляются Uз (модуль напряжения) и Pз (активная мощность). Искомыми явля-ются Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru (угол напряжения) и Q (реактивная мощность), либо вместо Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru могут быть составляющие напряжения - Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru

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

Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru (7)

где Uiз -заданный модуль напряжения в i-ом узле,

Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru - искомые составляющие напряжения в i-ом узле.

При этом должно выполняться условие:

Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru (8)

Для того чтобы это условие выполнялось на каждом шаге итерационно-го процесса, составляющие напряжения корректируют в соответствии со следующими формулами:

Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru (9)

Итерационный процесс заканчивается, когда выполняется условие (6) для всех узлов:

Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru

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

Сходимость метода Зейделя

Анализируя сходимость итерационного метода – рассматриваем, прежде всего, скорость сходимости (необходимое количество итераций) и характер сходимости (колебательная или экспоненциальная).

Экспериментально выявлен ряд факторов, влияющих на скорость сходи-мости. Сходимость улучшается при:

1) более точном задании начальных приближений по напряжению;

1)2) разгрузке сильно загруженных линий;

1)3) увеличении числа генерирующих узлов с заданным модулем напря-жения (узлы ФМ) и широким диапазоном изменения реактивной мощности;

1)4) увеличение числа контуров в сети;

1)5) увеличение числа связей БП с остальными узлами;

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

Сходимость метода Зейделя для СНАУ установившегося режима до-вольно медленная (для сети 100-200 узлов требуется около 300-500 итера-ций).

Сходимость может быть улучшена с помощью коэффициентов ускоре-ния (метод неполной релаксации).

Пусть Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru - значения переменной x на Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru итерациях. Тогда (к+1 )- е приближение можно определить по формуле:

Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru (10)

где Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru -коэффициент ускорения;

Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru приближения, вычисленные без учета Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru и с учетом Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru .

При Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru .

Значения коэффициента Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru выбирается в зависимости от характера сходимости итерационного процесса. Если сходимость монотонная, то:

Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru >0 (11)

Если сходимость монотонная, то принимают Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru >1. Использование ускоряющего коэффициента Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru на каждой итерации позволяет приблизить очередное значение к точному решению. Этим ускоряется сходимость. Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru

Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru

Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru Если характер сходимости колебательный, то выражение (11) будет меньше нуля и значение коэффициента принимают Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru <1 .

В этом случае коэффициент является демпфирующим – он позволяет уменьшить ам-плитуду отклонений приближений и умень-шить количество итераций.

Возможные причины отказа сходимости итерационного процесса

1. Не правильное задание исходных данных (параметры схемы, задан-ные параметры режима)

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

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

Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru

Пример:

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

Запишем для заданной схемы уравне-ния установившегося режима в форме ба-ланса токов вида:

Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru :

Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru

Так как балансирующий узел (0) связан только с 1-й узлом, то перено-сим соответствующий элемент 1-го уравнения в правую часть:

Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru

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

Решение систем нелинейных уравнений установившегося режима методом Зейделя - student2.ru

Наши рекомендации