
книги из ГПНТБ / Бондарь, Н. Г. Нелинейные стационарные колебания
.pdfколебаний, описываемых уравнением |
(1.54): |
|
|
|||
|
|
cos at |
— В |
cos (ш — 20) / |
+ |
|
, |
cos (со + 29) / |
•со2 |
• (со — 28)2 |
|||
} + A - ( l + - | - S c o s 2 0 / ) . |
(1.319) |
|||||
~г |
е а — (со+29)2 |
Сопоставляя решения (1.231) и (1.319), устанавливаем, что в не симметричных системах при несимметричном возбуждении возника ют суб- и ультрагармонические колебания такие же, как и в сим метричных системах. Следовательно, для определения амплитуд этих колебаний можно полностью воспользоваться результатами § 5 на стоящей главы, только 0 и В следует определять по формулам данного параграфа. Наличие несимметрии приводит к появлению ко лебаний с частотой 20, амплитуды которых, как видно из формулы (1.319), определяются уравнением
Л е Y a * + ~Т Hi e = |
. |
(1.320) |
Вычисления показывают, что имеет вместо неравенство сх^ ^> Это позволяет упростить (1.320) к виду
Л 2 0 |
д а - | | £ ^ . |
|
30 у а . |
1 ал?
-^-РЛ
(1.321)
Входящие сюда величины 0 и В определяются по приведенным в данном параграфе формулам в зависимости от амплитуды колеба ний с частотой со возбуждения. Обобщая эти результаты на случай колебаний при наличии вязкого трения, для решения уравнения (1.107) получаем выражение
У | / « * + 4 " № = -|р + a cos (at — р) + ас cos [(со — 20) t — рс ] —
— а у cos [(со -+- 20) г? —Р у ] + а2о cos (20с" — р2е),
где амплитуды колебаний и сдвиги фаз определяются по формулам, приведенным в § 5 данной главы. Кроме того,
<z28 = |
г |
206,5 |
\ |
|
Р29 = |
, |
4я0 |
|
' |
arctg |
л * - 3 0 * • |
||||
|
W - 3 e » ) * + 1 6 i i * e » |
|
V M |
S |
Сопоставляя полученные данные с результатами § 5 данной главы, видим, что в несимметричных системах при вязком тре нии возникают суб- и ультрагармонические колебания, амплитуды которых можно определять по формулам § 5. Кроме того, возникают колебания с частотой 20, амплитуды которых определяются из приближенного уравнения, полученного аналогично (1.321):
Л 2 9 |
да * |
= - |
2 6 6 * 5 |
. |
(1.322) |
|
Уа* |
|
/ а Л ( л 2 - 3 0 2 ) а + |
|
V |
Полагая в формуле (1.322) п = 0, как и следовало ожидать, при ходим к выражению (1.321).
100
В случае несимметричного возбуждения симметричных систем в приведенных выше формулах следует использовать равенства (1.291).
Более подробно1 рассмотрим стационарные колебания при турбу лентном сопротивлении. Используя замены (1.3) и формулы (1.133)
и (1.134), |
приведем уравнение (1.303) к виду, аналогичному (1.12): |
|||||
|
|
г" (Б) + 2 (е) = |
- Л - |
ек «-*) (F0 + F cos |
at). |
|
|
|
|
|
ф(<) |
|
|
Подставляя сюда формулу (1.227), приближенно получаем |
||||||
|
|
2 |
(г) + це) |
- е |
в (I + 2В cos 2Ш) |
|
~ |
|
(FQ + |
Fcosat)(\ |
|
„к (1-х) |
{F0 + Fcosat — |
е |
— 25 cos26/) = -=—g |
— 2BF0 cos 26i! — BF [cos (ю — 26) t + cos (со + 26) *]}.
Используя приближенное равенство (1.15), запишем это урав нение в виде
г"(е) + 2(е) = -^-е (Т ~*) \pQ + Fcos~е |
— 2BF0cos2е — |
—c o s ( 4 — 2 ) е - b c o s ( - r + 2 ) e }•
Частное решение этого уравнения, полученное аналогично из ложенному выше (см. § 3, 5, 7), имеет вид
2(e) = Л ° |
Ti{4 |
+ a c o s ( - f e _ p ) + |
a c cos[( - f — 2 ) |
+ |
|||||
|
|||||||||
+ |
aycos |
-g- + |
2) е — р у |
+ |
а2 е cos (2s — р2 е)}, |
(1.323) |
|||
где |
|
|
|
|
|
|
|
|
|
а — |
|
8F |
|
|
|
|
|
||
|
|
|
|
|
|
Л 2 я 2 |
|
||
|
|
|
|
|
|
+ |
16 |
|
|
|
р = |
arctg |
(4Л2 л2 — л2 ) со2 |
+ 82 it2 |
' |
|
|||
/ [ ( c o - 2 e ) 2 ( - i ^ n 2 _ i ) |
+ |
e2 |
+ |
16- АС*П?(со—26)4 |
|
||||
|
|
|
|
|
|
|
|
я а |
|
|
|
|
ААспк |
(со |
— |
28)2 |
|
|
|
|
рс --- arctg • (4Лс 2 л2 — я2 ) (со — 28)2 |
+ |
8 2 я 2 ' |
|
1 В § 5 влияние сопротивления не рассматривалось.
101
|
|
|
|
|
|
QFB |
|
|
|
|
|
|
|
|
||
V |
|
|
|
У . я 2 |
— 1 |
+ |
0 2 1 |
+16 |
У |
|
(со + |
20)1 |
||||
|
|
я 2 |
/ |
|
|
J |
|
|
я 2 |
|
|
|
|
|||
|
|
|
|
|
4Лу пл (ш-f 29)2 |
|
|
|
|
|
||||||
|
|
р у = arctg • (4Л2 я2 — я2 ) (со + |
2G)2 + |
0 2 я 3 |
' |
|
|
|
||||||||
а 2 е |
= |
|
|
|
|
Рге |
|
|
|
|
16Л2дЛЛ |
|
|
|||
0 ( 1 6 Л 2 2 0 л 2 |
—Зя 2 ) * |
a r |
c |
t g |
16AW— Зл2 |
|
|
|||||||||
Возвращаясь в решении (1.323) к старым переменным в соответ |
||||||||||||||||
ствии с заменами (1.3) и (1.15) и учитывая приближенные |
равенства |
|||||||||||||||
(1.134) и (1.143), для стационарных |
колебаний |
получаем выражение |
||||||||||||||
/ (х) = |
- А . _|_а |
cos (at — р) + |
ас |
cos [(со — 2В) t — рс ] |
+ |
|||||||||||
|
+ |
ау cos [(со + 28) t — ру ] + |
а2 9 |
cos (20t — р2о). |
|
|
||||||||||
Отсюда видно, что амплитудно-частотные характеристики |
суб- и |
|||||||||||||||
ультрагармонических |
|
колебаний |
определяются |
|
по |
формулам |
||||||||||
/ (Лс ) = ас ; |
/ (Лу ) = ау; |
/ (Л2в) = а2в. |
|
Подставляя сюда |
приве |
|||||||||||
денные выше выражения для а^, ау |
и а2 в |
и используя приближенное |
||||||||||||||
равенство / (Л) =а Л У~а, для амплитудных |
кривых |
получаем урав |
||||||||||||||
нения |
|
|
|
|
ЛС | У — |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
± |
|
|
|
|
|
9FS |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
/ |
|
со Т 20). 1 |
4 < |
У |
1 |
+ 0 2 |
|
+ 16 |
с, у |
л 2 (со ¥ |
20)4 |
|||||
|
|
|
||||||||||||||
|
|
^ 2 0 |
= =F |
^ " 5 |
л 3 |
|
|
|
. |
|
|
|
(1.324) |
|||
|
|
|
|
|
QVa.(\&A2№n> |
— Зл2 ) |
|
|
|
|
v |
|||||
Знак «минус» в выражении |
(со =р 29) соответствует |
амплитудам Л с |
||||||||||||||
субгармоник, |
знак «плюс»—амплитудам Л у |
ультрагармоник. Пола |
||||||||||||||
гая в формуле (1.324) п = |
0, как и следовало ожидать, приходим к |
|||||||||||||||
выражению |
(1.321). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Следует |
отметить, |
что |
если |
использовать |
разложение |
(1.240), |
то аналогично изложенному выше (см. § 5) можно обнаружить коле
бания с частотами | со ± |
49|; | со ± 691; | со ± 801; |
а также коле |
||
бания с частотами 2iB, |
I = |
1, 2, 3, 4, |
амплитуды которых зна |
|
чительно меньше амплитуд |
колебаний |
с частотами |
| со ± 20 \ и 26. |
|
Устойчивость колебаний. Сначала |
рассмотрим |
простейший слу |
чай отсутствия трения, т. е. когда колебания описываются уравне нием (1.54), где параметры и \ определяются по формуле (1.282).
Как следует из изложенного в § 5, неустойчивыми могут быть только колебания в мягких системах (р < 0), характеристика кото рых имеет ненулевые корни (см. рис. 9, е). Используя выражения (1.32) и (1.58) и заменяя 0 на 0 (1 + 2В), аналогично (L247) получа-
102
ем следующее уравнение амплитудно-частотных характеристик коле баний, близких к неустойчивости:
Как показано выше, формула (1.212) остается справедливой и для случая несимметричных колебаний мягких систем. Используя вы ражения (1.212) и (1.296), для критического состояния имеем
9 ( 1 + 2 5 ) = - ^ . |
= 2ц. |
Из третьей формулы (1.296) видно, что величина ц. выражается через корни уравнения энергии (1.295). Выразить ее аналитически через параметры системы невозможно. Однако, как показывают вы числения, хорошую сходимость дает весьма простое выражение
е(1 + 2В) = 2 | х » | / Г - ^ - , |
(1.326) |
аналогичное формуле (1.250) для симметричных систем. Причем не симметрия системы учитывается параметром а), определяемым из уравнения (1.287) и первого равенства (1.285). Точность приближенно го выражения (1.326) тем выше, чем меньше несимметрия, т. е. чем меньше \у^\.
Подставляя формулу (1.326) в равенство (1.325), после простых преобразований получаем уравнение критических состояний
F = |
(1.327) |
|
m[n |
Входящая в эту формулу критическая амплитуда Л к р равна ненулевым корням характеристики уравнения (1.281), которые свя заны с ненулевыми корнями характеристики уравнения (1.284) пер вой зависимостью (1.292), т. е.
Ар = У = У* — - gjjj-r •
Ненулевые корни характеристики уравнения (1.284), как извест но [13], имеют вид
У * = ± ТГРГ ( ^ v 2 + 4 a : | P | ± У*
Следовательно, для критических амплитуд получаем два выра жения:
А«р = "QTRT f"Г" ± l / V ; + 4a:|P|). |
(1.328) |
2IPI |
|
Из этих-двух значений в формулу (1.327) подставляется то, которое дает минимальную по абсолютному значению амплитуду возбужде ния F.
Заметим, что принимая в формуле (1.327) Л к р согласно выраже нию (1.243) и полагая а* = а и б* = 0, получаем, как следовало
103
ожидать, уравнение (1.251) критических состояний симметричной системы. Полагая в формуле (1.327) попеременно F = 0 и со = О, получаем координаты точек пересечения кривой / критических со стояний с осями координат (рис. 71):
< и Fi =
|
V а |
|
|
а . |
min |
Для определения |
геометрического |
места точек срыва |
колебаний |
||
продифференцируем |
уравнение (1.54) |
по |
Л и положим |
|
= 0. |
|
Имеем |
|
|
|
|
|
(/'0 + |
/ 0 ' ) ( 0 2 - с о 2 ) + |
2 ( / 0 - |
||
|
|
— б*) 60' = ± 200'F. (1.329) |
о.сек
Рис. 71. Кривые критических состоя ний несимметричной системы без тре ния.
Это уравнение должно разре шаться совместно с уравнением амплитудно-частотной характе ристики (1.54), которое запишем так:
(/0 — SJ (92 — со2) = ± 0 2 F . (1.330)
Формулы (1.329) и (1.330) представляют собой параметри ческие уравнения кривой / / кри тических состояний (см. рис. 71). Координаты точек пересечения кривой // критических состоя ний с координатными осями, най денные так же, как в § 6 данной главы, следующие:
(1.331)
min
Использование выражений (1.329) и (1.330) затруднительно из-за громоздкости вычислений. Поэтому построим приближенную форму лу, приняв линейное приближение для частоты свободных колеба ний 02 = а;. Тогда 0' = 0 и выражение (1.329) принимает вид f У о ; (а; — со2) = 0.
Поскольку а ; ^ 0 и а ; — со2 Ф 0, то /' = 0, откуда / = с = = const, т. е. линейное приближение для частоты 0 (так же, как и для симметричных систем) влечет за собой осреднение амплитудной функции. Теперь уравнение (1.330) принимает вид
(с Val - 6J (о! - со2) = ± <F. |
(1.332) |
104
Используя граничные условия со = О, F = Fц, получаем
Здесь использована формула (1.331). Теперь уравнение критических состояний (1.332) принимает вид
F = |
( ± - г а * / з Т р Т + б * |
1 — со |
(1.333) |
|
а* |
|
На рис. 71 кривые критических состояний построены по форму лам (1.327) и (1.133) для системы с параметрами а = 1 се/с- 2 ; у =
= 0,15 см~х • сект"1; В = —0,2 см~2 • сект2; 6„ = —0,25 см • сект2; F0 = 0,25 см • сект2.
Кружками (кривая /) и точками (кривая //) представлены ре зультаты решения 1 на ABM МН-7. Как видим, совпадение аналити ческих и машинных решений можно признать удовлетворительным, если учесть, что машинные решения для кривой / занижены на 10—
15%. Физический смысл кривых и порядок пользования |
рис. 71 |
такой же, как на рис. 59 для симметричных систем (см. § 6 |
данной |
главы). |
|
Далее рассмотрим устойчивость колебаний с учетом вязкого трения.
Пользуясь формулой (1.110), записываем уравнение амплитудно-
частотной характеристики |
в виде |
|
F = ± - gr (/в - |
6J V(n* - со2 + 82 )2 + 4га2 |
(1.334) |
Используя выражения (1.32) и (1.334) и заменяя частоту 9 на |
9 (1 + |
+ 25), аналогично (1.248) получаем следующее уравнение для амп литудно-частотных характеристик колебаний, близких к неустой чивости:
± Q2 (1 + 2 S ) 2 |
=FB(I |
+ 2 5 ) Л ] / " а # - ^ - | В | Л 2 - (1.335) |
- б , У [п2 |
— со2 + |
92 (1 +- 2В)2 ]2 -+ 4/г2со2. |
Подставляя сюда приближенное выражение (1.326), приходим к уравнению критических состояний (кривая /)
F= |
А, |
(±-^rV2a*-\V\A2Kp- |
|
||
|
|
. (1.336) |
1 Результаты решения |
на АВМ, приведенные в этом параграфе, получены |
В. С. Горбатовым.
105
В формулу (1.336) подставляется то значение (1.328) критической амплитуды /4к р , которое дает минимальное абсолютное значение F.
Исследуя выражение (1.336) на экстремум, определяем координа ты минимума кривой / критических состояний:
со - / 4 - « : - / . « |
|
F =nV2 ±AK?V2am-\^\Alp-y^- |
(1.337) |
|
min |
Заметим, что приближенное равенство (1.326) справедливо в случае отсутствия трения. Поэтому формулой (1.336) можно поль зоваться для малых сопротивлений; при больших сопротивлениях естественнее воспользоваться приближенным равенством
/г2 + 6 2 ( 1 + 2 В ) 2 ^ - ^ - .
Подставляя это выражение в формулу |
С1.335), |
получаем уравнение |
критических состояний для больших значений |
/г и со <С Уа; г: |
|
± - ^ 1 / 2 а , - | 6 | 4 р |
|
со2 +4п2 со2 |
V V а* |
|
|
|
|
(1.338) |
Для приближенного определения геометрического места точек срыва колебаний в соответствии с изложенным выше примем в
уравнении (1.334) линейное приближение для частоты 02 + |
пг = ai |
||||||||
и осредним |
амплитудную |
функцию / |
с = const. |
Имеем |
|||||
F |
= |
• (с У а, |
— бJ У (а: — со2)2 + |
4n2co2. |
(1.339) |
||||
|
|
а; |
|
|
|
|
|
|
|
Постоянную |
с найдем |
из |
граничных |
условий |
со = |
О, F — F/t. |
|||
Используя формулу (1.231), получаем такое же значение с, |
как для |
||||||||
случая отсутствия трения. Подставляя его в формулу |
(1.339), прихо |
||||||||
дим к следующему уравнению критических состояний |
(кривая / / ) : |
F =
(1.340)
1 Ограничение со < J^a* вытекает из сопоставления аналитических и машин ных решений.
106
Этим уравнением следует пользоваться для умеренных значе
ний п. При |
больших значениях п лучшую |
сходимость с машинными |
||
решениями |
дает |
уравнение. |
|
|
|
|
F = - 4 ( ± - ! - « , / V r p T + |
|
|
|
+ |
l / ( a : — со2)2 + 12шо |
- со)2 |
(1.341) |
На рис. 72 приведены построенные по формулам (1.336) — (1.341) графики критических состояний для системы с теми же параметрами,
|
|
|
|
|
|
|
|
а.сен |
|
|
1 |
/ |
о,сек' |
Рис. 72. |
Кривые |
критических |
состоя |
Рис. 73. |
Кривые |
критических |
состоя |
||||||
ний несимметричной системы при вяз |
ний несимметричной системы при тур |
||||||||||||
ком трении: |
|
|
|
|
булентном |
сопротивлении: |
|
||||||
|
|
|
|
|
|
|
|
|
|
||||
/ |
— |
л |
= |
0,05 сек-1; |
2 |
— л = |
0,20 |
сек— |
t — п — |
0,05 |
cm~1;2 |
— л = 0,20 |
см~1. |
3 |
— |
п |
= |
0,30 |
4 |
— л = |
0,50 сек~1 |
|
|
|
|
|
что и на рис. 71. Здесь же условными обозначениями представлены результаты решения на ABM МН-7. Как видно, совпадение резуль татов можно признать удовлетворительным. Порядок пользова ния графиками такой же, как и для симметричных систем (см. § 6 динной главы).
Рассмотрим влияние турбулентного сопротивления |
на |
устой |
||
чивость колебаний. |
Корректируя путем |
замены 0 на |
0 (1 + 2В), |
|
а на а.^ и F0 на |
амплитудно-частотную |
характеристику |
(1.308), |
|
для несимметричных колебаний мягких систем, близких |
к неустой |
чивости, имеем
107
0(1 + 2B)F |
|
|
6 (1 + 2В) ^ |
' |
Л2 Я2 д |
] / [ с о 2 ( ~ - л - 1 ) + 8* (1 + 2Bf |
+ 16 |
— ш4 |
Подставляя сюда приближенное выражение (1.326), получаем урав нение критических состояний (кривая /)
+ 1 6 |
(1.342) |
В эту формулу подставляется то значение (1.328) критической ам плитуды Л к р , которое дает наименьшее абсолютное значение F.
Параметрические уравнения геометрического места точек срыва колебаний (кривая // ) весьма громоздки и неудобны для использо вания, так как они содержат амплитуду колебаний. Вычисления по казывают, что хорошее совпадение результатов дает простое вы ражение 1
|
F |
= |
|
|
+ 8^ У |
( о ; - со2)2 + 6Д, (У< - ш)3> |
(1.343) |
где Ал |
абсолютная величина минимума характеристики |
системы, |
3|PI Imin"
На рис. 73 изображены построенные по формулам (1.342) и (1.343) кривые критических состояний для системы с теми же параметрами, что и на рис. 71, при различных значениях п. Дискретными ус ловными обозначениями представлены результаты решения на ABM МН-7. Как видно, совпадение результатов удовлетворитель но. Порядок пользования графиками такой же, как и для сим метричных систем.
1 Эмпирические формулы, приведенные в этом параграфе, подобраны В. С. Гор батовым.
Г Л А В А I I
БИГАРМОНИЧЕСКОЕ
ВОЗБУЖДЕНИЕ
§ 1. Симметричное |
возбуждение |
|
|
систем без трения |
|
|
|
Будем рассматривать |
(§ 1,2) |
возмущение в виде двух |
|
гармоник: F (t) = Fx cos coji + ^2 c o s |
F (0 = |
^1 s ' n W J + |
|
-f- F2 sin co2*. Как известно, эти |
функции |
будут |
периодическими |
только в случае кратности частот, т. е. со2 = р.©!, ц = 0, 1,2,3, ...
Так как только при периодическом возбуждении возможно появле ние стационарных колебаний, то в дальнейшем будем рассматривать случай кратности частот, т. е.
F (t) = Fx cos со* + F2 cos LICO*; |
F (t) = ^ s i n at + |
F2 sin p,co*. (II. 1) |
||||||||||||
Частота со называется основной (низшей) частотой возмущающей |
||||||||||||||
силы F (t). При р, = 0 первая формула (11-1) дает |
несимметричное |
|||||||||||||
гармоническое возбуждение (см. § 7 гл. I), вторая формула |
( I I . 1) — |
|||||||||||||
гармоническое возбуждение |
(см. § 1—6 гл. I). При |
L I = |
|
1 формулы |
||||||||||
( I I . 1) описывают |
гармоническое |
возбуждение |
с частотой |
со и ам |
||||||||||
плитудой Fx -f- F2. |
Несмотря на аналогичность |
выражений |
(II.1), |
|||||||||||
описание |
стационарных колебаний |
|
для них |
несколько |
различно. |
|||||||||
Поэтому |
будем |
рассматривать их раздельно |
и |
именовать |
соответ |
|||||||||
ственно |
четным |
(косинусоидальное) |
и нечетным (синусоидальное) |
|||||||||||
возбуждениями. |
|
|
|
|
|
|
|
|
|
|
|
|
||
Четное возбуждение. Рассмотрим |
[34] стационарные |
вынужден |
||||||||||||
ные колебания |
нелинейной системы, описываемые уравнением |
|||||||||||||
|
|
x(t) |
+ R(x) |
= F x |
cos со* + F2 cos LICO*. |
|
|
|
(II.2) |
|||||
Используя замены |
(1.3), преобразуем нелинейное уравнение |
(П.2) |
||||||||||||
к линейному, аналогичному (1.12), т. е. |
|
|
|
|
|
|
||||||||
|
г" (е) + |
г (е) = —^— (F± |
cos со* + F2 cos ixco*). |
|
|
|
||||||||
|
|
|
|
|
ф(9 |
|
|
|
|
|
|
|
|
|
Линеаризуя |
фазовую функцию в соответствии с выражением (1.13) |
|||||||||||||
и используя |
зависимости (1.14) и (1.15), получаем |
|
|
|
|
|||||||||
|
г" (е) + г(е) = -L |
(F± cos |
|
е + F2 cos ц -f- е). |
|
|
(П.З) |
109