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

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

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

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

Добавлен: 22.10.2024

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

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

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

позволяет упростить формулы (1.140) к виду

(Л) = е

 

А

_1_

= е " л / ( Л ) .

(1.143)

пА

[2 ( R (и)(u

du]-

I

 

 

о

2

'

 

 

 

 

 

Здесь использованы обозначения (1.21) и (1.22). Подставляя выраже­ ние (1.143) в формулу (1.138) и учитывая первое равенство (1.136), получаем уравнение (1.99), если п заменить на k. В случае колебаний с малыми амплитудами при малых сопротивлениях (пА <^ 1) условие (1.141) приближенно можно записать так:

со2

пгl ) + 02 = 0.

(1.144)

Подставляя (1.144) и (1.143) в (1.139), получаем уравнение

кривой

 

Л И Г -

< 1 Л 4 5 >

пересечение которой с кривой (1.144) приближенно дает максимум амплитудно-частотной характеристики. Заметим, что для An <^ 1 кривая (1.144) близка с скелетной кривой.

Простые характеристики. Рассмотрим стационарные колебания в случае линейной характеристики. Они описываются частным реше­

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

х (f) + пх2

(t) sgn х -f- ах

(t) = F cos at.

Исполь­

зуя вторую 1

формулу

(1.140), определяем

амплитудную

функцию

для х — А:

 

 

 

 

 

 

U (А) =

[21

аие2пи

du^

= -L / - f \ i n A (2пА - 1) +

1].

Подставляя это выражение в равенство (1.139), получаем уравнение амплитудно-частотной характеристики

±У^{2пА-1)+Ц =

= ±

г. .

арРпА

(1-146)

Здесь значение 9, как известно [7, 13], может быть определено по формуле

e - y ^ ( i - 3 i * L n M » ) .

Менее точное уравнение амплитудно-частотных характеристик по­ лучаем, используя приближенные равенства (1.143), где в соответ-

А

ствии с формулами (1.21) и (1.22) f(A) = [2 ] awdu] = А У~а. Под-

о

1 К этому же результату приводит и первая формула (1.140).

50


ставляя выражение (1.143) в уравнение (1.139), приближенно на­ ходим

 

АУ^

= ±

 

,

 

,

 

в Р

. '

 

 

 

(1-147)

 

 

 

 

 

 

 

 

 

 

 

+ 1

+ 16

 

 

 

 

Для оценки точности полученных результатов на рис. 43 построе­

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

при а =

1 сект2,

п =

= 0,5 см.—1. Сплошными линия­

 

 

 

 

 

 

 

ми показаны кривые,

построен­

 

 

 

 

 

 

 

ные по формуле (1.146), штрихо­

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

формуле

(1.147);

точками

пред­

 

 

 

 

 

 

 

ставлены

результаты

решения,

 

 

 

 

 

 

 

полученные на ЭЦВМ «Урал-3»1 .

 

 

 

 

 

 

 

Как

видно

из

рисунка,

совпа­

 

 

 

 

 

 

 

дение аналитических

и

машин­

 

 

 

 

 

 

 

ных

решений

можно

признать

 

 

 

 

 

 

 

удовлетворительным. Для

опре­

 

 

 

 

 

 

 

деления

амплитуд стационарных

 

 

 

 

 

 

 

колебаний

вблизи

от

резонанса

 

 

 

 

 

 

 

следует пользоваться

более

точ­

 

 

 

 

 

1,6 о.сек*

ным

выражением

 

(1.146). В ос­

Рис. 43. Амплитудно-частотные

харак­

тальных

случаях

допустимо ис­

пользование

более

простого

теристики линейной

системы при тур­

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

для различ­

выражения (1.147). Следует отме­

 

 

.—2

 

 

 

 

тить, что впервые

приближен­

ных F, см • сек

 

 

 

 

 

 

 

 

 

 

 

ное

решение

рассмотренной задачи

получено 2

в

работе

[71 ].

Далее рассмотрим колебания маятника (см. рис. 1), описываемые

уравнением •ф +

/г-ф2 sgn яр -[—у- sfn гр = М cos со/. Используя

вторую

формулу

(1.140),

определяем

 

 

 

 

 

 

 

 

 

 

 

 

 

а

 

 

, 1

 

4gn

e2nAsinA

+

 

 

 

 

 

 

I

sin ue~ du

 

 

 

 

 

 

 

 

+

2 )

 

 

 

 

 

 

 

 

 

 

 

 

 

1_

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

Подставляя это выражение в равенство (1.139), получаем

уравнение

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

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

 

 

 

 

 

 

 

 

 

4gn

 

e2nAsmA

+

1 (l-e2nAcosA)

 

 

 

 

 

 

Ц1 + 4 / г 2 )

 

 

 

 

 

 

 

 

 

 

 

 

=

±

 

 

 

 

 

 

QFenA

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

Результаты,

приведенные в данном параграфе, получены Н. Г. Новиковой

и Т. П. Лотаревой.

 

 

 

 

 

 

 

 

 

 

 

 

2

Использовался метод энергетического баланса [60].

 

 

 

 

4*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

51


где 0 можно определить по формуле

в = | / " - f [ 1 - 2n2 - -L (1 - 1 In") (sin А ) "

