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

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

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

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

Добавлен: 22.10.2024

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

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

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

Здесь

 

 

 

 

 

 

 

 

 

 

 

 

 

k\ = Vl-kl;

^

- _

^

=

/

l

f

U J

|

^ _

.

< i .

Формула (1.204) для частоты с учетом выражений (1.194) и (1.206)

приобретает известный [2, 13, 35, 36, 57] вид

 

 

 

 

 

 

0 =

-ъРЬг

= —

 

 

=

 

"

,

Voc + ВЛ2.

(1.207)

Для

мягких

систем

(а >

0;

В <

0)

выражения

(1.194)

можно

записать так:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

- Y

 

J l

'

^

<

1.

(1.208)

Формула (1.204) для частоты, с учетом выражений (1.208) при­

обретает известный [2, 13, 35, 36, 571 вид

 

 

 

 

 

 

 

 

е — w / a " ^ - | p , i 4 " '

 

 

 

( L 2 0 9 )

где й определяется из второго выражения (1.208).

Перейдем к рассмотрению фазовой функции во втором прибли­ жении. Подставим первую формулу (1.204) в правую часть прибли­

женного равенства (1.203). Имеем

 

Ф = 0 / + - J^ - sin0/cos0 i = 0 / + В sin 20^,

(1.210)

где

 

Выражение (1.210) подтверждается решениями на ЭЦВМ «Урал-3» [10, 13]. Максимальное значение амплитуды В будет иметь место

при k = 0 и К (k)

Подставляя это значение в формулу (1.211),

приближенно получаем

 

 

(I.2I2)

Здесь использована первая формула (1.202).

Выражение (1.212) справедливо для мягких систем > 0, В <; 0). Значение k определяется по второй формуле (1.208). Для жестких систем (а > 0, В > 0) значение k = ikQ становится мнимым. Поэтому, подставляя соотношение (1.205) в равенство (1.212), имеем

л

(1.213)

 

Здесь и k\ определяются по формулам (1.206).

Графики амплитуд В, построенные по формулам (1.212) и (1.213), изображены на рис. 50, где т = - t ^ - -

64


Рассмотрим системы с перескоком (а <с О, Р >• 0).

Для больших колебаний решение (1.198) остается справедливым, если поменять знак а на противоположный. Следовательно, будут справедливыми также формулы (1.206), (1.207), (1.210) и (1.211). В частности, в соответствии с выражением (1.42), меняя знак а на

противоположный в равенствах

(1.206) и (1.207), получаем извест­

ные [13] точные формулы для частоты больших

свободных

коле­

баний:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6 =

2К {k) У р Л 2

— | а | ;

k =

/ 2 ( P i 4 s - | « l )

(1.214)

Так как для

больших колебаний рЛ2

>. 2 | а | , то 1 ;> k ;>

j / l / 2 -

Для определения максимальной амплитуды В подставим в форму­

лу

(1.211)

К (У"1/2) «

1,85. Будем иметь В

4АГ(У1/2) те 0,425 а0.

Используя

первое

выражение

 

\В\

 

 

 

 

 

 

(1.202), получаем

 

 

 

 

Т<0

\

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

- )

В = 0,425

 

 

 

 

 

 

 

В>0\

 

 

 

 

 

 

 

 

0,10

 

 

 

 

 

 

 

 

 

 

 

(1.215)

 

от

 

 

 

 

 

 

Для малых колебаний си­

 

 

 

 

тУо~

 

стем с перескоком

воспользо­

 

 

 

 

 

 

 

 

 

 

 

8<0

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

зультатами

затруднительно,

 

 

 

 

/

 

г

\т\

поскольку

в

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

с

 

Рис. 50. График амплитуд В уточненной

условиями

(1.38),

а

именно:

 

 

фазовой функции систем с кубической ха­

| а |

< . РЛ2

<

2 |а |,

 

первая

 

рактеристикой.

 

 

 

формула (1.194)

дает

мнимое

 

 

 

 

 

 

 

 

значение параметра г|). Для того

чтобы избежать этого, будем ис­

кать решение уравнения (1.192) при начальных

условиях

 

в виде

 

 

 

 

t = 0,

X = Л,

х = 0

 

 

