Динамика жесткого тела

Исследование воздействия сил на недеформируемые тела

Движение каждого из компонентов паровой машины Бултона и Уатта (1784 г.) можно описать системой уравнений кинематики и кинетики.

В физической науке динамика динамика твердого тела изучает движение систем взаимосвязанных тел под действием внешних сил . Предположение о том, что тела являются жесткими (т.е. не деформируются под действием приложенных сил), упрощает анализ, сводя параметры, описывающие конфигурацию системы, к перемещению и вращению систем отсчета , прикрепленных к каждому телу. ^{[1] [2]} Сюда не входят тела, обладающие текучестью , высокой упругостью и пластичностью .

Динамика системы твердого тела описывается законами кинематики и применением второго закона Ньютона ( кинетики ) или его производной формы — лагранжевой механики . Решение этих уравнений движения дает описание положения, движения и ускорения отдельных компонентов системы и всей системы в целом как функции времени . Формулировка и решение динамики твердого тела является важным инструментом компьютерного моделирования механических систем .

Плоская динамика твердого тела

Если система частиц движется параллельно неподвижной плоскости, говорят, что система ограничена в плоском движении. В этом случае законы Ньютона (кинетика) для жесткой системы из N частиц Pi, _i = 1 ,..., N упрощаются, поскольку нет движения в направлении k . Определите результирующую силу и крутящий момент в контрольной точке R , чтобы получить

$${\displaystyle \mathbf {F} =\sum _{i=1}^{N}m_{i}\mathbf {A} _{i},\quad \mathbf {T} =\sum _{i=1}^{N}(\mathbf {r} _{i}-\mathbf {R} )\times m_{i}\mathbf {A} _{i},}$$

где r _i обозначает плоскую траекторию каждой частицы.

Кинематика твердого тела дает формулу ускорения частицы P _i через положение R и ускорение A опорной частицы , а также вектор угловой скорости ω и вектор углового ускорения α твердой системы частиц как ,

$${\displaystyle \mathbf {A} _{i}={\boldsymbol {\alpha }}\times (\mathbf {r} _{i}-\mathbf {R} )+{\boldsymbol {\omega }}\times ({\boldsymbol {\omega }}\times (\mathbf {r} _{i}-\mathbf {R} ))+\mathbf {A} .}$$

Для систем, ограниченных плоскостным движением, векторы угловой скорости и углового ускорения направлены вдоль k перпендикулярно плоскости движения, что упрощает это уравнение ускорения. В этом случае векторы ускорения можно упростить, введя единичные векторы e _i из опорной точки R в точку r _i и единичные векторы , так ${\textstyle \mathbf {t} _{i}=\mathbf {k} \times \mathbf {e} _{i}}$

$${\displaystyle \mathbf {A} _{i}=\alpha (\Delta r_{i}\mathbf {t} _{i})-\omega ^{2}(\Delta r_{i}\mathbf {e} _{i})+\mathbf {A} .}$$

Это дает результирующую силу, действующую на систему, как

$${\displaystyle \mathbf {F} =\alpha \sum _{i=1}^{N}m_{i}\left(\Delta r_{i}\mathbf {t} _{i}\right)-\omega ^{2}\sum _{i=1}^{N}m_{i}\left(\Delta r_{i}\mathbf {e} _{i}\right)+\left(\sum _{i=1}^{N}m_{i}\right)\mathbf {A} ,}$$

$${\displaystyle {\begin{aligned}\mathbf {T} ={}&\sum _{i=1}^{N}(m_{i}\Delta r_{i}\mathbf {e} _{i})\times \left(\alpha (\Delta r_{i}\mathbf {t} _{i})-\omega ^{2}(\Delta r_{i}\mathbf {e} _{i})+\mathbf {A} \right)\\{}={}&\left(\sum _{i=1}^{N}m_{i}\Delta r_{i}^{2}\right)\alpha \mathbf {k} +\left(\sum _{i=1}^{N}m_{i}\Delta r_{i}\mathbf {e} _{i}\right)\times \mathbf {A} ,\end{aligned}}}$$

