
ЭиУСУ / ЛИТЕРАТУРА_ЭиУСУ / Глазырин - МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ЭЛЕКТРОМЕХАНИЧЕСКИХ СИСТЕМ
.pdf
41
Рис. 17. Структурная схема замкнутой ЭМС с П-регулятором
Уравнение электрического равновесия ДПТ НВ –
Ud (t) = LДВ didt(t) + RДВ i(t) +c ω(t).
Уравнение механического равновесия ДПТ НВ –
с i(t) −MС = JДВ dωdt(t) .
Разрешая каждое уравнение относительно производных, запишем СДУ в нормальной форме Коши:
dUd (t) |
= |
|
1 |
|
|
kПР kР (UЗАД −kОС kТГ ω(t))−Ud (t) ; |
||||||
|
|
|
|
|
|
|||||||
dt |
|
T |
|
|||||||||
|
|
|
|
1 |
ПР |
|
|
|
||||
|
|
|
|
|
|
|
|
|
||||
di(t) |
|
|
|
|
|
|
||||||
|
dt |
= |
|
|
|
Ud (t) −RДВ i(t) −c ω(t) |
; |
|||||
L |
|
|
||||||||||
|
|
|
|
ДВ |
|
|
|
|
|
|
|
|
dω(t) = |
|
1 |
с i(t) −M |
. |
|
|||||||
JДВ |
|
|||||||||||
|
dt |
|
|
|
|
|
С |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
СДУ в матричном виде –
|
|
|
|
|
|
|
− |
|
1 |
|
|
|
|
0 |
|
|
− |
k |
ПР |
k |
Р |
k |
ОС |
k |
ТГ |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
TПР |
|
|
|
|
|
|
|
TПР |
|
|
|
|||||||||||||||
d |
U |
|
(t) |
|
|
|
|
|
RДВ |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
d |
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
c |
|
|
|
|
|
|||||||
|
|
i(t) |
|
= |
|
|
|
|
|
|
− |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
× |
||||
dt |
|
|
L |
|
|
L |
|
|
|
|
|
|
L |
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
ω(t) |
|
|
|
|
ДВ |
|
|
|
|
ДВ |
|
|
|
|
|
|
|
|
ДВ |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
c |
|
|
|
|
|
|
|
0 |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
JДВ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
kПР kР UЗАД |
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
U |
|
|
(t) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
d |
|
|
|
|
T |
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ПР |
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
× i(t) |
+ |
|
|
|
|
|
0 |
|
|
|
|
1(t). |
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
MС |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ω(t) |
|
|
|
|
|
− |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
JДВ |
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Здесь
Ud (t)
x(t) = i(t) – вектор переменных состояния;
ω(t)
|
− |
1 |
|
||
|
|
|
|||
TПР |
|||||
|
|
|
|||
|
|
|
1 |
|
|
A = |
|
|
|
|
|
|
L |
||||
|
|
||||
|
|
|
ДВ |
||
|
|
|
0 |
|
|
|
|
|
|
|
|
|
0 |
|
− |
k |
ПР |
k |
Р |
k |
ОС |
k |
ТГ |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
TПР |
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
RДВ |
|
|
|
|
|
c |
|
|
|
|
– матрица коэффициентов перед пе- |
||
− |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
L |
|
|
|
|
L |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
ДВ |
|
|
|
|
|
ДВ |
|
|
|
|
|
||
|
|
c |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
JДВ |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
ременными состояния;
42
kПР kР UЗАД |
|
|||||
|
|
|
|
|
|
|
TПР |
||||||
|
|
|||||
B = |
0 |
|
– вектор свободных членов СДУ. |
|||
|
− |
MС |
|
|
||
|
|
|
||||
|
||||||
|
|
J |
ДВ |
|
||
|
|
|
|
Из полученной математической модели видно, что помимо двух переменных состояния двигателя (тока и скорости) ЭМС имеет третью переменную состояния – выходное напряжение ШИП, обусловленную инерционностью преобразователя.
2.5.2. Модельзамкнутойэлектромеханическойсистемы сПИ-регулятором, двигателемпостоянноготока независимоговозбуждения исиловымпреобразователем, представленным пропорциональнымзвеном
Структурная схема ЭМС с ПИ-регулятором, ДПТ НВ и ШИП, представленным пропорциональным звеном, показана на рис. 18.
Втакомслучаеуравнениеэлектрическогоравновесиябудетиметьвид
kР (UЗАД −kОС kТГ ω(t))+UУ.И(t) kПР =
= LДВ didt(t) +RДВ i(t) +c ω(t),
где UУ. И(t) – напряжение управления интегральной частью регулятора. Уравнение механического равновесия останется без изменений:
сi(t) −MС = JДВ dωdt(t) .
Кэтим двум уравнениям прибавляется уравнение состояния
ПИ-регулятора |
|
|
|
|
|
|
|
|
dUУ.И(t) |
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
T |
|
=U |
|
|
−k |
|
k |
|
ω(t), |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
Р |
|
|
|
|
|
|
|
dt |
|
ЗАД |
|
|
ОС |
|
ТГ |
|
|||
где ТР – постоянная времени регулятора. |
|
|
|
|
||||||||||||||||||||||
СДУ в нормальной форме Коши – |
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
kР (UЗАД −kОС kТГ ω(t))+UУ.И |
||||||||||||||
di(t) |
|
|
|
1 |
|
|
|
|||||||||||||||||||
|
|
= |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
dt |
L |
×kПР −RДВ i(t) +c ω(t) |
||||||||||||||||||||||||
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
ДВ |
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
dω(t) |
= |
|
|
|
|
|
с i(t) −MС |
; |
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
dt |
|
|
JДВ |
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
dUУ.И(t) |
= |
|
|
|
1 |
|
|
U |
−k k |
|
ω(t) . |
|
|||||||||||||
|
T |
|
|
|
|
|||||||||||||||||||||
|
dt |
|
|
|
|
|
|
|
ЗАД |
ОС |
|
|
ТГ |
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
Р |
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(t) × ;
43

