Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Учебное пособие 1920

.pdf
Скачиваний:
14
Добавлен:
30.04.2022
Размер:
2.9 Mб
Скачать

 

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) получим с01. То-

гда совместное решение (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