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

Лекции по аналитической механике. Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
971 Кб
Скачать

Обобщенную силу Q1 найдем, определив работу веса на

возможном перемещении цилиндра, соответствующем возможному изменению δϕ угла ϕ. Эта работа равна

δ A = Pδ yC = Pδ[(R r) cosϕ] = −P(R r) sinϕδϕ.

 

Следова-

тельно, Q1 = −P(R r)sinϕ.

 

 

T

T = Q

 

Уравнение Лагранжа второго рода

d

в нашем

 

 

 

 

 

ɺ

∂ϕ

1

 

 

 

 

 

 

 

 

 

dt ∂ϕ

 

 

случае будет иметь вид

 

 

 

 

sinϕ = 0,

2

ɺɺ

 

 

ɺɺ

2

3m(R r)

ϕ = −mg(R r) sin ϕ,

или ϕ + k

 

где k 2 = 2g / 3(R r).

В случае малых колебаний sinϕ ≈ ϕ , и уравнение преобразуется к виду

ϕɺɺ+ k 2ϕ = 0.

Мы

получили дифференциальное уравнение свободных ко-

лебаний

с частотой k . Период этих колебаний T = 2π / k =

= 2π 3(R r) / 2g .

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

Задача 3.4. Два груза весом

P

и весом

P

, связанные меж-

 

1

 

2

 

ду собой пружиной жесткостью c,

могут двигаться без трения

по наклонной плоскости, составляющей угол α с горизонтом (рис. 3.4). Найти закон движения грузов, если в начальный момент времени их скорости равнялись нулю, а пружина имела длину 0.5l0 , где l0 – длина свободной пружины.

61

y

 

 

 

x1

x2

 

 

 

l

 

 

 

F1

 

 

P1

F2

 

 

 

 

 

 

P2

α

x

 

 

 

B

 

Рис. 3.4

 

 

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

К активным силам относятся внешние силы P1 и P2 – вес

грузов, а также внутренние силы упругости пружины F1

и F2

(см. рис. 3.4). Причем

F1

= F2

= cλ,

где

λ = l l0

удлинение

пружины.

 

степени свободы, s = 2.

 

 

 

Система имеет две

В

качестве

обобщенных координат возьмем координаты грузов

x1

и x2

в системе координат

Oxy.

Точку

O

выберем

так,

чтобы

в начальный момент времени t = 0, x1 (0) = 0. Тогда x2 (0) = 0.5l0 ,

а λ = x2 x1 l0 .

62

Кинетическая

энергия

системы равна

T = m v2

/ 2 +

 

 

 

 

 

 

 

1

1

 

2

ɺ2

/ 2

ɺ

2

/ 2 ,

где m1 = P1 / g, m2 = P2 / g

+ m2 v2

/ 2 = m1x1

+ m2 x

2

массы грузов.

Для нахождения обобщенных сил Q1 и Q2 дадим системе возможное перемещение, соответствующее возможным перемещениям грузов δ x1 и δ x2 . Тогда возможная работа активных сил на этих перемещениях

