
- •Учебное пособие
- •Предисловие
- •Список основных обозначений
- •Выдающиеся ученые, создатели основ механики
- •ВВЕДЕНИЕ
- •1. МАТЕМАТИЧЕСКИЕ СВЕДЕНИЯ
- •1.1. Матрицы и векторы
- •1.2. Конкатенация матриц и линейных функций
- •1.3. Перепроецирование векторов и тензоров второго ранга
- •1.4. Решение системы алгебраических линейных уравнений MATLAB
- •2. ДИНАМИКА МАТЕРИАЛЬНОЙ ТОЧКИ
- •2.1. Инерциальные системы отсчета
- •2.2. Поступательные (невращающиеся) система отсчета
- •2.4. Пример
- •2.5. Второй закон динамики в кинетостатической форме Даламбера
- •2.6. Закон динамики относительного движения
- •2.7. Уравнения относительного движения в форме Даламбера
- •2.8. Основной закон динамики относительного движения
- •2.9. Случай поступательного движения системы отсчета
- •2.10. Случай поступательной системы отсчета с началом на поверхности Земли
- •2.11. Прямолинейное движение материальной точки во вращающейся системе отсчета
- •3. ДИНАМИКА ТВЕРДОГО ТЕЛА
- •3.1. Геометрия масс твердого тела
- •3.2. Статические моменты массы и координаты центра масс тела
- •3.3. Осевые моменты инерции и радиусы инерции тела
- •3.4. Осевые моменты инерции площади поперечного сечения тела
- •3.5. Центробежные моменты инерции твердого тела
- •3.6. Перепроецирование вектора с применением матрицы поворота
- •3.7. Тензор инерции в точке твердого тела
- •3.8. Угловая скорость сферического движения тела
- •3.9. Кинетический момент сферического движения тела
- •3.10. Динамические уравнения Эйлера
- •3.11. Элементарная теория гироскопа
- •3.12. Динамика произвольного движения твердого тела
- •4.1. Сведения из статики
- •4.2. Момент силы и пары сил, расположенных в плоскости
- •4.3 Бивектор плоской системы сил
- •4.4. Уравнения равновесия плоской системы сил
- •4.5 Бивектор и M-функции пространственной системы сил
- •4.6. Матричное уравнение равновесия пространственной системы сил
- •4.7. Бивекторы реакций опор
- •4.8. Пример конкатенации матриц в матричном уравнении равновесия
- •4.9. Кинетостатическое матричное уравнение произвольного движения тела
- •5. ДИНАМИКА ГОЛОНОМНЫХ МЕХАНИЧЕСКИХ СИСТЕМ
- •5.1. Фазовый вектор системы (вектор состояния системы)
- •5.2. Кинетическая энергия твердого тела и голономной стационарной одностепенной системы
- •5.3. Мощность силы
- •5.4. Работа и потенциальная энергия
- •5.5. Обобщенная сила одностепенной голономной стационарной системы
- •5.6. Обобщенная сила нестационарной одностепенной системы
- •5.7. Работа обобщенной силы одностепенной стационарной системы
- •5.8. Идеальные связи и реакции связей
- •5.9. Уравнение Лагранжа для одностепенной стационарной голономной системы
- •5.10. Уравнение Лагранжа для голономной нестационарной одностепенной системы
- •6. ОБЩИЕ ТЕОРЕМЫ ДИНАМИКИ МЕХАНИЧЕСКОЙ СИСТЕМЫ
- •6.1. Уравнения Лагранжа для голономной стационарной системы
- •6.2.Уравнения равновесия стационарной системы в обобщенных координатах
- •6.3. Матричная форма уравнений Лагранжа
- •6.4. Матричное уравнение Гамильтона
- •6.5. Теорема о движении центра масс механической системы
- •6.6. Теорема об изменении количества движения механической системы
- •6.7. Теорема импульсов
- •6.8. Теорема об изменении кинетической энергии
- •6.9. Теорема об изменении кинетического момента механической системы
- •Приложение
- •Равновесие плоской системы сил
- •Равновесие пространственной системы сил
- •Динамика относительного движения материальной точки
- •Литература

