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

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

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

§ 4

А М О Р Т И З А Ц И Я Д И Н А М И Ч Е С К И Х С И С Т Е М

351

З а д ач а 4.

Найти оптимальную характеристику

такую, что

v4x, г ) е У

 

 

 

 

 

шах Л (г;0, F) =

min max J.2 (и, F),

 

 

re ф

нет геФ

 

 

max

(r-°, F ) ^ . R .

 

 

ГёФ

 

 

 

Здесь

V и I i — заданные положительные

величины.

Решению задач об оптимальной амортизации систем, рас­ считанных на класс внешних воздействий, посвящены ра­ боты L26G, 189, 51] и другие. Задачи 1—4 могут быть обобщены на системы со многими степенями свободы (см. [75, 109, 266] и другие).

Перечисленные выше задачи являются простыми при­ мерами, отражающими типичные особенности оптимиза­ ционных задач, связанных с проектированием амортиза­ ционных систем. Опи, конечно, не исчерпывают всех воз­ можных постановок. Читателю, интересующемуся пробле­ мами проектирования оптимальных амортизаторов в более широком аспекте, можно обратиться к монографиям [109, 207, 266]. Ниже рассматриваются в качестве примеров некоторые простые задачи оптимальной амортизации для систем с одной степенью свободы.

2. Оптимизация параметров упруго-демпфированных противоударных амортизаторов. В книге [97] показано, что амортизаторы целесообразно применять лишь в том случае, когда путь торможения (или разгона) объекта, содержащего амортизируемое тело, не больше хода амор­ тизатора. Поэтому особенный интерес представляет воп­ рос о защите при кратковременном (ударном) внешнем воздействии па корпус.

Если времл действия внешней силы много меньше ха­ рактерного времени амортизационной системы (например, периода ее колебании), то в некоторых случаях можно считать, что корпус подвергается мгповеппому удару и полагать Fit) = {}6U), где (} — постоянная (интенсивность удара), a 6U) - дельта-функция.

Рассмотрим механическую систему, представляющую собой абсолютно твердое тело, связанное с движущимся корпусом посредством упруго-демлфировапиого амортиза­

352 П Р И К Л А Д Н Ы Е З А Д А Ч И О П Т И М И З А Ц И И К О Л Е Б А Н И Й

(ГЛ . 8

тора со степенной характеристикой вида

 

и(х, х) = k\x\rsignх + c|.r|wsign х,

(8.4.7)

/• > 0, т > 0.

 

Здесь к > 0 н с ** 0 — постоянные коэффициенты демп­ фирования и жесткости соответственно, г и т — копстап-

ты. Предполагается, что в начальный момент

времени

t 0 система,

находящаяся

в положении равновесия

ж(0) = ,т(0) = 0,

подвергается

мгновенному удару

{Fit) =

= (М5Ш), в результате которого амортизируемое тело при­ обретает конечную относительную скорость (3. Движение амортизируемого тела относительно корпуса определяет­ ся дифференциальным уравнением с начальными усло­ виями

х + k\x\rsign х + сЫ '" sign х = 0,

 

 

 

(8.4.8)

 

я:С0) = 0, кО) =

р.

 

Введем следующие обозначения

 

 

I i m(k ,c)=

max |z(J)|,

 

(8.4.9)

 

te[o,«o

 

 

c) =

 

 

 

= max

\k\x (t) |T sign x (t) + c |x (t) |m sign x (t) |=

te[o,oo)

 

 

 

=

max |x(i)|.

(8.4.10)

 

 

teto.oo)

 

Ставится задача об определении оптимальных коэф­ фициентов демпфирования и жесткости, минимизирую­ щих максимум модуля перегрузки при ограниченном мак­ симуме модуля отклонения. Таким образом, требуется пайти параметры /со, с0 такие, что

Лг,ж (ко. Се) = min/;• "(*, С),

ft,c>o

(8.4.11)

11 (к01со)

Такая постановка соответствует задаче 2, описанной в п. 1. В рассматриваемом случае множество Y допусти­ мых характеристик представляет собой параметрическое

9 41

АМОРТИЗАЦИЯ ДИНАМИЧЕСКИХ

СИСТЕМ

353

семейство

функций,

а внешнее воздействие — ударного