где и — единичный вектор, перпендикулярный плоскости для всех частиц _Pi . ${\textstyle \mathbf {e} _{i}\times \mathbf {e} _{i}=0}$ ${\textstyle \mathbf {e} _{i}\times \mathbf {t} _{i}=\mathbf {k} }$

Используйте центр масс C в качестве опорной точки, чтобы эти уравнения законов Ньютона упростились до следующего вида:

$${\displaystyle \mathbf {F} =M\mathbf {A} ,\quad \mathbf {T} =I_{\textbf {C}}\alpha \mathbf {k} ,}$$

где M — полная масса, а I _C — момент инерции относительно оси, перпендикулярной движению жесткой системы и проходящей через центр масс.

Твердое тело в трех измерениях

Описания ориентации или отношения

Было разработано несколько методов описания ориентации твердого тела в трех измерениях. Они обобщены в следующих разделах.

углы Эйлера

Первая попытка представления ориентации приписывается Леонарду Эйлеру . Он представил себе три системы отсчета, которые могли вращаться одна вокруг другой, и понял, что, начав с фиксированной системы отсчета и выполнив три вращения, он может получить любую другую систему отсчета в пространстве (используя два вращения для фиксации вертикальной оси и еще одно для фиксации). зафиксируйте две другие оси). Значения этих трех поворотов называются углами Эйлера . Обычно используется для обозначения прецессии, нутации и внутреннего вращения. ${\displaystyle \psi }$ ${\displaystyle \theta }$ ${\displaystyle \phi }$

• 
  Схема углов Эйлера
• 
  Собственное вращение шара вокруг неподвижной оси.
• 
  Движение волчка по углам Эйлера.

Углы Тейта – Брайана

Углы Тейта – Брайана, еще один способ описания ориентации.

Это три угла, также известные как рыскание, тангаж и крен, углы навигации и углы кардана. Математически они представляют собой набор из шести возможностей внутри двенадцати возможных наборов углов Эйлера, причем этот порядок лучше всего использовать для описания ориентации транспортного средства, например самолета. В аэрокосмической технике их обычно называют углами Эйлера.

Вектор ориентации

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

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

Подобный метод, называемый представлением угла оси , описывает вращение или ориентацию с использованием единичного вектора, совмещенного с осью вращения, и отдельного значения для обозначения угла (см. рисунок).

Матрица ориентации

С введением матриц теоремы Эйлера были переписаны. Вращения описывались ортогональными матрицами , называемыми матрицами вращения или матрицами направленного косинуса. При использовании для представления ориентации матрицу вращения обычно называют матрицей ориентации или матрицей отношения.

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

Конфигурационное пространство несимметричного объекта в n -мерном пространстве — это SO ( n ) × Rn . Ориентацию можно визуализировать, прикрепив к объекту основу из касательных векторов . Направление, в котором указывает каждый вектор, определяет его ориентацию.

Кватернион ориентации

Другой способ описания вращения — использование кватернионов вращения , также называемых версорами. Они эквивалентны матрицам вращения и векторам вращения. Что касается векторов вращения, их легче преобразовывать в матрицы и обратно. При использовании для представления ориентации кватернионы вращения обычно называются кватернионами ориентации или кватернионами ориентации.

Второй закон Ньютона в трех измерениях

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

Ньютон сформулировал свой второй закон для частицы следующим образом: «Изменение движения объекта пропорционально приложенной силе и происходит в направлении прямой линии, по которой действует сила». ^[3] Поскольку Ньютон обычно называл «движением» частицы массу, умноженную на скорость, фраза «изменение движения» относится к произведению массы, умноженной на ускорение частицы, и поэтому этот закон обычно записывается как

$${\displaystyle \mathbf {F} =m\mathbf {a} ,}$$

под Fm—

Жесткая система частиц

