Файл: Бондарь, Н. Г. Нелинейные стационарные колебания.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 22.10.2024

Просмотров: 115

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Зависимость между

амплитудами устанавливается

формулой

(1.301): Ах = А + 0,84.

В табл. 1 даны также

относительные

ошибки вычислений по формулам (1.60) и (1.293).

 

 

Как видно из таблицы, для малых амплитуд колебаний (вблизи

А1 =

0) следует пользоваться формулой (1.293), для

умеренно боль­

ших

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

формулой

(1.60) и, наконец, для больших амплитуд следует пользоваться точ­ ными формулами (1.296) и (1.297).

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

При жесткой характеристике будет иметь место только один вещественный корень уравнения (1.287) и, следовательно, одно зна-

Т А Б Л И Ц А 1

Относительная ошибка

 

А,

 

е,

в,

в»-в, 100%

e

f f

-e,

100%

-0,5

0,34

0,74

0,78

0,62

-3,4

 

 

16,2

 

-0,75

0,09

0,76

0,75

0,77

1.3

 

 

-1,3

 

-0,84

0

0,78

0,73

0,78

6,4

 

 

0

 

-1,00

—0,16

0,77

0,70

0,75

9,1

 

 

2,6

 

-1,25

—0,41

0,74

0,64

0,55

13,5

 

 

25,6

 

-1,55

—0,71

0,61

0,58

0,09

4,9

 

 

85

 

чение а* >

0. Условия определения

частоты свободных

колебаний

по точным и приближенным формулам такие же, как для мягкой характеристики.

При полужесткой характеристике определение частот аналогич­ но определению их при мягкой характеристике.

В системах е перескоком (полумягкая характеристика) будут иметь место два ненулевых вещественных корня уравнения (1.287), для которых а* > 0. Они определяют собой два устойчивых положе­ ния равновесия. Для определения малых колебаний относительно этих положений можно пользоваться как точными, так и прибли­ женными формулами. При определении больших колебаний следует пользоваться точными формулами для корня уравнения (1.287), при котором а* <с 0. В этом случае можно также пользоваться фор­ мулами (1.63). Однако, поскольку вычисления приходится произво­ дить путем последовательных приближений, то предпочтительнее пользоваться точными формулами.

Ангармонические колебания \ Как показано выше, несимметрич­ ное возбуждение эквивалентно симметричному возбуждению не-

1 Ангармоническими колебаниями иногда называют [1] несимметричные ко­ лебания.

92


симметричной системы. Это позволяет воспользоваться полученны­ ми выше результатами.

Так, в случае отсутствия трения, используя формулу (1.58), получаем выражение для амплитудно-частотных характеристик

f [ - F -^шах\ =

7Г— ± -

\

min /

°

Совершенно аналогично в случае наличия вязкого трения, ис­ пользуя формулу (1.110), имеем

/ ( = М mto) 8 W — со2 + 92 )а + 4па со2

Для определения наибольших по модулю амплитуд колебаний можно воспользоваться формулой (1.111):

В частности, для стационарных колебаний, описываемых уравнением х -f- 2пх + ах + fix3 = F0 + F cos coif, используя формулы (1.32) и (1.291), получаем следующее выражение для амплитудно-частот­ ных характеристик:

S i n "

+

2 v

=-&-±

г

6 f

2 )2 + 4n2co2

(1-302)

 

9 =

