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

книги из ГПНТБ / Бондарь, Н. Г. Нелинейные стационарные колебания

.pdf
Скачиваний:
30
Добавлен:
20.10.2023
Размер:
7.67 Mб
Скачать

характеристике становится вертикальной. Следовательно, Для1 'оп­ ределения условий перескока колебаний необходимо продифферен­

цировать уравнение амплитудной кривой по Л и положить

= 0.

Поскольку точки срыва колебаний находятся вдали от критической амплитуды колебаний, то можно воспользоваться уравнением ам­ плитудно-частотной характеристики в виде (1.23), т. е.

 

f(A)[Q2(A) — со2] =FB(A).

(1.255)

Дифференцируя

выражение (1.255) по Л и полагая

= 0, полу­

чаем 1

/'(82 — со2 )+ 2/0 • 9' = F9\

(1.256)

 

Выражения

(1.255) и (1.256) представляют собой уравнения кри­

вой / / критического состояния в параметрической форме. Вычисле­ ния по ним производятся в следующем порядке. Для мягкой систе­ мы с заданными параметрами задаются амплитудой Л и, вычисляют / (Л), 8 (Л), f (А) и 9' (Л). Затем решают систему уравнений (1.255) и (1.256) и находят со и F. Если эти величины вещественны и поло­ жительны, то на графике наносится точка с координатами со и F. Следует заметить, что эти вычисления весьма громоздки. Оказывает­ ся, что результаты могут быть аппроксимированы простым выраже­ нием |

Это выражение получено из следующих соображений. Принимаем линейное приближение для частоты свободных колебаний 92 = а. Тогда£0' = 0 и уравнение (1.256) принимает вид f (a — со'2) = 0. Поскольку a —1 со2 Ф 0, то f = 0. Отсюда f = с = const, т. е. ли­ нейное приближение для частоты влечет за собой осреднение ампли­ тудной функции. Теперь уравнение (1.255) принимает вид

c(a — a2) = VaF.

(1.258)

Ниже будет показано, что кривая

/ /

пересекается

с осью орди­

нат в точке

 

 

 

со = 0; F - F u

-

^ y ^ .

d-259)

Используя эти граничные условия, находим

—§-«Vw- . (L260)

Подставляя выражение (1.260) в (1.258), получаем уравнение (1.257). На рис. 59 кривая / / построена по формуле (1.258) для тех же параметров, что и кривая /. Точками и кружками показаны

1 Штрихами обозначены производные по А.

80

результаты решения 1 на ABM МН-7. Как видим,

совпадение

ана­

литического и машинного решений удовлетворительно 2 .

 

 

 

Заметим, что кривые / и / / н а рис. 59 пересекаются. Это проис­

ходит потому, что кривая

/, как видно из выражения

(1.251), при

F

= 0 и со =

0 пересекается

с горизонтальной осью в

точке

со =

=

"j/' - y - и с

вертикальной

осью

в точке Ft =

|/"-гуг- •

Кри­

вая / / пересекается с

горизонтальной

осью в точке со = ] / а .

Этот

результат вытекает из следующих

соображений.

Полагая

в

урав­

нениях

(1.255) и (1.256)

F = 0,

получаем

 

 

 

 

 

9 = со;

2/09'=0;

299' = 0; _dQ2/dA = 0;

02 = const.

 

Поскольку для малых F при со =

У а (см. рис. 59) будут иметь место

малые линейные колебания, то 03

= а.

 

 

 

 

 

 

Определим точку пересечения

кривой

/ / с вертикальной

осью.

Полагая в формулах (1.255) и (1.256) со =

0, находим

 

 

 

 

 

 

 

/9 =

/?//;

/'92

+

2/00' =

FnQ'.

 

 

(1.261)

Исключая отсюда Fu,

получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/'0 +

/0' =

0.

 

 

 

(1.262)

 

При статическом воздействии

на

систему (со = 0,

х =

0,

х =

=

Л„)

уравнение (1.192)

принимает

вид

 

 

 

 

 

 

 

 

 

 

« A 0 - | B | A o W / 7 .

 

 

