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

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

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

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

Добавлен: 22.10.2024

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

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

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

При указанных значениях р, коэффициент взаимного влияния гар­

моник,

представляющий собой отношение вторых максимумов,

kzl

1.

Из

приведенных данных и рис. 78 можно заключить следующее.

1.Влияние высшей гармоники на основную весьма мало, и при определении второго максимума амплитудной кривой им практи­ чески можно пренебречь.

2.Влияние основной гармоники на высшую заметно увеличи­ вает первый максимум амплитудной кривой, и им пренебрегать нельзя.

3. Кратность и частот воз­ буждения практически не влия­ ет на взаимодействие гармоник.

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

Сначала следует построить скелетные кривые 0 (А) и — ,

Т А Б Л И Ц А 2

к 1!

V-

Ft

=

Ft=.

 

=

0,2

 

 

= 0,1

см • сек

*

см • сек

2

 

1,09

 

 

 

1,07

3

 

1,07

 

 

 

1,05

4

 

1,07

 

 

 

1,05

5

 

1,07

 

 

 

1,05

П р и м е ч а н и е .

 

Для

исех

значений

ftal =

1.

 

 

 

 

 

 

затем по формуле (1.101) найти максимумы амплитудных кривых монохроматических движений (возбуждение одной гармоникой). Полученный максимум для основного тона возбуждения может быть принят в качестве второго максимума искомой амплитудной кривой и отложен на скелетной кривой 0 {А). Полученный максимум для высшего тона следует умножить на коэффициент k14 и отложить

6

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

V-

Для частоты со основного тона возбуждения, близкой или пре­ вышающей частоту 0 свободных колебаний линейной системы, в качестве искомой амплитудной кривой может быть принята ам­ плитудно-частотная характеристика монохроматического движения при возбуждении с основной частотой. В зоне минимума и при малых частотах возбуждения искомая амплитудная кривая может быть получена как сумма амплитудно-частотных характеристик моно­ хроматических движений (штриховые линии с крестиками на рис. 78).

Аналогичные исследования были проведены для систем с пере­ скоком. На рис. 79 представлены амплитудно-частотные характерис­ тики, вычисленные 1 по формулам (11.42) с помощью ЭЦВМ «Про-

мшь» для

систем с2

перескоком с параметрами- 1

а =

1

сект2;

В =

см~

сект ;

п = 0,05 се/с ; FL = F2

=

F,

см

сект*

при различных значениях р,. Точками показаны результаты 2

реше­

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

Ь2 Результаты получены Н. Г. Новиковой и М. Р. Суровицкой.

123


На основании табл. 2

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

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

характеристиками монохроматических

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

Турбулентное сопротивление. Рассмотрим стационарные коле­

бания нелинейной системы, описываемые уравнением

 

 

х + ях2 • sgn х + R (х) = Fi cos (at -f- vx ) + F2 cos (pat - f v2 ).

(11.45)

Используя метод переменного масштаба (см. § 3 гл. I), преобразуем

нелинейное уравнение (11.45) в линейное:

 

 

2"(е) +

г(е) =

 

 

Fx cos ^ - | - е +

Vjj

+ F2 cos ^ - | - e + vs

.

(11.46)

Сопоставление уравнений (11.32)

и (11.46) показывает,

что для

частных решений уравнения (11.46) можно воспользоваться выра­

жениями

(11.33), заменив

показатель

степени

на

т|,

т. е.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

гх

= ахе

" ( т - х )

 

 

/с о

,

 

\

 

 

 

 

 

 

 

 

 

 

 

 

> cos

I -g- е + vx — рг J;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(П.47)

 

 

 

 

 

 

г2

 

 

 

 

* ( т ~ т )

 

 

/

со .

v2

\

 

 

 

 

 

 

 

 

= й,е

1

 

'

cos Ifi - у е +

— р 2 1 .

 

 

Здесь в соответствии с формулами (1.136) обозначено

 

 

 

 

 

 

 

 

 

 