44
Рис. 18. Структурная схема замкнутой ЭМС с ПИ-регулятором

СДУ в матричной форме –
d dt
|
|
|
|
− |
RДВ.Г |
||
|
|
|
|
||||
|
i(t) |
|
|
|
|
L |
|
|
|
ДВ |
|||||
|
ω(t) |
= |
|
|
c |
|
|
|
|
JЭ |
|||||
|
|
|
|
|
|
||
UУ.И(t) |
|
|
|
|
|
||
|
|
|
|
0 |
|
||
|
|
|
|
|
|
|
|
−kОС kТГ kР kПР +c LДВ
0
−kОС kТГ
TP
kР kПР UЗАД |
|
|
|||||
|
|
|
|
|
|
|
|
|
LДВ |
|
|||||
|
|
|
|
||||
|
|
|
MС |
|
|
|
|
+ |
− |
|
|
1(t). |
|||
JЭ |
|||||||
|
|
|
|
|
|||
|
UЗАД |
|
|
||||
|
|
|
|||||
|
|
|
|
|
|
|
|
|
TP |
|
|||||
|
|
|
|
k |
ПР |
|
|
|
|
|
|
||
L |
i(t) |
|||
|
||||
|
ДВ |
|||
0 |
|
|
||
|
ω(t) |
|||
|
|
|
|
|
|
|
|
UУ.И |
|
0 |
|
|
||
|
|
|
|
+ (t)
Здесь
|
|
|
|
i(t) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
ω(t) |
|
– вектор переменных состояния; |
||||||||||||||||||
x(t) = |
|
|
|||||||||||||||||||||
|
UУ.И(t) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
− |
|
RДВ.Г |
− |
k |
|
k |
k |
Р |
k |
ПР |
+c |
|
k |
ПР |
|
|
||||||
|
|
|
|
|
|
|
|
ОС |
|
ТГ |
|
|
|
|
|
|
|
||||||
|
|
|
|
LДВ |
|
|
|
|
|
|
LДВ |
|
|
|
|
LДВ |
|
|
|||||
|
|
|
|
c |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
A = |
|
|
|
|
|
|
|
|
0 |
|
|
|
|
0 |
|
– матрица коэффициентов пе- |
|||||||
|
|
|
JЭ |
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
kОС kТГ |
|
|
|
|
|
|
|
|||
|
|
|
0 |
|
|
|
|
|
|
− |
|
|
0 |
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
TP |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
ред переменными состояния; |
|
|
|
|
|
|
|||||||||||||||||
kР kПР UЗАД |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
LДВ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
MС |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
B = |
|
|
|
− |
|
|
|
– вектор свободных членов СДУ. |
|||||||||||||||
|
|
|
JЭ |
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
UЗАД |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
TP |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Третьей переменной состояния ЭМС в данном случае выступает напряжение управления интегральной частью ПИ-регулятора.
45