Если система из N частиц Pi _, i=1,..., N собрана в твердое тело, то второй закон Ньютона можно применить к каждой из частиц в теле. Если F _i — внешняя сила, приложенная к частице Pi _с массой m _i , то

$${\displaystyle \mathbf {F} _{i}+\sum _{j=1}^{N}\mathbf {F} _{ij}=m_{i}\mathbf {a} _{i},\quad i=1,\ldots ,N,}$$

F _ijji

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

Важным упрощением этих уравнений сил является введение результирующей силы и крутящего момента, действующих на жесткую систему. Эта результирующая сила и крутящий момент получаются путем выбора одной из частиц в системе в качестве опорной точки R , где каждая из внешних сил применяется с добавлением соответствующего крутящего момента. Результирующая сила F и крутящий момент T определяются по формулам:

$${\displaystyle \mathbf {F} =\sum _{i=1}^{N}\mathbf {F} _{i},\quad \mathbf {T} =\sum _{i=1}^{N}(\mathbf {R} _{i}-\mathbf {R} )\times \mathbf {F} _{i},}$$

R _iPi

Второй закон Ньютона для частицы сочетается с этими формулами для получения результирующей силы и крутящего момента:

$${\displaystyle \mathbf {F} =\sum _{i=1}^{N}m_{i}\mathbf {a} _{i},\quad \mathbf {T} =\sum _{i=1}^{N}(\mathbf {R} _{i}-\mathbf {R} )\times (m_{i}\mathbf {a} _{i}),}$$

F _ijiRa опорной частицы ,

$${\displaystyle \mathbf {a} _{i}=\alpha \times (\mathbf {R} _{i}-\mathbf {R} )+\omega \times (\omega \times (\mathbf {R} _{i}-\mathbf {R} ))+\mathbf {a} .}$$

Массовые свойства

Массовые свойства твердого тела представлены его центром масс и матрицей инерции . Выбрать опорную точку R так, чтобы она удовлетворяла условию

$${\displaystyle \sum _{i=1}^{N}m_{i}(\mathbf {R} _{i}-\mathbf {R} )=0,}$$

тогда он известен как центр масс системы.

Матрица инерции [I _R ] системы относительно опорной точки R определяется выражением

$${\displaystyle [I_{R}]=\sum _{i=1}^{N}m_{i}\left(\mathbf {I} \left(\mathbf {S} _{i}^{\textsf {T}}\mathbf {S} _{i}\right)-\mathbf {S} _{i}\mathbf {S} _{i}^{\textsf {T}}\right),}$$

где вектор-столбец R _i − R ; является его транспонированием и представляет собой единичную матрицу размером 3 на 3. ${\displaystyle \mathbf {S} _{i}}$ ${\displaystyle \mathbf {S} _{i}^{\textsf {T}}}$ ${\displaystyle \mathbf {I} }$

${\displaystyle \mathbf {S} _{i}^{\textsf {T}}\mathbf {S} _{i}}$является скалярным произведением самого себя, а является тензорным произведением самого себя. ${\displaystyle \mathbf {S} _{i}}$ ${\displaystyle \mathbf {S} _{i}\mathbf {S} _{i}^{\textsf {T}}}$ ${\displaystyle \mathbf {S} _{i}}$

Уравнения силы и момента

Используя центр масс и матрицу инерции, уравнения силы и крутящего момента для одного твердого тела принимают вид

$${\displaystyle \mathbf {F} =m\mathbf {a} ,\quad \mathbf {T} =[I_{R}]\alpha +\omega \times [I_{R}]\omega ,}$$

Динамика взаимосвязанной системы твердых тел B _i , j = 1, ..., M формулируется путем изоляции каждого твердого тела и введения сил взаимодействия. Результирующая внешних сил и сил взаимодействия, действующих на каждое тело, дает уравнения силы и момента

$${\displaystyle \mathbf {F} _{j}=m_{j}\mathbf {a} _{j},\quad \mathbf {T} _{j}=[I_{R}]_{j}\alpha _{j}+\omega _{j}\times [I_{R}]_{j}\omega _{j},\quad j=1,\ldots ,M.}$$

