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

ПриклТеорияКолеб

.pdf
Скачиваний:
26
Добавлен:
26.03.2015
Размер:
2.92 Mб
Скачать

 

 

 

 

 

 

24

α(2)2

= −

a11 k22

α(1)2 = −

a21

α(1)2 .

(19б)

a12

a22 k22

 

 

 

 

 

Принято использовать коэффициенты χi( j) форм колебаний, которые вводят-

ся

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

с

соотношением

αi( j) = χi( j)α1( j) . Ясно,

что

при этом

χ( j) =1, то есть в нашем случае χ(1) =1

и χ(2) =1. Два других коэффициента

1

 

 

 

 

 

 

1

1

 

 

 

χ(1)2 = −

a11 k12

 

,

χ(2)2

= −

a11 k22

.

 

(20а)

 

 

 

 

 

 

 

a12

 

 

 

a12

 

 

 

 

 

 

 

 

 

 

 

 

и частот

Подставив численные значения коэффициентов матрицы A

(13)

свободных колебаний k1,2 (17) в равенства (20а), получим:

 

 

χ(21)

= −0,65 ,

χ(22)

= 0,71.

 

 

 

 

 

(20б)

Первая и вторая формы колебаний построены на рис. 8. Из рисунка видно, что первое главное колебание с меньшей частотой k1 = 61,5 рад/с является симметричным (в том смысле, что колеса поворачиваются в разные стороны), причем амплитуда колебаний колеса 2 приблизительно в полтора раза меньше амплитуды колебаний колеса 1. Главное колебание с большей частотой k2 = 93,2 рад/с является антисимметричным, при этом амплитуда колебаний колеса 2 в 1,4 раза меньше амплитуды колебаний колеса 1.

Уравнения (18а), определяющие первое главное колебание, примут вид

ϕ(1)

= α(1)

sin(61,5t + δ ) ,

ϕ(1) = −0,65α(1)

sin(61,5t + δ ) .

(21а)

1

1

1

 

 

2

1

1

 

 

Уравнения (18б), определяющие второе главное колебание:

 

ϕ(2)

= α(2) sin(93,2t + δ

2

) ,

ϕ(1)

= 0,71α(2) sin(93,2t + δ

2

) .

(21б)

1

1

 

 

2

1

 

 

 

χ(1)

 

 

 

 

χ(2)

 

 

 

ϕ2

ϕ1

q

ϕ1

ϕ2

q

Первая форма колебаний

Вторая форма колебаний

Рис. 8

3.4.Закон движения механической системы

Общее решение системы уравнений движения (12б) определяется линейной суперпозицией собственных форм колебаний рассматриваемой механической системы (21а) и (21б):

 

 

 

 

 

 

 

 

25

 

 

 

 

 

 

 

 

 

 

 

 

ϕ

 

(t) = ϕ(1)

(2)

= α(1) sin(k t

) (2) sin(k

2

t

2

) ,

 

 

 

(22а)

1

1

1

1

 

1

1

1

 

 

 

 

 

 

 

 

 

 

ϕ

2

(t) = ϕ(1)

+ ϕ(2)

=μ α(1) sin(k t + δ ) + μ

2

α(2) sin(k

2

t + δ

2

) ,

 

(22б)

 

2

2

1

1

 

1

1

1

 

 

 

 

 

 

 

= χ(1) и

где введены обозначения для коэффициентов форм колебаний:

μ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

μ2 = χ(22) при условии, что χ1(1) =1 и χ1(2) =1.

Найдем из начальных условий постоянные интегрирования: амплитуды ко-

лебаний

α(1)

 

и α(2) ,

а также начальные фазы колебаний δ1 и δ2. Для этого

 

 

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

продифференцируем уравнения (22) по времени:

 

 

 

ϕ

(t) = k α(1) cos(k t + δ ) + k

2

α

(2) cos(k

2

t + δ

2

) ,

 

 

(23а)

1

1

1

 

1

1

1

 

 

 

 

 

 

 

 

 

ϕ

2

(t) k α(1) cos(k t + δ ) + μ

