Техническая постановка задачи расчета и анализа

Стр 1 из 4Следующая ⇒

СОДЕРЖАНИЕ

ВВЕДЕНИЕ................... 4

1.ПРОГРАММА КУРСА................ 4

1.1 Техническая постановка задачи расчета и анализа

установившихся режимов электрических систем...... 4

1.2 Уравнения состояния электрических систем...... 5

1.3 Методы решения уравнений состояния электрических систем. 5

1.4 Анализ статической устойчивости электрических систем... 6

2.ОСНОВЫ РАСЧЕТА И АНАЛИЗА УСТАНОВИВШИХСЯ РЕЖИМОВ

ЭЛЕКТРИЧЕСКИХ СИСТЕМ............... 7

2.1 Техническая постановка задачи.......... 7

2.2 Расчет установившегося режима с использованием

линейных математических моделей......... 9

2.3 Уравнения состояния ЭС............. 12

2.4 Пример расчета на основе линейной модели....... 15

2.5 Реализация расчета режима в среде Mathcad.... 19

2.6 Задание № 1 по контрольной работе......... 23

3.МАТЕМАТИЧЕСКИЕ МЕТОДЫ АНАЛИЗА СТАТИЧЕСКОЙ УСТОЙЧИВОСТИ

УСТАНОВИВШИХСЯ РЕЖИМОВ............... 23

3.1 Техническая постановка задачи...........23

3.2 Пример анализа статической устойчивости по корням

характеристического уравнения........... 24

3.3 Пример анализа устойчивости по критерию Гурвица.....28

3.4 Пример использования критерия Михайлова

для анализа статической устойчивости........ 31

3.5 Реализация задачи анализа устойчивости в среде Mathcad..33

3.6 Задание №2 для контрольной работы......... 35

Литература.................. 36

ПРИЛОЖЕНИЕ.................. 36

ВВЕДЕНИЕ

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

В курсе “Математические задачи электроэнергетики” будут представлены основные научные достижения в области управления таким сложным объектом как электроэнергетическая система (ЭЭС), рассмотрены возможности современной математики и вычислительной техники, позволяющие смоделировать реальные процессы, происходящие в ЭЭС.

В пособии представлены программа курса с ссылками на литературу, задание на выполнение контрольной работы для студентов заочного отделения, методические указания для проведения расчетов вручную и с использованием ЭВМ. Материал, изложенный в пособии может быть использован также студентами дневного отделения ФЭН при подготовке в практическим занятиям по курсу.

ПРОГРАММА КУРСА.

Техническая постановка задачи расчета и анализа

Установившихся режимов электрических систем.

Электрическая система (ЭС) как объект математического моделирования. Понятие режима работы ЭС. Виды режимов. Параметры режима функционирования ЭС. Общая характеристика разделов прикладной математики, используемых при решении задачи расчета установившихся режимов ЭС ([1]с.7-30, [2]c.5 - 6).

Уравнения состояния электрических систем

1.2.1 Понятие схемы замещения электрической системы. Схемы замещения источников энергии, потребителей и элементов электрической сети. Пример перехода от реальной схемы электрической системы к схеме замещения. Моделирование электрической сети с помощью направленного графа ([1]c.31-33, 37-38, [2]c.6-10, [3]c.109-111).

1.2.2. Использование матричных методов прикладной математики для моделирования процессов, происходящих в электрической системе. Основы матричной алгебры. Действия с матрицами. Виды матриц, используемых при расчете установившихся режимов. Матрицы инциденций первого и второго рода. Правила формирования матриц инциденций, исходя из структуры электрической сети, представленной в виде графа. Матрицы режимных параметров ([5]c.35-52, [1]c.38-40, [2]c.10-13, [3]c.114-118).

1.2.3 Виды уравнений состояния электрической системы. Представление в матричной форме основных законов электротехники: закона Ома, первого и второго законов Кирхгофа. Получение обобщенного уравнения состояния на основе двух законов Кирхгофа. Уравнения узловых напряжений. Структура и физический смысл элементов матрицы узловых проводимостей. Контурные уравнения состояния ЭС. Преимущества и недостатки различных форм представления уравнений состояния с учетом удобства реализации алгоритмов на ЭВМ ([1] c.40-42, c.48-56, [2] c.13-30, [3] c.118-129, [4]c.25-35).

Методы решения уравнений состояния электрических систем.

1.3.1. Математическая модель задачи расчета и анализа установившихся режимов ЭС. Способы представления параметров генераторных и нагрузочных узлов схемы замещения электрической системы. Классификация методов решения систем уравнений состояния ([1]c.70-74, [2]c.30-33, [4]c.8-14, c.35)