типа IFit) = рб(£)).

Будем считать в

дальнейшем,

что

начальная скорость [} и максимально допустимое отклоне­ ние D равняются единице. Это отвечает переходу в

(8.4.8)— (8.4.11)

к

безразмерным

переменным

по фор­

мулам

 

 

 

 

 

*' = 7Г'

V =

\Г 4'

k' = bD f~\

=

(8-4-12)

Отметим некоторые свойства системы (8.4.8), которые

понадобятся в дальнейшем.

 

 

С в о й с т в о

1. Максимум модуля отклонения аморти­

зируемого

тела

(8.4.9) достигается

в момент

£* первого

локального экстремума функции x(t). Это свойство явля­

ется простым

следствием диссипативиости системы (8.4.8).

С в о й с т в о

2. Максимум модуля перегрузки (8.4.10)

достигается на

отрезке 0 <£?^£*, где £ * — момент пер­

вого локального экстремума функции я(£). Действительно, в силу свойства 1 для любого момента

времени £ > £* выполняется неравенство |х (£) J^ |х (£*) |. Поскольку функция хШ непрерывна и £* — точка ее эк­ стремума, то для любого момента £j > £* существует мо­

мент £2 < £* такой, что

 

I*(£I)I = I*(£2)I.

(8.4.13)

В силу диссипативиости системы (8.4.8)

ее механиче­

ская энергия Е = хУ2 + сЫ т+1/(тп + 1) не возрастает с ростом времени и, следовательно, справедливо неравен­ ство

|i(£i)l^|i(£2)l. (8.4.14)

Решение задачи Коши (8.4.8) ж(£) неотрицательно

вместе со своей производной xit) на отрезке 0 ^ £ ^ £ *. Отсюда и из (8.4.13), (8.4.14) следует, что

|'iUi)| </с|г(£1)1г + с|ж(£1)|т ^А:а:г(£2) + са:т (£2) = |а?(£2)|.

Таким образом, для любого момента времени tx>►£+ справедливо соотношение

|i(£i)| < Ы£г)1 *Smax|z(£)l, <6Е0.

которое и доказывает свойство 2 .

23 ф. л. Черноусыю, Л. Д, Акуленко. В. Н. Соколов

354 ПРИКЛАДНЫЕ ЗАДАЧИ ОПТИМИЗАЦИИ КОЛЕБАНИЙ [ГЛ. 8

Приведем решение задачи оптимизации (8.4.11) для амортизаторов с линейной жесткостью (m. = 1 ) и с линей­ ным (г = 1 ) и квадратичным (г = 2 ) демпфированием. Амортизаторы указанных типов широко распространены

втехнике.

3.Амортизатор с квадратичным демпфированием.

Прежде, чем решать задачу определения оптимальных па­ раметров амортизатора с квадратичным демпфированием, решим следующую задачу о предельных возможностях защиты от ударных воздействий. Пусть движение систёмы описывается дифференциальным уравнением с на­ чальными условиями

я + н(£) = 0, .т(0) = 0, я(0) = 1.

(8.4.15)

Требуется найти кусочно непрерывную функцию u0it) такую, что

max

|и0(t) |=

min

max |и (t) |,

(8.4.16)

ieEo,°°)

|я (t) I ^

и

*=[0.00)

max

1..

 

 

<е[о,оо)

 

 

 

Справедливо следующее соотношение .

 

 

max |и0(t) |^ 0,5.

(8.4.17)

 

teEo.oo)

 

 

 

Допустим противное, т. е. что lifo(i)l<0,5 (и, следо­ вательно Uo(i)<0,5) для всех i s [0, °°). Тогда справед­ лива оценка решения уравнения (8.4.15)

 

t

 

 

 

t

 

 

а

z(i) = t -

j ( i - T ) u ( T ) d T > t - 0 , 5 J(* — * ) t fc = f —

 

о

x (i) >

 

0

 

 

 

 

max

max (i — i2/4) =

1.

 

 

 

*=[0,°o)

 

iS[0,oo)

 

 

 

 

Отсюда

вытекает,

что

если

I щ (i) I <

0,5

для

всех

