
книги из ГПНТБ / Бондарь, Н. Г. Нелинейные стационарные колебания
.pdfхарактеристике становится вертикальной. Следовательно, Для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 8а |
12(a) 2 |
' |
|
где Ах |
= у% (0) — максимальная амплитуда |
колебаний, |
связанная |
с амплитудой А — у (0) зависимостью, вытекающей из первой фор мулы (1.292),
Л 1 = л - Ж * |
( L 2 9 4 ) |
Для сопоставления точности формул (1.60) и (1.293) и определе ния пределов их применения необходимо иметь точную формулу.
89