1.3.2. Характеристика прямых методов решения уравнений состояния, представленных в виде систем линейных алгебраических уравнений. Алгоритм метода Гаусса с обратным ходом и без обратного хода. Анализ точности полученного решения по сумме невязок. Факторы, влияющие на точность решения систем уравнений методом Гаусса. Применение метода Гаусса для решения системы линейных уравнений узловых напряжений.([1] c.79-88, [2] c.34-36, [4] c.35-39).

1.3.3. Общая характеристика итерационных методов решения систем уравнений. Понятие итерационного процесса. Виды итерационных процессов. Критерии сходимости итерационных процессов. Метод простой итерации. Метод Зейделя и его применение для решения систем линейных уравнений узловых напряжений. Сравнительная характеристика итерационных методов с учетом их реализации на ЭВМ ([1] c.91-101, [2]c.37-38, [4] c.43-50).

1.3.4. Нелинейные уравнения установившегося режима. Представление нелинейных уравнений состояния в форме баланса токов и в форме баланса мощностей. Использование метода Зейделя для решения систем нелинейных уравнений состояния. Метод Ньютона. Графическая интерпретация метода Ньютона для функции одной переменной. Алгоритм метода Ньютона для функции многих переменных([4]c.57-60, c.64-76).

Анализ статической устойчивости электрических систем.

1.4.1.Понятие устойчивости применительно к электрической системе. Состояние устойчивого равновесия. Дифференциальные уравнения, описывающие изменение параметров системы во времени. Понятие свободной и вынужденной составляющих. Определение статической устойчивости как устойчивости в малом. Метод малых колебаний ([1]с.173-175, [2] c.39-41, [6]c.189-190).

1.4.2 Анализ статической устойчивости по корням характеристического уравнения. Характеристическое уравнение системы. Теорема Ляпунова. Графики изменения свободной составляющей в зависимости от корней характеристического уравнения. Пример анализа статической устойчивости системы " станция - шины бесконечной мощности" при подключении к узлу синхронного неявнополюсного генератора ([1]c.175-181, [2]c.41-44, [6]c.190-194).

1.4.3 Анализ статической устойчивости на основе алгебраических критериев. Понятие критерия устойчивости. Алгебраические критерии устойчивости. Необходимые и достаточные условия устойчивости системы. Критерий устойчивости Гурвица. Алгоритм составления определителя Гурвица по характеристическому уравнению системы. Примеры использования критерия Гурвица для характеристических уравнений 1-го-4-го порядков ([1]с.181-188, [2]c.44-47, [6]c.195-198).

1.4.4 Анализ статической устойчивости на основе частотных критериев. Частотные критерии анализа устойчивости. Принцип аргумента. Понятие годографа характеристического уравнения. Критерий устойчивости Михайлова. Примеры типовых годографов для устойчивых, неустойчивых систем. Граница колебательной и апериодической устойчивости. Пример использования критерия Михайлова для характеристического уравнения 3-го порядка ([1] c.188-196, [2] c.47-49).

ОСНОВЫ РАСЧЕТА И АНАЛИЗА УСТАНОВИВШИХСЯ РЕЖИМОВ

ЭЛЕКТРИЧЕСКИХ СИСТЕМ.

Расчет установившегося режима с использованием

Уравнения состояния ЭС.

Основные формы уравнений состояния ЭЭС: обобщенное уравнение состояния, уравнения узловых напряжений, контурные уравнения подробно описаны в [1].

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

Обобщенное уравнение состояния в матричной форме имеет вид:

, (2.1)

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

- матрица инцидениций 1-го рода, предназначена для описания структурных связей узлов и ветвей в расчетной схеме. ( - количество узлов, - количество ветвей);

Структура:

Правило формирования: каждый элемент матрицы , располагается на пересечении строки (номер узла) и столбца (номер строки), его значение определяется следующим образом:

-1, если ветвь входит в узел

= 1, если ветвь выходит из узла

0, если ветвь не соединена с узлом .

-произведение двух матриц:

- структурная матрица инциденций второго рода, отражающая связь ветвей в независимые контуры .

Структура:

Правило формирования:

-1, если направление ветки

противоположно направлению обхода контура

1, если направление ветки совпадает с

направлением обхода контура

0, если ветвь не входит в контур

- матрица сопротивлений ветвей.

- объединенная матрица свободных членов, включающая в себя:

- вектор задающих токов; - вектор ЭДС контуров.

При использовании обобщенного уравнения состояния расчет установившегося режима ЭЭС производится в следующем порядке: вначале определяются токи в ветвях схемы , а затем рассчитываются падения напряжения в ветвях , напряжения в узлах , потоки активной и реактивной мощностей и т.д. Пример расчета приведен ниже при описании реализации в среде Mathcad.

