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

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

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

Заметим, что при больших коэффициентах затухания естествен­ но вместо равенства (1.250) принять приближенное выражение

62 (1 +

2В)2 + / г 2 « - £ - .

(11.90)

Тогда уравнение (11.89)

критических состояний

упрощается 1

к виду

 

 

X ' У ( Т

 

2

- а

 

2

со2

+4п-шЧ

рАо2 )

+4п^*ш*

 

 

 

 

(11.91)

 

+ 4 / l 2 t 0 i +

( Т " Цаш* + 4л2 ц2 со2

Подставляя в эту формулу со = 0, находим координату точки пересечения кривой / критических состояний с осью F (см. рис. 84):

 

 

 

F ' =

2 ( 1 + р )

V

IPI "

 

 

Чтобы построить приближенную формулу для кривых

/ / крити­

ческих

состояний,

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

изложенной

выше

методикой

и подставим в уравнение (11.41) приближенные

выражения

 

 

О2

да/г2+

02 да а,.

/ « С

= const.

(11.92)

Будем

иметь

 

 

 

 

 

 

 

 

Р

С

У [(а - т'У» +

4п2 со2 ) [(а - ц2 со2 )2 +

4п»рЛо«Т

щ ggj

 

}Лх

р V (а — со2 )2 +

4/г2со2

+ V (а. — |x2coa)a + 4л2 ц2 со2

 

Поскольку

сопротивления

колебаниям не оказывают

влияния

на статические перемещения, точка пересечения кривой / /

критиче­

ских состояний с осью F должна совпадать с точкой Fu в случае от­

сутствия трения (п =

0). Подставляя в уравнение (11.93) граничные

условия со = 0 и F =

FJJ, а также используя формулу (11.85), по­

лучаем выражение (11.87). Подставляя его в равенство (11.93),

окончательно

имеем

 

 

 

р

2

-л Г

а

УЦа — со2 )2

+

2 со2 ] Ца - - р.2со2)2 +

4па ц2 со2 ]

~

3

V

3 | 0 |

р Y ( a — со2 )2

+

2 со2 + V(a — р.2со2)2

+ 4л2р.2со2

 

 

 

 

 

 

 

(11.94)

Полагая здесь я =

0, как и следовало ожидать, приходим к уравне­

нию (11.86).

 

 

 

 

 

Заметим, что уравнение (11.94) кривых 77 критических состоя­ ний, являющихся геометрическим местом точек срыва колебаний, справедливо для малых коэффициентов затухания. При больших коэффициентах затухания лучше соответствует машинным решени-

1 Формулой (II.9Г) следует пользоваться для со < У^сГ. Это ограничение вы­ текает из сравнения аналитических и машинных решений.

140

ям следующее выражение х :

 

 

 

 

3IPI

X

 

 

V[(а — со2 )2 + 4л2 со2

+ 12лсо (У а — со)2] X

 

 

X

[(а — |д.2со2)2

+

2 ц2 со2 +

12лцсо (Уа — |хш)2|

(11.95)

X

 

+

2 ш2 +

12лсо (Уа — со)2 +

 

р Via. — со2 )2

 

 

+ Via

— р.2со2)2

+

2 2 со2

+ 12лцсо (У а — р.со)2

 

 

Рассмотрим также влияние турбулентного сопротивления на

устойчивость стационарных

колебаний,

описываемых

уравнением

х + rixz- • sgn х + ах — | Р 1 я 3

=

 

Fx cos со/ 4 F-г, cos цсо/.

Корректируя

амплитудно-частотную характеристику (11.54) заменой Э на 0 (1 4

4 25) и используя

формулы

(1.243) и (1.250), получаем уравнение

кривой

/ критических

состояний:

 

 

 

 

 

 

 

 

 

 

 

PI X

 

 

 

 

а

/ 4

а

• Л 2 l j ' 4 16

 

2 со4 X

 

 

 

 

 

я 2 | Р |

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

ал 2

 

 

 

 

 

 

2

 

 

У [ - г + - ( -

л 2

 

•л 2 — 1

Н - 1 6 - — ^ ( 0 * 4

 

| Р I

 

 

 

 

,

о

 

 

/

4ал

2

УР

 

ал 2 .

 

+ У | - + ^ 2

( — - ! ) ]

4 1 6 ^ ^

Здесь

использованы обозначения

(11.78).

 

(11.96)

 

 

Кривая / / критических состояний может быть построена по урав­ нению, полученному аналогично уравнению (11.94) для вязкого

трения:

 

 

 

 

 

3IPI

X

а 4 со2

я 3

л3 — 1 +1 6

л 3 со4 X

\

/ J

я 2

X < а +

 

2

А1

ц2 со2

 

X

 

2

да

 

 

 

 

 

+ 16 —— л2 со4 +

+ | / | a + p A o 2

n 2 — l'

+ 1 6 - + - л 2 р . 4 с о 4

 

 

 

2

1 Эмпирические формулы, приведенные в этом параграфе, получены В. С.Гор­ батовым.

141

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

 

з

V з | р

i X

 

(а — со2 )2 +

6

 

• со)3

лсо

(а — ц 2 ш 2 ) 2 +

 

+ 6 у з | р

 

 

 

лр.ш

X

 

 

 

(11.97)

 

 

 

 

р | / " ( а -

со 2 ) 2 + 6 У

 