В случае нестационарной голономной системы реакция связи называется идеальной, если равна нулю её мощность при зафиксированном значении параметра t в связи.
5.9.Уравнение Лагранжа для одностепенной стационарной голономной системы
Рассмотрим голономную стационарную механическую систему, имеющую одну степень свободы, совершающую движение в инерциальной или поступательной системе отсчета. Движение механических устройств с твердыми звеньями в поступательной системе отсчета, движущейся с ускорением a0 эквивалентно их движению в ИСО, но при условии замены ускорения си-
лы тяжести g на измененное ускорение g0 = g −a0 , что эквивалентно присвоению силе тяжести каждого звена Gi значения Gi0 = Gi + Si при Si = −mi a0
Допустим, что мы выбрали обобщенную координату q и фазовый вектор [q&, q] . Пусть по формулам (5.1)-(5.6) определена кинетическая энергия меха-
низма в ИСО (либо в ПСО) как сумма кинетических энергий всех звеньев. Выразим все кинематические переменные через параметры q&, q , затем выне-
|
|
&2 |
. В результате получаем выражение вида |
||
сем за скобки общий множитель q |
|||||
с некоторым множителем a(q), называемым коэффициентом инерции: |
|||||
& |
2 |
где |
&2 |
. |
(5.35) |
2T = a(q)q |
|
a = 2T / q |
|
Цель перечисленных выше действий — получение выражения для коэффициента инерции a(q), называемого также приведенной массой mпр системы, если в качестве q принята декартова координата и приведенным моментом инерции Jпр механизма, если q есть угол. Приведенный коэффициент инерции может получиться постоянным или зависящим от обобщенной координаты. Обратный коэффициент инерции вычисляется по формуле b(q) = (a(q))−1 = q&2 / 2T
Далее допустим, что мы показали на рисунке “положительную фазу” го-
лономной системы [q > 0, q > 0], независимо от заданной начальной фазы |
||||||||
|
& |
|
|
|
|
|
|
|
& |
|
|
|
|
|
|
|
|
|
|
|
|
M |
|
j пар |
||
[q0 , q0 ], а также показали векторы приложенных сил Fi и моменты |
|
|
|
|||||
сил, |
скорости vi точек приложения сил и угловые скорости ω |
j пар сил твер- |
дых тел-звеньев (идеальные реакции можно не показывать). Затем пусть в мощности приложенных сил P = ΣFi vi + ΣM j ωj , выразили все кинемати-
ческие величины через фазовые переменные q, q и вынесли за скобки общий |
|||
множитель – обобщенную скорость |
& |
|
|
q , а выражение в скобках обозначили |
|||
|
|
& |
|
символом Q . В результате получили следующую математическую структуру |
|||
мощности: |
|
Q = P / q , |
(5.36) |
P =Q(q, q,t)q, |
|||
& |
& |
& |
|
где обобщенная сила Q в общем случае есть функция фазовых координат и времени. Таким образом, обобщенная сила численно равна “удельной мощ-
83
ности” системы сил, т.е. - мощности системы приложенных сил, отнесенной к обобщенной скорости. Обобщенная сила может получиться постоянной величиной, либо — функцией координаты q или фазового вектора [q&, q] , либо
— функцией “расширенного” фазового вектора [q&, q,t]. Последний, наиболее общий, вариант реализуется в случаях, когда имеются силы, явно зависящие от времени.
Теорема Лагранжа 1. Движение одностепенной голономной стационарной механической системы в ИСО (или ПСО с дополнительными силами) подчинено дифференциальному уравнению второго порядка, которое составляется по формуле:
d |
∂T |
|
− |
∂T |
= Q при |
& |
(5.37) |
|
|
|
& |
|
|
|
|||
|
∂q |
Q = P / q . |
|
|||||
dt |
∂q |
|
|
|
|
Таким образом, полная производная по времени от так называемого обобщенного импульса ∂T / ∂q&, за вычетом частной производной вида ∂T / ∂q , равна обобщенной силе Q .
Замечание 1. Действующие силы можно подразделять на потенциальные и непотенциальные и применять уравнение Лагранжа в форме
|
|
|
|
|
|
|
|
|
d ∂T |
∂T |
|
∂Π |
−Q |
нп |
= 0 |
(5.38) |
|||
|
|
|
|
|
|
|
|
|
|
|
& − |
∂q |
+ |
∂q |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
П(q) |
|
|
dt |
∂q |
|
|
|
|
|
|
|||||
Здесь |
— потенциальная |
энергия |
системы потенциальных |
сил, |
|||||||||||||||
нп |
|
|
нп & |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Q |
= P |
|
/ q — обобщенная сила системы непотенциальных сил. |
|
|
||||||||||||||
Иногда |
|
уравнения Лагранжа |
записывают |
в форме Q − LT = 0 , |
где |
||||||||||||||
L = |
|
d ∂ |
∂ |
— называется оператором Лагранжа и представляет собой |
|||||||||||||||
|
|
|
|
|
− |
|
|||||||||||||
|
|
& |
∂q |
||||||||||||||||
|
|
dt |
|
∂q |
|
|
|
|
|
|
|
|
|
|
|
|
совокупность математических действий над функцией T. Функцию S = -LT можно назвать обобщенной силой инерции системы. В таком случае можно сказать, что при движении механической системы, обобщенная сила системы приложенных сил уравновешивается обобщенной силой системы инерционных сил , Q + S = 0.
Замечание 2. Уравнению Лагранжа для стационарных голономных одно-
степенных систем можно придать вид |
= |
& |
(5.39) |
|
& |
= Q при Q |
|
P / q |
|
T / q& |
|
|
т.е. полная производная по времени от кинетической энергии одностепенной стационарной системы, отнесенная к обобщенной скорости равна обобщенной силе.
В результате выполнения математиеских действий в уравнении (5.37) или (5.39) получаем следующий более конкретный вид динамического уравнения.
Теорема 2. Математической моделью движения одностепенной голономной стационарной системы служит дифференциальное уравнение второго порядка
84

