Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

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

.pdf
Скачиваний:
7
Добавлен:
13.11.2023
Размер:
24.79 Mб
Скачать

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

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 / до \ а

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

Решение уравнений (53) должно удовлетворять-граничным условиям

 

V

г

и =

0 при Хг= 0;

 

(54)

 

=

 

у = —

= 0 ;

Ерел =

Р0 + Р( соз0* при * = /.

 

Исследование

устойчивости. Уравнения (53)

граничные

усло­

вия (54) будут удовлетворены,

если положить

 

 

_______________

Решение (55)

описывает установившиеся про-

о

м

дольные колебания; это

решение будем

назы-

 

рис, ю

вать- невозмущенным-. Исследуем устойчивость

 

 

этого решения.

Полагая

в уравнениях

(53)

и = иа-\- ы,. и — о0+ Ъ, где и и о — малые возмущения, и линеари­ зируя уравнение (53),относительно этих возмущений, придем к следую­ щим уравнениям в вариациях:.

Уравнения (56) подробно исследованы в книге [7]. Важнейший ре­ зультат, относящийся к главной области параметрического резонанса, состоит в следующем: помимо резонанса вблизи соотношения 0 = 2Й, возможно возбуждение поперечных колебаний 0 = ш/ (со/ — частота продольных колебаний).

Физическое истолкование результата. Вблизи 0 = со/ имеет место резонанс продольных колебаний, вследствие чего резко возрастает динамическая продольная сила в стержне. Поэтому жесткость стержня по отношению к поперечным колебаниям вблизи 0 = со; является перио­ дической функцией времени с большой амплитудой изменения. Глав­ ные области параметрического возбуждения при отсутствии демпфиро­ вания показаны на рис. 10. При больших значениях коэффициента возбуждения р области сливаются.

Дальнейшие подробности можно найти в книге [7]. Там же рас­ смотрен в уточненной постановке ряд других задач (параметрические колебания криволинейных стержней, балок, рам, пластин и т. д.).

Учет нелинейнык факторов

367

УЧЕТ НЕЛИНЕЙНЫХ ФАКТОРОВ ПРИ ПАРАМЕТРИЧЕСКИХ КОЛЕБАНИЯХ

УПРУГИХ СИСТЕМ

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

Простейшая нелинейная задача. Нелинейные пара­ метрические колебания упругих систем рассматривались впервые в работах [4, 101 и позднее в статьях [6, 16, 25, 28, 29]. Приведем некоторые сведения, относя­ щиеся к простейшей задаче о колебаниях опертого стержня, сжатого периодической продольной силой. Пусть на подвижной опоре стержня (рлс. 11) имеется линейная упругая связь, препятствующая продольным перемещениям, продольный линейный демпфер и сосре­ доточенная масса. Полагая, что

V {х, 0 = / (0 51П - р

*(х,и

получим, что колебания системы приближенно описы­ ваются нелинейным дифференциальным уравнением

• ж + 2е ■§■+ 22 (* ~ 2* с05 е'> / +

+ 1> (/■

Й .

Щ = о

(57)

<И'

АР )

 

В уравнении (57) использованы те же обозначения, что и в урав­ нении (48); кроме того, в это уравнение входит нелинейная функция

* = V/3 + 2е„^ 1 + 2 ,/ [ / - $ + ( 4 ) 2] !

(58)

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

здесь у. ро и х — некоторые постоянные коэффициенты. Первый член в формуле (58) учитывает влияние продольной упругой связи, второй — влияние продольного демпфера, третий — влияние сосредоточенной массы, закрепленной на подвижном конце.

Графики зависимости амплитуды а0 от частоты возбуждения 9 схематически показаны на рис. 12, а—в в окрестности главного параметрического резонанса.

Рассмотрены три случая:

а) у > 0, е0 = к = 0; б) е0 > 0, у —х = 0, в) и > 0, у = еп — 0*

Штриховыми линиями обозначены неустойчивые решения (подробнее см. работы [7,10]). Там же дано решение различных нелинейных задач о параметрических колебаниях стержней, стержневых систем, пласти­ нок. Параметрические колебания упругих систем в нелинейной поста­

новке рассмотрены также

в работах [16,

25, 27].

 

ЛИТЕРАТУРА

 

 

1. Б е й л и н Е. А.,

Д ж а н е л и д з е

Г. Ю.

Обзор работ по дина­

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

математика и механика.

Т.6. Вып. 5. 1952.

2.Б е л я е в Н. М. Устойчивость призматических стержней под дей­

ствием переменных продольных сил. Сб. «Инженерные сооружения и строи­

тельная механика».

Изд-во

«Путь», 1924.

 

3. Б о д н с р

В. А. Устойчивость пластин под действием продольных

периодических сил.

Прикладная математика и механика. Т. 2,

№ 1, 1938.

1 Б о л о т и н

В. В.

О поперечных колебаниях стержней,

вызываемых

периодическими продольными силами. Сб. «Поперечные колебания и крити­ ческие скорости». Вып. I. Изд-во АН СССР, 1951.

5. Б о л о т н н

В. В.

О

параметрическом возбуждении поперечных

колебаний. Сб. «Поперечные

колебания и критические скорости». Вып. 2.

Изд-во АН СССР. 1953.

 

 

6. Б о л о т и н

В. В. О взаимодействии вынужденных и параметрически

возбуждаемых колебаний. «Изв.

АН СССР*. ОТН, № 4, 1954.

7.Б о л о т и н В. В. Динамическая устойчивость упругих систем. Гостехиздат, 1956.

8.Б о л о т и н В. В. Неконсервативные задачи теории упругой у с т о й ­ чивости. Физматгиз, 1961.

9.

В о л ь м и р А. С.

Устойчивость упругих систем. Физматгиз,

1963.

10.

Г о л ь де

н б л а т

И.

И.

Динамическая устойчивость сооружений.

Стройнздат, 1948.

и б л а т

И. И..

С и з о в

А. М.

Справочник по расчету

11.

Г о л ь д е

строительных конструкций

на

устойчивость и

колебания. Стройнздат,

1952.

12.

Д ж а н е л и д з е

Г.

Ю.,

Р а д ц и г

М. А.

Динамическая устой­

чивость

кольца под действием

нормальных периодических сил. П р и к л а д н а я

математика и механика. Т.

4. Вып. 5—6. 1940.

 

 

 

13.Д ж а и е л и д з с Г. 10. Устойчивость упругих систем при динами­

ческих нагрузках. Сб. «Проблемы устойчивости в строительной механике*. Стройнздат, 1965.

14.К р е й н М. Г. Основные положения теории Л — зон устойчивости канонической системы линейных дифференциальных уравнений с периоди­

ческими коэффициентами. Сб. памяти А. А. Андронова. Изд-во АН СССР,

1955.

 

15. М а к - Л а х л а и

Н. В. Теория и приложения функций

Матье.

Изд-во ИЛ, 1953.

 

Г.

В. О динамической устойчивости пологих

упру­

гих

16. М и ш е н к о в

оболочек. «Инженерный

журнал». 1961. Т. I. Вып. 2.

 

 

 

17. Ст р е тт М. Д.

Функции Ляме, Матье и родственные им в физике

и технике. ДНТВУ,

1935.

Динамическая устойчивость авиационных

кон­

 

18. Ч е л о м е й

В. И.

струкций. Изд-по «Аэрофлот», 1939.

устойчивости

упру­

гих

19. Ч е л о м е й

В. Н.

О возможности повышения

систем при помощи

вибраций. Доклады АН СССР.

Т. ПО,

3,

1956,

Литература

369

20. Я к у б о в и ч В. А. Замечание к некоторым работам

по системам

линейных дифференциальных уравнений с периодическими коэффициентами. Прикладная математика и механика. Т. 21. Вып. 5, 1957.

21. Я к у б о в и ч В. А. Метод малого параметра для канонических систем с периодическими коэффициентами. Прикладная математика и меха­ ника. Т. 23. Вып. I, 1959.

22.Е V а п • I V а п о V з к I К. М. Оп 1ке рагатсЫс гезропзс о! 5(гис1игез. Арр1. МесЬ. КеУ., уо1. 18, п. 9. 1965.

23.М е Ш е г Е. А11дсте1пе ТНеопо с!ег 8(аЫШа1 еггугилдепег 5сН\у1п-

дипдеп ектзИзсЬег Кбгрег. 1пдеп1еиг—А г с Ы у , Вб. 17, п. 6, 1949.

24.М е I 11 е г Е. Оупапйс ЬискИпд. НапбЬоок о! Елд1пеегшд МесЬа- п!сз (еб. Ьу \у. Р1йдде), Мс Ога\у-Н111, №\у-Уогк, 1962.

25.Р 1з г с г е к К. Рагатс1пс сотЫлаНоп гсзопапсе (о! 1Нс зесолб

к1пб) 1п попПпеаг зуз1ешз. Рогргалуу 1пгк., 1. 8, 1960.

26.

5 с к ш 1<11

О. ЙЬег сПе

В1еде$сЬумпдипдеп без де1епк!д де1адег1еп

ах1а! ри1з1сгелс1

Ье1аз1е1еп

31аЬс&.

Ма1ет. №скпсН1еп. Вб. 23, п. 2, 1961.

27.

3 с Н т

1 б 1

О. ПЬсг б1е

(2исгзс1г\у1пдипдеп вск^асН уогдскг0тт1ег

51аЬе Ье1 ехгсп1п5ск

апдгеИелбег

ри|51сгспбег

ЬапдзЬе1а$1ипд. Ма1еш. Иаск-

г1сЫел., Вб. 23, п. 4—5, 1961.

 

е г

Р.

1пз1аЬШШеп дебатрГ-

28.

5 с Нт

1 б 1

О.,

V/ е I б е п к я т т

1ег гкеоПпеагег 5сЬ\ушдипдеп. Ма111ет. №скг1сЫеп,

Вб. 23, N 4—5, 1961.

29.

Ш е I б е п На т т

е г Р.

ЭДсМПлеаге В1едезс1т1пдипдеп без ах1а1—

риЫегепб Ье1аз1е1сп

51аЬез. 1пдсп1еиг—АгсЫу,

Вб. 20, N 5, 1952.

Глава 7

КОЛЕБАНИЯ ПЛАСТИНОК

ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ КОЛЕБАНИЙ ПЛАСТИНОК

Уравнение изгиба пластинки в ортогональной системе координат для динамического случая. Пусть срединная поверхность изотропной пластинки постоянной толщины Н отнесена к ортогональной криво­ линейной системе координат хгх2 (рис. 1). Дифференциальное уравнение малых поперечных колеба­ ний в рамках гипотез Кирхго­

фа-Лява будет

ОДДш + рЛ т = <?; (1)

здесь ш—прогиб пластинки; 0 =

ЕН.3

— 12 (1 _ ^2) ~ Цилиндрическая

жесткость; р — плотность мате­ Рнс. 1 риала; ц (х1( х2, 0 — интенсив­ ность поперечной нагрузки; Е

модуль упругости; V — коэффициент Пуассона;

Н гН Л дхЛ Н ! дх,)

обобщенный оператор Лапласа; Яц Я 2 — параметры Ламе.

Моменты Мп, М22. М12 и

перерезывающие силы

выража­

ются через прогиб

пластины

ш следующими

формулами:

ч[

1 а

дш

дНх

дьи