- •Методическое пособие к курсовой работе по дисциплине «Математические задачи энергетики»
- •Содержание
- •Уравнения установившихся режимов электрических систем
- •1.1 Понятие о режимах электрических систем и схемах замещения
- •1.2 Аналитическое представление информации о конфигурации электрической сети с помощью матриц инциденций и матричное выражение законов Кирхгофа
- •Первая матрица инциденций «узлы-ветви» и ее применение для записи 1-го закона Кирхгофа
- •Вопросы для самопроверки:
- •1.2.2 Вторая матрица инциденций «ветви-контуры» и матричная запись второго закона Кирхгофа
- •1.2.3 Запись уравнений состояния сети по законам Кирхгофа
- •1.3 Метод уравнений узловых напряжений
- •1.3.1 Вывод узловых уравнений
- •Здесь [м]т – транспонированная 1-я матрица инциденций,
- •1.3.2 Определение матрицы узловых проводимостей и ее характеристика
- •1.4 Контурные уравнения установившихся режимов электрических систем
- •Запись уравнений состояния сети с помощью матриц обобщенных параметров.
- •Вопросы для самопроверки
- •1.6 Расчёт режима электрической сети с использованием матрицы коэффициентов распределения
- •Расчётные токи в узлах сети можно определить как:
- •2. Методы решения уравнений установившихся режимов электрических систем
- •2.1 Итерационные методы решения систем уравнений
- •2.2 Критерии сходимости итерации и анализ их выполнения для узловых уравнений установившихся режимов
- •2.2.1 Теорема сходимости итерации
- •2.2.2 Факторы, влияющие на сходимость итерации для узловых уравнений установившихся режимов
- •2.2.3 Критерии и анализ сходимости итерации для нелинейных систем узловых уравнений установившихся режимов
- •2.3 Решение уравнений узловых напряжений итерационными методами
- •2.3.1 Решение уравнений узловых напряжений в форме баланса токов
- •2.3.2 Обращенная форма уравнений узловых напряжений и их анализ
- •2.4 Применение метода Ньютона для решения для нахождения корней уравнений установившихся режимов
- •2.4.1 Обоснование метода Ньютона для решения нелинейного уравнения
- •2.4.2 Применение метода Ньютона для систем нелинейных уравнений
- •2.4.3 Решение нелинейных узловых уравнений методом Ньютона.
- •III. Задание на курсовую работу
- •Содержание расчетно-пояснительной записки (перечень подлежащих разработке вопросов)
- •Перечень графического материала (в виде компьютерных рисунков в формате а4)
- •IV. Примеры для выполнения разделов курсовой работы
- •С бу оставляем граф-схему замещения электрической сети и нумеруем её ветви и узлы (ребра и вершины) в соответствии с принципом ярусности:
- •Составление элементарных матриц параметров режима [pу], параметров сети [dZв],[dYв] и матриц соединений [м] и [n].
- •Расчёт матрицы узловых проводимостей [Yy] и матрицы контурных сопротивлений [Zk]. Расчёт матрицы узловых проводимостей [Yy] (См) (без учёта балансирующего узла) производим по формуле:
- •Расчет режима электрической сети по узловым уравнениям путем обращения матрицы узловых проводимостей
- •Расчет режима электрической сети на основе линейных контурных уравнений
- •Решение нелинейных обращенных узловых уравнений с матрицей- методом простой итерации.
- •Затем по выражению (8) проверяется точность вычислений:
- •Пример расчета:
- •Третья итерация:
- •Пример расчета:
- •В общем виде итерационный процесс можно записать в виде
- •Третья итерация:
- •Итерационный процесс закончен!
- •Заключение Литература
Составление элементарных матриц параметров режима [pу], параметров сети [dZв],[dYв] и матриц соединений [м] и [n].
Приложим рассчитанные нагрузки в узлах сети и запишем матрицу независимых параметров режима Ру (МВ
ТМ(МВт)
(1)
Составим диагональную матрицу сопротивлений ветвей [dZв] (Ом) по также уже известным сопротивлениям:
(Ом)
(2)
Найдём матрицу проводимостей ветвей [dYв] (См) по формуле:
[dYв]=[dZв]-1
С(Ом)
(3)
Первая матрица инциденций (матрица соединений ветвей в узлах):
MΣ
=
,
где [M] —1 матрица соединений для ветвей дерева;
[M] 1 матрица соединений для хорд.
Вторая матрица инциденций (матрица контуры-ветви):
N = [N N],
где [N] — блок матрицы N для ветвей дерева;
[N] — блок матрицы N для хорд.
Рис. 2
Расчёт матрицы узловых проводимостей [Yy] и матрицы контурных сопротивлений [Zk]. Расчёт матрицы узловых проводимостей [Yy] (См) (без учёта балансирующего узла) производим по формуле:
[Yy] =[M]*[dZB]-1*[M]T
(5)
Матрица узловых проводимостей, вычисленная с помощью программы Mathcad, совпадает с матрицей, вычисленной вручную.
Матрица узловых проводимостей [YΣ] (с учетом строки для балансирующего узла) определяется по формуле:
[YΣ] =[MΣ]*[dZB]-1*[MΣ]T
(См)
YΣ=
Матрица [YΣ] является вырожденной матрицей, т. е. обратной для неё не существует. Это подтверждается тем, что при суммировании элементов строк YΣ получается нулевая строка, и, следовательно, определитель этой матрицы, вычисленный по теореме разложения определителя по элементам строки (столбца), обращается в 0, т.е. det YΣ = 0
Расчёт матрицы контурных сопротивлений [ZК] (Ом) производим по формуле:
[Zк]=[N]*[dZВ]*[N]Т;
(Ом)
Как и ожидалось, получили симметрическую матрицу 3-го порядка.
В матрице Zк получили два отрицательных элемента (взаимные сопротивления 1-ого и 3-ого контуров). Знак «—» объясняется тем, что направления обхода контуров 1 и 3, совпадающие с направлениями хорд I и III, противоположны.
Диагональные элементы матрицы Zii (40.96, 53.76, 51.2) имеют преобладающие в своих строках и столбцах значения, поскольку представляют собой полные сопротивления контуров, а побочные элементы равны только суммам сопротивлений ветвей дерева, общих для i-го и j-го контуров. Поскольку каждый контур включает в себя и сопротивление замыкающей его хорды, то Zii > Zij. Это можно увидеть на схеме рис.
Заметим, что для матрицы контурных сопротивлений Zк отсутствует прямая числовая связь между диагональными и суммой соответствующих побочных элементов строки или столбца, свойственная матрице YУ, но в ней отражено реально существующее для схем соотношение Zii > Zij. При другом способе формирования системы независимых контуров может быть обеспечено соотношение:
Zii > Zij и
Zii > Zij(cумма)
2. Расчёт режима электрической сети при задании нагрузок в токах приведем на примере работы студента Кабанова Павла
Рис.1
Исходные данные для этой схемы следующие:
Li=35 км;
Si=S0(1+0,1*i) МВа;
S0=24 МВа;
Uбаз=1,05Uном ; Uном=110 кВ;
X0=0,4 Ом/км;
L5=0,65Li