Формулировка Ньютона дает 6 M уравнений, которые определяют динамику системы M твердых тел. ^[4]

Вращение в трех измерениях

Вращающийся объект, независимо от того, находится ли он под действием крутящих моментов или нет, может проявлять поведение прецессии и нутации . Фундаментальным уравнением, описывающим поведение вращающегося твердого тела, является уравнение движения Эйлера :

$${\displaystyle {\boldsymbol {\tau }}={\frac {D\mathbf {L} }{Dt}}={\frac {d\mathbf {L} }{dt}}+{\boldsymbol {\omega }}\times \mathbf {L} ={\frac {d(I{\boldsymbol {\omega }})}{dt}}+{\boldsymbol {\omega }}\times {I{\boldsymbol {\omega }}}=I{\boldsymbol {\alpha }}+{\boldsymbol {\omega }}\times {I{\boldsymbol {\omega }}}}$$

псевдовекторы τLкрутящие моментыугловой моментIмомент инерцииωα

Решение этого уравнения при отсутствии приложенного крутящего момента обсуждается в статьях Уравнение движения Эйлера и Эллипсоид Пуансо .

Из уравнения Эйлера следует, что крутящий момент τ , приложенный перпендикулярно оси вращения и, следовательно, перпендикулярно L , приводит к вращению вокруг оси, перпендикулярной как τ , так и L. Это движение называется прецессией . Угловая скорость прецессии Ω _P определяется векторным произведением : ^{[ необходима ссылка ]}

$${\displaystyle {\boldsymbol {\tau }}={\boldsymbol {\Omega }}_{\mathrm {P} }\times \mathbf {L} .}$$

Прецессия гироскопа

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

При постоянном крутящем моменте величины τ скорость прецессии Ω _P обратно пропорциональна L , величине его углового момента:

$${\displaystyle \tau ={\mathit {\Omega }}_{\mathrm {P} }L\sin \theta ,}$$

θΩ _PL.

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

Виртуальная работа сил, действующих на твердое тело.

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

Виртуальную работу сил, действующих в различных точках одного твердого тела, можно рассчитать, используя скорости их точек приложения и результирующие силу и крутящий момент . Чтобы убедиться в этом, пусть силы F ₁ , F ₂ ... F _n действуют на точки R ₁ , R ₂ ... R _n в твердом теле.

Траектории R _i , i = 1, ..., n определяются движением твердого тела. Скорость точек R _i вдоль их траекторий равна

$${\displaystyle \mathbf {V} _{i}={\boldsymbol {\omega }}\times (\mathbf {R} _{i}-\mathbf {R} )+\mathbf {V} ,}$$

ω

Виртуальная работа

Работа вычисляется как скалярное произведение каждой силы со смещением точки ее контакта.

$${\displaystyle \delta W=\sum _{i=1}^{n}\mathbf {F} _{i}\cdot \delta \mathbf {r} _{i}.}$$

обобщенных координат q _jj = 1, ..., mδ r _i

$${\displaystyle \delta \mathbf {r} _{i}=\sum _{j=1}^{m}{\frac {\partial \mathbf {r} _{i}}{\partial q_{j}}}\delta q_{j}=\sum _{j=1}^{m}{\frac {\partial \mathbf {V} _{i}}{\partial {\dot {q}}_{j}}}\delta q_{j}.}$$

$${\displaystyle \delta W=\mathbf {F} _{1}\cdot \left(\sum _{j=1}^{m}{\frac {\partial \mathbf {V} _{1}}{\partial {\dot {q}}_{j}}}\delta q_{j}\right)+\dots +\mathbf {F} _{n}\cdot \left(\sum _{j=1}^{m}{\frac {\partial \mathbf {V} _{n}}{\partial {\dot {q}}_{j}}}\delta q_{j}\right)}$$

или собирая коэффициенты δq _j

