
книги из ГПНТБ / Бондарь, Н. Г. Нелинейные стационарные колебания
.pdfЗаметим, что при больших коэффициентах затухания естествен но вместо равенства (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 |
+ |
4л2 со2 ] Ца - - р.2со2)2 + |
4па ц2 со2 ] |
~ |
3 |
V |
3 | 0 | |
р Y ( a — со2 )2 |
+ |
4л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 |
+ |
4л2 ц2 со2 + |
12лцсо (Уа — |хш)2| |
(11.95) |
|||
X |
|
+ |
4л2 ш2 + |
12лсо (Уа — со)2 + |
|
|||
р Via. — со2 )2 |
|
|
||||||
+ Via |
— р.2со2)2 |
+ |
4л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 |
|
|
|
|
4а |
|
|
|
ал 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 < а + |
|
4Л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 + |
4л2 ц2 со2 |
|
|
Эта формула получена аналогично уравнению (11.94). |
|
|||
При больших |
значениях |
коэффициентов |
затухания лучшую |
|
сходимость с машинными решениями дает формула, |
аналогичная |
|||
выражению (11.95): |
|
|
|
|
( ± Т - а У ш |
- б * ) х |
У |
[(а* — со2 )2 + |
4л2 со2 + |
12лш ( К а * — ш)2 |
] |
|
|||
X |
[(а* — ц2 со2 )2 |
+ 4л2 ц2 со2 |
+ |
12л|ЛСО {Уа — (ш)2 ] |
(11.103) |
|||
* |
|
|
|
— |
* |
|
||
Р У (а |
— со2 )2 + |
4л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
|
|
4А |
кр Я2 —1 |
|
|
А2 |
|
|
|
|
|
|
|
|
|
+ |
16—?£-п*а>*\ X |
|
|
||||
|
X |
•|Х2 С02 |
4А кр |
Л 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. К оценке влияния взаимо~ |
имодействия гармоник на устой |
действия гармоник на устойчивость |
|
чивость симметричных колебаний |
несимметричных колебаний без тре |
|
при турбулентном |
сопротивле |
ния. |
нии.