(1.263)

Дифференцируя это равенство по Л0 , имеем

 

 

 

 

 

 

 

 

 

а 3 |Р \А\ = F'u.

 

 

(1.264)

В то же время, дифференцируя

по А0 первое выражение

(1.261),

находим f'0

-f- /8' = F/j.

Сопоставляя это равенство с формулой

(1.262), определяем F'u =

0. Далее, из выражения (1.263) получаем

А0

=

" j / " 3 "p j

. Подставляя это значение в равенство (1.263), по­

лучаем формулу

(1.259).

 

 

 

 

 

 

 

 

 

 

 

На рис. 59 точки, лежащие выше сплошных линий,

характери­

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

которых реализуются

нерезонансные

колебания

(см. § 2 данной

главы). Штриховая часть кривой / / характеризует

параметры

воз­

буждения, при которых происходит перескок с нерезонансных

коле­

баний на устойчивые

резонансные

(см. кривые

/ / / на рис. 60).

1 Решения на АВМ, приведенные в этом параграфе, получены В. С. Горбато­

вым. Для кривой / / здесь и далее опущены результаты, полученные вблизи

резо-

нансов субгармонических колебаний, т. е. при со = 6/3, 6/5...

 

2 Следует иметь в виду, что погрешность вычислений на АВМ МН-7 может достигать 10—15%. Причем особенность решения задач устойчивости на АВМ МН-7 состоит в том, что результат для кривой / всегда занижен.

6 4-5

81

Штриховая часть левой ветви кривой / характеризует параметры возбуждения, при которых резонансные стационарные колебания будут неустойчивыми.

Критерии

неустойчивости при вязком

трении. Подставляя вы­

ражения (1.243) и (1.250)

в формулу. (1.248), получаем уравнение

 

 

критического

состояния при наличии

^

/

вязкого трения

 

] / " | п 2 — со2 + - | - j 2 + 4л2 соа

(1.265)

Поскольку (/г2— со2 + a/2)a +4rt2 coV= Ф 0, то, как видно из уравнения (1.265), кривая / критических состояний не пересекает горизонтальную ось (рис. 61).

Рис. 61. К построению графика

Исследуем

выражение

(1.265) на

критических состояний при вяз­

минимум. Дифференцируя его по со

ком трении.

 

dF

п

 

 

 

 

 

и приравнивая

=

0, находим ко­

ординату минимального значения со0 = Vra/2 п2.

Подставляя

это выражение в формулу (1.265), получаем минимальное значе­ ние критической амплитуды возбуждения (см. рис. 61):

2

Заметим, что приближенное равенство (1.250) справедливо для случая отсутствия трения. Поэтому формулой (1.265) допустимо пользоваться для малых сопротивлений. При больших сопротивле­ ниях естественнее воспользоваться приближенным равенством п2 + + 92 (1 + 2В)2 = а/2. Подставляя это выражение, а также равен­ ство (1.243) в формулу (1.248), получаем следующее уравнение

критического состояния

F

, со< / а

(1.266)

Перейдем к рассмотрению кривой / / графика критических со­ стояний. Как показано выше, кривая / / представляет собой геомет­ рическое место точек срыва колебаний. Эти точки, естественно, лежат на амплитудно-частотной характеристике для умеренных амплитуд, которая описывается уравнением (1.99), т. е.

f У{пг

— со2 + 02 )2 + 4я2со2 = QF.

(1.267)

1 Ограничение со < У

вытекает из сопоставления аналитических и машин­

ных решений (см. ниже).

82

Далее поступаем аналогично изложенному выше, а именно: принимаем линейное приближение для частоты свободных колеба­ ний и осредняем амплитудную функцию, т. е. 02 л? п2 + 02 => а; f да с = const. Подставляя эти приближения в равенство (1.267), получаем

F = -у=- У (а — со2)2 + 4/г2со2.

(1.268)

Постоянную с найдем из следующих соображений. Точка пере­ сечения кривой / / (рис. 61) с вертикальной осью F по-прежнему определяется координатами (1.260), так как при статическом воз­ действии (со = 0) сопротивления не оказывают влияния на переме­ щения. Поэтому, подставляя в равенство (1.268) граничные условия (1.260), находим