Общий вид уравнения узловых напряжений [1, 2]:

, (2.2)

где - матрица узловых проводимостей;

- матрица проводимостей ветвей, обратная матрице сопротивлений ветвей ;

- матрица узловых напряжений;

- базисное напряжение балансирующего узла.

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

, (2.3)

Матрица узловых проводимостей рассчитывается по формуле:

, (2.4)

где

- транспонированная матрица инциденций первого рода;

— матрица узловых проводимостей.

Структура определяется физическим смыслом ее элементов:

· на главной диагонали расположены собственные проводимости узлов , равные сумме проводимостей ветвей, соединенных с узлом ;

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

Матрица является симметричной и слабо заполненной, т.е. содержит большое число нулевых элементов. Эти свойства позволяют реализовать на ЭВМ эффективные алгоритмы расчета режимов с учетом слабой заполненности массивов.

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

Пример расчета на основе линейной модели.

Условие задачи: Для расчетной схемы, представленной на рис. 2.5 записать матричное уравнение узловых напряжений и рассчитать значения узловых напряжений методом Гаусса.

Исходные данные:

- сопротивления

ветвей;

- задающие токи, моделирующие

подключение нагрузки.

Расчет начинается с формирования уравнения состояния по расчетной схеме:

1. Составим матрицу инциденций 1-го рода.

1 2 3 4 5 6

При правильном составлении матрицы М строка, соответствующая балансирующему узлу, дополняет каждый столбец до нуля.

2.Составим транспонированную матрицу

3.Определяем матрицу узловых проводимостей

В матричной форме уравнение узловых напряжений имеет вид:

(2.5)

5. Перейдем к системе уравнений:

(2.6)

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

1. Решая систему уравнений вида (1.12), определяются значения узловых напряжений . Произведем расчет с помощью метода Гаусса.

Прямой ход Гаусса состоит из однотипных шагов, связанных с формированием из матрицы коэффициентов верхней треугольной матрицы.

Шаг 1. Получим первое ключевое уравнение, для чего разделим первое уравнение системы (2.5) на коэффициент при , а затем исключим из всех уравнений, расположенных ниже ключевого.

(2.7)

Шаг 2. Принимаем за ключевое второе уравнение (разделим все коэффициенты на ) и исключим из уравнений ниже ключевого.

Преобразованная система, начиная с ключевого уравнения имеет вид:

(2.8)

Шаг 3. Принимаем за ключевое третье уравнение и исключаем из всех уравнений ниже ключевого, преобразованная система, начиная с ключевого уравнения имеет вид:

(2.9)

Шаг 4. Выбираем четвертое ключевое уравнение:

(2.10)

Обратный ход Гаусса:

Анализ точности расчета: Производится расчет невязок по исходной системе уравнений:

(2.11)

2. Из уравнения связи параметров режима [ 1 ] находятся падения напряжений в ветвях

. (2.12)

3. Из уравнения закона Ома (1.1) определяются токи в ветвях схемы

. (2.13)

4.По известным значениям и определяются остальные параметры режима и т.д.

Лабораторная работа N 1

Матрица инцидений 1-го рода

Сопротивления ветвей

Произведение матриц

Матрица коэффициентов

Лабораторная работа N 2

Задание № 1 по контрольной работе.

Задана расчетная схема электрической системы, представленная в виде направленного графа, который содержит 5 узлов, 7 ветвей и 2 независимых контура (по вариантам в приложении 1). Направление ветвей и независимых контуров может быть задано произвольно. Для указанной схемы ЭС необходимо рассчитать параметры установившегося режима. В связи с этим, требуется выполнить следующие пункты задания: 1.Составить обобщенное уравнение состояния на основе первого и второго законов Кирхгофа, записать этого уравнение в матричной форме и в виде системы уравнений;

2.Вычислить матрицу узловых проводимостей и записать уравнение узловых напряжений в матричной форме и в виде системы уравнений;

3. Рассчитать узловые напряжения и токи в ветвях с использованием метода Гаусса с обратным ходом. Оценить точность полученных результатов. Исходные данные, необходимые для проведения расчетов, приведены в таблице 1 приложения

УСТАНОВИВШИХСЯ РЕЖИМОВ.

Пример анализа статической устойчивости по корням

Пример анализа устойчивости по критерию Гурвица.

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

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

Рассмотрим применение алгебраического критерия Гурвица для анализа статической устойчивости простейшей электрической системы: станция - шины бесконечной мощности, рассмотренной в разд. 3.2. При этом учтем не только демпфирующие моменты, но и переходные процессы в обмотке возбуждения генератора.