i s [О, «О, то неравенство

(8.4.16)

не выполняется. Полу­

ченное противоречие

доказывает

неравенство

(8.4.17).

Можно проверить, что управление

 

 

 

 

, . , ( 0 , 5 ,

 

 

 

 

 

 

 

UoU

И *):

|i>(«)|<0,5,

|*(<}|< 1 ,

« >

2

(8.4.18)

§ 41

АМОРТИЗАЦИЯ

ДИНАМИЧЕСКИХ СИСТЕМ

355

приводит к соотношениям

 

 

 

max

K ( t ) 1 =

0,5,

max |*(*)| = я (2) = 1.

(8.4.19)

ie [o ,o o )

 

te E o .o o )

 

 

 

Выражепис (8.4.18)

означает, что.на

отрезке 0=^ f^ 2

до достижения

первого

нуля

скорости

x(t) управление

uyit)

постоянно,

а

при

t > 2

оно может

быть любой ку-

сочно-непрерывпой

функцией, удовлетворяющей

указан­

ным в (8.4.18) неравенствам. Из (8.4.17), (8.4.19) выте­ кает, что любая функция н0(£) вида (8.4.18) является оп­ тимальным управлением, решающим задачу (8.4.15), (8.4.16). Методом от противного доказывается, что на от­ резке 0 t ^ 2 оптимальное управление определяется единственным образом.

Соотношения (8.4.19) определяют предельные возмож­ ности защиты амортизируемого тела от ударных воздей­ ствии. Нельзя добиться значения максимальной перегруз­ ки, меньшего 0,5, в задаче (8.4.15), (8.4.16).

Покажем теперь, что можно так выбрать параметры амортизатора с линейной жесткостью и квадратичным демпфированием, что минимум максимальной нагрузки

будет

равен 0,5, т. е. что при соответствующем выборе

параметров амортизатор указанного тппа обеспечивает

абсолютный минимум перегрузки.

 

 

Если управление u0(t) имеет вид (8.4.18), то на отрез­

ке 0 < t ^

2 имеем

 

 

 

 

x(t) = t - m , k t) =

l-t/ 2.

(8.4.20)

После

подстановки выражении

(8.4.20)

в уравнение

(8 .4 .8 )

при т = 1 , г = 2 левая часть этого уравнения при­

мет вид

 

 

 

- т + * ( * - т ) + ^ - - г ) - ( * - т ) + < * - с>( т - ‘ ) '

0 < * < 2 .

Это выражение

тождественно обращается в нуль тог­

да и только тогда,

когда к — с = 0,5. Таким образом, до­

казано, что

при

к — с —0,5

выполняются

-соотношения

|я(*)| =

0 ,5 ,

 

0 < * < 2 ,

max И*)| =

s (2 ) - 1 .

 

 

 

 

te[o,a]

 

23*

356 ПРИКЛАДНЫЕ ЗАДАЧИ ОПТИМИЗАЦИЙ КОЛЕБАНИЙ [ГЛ. 8

Отсюда и из отмеченных выше свойств 1 и 2 системы (8.4.8) (здесь имеем £* = 2 ) вытекает, что для амортиза­ тора с квадратичным демпфированием и линейной жест­ костью оптимальные параметры равны fco = c0 = O,5. Пе­ реходя к исходным размерным переменным по формулам (8.4.12), получим зависимость оптимальных параметров и соответствующего значения функционала от величии Р п D

*„ = 0,5-1, с0 = 0,5 -£ ., II'1(7с0, с„) = 0,5

(8.4.21)

Итак, амортизатор с линейной жесткостью и квадра­ тичным демпфированием обеспечивает абсолютный мини­ мум максимальной перегрузки. Отметим, что в классе упруго-демпфированных амортизаторов с характеристи­ кой вида (8.4.7) только амортизатор с квадратичным демп­ фированием (г = 2 ) и линейной жесткостью = 1 ) обла­ дает таким свойством. Действительно, как было отмечено выше, оптимальное управление и0(£), доставляющее пре­ дельно возможное качество амортизации, определяется единственным образом на отрезке 0 < £ ^ 2 (см. (8.4.18)). Поэтому для того, чтобы амортизатор с характеристикой вида (8.4.7) обеспечивал абсолютный минимум перегруз­ ки, необходимо, чтобы выполнялось тождество