$${\displaystyle \delta W=\left(\sum _{i=1}^{n}\mathbf {F} _{i}\cdot {\frac {\partial \mathbf {V} _{i}}{\partial {\dot {q}}_{1}}}\right)\delta q_{1}+\dots +\left(\sum _{1=1}^{n}\mathbf {F} _{i}\cdot {\frac {\partial \mathbf {V} _{i}}{\partial {\dot {q}}_{m}}}\right)\delta q_{m}.}$$

Обобщенные силы

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

$${\displaystyle \delta W=\left(\sum _{i=1}^{n}\mathbf {F} _{i}\cdot {\frac {\partial \mathbf {V} _{i}}{\partial {\dot {q}}}}\right)\delta q=\left(\sum _{i=1}^{n}\mathbf {F} _{i}\cdot {\frac {\partial ({\boldsymbol {\omega }}\times (\mathbf {R} _{i}-\mathbf {R} )+\mathbf {V} )}{\partial {\dot {q}}}}\right)\delta q.}$$

Введем результирующую силу F и крутящий момент T , чтобы это уравнение приняло вид

$${\displaystyle \delta W=\left(\mathbf {F} \cdot {\frac {\partial \mathbf {V} }{\partial {\dot {q}}}}+\mathbf {T} \cdot {\frac {\partial {\boldsymbol {\omega }}}{\partial {\dot {q}}}}\right)\delta q.}$$

Величина Q, определяемая формулой

$${\displaystyle Q=\mathbf {F} \cdot {\frac {\partial \mathbf {V} }{\partial {\dot {q}}}}+\mathbf {T} \cdot {\frac {\partial {\boldsymbol {\omega }}}{\partial {\dot {q}}}},}$$

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

$${\displaystyle \delta W=\sum _{j=1}^{m}Q_{j}\delta q_{j},}$$

$${\displaystyle Q_{j}=\mathbf {F} \cdot {\frac {\partial \mathbf {V} }{\partial {\dot {q}}_{j}}}+\mathbf {T} \cdot {\frac {\partial {\boldsymbol {\omega }}}{\partial {\dot {q}}_{j}}},\quad j=1,\ldots ,m.}$$

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

$${\displaystyle Q_{j}=-{\frac {\partial V}{\partial q_{j}}},\quad j=1,\ldots ,m.}$$

Форма принципа виртуальной работы Даламбера.

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

Статическое равновесие

Статическое равновесие твердых тел механической системы определяется условием, что виртуальная работа приложенных сил равна нулю при любом виртуальном перемещении системы. Это известно как принцип виртуальной работы. ^[5] Это эквивалентно требованию, чтобы обобщенные силы для любого виртуального смещения были равны нулю, то есть Q _i =0.

Пусть механическая система построена из n твердых тел B _i , i = 1, ..., n , и пусть результирующая приложенных сил к каждому телу представляет собой пары сила- момент F _i и Ti _, i = 1, ..., н . Обратите внимание, что эти приложенные силы не включают силы реакции в местах соединения тел. Наконец, предположим, что скорость Vi и угловые скорости ω _i , _i = 1, ..., n для каждого твердого тела определяются одной обобщенной координатой q. Говорят, что такая система твердых тел имеет одну степень свободы .

Виртуальная работа сил и моментов Fi и Ti , приложенная к этой системе _{с одной степенью свободы ,} определяется выражением

$${\displaystyle \delta W=\sum _{i=1}^{n}\left(\mathbf {F} _{i}\cdot {\frac {\partial \mathbf {V} _{i}}{\partial {\dot {q}}}}+\mathbf {T} _{i}\cdot {\frac {\partial {\boldsymbol {\omega }}_{i}}{\partial {\dot {q}}}}\right)\delta q=Q\delta q,}$$

$${\displaystyle Q=\sum _{i=1}^{n}\left(\mathbf {F} _{i}\cdot {\frac {\partial \mathbf {V} _{i}}{\partial {\dot {q}}}}+\mathbf {T} _{i}\cdot {\frac {\partial {\boldsymbol {\omega }}_{i}}{\partial {\dot {q}}}}\right),}$$

