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

книги из ГПНТБ / Колебания и устойчивость упругих систем машин и приборов

..pdf
Скачиваний:
6
Добавлен:
24.10.2023
Размер:
8.97 Mб
Скачать

2.ТЕОРЕТИЧЕСКОЕ РАССМОТРЕНИЕ ЗАДАЧИ О

САМО ВОЗБУЖДАЮЩЕМСЯ СУБГАРМОНИЧЕСКОМ РЕЗОНАНСЕ .

Из резонансных соотношений (4.27) рассмотрим первое

Ш п

——— =2, соответствующее точному делению частоты внешне-

' еошо

со возбуждения в два раза. Ограничим рассмотрение слу­ чаем, 'когда характеристика R(x) представлена в виде куби­ ческого полинома, так как в дальнейшем, исходя из оценки влияния неучитываемых высших степеней R{x) на демультипликацнонный резонанс, покажем, что это влияние не ка­ чественного, а только количественного характера. Положив со= 2ш1, движение ооцнллятор-а опишем уравнением

dt2

+25 ^ + w5(e0+ e 1x + s2x2)A'= /0 sin 2(üJ,

(4.28)

dt

 

где 25—коэффициент затухания,

 

F0—амплитуда внешней силы.

Зная характер изменения параметра из (4.21), решение (4.28) (можно сразу принять в форме

х= Хsin 2(іѵ^+а sin u j + b cos

zD

(4.29)

и, подставив его в (4.28), найти а, b и 20,

воспользовавшись

одним из приближенных методов нелинейной механики, на­

пример

методом Ритца — Галеірікінна или гармонического ба­

ланса.

В решении (4.29), по-прежнему называемом «глав­

ным»,

амплитуда вынужденного решения Я (.амплитуда на­

качки) определяется приближенно из линеаризованного урав­

нения (4.28)

 

 

4 т г -

+ < еоИ=А) sin Зшр1,

(4-30)

dt2

 

 

в предположении, что

затухание мало, а нелинейность

R(x)

не очень велика. Тогда из (4.30)

с точностью до величин по­

рядка

ві, 62 имеем

 

 

 

Wo/*o«>3 - 4ы2.

(4.31)

Найденные после подстановки (4.29) в

(4.28) решения, пос­

кольку

они действительны, будут

приближенными решения-

2 2 9

іми задачи, и дальнейшие выкладки к этому и сводятся. Представляется, однако, более желательным процесс подста­ новки (4.29) в (4.28) выполнить в два этана, используя пе­ реход с помощью преобразования

х=Х sin 2шх; ~ z0+ z

(4-32)

ж уравнению возмущенного движения

d-z

dz

е о Н~ —

s.2/l2— 2 г 1г0 + 3 е 2г^-1-

---- +2о — +

о)-

dP

dt

 

 

 

 

+ 2 (bjl—3e,zu) Xsin 2шxt — — e.,X2cos 4ш1і

z-f

 

 

-\-Y[z,

X, z0,

t) = 0,

(4,33)

где Y(z, X, z0, t) — нелинейная функция, учитывающая чле­ ны высших порядков. Линейная часть уравнения (4.33), пред­ ставляющая собой уравнение в вариациях типа Хилла, отра­ жает характер модуляции жесткости стационарными колеба­ ниями вынужденного режима

,\:1 = Хsin 2й)^—20.

(4.34)

В частности, первая гармоника

2 ( — Зе2 z0) X sin 2ш-^

определяет величину приращения крутизны характеристики за период колебания X sin 2w1t—z0. Глубина модуляции парамет­ ра определяется разностью между максимальным МшІ[е0- -2гхг0+

+ 3£22o-(-2 (£j— Зе2г0)] X и минимальным

Мщ [s0— 2s1z0 -f 3e2z2

—2(s1—3s2z0)]X значениями жесткости,

отнесённой к среднему

ее

значению. Амплитуду 2(ег—3s2z0)X первой гармоники при изме­ нении параметра будем называть коэффициентом возбуждения

в первой области неустойчивости вынужденных колебаний, так как в первом приближении именно эта гармоническая составля­

ющая создает комбинационное

воздействие

[73] на систему при

частоте

 

 

 

 

ш1= 2 й ) 1-- (і)1

( 4 .3 5 )

субгарімоники

z = a sin w j

4- b coscö^ .

Учитывая, что

щи одна из высших гармоник, имеющих место в спектре ко­ лебаний контура, за -исключением .составляющей с частотой ші, не попадает в полосу пропускания контура, будем счи­ тать их амплитуды імнлыми.

230

Для ■решения (4.33)

воспользуемся

методом Ван-дер-

Поля, полагая

 

 

z = a{t) sin

cos (nxt,

(4.36)

где a(t), b{t) — функции, медленно меняющиеся во времени з период установления и постоянные для стационарных решений. Медленность изменения а, b, z0 означает, что

da

d-a

_ da

< a ,

----<

— ,

dt

dt2

dt

db

d-b

db

< b,

---

dt

dt2

dt

Подставим (4.36) в (4.33) и осуществим ряд стандартных операций по его укорачиванию. Прежде всего отбросим нерезонаноиые гармоники, затем, пользуясь малостью затуха­ ния ' 25, исключим из уравнений члены второго ,и более вы­

сокого порядка малости, т. >е.

25а

и 25Ь,

и вторые произ-

^

d2a

d2b

Приравняв коэф­

водные амплитуды еѵбгармоники----, —

 

dt2

dt2

 

фициенты при sin , cos

и постоянной нулю, получим

следующие два укороченных уравнения установления колеба­ ний:

 

 

 

 

 

 

éjX—3e2Xz0-j- ~H- х- I

a +

+

( — 1 —?+e0“t~ — г^ 2—

 

— e2^ 2+ 3e2zj'SI b

(4.37)

 

. . . .

db

1

'1—£ + so~b -g- èa^2—^Zg-b

 

ф(а,

b)S W l _ = _

+

e2A2+ 3 e 2Z5)

a + (sjX—3e2Xz0-

25ü),

 

(4.38)

0),T

 

 

 

 

 

 

 

 

 

 

и одно уравнение связи амплитуд

субгармоиики

А —У a2+62 ,

накачки X

и смещения z0, вида

 

 

 

 

 

 

 

£2zo~'éiZ<H~ (ео +

-f- Y Ч

20=

 

 

 

 

 

SiX2 ,

М 2

— sАаЬ.

 

(4.39)

 

 

 

 

— +

—±

+

 

 

 

 

 

2

2

 

2

 

 

231

В

уравнениях

(4.37)

и

(4.38)

через

$

обозначена

расстрой­

ка

половины

частоты

.внешнего воздействия

относительно

(собственной частоты

ш0, т. е.

 

 

 

 

 

 

 

 

 

W“

 

 

 

 

(4.40)

 

 

 

 

 

 

 

 

 

Стационарные решения задачи получим из

(4.37)

и (4.38)

как решения системы при

 

 

 

 

 

 

 

= 0;

ё±

,

0,

 

 

то

есть

dt

 

dt.

 

 

 

 

 

 

 

 

 

 

 

 

 

^£jX—3s2Xz0+

1 j

а - f- ^— 1 — ^- Ь + -2_e2^2—^eizo_b

 

 

+

M 3+

3e22sj

 

0,

 

(4.41)

^1 — (H-s0+ e2^2—- £izo~b — e2^ 2_b 3e2Zoj a-\-

^ l - 3 e 2Xz0- ~ ^ j b = 0 .

(4.42)

Одним из возможных решений системы (4.41)

и (4.42) явля­

ется а—Ь = 0- или лу= Хsin 2w^—z0, при котором парамет­

рическая генерация отсутствует. Таким образом, ввиду того,

что уравнения (4.41)

и (4.42)

всегда

удовлетворяются

при

а = &= 0, возможным решением задачи

является

в с е г д а

ее

в ы н у ж д е н н о е

р е ше в ие ,

определяемое

X sin 2ш^—г0.

Нетривиальное решение [93] получим, приравняв нулю

детерминант коэффициентов при

а и Ь, из которого

 

А2— —— Г^-(-1—£0—— £2X2-j-2e1z0—Зе22^:

 

^£2

L

^

 

 

 

 

± у ^ ( в 1Х -

3s2Xz0)2-

— (!;+1) ■

(4.43)

Динамическое смещение z0 должно удовлетворять уравнению (4.39) , которое упрощается, если положить Ь = 0, т. е. считат влияние затухания на z0 несущественным. Тогда ab = 0 « из (4.39) поломаем

232

3

3

г0= ±L (Х*+Л2). (4.44)

e220 ~ Sl Zo+ ( е 0 4 - — 82^2 4"

£2 ^ 2

Представляет 'Интерес случай, когда характеристика R{x) системы до смещения, т. с. R(x) соответствующей симмет­ ричной системы, выражена -кубическим полиномом R — K(x+ -f ß+). Тогда уравнения (4.43) и (4.44) удается выразить в безразмерной форме. В этом случае

 

£0= l-}-3j3A2; е, —3ßA;

ea=ß.

(4.45)

Подставив (4.45)

в (4.43) и (4.44) и умножив

все парамет­

ры на )/ß, .получим

 

 

 

 

 

Л2=

5—3(д—z0)3— ~

х2±

 

 

 

3

 

 

 

 

 

 

- 2~о)2- 4 - 4 іе+ і)

(4.46)

 

 

 

 

шо

 

 

—о

2г0[1+ 2„-ТЗД (Д—z0) ]—3 (Д—z0)X2

(4.47)

ДІ

_ _

---------------- - —

-------------------

 

 

 

3(Д

z0)

 

 

Фаза колебаний

ср

определится

из

 

 

tgCf,= b_

а

EjX—3e2Xz0 —{- ——( ^ + l )

 

___________ _________

(4.48)

 

 

£хХ— З е 2Хг0 —

(5 + 1)

 

 

ші

 

§ 1. ГРАНИЦЫ ОБЛАСТИ СУЩЕСТВОВАНИЯ «ГЛАВНОГО»

РЕШЕНИЯ

Исходя -из того, что A2 — a2jrb2— существенно положи­ тельная величина, можем заключить из (4.43), что парамет­ рический эффект возможен только в системе с несимметрич­

ной характеристикой (ерфЩ), так как если еі = 0, то

z0=0 и

значение кормя — мнимое. Далее, из (4.43) видно, что

ампли­

туда X играет в явлении резонанса двойную роль: она стиму­ лирует возникновение параметрического эффекта (множитель в выражении под корнем) и 'одновременно как бы гасит его, а, точнее, вызывает сдвиг резонансной кривой вправо, в сто­ рону больших расстроек 5 (множитель при е2). Стимули-

233

■рующее влияние Л, как легко видеть, пропорционально пер­ вой степени, а сдвигающее — квадрату X. По этой причине сдвигающий эффект при малых X выражен слабо и затем быстро возрастает три увеличении амплитуды внешнего воз­ действия. Главное решение А2> 0 возможно лишь при ус­ ловии

(£іХ—Зв,Хг0)2 >

(5+1).

(4.49)

Для частоты возбуждения 2ш0ф é0, условие (4.49) принимает вид (при £ = £0—1)

Хгаіп = -----

25У1е»

(4.50)

Если условие (4.49) выполнено, то существует конечная об­ ласть расстройки S, внутри шторой 'имеет место параметри­ ческое возбуждение. Границы этой области определяются из ус­ ловия: Л2 = 0. Тогда .из (4.43) имеем

3

с-{-1—

®2^" “Ь ^£і^о — Зе325 Нд

 

±

X2 — 4 4 (? + 1) = 0 .

(4 -5 1 )

Г

й>5

 

Решая это уравнение относительно I, получим предельные границы области существования главного решения; разность определяет относительную величину частотного ин­ тервала, внутри которого параметрический эффект возмо­

жен. В безразмерной форме эти границы имеют :вид

_

 

о

__

ара

 

£і= 3(Д—z0)2 +

— Яа— ———

 

 

 

2

 

 

 

р / з ( Д - г~Г‘ ( з Р - * £ )

- 4

J

( 1 + 4 ^ 5 ) ,

(4.62)

— —

9

 

9п2

 

Ь = 3 (А - г0)2 + 4

X2-

f i- +

 

 

2

 

 

 

 

+ у Г

 

 

Ȁ

+

(4.63)

 

 

 

 

 

234

о 5х=2 ] / 3(Д

“o' / (4.54)

Величина г0, входящая в вышеприведенные формулы, долж­ на удовлетворять уравнению (4.44), если положить в пос­ леднем Л2= 0; в безразмерном виде это уравнение имеет вид

3(Д—z0)X-= 2z„ [ 1 +Zq+ ЗД (Д—z0) ].

(4.55)

§ 2. ИССЛЕДОВАНИЕ УСТОЙЧИВОСТИ ПЕРИОДИЧЕСКИХ РЕШЕНИЙ

Для того чтобы «айденные периодические решения могли реально 'существовать, необходимо, чтобы они были устой­ чивыми. Устойчивость будем понимать в смысле Ляпунова, т. е. состояние системы будем считать устойчивым, если под действием достаточно малого возмущения ее координаты не изменяются заметным образом, .находясь в окрестности неівозмущенного состояния. В случае стремления координат системы возвратиться в исходное положение, существовавшее к моменту приложения возмущения, будем говорить об их асимптотической устойчивости. Как известно, необходимым признаком устойчивости является наличие отрицательных ве­ щественных частей всех корней характеристического урав­

нения.

Исследуем устойчивость решений, удовлетворяющих

уравнениям (4.41) и

(4.42). Пусть существуют действитель­

ные положительные

решения а0, Ь0> г0 системы (4.41) и

4.42). Чтобы исследовать их устойчивость, рассмотрим пове­ дение системы при слепка возмущенных значениях амплитуд:

а = а о+Р> Ь — Ь0-\~ѵ].

(4.56)

іПотставим (4.56) в (4.37) и (4.38) и разложим полученные функции в ряды по малым возмущениям; приняв во внима­

ние, что а0, Ь0, 20 являются решениями невозмущеиной сис­ темы, получим

(4.57)

235

Будем предполагать, что в период установления колебаний, а также в период движения системы под действием возмуще­

ния, связь

между

амплитуда ми субгармоники А,

накачки X

и смещением г0 такая

же, как и

в периодическом решении,

т. е. подчиняется

уравнению (4.44). В системе

(4.57)

для

іеокращения

обозначено

 

 

 

 

 

/<ЭФ\

__ /б?Ф\

 

ідФ\

_

/дФ\

 

 

\ d

ß

j 0

U

a j ^

=ob )00

ІД)a=Ь aa""'

 

 

 

 

 

b - b Q

 

 

ь = Ь й

 

 

а точкам

соответствуют нелинейные

члены разложения

по

степеням возмущений.

 

 

 

 

 

 

Уравнения в вариациях

(4.57)

определяют закон изме­

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

ределяется линейными

членами

dp

 

It

(4.58)

dt]

 

I t

 

Уравнения (4.58) будут удовлетворены, если решения их при­ нять в виде

р = р0еа/; 1? = р 0 еа/.

(4.59)

Подстановка (4.59) в (4.58) дает характеристическое уравнение

(4.60)

/гіф

— а

т?=0.

 

/Для совместности этих уравнений необходимо, чтобы

' £ )

10

( f

)

 

\да

\âb

10

(4-61)

 

 

 

= 0.

/д± ) ,

Ж _ «

\да ! 0

\дЬ 10

(Решая уравнение (4.61), найдем

23S

1

а =.—

2

(4.62)

Чтобы возмущения р и у со временем затухали, необходимо, чтобы оба значения а были отрицательными или имели отри­ цательную вещественную часть. Из (4.62) видно, что для ■этого должны быть выполнены два условия

(4.63)

 

 

 

 

 

 

 

(4.64)

Из

уравнений

(4.37) и

(4.38) имеем

 

 

 

(дФ\

= _

1

3e,Xz0+

4- — e.flb

 

\да ) 0

 

 

2

 

 

 

“ о

2

 

 

(дФ'

 

 

— (§+1) + S0+ — S2^2~ 2e1Z0 +

 

 

 

 

 

 

 

+

— e2/424-3e2Zo-!- — e2b2

(4.65)

 

 

 

 

4

2

 

 

l d± ]

= _ i

— (£ + 1 ) + eo + —

2eizo +

 

\d a )0

2

 

 

 

 

 

 

 

4- — е2Л2+ Зг22о+ — s2a2

 

 

 

 

 

4

2

 

 

 

(d±

 

1

EjX-

25t!)!

 

 

 

\db

 

2

-3c2Xz0

 

 

 

 

 

 

 

 

Если оба значения Rea, в (4.62) отрицательны,

то исследуе­

мое

движение

устойчиво. Первое

условие

устойчивости

(4.63) с учетом

(4.65)

дает

 

 

237

 

1 <'

SjA — öZoA^q-p

2 S w l

+

 

1 о 1 — I

 

 

 

 

CO"

 

+ — (еА ~Зеглг0-

iOÜ),

2otO!

< 0 ,

(4.66)

 

tos

CO«

 

 

что выполняется, если ö>0; это условие соответствует тре­ бованию, чтобы в системе не было самовозбуждения коле­ баний (автоколебаний). Из второго условия (4.64) получаем

'/дФ\ / < М Л _

/дФ\ Л ? Ф \ _ _

1

EjA— 3e.Xz0 +

2

+

W / o W / o

l ^ / o W / o

2

 

 

 

Ü>n

3

\

2

 

2o(Oi

3

,

-f-

+ -у

ЧаЪ ) X — UX—3e2Xz0 —

----- -f- — B„cib I

 

 

 

 

ton

2

 

 

 

— ( 5 + 1) + eo ~b — - 3^2~ 2 £120-1- —

е2Л3 + Зе22^+

e2^2

X

—(ë “h l) + £o

+

V*-2 - 2еіго +

e2^ 2+

 

 

+ 3 e22 „ -f-

£„tr

> 0.

 

 

(4.67)

Рассмотрим прежде всего вопрос об устойчивости началь­ ной формы равновесия, для которой А — а=Ь 0 (отсутст­ вие параметрических колебаний); из этих условий получает­ ся, что равновесие системы устойчиво, если

1.8 > 0 ;

2.

4 -І_ (£-j-l) — (еіА—Зе.Д20)2

+ 4

4

-

(?+1) + ®о +

 

4

. шо

J

 

 

 

 

 

+

2sl'£0~b 3e22J

> 0 .

(4.68)

Неравенство (4.68), как нетрудно видеть, дает для границ области возбуждения уже полученные 'нами результаты (4.52) и (4.53), разумеется, с учётам знака неравенства.

Тактам образом, область самовозбуждения субгармошгкм Ъ.—Іі совпадает с диапазоном (4.68) неустойчивости стацио­ нарных колебаний вынужденного режима. Устойчивым здесь оказывается «главное» решение (4.29), а вынужденные коле238

Соседние файлы в папке книги из ГПНТБ