&& |
|
|
&2 |
′ |
/(2a) = 0 |
при |
a′ = ∂a(q) |
≡ da(q) |
(5.40) |
|
q |
−Q / a + q |
a |
∂q |
dq |
|
|||||
|
|
|
|
|
|
|
|
|||
или |
|
|
|
|
|
при b = a−1 , b′ = db(q) |
(5.41) |
|||
q&&−bQ − q&2b′/(2b) = 0 |
||||||||||
|
||||||||||
|
|
|
|
|
|
|
|
dq |
|
|
|
& |
|
|
& |
& |
|
|
|
|
|
Здесь Q(q, q,t) = P(q, q,t) / q - обобщенная сила механической системы, |
||||||||||
&2 |
- приведенный коэффициент инерции механической системы, |
|||||||||
a(q) = 2T / q |
||||||||||
|
& |
2 |
/ 2T - обратный коэффициент инерции механической сис- |
|||||||
b(q) =1/ a(q) = q |
|
темы. Отметим, что наиболее предпочтительным при практическом применении является уравнение (5.41), но целесообразно знать и классические формы уравнения Лагранжа (5.37), (5.38).
Кроме уравнения (5.40) или (5.41), в математическую модель включают начальные условия, начальное состояние механической системы, задаваемое отдельно в виде двух конкретных начальных значений фазовых переменных:
q(0) = q0 , |
q(0) = q0 |
, или [q(0), q(0)] =[q0 , q0 ] |
(5.42) |
||
& |
& |
|
& |
& |
|
Вместо начальных условий иногда присоединяют граничные условия вида q(t1) = q1, q(t2 ) = q2 , своеобразным граничным условием может являться
условие периодичности движения.
Теорема 3. Движение одностепенной голономной стационарной механи-
ческой системы |
в случае постоянного обратного коэффициента инерции |
|||||
b=const подчинено дифференциальному уравнению второго порядка вида |
||||||
&& |
при b = a |
−1 |
&2 |
/ 2T = const , |
& |
(5.43) |
q −bQ = 0 |
|
= q |
Q = P / q |
|
Уравнение (5.43) имеет простой смысл: обобщенное ускорение системы во все моменты времени пропорционально обобщенной силе с коэффициентом пропорциональности равным обратному коэффициенту инерции системы.
В случае потенциальных сил обобщенное ускорение системы пропорционально производной от потенциальной энергии по обобщенной координате:
&& |
+b |
dΠ(q) |
= 0 |
(5.44) |
dq |
|
|||
q |
|
Теорема 4. О равновесии голономной одностепенной системы. Необходимым условием равновесия одностепенной механической системы является равенство нулю обобщенной силы механической системы:
Q(q&, q) q&=0 = 0 (5.45)
Итак, должна равняться нулю обобщенная сила Q механизма, в которой сле-
дует приравнять нулю обобщенную скорость, например, не учитывать диссипативные силы, обращающиеся в ноль при q& = 0. В случае консервативной механической системы, когда все приложенные активные силы потенциальны, а связи - идеальны, необходимым условием равновесия является выполнение функционального уравнения, получаемого из условия
85