Если механическая система определяется m обобщенными координатами, q _j , j = 1, ..., m , то система имеет m степеней свободы, а виртуальная работа определяется выражением

$${\displaystyle \delta W=\sum _{j=1}^{m}Q_{j}\delta q_{j},}$$

$${\displaystyle Q_{j}=\sum _{i=1}^{n}\left(\mathbf {F} _{i}\cdot {\frac {\partial \mathbf {V} _{i}}{\partial {\dot {q}}_{j}}}+\mathbf {T} _{i}\cdot {\frac {\partial {\boldsymbol {\omega }}_{i}}{\partial {\dot {q}}_{j}}}\right),\quad j=1,\ldots ,m.}$$

q _j

$${\displaystyle Q_{j}=0,\quad j=1,\ldots ,m.}$$

Эти m уравнений определяют статическое равновесие системы твердых тел.

Обобщенные силы инерции

Рассмотрим одно твердое тело, которое движется под действием равнодействующей силы F и крутящего момента T , с одной степенью свободы, определяемой обобщенной координатой q . Предположим, что точкой отсчета результирующей силы и крутящего момента является центр масс тела, тогда обобщенная сила инерции Q* , связанная с обобщенной координатой q, определяется выражением

$${\displaystyle Q^{*}=-(M\mathbf {A} )\cdot {\frac {\partial \mathbf {V} }{\partial {\dot {q}}}}-\left([I_{R}]{\boldsymbol {\alpha }}+{\boldsymbol {\omega }}\times [I_{R}]{\boldsymbol {\omega }}\right)\cdot {\frac {\partial {\boldsymbol {\omega }}}{\partial {\dot {q}}}}.}$$

Эту силу инерции можно вычислить по кинетической энергии твердого тела:

$${\displaystyle T={\tfrac {1}{2}}M\mathbf {V} \cdot \mathbf {V} +{\tfrac {1}{2}}{\boldsymbol {\omega }}\cdot [I_{R}]{\boldsymbol {\omega }},}$$

$${\displaystyle Q^{*}=-\left({\frac {d}{dt}}{\frac {\partial T}{\partial {\dot {q}}}}-{\frac {\partial T}{\partial q}}\right).}$$

Система n твердых тел с m обобщенными координатами имеет кинетическую энергию

$${\displaystyle T=\sum _{i=1}^{n}\left({\tfrac {1}{2}}M\mathbf {V} _{i}\cdot \mathbf {V} _{i}+{\tfrac {1}{2}}{\boldsymbol {\omega }}_{i}\cdot [I_{R}]{\boldsymbol {\omega }}_{i}\right),}$$

^[6]

$${\displaystyle Q_{j}^{*}=-\left({\frac {d}{dt}}{\frac {\partial T}{\partial {\dot {q}}_{j}}}-{\frac {\partial T}{\partial q_{j}}}\right),\quad j=1,\ldots ,m.}$$

Динамическое равновесие

Форма принципа виртуальной работы Даламбера гласит, что система твердых тел находится в динамическом равновесии, когда виртуальная работа суммы приложенных сил и сил инерции равна нулю при любом виртуальном перемещении системы. Таким образом, для динамического равновесия системы n твердых тел с m обобщенными координатами необходимо, чтобы

$${\displaystyle \delta W=\left(Q_{1}+Q_{1}^{*}\right)\delta q_{1}+\dots +\left(Q_{m}+Q_{m}^{*}\right)\delta q_{m}=0,}$$

δq _jm

$${\displaystyle Q_{j}+Q_{j}^{*}=0,\quad j=1,\ldots ,m,}$$

$${\displaystyle {\frac {d}{dt}}{\frac {\partial T}{\partial {\dot {q}}_{j}}}-{\frac {\partial T}{\partial q_{j}}}=Q_{j},\quad j=1,\ldots ,m.}$$

Уравнения Лагранжа

Если обобщенные силы Q _j выводятся из потенциальной энергии V ( q ₁ , ..., q _m ) , то эти уравнения движения принимают вид

