Материал: Прочность, устойчивость, колебания. Т. 3

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

Области неустойчивости для систем уравнение

361

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

 

0 = ^ - (/1= 1, 2,.

. );

(38)

здесь

— частоты собственных колебаний системы, загруженной ста­

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

 

 

 

 

\Е — аА — Й2С | = 0.

 

 

 

 

 

 

(39)

Кроме того, уравнение (37) позволяет

выделить две группы

«подо­

зрительных»

частот. Первую

группу

частот определяют

по

формуле

 

 

 

д _ _ 0 / + Оа

(/1 =

1, 2,.

!ф к),

 

 

 

 

(40)

 

 

 

 

п

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

вторую группу частот

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0 / —

(п= 1, 2 .

 

 

 

 

 

 

(41)

 

 

 

п

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из

более

строгого

анализа

выражений

[14,

20,

21]

следует,

что

в случае гамильтоновых систем * области неустойчивости

 

образуются

лишь вблизи

частот, определяемых

по формуле (40). Если

система

является негамильтоновой, то возможны области

неустойчивости,

располагающиеся вблизи второй группы частот.

 

 

 

 

 

 

 

Области

неустойчивости, лежащие вблизи частот, соответствующих

формуле (38), называют основными областями

[5];

остальные

обла­

сти — комбинационные. Последнее название подчеркивает

тот

факт,

что колебания внутри этих областей осуществляются главным

образом

за счет взаимодействия какой-либо пары форм

колебаний.

Это

непо­

средственно следует из формул (40) и (41),

в которые

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

форме входят две собственные частоты й ; и

Аналогично можно го­

ворить об основных и комбинационных параметрических

 

резонансах.

 

Метод построения границ основных областей неустойчивости [7].

Границы

областей,

которым

в формуле (38) отвечают значения

п =

=

1, 3,

 

определяют из

уравнения

 

 

 

 

 

 

 

 

Е - а А ± 1 р я --\-№ С

 

 

 

 

О

 

 

 

 

 

 

 

5 -1»

г - < м — | е * с

-

4 -№

 

 

:

=о.

 

 

 

О

 

В

 

Е - а А - ^ &

‘С'.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

'

(42)

* Механическую систему с конечным числом степеней свободы называю! гамильтоновой, если уравнения ее движения могут быть представлены в кано­ нической форме Гамильтона. К нсгамильтоновым системам приходим, напри­

мер, рассматривая колебания стержня, нагруженного следящими силами [8]

362 Параметрические колебания упругих систем

Определитель, стоящий в левой части, следует понимать в том смысле, что на месте каждого его элемента стоит матрица порядка п, где п — число членов ряда (11). Например, удерживая в определителе (42) одни лишь выписанные члены, получим определитель порядка 3п. Границы главной области неустойчивости (п = 1) в первом прибли­ жении могут быть найдены из уравнения 17]

|е аА ± 4-РД---^-еас|=0.

(43)

Границы областей, которым в формуле (38) отвечают значения п —

= 2, 4, . . ., находят из уравнений, аналогичных выражению

(42).

Эти

уравнения неудобны для аналитических вычислений, поскольку

их

решение требует развертывания определителей высокого порядка

и отыскания корней алгебраических уравнений высоких степеней. Эти операции, однако, не представляют затруднений для электронных цифровых машин.

Применение метода малого параметра. Для аналитических вычисле­

ний более удобны формулы,

основанные на методе малого параметра

[21,

23].

Преобразуем

уравнение

(36)

к

главным

осям

матрицы

С ^ ( Е - а А ) :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

Р Я Л

СЮ 0/ = 0;

 

 

 

(44)

здесь / г — преобразованный вектор / ;

матрица С^ 1 имеет вид

 

 

 

 

 

 

0

 

 

0

 

 

 

 

 

 

 

 

0

&\

 

0

 

 

 

 

 

 

 

 

0

0

 

 

 

_

 

 

 

 

 

 

 

-

 

 

 

 

 

 

 

а матрицу И вычисляют применением преобразования

подобия к ма*

трице С_1Д. Иначе,

И —

 

 

 

 

 

 

(45)

 

 

 

 

 

 

 

 

 

где V — матрица, составленная по столбцам из собственных векторов

матрицы

С-1 аА).

 

1. Если

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

то

Ограничимся случаем п =

границы главных областей неустойчивости могут быть найдены

по

формуле

[28]

 

 

 

 

 

 

 

 

 

 

 

 

 

0 .= 2 О / ±

- |- .- |Щ ,

 

 

 

(46)

где

Л/д— элементы матрицы

Н. Границы

комбинационных

областей

определяют по формуле

 

 

 

 

 

 

 

 

 

 

 

 

0, =

П/ +

а * ±

А

 

]

/ » .

 

 

 

(47)

Если

диагональные

элементы матрицы

Нвелики

по

сравнению

с побочными элементами, то

ширина

комбинационных

областей

не­

устойчивости будет мала по сравнению с шириной главных

областей.

Влияние демпфирования на области неустойчивости

363

