Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Binder1[1].pdf
Скачиваний:
239
Добавлен:
21.03.2016
Размер:
7.46 Mб
Скачать

В уравнении Лагранжа

 

d

T

 

dT

Q = 0

(5.1)

 

 

 

q

 

dq

 

 

 

dt

 

 

 

 

 

присваиваем обозначения q ϕ,

q ω. Получаем уравнение в новых обо-

значениях

 

 

 

 

 

 

 

 

d T

 

dT

Q = 0

,

 

 

 

 

 

 

dϕ

 

 

 

ω

 

 

dt

 

 

 

 

 

где кинетическая энергия T должна быть найдена в виде функций фазовых переменных, T =T (ω,ϕ), а обобщенная сила Q – в виде функций фазовой ко-

ординаты Q =(ϕ) (в более сложных задачах имеем Q(ϕ,ω,t).

В данной теме двойная кинетическая энергия механизма приводится к следующей структуре:

2T = Jω2 ,

(5.2)

где по рисунку 5.1 определяем постоянный коэффициент, именуемый приведенным моментом инерции J механизма. Он является угловой мерой инертности механизма.

5.1.1 Потенциальная энергия и обобщенная сила механизма

Обобщенная сила Q представляет собой частное от деления мощности системы приложенных сил на обобщенную скорость ϕ =ω.

Q =

P

(5.3)

ω

 

 

Мощность в общем случае состоит из скалярных произведений сил на скорости их точек приложения, а также из моментов приложенных пар сил, скалярно умноженных на угловые скорости соответствующих звеньев:

 

 

P =

 

 

 

 

 

 

 

 

 

 

 

Fi

vi

 

 

j ωj

 

+ M

(5.4)

Здесь

 

 

 

 

 

 

vi

- скорость точки приложения силы Fi

, ω

j

- угловая скорость

тела, к которому приложена пара сил с моментом M j . Если все силы потен-

циальны и найдена их потенциальная энергия П(ϕ), то можно вместо формулы (5.4) применить другую формулу:

Q = −

П(ϕ)

(5.5)

ϕ

 

 

5.1.2.Математическая модель движения механизма

Вданных задачах коэффициент инерции J получается постоянным, не зависящим от угловой обобщенной координаты φ, т.е. J = const . Поэтому уравнение Лагранжа (5.1) упрощается, принимая вид

Jϕ Q =0

(5.6)

или

75

ϕ bQ =0 при b = J 1

(5.7)

т.е. угловое ускорение ϕ =ω прямо пропорционально обобщенной силе с не-

которым постоянным коэффициентом пропорциональности b , называемым

обратным приведенным моментом инерции механизма.

Таким образом, составление дифференциального уравнения движения второго порядка вида (5.6) сводится к вычислению коэффициента b и обоб-

щенной силы Q по формулам:

 

 

 

 

 

 

J =

2T

,

b =

ω2

, Q =

P

(5.8)

ω2

2T

ω

 

 

 

 

 

Математическая модель движения состоит из однородного дифференциального уравнения (ОДУ) вида (5.6) и начальных условий к нему, т.е. начального состояния механизма:

ω =ω0 , ϕ =ϕ0

(5.9)

В данном случае по условию задачи имеем ω0 =0, а также ϕ0 = 0, если угол ϕ отсчитывать от такого положения, при котором пружина была не деформирована.

Замечание 1 В более сложном случае, когда имеем переменный коэффициент инерции J = J (ϕ), в уравнении (5.6) добавляется нелинейное слагаемое:

Jϕ +

1 J

ϕ2 Q = 0,

ϕ bQ +

1 b

J

ϕ2

= 0, b =1/ J .

 

 

 

2 ϕ

ϕ

 

 

 

2

 

 

Замечание 2 Динамическое уравнение одностепенной системы со стационарными связями можно также получать из дифференциального урав-

нения энергии вида T = P , где P - мощность системы приложенных сил, точнее из уравнения T / q =Q .

5.2 Задание 5. Линейные колебания одностепенной механической системы

Кинематика

Изобразим механизм (рис. 5.1) в движении, т.е. в «положительной» фазе, когда обобщенная координата ϕ > 0 обобщенная скорость ϕ > 0 . Выразим

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

ниже.

Из рис. 5.1, а также условий нерастяжимости и непроскальзываемости троса, получаем кинематические выражения:

76

v =ωr2 ,

s = vdt =ϕr2 ,

ω3

= v3

=

v

=

r2

ω,

ϕ3

=

r2

ϕ

 

 

r3

r3

r3

 

 

 

r3

 

 

 

 

 

 

Двойная кинетическая энергия системы равна сумме кинетических энергий: поступательного движения груза 2T1 = m1v2 = m1r22ω2 ; вращательного

движения блока 2T2 = m2i22ω2 ; плоского движения колеса, состоящего из по-

ступательного движения вместе с центром масс и вращательного движения колеса вокруг центра масс:

2T3 = m3v32 + m3i32ω2 = m3 (r22 +i32 )ω2

Двойная кинетическая энергия механизма равна:

2T = 2(T1 +T2 +T3 ) =(m1r22 + m2i22 + m3 (r22 +i32 ))ω2

Приведенный момент инерции механизма есть множитель при квадрате

обобщенной (угловой) скорости, именно,J =

2T

 

 

ω2

 

J = m1r22 + m2i22 + m3 (r22 +i32 ),b = J 1

(5.10)

Из последнего выражения видим, что J = const , тогда уравнение Лагранжа будет иметь упрощенный вид (5.6)- (5.7).

Определение обобщенной силы методом мощности

Изобразим на рисунке 5.1 приложенные силы тяжести звеньев G1 , G2 , G3 , упругую силу пружины F =−csi , где s – деформация пружины, а также

реакции связей N2 , N3 , которые не имеют мощности (связи являются иде-

альными) и в дальнейшем не принимаются во внимание. Вычислим мощность системы сил:

P =G1 v1 +G2 O +G3 v3 + F v3 +(N2 O + N3 O),

отсюда

P =G1v +G3v cos π2 + β csv = (m1gr2 m3g sin βr2 cr2ϕ)ω = (mg cϕ)r2ω

при m = m1 m3 sin β .

Коэффициент Q = ωP в выражении мощности при обобщенной скорости

представляет собой обобщенную силу системы приложенных к механизму сил, т.е. меру эффективности сил.

Q =(mg cϕ)r2 .

(5.11)

Здесь мы не использовали понятие потенциальной энергии сил тяжести и упругости.

77

Определение обобщенной силы методом потенциальной энергии

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

работе всех действующих на механическую систему потенциальных сил при перемещении системы из произвольного положения q в «нулевое» положе-

ние q = 0, либо - в какое-либо другое выбранное положение q = q0 . Потенци-

альная энергия характеризует способность сил совершать работу на таком перемещении. Ее можно определять с точностью до аддитивной константы, т.е. можно отбрасывать или добавлять любую константу. Это объясняется тем, что в уравнение Лагранжа потенциальная энергия входит лишь в виде производной ∂Π / q , а при дифференцировании константы исчезают.

Переведем условно механизм из произвольного положения в «нулевое», в котором пружина не деформирована, и подсчитаем работу приложенных сил, которая принимается за потенциальную энергию системы сил в произвольном положении. Работа силы тяжести равна произведению этой силы на падение высоты, для груза работа отрицательна, так как груз подни-

мается. Потенциальная энергия пружины равна

cs

2 / 2 , где s = r ϕ – дефор-

мация пружины.

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

Потенциальная энергия всей системы равна

 

 

 

П = A

= −G s +G ssin β +

cs2

= −Gs

+

cs2

, G =G G sin β .

so

1

3

 

 

2

 

 

2

1 3

 

 

 

1 cr2ϕ2

 

 

 

или окончательно П = −mgrϕ +

. Дифференцируя это выражение,

 

 

2

2

2

 

 

 

 

 

 

находим выражение обобщенной силы

 

 

 

 

 

 

 

 

Q = −

П(ϕ)

= (mg cϕ)r

 

(0.1)

 

 

 

 

 

 

 

 

ϕ

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

.

Результат совпадает с выражением для Q , полученным методом мощно-

сти (5.11).

Определение положения равновесия и решение дифференциального уравнения движения системы

Подставим величину (5.11) в уравнение (5.7). Получим

ϕ +bcr2ϕ =(m1 m2 sin β)r2 g

или, вводя краткие обозначения,

 

ϕ +k2ϕ = h

(5.12),

где k = bcr2 , h =(m1 m2 sin β)r2 gb .

Будем считать, что m1 m2 sin β , т.е. h 0.

К линейному неоднородному уравнению второго порядка (5.12) присоединяем заданные начальные условия, характеризующие предысторию движения:

78

ст

k2

Выполним замену переменной

 

ϕ =ϕст +ψ ψ =ϕ ϕст

(5.13)

ω(0) =0

ϕ(0) =0

Получаем математическую модель движения в виде начальной задачи

Коши (5.12)-(5.13).

 

 

Статическое положение равновесия ϕ =ϕст. находим

как постоянное

решение неоднородного уравнения (5.12), а именно,

 

ϕ

= h

 

т.е. переносим начало отсчета координаты в статическое положение ϕст . Получим однородное уравнение для нового угла поворота ψ , который

отсчитывается от положения равновесия. Получаем начальную задачу для

угла поворота ψ :

 

=ψ (0)

 

 

 

 

 

 

ψ + k2ψ =0 при ψ

0

= −ϕ

ст.

, ψ

0

=ω =0

(5.14)

 

 

 

 

0

 

Решение задачи Коши для ОДУ вида (5.14) имеет вид кинематического

уравнения гармонических колебаний

 

 

 

 

 

 

ψ =ψ0 cos(kt)+ w0

sin(kt)= Acos(kt α) ,

 

 

 

k

 

 

 

 

 

 

где

2

A = ψ02 + wk0 - амплитуда колебаний

α=sgn(w0 )arccos ψ0 – начальная фаза колебаний

A

k– циклическая частота колебаний

τ= 2kπ - период колебаний, ν = 2kπ (Гц) – частота колебаний.

При заданных начальных условиях получаем

ψ = −kh2 coskt = kh2 cos(kt h)

Исходные данные для выполнения самостоятельной работы

Необходимо построить и решить динамическое уравнение колебаний одностепенной механической системы. Исходная схема механической системы для самостоятельного выполнения задания приведены в таблице 5.1. в соответствии с номером варианта. Жесткость с пружины, внутренние и наружные радиусы ri , Ri блоков, геометрический размер l стержней, массы mi всех

звеньев, радиусы инерции i блоков считать известными, не отмеченные номерами соединительные стержни считать невесомыми. Задачу решать в символьном виде.

79

 

Таблица 5.1

1

2

Дано: m1, m2, m3; R3 r3, i3, l; c.

Дано: m1, m2, m3; r2, i2,l; c.

3

4

Дано: m1, m2, m3, m4 ; r2, i2, r4, i4, l ; c. Дано: m1, m2, m3; r2, i2, R3, r3, i3; c.

5

6

Дано: m1, m2, m4; r2, i2, r4, i4; c. Дано: m1, m2, m3; r2, i2, R3, r3, i3; c.

80

 

Продолжение таблицы 5.1

7

8

Дано: m1, m2, m3; r2, i2, R3, r3, i3; c. Дано: m1, m2, m4; r2, i2, r4, i4; c.

9

10

Дано: m1, m2, m3; r2, i2, l; c

Дано: m1, m2; r2, i2, l; c.

11

12

Дано: m1, m2, m3; R3, r3, i3, l; c

Дано: m1, m2, m3, m4; r2, i2, r4, i4, l ; c.

81

 

Продолжение таблицы 5.1

 

 

13

14

Дано: m1, m2, m4; r4, i4, l ; c.

Дано: m1, m2, m3; r2, i2, l; c.

15

16

Дано: m1, m2, m3, m4; r2, i2, l; c.

Дано: m1, m2, m3; r2, i2, l; c.

17

18

Дано: m1, m2, m4; r2, i2, r4, i4; c. Дано: m2, m4; r2, i2, r4, i4; c.

82

 

Окончание таблицы 5.1

 

 

19

20

Дано: m1, m2, m3; R3, r3, i3, l ; c.

Дано: m1, m2, m3; r2, i2, l; c.

21

22

Дано: m1, m2, m4; r2, i2, r4, i4; c.

Дано: m1, m2, m3; r2, i2, r1, i1, l; c.

23

24

Дано: m1, m2, m3; r2, i2, r1, i1, l; c.

Дано: m1, m2, m3; r2, i2, r3, i3, l; c.

83

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]