[/а

-

со)2 лсо +

f | / ~ ( a - p . 2 c u 2 ) 2 - f - 6 У

-г

• CV^a

— М-м)3 лцсо

Для оценки точности полученных результатов на рис. 85 постро­ ены кривые критических состояний для системы с параметрами а =

= 1 сек-2;

| В | = 0,2 см~2

• сект2;

р. = 2; р = 1 при различных

значениях

п по формулам

(11.89),

(11.91), (11.94) и (11.95), а на

рис. 86 — по формулам (11.96) и (11.97). Здесь же приведены резуль­ таты решений на ABM МН-7. Как видно, совпадение аналитических

 

 

Z

асек*

 

/

 

Z а.се К1

Рис.

85.

Кривые критических

Рис. 86. Кривые критических со­

состояний

симметричной

систе­

стояний

симметричной

системы

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

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

моническом возбуждении:

 

и бигармоническом возбуждении:

/ = г л = 0 , 0 5

сек~1; 2—л=0,20

сек—1;

л = 0 , 0 5 см~1;

2—п=0,20

см—1;

В =

л = 0,50 с е к - 1 .

 

3 _ л =

0,40

см~1.

 

142

и машинных результатов можно признать удовлетворительным *•. Порядок пользования графиками рис. 85 и 86 такой же, как и гра­ фиками рис. 71.

Несимметричное возбуждение. Рассмотрим случай отсутствия трения, т. е. когда колебания описываются уравнением (11.57). В § 3 данной главы получена амплитудно-частотная характеристика

для колебаний

с

умеренно

большими

амплитудами. Заменяя 0

на

0 (1 +

25), получаем характеристику

для колебаний,

близких

к

неустойчивости:

 

 

 

 

 

 

ч - у П / а

1

161 Л 2 - -

6*

I

9 ( 1 + 2 B ) F 1

 

 

1

F а *

2

1^1

9(1 + 2В) ^

92 (1+2S)2 - со3

^

 

 

 

 

.

9(1+25)/ ^

 

 

 

^92(1 + 2В)3 — рАо2 -

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

F =

. — со2

2 ш2

X

(11.98)

+

min

Здесь использованы обозначения (1.282) и (11.78). Параметр а"опре­ деляется из уравнения (1.287) и первого равенства (1.285).

Из двух значений (1.328) критической амплитуды Л к р в формулу (II.98) подставляется то, которое дает минимальную по абсолютному значению амплитуду возбуждения F.

Исследование выражения (11.98) показывает, что кривая / критических состояний (рис. 87) пересекается с осями в точках

скоординатами

F,=

± J ^ V 2

a , - \ p - \ A $ p - 2 - A

1 + Р

 

к а

а .

Кривая I I критических состояний от срыва колебаний может быть построена путем совместного решения системы уравнения амплитуд-

1 Следует иметь в виду, что ошибка решений на ABM МН-7 может составить 10—15%.

143

но-частотной характеристики и ее производной по Л, где принято

J®. - О- dA