dΠ(q) |
= 0 |
(5.46) |
dq |
|
|
|
|
т.е. - равенство нулю производной от потенциальной энергии по обобщенной координате. Из функционального (линейного или нелинейного) уравнения (5.46) можно найти положение статического равновесия системы в виде q = qСТ = const . В задачах статики закрепленного тела можно поочередно ос-
вободить тело от той или иной связи, придавая ему степень свободы, вводя одну обобщенную координату и вычислять одну реакцию отброшенной связи через уравнение вида (5.46)
Пример 1.
y1 30 |
|
|
|
c |
y |
|
F |
j1 |
|
|
|
|
|
N |
G |
O |
l |
x.0 |
. |
|
i |
x0 |
x=vx |
|
|
x x
Рис. 17 Свободные колебания одностепенной системы.
Тележка на четырех колесах удерживается пружиной жесткости c на наклонной плоскости с наклоном 30° (рис. 17). Масса корпуса тележки m1, масса каждого колеса m2, радиус инерции колеса относительно его оси равен ρ. В начальное мгновение тележка занимала такое положение, при котором пружина растянута на величину l , и имела начальную скоростью v0, направленную по склону. Определить дальнейшее движение тележки. Трением пренебречь.
Решение. Механизм состоит из пяти подвижных звеньев: корпуса, четырех колес и невесомой пружины. В качестве обобщенной координаты возьмем растяжение пружины x, от нерастянутого положения конца пружины.
Покажем начальное состояние механизма [x0 |
= l, x0 |
= v0 ] и положительное |
|
произвольное состояние [ x > 0 |
, x > |
& |
Приложенные силы: |
0]. |
|||
& |
|
|
|
G = mg, F = −cxi - сила тяжести механизма и упругая сила. Реакция N опо-
ры на действия колес, показаная пунктиром, включает нормальную реакцию и силу сцепления колес с дорогой. В дальнейшем она не будет принимаеться во внимание, поскольку не имеет мощности (точка приложения силы N не имеет скорости).
Необходимо выразить все кинематические параметры через вектор состояния [x&, x] . В данном случае выразим угловую скорость и угол поворота
каждого колеса: ω = x&/ r , φ = x/r. Двойная кинетическая энергия механизма определяется по формулам поступательного движения корпуса и плоского движения четырех колес:
86

2T = 2Tкорп +8Tколес = m1x&2 +4(m2 x&2 + m2 ρ2 x&2 / r2 ) = mпрx&2 = x&2 / b,
при mпр = m1 + 4m2 (1 + ρ2/r2), b = 1/ mпр .
Таким образом найдена приведенная масса механизма mпр и обратная приведенная масса b, они получились постоянными. Перейдем к определению мощности упругой силы F = −cxi и общей силы тяжести G = −mgj при m = m1 +4m2. Получаем
P = F v +G v = Fx vx + Gx vx = (−cx + mg 3 / 2)x&
где m = m1 +4m2. Отсюда находим обобщенную силу устройства
Q = P / x& = −cx + mg 32
Поскольку в данной задаче приведенная масса mпр механизма постоянна, то применяем теорему Лагранжа 3, в виде уравнения (5.43) при обозначении x=q. Получаем:
|
|
x − bQ = 0; |
x + bcx − bmg 3 / 2 = 0. |
||
|
|
&& |
|
&& |
|
Окончательно |
|
|
|
||
&& |
2 |
x −h = 0 при k = |
bc = |
c mпр , |
h = bmg 3 / 2 = mg 3 /(2mпр) |
x + k |
|
||||
Начальное состояние системы: |
[x(0) = v0 , |
x(0) = l] |
|||
|
|
|
& |
|
|
Постоянное частное решение дифференциального уравнения xч = h / k2 определяет положение равновесия системы. Общее решение имеет вид
x = xч +C1 cos kt + C2 sin kt.
Подставив начальные значения в общее решение и в выражение производной от него, находим
b = xч + C1 , v0 = C2k.
Отсюда C1 = b - x2 , C2 = v0 / k Окончательный вид решения:
x = (b - h / k2 ) cos kt + (v0 / k ) sin kt + h / k2
Таким образом, система совершает гармонические колебания с периодом τ = 2π/k вокруг положения статического равновесия xч = h/k2.
Замечание Мы отсчитывали координату x от нерастянутого положения пружины. Если отсчитывать координату от положения статического равновесия, то динамическое уравнение получается без постоянного слагаемого xч .
Действительно, в таком случае текущая деформация пружины состоит из деформации f в положении равновесия и координаты x. Упругая сила
имеет вид F = −c( f + x)i .
& |
|
|
& |
3 |
mg . |
|
|
|
|
||||
|
|
|
2 |
|||
Обобщенная сила: Q = P / x = (F v +G v ) / x = −cf −cx + |
Мы взяли начало координат x = 0 в положении равновесия, т.е. в положении, в котором обобщенная сила равна нулю. Приравнивая нулю обобщенную силу Q и подставляя x=0, находим величину статической деформации пружины:
87