В этом случае характеристическое уравнение будет иметь третий порядок [4]

(3.10)

где - постоянная инерции генератора; -переходная постоянная времени генератора по продольной оси; - коэффициент демпфирования.

Значение коэффициента вычисляется по (3.4), а для определения используется выражение из [4]:

(3.11)

где - переходное реактивное сопротивление генератора по продольной оси;

Переходная постоянная времени генератора рассчитывается из выражения

(3.12)

где - постоянная времени обмотки возбуждения синхронной машины при разомкнутой обмотке статора. Произведем расчет коэффициентов характеристического уравнения (3.10), используя исходные данные примера разд. 3.2 и дополнительные справочные характеристики генератора в относительных единицах [7]:

= 7.26 = 0.172.

Находим по (3.12)

;

тогда

Найдем значение коэффициента по (3.11)

тогда .

Запишем характеристическое уравнение (3.10) с учетом значений коэффициентов:

(3.13)

Для использования Гурвица составим определитель Гурвица по следующим правилам:

- по главной диагонали располагаются коэффициенты уравнения (3.10) в порядке возрастания индексов, начиная с ;

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

- все недостающие коэффициенты заменяются нулями.

Определитель Гурвица для характеристического уравнения (3.13) имеет вид:

Выделим миноры относительно главной диагонали и применим критерий Гурвица: для устойчивости системы необходимо и достаточно, чтобы при > 0 все главные диагональные миноры определителя Гурвица были положительны.

= 3824.5> 0;

= = 448.2> 0;

;

Таким образом, рассмотренная электрическая система является статически устойчивой.

Пример использования критерия Михайлова

Лабораторная работа N 3

Given

Задание №2 для контрольной работы.

1. Проанализировать устойчивость системы по корням характеристического уравнения вида (3.2), считая, что генератор подключен к 4 узлу расчетной схемы. Параметры генератора одинаковы для всех вариантов и приведены в методических указаниях, кроме угла , заданного в приложении по вариантам.

2. Определить устойчивость системы на основе характеристического уравнения (3.10) по критерию Гурвица.

3. Провести анализ устойчивости системы по критерию Михайлова.

4. Подготовиться к выполнению лабораторной работы по разд. 3.5

Литература

1. Веников В.А. Электрические системы. Математические задачи энергетики: Учебник для студентов вузов.- М.: Высшая школа, 1981, 288 с.

2. Любченко В.Я., Манусов В.З. Физико-математические основы электроэнергетики: Учеб. пособие; В 2-ч., часть 1/Новосиб.гос.техн.унт. -Новосибирск, 1994. - 58 с.

3.Мельников Н.А.Электрические сети и системы.- М.: Энергия, 1969.- 456 с.

4. Идельчик В.И. Расчеты и оптимизация режимов электрических сетей и систем.- М.: Энергоатомиздат, 1988.- 288 с.

5. Данилина Н.И., Дубровская Н.С. и др. Численные методы.- М.: Высшая школа, 1976.- 368 с.

6. Веников В.А. Переходные электромеханические процессы в электрических системах.- М.: Высшая школа, 1970.- 472 с.

ПРИЛОЖЕНИЕ

Исходные данные для задания №1

№	вар	Сопротивления ветвей, Ом	Задающие токи, кА
вар	схемы	Z1	Z2	Z3	Z4	Z5	Z6	Z7	J1	J2	J3	J4	J5
		0.2	0.3	0.4	0.7	0.8	0.3	0.5
		0.1	0.5	0.7	0.3	0.6	0.4	0.3
		0.3	0.4	0.8	0.9	0.5	0.7	0.6
		0.2	0.5	0.7	0.9	0.5	0.7	0.3
		0.2	0.4	0.3	0.5	0.3	0.6	0.4
		0.5	0.3	0.6	0.4	0.9	0.7	0.8
		0.2	0.3	0.4	0.8	0.3	0.5	0.8
		0.3	0.5	0.4	0.5	0.6	0.9	0.8
		0.1	0.5	0.3	0.5	0.4	0.1	0.2
		0.2	0.3	0.5	0.6	0.8	0.2	0.3

Исходные данные для задания №2

Параметры генератора заданы для всех вариантов в разделе 3.2

№ вар	угол d	№ вар	угол d
	p/6		p/4
	p/3		p/6
	p/4		p/3
	p/6		p/4
	p/3		p/6

Варианты расчетных схем для практических занятий

Вариант 1

Вариант 2

Вариант 3

Вариант 4

Вариант 5