Примером, где это условие не выполнено, может служить задача об устойчивости плоской формы изгиба (см. 354—355). Пусть в уравне­ ниях (20) и (22) М (0 = рсоз 0/. Тогда уравнение (22) имеет вид (44), где С0 = С,

•Ч 9 1-а>1

И = С”1А =

Ык

Ширина главных областей имеет в этом случае порядок р2. Границы комбинационной области согласно формуле (47) определяют как

где Р* — критическое зи чеиие параметра для соответствующе;’ ческой задачи:

ВЛИЯНИЕ ДЕМПФИРОВАНИЯ НА ГРАНИЦЫ ОБЛАСТЕЙ НЕУСТОЙЧИВОСТИ

Уравнения для границ областей неустойчивости. Ограничимся про­ стейшим случаем уравнения Матье (26) с диссипативным членом

^ 1

+ 2 * 4 - + П2 (1 — со5 00 = 0;

(48)

л 2

+ г а

 

здесь е — коэффициент демпфирования. Области неустойчивости попрежнему располагаются вблизи частот, определяемых по формуле (2/). При п = 1,3, границы областей неустойчивости находят из урав­ нения [7]

.

902

 

ц

0

-

А

. Л

 

4Й2

 

 

 

 

я

— И

, ,

02

Д .

о

 

0

 

 

я ' 2Й

 

0

 

_д___ о_

,

02

 

= 0. (49)

 

_ р

 

л

20

1 11

4ЕГ-

А

39

 

0

- И

 

 

902

я ’

20

 

 

 

 

 

 

 

 

402

 

 

 

 

 

 

 

364

Параметрические колебания упругих систем

При п =

2, 4,

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

 

 

— Р

 

0

0

А

20

1

 

л ’

й

 

 

 

- Р

 

1-----

 

0

до

0

 

 

 

лС2

 

 

 

12*

 

 

 

0

 

 

 

1

— Р

0

 

0

 

Д

0

- 2р

1 - ^ 1

— Р

 

 

л

Ф

 

 

 

 

й 2

 

 

Д

20

0

 

0

— р

1 _ 4 2 1

л

П

 

 

 

 

 

о*

Критические значения коэффициента возбуждения. Наиболее суще­ ственный факт состоит в том, что уравнение (48) имеет возрастающее решение лишь при достаточно больших значениях параметра р. Так, для главной области неустойчивости из уравнения (49) получаем

формулу

О*

’2а] /

, ±

/ »

4М У

 

 

 

 

 

(51)

где

Д — декремент

затухания

собствен-

пых

колебаний,

 

 

 

 

 

Д = 2 л е

 

Из

формулы

(51)

видно, что вблизи

0= й

 

неустойчивость возможна лишь при

р >

р*, где р+ — критическое

значение

коэффициента возбуждения,

 

Д

Вообще можно показать [7], что для я-й области параметрического возбуждения критическое значение р* по порядку величины равно

 

_1_

 

^ ~ (4 )

л< = | ’ 2>-

>•

Следовательно, демпфирование подавляет побочные параметрические резонансы тем сильнее, чем больше число /г. Это показано на рис. 9, где построены области неустойчивости для случая, когда Д = 0,1.

Влияние перемещений в невозмущенном состоянии

365

Влияние демпфирования в общем случае. Если параметрические ко­ лебания описываются системой дифференциальных уравнений с периоди­ ческими коэффициентами, то расчет областей неустойчивости при на­ личии демпфирования усложняется. Некоторые приближенные методы указаны в книге [71. Метод малого параметра дает для определения границ комбинационных областей следующую приближенную фор­ мулу [28]:

О = й /

±

6/7 +

екк

Е А А /

 

(52)

 

 

2

 

 

 

 

здесь ^ — элементы

матрицы

демпфирования

после

при­

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

 

 

 

 

 

 

 

 

+ ( Е - а Д

- Р с°8 9 ^ ) / = 0

 

 

к главным осям матрицы

С-1 (Е — аА). Из формулы (52) видно,

что

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

ВЛИЯНИЕ ПЕРЕМЕЩЕНИЙ В НЕВОЗМУЩЕННОМ СОСТОЯНИИ

Постановка задачи. Рассмотренные выше задачи параметрических колебаний можно трактовать как задачи об устойчивости некоторых

режимов

установившихся вынужденных

 

колебаний.

Поясним

это

на примере задач, показанных на рис. 1.

В случае,

показанном

на

рис. 1, а,

роль невозмущенного движения

играют продольные колеба­

ния стержня, в случае рис. 1, 6 — радиальные колебания кольца, в случае 1, в — колебания пластинки в своей плоскости и т. д. Однако весь предыдущий анализ базировался на предположении, что переме­ щения в невозмущенном состоянии пренебрежимо малы. Рассмотрим уточненную постановку задачи для случая упругого стержня, сжимае­ мого периодической продольной силон (рис. 3).

Пусть и (х, I) — продольное перемещение точек; V(*, 0 — попереч­ ное перемещение точек, принадлежащих оси стержня; ЕР — жесткость сечения при растяжении-сжатии. С учетом наиболее существенных не­

линейных

членов

уравнения

совместных

продольных и поперечных

колебаний

имеют

вид

[6, 7,

 

13]

 

 

 

 

 

 

- ^ 4

 

г + " ‘ 4 т г = 0;

 

 

 

 

 

дх

1

д^^

 

 

 

 

д1и

сс

 

д

(

дь \

д2у

(53)

 

Е* ~ д ^ ~ ЕР дх ( 8х

дх )

+ '" д(2 ~ °

 

 

Через гх обозн чена продольная деформация оси стержня, т. е.

ди . 1 / до \ а

Источник: https://studfile.net/preview/19992756/