Рассмотрим механическую систему из n материальных точек Mk (k=1,2,...,n), ограниченную l идеальными голономными связями. Ее положение в пространстве определим независимыми обобщенными координатами qi (i=1,2,...,s), число которых соответствует числу степеней свободы s механической системы. Обратимся к общему уравнению динамики (25) и запишем его в обобщенных координатах. Работа активных сил на возможных перемещениях, согласно (38) и (40), будет
δA=∑si=1(∑nk=1→Fk∂ˉrk∂qi)δqi=∑si=1Qi⋅δqi. (46)
Преобразуем оставшиеся слагаемые
n∑k=1mk−→Wk⋅δ→rk=n∑k=1mkd→Vkdts∑i=1∂→rk∂qiδqi=
=∑si=1(∑nk=1mkd→Vkdt⋅∂→rk∂qi)δqi . (47)
Выражение в скобках представим так
n∑k=1mkd→Vkdt⋅∂→rk∂qi=ddtn∑k=1mk→Vk∂→Vk∂˙qi−n∑k=1mk→Vk∂→Vk∂qi=
=ddt⋅∂∂˙qi∑nk=1mkV2k2−∂∂qi∑nk=1mkV2k2=ddt(∂T∂˙qi)−∂T∂qi , (48)
где T=∑nk=1mkV2k2 — кинетическая энергия системы, ˙qi=dqidt — i−ая обобщенная скорость.
Таким образом, общее уравнение динамики с учетом (48), (46) и (16) можно представить в виде
∑si=1(Qi−ddt∂T∂˙qi+∂T∂qi)δqi=0 .
Поскольку δqi независимы, равенство выполняется, когда
ddt(∂T∂˙qi)−∂T∂qi=Qi ; i=1,2,...,s. (49)
Мы получили систему обыкновенных дифференциальных уравнений второго порядка относительно s неизвестных функций qi(t), определяющих закон движения механической системы. Число уравнений равно числу степеней свободы; форма их записи не зависит от конкретного выбора обобщенных координат; они не содержат реакций идеальных голономных связей. Уравнения (49) называются уравнениями Лагранжа второго рода.
ПРИМЕР 10. Используя Уравнения Лагранжа второго рода, составить дифференциальное уравнение движения первого груза механической системы из примера 7.
РЕШЕНИЕ. На рис.13 изображена механическая система с действующими на нее внешними силами.
Выберем за обобщенную координату линейное перемещение первого груза. Уравнения кинематических связей останутся прежними:
V1=ω2⋅R или ˙s1=˙ϕ2⋅R;
ω2⋅r=ω32r3 или ˙ϕ2⋅r=2˙ϕ3r3;
V3=ω3⋅r3 или ˙s3=˙ϕ3⋅r3.
Составим выражение для кинетической энергии механической системы:
T=Tпост1+Tвр2+Tпост3+Tвр3=P1˙s212g+I2˙ϕ222+P3˙s232g+I3˙ϕ232 .
При записи учтено, что мгновенный центр скоростей диска, катящегося без скольжения, расположен в точке его соприкосновения с плоскостью.
Вынесем за скобки квадрат обобщенной скорости
T=˙s212(P1g+IR2+3P3r28gR2).
При получении результата учтено, что осевой момент инерции однородного диска I3=P3r232g.
Частная производная от кинетической энергии по обобщенной координате будет равна нулю, так как полученное выражение от обобщенной координаты не зависит.
Частная производная от кинетической энергии по обобщенной скорости будет
∂T∂˙s1=˙s1(P1g+IR2+3P3r28gR2) . Производная по времени от частной производной кинетической энергии по обобщенной скорости будет
ddt(∂T∂˙s1)=¨s1(P1g+IR2+3P3r28gR2) .
Для вычисления обобщенной силы дадим системе возможное перемещение. Составим сумму работ всех внешних сил на соответствующих возможных перемещениях и вынесем возможное перемещение по обобщенной координате s1 за скобки:
δA=P1δs1−Fуδs1−P3sinαδs3−Fсопрδs3==δs1[P1−c(s1+Δ)−P3rRsinα−b˙s1r2R2]
Выражение в квадратных скобках и есть обобщенная сила для возможного перемещения по выбранной обобщенной координате.
После подстановки в (49) выражений для частных производных и обобщенной силы, а так же учета условия равновесия (см. пример 17 из лекции 6), окончательно имеем ¨s1(P1g+IR2+3P3r28gR2)+˙s1br2R2+s1⋅c=0
Полученное уравнение совпадает с результатом решения примера 7.
ПРИМЕР 11 (задача 48.26 из [2]). Решим задачу из примера 8 (см. рисунок 14), используя уравнения Лагранжа второго рода.
РЕШЕНИЕ. Рассматриваемая механическая система обладает двумя степенями свободы. Выберем, как и при решении примера 8, в качестве независимых обобщенных координат горизонтальное перемещение x груза А и вертикальный подъем y груза В. Тогда, в силу нерастяжимости нити, опускание по вертикали z центра блока Д (и груза К) будет зависеть от выбранных обобщенных координат.
Уравнение кинематической связи были получены при решении примера 6:
VK=VA+VB2=˙x+˙y2=˙z .
Составим выражение для кинетической энергии системы
T=TA+TB+TK=MV2A2+M1V2K2+MV2B2==(M1/4+M)˙x22+M1˙x˙y/4+(M1/4+M)˙y22
и возьмем соответствующие производные:
∂T∂x=∂T∂y=0 ; ∂T∂˙x=(M1/4+M)˙x+M1˙y/4 ;
ddt(∂T∂˙x)=(M1/4+M)¨x+M1¨y/4 ; ∂T∂˙y=(M1/4+M)˙y+M1˙x/4 ;
ddt(∂T∂˙y)=(M1/4+M)¨y+M1¨x/4 ;
Выражения для обобщенных сил были получены при решении примера 6. Изменив индексы масс, (так как в задачах 46.21 и 48.26 из [5] они не совпадают) имеем:
Qx=δAxδx=−fMg+M1g/2 ;
Qy=δAyδy=−Mg+M1g/2 .
Теперь составим систему из двух дифференциальных уравнений второго порядка
(M1/4+M)¨x+M1¨y/4=−fMg+M1g/2 ;
M1¨x/4+(M1/4+M)¨y=−Mg+M1g/2 .
Сложим полученные уравнения и учтем результат дифференцирования уравнения кинематических связей ¨z=¨x+¨y2. Окончательно получим
¨z=gM1−M(1+f)2M1+M .
Очевидно, что груз К будет двигаться вниз, если числитель полученного выражения положителен, т.е. при M1>M(1+f) .
Уравнение Лагранжа второго рода для консервативных систем
В случае действия на систему потенциальных сил, обобщенные силы могут быть найдены как
Qi=−∂П∂qi , i=1,2,...,s (50)
где П=П(q1,q2,...,qs) — потенциальная энергия системы.
Тогда уравнения Лагранжа принимают вид
ddt(∂T∂˙qi)−∂T∂qi=−∂П∂qi, i=1,2,...,s (51)
Если на систему действуют силы, имеющие потенциал и неконсервативные силы, то получим уравнения Лагранжа в виде
ddt(∂T∂˙qi)−∂T∂qi=−∂П∂qi+QHKi , i=1,2,...,s (52)
Уравнения (51) можно преобразовать путем введения функции Лагранжа L=T−П, называемой кинетическим потенциалом.
Выразим кинетическую энергию через кинетический потенциал и возьмем соответствующие производные из (51). Тогда, с учетом зависимости потенциальной энергии только от обобщенных координат, получим:
T=L+П ; ∂T∂qi=∂L∂qi+∂П∂qi ; ∂T∂˙qi=∂L∂˙qi.
Окончательно имеем ddt(∂L∂˙qi)−∂L∂qi=0, i=1,2,...,s (53)
Уравнения (51) или (53) называются уравнениями Лагранжа для консервативных систем.
ПРИМЕР 12. Составить уравнения движения двойного математического маятника, изображенного на рис. 15.
РЕШЕНИЕ. Составим выражение для кинетической энергии двойного математического маятника:
T=T1+T2=P1V212g+P2V222g .
Скорость второй тяжелой точки →V2=→V1+→V21, где V1=˙ϕ1l1 — скорость переносного движения с первой тяжелой точкой, а V21=(˙ϕ2−˙ϕ1)l2— относительного вращения вокруг первой точки. Тогда V22=V21+V221+2V1V21cos(ϕ2−ϕ1).
Выбрав, как и ранее, в качестве независимых обобщенных координат углы с вертикалью ϕ1 и ϕ2, выражение для кинетической энергии можно привести к виду:
T=12g[(P1+P2)l21˙ϕ21+2P2l1l2˙ϕ1˙ϕ2cos(ϕ2−ϕ1)+P2l22˙ϕ22].
Вычислим соответствующие производные:
∂T∂ϕ1=1gP2l1l2˙ϕ1˙ϕ2sin(ϕ2−ϕ1); ∂T∂ϕ2=−1gP2l1l2˙ϕ1˙ϕ2sin(ϕ2−ϕ1) ;
∂T∂˙ϕ1=1g[(P1+P2)l21˙ϕ1+P2l1l2˙ϕ2cos(ϕ2−ϕ1)] ;
∂T∂˙ϕ2=1g[P2l1l2˙ϕ1cos(ϕ2−ϕ1)+P2l22˙ϕ2] ;
ddt(∂T∂˙ϕ1)=1g[(P1+P2)l21¨ϕ1+P2l1l2¨ϕ2cos(ϕ2−ϕ1)−−P2l1l2˙ϕ2(˙ϕ2−˙ϕ1)sin(ϕ2−ϕ1)]
ddt(∂T∂˙ϕ2)=1g[P2l1l2¨ϕ1cos(ϕ2−ϕ1)−−P2l1l2˙ϕ1(˙ϕ2−˙ϕ1)sin(ϕ2−ϕ1)+P2l22¨ϕ2]
Для системы, изображенной на рис.15, ранее (см. пример 9) было составлено выражение для потенциальной энергии и найдены обобщенные силы, соответствующие обобщенным координатам ϕ1 и ϕ2:
П=− P1l1cosϕ1−P2(l1cosϕ1+l2cosϕ2)=
=− (P1+P2)l1cosϕ1−P2l2cosϕ2 ;
Q1=−∂П∂ϕ1=−(P1+P2)l1sinϕ1 ;
Q2=−∂П∂ϕ2=−P2l2sinϕ2 .
Подставив полученные выражения в (51), получим систему из двух дифференциальных уравнений второго порядка, описывающих движение двойного математического маятника:
(P1+P2)l21¨ϕ1+P2l1l2¨ϕ2cos(ϕ2−ϕ1)+P2l1l2˙ϕ22sin(ϕ2−ϕ1)==−g(P1+P2)l1sinϕ1
P2l1l2¨ϕ1cos(ϕ2−ϕ1)+P2l1l2˙ϕ21sin(ϕ2−ϕ1)+P2l22¨ϕ2==−gP2l2sinϕ2
Циклические координаты, циклические интегралы
Обобщенные координаты, не входящие в явном виде в кинетический потенциал L, называются циклическими координатами.
Предположим, что среди s обобщенных координат q1;q2;...;qk;k<s являются циклическими.
Тогда, в соответствии с определением циклических координат, производные от кинетического потенциала по этим координатам равны нулю:
∂L∂qj=0 ; j=1,2,...,k.
В этом случае уравнения (53) принимают вид
ddt(∂L∂˙qj)=0 ; j=1,2,...,k.
Откуда ∂L∂˙qj=Cj=const; j=1,2,...,k. (54)
Равенства (54) называются циклическими интегралами.
Как видно, циклические координаты позволяют сразу записывать первые интегралы, что значительно упрощает решение уравнений Лагранжа. По существу уравнения (54) представляют собой соотношения между кинематическими характеристиками движения, а учет каждого такого соотношения позволяет понизить на единицу порядок системы дифференциальных уравнений Лагранжа. В итоге вместо решения системы дифференциальных уравнений порядка 2s при k циклических координатах требуется решать систему дифференциальных уравнений порядка 2s−k. Поэтому обобщенные координаты следует выбирать так, чтобы возможно большее их число было циклическими.
ПРИМЕР 13 (задача 48.28 из [5]). Призма А массы m скользит по гладкой боковой грани призмы В массы m1, образующей угол α с горизонтом (см. рис.16.а).
Определить ускорение призмы В. Трением между призмой В и горизонтальной плоскостью пренебречь.
РЕШЕНИЕ. Система обладает двумя степенями свободы, так как для остановки всех тел системы необходимо зафиксировать две координаты, например перемещение s призмы А относительно призмы В и горизонтальное перемещение x самой призмы В.
Именно эти перемещения выберем в качестве независимых обобщенных координат. Составляя выражение для кинетической энергии учтем, что призма А участвует в двух движениях – переносном с призмой В и движении относительно нее (см. рис.16.б).
Тогда
T=TB+TA=12m˙x2+12m1(˙x2+˙s2−2˙x˙scosα)=
=12(m1+m)˙x2+12m1˙s2−m1˙x˙scosα .
Составим выражение для потенциальной энергии (ее значение в начальном положении принято за нулевое):
П=−m1gssinα .
Подставим выражения энергий в функцию Лагранжа
L=T−П=12(m1+m)˙x2+12m1˙s2−m1˙x˙scosα+m1gssinα .
Выполнив ряд формальных действий, довольно легко получить систему из двух дифференциальных уравнений Лагранжа, описывающих движение тел механической системы
(m+m1)¨x−m1¨scosα=0 ;
¨s−¨xcosα=gsinα .
Решив систему уравнений, найдем ускорение призмы В:
¨x=mgsin2α2(m1+msin2α) .
Заметим, что обобщенная координата x является циклической, так как отсутствует в выражении кинетического потенциала. В таком случае мы могли без составления уравнений Лагранжа сразу записать выражение
∂L∂˙x=(m1+m)˙x−m1˙scosα=const,
совпадающее с первым интегралом первого дифференциального уравнения. Полученное соотношение показывает, что абсцисса центра масс механической системы, в зависимости от начальных условий, либо сохраняет свою величину, либо изменяется пропорционально времени.
Особенности применения уравнений Лагранжа второго рода для систем с неидеальными связями
Уравнения Лагранжа получены для механической системы с идеальными, удерживающими и голономными связями. Это не означает, что их нельзя использовать для систем с неидеальными неудерживающими связями.
Например, если связи, наложенные на механическую систему неидеальны, их неизвестные реакции будут входить в выражения обобщенных сил. Очевидно, что наличие дополнительных неизвестных требует составления соответствующего числа дополнительных уравнений.
Именно так и было сделано при решении примера 6, где сила сухого трения вычислялась предварительно как произведение коэффициента трения на силу нормального давления; в свою очередь, сила нормального давления определялась из равенства нулю проекции на вертикаль главного вектора сил, действующих на груз А , то есть F=fN=fMg.
В более сложных случаях в процессе движения реакции неидеальных связей могут изменяться. Для составления дополнительных уравнений обычно используют общие теоремы механики. Сказанное проиллюстрируем примером.
ПРИМЕР 14 (пример из [1], п.19.4). Ползун А линейки эллипсографа АВ скользит по шероховатой вертикальной стене, а ползун В – по гладкой горизонтальной поверхности (рис.17).
В этом случае связь, наложенная на конец А линейки эллипсографа — удерживающая неидеальная, а связь, наложенная на конец В – удерживающая идеальная. Линейка имеет одну степень свободы, ее положение может быть задано, например, углом ϕ с вертикалью. Требуется получить дифференциальное уравнение движения линейки эллипсографа.
Составим выражение для кинетической энергии линейки:
T=12mV2C+12IC˙ϕ2=16ml2˙ϕ2 .
При выводе учтено, что точка С – центр масс линейки эллипсографа; скорость центра масс V2C=˙x2C+˙y2C,
где xC=12lsinϕ; yC=12lcosϕ, а осевой момент инерции линейки при вращении вокруг центра масс IC=112ml2,
Составим выражение для потенциальной энергии силы веса, приняв горизонтальное положение линейки за нулевое:
П= 12mglcosϕ.
Непотенциальная сила трения FA=fNA, где f — коэффициент трения скольжения (предполагается, что ползун А движется вниз, поэтому на рисунке сила трения направлена вверх). Для составления уравнения Лагранжа нужно определить обобщенную силу, соответствующую углу поворота, с учетом непотенциальной силы трения FA :
δA=δAПОТ+δAНЕПОТ=δAПОТ−FAδyA .
Поскольку yA=lcosϕ, то δyA=−lsinϕ⋅δϕ.
В таком случае Qϕ=−∂П∂ϕ−FAlsinϕ=l(12mg−fNA)sinϕ .
Возьмем соответствующие производные от кинетической энергии, выражение для обобщенной силы и составим уравнение Лагранжа второго рода:
13ml¨ϕ=(12mg−fNA)sinϕ.
Это уравнение не может быть проинтегрировано, так как величина NA неизвестна. Для ее определения используем теорему о движении центра масс механической системы, спроецировав соответствующее векторное уравнение на ось абсцисс: m¨xC=NA.
Поскольку xC=12lsinϕ, то ¨xC=l2(¨ϕcosϕ−˙ϕ2sinϕ).
Отсюда NA=ml2(¨ϕcosϕ−˙ϕ2sinϕ).
Внесем это выражение в уравнение Лагранжа. После преобразований, получим:
(1+34fsin2ϕ)¨ϕ=32(gl+f˙ϕ2sinϕ)sinϕ.
Это уравнение может быть проинтегрировано приближенными методами, однако дальнейшее решение выходит за рамки настоящего пособия.
Особенности применения уравнений Лагранжа второго рода для систем с неудерживающими связями
Ели связи, наложенные на механическую систему неудерживающие, может иметь место ситуация, когда связь в некоторый момент перестает действовать (либо наоборот – накладывается). Силовые схемы на этапах движения с действующей связью и без нее, различны и, как следствие, различны на этих этапах дифференциальные уравнения, описывающие движение механической системы.
Характерным примером является спуск с горы лыжника, который может оторваться от поверхности горы вследствие резкого изменения ее профиля (трамплин) либо достаточно высокой скорости на выпуклости. В последнем случае точка отрыва от связи (начало этапа полета) определяется равенством нулю силы нормального давления лыжника. Конечные положение и скорость на этапе скольжения есть начальные положение и скорость для этапа полета.
Точкой окончания этапа полета является точка пересечения траектории полета и профиля поверхности горы. Если в этой точке касательные к траектории движения и к профилю поверхности окажутся достаточно близкими, лыжник приземлится «мягко» (иногда говорят, что он хорошо «вписался»), если нет – будет иметь место механический удар.
Таким образом, решение задачи о движении механической системы с неудерживающими связями может представлять собой ряд этапов, для каждого из которых составляется своя силовая схема и выводятся свои дифференциальные уравнения; при этом конечные кинематические характеристики предыдущего этапа являются начальными для последующего.
На протяжении рассматриваемого этапа осуществляется контроль величины параметра, определяющего переход к следующему этапу движения, на котором либо какие–то связи прекращают действовать либо начинают действовать дополнительные связи. Сказанное проиллюстрируем примером.
ПРИМЕР 15 (пример из [1], п.19.4). Стержень АВ скользит без трения по сторонам прямого угла (рис. 18).
До тех пор, пока стержень опирается на стороны угла, связи не нарушены и уравнение движения стержня может быть получено из уравнения движения линейки эллипсографа при f=0, т.е.
¨ϕ=32glsinϕ. Это уравнение может быть представлено в виде
d(˙ϕ22)=−3g2ld(cosϕ) и проинтегрировано при начальных условиях ϕ=ϕ0, ˙ϕ=0 (движение наклоненного стержня из состояния покоя).
Тогда первый интеграл будет ˙ϕ2=3gl(cosϕ0−cosϕ).
Однако сразу использовать эти уравнения нельзя, так как не известно, произойдет или нет отрыв стержня от вертикальной стенки. Если отрыв может иметь место, то надо определить угол ϕ1, соответствующий этому моменту. Очевидно, что сила NAнормального давления на стену в момент отрыва должна стать равной нулю, то есть
NA=ml2(¨ϕ1cosϕ1−˙ϕ21sinϕ1)=0.
Выразим ¨ϕ и ˙ϕ2 из дифференциального уравнения движения и его первого интеграла и подставим в выражение для силы нормального давления. Тогда
NA=34mgsinϕ1(3cosϕ1−2cosϕ0)=0,
Отсюда cosϕ1=23cosϕ0. Таким образом, полученное уравнение Лагранжа определяет движение не при всех ϕ, как это было в случае удерживающих связей, а только если ϕ0≤ϕ≤ϕ1.
С момента отрыва точки А от вертикальной стены стержень будет иметь не одну, а две степени свободы (по предположению движение происходит в вертикальной плоскости). Если нас интересует дальнейшее движение стержня, то, как уже говорилось выше, для второго этапа движения следует изобразить стержень АВ с действующими на него силами (см. рис.18 при отсутствии реакции →NA), получить каким-либо методом дифференциальные уравнения движения стержня и проинтегрировать их. При этом начальными условиями для второго этапа будут найденный угол ϕ1, при котором происходит отрыв точки А стержня от вертикальной стены, и угловая скорость ˙ϕ1, которой обладал стержень в этот момент времени.
Следует отметить, что единообразие в составлении уравнений Лагранжа, с одной стороны являясь его достоинством, с другой стороны не позволяет, в некоторых случаях, произвести анализ действия сил, и, как следствие, поведения механической системы, что для инженера оказывается серьезным недостатком. Избежать этого позволяет комбинирование уравнений Лагранжа с общими теоремами механики (как это было сделано выше).