Учебное пособие 1920
.pdf
|
T |
|
2 |
]dt min, |
|
|
|
|
|||
[(J Mc ) |
|
||||
|
0 |
|
|
|
(2.7) |
|
|
|
|
|
|
T |
|
(0) 0; |
(T) 0. |
||
(t)dt ; |
|||||
0 |
|
|
|
|
|
Задача (2.7) представляет собою изопериметрическую задачу с вспомогательной функцией F и закреплёнными концами. В ней необходимо найти две неизвестные функции ω(t) и(t), причём в данном частном случае в силу изопериметрического условия (2.3) справедливо (t)= = const. Для определения ω(t) и имеем два уравнения – уравнение Эйлера
|
|
|
|
|
|
|
|
|
|
|
|
|
|
F |
|
|
|
d F |
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dt |
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
T |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
и уравнение связи (t)dt . |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Подставляя (2.6) в уравнение Эйлера, получаем |
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
d |
|
(2J |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
(2.8) |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
dt |
|
|
2JMc ) 2J |
0. |
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Решение дифференциального уравнения (2.8) имеет вид |
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t2 c t c . |
|
|
|
|
|
(2.9) |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4J2 |
|
|
|
1 |
0 |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
Подставляя (2.9) в (2.3), получим |
|
|
|
|
|
|
|
||||||||||||||||||||||||
T |
|
|
|
2 |
|
|
|
|
t3 |
|
|
c t |
2 |
|
|
|
|
T |
T3 |
|
cT2 |
|
|
(2.10) |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
t |
|
c1t c0 |
dt |
|
|
|
|
|
|
|
|
|
|
|
|
c0t |
|
|
|
|
|
|
c0T |
. |
|
||||
4J |
2 |
|
|
12J |
2 |
2 |
|
|
12J |
2 |
2 |
|
|||||||||||||||||||||
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Положив ω(0)=ω1 и ω(Т)=ω2, из (2.9) получим с0=ω1. То-
гда совместное решение (2.9) для t=T и (2.10) даёт:
12J2 ( 1 T23)T 2 ; c1 6 2(2T21 2)T .
Для условий (2.7) ω(0)=0 и ω(Т)=0 имеем:
|
24J 2 |
|
; c |
6 |
, |
T3 |
|
T2 |
|||
|
|
1 |
|
130
и с учётом (2.9) окончательно: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
(t) |
6 |
t2 |
|
6 |
t , |
|
|
|
|
|
|
(2.11) |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
и |
|
|
|
|
|
|
|
|
|
|
T3 |
|
|
T2 |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
M |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
i(t) |
J d |
|
|
c |
1 |
|
12J |
t |
|
|
6J |
Mc |
|
(2.12) |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
||||||||||
cm dt |
|
|
cm |
|
|
T3 |
T2 |
|||||||||||||||||||||
|
|
|
|
|
|
|
cm |
|
|
|
|
|
|
|
|
|||||||||||||
Для определения вида найденного экстремума проверяем |
||||||||||||||||||||||||||||
|
|
|
|
|
2 |
F |
|
|
|
2J |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||
условие Лежандра: |
|
|
|
|
2JMc |
2J 2 |
0 |
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таким образом, найденное решение (2.11) и (2.12) обеспечивает функционалу Q минимум.
|
i |
10s+70 |
1 |
Signal 1 |
2.5 |
s |
0.05s+0.06 |
Signal 2
|
|
omega |
0.01 |
1/s |
1/s |
2.5 |
|
fi |
|
|
Рис. 2.1. Имитационная модель системы оптимального управления поворотом платформы
Для иллюстрации полученного решения составим имита-
ционную модель объекта (рис. 2.1) с параметрами L=0,05 Гн;
R=0,06 Ом; J=100 кг м2; U=400 В; M=300 Н м; се=сm=2,5; iном=200 А; ωном=1500 об/мин. Функция момента сопротивле-
ния формируется |
в виде момента сухого трения |
Mc M sign( ). |
Двигатель снабжён пропорционально- |
интегрирующим регулятором с передаточной функцией W(s):
W(s) 10s 70 , в совокупности с которым образует следя- s
щую систему регулирования тока двигателя по желаемому закону, в частности, вида (2.12). Графики изменения тока i(t) при оптимальном управлении (2.12), частоты вращения ω(t) и угла
131
поворота (t) при заданном значении последнего 100 радиан за время Т=5 с, представлены на рис. 2.2а,б.
1000 |
Оптимальное |
|
управление |
||
i, A |
||
|
||
500 |
i(t) |
|
|
||
0 |
|
|
-500 |
|
|
-850 |
|
|
-10000 1 2 3 4 5 6 t, с |
Рис. 2.2а. Изменение тока при оптимальном управлении (2.12)
Потери энергии при полученном оптимальном управлении (2.12) для Т=5 с, М=300 Н м и =100 рад составляют
T |
(t)dt |
T R |
12J |
|
6J |
2 |
||||
P R i2 |
|
c2 |
|
|
T3 |
t |
M dt |
|||
|
|
|
|
|
|
T2 |
|
|||
0 |
|
|
0 |
m |
|
|
|
|
|
|
R12J2 2 |
T4M2 |
J 100 |
|
96,48кДж. |
||||||
|
|
2 |
|
3 |
|
|
100 |
|
|
|
|
|
cmT |
|
|
|
T 5 |
|
|
|
|
|
|
|
|
|
|
|
cm 2,5 |
|
|
|
|
|
|
|
|
|
|
R 0,06 |
|
|
|
, рад |
|
|
|
|
|
|
|
|
|
|
1 0 0 |
|
|
|
|
|
|
|
|
|
|
ω, с-1 |
|
|
|
|
(t) |
|
|
|
||
5 0 |
|
|
|
|
|
|
|
|
|
|
3 0 |
|
|
|
|
|
|
|
|
ω(t) |
|
0 |
|
|
|
|
|
|
|
|
|
6 t, c |
0 |
1 |
|
2 |
|
3 |
4 |
|
5 |
Рис. 2.2б. Изменение частоты вращения ω(t) и угла поворота (t) при оптимальном управлении
132
Для сравнения рассмотрим изменения тока при оптимальном управлении (2.12) с тривиальным решением задачи, при котором половину времени поворота платформы её частота вращения линейно возрастает до значения ωmax, а вторую половину – снижается до нуля.
В этом случае для того чтобы за время T повернуть платформу на угол , потребуется выполнение равенства
|
|
|
T |
max |
dt |
maxT |
. |
||||||||||||
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
2 |
|
|
2 |
|
|
|
|
||
|
|
|
|
2 |
|
0 |
|
|
|
|
|
|
|
|
|
|
|||
Отсюда |
max |
|
|
и далее |
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
T |
|
|
|
|
|
|
|
|
|
T |
|
|
|
||
|
|
|
4 |
t, |
при 0 t |
|
; |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
2 |
2 |
|
|||||||||||||
(t) |
T |
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
4 |
|
|
|
|
|
|
T |
|
|
|
||||||
|
|
|
|
|
|
|
(t T ), |
при |
|
|
|
t T . |
|||||||
|
|
|
|
|
T |
2 |
|
2 |
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
С учётом (2.4) получаем
|
|
i(t) |
J |
|
d |
|
Mc |
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
cm |
|
|
|
|
dt |
|
|
cm |
|
||||||||
|
4J |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
T |
|
||||||
|
|
|
Mc |
|
|
|
, |
|
|
при0 t |
|
; |
(2.13) |
|||||||||
|
2 |
|
|
|
|
|
||||||||||||||||
|
T |
|
|
cm |
2 |
|
|
|||||||||||||||
|
|
|
|
|
|
|||||||||||||||||
|
|
4J |
|
|
|
|
1 |
|
|
|
|
T |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
M |
|
|
|
|
|
|
|
|
, при |
|
t T. |
|
||||||
|
|
|
c |
|
|
|
|
|
||||||||||||||
|
2 |
|
|
|
||||||||||||||||||
|
|
T |
|
|
cm |
2 |
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
Графики изменения тока i(t) при неоптимальном управлении (2.13), частоты вращения ω(t) и угла поворота (t) при тех же заданных значениях =100 радиан и Т=5 с представлены на рис. 2.3а,б.
133
1000 |
Неоптимальное |
|
|
750 |
управление |
500 i(t)
i, A
0
-500 0 1 2 3 4 5 6t, c
Рис. 2.3а. Изменение тока при неоптимальном управлении (2.13)
, рад |
|
|
|
|
|
|
|
||
1 0 |
|
0 |
|
|
(t) |
|
|
|
|
ω, с-1 |
|
|
|
|
|
|
|
||
5 0 |
|
|
|
ω(t) |
|
|
|
||
4 |
0 |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
0 |
|
0 |
1 |
2 |
3 |
4 |
5 |
6 |
t, c |
|
|
Рис. 2.3б. Изменение частоты вращения ω(t)
и угла поворота (t) при неоптимальном управлении
Потери энергии при неоптимальном управлении (2.13)
|
|
|
|
|
|
|
T |
|
|
|
2 |
|
|
|
|
|
|
|
|
2 R |
|
|
|
T |
|
2 |
|
|
2 |
4J |
|
|
R |
|
|
T |
|
4J |
|
||||||
P |
|
i |
|
(t)Rdt |
|
|
|
M |
|
|
dt |
|
|
|
|
M |
|
dt |
||||
|
|
c2 |
|
c2 |
||||||||||||||||||
|
|
|
|
|
T2 |
|
|
|
T2 |
|
|
|||||||||||
|
0 |
|
|
|
|
0 |
|
|
|
|
|
m |
|
|
T |
|
|
|
|
m |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
16J2 2 T4M 2 |
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
R |
|
|
|
|
|
|
|
|
J 100 |
127,2 кДж. |
|
|
||||||
|
|
|
|
|
|
|
2 |
|
3 |
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
100 |
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
cmT |
|
|
T 5 |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
cm 2,5 R 0,06
134
Значение потерь при неоптимальном управлении (2.13) превышает найденные потери при управлении (2.12) на 30,72 кДж, причём эта разница растёт пропорционально второй степени всех основных рабочих характеристик объекта: угла поворота , момента инерции J и момента сухого трения М.
Найденные выражения для экстремалей (2.11) и (2.12) позволяют реализовать управление электроприводом не только по току i в соответствии с (2.12), но и по частоте вращения ω, подавая на вход системы регулирования задающее воздействие вида (2.11).
omega1 |
0.05s2+0.06s+0.0625 |
i |
|
|
omega |
0.0001s2+0.001s+0.025 |
|
|
|||
|
|
|
|
||
|
|
|
|
|
|
f(u) |
50 |
1 |
0.01 |
1/s |
1/s |
2.5 |
|||||
|
|
0.05s+0.06 |
|
|
|
|
|
|
2.5 |
|
fi |
|
|
|
|
|
Рис. 2.4. Имитационная модель системы оптимального управления поворотом платформы с обратной связью по частоте вращения
На рис. 2.4 представлена имитационная модель такой системы регулирования с обратной связью по частоте вращения.
i, A1500 |
|
|
|
|
|
|
1000 |
|
i(t) по (2.12) |
|
|
i(t) |
|
50 0 |
|
|
|
|||
|
|
|
|
|
||
0 |
|
|
|
|
|
|
-500 |
|
|
|
|
|
|
-1 0000 |
1 |
2 |
3 |
4 |
5 |
6 t, c |
Рис. 2.5а. Изменение тока при оптимальном управлении частотой вращения (2.11)
135
100 |
|
|
|
|
|
|
|
, рад |
|
|
(t) |
|
|
|
|
ω, с-1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
50 |
|
|
|
|
ω(t) |
|
|
30 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
0 |
1 |
2 |
3 |
4 |
5 |
6 |
t, c |
Рис. 2.5б. Изменение частоты вращения ω(t) и угла поворота (t) при оптимальном управлении
по экстремали (2.11)
На рис. 2.5а,б представлены полученные на этой модели фактические законы изменения тока, частоты вращения и угла поворота платформы.
Анализ этих рисунков показывает, что кинематические зависимости частоты вращения и угла поворота практически не отличаются от расчётных, однако отклонение функции тока от экстремали (2.12) уже достаточно заметно и снижает отмеченную выше эффективность оптимального управления. Из этого следует, что непосредственное управление физической величиной, входящей в выражение критерия оптимизации задачи, является более предпочтительным по сравнению с математически эквивалентными, но опосредованными вариантами управления. Важно отметить, что оптимальное решение (2.12) задачи найдено без учёта технических ограничений на максимально допустимое значение тока i(t), иными словами: рассмотренная задача решена без учёта ограничений на управляющее воздействие.
Задача 2. Управление электроприводом подъёмного крана с гибкой связью его подвижной тележки с перемещаемым по горизонтали грузом должно предусматривать демпфи-
136
рование возникающих при этом раскачиваний (колебаний) груза (рис. 2.6).
m1 Fтр F
μ x
L
v
m2
Рис. 2.6. Иллюстрация постановки задачи
На рис. 2.6 обозначены: m1, m2 – массы тележки и груза соответственно, кг; L – длина подвеса груза, м; х – координата горизонтального перемещения каретки, м; – угловое отклонение подвешенного груза от вертикали, рад; μ – коэффициент сухого трения тележки (вязким трением в силу малых значений скорости перемещения подвижных частей системы пренебрегаем); Fтр – сила сухого трения, действующая на каретку, Н.
Ставится задача: найти закон изменения силы F, прикладываемой со стороны электропривода к тележке, обеспечивающий за заданное время Т перемещение груза на расстояние S без раскачиваний, т. е. проекция абсолютной скорости v груза на неподвижную ось х не должна менять своего знака.
Рассматриваемый в задаче объект представим схемой математического маятника со свободной и подвижной по горизонтали точкой подвеса (рис. 2.7), на котором (t)= (t) – частота вращения груза в момент времени t, с-1; Fн – сила на-
тяжения троса, Н; g – ускорение свободного падения, g=9,81 м/с2.
137
Fтр m1 F
Fн
L
Fн
х
m2
Lv
m2g
Рис. 2.7. Расчётная схема объекта “тележка-груз”
Для описания объекта составим уравнения его движения в координатах х и в виде уравнений Лагранжа:
|
E |
|
d |
|
E |
F F ; |
|
|
|
|
|
||||
|
|
|
|
|
|
тр |
|
|
x |
|
|
|
|
||
|
dt x |
|
(2.14) |
||||
E |
d E |
|
|||||
|
m2gLsin , |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
dt |
|
|
где E m1x2 m2v2 – кинетическая энергия системы.
2 2
Получая абсолютную скорость v груза в виде
v2 x2 ( L)2 2x Lcos x2 L2 2 2Lx cos ,
придём к развёрнутому выражению для Е:
|
m1 |
m2 |
2 |
|
m2L2 |
|
2 |
|
|
|
|
||
E |
|
|
x |
|
|
|
|
|
|
m2Lx cos . |
(2.15) |
||
|
2 |
|
2 |
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||
Подставив (2.15) в (2.14), получаем: |
|
||||||||||||
|
|
|
|
|
|
|
|
2 |
sin ) F Fтр |
; |
|||
(m1 m2)x m2L( cos |
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
(2.16) |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
xcos L gsin 0. |
|
|
|
|
|
|
Силу сухого трения Fтр определим в виде Fтр Fн cos ,
где сила натяжения Fн: |
2 |
|
|
|
L. |
(2.17) |
|
Fн m2gcos m2xsin m2 |
138
Система уравнений (2.16), (2.17) полностью описывает движение рассматриваемого объекта.
Имитационная модель, составленная по этим уравнениям в пакете MatLab, представлена на рис. 2.8 для параметров: m1=100 кг; m2=200 кг; L =5 м; Т=5 с; S=10 м; μ=0,01.
Поведение тележки и груза при тривиальном управлении таким объектом (включение привода на постоянный номинальный режим и его выключение при достижении кареткой заданной координаты х) показано на рис. 2.9 и хорошо иллюстрирует колебательность системы “тележка-груз”, причём при отсутствии активного торможения после отключения привода тележка продолжает двигаться по инерции, постепенно теряя избыточно запасённую энергию системы на преодоление силы трения. Для дальнейшей работы с моделью (2.16), (2.17) линеаризуем её для диапазона малых значений и учтём некоторые упрощающие условия:
|
2 |
|
|
|
|
|
|
2 |
L. |
cos |
sin ;cos 1; sin ;gcos xsin |
||||||||
|
Получим |
|
|
|
|
|
|
|
|
|
|
(m |
|
|
|
|
F m g; |
|
|
|
|
m )x |
m L |
|
|
||||
|
|
1 |
2 |
|
2 |
|
2 |
(2.18) |
|
|
|
x L g 0, |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
или, выразив в явном виде x : |
|
|
(2.19) |
||||||
|
|
|
|
|
|
|
|||
|
|
|
x L g , |
||||||
подставив (2.19) в первое уравнение из (2.18): |
|
|
|||||||
|
|
|
|
|
m2)g F m2g . |
(2.20) |
|||
|
|
m1L (m1 |
|||||||
|
Поскольку слагаемое |
μm2g не |
зависит от |
времени, то |
правую часть равенства (2.20) в дальнейшем будем обозначать как одну функцию F (t)=F m2g. Имитационная модель линеаризованной системы представлена на рис. 2.10.
Перейдём к нахождению оптимального управления F (t). В качестве критерия оптимальности, обеспечивающего наименьшую колебательность системы, выберем используемый в теории управления интегральный квадратичный критерий
139