(/0 К) С02 - й 2 ) (02 — I-*2"2) = ^02 1 0 2 ы 2 1 + 1 0 2 — ^2ю21); (11.99)

(/'0 + /0') (02 — со2) (03 — uAo2) + 200' (/0 — бJ [202 — (1 + ii2 ) со2] =

= 2FQQ' (р 1202 — со21 + 1203 — рАо21).

Решая параметрические уравнения (11.99), аналогично изложен­ ному выше находим координаты пересечения кривой / / с осями:

 

•та3 | Р |

min 1

С0Ч

_

со4

= Ya\-

Построение кривой / / критичес­ ких состояний с использованием па­ раметрических уравнений (П. 99) весьма сложно из-за громоздкости вычислений. Результаты этих расче­ тов хорошо аппроксимируются прос­ тым выражением,'полученным анало­ гично формуле (11.86):

 

 

 

 

 

F

=

а .

 

3|Р1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

о,сек

 

 

•б*

X

 

 

 

 

 

 

 

 

 

 

 

Рис. 87.

Кривые критических со­

 

 

(а — ш2) (а — |х2ш2)

 

стояний

несимметричной

системы

X

 

 

без трения

при

бигармоническом

 

+

К

• и,'-ш2

 

возбуждении.

 

 

 

 

min

-

 

 

 

 

 

 

 

 

(11.100)

Для оценки точности полученных результатов на рис. 87 по фор­

мулам (11.98) и (11.100) построены

графики критических состояний

для системы с параметрами а =

1 се/с- 1 ; |В\ =

0,2 см—2

• сек-2;

у = 0,15 см—1

• сек~2;

б0 F0

=

— 0,5 см • сек—2; р = 1; |л = 2.

Кружками

(кривая

/) и точками

(кривая II)

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

зультаты

решений на

ABM МН-7.

Соответствие

аналитических

и машинных результатов можно признать хорошим. Порядок поль­ зования графиками рис. 87 такой же, как для графиков рис. 84.

Далее рассмотрим устойчивость стационарных колебаний при наличии вязкого трения, т. е. когда колебания описываются урав­ нениями, (11.60) или (11.61). Амплитудно-частотная характеристика в этом случае определяется формулой (11.63). Заменив в ней 0

1 4 4

на 0 (1 + 2В) и используя выражения (1.326) и (11.78), получаем уравнение кривой / критических состояний:

 

F

=

 

 

 

 

а

 

У

1Г/ а

у

X

 

+ л г _Ш 2 1 + 4„2Ш2

+ л2 — |л2ю2J + 4л2ц,2со2

 

 

 

 

 

] /

" +

я 2 -

со 2 ) 2 + 4л2 со2 + ] / ^ ( " ^ - + п " ~ l ^ 2 » 2 ) 2 + 4л*ц»ш»

 

 

 

 

 

(11.101)

Из двух

значений

(1.328)

критической амплитуды Л к р в формулу

(11.101) подставляется то, которое дает минимальную по абсолют­ ному значению амплитуду возбуждения F. При больших коэффи­ циентах затухания и со <; |/"^*следует использовать приближенное

выражение (11.90), т. е. в круглых скобках формулы (11.101) сле­

дует положить п = 0.

 

 

 

Для построения

кривой / /

критических состояний

при умерен­

ных сопротивлениях можно пользоваться выражением

 

F =

( ± 4 а * / з Щ - - б * ) Х

 

 

 

У [[а] — со2 )2 -|- 4л2 ш2 ] [(а* — р.2со2)2 + 4л2 р.2 со2 ]

(11.102)

X

 

 

 

р ] / ( а * — со2 )2 + 4л2 со2 +

]/"(а] — ц2 со2 )2 +

2 ц2 со2

 

Эта формула получена аналогично уравнению (11.94).

 

При больших

значениях

коэффициентов

затухания лучшую

сходимость с машинными решениями дает формула,

аналогичная

выражению (11.95):

 

 

 

 

( ± Т - а У ш

- б * ) х

У

[(а* — со2 )2 +

2 со2 +

12лш ( К а * — ш)2

]

 

X

[(а* — ц2 со2 )2

+ 4л2 ц2 со2

+

12л|ЛСО {Уа — (ш)2 ]

(11.103)

*

 

 

 

*

 

Р У (а

со2 )2 +

2 со2 +

12лш (У а

— со)2

+

 

+ У(а'

— ц.2 ша )

г + 4л2 ц2 ш2 +

12лцсо ( ] / а * — jico)2

min

 

 

 

 

 

 

 

 

10 4-5

145

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

на рис. 88

по формулам (11.101) — (11.103) построены графики

критических

состояний для системы с теми же параметрами, что и на рис. 87 при различных значениях п. Здесь же приведены результаты реше­ ния на АВМ МН-7. Как видим, соответствие аналитических и машинных результатов можно признать удовлетворительным.

Рассмотрим устойчивость ста­ ционарных колебаний при тур­ булентном сопротивлении, т. е. когда колебания описываются уравнением (11.66) или (11.67). Амплитудно-частотная характе­ ристика для этого случая опре­ деляется выражениями (11.70) и (11.72). Заменив в них 0 на

Рис. 88. Кривые критических состоя­ ний несимметричной системы при вяз­ ком трении и бигармоническом возбуж­ дении:

/ ~ п

=

0,05

с е / с - 1 ; 2 — л = 0,20 сек~

S — п

•=

0,50

« к - 1 .

0 (1 + 2В) и используя формулы (1.326) и (11.78), получаем урав­ нение кривой / критических состояний:

F = ± _ ^ ] / 2 a , - | P | 4 p - 2 J ± \ X

 

 

кр Я2 1

 

 

А2

 

 

 

 

 

 

 

 

+

16—?£-п*а>*\ X

 

 

 