3. АНАЛИЗ ВЫХОДНЫХ ПРОЦЕССОВ ЭЛЕКТРОМЕХАНИЧЕСКИХ СИСТЕМ С ПРИМЕНЕНИЕМ КЛАССИЧЕСКИХ СПОСОБОВ РЕШЕНИЯ ЗАДАЧИ КОШИ
3.1. Решение дифференциальных уравнений классическим способом на примере RL-, RC- иRLC-цепей
3.1.1. Нахождениерешениядифференциальногоуравнения напримереRL-цепи
Схема коммутации RL-цепи на источник постоянного напряжения представлена на рис. 19.
Рис. 19. Схема коммутации RL -цепи на источник постоянного напряжения
Процессы, протекающие в цепи при замыкании ключа, описываются дифференциальным уравнением 1-го порядка, составленным по второму закону Кирхгофа:
E 1(t) =i(t) R + L didt(t) .
Характеристическое уравнение получается из однородного ДУ для переменной x(t) путем его алгебраизации, которая заключается в заме-
не оператора дифференцирования на переменную λ. При этом степени новых переменных λ равны порядку соответствующих производных. После алгебраизации переменная x(t) выносится за скобки, а полином,
находящийся в скобках, после приравнивания его к нулю будет представлять собой характеристическое уравнение.
Однородное дифференциальное уравнение – i(t) R + L didt(t) = 0.
46

Определим корень характеристического уравнения
R +L λ =0;
λ = − RL .
Общее решение i0 (t) однородного уравнения –
i0 (t) = N eλ t = N e−RL t .
Для нахождения частного решения неоднородного ДУ подставим в исходное уравнение значение t = ∞. Тогда оно будет выглядеть сле-
дующим образом:
E =iЧ R; iЧ = ER .
Общее решение неоднородного уравнения – i(t) =iЧ +i0 (t) = ER + N e−RL t .
Решим задачу Коши с нулевыми начальными условиями. Для этого подставим в общее решение неоднородного ДУ значение t = 0. Получим
i(0) = ER + N = 0,
откуда постоянная интегрирования
N = − ER .
λ
Рис. 20. Изменение тока в RL-цепи при коммутации на источник постоянного напряжения
47

В итоге временная зависимость тока в RL -цепи при коммутации ее на источник постоянного напряжения
i(t) = ER − ER e−RL t .
Графическое решение ДУ представлено на рис. 20.
На графике обозначена также постоянная времени цепи Т, обратная модулю корня характеристического уравнения T =1/ λ = RL .
3.1.2. Нахождениерешениядифференциальногоуравнения напримереRC-цепи
Схема коммутации RC-цепи на источник постоянного напряжения представлена на рис. 21.
iC(t) R
E·1(t) |
C |
UC(t)
Рис. 21. Схема коммутации RC-цепи на источник постоянного напряжения Е
Дифференциальное уравнение, составленное по второму закону Кирхгофа, описывает процессы в цепи после замыкания ключа:
E 1(t) =iC (t) R +UC (t).
Учитывая, что i (t) =C |
dUC (t) |
, это уравнение можно записать |
|||||
|
|
||||||
C |
|
|
dt |
|
|
|
|
в следующем виде: |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
E 1(t) = RC |
dUC (t) |
+UC (t). |
||||
|
|
||||||
|
|
|
|
|
|
dt |
|
Найдем общее решение однородного уравнения: |
|||||||
|
RC |
dUC (t) |
+UC (t) =0. |
||||
|
|
||||||
|
|
|
dt |
|
|
|
48