$${\displaystyle {\frac {d}{dt}}{\frac {\partial T}{\partial {\dot {q}}_{j}}}-{\frac {\partial T}{\partial q_{j}}}=-{\frac {\partial V}{\partial q_{j}}},\quad j=1,\ldots ,m.}$$

В этом случае введите лагранжиан L = T − V , чтобы эти уравнения движения приняли вид

$${\displaystyle {\frac {d}{dt}}{\frac {\partial L}{\partial {\dot {q}}_{j}}}-{\frac {\partial L}{\partial q_{j}}}=0\quad j=1,\ldots ,m.}$$

уравнения движения Лагранжа

Линейный и угловой момент

Система частиц

Линейный и угловой момент жесткой системы частиц формулируется путем измерения положения и скорости частиц относительно центра масс. Пусть система частиц P _i , i = 1, ..., n расположена в координатах r _i и скоростях v _i . Выберите опорную точку R и вычислите векторы относительного положения и скорости,

$${\displaystyle \mathbf {r} _{i}=\left(\mathbf {r} _{i}-\mathbf {R} \right)+\mathbf {R} ,\quad \mathbf {v} _{i}={\frac {d}{dt}}(\mathbf {r} _{i}-\mathbf {R} )+\mathbf {V} .}$$

Полные векторы линейного и углового момента относительно опорной точки R равны

$${\displaystyle \mathbf {p} ={\frac {d}{dt}}\left(\sum _{i=1}^{n}m_{i}\left(\mathbf {r} _{i}-\mathbf {R} \right)\right)+\left(\sum _{i=1}^{n}m_{i}\right)\mathbf {V} ,}$$

$${\displaystyle \mathbf {L} =\sum _{i=1}^{n}m_{i}\left(\mathbf {r} _{i}-\mathbf {R} \right)\times {\frac {d}{dt}}\left(\mathbf {r} _{i}-\mathbf {R} \right)+\left(\sum _{i=1}^{n}m_{i}\left(\mathbf {r} _{i}-\mathbf {R} \right)\right)\times \mathbf {V} .}$$

Если R выбран в качестве центра масс, эти уравнения упрощаются до

$${\displaystyle \mathbf {p} =M\mathbf {V} ,\quad \mathbf {L} =\sum _{i=1}^{n}m_{i}\left(\mathbf {r} _{i}-\mathbf {R} \right)\times {\frac {d}{dt}}\left(\mathbf {r} _{i}-\mathbf {R} \right).}$$

Жесткая система частиц

Чтобы специализировать эти формулы для твердого тела, предположим , что частицы жестко связаны друг с другом, поэтому Pi _, i=1,...,n расположены по координатам r _i и скоростям vi _. Выберите опорную точку R и вычислите векторы относительного положения и скорости,

$${\displaystyle \mathbf {r} _{i}=(\mathbf {r} _{i}-\mathbf {R} )+\mathbf {R} ,\quad \mathbf {v} _{i}=\omega \times (\mathbf {r} _{i}-\mathbf {R} )+\mathbf {V} ,}$$

^{[7] [8] [9]}

Линейный момент и угловой момент этой твердой системы, измеренные относительно центра масс R, равны

$${\displaystyle \mathbf {p} =\left(\sum _{i=1}^{n}m_{i}\right)\mathbf {V} ,\quad \mathbf {L} =\sum _{i=1}^{n}m_{i}(\mathbf {r} _{i}-\mathbf {R} )\times \mathbf {v} _{i}=\sum _{i=1}^{n}m_{i}(\mathbf {r} _{i}-\mathbf {R} )\times (\omega \times (\mathbf {r} _{i}-\mathbf {R} )).}$$

Эти уравнения упрощаются и становятся

$${\displaystyle \mathbf {p} =M\mathbf {V} ,\quad \mathbf {L} =[I_{R}]\omega ,}$$

_]момента инерции,