В]\

 

 

 

 

 

 

.

 

 

,

_

2feco2feco

.

Q

l

~

/(ft

2

- co

2

+

6

2

)

2

+

2

 

2

'

 

 

g P

l _

2

- со2 2 +~82 ~'

 

 

 

 

 

 

4A o)

 

 

 

 

 

ft

-co +>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(11.48)

о2

 

=

л

-

ц2

со2

 

 

 

 

 

 

 

 

;

 

'

tg р2

=

^—

6а

 

 

 

]/"(ft3

+ 02 )2 + 4ft2|A2co2

 

 

 

& — И2 ®2 +

Поскольку

 

уравнение

(11.46)

 

линейное,

то, суммируя

частные

решения

(11.47),

имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

г = е ( е

) ах c o s ^ e + vx + < 3 2 c o s ^ [ x - | - 8 + v 2 — p 2 j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(11.49)

Возвращаясь

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

заменами

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

/* (•*) = е * р (пх • sgn х) х cos (at + vx

+

+ a2 cos (pat + v2 — p2 )].

(11.50)

125


Далее рассмотрим

частный

случай, когда

vx да рг

и v2

да р2 .

В этом случае решение упрощается к виду

 

 

 

 

 

/* (х) = е х р (пх • sgn х) cos со/ + а2

cos ixot).

 

(II.51)

Заметим, что решения

(П.5) и (11.51) аналогичны. Следователь­

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

выра­

 

 

 

 

жением

(11.12). При

этом обе

 

 

 

 

гармоники одновременно прини­

 

 

 

 

мают

максимальные

значения

 

 

 

 

cosco^ = 1 и

cosfj-cof

=

1.

Под­

 

 

 

 

ставляя

эти

значения,

а

также

 

 

 

 

х = А и sgn х — 1

в

решение

 

 

 

 

(11.51), получаем выражение для

 

 

 

 

амплитудной

кривой

 

 

 

 

 

 

 

/*И) =

е ^ ( а 1 +

а2 ).

(11.52)

 

 

 

 

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

 

 

 

 

теристики жесткой системы при тур­

0,5

1,0

1,5

о.сен

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

нии частот возбуждения р. = 2.

Здесь амплитудная функция [А) определяется по формулам (1.140). В § 3, 7 гл. I и § 3 данной главы показано, что можно пользо­ ваться приближенным выражением (1.143). Подставляя его, а

также формулы (1.133) и (11.48) в равенство (11.52), имеем

 

 

 

6F,

 

 

 

+

] / ] с о * (

 

Л2

 

2

Л2

я^со-

4

1ГЗ

 

 

л-

 

 

 

 

 

+ "|/"|ц я ш а ^4 - 1 1

 

 

Л2

 

(11.53)

л 2 — 1 j + 92

 

 

2 "

 

 

 

 

 

 

 

 

 

В качестве примера

 

рассмотрим

кубическую

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

(1.31). В соответствии с формулами (1.32) и (11.53) для амплитудной характеристики получаем выражение

+

9^3

+

 

 

(11.54)

" | / | р А о 2 ^ 4 - ^ - я 2 — l j + 0а + 16 —— я2 р.4 со4

 

На рис.-80 изображены амплитудно-частотные характеристики, вычисленные1 по формуле (11.54) на ЭЦВМ «Наири-С» для систе­ мы с параметрами а = 1 сек*2; J3 = 1 см~2 - сект-2; п = 0,1 см.—1. Точками показаны решения уравнения (11.45), полученные2 на ЭЦВМ «Урал-3». Как видим, совпадение аналитических и машин­ ных решений можно признать удовлетворительным.

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

126


§ 3. Несимметричное возбуждение несимметричных систем

Рассмотрим часто встречающийся случай, когда воз­ мущение имеет постоянную составляющую, т. е.

 

F

if) =

Fq

+

F i

cos

со* + Fa

cos

LICO*;

Ц

^

 

F (*) =

F0

+

Ft

sin

со* -f- Fa

sin

LICO*.

 

 

Как

показано ранее (см. § 7 гл. I), наличие

постоянной составляю­

щей

возбуждения

симметричной

системы

переводит ее в

несим­

метричную систему. Поэтому в данном параграфе

рассматриваются

несимметричные системы и, следовательно, частота свободных ко­ лебаний 0 должна определиться в соответствии с указаниями § 7 гл. I .

Системы без трения. Сначала рассмотрим случай несимметрич­ ной кубической характеристики системы без трения при четном

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

т. е. когда

колебания

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

х + б0

+ OLX 4- ух2

+ (Зх3 = Fn

-J- Ft cos со* -f- F2 cos LICO*. (11.56)

Если воспользоваться заменой (1.278) и формулой (1.280), то урав­ нение (11.56) приводится к виду

У + а*У + $У3 = б* + F\ cos со* + F 2 cos р.со*,

(11.57)

где параметры а„. и б,,, определяются по формулам (1.282).

Путем

замен (1.3), аналогично (1.56)

приведем уравнение

(11.57)

к виду

г»(в) + г(е) = - g - (6 * +

^iCOS - J - е + F 2 c o s р

е

 

Частное решение этого уравнения, найденное аналогично (1.57), имеет вид

 

о г»

со

ш

 

Вг

г COS —— Б

8г COS Ц, —— 8

г ~ ~ е ~ '

е а — ш 2

' е а — ц2 ш2 *

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

# f/ a* +

P# - — + - § 2 — + ea -u.2 coa • ( 1 L 5 8 >

Поскольку член

в решении

(11.58) не зависит от времени *,

то условия максимума выражений

(П.5) и (11.58) ничем не отлича­

ются, т. е. они определяются формулой (11.12). Используя ее и под­

ставляя в решение

(11.58)

 

 

cos со* =

± 1; cos LICO* = ± 1; у =

Amax,

(11.59)

 

 

min

 

127


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

Т Л ™ п Va* + " 2 ~ Р Л 2 = "Т ~±

в*-а*

* б 3 - ^ с о 2

Здесь верхние знаки соответствуют амплитудной

кривой

для ЛШ ах,

нижние — для Л т т .

определяемого второй форму­

В случае нечетного возбуждения,

лой (11.55), совершенно аналогичные выкладки

свидетельствуют

о

том, что можно воспользоваться результатами §

1 гл. I I , добавляя

к

правой части решений уравнений и выражений

для амплитудно-

частотных характеристик постоянный член б^/б.

В частности, для случая р, = 2 в соответствии с формулой (11.28) амплитудно-частотная характеристика имеет вид

Значения входящих сюда величин приведены в § 1 данной главы.

Учет вязкого

трения.

Рассмотрим

влияние вязкого

трения

на стационарные

колебания,

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

 

х + 2пх + б0 + ax + ух2

-f- Br* = F0 + Fx cos (at -\- vx) +

 

+

F2

cos (\xai +

v2 ).

(11.60)

Если принять замену (1.278) и воспользоваться формулами (1.280) то уравнение (IT.60) приводится к виду

у + 2пу + а*у + Ру3 = 6*+ F t cos (at + + F2 cos (\iat + v2 ). (11.61)

Здесь использованы обозначения (1.282).

Путем замен (1.86), аналогично (1.108), приведем уравнение (11.61) к уравнению с постоянными коэффициентами

п

г»(е) + г(е) = - 1 - е Т Е б* + Fi cos J -J- e + vx J + F2 cos fji e + v2 j

Частное решение этого уравнения, найденное аналогично (1.109), имеет вид

 

 

 

 

п

 

 

 

г = -f- е + гх + г2 ,

где 2(, i =

1,2, определяются

по формулам (11.33), используя кото­

рые

имеем

 

 

 

2 =

6

+ ахcos

е +

v, P l J + а2cos J|i -J- s + v2 — p2 j j ,

128