![](/user_photo/_userpic.png)
книги / Управление колебаниями
..pdf§4 ] |
Р А З Г О Н М А Я Т Н И К А П Е Р Е М Е Н Н О Й Д Л И Н Ы |
зи |
систему из состояния покоя при i = 0 в состояние дви жения со скоростью 1> без колебаний. Краевые условия имеют вид
ф(0) = ф(0) = у(0), <р(Т) = ф(Г) = 0, и(Т) = v0. (7.4.3)
В задаче торможения начальный и конечный момен ты времени О, Т в (7.4.3) следует поменять местами.
Введем параметр и и безразмерные переменные
и сделаем замену (7.4.4) в соотношениях (7.4.1)—(7.4.3). Опуская штрихи, получим для задачи разгона
(1 ± цОф ± 2шр + ф = v, 0 ^ у < 1,
(7.4.5)
Ф(0) = <р(0) = и(0) = ф (Г) = ф (Я = 0, v(T) = 1.
Взадаче торможения нужно поменять местами 0 и Г
вкраевых условиях (7.4.5).
Ниже будут построены режимы разгона и торможе ния в случаях подъема и опускания груза. Оптимальность этих режимов не проверяется, однако построенные режи мы обладают минимально возможным числом точек пере ключения п переходят в оптимальные по быстродействию (Т min) в случае маятника постоянной длины (и = 0).
2. Построение режимов разгона и торможения. Опти мальный разгон маятника постоянной длины (и = 0) с
ограничением |
(7.4.5) построен в- п. |
2 § 2 п состоит из |
двух участков |
постоянства скорости. |
При и Ф 0 управле |
ние в задаче разгона будем искать |
в аналогичном виде, |
|
а именно |
|
|
|
|
(7-4.0) |
Здесь £* — постоянная. Спачала построим разгон при опускании груза (знак «+» в (7.4.5)). Уравнение (7.4.5) питегрпруется в бесселевых функциях на интервалах [0, tif) п О*, 74. Сопрягая эти решения в точке t — t^u удовлетворял граничным условиям (7.4.5) с учетом соот ношений (7.4.0), получим всего шесть условий для четы
312 |
Р А З Г О Н К О Л Е Б А Т Е Л Ь Н Ы Х С И С Т Е М |
1ГЛ . 7 |
|
|
рех постоянных интегрирования (на двух интервалах) и для двух параметров **, Т. После преобразований эти ус ловия приводятся к системе двух трансцендентных урав нений относительно параметров р, q
• U M N ^ p )- ЛГ,(гШу>Шд) + В Д Л М =
|
*= l N d q ) U p ) - М я Ш р Ш М , |
|||||
Ш Ш р ) - /| (гШ р )Ш д ) = |
|
|
|
|
|
|
|
= [N,(pWq) - W q W p V X M . |
а |
(7.4.7) |
|||
Здесь введены обозначения |
|
|
|
|
|
|
|
ms) = Jo(s)Ni{s)-Ms)No(s)t |
|
|
|||
г = 2/н, |
р = 2 / 1 + и**/»* 5 = |
2 у Т + И г /м , |
(7.4.8) |
|||
. /о, 7i и TVo, N1 — функции Бесселя и |
Неймана |
соответ |
||||
|
ствующих |
порядков. |
Значе |
|||
|
н и я м и ^ |
в зависимости от и |
||||
|
находились численно на ЭВМ |
|||||
|
путем |
решения |
|
системы |
||
|
(7.4.7) |
относительно парамет |
||||
|
ров р, q. При этом отбира |
|||||
|
лись корпи, соответствующие |
|||||
|
наименьшему Т. |
представле |
||||
|
На |
рис. 7.14 |
||||
|
ны |
результаты |
расчетов — |
|||
|
графики зависимостей £*(и) |
|||||
|
(кривая 7), ТЫ) (кривая 2) |
|||||
и t*(w) = Г М — f*(w) (кривая 3). |
При к = 0 |
в |
соответ |
|||
ствии с |
(7.2.15) имеем £* —t* = |
jt/3, |
Т = 2я/3. |
|
||
Построим реяшм торможения при опускании груза. |
||||||
Сделаем |
в соотношениях (7.4.5) |
|
замену v -*■ 1 — v, ф -*■ |
-н— ф. Уравнение и ограничения при этом не изменятся, а краевые условия (7.4.5), соответствующие разгону, пе рейдут в краевые условия для торможения. Поэтому ре жим торможения получим из (7.4.6), если поменять мес
тами 0 и 1, а именно |
при |
0 < t < t |
|
|
|
[О |
t > T . |
(7.4.9) |
|||
i |
при |
* * < * < Г , |
|||
|
|
/Зависимости t* Ы) и ТЫ) — те же, что и в задаче раз гона (см. рис. 7.14),
§ 4] |
РЛЗГОЫ МАЯТНИКА ПЕРЕМЕННОЙ ДЛИНЫ |
313 |
Обратимся к построению разгона и торможения маят ника в случае подъема груза (знак «—» в (7.4.5)). Вы полним в уравнении (7.4.5) лилейную замену перемен ных t== — Aiti + А2, ср — ^43ф1, v = 1 — и\. Подберем по стоянные Ait А2, А3 так, чтобы уравнение (7.4.5) со зна ком «—» в исходных переменных перешло в уравнение (7.4:5) со знаком «+» в новых переменных ti%qpi, v\. Ока зывается, что этого можно добиться, если замену пере менных осуществить по формулам
ut |
U2 — VT |
U,ф, |
« = - —* + |
- ^ 5-,. |
ф = - р , v = l - v v (7.4.10) |
Здесь и — заданное значение параметра (7.4.4), щ — соответствующий параметр уравнения (7.4.5) со знаком «+» после преобразования (7.4.10).
Чтобы граничные условия (7.4.5) не нарушались при замене (7.4.10), потребуем взаимного соответствия момен тов времени
t = 0 ^ t 1= Т(щ), t = Г (и) «=> = 0. (7.4.11)
Здесь Т(щ) — время разгона при опускании груза для параметра щ (рис. 7.14), а Т'(и) — искомое время раз гона при подъеме груза для параметра и. Из соотноше ний (7.4.10), (7.4.11) получим условия
QЫ = |
Г (и ) = - |
/ i |
+ V K ) |
|
(7.4.12) |
Первое равенство (7.4.12) служит для определения параметра щ по заданному и. Па рис. 7.14 построен гра фик Q(u) (кривая 4). Эта зависимость монотонна, поэтому уравнение Q{u\) = и в (7.4.12) определяет единственное щ > 0.
Второе соотношение (7.4.12) выражает время разгона Т'(и) при подъеме груза. Таким образом, расчет режима разгона сведен к рассчитанному выше разгону при опус кании груза, задаваемому формулами (7.4.6) — (7.4.8) и кривыми рис. 7.14. При этом по заданпому и нужно оп ределить щ. Пересчет осуществляется по формулам (7.4.10) — (7.4.12) и рис. 7.14. Торможение при подъеме груза рассчитывается аналогично режиму торможения (7.4.9) -при опускании груза.
|
Г Л А В А 8 |
НЕКОТОРЫЕ ПРИКЛАДНЫЕ ЗАДАЧИ |
УПРАВЛЕНИИ |
И ОПТИМИЗАЦИИ |
НОЛЕБАНИЙ |
Вданной главе рассмотрен ряд прикладных задач управления и оптимизации для механических колебатель ных систем.
В§§ 1—3 исследуются управляемые колебательные
системы маятникового типа; этот материал по своему содержанию примыкает к главам 6, 7. В § 1 построены способы перемещения маятника с переменной длиной подвеса в вертикальной плоскости. В качестве составных элементов движения используются оптимальные и квазиоптимальные режимы перемещения, найденные в гла ве 6, и режимы разгона и торможения из главы 7. По строенные способы управления позволяют путем совме щения вертикального и горизонтального движений пере местить качающийся груз из одного состояния покоя в другое. § 2 посвящен вопросам управления грузоподъем ными машинами. Здесь обсуждаются механические мо дели таких машин и постановки задач оптимального управления. Отмечено, что реяшмы, построенные в гла вах 6, 7 и в § 1 главы 8, могут быть использованы для автоматизации управления грузоподъемных машин с целью сокращения времени их рабочего цикла. Помимо задач быстродействия, представляют интерес также за дачи минимизации энергетических потерь в электродви гателе. Дано решение таких задач в случае перемещения н разгона груза. В § 3 исследуются некоторые вопросы уп равления системой нескольких маятников путем переме
щения тела, несущего их точки подвеса.
§§ 4—5 посвящены задачам оптимальной амортиза ции колебательных систем. В § 4 приведен краткий об8ор проблем оптимального выбора параметров амортиза ционных систем и дано решение некоторых характерных задач этого класса в случае ударных воздействий. В § 5 исследуется задача оптимальной амортизации роторной системы в процессе ее раскрутки. В отличие от § 4, на-
§ 1] |
ПЕРЕМЕЩЕНИЕ МЛЯТНПКЛ В ПЛОСКОСТИ |
315 |
раметры системы (жесткость амортизации) могут изме няться в процессе движения. Указаны оптимальные за коны изменения параметров, позволяющие снизить неже лательные резонансные эффекты.
Отмстим, что много задач управления колебательны ми системами возникает в области динампки и управле ния роботами и манипуляторами. Эта проблематика, пред ставляющая собой быстро развивающуюся самостоятель ную сферу исследований, в дайной кпиге пе затрагивается (см., например, [161—163, 178, 41, 42, 61, 127, 1281).
§1. Перемещение маятника переменной длины
ввертикальной плоскости
1.Постановка задачи. Рассмотрим задачу о перемеще нии маятника переменной длины в вертикальной плоскос
ти из одного состояния покоя в другое. Эта задача важна
всвязи с управлением грузоподъемными машинами типа мостовых шш козловых крапов и перегружателей. Изложе ние следует работе [141].
Математический маятник (рис. 6.1) может перемещать ся вдоль горизонтальной оси х со скоростью vit), а его длина Lit) — изменяться со скоростью uit).
Требуется переместить маятник на заданное расстоя ние и высоту из состояния покоя в состояние покоя. Го ризонтальное перемещение должно равняться а, а началь ная и конечная длины подвеса равны соответственно L0 и L\. Граничные условия в обозначениях глав 6, 7 запишутся
ввиде
я(0) = у(0) = ф(0) = ф(0) = 0, L(0) = Lo,
(8. 1 . 1 )
xiT) = a, viT) = ф(Я = ср(Г) = О, LiT) = Ь\.
Без ограничения общности будем считать, что маятник опускается: L\ ^ Lo. В противном случае можно обратить время и поменять местами начальную и конечную точки. Кроме того, за счет выбора направления оси х можно счи-, тать, что а > 0. На управляющие функции uit) и vit) на ложены ограничения
0 uil) ^ и0, 0 vit) s? vo, |
(8.1.2) |
где u0и Vo— заданные постоянные.
316 |
П Р И К Л А Д Н Ы Е З А Д А Ч И О П Т И М И З А Ц И И К О Л Е Б А Н И И |
[Г Л . 8 |
Задача состоит в отыскании управляющих функции и и v, удовлетворяющих ограничениям (8.1.2) и реализую щих для' рассматриваемой системы граничные условия (8.1.1). Длительность процесса Т должна быть равна или близка к времени оптимального быстродействия. Постав ленная задача практически важна в связи с управлением грузоподъемными машинами.
Введем безразмерные константы d, А, с и перейдем к безразмерным (штрихованным) переменным х , cp', L
и', v' по формулам
■ |
х = гфоТ1*)"1, |
ф' = (рV Lxg/ь’о, |
I/ = LILU |
||||
|
t' = t/T#, |
Т * = У Ш , |
Т’ = |
т/т* , |
|
||
|
V |
= vlv0, |
и = u/YLrf, |
(8.1.3) |
|||
|
d = a (РоГ#)"1, |
h = |
LJLX, |
с = |
uJYT^g. □ |
||
|
Далее предполагается, что L\ > |
0. Если же |
L\ — 0, то |
||||
Lo ^ L\ = 0 и решение |
задачи оптимального |
быстродей |
ствия очевидно: маятник пмеет нулевую длину и переме щается с максимальной скоростью Vo вдоль оси х.
В переменных (8.3.1) (далее штрихи всюду опущены) уравнения движения в случае малых колебаний и соотно
шения (8.1.1), (8.1.2) примут вид |
|
|
|
■ |
2лр + 2£ф + ф = х, х = |
v, |
L = и,, |
|
х(0) = х{0) = ф(0) =ф (0) = О, |
||
|
£(0) = h ^ l, ( X |
к < |
с, |
|
x(T) = d, ^ Г )= ф (Я |
= ф(Л = 0 , |
|
|
UT) = 1, C X z X l . |
□ (8.1.4) |
Соотношения (8.1.4) содержат три постоянных пара метра d, Л, с, которые связаны с исходными константами формулами (8.1.3).
Точное аналитическое решение задачи оптимального быстродействия (Т -*• min) представляет значительные трудности в силу нелинейности системы уравнений (8.1.4) и ее сравнительно высокой размерности. Реализация на практике законов управления, которые могут быть рассчи таны на ЭВМ, также будет представлять известные труд-
ПЕРЕМЕЩЕНИЕ МАЯТНИКА В ПЛОСКОСТИ |
317 |
ности. Оптимальные управления будут содержать много точек переключения и, кроме того, зависеть от трех па раметров d, h, с. Поэтому ниже предлагаются некоторые достаточно простые (квазиоптпмальпыс) способы управ ления с небольшим числом переключений, удовлетворяю щие (8.1.4). Этп управления построены па основе соче тания оптимальных законов, найденных в главах 6,7 для более простых случаев постоянной длины маятника. Здесь предполагается (см. (8.1.4)), что скорость подвеса ограни чена и может изменяться практически мгновенно. Анало гично при других ограничениях могут быть использованы другие реяшмы глав 6, 7, в частности, при совместных ог раничениях па скорость и ускорение — результаты § 3 главы 7.
2. Простейшие типы движения. Укажем те простейшие движения, из которых будет составлено решение исходной задачи управления. Время t всюду отсчитывается от нача ла соответствующего режима.
A. Опускание груза с максимальной скоростью при от сутствии колебаний п перемещения точки подвеса:
Ф ев 0, х = const, и — с.
B. Оптимальный по быстродействию разгон маятника из состояния покоя до поступательного движения при пос тоянной длине подвеса L —const. В момент окончания
разгона |
накладываются |
условия |
ф(7,в) = ф(7,в) = О, |
||
v(TB) = 1. Это движение согласно (7.2.15) имеет вид |
|||||
|
fl |
при |
0 < t < |
ll2.TB, Гб = 2/зЛ]/Х, |
|
v(t) = |
\ |
|
ЧцТв < |
t < Тв, |
(8.1.5) |
|
lO |
при |
v ('Гв) = 1. |
С. Оптимальное но быстродействию торможение (до покоя) поступательно движущегося маятника при посто янной длине подвеса. Этот режим аналогичен В и задается в виде
0 |
при 0 < .t < .llzTB, |
Тв — |
L, |
v(t) = |
при V2Т в < * < Т в , |
|
(8. 1. 6) |
1 |
V ( T B ) = |
0 . |
D.Горизонтальное перемещение системы без колеба
ний с опусканием груза <р = 0, v —const, и = uit), 0 |
^ |
318 |
П Р И К Л А Д Н Ы Е З А Д А Ч И О П Т И М И З А Ц И И К О Л Е Б А Н И Й [Г Л . 8 |
< 1. Здесь скорость u(t)' может быть произвольной функ цией (0 < и < с).
Е.Оптимальное по быстродействию перемещение груза
из точки в точку прп постоянной |
длине |
подвеса |
(L = |
= const). Это движение состоит из |
2 к + 3 |
участков |
пос |
тоянства скорости и имеет вид (см. формулы (6.2.11), (6.2.31), (6.2.32), где надо положить и = 1, Y = 0)
v(t) = 1 |
при |
i—1 |
< |
* < |
i |
bi |
i = |
1, 3, |
. . . , |
2к + 3, |
|
S |
2 |
||||||||||
|
|
|
i=i |
|
|
J=i |
|
|
|
|
|
|
|
|
i—1 |
Ь< |
t < |
г |
|
|
|
|
|
v (£) = 0 |
при 2 |
S |
|
i = |
2, 4, |
. . . , |
2/c -|- 2, |
||||
|
|
|
j=i |
|
|
;=i |
|
|
|
|
|
|
|
|
|
i;(0) = V (TE ) — 0. |
|
|
|
||||
Целое число к и длины интервалов t,■в принятых обоз-, |
|||||||||||
начениях определяются формулами |
|
|
|
||||||||
|
= |
|
Ъ = - 4 ъ ' |
“ kW = |
«csin[5!£W2).|, |
||||||
t i = |
^гл+з = (т/2 |
a ft) У" |
^2— ^4 = |
•••= h h + 2 — 2aft'|/rX, |
|||||||
|
|
|
|
|
|
|
|
|
|
|
(8.1.7) |
|
|
*3 = |
^Б = |
• • • —^2Ь+1 = |
2 (Л — Gt/jyZ, |
|
|||||
где |
[...] — целая часть числа, |
а |
т(Ь) — единственный ко |
||||||||
рень трансцендентного уравнения |
|
|
|
||||||||
|
6 = 2лк + х - |
2(к + 1)а*(т), |
т е= [0, 2л). (8.1.8) |
||||||||
|
Время быстродействия Т Е задается равенством |
||||||||||
|
|
|
ГЕ= |
2А+3 |
(2л/с + x)VZ. |
|
(8.1.9) |
||||
|
|
|
2 |
tj = |
|
||||||
|
|
|
|
3=1 |
|
|
|
|
|
|
|
В частном случае, когда Ь/(2л) — целое число, ско рость точки подвеса в задаче оптимального быстродей ствия постоянна и вместо формул (8.1.7) — (8.1.9) имеем
|
v = 1 при 0 < t < |
Те, |
TE = d. |
(8.1.10) |
|
F. |
Квазиоптпмальпое по |
быстродействию |
перемеще |
||
ние груза из точки в точку при постоянной длине под |
|||||
веса |
L = const. Это движение, построенное в § 2 |
главы |
|||
6 (см. формулы (6.3.15), (6.3.16)), |
отличается |
от |
опти- |
§ И ПЕРЕМЕЩЕНИЕ МАЯТНИКА В ПЛОСКОСТИ 319
мального том, что содержит три участка постоянства ско рости. Оно задается соотношениями
|
|
J1 при |
0 < |
%< |
ilf |
tL+ U< |
t < |
TF, |
|
|||
" |
u |
\0 при |
tAc |
t < |
t* + |
*8, |
* = |
0, |
t = |
Tjr, |
||
|
|
|
|
_ |
TF= |
3 |
tj, |
к = |
[6/(2я)1, |
|||
|
|
tL= i3 -I- 2nkVL, |
S |
|||||||||
|
| = b — 2я/.-, |
ii = 0, |
|
2V = |
2я/с 1 /Z/ |
при |
t = |
0, |
||||
т| = |
я -f |/2, |
= (2яА: + |
я + |
V2)VL при 0 < ^ < 2л1 |
||||||||
|
|
i2 = |
0, |
i3 = 0 |
при 0 < т | < я , |
|
|
|||||
t., = |
(2я — rj)}/ L, |
t3 = |
(rj — n)YL при |
я |
г) < |
2я. |
||||||
|
|
|
|
|
|
|
|
|
|
|
□ |
(8.1.11) |
Из перечисленных движении можно построить три простых способа управления, решающих поставленную задачу о перемещении груза, а именно движения: АЕА, AFA, ABDCA.
Режимы АЕА и AFA содержит один свободный па раметр — длину L па участках Е и F соответственно. В режиме ABDCA имеются два свободные параметра — по стоянные длины Li и Ь2 при движепиях В, С. Эти пара метры естественно выбрать так, чтобы минимизировать суммарное время перемещения.
3. Построение управления в виде оптимального соче тания исходных режимов. При реализации режима АЕА суммарное время перемещения Т\ согласно соотношениям (8.1.4), (8.1.0) задастся формулой Т\ = (1 — h)/c + ТЕШ. Параметр L найдем из условия minTi по L при ограниче нии h ^ L < 1.
Первое слагаемое в формуле для Т\ не зависит от L, поэтому достаточно пайтн минимум ТЕиз (8.1.9). ПреОбразуем ТЕк виду Те = сЩЬ), где использованы обозначения (8.1.3), (8.1.7) и
il)(b)= iM b)/6, i|)* (6) |
= 2я/с(Ь) + т(Ь) .(&>0). (8,1.12) |
|
Из выражения (8.1.7) |
для b и неравенств |
h < L ^ i |
получим ограничения па b вида |
Та |
|
ким образом, исходная задача минимизации |
Т\ сведена |
в терминах новой переменной |
Ъ к определению minij>(&) |
по Ьпри ограничениях d<1 b^ |
d]/Tl. |
320 ПРИКЛАДНЫЕ ЗАДАЧИ ОПТИМИЗАЦИИ КОЛЕБЛИИЙ [Г Л . 8
Функция ф*(6) из (8.1.12) задаст время оптимального быстродействия как функцию расстояния. Эта функция изу чена в § 2 главы 6 (где она обозначена через ТЫ)). Функ ция ф*(&) строго возрастает, и справедливы соотношения
ф *(Ь)^Ь при 6 > 0 , |
|
|
|
|||
ф* (6) = |
л [1 + Ы{2я)] при |
Ъ< |
2я, |
(8.1.13) |
||
ф* ( 2 J U ) = |
2ni при |
i = |
1 ,2 , .. . , |
lim ф* (Ъ)/Ь= 1. |
||
|
|
|
|
Ь-»оо |
|
|
На основании свойств (8.1.13) и соотношений (8.1.12) |
||||||
имеем (см. рис. 8.1, сплошная кривая 1) |
|
|||||
ф (Ь)> 1 при & > 0, |
ф(6) = |
1/2 + |
я/6 |
при Ь < 2я, |
||
|
|
|
|
|
|
(8.1.14) |
lim ф(Ь) = 1, |
ф (bt) = |
1, |
= |
2JU, |
i = |
1, 2, 3, .. . |
На интервалах b{ < |
6 < 6i+j функция ф(Ь) имеет един |
|||||
ственный внутренний максимум. Отсюда |
и из неравенств |
ФУ
W
IJr
1,0
Рис. 8.1.
d^b^d/YIn вытекает, что при d < 2яУАимеем Ъ< 2л, и
пппф(&) достигается в точке |
b — dfh. Это соответствует |
режиму управления с тремя |
(к = 0) участками постоян |
ства скорости v(t). __ |
|
Если же d> 2nYh, то оптимальное управление v(t) мо жет содержать один либо более трех участков постоянства скорости. Если при некотором i > 1 выполняется включе ние d/j/H], то на отрезке [d, d/Yh\ имеется хотя бы одна точка Ь<, где достигается абсолютный минимум функ ции ф(Ь), равный ф(Ь<) = 1 (см. рис. 8.1). В качестве оптимального значения Ъ здесь следует выбрать любое