X

•|Х2 С02

кр

Л 2

1

+ 16 — i £ _ n iu«o> 4

X

 

 

 

 

 

 

 

я 2

 

 

 

У14-

 

 

 

 

 

 

 

 

 

 

 

 

 

л 2 — 1

 

 

я 2

ш4

+

 

 

 

+ С02!

4 А К Р

 

+ 16—!2_ л2

 

 

+

-|12 С02

4А:кр

Л 2

1

+16 -

- р - я 2

ц 4

с о

4

 

 

Я2

 

 

 

 

я 2

 

 

|min

(11.104)

146

Из двух значений (1.328) критической

амплитуды Л к р

в

формулу

(11.104) подставляется то, которое дает минимальную по абсолют­

ному значению возбуждения

 

F.

 

 

 

 

 

 

Для построения кривой / / критических состояний можно пользо­

ваться выражением, аналогичным уравнению (11.97):

 

 

F =

( ± 4 A * / 3 ¥ T ~ 6

* ) x

 

 

 

 

 

 

] / к а * — со2 )2 +

6Л,пш(]/а* — со)3] [(а* — ц2 ш2 )2 +

б/^лцш)Уа] — цш)3]

 

р 1/"(а* - со2 )2

+ 6Л*лй) (]/"«* — со)3 +

 

 

+

~

И2 ш2 )2

+

6Л»пцсо ( ] / ^ _

1^)3

 

 

 

 

 

 

 

 

 

 

 

 

(11.105)

где Л.,. абсолютная

величина

 

координаты

экстремальной точки

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

системы,

 

 

 

 

 

 

 

 

 

А*. —

1

• ( ?

 

± ^ . 2

+ За;|В|)

i n

(П.106)

 

3 | В |

V I"*

'

Г. I

1 Г I/ | m

Входящие сюда параметры определяются из выражений (1.285) — (1.287).

Для оценки точности полученных результатов на рис. 89 по фор­ мулам (11.104) и (11.105) построены графики критических состоя­

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

1 сек—1

; | В \ = 0,2 см~2

х

% сек~2; у =

0,15 см~1 • сек—2; 80

— F0 =

0,5

см • сект2; р

=

= I," [i = 2.

Здесь же представлены

результаты

решений на АВМ

МН-7. Как видно, совпадение результатов можно признать удовлет­ ворительным.

Взаимодействие гармоник. Оценим влияние взаимодействия гармоник на устойчивость стационарных колебаний. На рис. 90— 95 выполнено сопоставление графиков критических состояний при бигармоническом возбуждении (сплошные линии) и при моно­ гармоническом возбуждении (штрих-пунктирные линии). На этих же рисунках приведены результаты решения на АВМ МН-7. Причем моногармоническому возбуждению соответствуют кружки (кривая Г)

и

точки (кривая II), а бигармоническому — крестики (кривая II)

и штрихи в кружках (кривая I). Рис. 90—92 соответствуют симмет­

ричным системам, а рис. 93—95 — несимметричным.

Рис. 90 представляет собой совмещение рис. 59 и 84. Рис. 91

построен с использованием данных рис. 62 и 85. Рис.

92 построен

с использованием данных рис. 64 и 86 и соответствует

случаю тур­

булентного сопротивления с коэффициентом п =

0,2

см—1 . Рис. 93

представляет собой

совмещение рис. 71 и 87 и

отражает случай

отсутствия трения в

несимметричной системе. Рис.

94 построен с

ю*

 

 

147

Рис. 89.

Графики

критических

состоя­

Рис. 90.

К оценке

влияния взаимо­

ний несимметричной системы при тур­

действия

гармоник

на устойчивость

булентном сопротивлении, и бигармони­

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

ческом возбуждении:'

 

 

 

ния. :

 

 

 

/ _ п =

0.05 с м ~

2 —

п =

0.20

см.—1.

 

 

 

 

 

сисе/г''

 

 

 

i 1

 

 

У

 

 

 

 

 

 

1

 

 

 

 

 

 

i

1

 

 

11

 

t,

 

 

1

 

 

 

 

//.

 

 

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

 

 

 

 

 

 

 

1

1/f

I/

 

 

 

 

 

 

 

1

 

 

 

/

 

 

 

 

 

 

 

 

 

 

 

4

2 o,ce/f'

Рис. 91. К оценке влияния взаимодействия гармоник на устой­ чивость симметричных колебаний при вязком трении:

•а — п = 0,2 сек—1; б — п = 0,5 сек—1.

Рис. 92. К оценке

влияния вза­

Рис. 93. К оценке влияния взаимо~

имодействия гармоник на устой­

действия гармоник на устойчивость

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

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

при турбулентном

сопротивле­

ния.

нии.

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