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

книги / Управление колебаниями

..pdf
Скачиваний:
2
Добавлен:
12.11.2023
Размер:
15.93 Mб
Скачать

§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). В качестве оптимального значения Ъ здесь следует выбрать любое