2

k

2

α(2) cos(k

2

t + δ

2

) .

(23б)

 

1 1

1

1

1

 

 

 

1

 

 

 

 

 

 

Уравнения (22) и (23) справедливы для любого момента времени t. Следовательно, при t = 0:

ϕ

 

= α(1) sin δ +α(2) sin δ

2

,

 

 

10

1

 

1

1

 

 

 

 

(24а)

 

 

= μ α(1) sin δ +μ

 

α(2) sin δ

 

ϕ

20

2

2

,

 

1

1

 

1

1

 

 

 

ϕ

 

= k α(1) cosδ +k

2

α(2) cosδ

2

,

 

 

 

 

 

 

10

1

1

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

(24б)

ϕ

 

= μ k α(1) cosδ +μ

 

k

 

α(2) cosδ

 

.

 

 

20

2

2

2

 

 

 

 

1 1

 

1

1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

Решая систему алгебраических уравнений (24), находим

tg δ = k

μ2ϕ10 − ϕ20

,

 

 

tg δ

2

= k

 

μ1ϕ10 − ϕ20

,

 

1 1 μ

2

ϕ − ϕ

20

 

 

 

 

 

 

 

 

 

 

2 μ ϕ − ϕ

20

 

 

 

 

 

10

 

 

 

 

 

 

 

 

 

 

 

 

 

1

10

 

 

(1)

=

(μ2ϕ10 − ϕ20 )2 + k12 (μ2ϕ10

− ϕ20 )2

 

(25)

α1

 

 

 

 

 

(μ

 

−μ

2

)2

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2)

 

(μ1ϕ10 − ϕ20 )2 + k22 (μ1ϕ10

− ϕ20 )2

 

 

α1

=

 

 

 

 

 

(μ

 

−μ

2

)2

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя в выражения (25) численные значения входящих в них величин, получаем

α(1)

= 0,00769 рад, α(2) = 0,00736 рад,

(26)

1

1

δ1 = arctg(8,65) = −1,46 рад, δ2 = arctg(14,31) =1,50 рад.

 

В итоге уравнения движения (22) механической системы примут вид

 

ϕ1(t) = 0,00769sin(61,5t 1,46) + 0,00736sin(93,2t +1,50) ,

(27а)

ϕ2 (t) = −0,00481sin(61,5t 1,46) + 0,00523sin(93,2t +1,50) .

(27б)

3.5.

Исследование вынужденных колебаний

 

С учетом обобщенных возмущающих сил

 

 

26

(28)

Qвозм = F sin ωt ,

где F — вектор-столбец амплитуд возмущающих сил, уравнения Лагранжа

(1) можно записать в следующем виде:

 

 

q+ A

q= f sin ωt ,

В нашем случае

M 0

 

F

=

 

.

 

 

0

 

 

 

 

1

(29)

f

= M

F .

 

 

 

(30)

Тогда система неоднородных дифференциальных уравнений (29) примет вид

ϕ1 + a11ϕ1 + a12ϕ2 = m111M 0 sin ωt,

(31)

ϕ2 + a21ϕ1 + a22ϕ2 = 0.

 

Решение системы уравнений (31) состоит из двух слагаемых, одно из которых является общим решением соответствующей однородной системы уравнений и для рассматриваемой механической системы определяется формулами (18). Второе слагаемое является частным решением неоднородной системы уравнений (29), которое ищем в виде гармонических колебаний с частотой ω:

ϕ1вын = h1 sin ωt,

(32)

ϕ2вын = h2 sin ωt.

 

Подстановка решения (32) в уравнения (31) приводит к системе алгебраических уравнений для амплитуд вынужденных колебаний h1 и h2:

(a

− ω2 )h + a h = m1M

0

,

 

11

1

12

2

11

 

(33)

(a

 

− ω2 )h + a

 

h = 0,

 

 

22

21

 

 

 

 

2

 

1

 

 

 

 

решение которой находим по методу Крамера:

h =

M 0

 

a22 − ω2

,

h

= −

M 0

 

a21

,

(34)

 

 

 

 