δAakt = (P sinα + F x + (P sinα − F x

2

=

1

 

1

1

2

 

2

 

 

= (P sinα + cλ)δx + (P sinα − cλ)δx

.

 

 

1

 

 

1

2

 

2

 

 

 

Следовательно,

Q = P sinα + cλ,

Q = P sinα − cλ.

 

1

1

 

 

2

2

 

 

 

Теперь составим уравнения Лагранжа второго рода (3.12):

m1 x1 = m1 g sinα + cλ,

m2 x2 = m2 g sinα − cλ.

(3.13)

ɺɺ

ɺɺ

 

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

Нам осталось теперь проинтегрировать уравнения (3.13)

при начальных данных

x1 (0) = 0, x2

(0) = 0.5l0

, x1

(0) = x2 (0) = 0.

 

 

 

ɺ

ɺ

Для этого умножим первое из уравнений (3.13) на m2 , а второе

на m1 и сложим их. Тогда получим уравнение m1m2 (x2

x1 ) =

 

ɺɺ

ɺɺ

= −c(m1 + m2 )λ . Используя

равенство λ = x2 x1 l0 ,

перепи-

шем его в виде λɺɺ+ k 2 λ = 0,

где k 2 = c(m1 + m2 ) / m1m2 .

Его ре-

шение мы знаем: λ = C1 sin kt + C2 cos kt. Из начальных условий находим, что C1 = 0,C2 = −0.5l0 . Следовательно, мы получили, что

 

x2 x1 = l0 − 0.5l0 coskt.

(3.14)

m1 x1

Второе соотношение получим, сложив уравнения (3.13):

+ m2 x2 = (m1 + m2 )g sinα. Перепишем последнее равенство

ɺɺ

ɺɺ

 

63

в виде

d 2

(m x + m x ) = (m + m )g sinα и дважды проинтегри-

dt 2

 

1

1

2

2

1

2

руем с учетом начальных данных. Тогда

(m1 x1 + m2 x2 ) = (m1 + m2 ) g sin α × t 2 / 2 + m2 l0 / 2. (3.15)

Из уравнений (3.14) и (3.15) находим закон движения грузов: x1 = g sin α × t 2 / 2 - 0.5m2 l0 [1 - cos kt] /(m1 + m2 ),

x2 = g sin α × t 2 / 2 + 0.5m1l0 [1 - cos kt] /(m1 + m2 ) + 0.5l0 .

Видим, что движение грузов слагается из двух движений – равномерно ускоренного движения с ускорением g sin α , обу-

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

3.4.Уравнения Лагранжа второго рода

вслучае консервативных систем

Напомним, что под консервативной системой мы понимаем механическую систему, в которой все действующие на нее си-

лы Fk консервативны,

 

т.е. имеют потенциал. При этом связь

между Fk

и потенциальной энергией

П системы дается фор-

мулами

 

 

П

 

 

 

 

 

П

 

 

П

 

 

 

F

= −

,

F

= −

,

F = −

.

(3.16)

 

kx

 

xk

 

ky

 

 

yk

kz

zk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Покажем,

что

в

 

этом случае

обобщенные силы Qj ,

j = 1,..., s,

связаны

с потенциальной

энергией системы

фор-

мулами

 

 

 

 

 

 

 

П

 

 

 

 

 

 

 

 

 

 

 

Q

 

= −

, j = 1,..., s,

 

 

(3.17)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

q j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т.е. в случае консервативной системы обобщенная сила, соответствующая обобщенной координате q j , равна взятой со

знаком «минус» частной производной от потенциальной энергии механической системы по этой координате.

64

Действительно,

 

 

 

 

 

rk

 

 

 

 

 

xk

 

 

 

yk

 

 

zk

 

 

Q

j

=

F ×

 

 

=

F

 

 

 

+ F

 

 

 

 

+ F

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

k

q j

 

 

kx

q j

 

ky

q j

 

kz

 

 

 

 

 

 

k

 

 

 

k

 

 

 

 

 

 

 

q j

 

 

 

 

 

 

П xk

 

П yk

 

П zk

 

 

 

 

 

 

 

-

 

 

 

 

 

+

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q j

yk q j

 

 

 

 

 

.

 

 

 

 

 

 

 

k xk

 

 

zk q j

 

 

 

Так как потенциальная энергия системы является функцией декартовых координат точек системы, П = П(xk , yk , zk ), а последние, в свою очередь, являются функциями ее обобщенных коорди-

нат,

xk = xk (q1,...,qs ,t), yk = yk (q1,...,qs ,t),

zk = zk (q1,..., qs ,t), то,

по

 

правилам

 

 

дифференцирования

 

сложной

функции,

 

П

 

П

 

x

 

 

 

П y

k

 

П z

k

 

 

 

 

 

 

=

 

k

+

 

 

+

 

 

 

 

,

откуда и следует форму-

 

q

 

x q

 

y q

 

z

 

 

 

 

 

j

 

j

 

 

j

 

k

 

q

 

 

 

 

 

 

 

 

k

k

 

 

 

 

 

 

k

 

 

 

 

 

 

 

j

 

 

 

ла (3.17).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Уравнения Лагранжа второго рода для консервативных сис-

тем тогда принимают вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

T

T

= −

 

П

,

j = 1,..., s .

(3.18)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt q

 

 

q j

 

 

 

 

q j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ɺ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если мы введем функцию Лагранжа L = T П и учтем,

что потенциальная энергия

П не зависит от обобщенных ско-

ростей, то уравнения (3.18) можно переписать в виде

 

 

d

L

L

= 0 , j = 1,..., s .

(3.19)

 

 

 

 

dt q

q j

 

 

 

ɺ

 

 

 

Уравнения (3.19) называются уравнениями Лагранжа второго рода для консервативных систем. Они широко использу-

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

Задача 3.5. Определить период свободных колебаний груза весом G , лежащего на двухопорной балке AB (рис. 3.5). Рас-

65

стояния от груза до опор балки равны a и b . Модуль упругости материала балки равен E , момент инерции поперечного сечения J . Весом балки пренебречь.

λst

 

F

A

0

B

 

 

x

 

G

 

 

x

a

 

b

Рис. 3.5

Решение. В качестве механической системы рассмотрим груз вместе с балкой. Связями для нее будут неподвижный шарнир A и подвижный шарнир B . Пренебрегая трением в шарнирах и считая шарнир B удерживающим, видим, что все условия для применимости уравнений Лагранжа второго рода выполнены.

Активными силами, действующими на систему, будут вес груза G и сила упругости балки F . Обе силы являются консервативными, поэтому уравнение (3.19) применимо.

Система имеет одну степень свободы. За обобщенную координату примем перемещение груза x , отсчитываемое от положения его равновесия.

Кинетическая энергия системы будет равна кинетической

энергии груза (весом балки пренебрегаем): T = mv2 / 2 = mxɺ2 / 2, m = G / g . Потенциальная энергия определяется как сумма по-

тенциальной энергии веса груза и потенциальной энергии силы упругости балки:

П = −Gx + [c

st

+ x) 2 / 2 − cλ2

/ 2] ,

(3.20)

 

st

 

 

66

где с – жесткость балки, а λst – ее статический прогиб. Так как в состоянии покоя G = cλst , то из (3.20) следует, что П = cx2 / 2.

Тогда функция Лагранжа L = T П = mxɺ2 / 2 − cx2 / 2, и уравнение (3.19) преобразуется к виду

ɺɺ

+ k

2

x = 0 ,

(3.21)

x

 

где k 2 = c / m. Уравнение (3.21) нам хорошо знакомо – это дифференциальное уравнение свободных колебаний груза с частотой k . Период этих колебаний TG = 2π / k. Воспользовавшись

еще формулой c = 3EJ (a + b) / a2b2 , известной из курса сопротивления материалов, окончательно получим, что

TG = 2π abG / 3EJg(a + b).

Задача 3.6 (колебания эллиптического маятника). Пол-

зун A массой mA может перемещаться вдоль гладких горизонтальных направляющих (рис. 3.6).

 

x

 

A

O

x

O1

 

B v2l φ

v2r

Рис. 3.6

67

Маятник B имеет массу mB , момент инерции J относительно точки подвеса O1 . Расстояние от точки подвеса до центра масс O2 равно h. Составить дифференциальные уравнения

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

Решение. В качестве механической системы рассмотрим оба тела вместе – ползун и маятник. Она будет подчинена действию двух связей (внешней связью для нее будут гладкие на-

правляющие, а внутренней – идеальный шарнир O1 ). Видим,

что обе связи таковы, что применять уравнения Лагранжа второго рода можно.

Система имеет две степени свободы, s = 2 . В качестве обобщенных координат возьмем перемещение x ползуна и угол ϕ маятника с вертикалью.

Ползун движется поступательно, поэтому его кинетическая энергия TA = mA v2A / 2 = mA xɺ2 / 2. Маятник совершает плоское движение, и для вычисления его кинетической энергии воспользуемся теоремой Кенига: TB = mB v22 / 2 + J2ω2 / 2. Здесь

 

ɺ

 

– скорость центра масс O2 маятника, а J2 – его

ω = ϕ , v2

 

момент инерции относительно центра масс O2 . Воспользуемся

еще кинематической формулой v2 = v2r

+ v2e

= v2r + v A . То-

гда

v2 = (v

A

+ ωh cosϕ)2 + (ωh sinϕ)2 . Учитывая теперь форму-

 

2

 

 

 

 

 

 

 

 

лу

 

 

 

 

2

,

получим,

что

ɺ2

/ 2 +

Гюйгенса J = J2 + mB h

TB = mB x

 

ɺ ɺ

 

ɺ

2

/ 2. Следовательно,

 

 

 

+ mB hxϕ cosϕ + Jϕ

 

 

 

 

T = TA + TB = (mA + mB )xɺ2 / 2 + mB hxɺϕɺcosϕ + Jϕɺ2 / 2.

68

Так как вес ползуна при его горизонтальном перемещении работы не совершает, то потенциальная энергия системы будет равна потенциальной энергии веса маятника, т.е.

П = −mB gh cosϕ.

Уравнения (3.19) для нашей задачи примут вид

d

 

[(mA

ɺ

ɺ

 

 

 

 

+ mB )x + mB hϕ cosϕ ] = 0,

dt

 

 

(3.22)

 

ɺɺ

 

ɺɺ

 

Jϕ

+ mB hx cosϕ + mB hg sin ϕ = 0.

Заметим, что первое из уравнений (3.22) мы могли бы получить, воспользовавшись теоремой о движении центра масс

системы в направлении оси x.

 

 

 

 

Интегрируя его, найдем,

с учетом

начальных

условий

ɺ

ɺ

что

(mA

ɺ

ɺ

или

x(0)

= ϕ(0) = 0,

+ mB )x + mBϕ cosϕ = 0,

d [(mA + mB )x + mB h sin ϕ ]= 0. Следовательно, dt

(mA + mB ) x + mB h sin ϕ = C.

Пусть при t = 0 маятник отклонен на угол ϕ(0) = ϕ0 . Воспользуемся произвольностью выбора начала O оси x так, что-

бы (mA + mB )x(0) + mB h sin ϕ0 = 0. Тогда C = 0,

и мы получаем

связь между величинами x и ϕ :

 

 

 

(mA + mB )x + mB hsinϕ = 0.

(3.23)

Второе дифференциальное уравнение (3.22) рассмотрим с

учетом

малости

угла отклонения маятника

от вертикали:

ϕ << 1,

sin ϕ ≈ ϕ,

cos ≈ 1 . Тогда из уравнения (3.23) будем иметь

x = −mB hsinϕ /(mA + mB ). Дифференцируя это выражение дваж-

69

ды по времени и подставляя значение ɺɺx во второе уравнение (3.22), получим, с учетом малости угла ϕ :

2

2

ɺɺ

(3.24)

 

 

(J h mB /(mA + mB ))ϕ + mB hgϕ = 0.

Так как

Jh2 mB2 /(mA + mB ) = J2 + mB h2 h2 mB2 /(mA + mB ) =

=J2 + mA mB h2 /(mA + mB ) > 0,

то

уравнение (3.24)

описывает гармонические

колебания

ϕ + k ϕ = 0 с частотой

k = mB hg /(J h

mB /(mA + mB )). Его ре-

ɺɺ

2

2

2

ϕ(0) = 0

 

шение при начальных данных ϕ(0) = ϕ0 ,

имеет вид

 

 

 

 

ɺ

 

ϕ = ϕ0 cos kt.

 

 

 

 

 

При движении системы центр тяжести маятника O2 (x2 , y2 )

движется по эллипсу:

 

 

 

 

 

x2 = x + hsinϕ = (mA /(mA + mB ))hsinϕ,

y2 = −h cosϕ,

поэтому маятник и называют эллиптическим.

Вопросы для самоконтроля

1.Какое уравнение называют общим уравнением динамики?

2.В чем состоит принцип ДаламбераЛагранжа? Почему его так называют?

3.При каких связях, действующих на механическую систему, можно применять принцип ДаламбераЛагранжа?

4.Как можно применять принцип ДаламбераЛагранжа, если среди связей есть связи с трением? А как его применять для определения реакции идеальной связи?

5.Можно ли применять принцип ДаламбераЛагранжа, если среди связей есть неудерживающие?

6.Какие величины называют обобщенными силами инерции?

70

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