- •Методическое пособие к курсовой работе по дисциплине «Математические задачи энергетики»
- •Содержание
- •Уравнения установившихся режимов электрических систем
- •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) проверяется точность вычислений:
- •Пример расчета:
- •Третья итерация:
- •Пример расчета:
- •В общем виде итерационный процесс можно записать в виде
- •Третья итерация:
- •Итерационный процесс закончен!
- •Заключение Литература
Запись уравнений состояния сети с помощью матриц обобщенных параметров.
Обратимся к уравнениям состояния сети (54), (55) по законам Кирхгофа.
Представим матрицу коэффициентов системы [A-1] в виде блочной матрицы с размерностью блоков по числу узлов n и числу контуров k
(60)
Тогда
(61)
Здесь — квадратная матрица, называемая матрицей входных и взаимных проводимостей ветвей схемы. Её элементы yij определяют величину и фазу тока в i-ой ветви от действия ЭДС j-ой ветви и называются взаимными проводимостями, а элементы уii определяют величину и фазу тока в i-ой ветви от действия ЭДС этой же ветви и называются собственными или входными проводимостями.
При отсутствии ЭДС ветвей ([Eв] = 0), выражение (61) обращается в
, (62)
откуда наглядно виден смысл матрицы [С] и её элементов.
[C] — матрица порядка , называется матрицей коэффициентов распределения задающих токов узлов по ветвям сети. Её произвольный элемент cij представляет собой долю от тока j-го узла в токе i-ой ветви:
, (63)
где cij — элемент матрицы коэффициентов распределения [C];
Ii — ток i-ой ветви;
Jj — задающий ток j-го узла.
Матрицы [С], [D] и, следовательно — [Yв], вычисляются путем обращения матрицы [А] с помощью разбиения на блоки, и представляют собой линейные комбинации блоков матрицы [А], показанных в (53). При этом [С] и [Yв] могут быть выражены как на основе узловой модели сети:
, (64)
где — обратная к матрице узловых собственных и взаимных проводимостей,
так и на основе контурной модели:
(65)
где — обратная к матрице контурных сопротивлений.
Выражения (64), (65) показывают, что процедура нахождения [С] достаточно громоздкая, но вычисленная один раз, эта матрица позволяет вести многократные расчёты режима по выражениям (55), (61) вручную или на ЭВМ с высоким быстродействием. После нахождения токов ветвей остальные параметры режима рассчитываются по известным формулам.
Правильность выражений (64), (65) и результатов конкретных вычислений матрицы [С] по этим выражениям можно проверить по соотношению:
, (66)
которое получается если в 1-ый закон Кирхгофа подставить вектор токов ветвей из (55):
, а следовательно
(67)
Логика выражения (60) наглядно видна, если в выражение (55) подставить [С] из (60):
(68)
[U]
[Uв]
[I]=[I]
тождество
Аналогично можно показать логику выражения для матрицы [С] на основе контурной модели сети:
Перемножим в правой части и сгруппируем сомножители скобками
[Iα]
[ΔUα]
[Iα΄]
[Iв]
[Ек]
[Iβ]
[Iв΄΄]=[Iα΄΄
Iβ]Т
[Iв]
[Iв]
=
Здесь [Iα΄] – составляющая токов в дереве сети, обусловленная задающими токами при отсутствии хорд;
[Iα΄΄] – составляющая токов в дереве сети, вызванная замыканием хорд.
[Iα]=[Iα΄]+[Iα΄΄]; (69)
[Iв]= [Iα΄]+[Iв΄΄] (70)