![](/user_photo/_userpic.png)
книги / Управление колебаниями
..pdf§ 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) |
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, |