$${\displaystyle [I_{R}]=-\sum _{i=1}^{n}m_{i}[r_{i}-R][r_{i}-R],}$$

₋r _iR.

Приложения

• Для анализа робототехнических систем
• Для биомеханического анализа животных, людей или гуманоидных систем.
• Для анализа космических объектов
• Для понимания странных движений твердых тел. ^[10]
• За проектирование и разработку датчиков, основанных на динамике, таких как гироскопические датчики.
• За проектирование и разработку различных приложений для повышения устойчивости автомобилей.
• Для улучшения графики видеоигр, в которых используются твердые тела.

Смотрите также

• Аналитическая механика
• Аналитическая динамика
• Вариационное исчисление
• Классическая механика
• Динамика (механика)
• История классической механики
• Лагранжева механика
• лагранжиан
• гамильтонова механика
• Жесткое тело
• Жесткий ротор
• Динамика мягкого тела
• Многотельная система
• Полходе
• Герпоход
• Прецессия
• Эллипсоид Пуансо
• Гироскоп
• Физический движок
• Физический процессор
• Уровень абстракции физики – унифицированный симулятор многотельных тел.
• RigidChips — японский симулятор твердого тела.
• Уравнение Эйлера

Рекомендации

1. Б. Пол, Кинематика и динамика планарных машин, Прентис-Холл, Нью-Джерси, 1979.
2. Л.В. Цай, Анализ роботов: механика последовательных и параллельных манипуляторов, Джон-Уайли, Нью-Йорк, 1999.
3. Британская энциклопедия, Законы движения Ньютона.
4. К. Дж. Уолдрон и Г. Л. Кинзель, Кинематика и динамика и проектирование машин, 2-е изд., John Wiley and Sons, 2004.
5. Торби, Брюс (1984). «Энергетические методы». Advanced Dynamics для инженеров . Серия HRW в области машиностроения. Соединенные Штаты Америки: Издательство CBS College Publishing. ISBN 0-03-063366-4.
6. Т.Р. Кейн и Д.А. Левинсон, Динамика, теория и приложения, McGraw-Hill, Нью-Йорк, 2005.
7. Мэрион, JB; Торнтон, Северная Каролина (1995). Классическая динамика систем и частиц (4-е изд.). Томсон. ISBN 0-03-097302-3..
8. Саймон, КР (1971). Механика (3-е изд.). Аддисон-Уэсли. ISBN 0-201-07392-7..
9. Тененбаум, РА (2004). Основы прикладной динамики. Спрингер. ISBN 0-387-00887-Х..
10. Гомес, RW; Эрнандес-Гомес, Джей-Джей; Маркина, В. (25 июля 2012 г.). «Прыгающий цилиндр на наклонной плоскости». Евро. Дж. Физ . ИОП. 33 (5): 1359–1365. arXiv : 1204.0600 . Бибкод : 2012EJPh...33.1359G. дои : 10.1088/0143-0807/33/5/1359. S2CID  55442794 . Проверено 25 апреля 2016 г.

дальнейшее чтение

• Э. Лейманис (1965). Общая задача о движении связанных твердых тел вокруг неподвижной точки. ( Спрингер , Нью-Йорк).
• ВБ Херд (2006). Механика твердого тела: математика, физика и приложения. ( Вайли-ВЧ ).

Внешние ссылки

• Информация Криса Хекера о динамике твердого тела, заархивированная 12 марта 2007 г. в Wayback Machine.
• Физически обоснованное моделирование: принципы и практика
• База знаний DigitalRune. Архивировано 20 ноября 2008 г. на Wayback Machine. Содержит магистерскую диссертацию и коллекцию ресурсов по динамике твердого тела.
• Ф. Кляйн, «Заметка о связи линейной геометрии и механики твердых тел» (английский перевод)
• Ф. Кляйн, «О теории винтов сэра Роберта Болла» (английский перевод)
• Э. Коттон, «Применение геометрии Кэли к геометрическому исследованию смещения твердого тела вокруг фиксированной точки» (английский перевод)