Определим корень характеристического уравнения
RC λ +1=0;
λ = − RC1 .
Общеерешение UC0 (t) однородногодифференциальногоуравнения–
UC0 (t) = N eλt = N e−RCt .
Для нахождения частного решения неоднородного ДУ подставим в исходное уравнение значение t = ∞. Тогда оно будет выглядеть следующим образом:
UC.ЧАСТН = E.
Общее решение неоднородного уравнения –
− t
UC (t) =UC.ЧАСТН +UC0 (t) = E + N e RC .
Решим задачу Коши с нулевыми начальными условиями. Для этого подставим в общее решение неоднородного ДУ значение t = 0. Получим
UC (0) = E + N =0,
откуда постоянная интегрирования
N= −E.
Витоге временная зависимость тока в RC-цепи при коммутации ее на источник постоянного напряжения
− t
UC (t) = E −E e RC .
Графическая зависимость напряжения на конденсаторе от времени представлена на рис. 22.
λ
Рис. 22. Изменение напряжения на конденсаторе при коммутации на источник постоянного напряжения
49

На графике обозначена также постоянная времени цепи Т, обратная модулю корня характеристического уравнения T =1/ λ = R C.
3.1.3. Нахождениерешениясистемыдифференциальных уравненийнапримереRLC-цепи
Коммутация RLC-цепи на источник постоянного напряжения представлена на рис. 23.
Рис. 23. Схема коммутации RLС-цепи на источник постоянного напряжения
Составим СДУ по первому и второму законам Кирхгофа:
|
E 1(t) |
=i(t) R + L di(t) +UC (t); |
|
|||||||
|
|
|
|
|
|
|
|
dt |
|
|
|
|
|
|
dUC |
(t) |
|
|
|||
|
|
|
|
|
|
|||||
|
i(t) =C |
|
|
|
|
. |
|
|||
|
|
|
|
|
|
|||||
|
|
|
|
|
|
dt |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Представим СДУ в нормальной форме Коши: |
|
|||||||||
|
|
|
1 |
|
|
|
|
|
|
|
di(t) = |
|
E 1(t) −i(t) R −UC (t) |
; |
|||||||
|
|
|||||||||
|
dt |
L |
|
|
|
|
|
|
||
|
dUC (t) |
= |
|
1 |
i(t). |
|
|
|||
dt |
|
C |
|
|
||||||
|
|
|
|
|
|
|
|
Запишем СДУ в матричной форме:
|
|
|
|
|
|
|
|
R |
|
d |
|
i(t) |
|
= |
− L |
||||
|
|
|
|
||||||
dt U |
|
(t) |
|
1 |
|
||||
|
|
|
C |
|
|
|
|
|
|
|
|
|
|
C |
|||||
|
|
|
|
|
|
|
|
− |
1 |
|
i(t) |
|
|
E |
|
|
|
|
|||||||
|
L |
+ |
L |
1(t). |
||||
|
|
|
U |
(t) |
|
|
|
|
0 |
|
C |
0 |
|
||||
|
|
|
||||||
|
|
|
|
|
|
|
|
|
Здесь
|
− R |
− |
1 |
|
|
||
|
|
|
|||||
A = |
|
L |
|
L |
– матрица коэффициентов перед переменными состоя- |
||
|
1 |
|
0 |
|
|||
|
|
|
|
|
|
||
|
C |
|
|||||
|
|
|
|
|
|
ния RLC-цепи;
50