(1.216)

 

 

 

 

X = A cn (т|)/, k) =

Л сп и.

 

 

(1.217)

 

 

 

 

 

 

 

 

Подставляя решение (1.217) в уравнение (1.192) и группируя

члены, имеем Л [— ip* (1 — 2k2

sn2 и) — \а\ + рЛ2

сп2 и] сп и = 0.

Используя известную зависимость сп2 и = 1 sn2

и и приравнивая

нулю квадратную

скобку,

получаем

 

 

 

 

 

 

 

(2А*ф« рЛ2 ) sn2 « - f (рЛ2 — | а | — ip2) = 0.

 

Это

уравнение

удовлетворяется

при условиях

21гЩ>2 — рЛ2 = 0;

РЛ2

|ос | — ip2

=

0,

откуда

находим

 

 

 

 

 

 

 

 

ф = У"рЛ2 | о | ,

А2

=

 

рМ2

 

 

(1.218)

 

 

 

 

2 ( р Л 2 - | с с |) '

 

 

 

 

 

 

 

 

 

 

 

 

 

В соответствии с формулами (1.196) и (1.45) решение уравнения (1.192) методом переменного масштаба для малых колебаний при

Б 4-5

65


начальных условиях (1.216) имеет вид

Будем рассматривать случай реализации максимальных ампли­

туд, т. е. когда Л да 1|/~2

.

Тогда

 

У 2

/ 2 Р

Г 2

у 2 р

Следовательно, решение (1.219) приближенно можно записать

так:

= . 1/2Л* [1 + cos cp (г)] = Л 1 Icos 1/2 ф (01 *. Отсюда

 

•xssi4cos-j-q>(')-

(1.220)

Сопоставляя решения (1.217) и (1.220), имеем

I cos-i-tpCO^criCipf. k) = cnu, т. е. - | - = amu. . .(1.221)

Поскольку для малых колебаний | а | <; рЛа <;2 | а |, вторая фор­ мула (1.214) дает k >• 1. Поэтому воспользуемся известными [23, 55]

свойствами

эллиптических

функций:

 

 

 

 

, ...,

 

 

 

 

 

сп (и, k) = dn [ku,

- | - j =

 

 

 

 

d JB_

 

 

 

. .

 

 

2я2я

 

V у

qn

 

nnu

 

 

 

2

 

 

я

 

 

 

 

 

 

 

dku

 

" « ( 4 f

1

* ( + )

 

\

1

C Q

S .

 

 

 

 

 

 

 

 

i " ^ " " " > f j y

Интегрируя,

получаем

 

 

 

со

 

 

 

 

 

 

Ф

o m

 

 

 

i

 

 

 

 

nnku

 

 

 

 

 

о V

1

 

 

 

^

=

И

П

В =

^ 7 Х Т +

2

2 - • T T ? ^

s i n

 

 

 

(±)

 

 

 

 

 

*(±-)'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Сохраняя здесь один член ряда и

подставляя

и = ip/,

получаем

приближенное равенство ф (t)

= 8г + В sin6/, где

 

 

 

-

e

= - r V

 

5 =

- г т ^ '

 

^ =e x p ( - f - ) -

(L222 >

Принимая во внимание формулы (1.218) и вводя обозначение

 

 

 

 

 

^ = 4 -

=

т ^ - | / 2 ( р Л 2 - | а | ) ,

 

(1.223)

в соответствии с первым выражением (1.222) получаем точную фор­ мулу для частоты малых колебаний систем с перескоком:

6 = VIК (*.) .

(1.224)

66


Выражение для амплитуды В, в соответствии с последними двумя» формулами (1.222), можно представить так:

 

В =

 

 

 

 

 

 

 

= 2sech

я/С (А.)

 

 

 

 

 

 

 

 

/С (ft,)

 

 

 

 

 

+ 4

 

е х р ( — J + exp(

—J

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.225)

Здесь использованы известные [23, 55] соотношения

-Y

 

 

К'(К)

= К(К);

- ^

= - ^ j / 2 | a [ - p \ 4 2 .

(1-226)

 

Графики амплитуд Б для систем с перескоком, построенные по

формулам

 

(1.215)

и (1.225), приведены на рис, 51, где обозначено

 

 

 

 

 

 

 

 

Л2 .

 

 

 

 

 

Основной

тон. Дифференци

 

 

 

 

руя

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

равенство

°>г

 

 

 

(1.210)

по

времени

t, получаем

 

 

 

 

 

 

ср (0 =

6(1 +2Bcos2'e/).

0,1

 

 

 

 

 

 

 

 

 

 

 

(1.227)

0

 

 

 

 

Подставляя эту формулу в урав-

 

 

 

 

нение'(1.12),

имеем

 

 

 

 

 

 

« . . .

-

,

 

 

Fcased

Рис. 51. График амплитуд В 'уточнен­

z

~

т

г

К

* )

 

ной фазовой функциисистем- с

пере­

 

 

в(1 + 2Вcos2№) *

 

 

 

скоком.

(1.228)

Далее будем рассматривать умеренно нелинейные системы, для которых <^ 1. Уравнение (1.228) можно представить так:

г" (е) + г (е) = cos со* (1 — 25 cos 2Ы) =

={cos со* — В [cos (со — 20) ^ -f- cos (со + 20) t]}.

Используя

приближенное равенство

(1.15), запишем это уравнение

в виде

 

 

 

 

г" (е) + г (е) = -J- {cos - f - в -

В [cos

- 2) е + cos

+ г) ej J.

 

 

 

 

(1.229)

Частотное решение этого уравнения будем искать в форме

г (е) =

Сх cos -J- е + Са cos

— 2J е + Cs cos

+ 2J 8. (1.230)

Подставляя это решение в уравнение (1.229) и приравнивая коэффи­ циенты при членах, содержащих косинусы с одинаковыми частотами,

б*

67


определяем

 

 

 

С2 =

FB

 

 

 

1 —

I

е2

1

 

 

 

С3

=

 

- -

 

 

 

 

е 1

 

Теперь частотное решение (1.230) можно записать так:

со cos-g- е

г(е) = 9^

6» — со»"

г

/

со

\

/

со

\

— В

о з ( т

- 2 ) в

« ( - ё -

+ 2 ) .

В2

— (со — 26)2

+ б3

— (со + 29)а

 

Переходя здесь к старым переменным в соответствии с формулами (1.3), (1.15) и (1.32), находим решение для стационарных колеба­ ний:

cos со<.

 

"

cos (со — 29) I

+

* | Д с + - | - л : 2 = 8F{- еа

со2

— 5

 

б 2 — (со - 20)2

cos (со +

26) t

 

 

(1.231)

+ б 2 - • (со +

26)а

 

 

 

 

по форму-

Здесь В определяется по формулам (1.212) и (1.213), а 9

ле (1.33) или (1.207) и (1.209).

 

 

 

 

 

Как видно из решения (1.231), помимо гармонических

колебаний

С частотой возмущения со имеют место колебания с большей частотой

(со + 29), которые условимся

называть ультрагармоническими 1 ,

и колебания с меньшей частотой

(со 20), которые будем называть

субгармоническими. Последние

удовлетворяют

условию

|со 2 0 | < с о .

(1.232)

Исследуем закономерности изменения этих колебаний в зависимости от соотношения частот

Сначала

 

о

 

 

 

 

 

 

рассмотрим ультрагармонические колебания. Сохраняя

в решении- (1.231) только член с частотой

(со+20)

и полагая х = Ау\

cos (со +

20) t = ± 1,

получаем

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

ристики

для

ультрагармонических

колебаний:

 

 

 

 

Ау ]/«

+ •

А2У =

±

 

QFB

 

(1.233)

 

 

б 2

—(со +

26)2

 

 

 

 

 

 

 

Порядок вычисления по этой формуле таков. По заданной часто­ те со, пользуясь уравнением (1.34), находят А. Затем по формулам (1.212), (1.213)' и (1.33) или (1.207), (1.209) вычисляют В и 0 и, на­ конец, из выражения (1.233) определяют амплитуду Ау.

На рис. 52 приведены амплитудные кривые ультрагармонических колебаний при а = 1 сект2. Как видно из рис. 52, амплитуды ультра-

1 Принятая здесь терминология несколько отличается от предложенной ранеё"! [57]. •

68