1

m11

 

2

m11

 

 

 

 

 

 

 

 

где определитель системы уравнений

= (a11 − ω2 )(a22 − ω2 ) a12a21 = (ω2 k12 )(ω2 k22 ) .

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

Тогда для рассматриваемой механической системы в соответствии с (32) для вынужденных колебаний запишем

ϕ1вын = H11M 0 sin ωt , ϕ2вын = H21M 0 sin ωt , (35)

где введены коэффициенты гармонического влияния:

 

 

 

 

 

 

 

 

 

 

27

H11

=

 

(a22 −ω2 )m111

 

 

,

 

(36а)

(ω2 k 2 )(ω2 k

2 )

 

 

 

 

1

 

2

 

 

 

 

 

 

 

a

21

m1

 

 

 

 

 

H21 = −

 

11

 

 

 

.

(36б)

(ω2 k

2 )(ω2 k

2 )

 

 

 

1

 

 

2

 

 

 

Из системы уравнений (31) найдем значения углов поворота дисков при статическом действии момента, равного амплитуде возмущающего воздействия (M = M0). При этом sinωt = 1, а ϕ1 = ϕ2 = 0 . В результате получаем

ϕ

=

 

a22

 

M

0

,

ϕ

2сста

= −

a21

ϕ

.

(37)

a a

 

 

 

1стат

 

22

a a

21

 

 

 

 

a

22

1сста

 

 

 

 

11

12

 

 

 

 

 

 

 

 

 

 

С помощью (37) перепишем уравнение (35) через статические значения углов поворота:

ϕ1вын = Φ1(ω)ϕ1стат sin ωt ,

ϕ2вын = Φ2 (ω)ϕ2стат sin ωt ,

(38)

где

(a11a22 a12a21)m111 a22 −ω2

 

 

Φ1(ω) =

 

 

 

 

 

 

 

 

 

 

,

 

(39а)

 

(ω2 k

2 )(ω2

k 2 )

 

 

 

 

a22

 

 

 

1

 

2

 

 

 

 

 

 

 

 

(a a

a a

21

)m1

 

 

 

 

Φ2 (ω) =

11 22

12

11

.

 

 

 

(39б)

 

 

 

 

 

 

 

 

 

 

(ω2 k 2 )(ω2 k 2 )

 

 

 

 

 

 

 

1

 

2

 

 

 

 

 

Учитывая численные значения элементов матрицы

 

A и матрицы коэффици-

ентов инерции M , для коэффициентов динамичности получим

K

=| Φ (ω) |=

 

1618 (6340 − ω2 )

,

(40а)

 

 

 

1дин

1

(ω2 3781)(ω2 8684)

 

 

 

 

 

 

 

 

K2дин =| Φ2 (ω) |=

 

10260700

 

 

.

(40б)

 

 

 

 

(ω2 3781)(ω2 8684)

 

 

 

 

 

 

 

Рассмотрим поведение Φ1(ω) и Φ2(ω) для трех интервалов изменения частоты возмущающего момента: а) 0 ≤ ω< k1; б) k1 < ω< k2 ; в) k2 < ω< ∞ (см.

рис. 9).

а) 0 ≤ ω< 61,5 .

В этом случае знаменатель в формулах (39) положителен. Следовательно, Φ1 > 0 и Φ2 > 0. С ростом частоты ω амплитуды колебаний обоих колес возрастают, причем их колебания происходят в одной фазе с изменением возмущающего момента. При ω = k1 в системе наступает первый резонанс и функции Φ1(ω) и Φ2(ω) претерпевают бесконечный разрыв. При действии момента сил резонансной частоты амплитуды колебаний возрастают пропорционально времени.

28

б) 61,5 < ω<93,2 .

В этом случае знаменатель в формулах (39) отрицателен. Следовательно, на рассматриваемом интервале частот Φ2 < 0, то есть, колебания колеса 2 про-

исходят в противофазе с изменением возмущающего момента.

Для определения знака функции Φ1(ω) нужно рассмотреть числитель форму-