* ( 1 _ т ) Г + с (‘ - т Г " т

-

<8-4-22)

Подставляя в (8.4.22) значения

t = 0 и t =

2, получим,

что коэффициенты демпфирования и жесткости должны быть равны к = с = 0,5. Дифференцируя тождество (8.4.22) при к = с = 0,5, получим новое тождество

которое в случае

г > 0 , пг> 0

справедливо только при

г —2 , т—1.

 

амортизатор. Решив при

m =r= 1 зада­

4.

Линейный

чу Коши

(8.4.8)

и

последовав

на максимум

функции

(£)1 и 1я(£)1 при помощи обычных методов анализа, по-

§ «Г

АМОРТИЗАЦИЯ ДИНАМИЧЕСКИХ СИСТЕМ

357

лучим выражения для функционалов (8.4.9), (8.4.10)

 

1

f Ic-Vk* — 4с\ гУл2- 4С f 2 — 4 0

0,

 

У с

\ Л + / ^

I

 

 

 

 

 

_ 4с J

 

 

 

 

И г(Ьс) = . 2{(кеГ\

 

 

 

к2— 4с =

О,

 

 

 

 

ft

.

V 4с -

ft2

 

 

17г

е х р ( - V~4^W

arctg------г-----

 

 

 

 

 

 

 

2 — 4с<С0,

 

 

 

(8.4.23)

 

7с,

 

 

 

 

 

 

 

 

 

7с8 — с >

О,

 

 

е х р// —

*

X (/с, с) [/с cos х (/с, с) -

 

1\'х{к, с) =

У 4с — fts

 

 

 

 

 

+

 

sin х (7с, с)],

 

к2 — с С О,

 

 

У 4с — Л2

- к2)У 4с-

 

 

 

 

X (Тс, с) =

arctg

к2

(8.4.24)

 

 

 

 

 

 

 

Формулы (8.4.23), (8.4.24) получены и исследованы в

работе

[491. Доказано, что

grad

(к, с) Ф 0

всюду в

области

/с^О, с^ О .

Следовательно, оптимальные пара­

метры находятся на границе области своих допустимых

значений, т. е. на кривой у =

{&, с : /^ д (/с, с) =

1 }.

Отметим некоторые свойства кривой у. Дифференци­

рованием функции

II’1(к, с)

и последующим анализом

доказывается, что

 

fot’1 (к, с)

 

 

d l1'1 (ft, с)

 

 

дк

< о,

дс

< 0 .

(8.4.25)

Неравенства (8.4.25) имеют простое механическое ис­ толкование. Они означают, что с увеличением как коэф­ фициента демпфирования, так и жесткости амортизатора, максимум модуля отклонения амортизируемого тела уменьшается. Из (8.4.23) следует, что

II'1(к, с)

у - , / } д (к, с)

j .

(8.4.26)

Рис. 8.12.

358 ПР1ШЛАДЦЫЕ ЗАДАЧИ ОПТИМИЗАЦИИ КОЛЕБАН ИЙ {ГЛ. 8

В силу теоремы о неявных функциях из (8.4.25) вы­ текает, что кривая Y представляет собой график монотон­ но убывающей функции с {к), и в силу (8.4.26) выполня­ ются соотношения

lim с (к) = 1 ,

Ига. с (к) =

0 .

 

 

 

 

ft-*0

h~*1

 

 

 

 

 

 

 

На рис. 8 .1 2

сплошными линиями изображены линии

уровня функции

’1 (^, с) , а штриховой — кривая 4 . За­

 

 

дача

поиска

оптимальных

 

 

параметров

сводится к по­

 

 

иску

минимума

функции

 

 

1\л (к, с)

на

кривой

у.

 

 

В

результате

численного

 

 

решения на ЭВМ получе­

 

 

ны

 

следующие

значения

 

 

оптимальных

параметров

 

 

и

соответствующего

им

 

 

минимума

максимальной

 

 

перегрузки

 

 

 

 

 

 

/с0«

0,481,

 

с0 «0 ,3 6 1,

 

 

 

 

(*..<>.)* 0,521.

 

Переходя к исходным размерным переменным по формулам (8.4.12), получим выражения, определяющие

зависимость оптимальных параметров линейного аморти­ затора ц функционала от начальной скорости р н макси­ мально допустимого отклонения D

0,481 с0* 0 , 3 6 1 0 , 5 2 1 $ .

(8.4.27)

Сравним значение критерия качества амортизации, со­ ответствующего оптимальному линейному амортизатору,

спредельно возможным, обеспечиваемым амортизатором

слинейной жесткостью и квадратичным демпфированием. Из (8.4.21), (8.4.27) следует, что

^

f a °.) .. 0,04,

Пд (‘ „ .д

§ 5]

УПРАВЛЯЕМ АЯ - АМОРТИЗАЦИЯ РОТОРА

359

т. е. линейный амортизатор с оптимальными параметрами обеспечивает хорошее качество амортизации, отличаю­ щееся от предельно возможного всего па 4%. Некоторые задачи амортизации крутильных колебаний в случае удар­ ного воздействия, близкие по постановке к изложенным выше, решены в статье [52].

§5. Об управляемой амортизации ротора

1.Уравнения движения и квазнстацпонарное прибли­ жение. Рассматривается механическая система, состоя­ щая из несбалансированного ротора, вращающегося во­

круг оси, жестко связанной с корпусом. Корпус укреплен на неподвижном основании посредством вязко-упругого амортизатора с линейной характеристикой и может пере­ мещаться поступательно вдоль горизонтальной оси х (рис. 8.13). Пусть М — масса корпуса, т — масса ротора,I

I— расстояние от оси вращения до центра инерцпи рото­ ра; & > 0 , с > 0 — коэффициенты демпфирования и жест­ кости амортизатора соответственно; х — отклонение кор­ пуса от положения равновесия, <р — угол поворота ротора вокруг оси вращения. Тогда имеем уравнение движения центра масс системы

+ М)х + кх + cx = ml(<p2 sin ф — фcos ф).

(8.5.1)

Изменение угловой переменной ф в процессе раскрут­

ки зададим уравнением и начальными условиями

 

Ф = со, со = в/ ( £ ) , ф(0) = ф°, со'(0) = со0,

(8.5.2)

где ф°, ю° — постоянные, в < 1 — малый параметр. Функ­ ция / и ее производная df/dt предполагаются непрерыв­ ными и ограниченными для всех £ е [0, «>). Конкретный

360 ПРИКЛАДНЫЕ ЗАДАЧИ ОПТИМИЗАЦИИ КОЛЕБАНИЙ (РЛ. 8

их вид далее не используется; предполагается лишь, что

угловая скорость © в

силу

(8.5.2) монотонно

возрастает

от а 0 до некоторой скорости установившегося

вращения

ю1, т. е. © а 1 при £

 

Малость числового параметра

е означает, что время, за которое существенно изменяет­ ся угловая скорость вращения ротора, много больше пе­ риода установившегося вращения 2 л/© 1. Приведенная модель описывает подвижные узлы машин и механизмов, широко распространенных в технике. Отметим, что дина­ мика колебательных систем, содержащих неуравновешен­ ный ротор, исследовалась в ряде работ (см., например, [111, 144, 203, 222, 240, 255, 256] и другие).

Перейдем в

(8.5.1), (8.5.2) к безразмерным величинам

 

 

т + м

*

М =

-|- М) 0 1

+

М) (со1 )3

 

 

 

(8.5.3)

« " - г

(ш0)' = - £

-

 

После замепы (8.5.3) получим (штрихи далее опус­ каем)

х + кх + сх == qr sm ф — ф cos ф,

(8.5.4)

Ф = ©, © = е/(£), ф(0) = ф°, ©(0) = ©°.

Решение уравнения (8.5.4) для х представим в виде суммы свободных и вынужденных колебаний

x(t) = CiGxit) + C2G2{t) + уit).

(8.5.5)

Здесь 6?i, G2 — функции, отвечающие решению одно­ родного уравнения (8.5.4), Си С2— постоянные интегри­ рования. Вынужденные колебания у(£) в случав медлен­ но изменяющейся угловой скорости приближенно пред­ ставляются в виде

л у (i) =

U (с, k, z) sin (ф + а),

U{с, к, z) =

z [(с — z f + k2z]~112,

sin а = B/TJ,

COS а = A/Ut z = ©2,