Если положим здесь я = 0, то получим выражение (1.28). Прибли­ женное уравнение амплитудных кривых получаем, используя равен­ ство (1.27):

QF

Рассмотрим колебания осциллятора с кубической характеристи­

кой,

описываемые уравнением х + их2 sgn х + ссх + fix3

= F cos Ы,

а >

О, Р g 0. Используя вторую1

формулу

(1.140),

определяем

 

a

j

_

 

 

U(A)=[2\(au+fiu3)e2nudu]2

 

=

 

- Т Г

{а {2пА

+

 

[А3

+ -ir Р'гА -

+

 

 

 

+

_зр_

 

(1.148)

 

 

 

2п2

 

Полагая

здесь 6 = 0 ,

получаем

выражение,

стоящее в левой

части

уравнения (1.146). Подставляя (1.148) в (1.139), находим уравнение амплитудно-частотной характеристики, которое оказывается весьма громоздким. Его можно упростить (за счет точности), если восполь­ зоваться равенством (1.32). Имеем

+ 4 ВЛ2 = ±

QF

1 + 6 = + 16 Л2 л2

 

(1.149)

Здесь 0, как известно 113], можно определить по формуле

e - V ^ [ i + - T - ( 3 i r - 2 п \ -

При п = 0 это выражение превращается в известную [13, 35, 57, 69] формулу Дуффинга

0

3

р

(1.150)

 

 

которая несколько уступает в точности выражению (1.33). Формулой (1.150) можно также пользоваться при условии 2гС- <^

<^3р7а. Если в (1.150) положить 0 = со и разрешить относительно

1 К этому же результату приводит использование "первой формулы (1.140).

62


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

Л = | / -

^ -

(1.151)

43

Р

 

Предполагается, что член в скобках формулы (1.150) мал по сравне­ нию с единицей.

Уравнение (1.145) кривой,

ограничивающей амплитудно-частот­

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

A ] / а +

_!_М2 -

n Q F

 

 

 

2 И

4Ллш2

Возводя это выражение в квадрат и используя формулу (1.150), приближенно имеем

Н1-&А*)--ШЯ--

(U52)

Для оценки точности полученных результатов сопоставим их с дан­ ными расчетов, проведенных одним из классических методов [60]. Исследования [18] приводят к следующим результатам.

Уравнение скелетной кривой

V1

Р "

(1.153)

Уравнение амплитудно-частотной характеристики

 

аЛ + 6Л3 =

1 -jJ-coM4 + Ао2 .

(1.154)

Уравнение ограничивающей

кривой

 

 

 

 

(1.155)

Амплитудные кривые, построенные по формулам (1.149) и (1.154)

при а =

20 сек-2, 6 = 4 см~2 • сект2, п = 0,1 см~1, F = 2,5 см X

X сек-2,

приведены на рис. 44. Как видим, совпадение кривых мож­

но признать удовлетворительным.

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

Сложные характеристики. Рассмотрим стационарные колебания

систем с перескоком. Они описываются частным

решением уравне­

ния

х + пх • sgn х — ах + 6л:3 = F cos at,

а > 0,

6 >• 0. Подстав­

ляя

выражения (1.42) и (1.43) в формулу

(1.139), получаем урав­

нение амплитудно-частотных характеристик для больших колеба­ ний, т. е. когда выполняется условие (1.39):

QF

(1.156)

63


Здесь 0, как известно [13, 32], определяется по формуле

 

9 = 0*(!

"9*

 

(1.157)

 

2УЦ

)

где 0„. = &6 находится

из выражения (1.41).

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

Уравнение (1.145) кривой, ограничивающей

ные характеристики,

в рассматриваемом случае, принимает вид

л ] / " 4 - р л а -

а

4Ллш2

'

Подставляя выражения (1.45) и (1.143) в формулу (1.139), полу­ чаем уравнение амплитудно-частотной характеристики для малых

а,сек

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

Рис. 45. Сопоставление аналитических

системы с кубической

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

и машинных

решений для амплитуд­

при турбулентном сопротивлении, по­

ной кривой системы с кубической

ха­

строенные по формулам:

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

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

со­

/ — (I . M9); 2 — (1.154);

3 — (I-151); 4 —

противлении.

 

 

(1.153); 5 — (I.155);ff —

(1.152).

 

 

 

колебаний, т. е. когда выполняется условие (1.38):

л 2 -

 

 

QF

К2р

. [ ш 2

( т л з - 1 ) + е Т + 1 6 4 ^ ш 4

 

 

 

(1.158) Здесь, аналогично изложенному выше, 0 можно определить по фор­

муле (1.157), где0„. =

0М находится из выражения (1.40).

На

рис. 46 приведены амплитудные кривые,

построенные при

а =

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

по формуле

(1.156) для

F =

1; 2; 4 см • сект2

и по формуле (1.158) для F =

0,25 см • сект2.

Кривые

для

малых

колебаний

практически

не

отличаются от

амплитудных

кривых тех

же параметров

в системах

без трения

(см. рис. 7) и в системах

с вязким трением

(см. рис. 29

и 30). Как

видим,

из рис. 46, совпадение аналитических и машинных данных,

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

можно

признать

хорошим.

 

54