лы (39а). Если k1 < ω< a22 ( 61,5 < ω< 79,3), то Φ1 < 0 и колебания колеса 1, так же как и колеса 2, происходят в противофазе с изменением возмущающего

Φ1

ω, с–1

ω = k1

 

ω = k2

 

 

 

Φ2

ω, с–1

ω = k1

 

ω = k2

 

 

 

29

Рис. 9

момента. Если же a22 < ω< k2 ( 79,3 < ω<93,2 ), то Φ1 > 0 и колебания колеса 1 происходят синфазно с изменением M. При ω= a22 = 79,3 с–1, Φ1 = 0.

Вэтом случае колесо 2 можно рассматривать как гаситель колебаний колеса

1.Явление отсутствия вынужденных колебаний колеса при действии на него

возмущающего момента частоты ω= a22 называется антирезонансом, а

само значение частоты – антирезонансной частотой. При ω = k2 = 93,2 в системе наступает второй резонанс.

в) 93,2 < ω< ∞.

В этом случае знаменатель в формулах (39) положителен. Так как k2 > a22 ,

то Φ1 < 0, а Φ2 > 0. С ростом частоты ω амплитуды вынужденных колебаний обоих колес уменьшаются, стремясь к нулю. Колебания колеса 1 происходят в противофазе, а колеса 2 — в одной фазе с изменением возмущающего момента.

30

ЛИТЕРАТУРА

1.Андронов А.А., Витт А.А., Хайкин С.Э. Теория колебаний.–М: Наука, 1981.

2.Бабаков И.М. Теория колебаний. – М: Наука, 1965.

3.Бидерман В.Л. Теория механических колебаний. – М.: Высш. шк., 1980.

4.Бицено К.Б., Граммель Р. Техническая динамика.–М.: ГИТТЛ, 1952.–Т.2.

5.Бутенин Н.В., Неймарк Ю.И., Фуфаев Н.А. Введение в теорию нелинейных колебаний. – М.: Наука, 1987.

6.Вибрации в технике. Справочник: В 6 т./ Ред. совет: В.Н. Челомей (предс.). – М.: Машиностроение, 1978 – 1982.

7.Вихренко В.С. Линейные колебания в системах с дискретными параметрами. – Минск: БТИ, 1993.

8.Вихренко В.С. Устойчивость и нелинейные колебания.–Минск: БТИ, 1993.

9.Вихренко В.С. Численные методы в теории линейных колебаний. – Минск: БГТУ, 1995.

10.Вихренко В.С. Колебания упругих балок и пластин.– Минск: БГТУ, 1996.

11.Пановко Я.Г. Введение в теорию механических колебаний. – М.: Наука, 1980.

12.Пановко Я.Г., Губанова И.И. Устойчивость и колебания упругих систем.

– М.: Наука, 1979.

13.Потсон Т., Стюарт Н. Теория катастроф и ее приложения.–М.: Мир, 1980.

14.Сборник заданий для курсовых работ по теоретической механике / Под ред. А.А. Яблонского. – М.: Высш. шк., 1985.

15.Тимошенко С.П. Колебания в инженерном деле. – М.: Физматгиз, 1963.

 

СОДЕРЖАНИЕ

 

Введение ...................................................................................................................

3

Программа курса......................................................................................................

4

Задача 1. Исследование устойчивости положений равновесия механической

системы и построение диаграммы устойчивости.................................................

5

Задача 2. Интегрирование уравнений движения системы с кусочно-линейной

характеристикой методом припасовывания (сшивки решений).......................

12

2.1

Составление дифференциального уравнения движения .............................

12

2.2

Интегрирование уравнения движения по методу припасовывания...........

14

Задача 3. Исследование свободных и вынужденных колебаний механической

системы с двумя степенями свободы...................................................................

21

3.1. Составление уравнений движения................................................................

23

3.2

Нахождение собственных частот и форм колебаний...................................

26

3.3

Формы главных колебаний.............................................................................

26

3.4

Закон движения механической системы.......................................................

28

3.5

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

29

Литература..............................................................................................................

37

31

32

33