- •СОДЕРЖАНИЕ
- •Введение
- •ГЛАВА 1. СТАТИКА ТВЕРДОГО ТЕЛА
- •ГЛАВА 4. ДИНАМИКА МЕХАНИЧЕСКОЙ СИСТЕМЫ
- •5.1.1 Потенциальная энергия и обобщенная сила механизма
- •5.1.2. Математическая модель движения механизма
- •9. Н.В. Бутенин, Я.Л. Лунц, Д.Р. Меркин. Курс теоретической механики. Статика и кинематика. Динамика: в 2 т. — СПб.: Лань, 2004 .— 730 с.
- •КАФЕДРА ТЕОРЕТИЧЕСКОЙ И ПРИКЛАДНОЙ МЕХАНИКИ
В уравнении Лагранжа
|
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