/ ( / t 2

со2 + 9

 

В приведенных выше формулах верхний знак соответствует макси­ мальной амплитуде Л т а х , нижний — минимальной амплитуде Л т ш, а частота свободных колебаний 8 вычисляется в соответствии с ука­ заниями предыдущего раздела.

На рис. 68 изображены амплитудно-частотные характеристики, построенные по формулам (1.302) для системы с параметрами а = = 1 сек-2; 6 = 0,6 см-2 • сект2, п = 0,1 сект1.

Здесь же точками представлены результаты решения 1 на ЭЦВМ «Урал-3». Как видим, совпадение аналитических и машинных ре­ зультатов можно признать удовлетворительным.

Более подробно2 рассмотрим стационарные колебания при тур­

булентном

сопротивлении, которые

описываются уравнением

 

х + nx2sgnx-f

R (х) = F0 + F cos со*.

(1.303)

Используя

(1.3) и формулы (1.14), (1.15), (1.133) и (1.134), преобра­

зуем уравнение (1.303) к виду

 

 

 

 

г" (8).+2 (е) = 4 - е

(Т _Т)

(F0 + F cos - f в).

 

1Результаты получены Н. Г. Новиковой и В. В. Потаповой.

аЭтот вопрос необходимо рассмотреть более подробно, так как в § 3 данной главы не рассматривались несимметричные системы.

93


Частное решение этого уравнения, полученное аналогично.изло­ женному выше (см. § 2 и 3) и определяющее собой стационарные колебания, имеет вид

г(е) = 4 - + a c o s ( - f e - p )

где а и р определяются по формулам (1.136). Возвращаясь к. старым переменным в соответствии с заменами (1.3) и (1.15) и учитывая

Рис. 68. Амплитудно-частотные характеристики несимметричных колебаний при вязком трении для различных F, см • сек~2:

а — в о з б у ж д е н и е F + F cos at; б — в о з б у ж д е н и е F + Fcos 2a>t.

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

-д2- + acos(co* — р) ехр {пх • sgn*).

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

х = + Л т а х ; cos (со*— р) = ± 1 ; sgn х =

1, полу-

чаем уравнения

min

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

 

амплитудно-частотных

 

 

Fn

min

(1.304)

94

Здесь амплитудная функция определяется по формуле (1.130). Используя приближенное равенство (1.143), упростим уравнение амплитудной кривой к виду

/ ( = F 4 n « ) = - § - ± a .

(1.305)

Рассмотрим два примера.

Исследуем простейший случай линейной характерцстики, т. е. когда в уравнении (1.303) R (х) = ах. Используя равенства (1.32),

 

 

а

 

 

б

 

 

Рис. 69. Амплитудно-частотные характеристики несимметричных коле­

баний линейной системы

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

сопротивлении для различ­

ных F,

см • сек~2:

 

 

 

 

 

а — в о з б у ж д е н и е F + Fcos cut; б

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

F + F cos 2a>t.

 

(1.136), (1.146) и (1.147), по формулам (1.304) и (1.305) получаем сле­

дующие выражения для амплитудно-частотных

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

более точное

 

 

 

 

 

{ F

 

[1 +

^2/гДпах — l j exp jq=2n^m aX jj

=

 

 

QF

~\

 

 

-if" ±

г

 

 

 

 

* е х Р (=F пАтйх\;

менее

точное

 

 

 

 

(1.306)

 

 

 

 

 

 

 

(1.307)

На рис. 69 изображены

амплитудно-частотные характеристики,

вычисленные 1

по формулам

(1.306) (сплошные

линии)

и (1.307)

1 Вычисления

выполнены Л. И. Кожуховой.

 

 

95


(штриховые линии) на ЭЦВМ «Промшь» для системы с параметрами а = 1 сект-2; п = 0,5 см~х. Здесь же точками представлены резуль­ таты решения 1 уравнения на ЭЦВМ «Урал-3». Из рисунка видно, что, как и следовало ожидать, более точная формула (1.306) лучше соответствует машинным решениям. Однако менее точная и более

Рис. 70. Амплитудно-частотные характеристики несимметричных колебаний нели­ нейной системы при турбулентном сопротивлении для различных F, см • сек"2:

а — в о з б у ж д е н и е F + F cos at; 6 — в о з б у ж д е н и е F + F cos 2at.

простая формула (1.307) удовлетворительно соответствует машинным решениям. Это дает основание рекомендовать ее для более сложных характеристик.

Рассмотрим случай кубической характеристики R (х) = ах + + Р*3 . Используя равенства (1.32), (1.136) и (1.149), по формуле (1.305) получаем следующее приближенное выражение для ампли­ тудно-частотных характеристик:

min '

z

 

QF

 

 

. (1.308)

в ^

2

|

Л 2 я 2

 

2 " - П - 1 +

+ 16

—со"

На рис. 70 изображены амплитудно-частотные характеристики,

1 Результаты получены Н. Г. Новиковой и Л. И. Кожуховой.

96

вычисленные 1

по формуле (1.308) для системы с параметрами а =

= 1 сек-1; В =

1 см~2 • сек—2; п = 0,1 см-1.

Точками представлены результаты решения 2 на ЭЦВМ «Урал-3». Как видим, совпадение аналитических и машинных решений удов­ летворительно.

Суб- и ультрагармонические колебания. Для выяснения законо­ мерностей стационарных колебаний с частотой, отличной от частоты возбуждения со, необходимо прежде уточнить фазовую функцию с несимметричной характеристикой. Для этого найдем приближенное

решение уравнения

(1.281). Будем искать его в виде

 

 

 

 

х = N + А сп и,

 

 

 

(1.309)

удовлетворяющем начальным

условиям

 

 

 

 

 

 

* = 0;

x = N + А;

х = 0.

 

 

(1.310)

Подставляя решение (1.309) в уравнение

(1.281) и

группируя

члены,

имеем [(2fe2ip2 — РЛ2 ) sn2

и + (ВЛ2

+ а* — ф2

+

3pW2 +

+ ЗРУУЛ сп и) ] X Л сп и + (aN + $N3

=

0. Это

равенство

может

быть удовлетворено,

если

величины

в

круглых

скобках

равны нулю, т. е.

 

 

 

 

 

 

 

 

 

 

2AV

— М ' = 0;

cVV + pW3 — 6* = 0;

 

(1.311)

 

РЛ2 +

а* — г|>2 +

З Р # 2

+ ЗРЛМ сп и = 0.

 

(1.312)

Равенства (1.311) могут быть удовлетворены точно; причем корни второго уравнения (1.311) определяют нулевые точки харак­ теристики уравнения (1.281). Равенство (1.312) может быть удовле­ творено приближенно, если заменить в нем функцию сп и некоторой постоянной г%. Введя обозначение пх = 3$N2 + 3$ЫАг^, запишем приближенно равенство (1.312) так: РЛ2 2 + «i ~ 0. Раз­ решая это уравнение совместно с первым равенством (1.311), при­ ближенно получаем

+

+ а * + п » ^ w + w -

( Ш З )

Как известно [13], приближенное решение уравнения (1.281) методом переменного масштаба при начальных условиях (1.310) имеет вид

 

(N + А)

+ - i - Р (N + Л)2 - J

С О 8 ф ( 0 * 4 - -

1

Вычисления проведены Л. И. Кожуховой на ЭЦВМ «Наири-О

2

Результаты получены Н. Г. Новиковой и Л. И. Кожуховой.

7 4-5

 

 

97


Осреднив в левой части этого равенства х2 да пг

(N + А)2,

запишем

его приближенно

так:

 

 

 

 

 

 

х =

 

 

 

(N + A) ] Д ^ +

-1-6 (ЛА +

/ 4 ) 2 -

^-faffl

+ Л)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б

 

 

 

 

(1.314)

 

 

 

0 coscp(0 +

-7n

 

 

 

 

 

 

 

 

Сопоставив

равенства (1.309)

и (1.314), приравняем амплитуды

и свободные

члены:

 

б,

 

 

 

 

 

 

N =

 

 

 

 

(1.315)

 

 

У а, + -L Р/ц (;V +

А)2

 

А =

(N +

А) У а,

+ 4~ Р

+

'

 

 

 

 

 

 

 

 

Л)2

 

 

Разрешая эту систему, получаем

 

 

 

 

 

п0 =

1;

a * - j - - f В(ЛГ +

Л)*

1.

(1.316)

Из последнего равенства можно найти коэффициент п, если вос­ пользоваться какой-либо из формул для частоты Э. В частности, для а* + пг > 0, как видно из формул (1.313), можно воспользоваться известной [7, ГО, 13, 35, 36, 57, 69] формулой Дуффинга

0

=

+

Р

2

 

(1.317)

 

 

 

«» + Ч

 

 

 

Подставляя это выражение

во вторую

формулу

(1.316),

находим

a* +

«i =

б 2

-

-

f И ' -

(1.318)

 

W2 а* + — P ( W + .4)2

Сопоставляя равенства (1.314) и (1.309) при условиях (1.315), приходим к выражению (1.219), откуда следуют приближенные равенства (1.201) и (1.203). Следовательно, остаются справедливыми формулы (1.210), (1.211) и (1.212). Последние две формулы могут иметь модификации в зависимости от типа характеристики и, сле­ довательно, модуля полного эллиптического интеграла первого рода, который может оказаться мнимым или большим единицы. В част­ ности, для систем с небольшой несимметрией будут справедли­

вы следующие

формулы: (1.212) — для мягких

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

(см. рис. 9, в),

(1.213) —для жестких характеристик

(см. рис. 9, а),

98