Теперь выражение (1.268) принимает вид

з |ТГ [ ( а — 0 ) 2 ) 2 + 4 " , ( D ' 1 '

( 1 , 2 6 9 )

Полученная формула описывает кривые / / критических состоя­ ний. Поскольку при ее выводе использовались выражения, справед­ ливые при малом трении, этому условию соответствует и формула (1.269).

Для больших сопротивлений лучше соответствует машинным ре­ шениям следующее выражение х :

F = X У ~зЩ1{а ~ 0)2)2 + 1 2 т о ° ^ ~Сй)2]-

( 1 - 2 7 0 >

Как видно из формулы (1.265), кривая / критических состояний

(см. рис. 61) пересекает вертикальную ось в точке

 

F> = (тг+п2)Уж > -fа У ШТ

= Р ц -

Подставляя значение со0 в формулу (1.269), имеем

 

f — 4 - V r

w [ - T - + ^ a - ^ ] -

При условии Fml„ < Fa

кривые I к II критических состояний

имеют две точки пересечения и образуют область

(на рис. 61 она

заштрихована), характеризующую такие параметры возбуждения,

при которых устойчивость стационарных колебаний

(аналогич­

но тому, как это было при п = 0 на рис. 59) зависит от

начальных

условий.

 

Для оценки точности полученных результатов было проведено

исследование на ABM МН-7 устойчивости стационарных

колебаний

1 Эмпирические формулы, приведенные в этом параграфе,

подобраны

В. С. Горбатовым.

 

6*

83

системы с параметрами а =

1 сект2

и | Р | =

0,2 см.—2 • сект2 для раз­

личных

значений я. Результаты

исследования

представлены на

рис. 26

условными обозначениями. Здесь же показаны аналитичес­

кие решения, полученные

по формулам

(1.265),

(1.266), (1.269)

Рис.

62. Кривые

критических

Рис. 63. Кривые критических состояний

состояний

симметричных коле-

при вязком трении в безразмерных коорди-

баний

при вязком

трении:

 

натах.

I

п

= 0,05

сек—1:

 

2 — л

=

 

=

0,2

сек—1;

3 -

л = 0 , 3

сек~2:

4—

 

а

- 0,5 сек-1:

5

— л = 0 , 7

сек—1.

 

 

и (1.270). Как видно, совпадение результатов оказывается удовле­ творительным.

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

В этом случае уравнения критического состояния (1.265), (1.266), (1.269) и (1.270) принимают следующий вид:

для кривых / при любых значениях п и со ;> у~~а

при больших значениях п и со < у~а

F* = К ( 1 - с о 2 ) 2 + 4«2со2; для кривых / / при малых значениях п

^ - 4 У х [ ( 1 - - Н ' + я И .

84

при больших значениях п

/ 7 * = 4 - / т - [ ( 1 - т - ( й ! ) 2 + 6 ^ ( 1

1

\ 2

У 2 *

 

Построенные по этим формулам кривые критических состояний приведены на рис. 63.

Критерии неустойчивости при турбулентном сопротивлении.

Подставляя выражения (1.243) и (1.250) в (1.249), получаем уравне­ ние критического состояния при наличии турбулентного сопротив­ ления:

 

 

 

 

 

(1.272)

а 2

( .4 ," , • я 2 l) + — ? + 16 — а

/ .

I я 2 1 Р I

/

2 J т

я 2

I P I

Исследуем выражение (1.272) на минимум. Дифференцируя его

по со и приравнивая

= 0, находим координаты минимального

значения кривой /

(см. рис. 61):

 

 

 

 

 

 

 

I P I

1 + Л 2 '

где

 

 

 

 

 

 

 

2 / а

 

(1.273)

 

 

 

 

 

 

*

я / |

Р |

 

 

Полагая в формуле (1.272) со = 0, получаем Fj — а/2 У а/\ В |. Перейдем к рассмотрению кривой ТУ критических состояний.

Используем уравнение (1.149) для амплитудно-частотной характе­ ристики при умеренных амплитудах, которую запишем так:

Принимая линейное приближение для частоты свободных коле­ баний (Э да "j/a) и оередняя амплитудную функцию (/ да с), имеем

2 + 16

со* = F у а. (1.274)

Постоянная с, найденная с использованием начальных условий (1.260), оказалась такой же, как и для случая вязкого трения. Под­ ставляя значение с в равенство (1.274), получаем уравнение критиче­ ских состояний в виде

+ 1 6 — 5 - СО

(1.275)

Этим уравнением неудобно пользоваться, так как оно включает в се­ бя амплитуду колебаний А. Вычисления показывают, что уравнение

85

(1.275) хорошо аппроксимируется

выражением

 

 

 

276)

На рис. 64 представлено сопоставление аналитических

данных

(кривые /—3) с полученными на ABM МИ-7 для системы с парамет-

6"

п-LO-

1

 

0.6-

 

 

 

 

1/

 

 

 

г

//////

 

 

Z сз.сат'

 

 

 

 

 

Рис. 64. Кривые критических со-

Рис. 65.

Кривые критических состояний при

стояний симметричных колебаний

турбулентном сопротивлении в безразмерных

при турбулентном сопротивлении:

координатах.

/

л = 0 , 0 5 см~2—л=0,2

см~

 

 

3

— л •= 0,4 см~

 

 

 

рами а = 1 сект-1

и | В | =

0,2 смг2

• сек - 2 при различных значе­

ниях л. Как видно из рис. 64, совпадение результатов можно при­

знать удовлетворительным.

Практически удобнее пользоваться графиками критических со­ стояний, представленными в безразмерных координатах (1.271), за исключением безразмерного параметра сопротивления, который следует определять по формуле (1.273). В этом случае уравнения критического состояния (1.272) и (1.276) принимают следующий вид:

для кривых /

для кривых / /

со2 ^

86

Построенные по этим формулам кривые критических состояний при­ ведены на рис. 65. Порядок пользования графиками такой же, как и для случая вязкого трения (см. рис. 63). Необходимо отметить, что выполненные в этом параграфе исследования устойчивости стационарных колебаний уравнения Дуффинга тривиально обоб­ щаются на случай любой мягкой системы путем замены истинной характеристики кубической. Параметры а и В определяются из условий равенства характеристик и их первых производных для ненулевого корня истинной характеристики.

§ 7. Несимметричное возбуждение

Рассмотрим часто встречающийся в приложениях слу­ чай несимметричного возбуждения, когда возмущение имеет постоян­ ную составляющую, т. е. F (f) = F0 + F cos at.

Одно из принципиальных различий линейных и нелинейных систем состоит в том, что загружение линейной системы постоянной силой не влияет на частоту свободных колебаний, а загружение не-

Рис. 66. Определение частоты свободных колебаний загруженных систем:

а — линейная система: б — нелинейная система.

линейной системы может существенно изменить частоту свободных колебаний (рис. 66). Как известно, частота малых колебаний опре­ деляется квадратным корнем из тангенса угла наклона касательной к характеристике в нулевой точке. Поэтому для линейной системы 0Х = у tg г|) = 02 , а для нелинейной 0, = ] A g ^ i ф V^tg tpa = 02 .

Таким образом, загружение нелинейной симметричной системы постоянной силой переводит ее в несимметричную систему. Следо­ вательно, несимметричное возбуждение симметричной системы экви­ валентно симметричному возбуждению (F = cos at) несимметрич­ ной системы.

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

87

Частоты свободных колебаний загруженных систем. Определим частоту свободных колебаний нелинейной системы, описываемых уравнением

 

x+R(x)

 

= F0,

R (х) =

S0 +

ax + ух2 + бх3 .

(1.277)

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

(1.52), т. е.

 

 

 

получаем

 

 

х = у — х0,

х0

=

const,

 

 

(1.278)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

у + (а -

2ух0+

3*§Р) у + (у-

ЗрХ) у2

+ рУ = F0

+ ах0

- ух* +

 

 

 

 

 

 

+ Р * § - 6 0 .

 

 

 

 

А-279)

Если принять

зависимость (1.53)

 

 

 

 

 

 

 

 

 

 

 

 

А-О = ^ Г .

 

 

 

 

(1-280)

то уравнение (1.277)

преобразуется

к

виду

 

 

 

где

 

 

 

 

# + °У/ + рУ =

6„

 

 

(1.281)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а * = а ~

ж

;

8*=F°-8°+w(a~~9~)-

 

 

 

( L 2 8 2 )

Если

принять

 

зависимость

 

 

 

 

 

 

 

 

 

 

F 0

+

a r 0 - Y * g +

p*g--60 = 0,

 

(1.283)

то уравнение (1.277) преобразуется

к виду

 

 

 

Здесь

 

 

У* + <У* + Ч*У1 + Ш =

0.

 

(1.284)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а: =

а - 2 у * 0

+

3*2р\-

y* = y-2>$xQ,

(1.285)

где х0 — корни уравнения (1.283). Подстановка

 

 

 

 

 

 

 

 

*о = *о + ^ Г

 

 

 

 

< L 2 8 6 )

приводит уравнение (1.283) к виду

 

 

 

 

 

 

 

 

 

 

«)3

+ - ^ 4 + -j- = °-

 

 

a - 2 8 7 )

В соответствии с зависимостью между

корнями

алгебраических

уравнений и их коэффициентами

корни

уравнения (1.287) будут

равны по величине и противоположны по

знаку

корням

характе­

ристики

уравнения

(1.281). Анализ

характеристики

уравнения

(1.281) приводит к следующим

результатам.

 

 

 

При условиях 1

 

 

 

 

 

 

 

 

 

 

 

 

1 Вторые условия получены ранее [13] и приведены в § 1 данной главы.

88

имеет место один вещественный корень и, следовательно, характе­ ристика является жесткой (см. рис. 9, а).

При условиях

- ? - + w = 0

и л и с с : > 0 - Р > 0 ' 4 « # = tf < 1 2 8 9 )

имеют место два вещественных корня и, следовательно, характери­ стика является полужесткой (см. рис. 9, б).

При условии

имеют место три вещественных корня. Характеристика будет мягкой (см. рис. 9, в), если а* > 0 и р < 0, т. е. когда наибольший и наи­ меньший ненулевые корни имеют разные знаки, или полумягкой

(система с перескоком, см. рис. 9, г), если а* > О, В >

0, 4а*В

<

<

у2, Т. е. когда максимальный и минимальный корни имеют оди­

наковые знаки.

 

 

 

 

 

 

 

Из проведенного анализа вытекает, что в равенства (1.285) под­

ставляется только тот корень уравнения (1.287), для которого

>•

>

0. Заметим, что для частного

случая

симметричных

систем, т. е.

когда 8„ = 0 и у = 0 или у^. =

0 имеют место равенства

 

 

 

а; =

=

а;

6* =

^ .

(1.291)

 

Установим зависимость

между

параметрами уравнений (1.281)

и (1.284). Из сопоставления выражений (1.277) и (1.283) по формулам (1.278), (1.280) и (1.282) получаем

v* = y—w>

a * = « : - | r ; 8* = / ? о + ж ( а *

~ " Ж " )

 

 

(1.292)

Выше (см. § 1) приведены формулы (1.60) и (1.61) для определе­

ния частот свободных колебаний в случае а„. > 0.

 

В литературе [40] для уравнения (1.284) приводится

следующая

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

баний

в случае а* > 0:

 

 

 

 

 

 

 

2

 

 

 

0 = "|/с?

i + 1-JL.

! ! ^ _ ) л 2

(1.293)

 

 

1

12(a) 2

'

 

где Ах

= у% (0) — максимальная амплитуда

колебаний,

связанная

с амплитудой А — у (0) зависимостью, вытекающей из первой фор­ мулы (1.292),

Л 1 = л - Ж *

( L 2 9 4 )

Для сопоставления точности формул (1.60) и (1.293) и определе­ ния пределов их применения необходимо иметь точную формулу.

89

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