A7TZrRHEJI
.pdfМИНОБРНАУКИ РОССИИ
–––––––––––––––––––––––––––––––––––––––––––––––––––––––
Санкт-Петербургский государственный электротехнический университет «ЛЭТИ» им. В. И. Ульянова (Ленина)
–––––––––––––––––––––––––––––––––––––––––––––
П. А. КОШЕЛЕВ С. В. ПАРАМОНОВ
КОМПЬЮТЕРНО-ИНФОРМАЦИОННЫЕ ТЕХНОЛОГИИ АНАЛИЗА И СИНТЕЗА ЭЛЕКТРОТЕХНИЧЕСКИХ КОМПЛЕКСОВ (ЭТК)
Часть 2
Нелинейные и нестационарные динамические системы
Учебное пособие
Санкт-Петербург Издательство СПбГЭТУ «ЛЭТИ»
2016
1
УДК 621.3.011.71: 621.3.011.73(07) ББК З29-5с11я7
К76
Кошелев П. А., Парамонов С. В.
К76 Компьютерно-информационные технологии анализа и синтеза электротехнических комплексов (ЭТК): в 2 ч. Ч. 2. Нелинейные и нестационарные динамические системы: учеб. пособие. СПб.: Изд-во СПбГЭТУ
«ЛЭТИ», 2016. 28 с.
ISBN 978-5-7629-2033-9 (ч. 2)
ISBN 978-5-7629-2029-2
В учебном пособии описаны методы анализа, синтеза и оптимизации параметров элементов электротехнических устройств и комплексов.
Методы и алгоритмы основаны на фундаментальных законах электротехники и реализованы в современных пакетах программ символьных и численных вычислений.
УДК 621.3.011.71: 621.3.011.73(07) ББК З29-5с11я7
Рецензенты: кафедра электротехники и прецизионных электромеханических систем ФГАОУ ВПО «Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики» (Университет ИТМО); канд. техн. наук В. С. Федорова (ПГУПС Императора Александра I).
Утверждено редакционно-издательским советом университета
в качестве учебного пособия
ISBN 978-5-7629-2033-9 (ч. 2) |
|
ISBN 978-5-7629-2029-2 |
© СПбГЭТУ «ЛЭТИ», 2016 |
2
Введение
Вреальной динамической системе, в том числе, в электротехническом устройстве, присутствуют нелинейные и нестационарные элементы. Они вносят зависимость параметров схемы от токов и (или) напряжений, от момента времени, заданного алгоритмом действия устройства или вследствие естественного хода процесса и приводят к непостоянству коэффициентов дифференциальных уравнений (ДУ).
Частный случай нелинейности (его можно назвать предельным) – нестационарность, вызванная наличием коммутаторов. Это типично для преобразователей электроэнергии, которые присутствуют практически в каждом электротехническом комплексе.
Вследствие коммутаций изменяться могут не только коэффициенты, но количество и содержание уравнений, которые описывают нелинейную динамическую систему (НДС).
При исследовании и моделировании НДС иногда используют искусственный прием, сущность которого заключается в заранее заложенной избыточности математического описания задачи.
Взависимости от применяемого аппарата исследования, система ДУ или топология математической модели (ММ) предусматривают все ожидаемые состояния НДС.
Возрастание порядка матриц коэффициентов ДУ неизбежно, однако, при использовании современных аппаратных средств и численно-аналитических программ, таких, как MAPLE, MATHEMATICA, MATLAB в большинстве практических задачвремярасчета и требуемый ресурспамяти вполнеприемлемы.
3
1. ОБЩИЕ СВОЙСТВА ЗАДАЧ АНАЛИЗА НДС
Один из распространенных методов анализа нелинейных систем основан на принципе припасовывания.
Сущность этого принципа в том, что система считается линейной на достаточно малом временнόм интервале анализа.
Если в процессе функционирования устройства количественно или качественно изменяется его структура а, соответственно, при моделировании реальной ситуации изменяется какой-либо коэффициент или группа коэффициентов ДУ, вектор времени анализа обнуляется и дальнейший расчет ведется с начальными условиями, сохранившимися в предыстории задачи, но с новой совокупностью коэффициентов.
Начальные условия формируются в соответствии с законами коммутации только для реактивных элементов (индуктивностей, емкостей).
Результат расчета на очередном шаге присоединяется к предшествующим шагам.
Физически начальные условия отражают процесс запасания электромагнитной энергии в реактивных элементах за счет токов индуктивностей и напряжений емкостей.
После изменения какого-либо параметра цепи предшествующие значения токов индуктивностей и напряжений на емкостях распространяются на следующий интервал расчета, т. е. в момент перехода на следующий шаг считаются постоянными.
Следует различать 3 случая.
●Моменты коммутации заранее известны. Примеры: – автономные инверторы с ШИМ управлением, скачкообразное изменение сопротивления нагрузки при выполнении определенных условий и т. п. Это наиболее «простой» случай.
●Моменты коммутации заранее неизвестны и определяются мгновенными значениями или знаком сигнала текущего процесса, в частности, полярностью напряжения на участке анод-катод неуправляемого вентиля.
●Система содержит «гладкие» нелинейности. В любой момент коэффициенты ДУ «ожидают» переопределения. Это наиболее сложный случай с точки зрения оценки ошибки округления и адекватности моделирования процесса рассеяния запасенной энергии.
Приемлемым средством контроля может стать расчет в альтернативных
вычислительных средах, моделирование в примененном пакете программ с вариацией данных.
4
Применительно к методу узловых потенциалов (МУП) начальные условия отображаются дополнительными источниками тока, а для метода контурных токов (МКТ) – источниками ЭДС или (в операторной форме) их изображениями по формуле прямого преобразования Лапласа.
2. УЧЕТ НАЧАЛЬНЫХ УСЛОВИЙ (МЕТОД УЗЛОВЫХ ПОТЕНЦИАЛОВ)
Если в цепи некоторые двухполюсники являются реактивными, параллельно им включаются источники тока.
Ток емкости
ic C |
dUc |
(2.1) |
|
dt |
|||
|
|
Изображение тока емкости, распространенного на следующий шаг:
i |
(p) |
pCUc0 |
CU |
c0 |
(2.2) |
|
|||||
c |
|
p |
|
||
|
|
|
|
Ток индуктивности и операторное изображение тока дополнительного источника:
i |
i |
const; i |
L |
( p) iL0 . |
(2.3) |
L |
L0 |
|
p |
|
|
|
|
|
|
|
3. УЧЕТ НАЧАЛЬНЫХ УСЛОВИЙ (МЕТОД КОНТУРНЫХ ТОКОВ)
Для учета начальных условий в уравнения линейной системы добавляется слагаемое, представляющее собой сумму изображений начальных условий.
Матрица сопротивлений [Z] при этом не изменяется, а в столбцы вектора задающих воздействий [J] включаются изображения токов емкостей и индуктивностей, причем номер строки вектора [J] и номер узла, к которому подключен реактивный элемент, совпадают.
Начальные условия изображаются источниками ЭДС или (в операторной форме) их изображениями.
Расчетная схема цепи питания плазмотрона, в которой учтены начальные условия, приведена на рис. 3.1.
Для выбранных направлений контурных токов (в этом случае «по часовой стрелке») напряжения источников ЭДС, отражающих начальные условия, записываются следующим образом.
5
|
R1 L1 |
L1i10 |
L3 |
L3ii20 |
|
|
|
|
|
– |
+ |
i2 |
– |
+ |
+ |
|
|
|
|
|
|
С |
|
|
||
|
|
L2 |
|
|
C |
|
|
R3 |
|
ii1 |
|
U |
//p |
|
i |
||
U |
|
–– |
c0 |
– |
i3 |
R3 |
||
|
|
|
|
|
|
|||
|
i |
– i ) |
|
R2 |
|
|
||
|
L2(i10-i20) |
+ |
|
|
|
|
||
|
|
|
|
|
|
|
|
Рис. 3.1
Напряжения на индуктивностях в естественном виде и в операторной форме (3.1):
UL (tm) L di |
; UL ( p) pLi( p) |
Li0. |
(3.1) |
dt |
p |
|
|
В конце исследуемого интервала расчета tm соответствующий ток уменьшается, а напряжение считается постоянным вплоть до начала следующего интервала расчета, поэтому производной присваивается знак (–).
Здесь i0 – значение тока индуктивности в конце интервала tm. Напряжение на конденсаторах и его изображение распространяются на
следующий интервал со своим знаком (3.2):
|
1 t |
|
|
|
|
Uc0 |
|
||
|
|
|
|
||||||
|
|
|
|
|
|
|
|||
Uc (t) C |
ic (t)dt ; |
p . |
(3.2) |
||||||
Uc (p) L Uc (t) |
|||||||||
|
|
0 |
|
|
|
|
|
|
Uc0 – значение напряжения на конденсаторе в конце интервала расчета
на текущем шаге.
Матрица [Z] не зависит от наличия начальных условий. Изменяется только вектор задающих источников [Е] (рис. 3.2)
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
2 |
|
|
|
3 |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
R1 + pL1 + pL2 |
|
|
–pL2 |
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
1 |
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
0 |
|
|||||||||||||||||||||
[Z] = |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
–pL2 |
|
pL2 + pL3 + 1/(pC) + R2 |
|
|||||||||||||||||||||
|
2 |
|
|
|
|
|
–(1/p/C + R2) |
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
–(1/p/C + R2) |
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1/p/C + R2 + R3 |
|||||||||
|
|
3 |
0 |
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
1 |
|
|
|
|
(U/p) |
L1i10 |
L2(i10 |
i20) |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
[E] = |
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
L3i20 + Uc0/p + L2(i10 i20) |
|||||||||||||||
|
|
|
2 |
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Uc0/p |
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 3.2
6
Это справедливо и для МУП (для матрицы [Y] и вектора [I01]).
На рис. 3.3, а – схема замещения для МКТ, на рис. 3.3, б – для МУП. Отметим, что в расчетной схеме для МКТ последовательно с источни-
ками ЭДС и конденсаторами должны быть включены сопротивления, а для МУП параллельно подключены проводимости. Если в исходной схеме этих элементов нет, то они вводятся искусственно (g4 на рис. 3.3, б).
R1 L
U
G
i1
L· i10
a
C
2R
uc0/p
|
|
|
|
|
|
|
I0 |
|
|
||
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
3R |
|
|
|
|
|
|
|
|
|
|
|
|
|
G |
||
|
|
|
|
|
|
|
|
|
|||
|
i2 |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 3.3 |
|
|
U1 L
1g
0
i10/pU2 C |
C · uc0 U3 |
|
2g |
g4 |
3g |
б |
|
|
Матрица операторных сопротивлений Z и вектор воздействий E1 для нелинейной стационарной схемы (МКТ), по рис. 3.3, а – на рис. 3.4.
|
R1 pL R2 |
|
R2 |
U |
|
|
UCO |
||||
Z : |
|
|
|
1 |
|
L110, |
|
||||
|
R2 |
R2 R3 |
|
E : |
p |
p |
|
||||
|
|
pC |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 3.4
Матрица операторных проводимостей Y схемы (МУП) по рис. 3.3, б – на рис. 3.5.
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
1 |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
Lp |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
R1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Lp |
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
pC |
|
|
|
|
||||
Y : |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
R4 |
|
|
|
|
||||||||
|
|
1 |
|
|
1 |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
1 |
pC |
|
1 |
|
|
|||||||||
|
|
|
Lp |
|
|
pC |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Lp |
R2 |
R4 |
|||||||||||||||||
|
|
Lp |
|
R2 |
R4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
pC |
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
R4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
pC |
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
R3 |
|
R4 |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 3.5
7
Вектор воздействий I01 для нелинейной (нестационарной) схемы (МУП) на рис. 3.6.
|
|
U |
|
|
i10 |
|
|
|
i10 |
C uC0 |
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
p |
|
|
|
p |
|
|
|
|
C uC0 |
|
|
|
||||||||||
R1p |
|
|
|
|
|
|
|
|
||||||||||||||||||
I 01: |
|
|
|
|
|
|
|
|
|
, |
|
|
|
|
|
|
|
|
, |
|
|
|
|
|
|
|
|
1 |
|
|
|
|
1 |
|
|
1 |
1 |
|
|
1 |
|
|
1 |
|
1 |
|
|||||||
|
|
|
|
|
|
|
|
|
pC |
|
|
|
pC |
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
R1 |
|
|
Lp |
|
|
|
Lp |
R1 |
R4 |
|
|
|
R3 |
R4 |
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 3.6
Результаты моделирования процесса при переключении сопротивления нагрузки R3 по рис. 3.3, а – на рис. 3.7. В начале R3 = 10 Ом, затем при t = 0,35 ms оно становится 100 Ом, в момент 0,6 ms возвращается к прежнему значению 10 Ом.
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Напряжение на конденсаторе, В |
|
|||||||||||||||||||||
|
4 |
|
Нелинейнаясистема. Напряжение нагрузки, В |
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
10 |
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
3,5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
9 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
8 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
2,5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
7 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
1,5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
0,5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
0,5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
–1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
0 |
|
20 |
40 |
60 |
80 |
100 |
120 |
140 |
160 |
|
180 |
200 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
20 |
40 |
60 |
80 |
100 |
120 |
140 |
160 |
|
180 |
200 |
|
|
|||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
№ шага |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
№ шага |
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 3.7
Структура модели схемы по рис. 3.3 в системе MATLAB, в пакете SIMULINK POWER SYSTEMS (SPS) показана на рис. 3.8, результаты моделирования на рис. 3.9.
8
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
R1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
L |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
Voltage Measurement1 |
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
C 1e-5 |
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
1e-3 |
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Scope1 |
||||||||
|
|
|
|
|
|
|||||||||||||||||||
1 |
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
10 V |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
Series RLC Branch 1 |
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
Series RLC Branch |
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
Voltage Measurement |
|||||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
Controlled Voltage Source |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
10 |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
Ideal Switch2 |
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Step |
|
|
|
|
|
|
|
|
|
|
|
|
|
R3 |
|
|
|
|
|
|
|
|
|
|
|
Timer1 |
|
|
|
||
|
|
|
|
|
|
||
|
|
|
|
|
|
90 |
|
|
|
|
|
|
|
|
|
|
|
|
R2 10 |
|
|
||
|
|
|
|
|
|
Рис. 3.8
Напряжение нагрузки, В Напряжение на конденсаторе, В
10
9
8
7
6
5
4
3
2
1
0
3,5
3
2,5
2
1,5
1
0,5
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0,1 |
|
0,2 |
|
0,3 |
|
0,4 |
|
0,5 |
|
0,6 |
|
0,7 |
|
0,8 |
|
|
0,9 |
|
1 |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
Time |
offset: 0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Время |
, с · 10–3 |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 3.9
4.МОДЕЛИРОВАНИЕ ВЕНТИЛЬНЫХ ПРЕОБРАЗОВАТЕЛЕЙ
СНЕОПРЕДЕЛЕННЫМИ МОМЕНТАМИ КОММУТАЦИИ
4.1.Анализ 3-фазного управляемого выпрямителя (УВ, мостовая схема)
Расчетная схема преобразователя показана на рис. 4.1.
9
Управляемые вентили (тиристоры) представлены R – моделями. Сопротивление открытого вентиля rl ≈ 0, сопротивление запертого вентиля достаточно велико: rm ≈ 105 Ом, примерно соответствует среднему реальному значению для низкочастотных приборов.
В соответствии с алгоритмом R = rl (вентиль открыт), если его анодное напряжение Ua положительно и пройден интервал задержки включения, соответствующий в угловых единицах alpha в градусах.
Запирание вентиля (R = rm) происходит при изменении полярности Ua. Начальные условия переопределяются (см. метод припасовывания) на каждом шаге вычислений. Это позволяет имитировать динамические процессы практически при всех реальных ситуациях, в т. ч. изменении параметров схемы, аварийных режимах и др., что является преимуществом по сравнению
с использованием имитационных блоков SPS.
|
|
|
|
|
|
|
|
|
|
|
|
|
Ld |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Ld |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
R |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
R |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
r1c |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
r1c |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
r2c |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
r1b |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
r1b |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
r2c |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
r2b |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
r2b |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
i2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i5 |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
r1a |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
r1a |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
r2a |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
r2a |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
ea |
|
|
eb |
|
|
|
ec |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ea |
|
|
eb |
|
|
|
|
ec |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
L |
|
|
|
|
|
|
|
|
|
L |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
L |
|
|
L |
|
|
L |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
L |
|
|
|
|
|
|||||||||||||||||
|
i1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
i3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i6 |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Риc. 4.1
Поскольку в 3-фазной схеме УВ в рассматриваемой модели 1 вентиль анодной группы может работать поочередно с 2 вентилями катодной группы и наоборот, например, открытыми могут быть пары r1a, r2b или r1a, r2c, в каждой фазе сети предусмотрено протекание 2 токов: соответственно i1, i4 (рис. 4.1).
10