|
X1 |
|
|
X1 |
|
|
|
|
Y |
|
A |
Y |
; |
|
|
|
|
1 |
|
2 |
|
2 |
|
|
|
|
|
|
|
|
|
z |
|
|
|
ZK |
|
|
|
|
|
|
|
1 |
|
0 |
|
|
0 |
|
A |
|
|
0 |
|
cos |
|
sin |
. |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
sin |
|
cos |
|
|
|
|
|
|
|
Рис. 9.4. Переход от промежуточной системы 0X1Y1ZK к системе 0X KY2z
Здесь M 2 – проекция точки M в плоскости 0Y1ZK .
3. Третий поворот вокруг оси 0z переводит промежуточную систему координат 0X1Y2 z в основную с помощью матрицы A3, рис. 9.5:
|
X1 |
|
|
|
x |
|
|
|
|
|
Y |
|
A |
y |
|
; |
|
|
|
|
2 |
|
|
3 |
|
|
|
|
|
|
|
z |
|
|
|
|
|
|
|
|
|
|
|
|
|
z |
|
|
|
|
|
cos |
sin |
0 |
|
A |
sin |
cos |
|
0 |
. |
3 |
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
0 |
|
|
1 |
|
|
|
|
|
|
|
|
|
Рис. 9.5. Третий поворот вокруг оси 0z
Тогда матрицу A перехода от системы 0xyz к системе 0X KYK ZK представим в виде
A A1A2 A3 ,
а ее элементы имеют вид
a11 cos cos sin sin cos ; |
|
a12 |
cos sin sin cos cos ; |
|
|
a13 sin sin ; |
|
a21 sin cos cos sin cos ; |
|
a22 |
sin sin cos cos cos ; |
(9.3) |
|
a23 cos sin ; |
|
a31 sin sin ; a32 cos sin ; a33 cos .
Как следует из (9.3), при 0, линия узлов не определена
и углы , не определены, а определена только их сумма |
. |
|
261 |
Следовательно, при движениях твердого тела, когда ось 0z мало отклоняется от 0ZK , целесообразно применять другие углы, опре-
деляющие ориентацию тела в пространстве.
Отметим также, что в силу того, что произведение матриц A A1A2 A3 некоммутативно, то и конечные повороты твердого тела
некоммуникативны. Это означает, что в общем случае ориентация твердого тела (мини-робота), получаемая им в результате двух последовательных конечных поворотов, зависит от порядка выполнения этих поворотов. Для мини-роботов, перемещающихся в автономном режиме, алгоритм подачи управляющих воздействий должен разрабатываться с учетом этого.
Отметим некоторые свойства матрицы A, следующие из ее ортогональности:
1.A 1 AT, индекс T означает транспонирование;
2.AAT E – единичная матрица;
3.det A 2 1.
Если зависимость A t от времени непрерывная, то знак A o 1 сохраняется в течение движения.
9.1.1. Мгновенное кинематическое состояние твердого тела
Рассмотрим свободное движение твердого тела, скорости точек которого характеризуются вектором
V V r1 t r1 t ,
зависящим от радиуса-вектора r1 t в неподвижной системе коорди-
нат времени t .
Если в момент времени t
V r1 t Vo t const,
то тело совершает мгновенное поступательное движение в мо-
мент времени t , а при V0 t 0 находится в мгновенном покое. Очевидно, что ускорения
W V r1 t r1 t
могут быть разными в момент времени t в разных точках тела. Если в данный момент времени t в теле (мини-роботе) имеются
две неподвижные точки, то тело совершает мгновенное вращательное движение вокруг прямой, проходящей через эти точки, называемой мгновенная ось вращения. Мгновенная ось вращения меняет свое положение в неподвижной и подвижной системах координат при конечном движении тела.
Если в данный момент времени тело (мини-робот) участвует в двух мгновенных движениях: мгновенном поступательном вдоль некоторой оси имгновенном вращательном вокруг этойоси, то тело (мини-робот) совершает мгновенное винтовое движение, которое является самым общиммгновенным движением свободного твердого тела.
9.1.2.Скорость и ускорение точек твердого тела
вобщем случае движения
Дифференцируя (9.1) по времени t , получим
|
|
|
|
|
|
|
|
|
|
r , |
(9.4) |
|
|
V |
V |
r V |
0 |
0 |
|
|
|
|
|
|
|
|
|
|
где V0 – скорость полюса 0; |
|
|
|
|
|
– угловая скорость вращения системы 0xyz |
вокруг мгновен- |
ной оси вращения, проходящей через 0 .
Формула (9.4) справедлива в силу формулы Эйлера
r r.
Ускорение произвольной точки M твердого тела (мини-робота) находится дифференцированием по времени t формулы (9.4):
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
W0 r r W0 W Wn , (9.5) |
W V0 |
r r |
где W0 – ускорение точки 0;
W r – вращательное ускорение; Wn – осестремительное ускорение;
– угловое ускорение.
9.1.3.Угловая скорость в подвижной и неподвижной системах координат
Вектор мгновенной угловой скорости , направленный вдоль мгновенной оси вращения, имеет компоненты x , y , z , кото-
рые являются проекциями на оси подвижной системы координат 0xyz . С другой стороны, можно представить в виде суммы трех
|
|
|
|
|
угловых скоростей вращения вокругосей 0ZK ,0N ,0z : |
|
|
|
|
(9.6) |
. |
|
Записывая (9.6) покомпонентно, получим кинематические
формулы Эйлера
x sin sin cos ;y sin cos sin ;
z cos .
9.1.4. Дискретная аппроксимация твердого тела. Геометрия масс твердого тела и эквивалентной системы материальных точек
При решении многих конкретных задач удобно заменять сплошное твердое тело эквивалентной системой материальных точек,
рис. 9.6.
Рис. 9.6. Замена сплошного твердого тела эквивалентной системой материальных точек
Разбиваем тело на конечные объемы (элементы), массы которых
сосредоточим |
в |
точках, радиусы-векторы которых |
обозначим |
r1 m m |
|
, |
r m – вектор, соединяющий точку 0 и M в системе |
1,n |
0X Y Z , m |
– вектор |
|
|
в системе 0xyz. |
|
OM |
|
1 1 |
1 |
|
|
|
|
|
|
|
|
|
|
|
Центр масс rc |
твердого тела вычисляется по формуле |
|
|
|
|
|
|
r |
|
Vrdm |
|
Vr r dr |
, |
(9.7) |
|
|
|
|
|
|
|
|
|
|
|
|
c |
|
Vdm |
|
V r d r |
|
где r ddm ; m – масса тела;
V– объем;
– плотность.
265
Проекцируя (9.5) на оси координат, получим три уравнения. Центр масс эквивалентной системы материальных точек вычис-
ляется по формуле
(9.8)
n
m m ,
1
где m – масса -го элемента объема; m – масса всей системы.
При n разность между rc , вычисляемому по формуле (9.8), стремится к rc , вычисляемому по формуле (9.7).
Центр масс с может быть выбран в качестве полюса 0, что будет использовано дальше.
Моменты инерции сплошного твердого тела относительно осей x, y, z вычисляются по формулам
Ix y2 z2 dm, |
I y x2 z2 dm, |
Iz x2 y2 dm, |
(9.9) |
V |
V |
|
V |
|
а дискретной системы материальных точек – по формулам |
|
|
n |
; |
|
|
|
Ix m y2 z2 |
|
|
|
1 |
|
|
|
|
n |
; |
|
|
I y m x2 z2 |
(9.10) |
|
1 |
|
|
|
|
n |
. |
|
|
|
Iz m x2 y2 |
|
|
|
1 |
|
|
|
Центробежные моменты твердого тела вычисляются – по формулам
Ixy xydm; |
Ixz xzdm; |
I yz yzdm, |
(9.11) |
|
V |
V |
V |
|
а для дискретной системы материальных точек – по формулам |
|
n |
Ixz |
n |
n |
|
Ixy m x y ; |
m x z ; |
I yz m y z . |
(9.12) |
1 |
|
1 |
1 |
|
Величины Iij i, j x, y, |
z зависят от выбора тоски 0 и ориента- |
ции осей координат x, y, z. При n формулы (9.10), (9.12) пе-
реходят в формулы (9.9), (9.11) соответственно.
Компоненты Iij образуют симметричный тензор второго ранга
(тензор инерции), который характеризует сопротивление твердого тела вращательным движением:
|
Ix |
Ixy |
Ixz |
|
Ixy |
I y |
|
|
|
I |
I yz . |
|
I |
|
I |
|
I |
|
|
|
xz |
yz |
z |
|
|
|
|
|
|
Моменты инерции относительно любых параллельных осей вы-
числяются по формуле |
|
|
Iu I m du2 d 2 , |
(9.13) |
где Iu |
– момент инерции относительно оси u; |
|
I |
– момент инерции относительно параллельной оси ; |
|
du , d – расстояния осей u, от оси, проходящей через центр масс.
9.2. Основные законы динамики твердого тела (мини-робота)
В динамике твердого тела основные законы представляют собой уравнения баланса относительно основных динамических величин мехатронной системы, какими являются:
1)векторные величины: импульсы (количество движения), момент количества движения (момент импульса), силы, моменты сил;
2)скалярные величины: кинетическая и потенциальная энергия, работа, диссипация энергии, энтропия.
Предварительно напомним некоторые определения. Количество движения (импульс) твердого тела – вектор Q, вы-
числяется по формуле
Импульс дискретной мехатронной системы вычисляется по формуле
n
Q m V .
1
Кинетический момент K (момент импульса, момент количества движения) твердого тела (мини-робота) относительно точки A вычисляется по формуле
KA AQdr,
V
где A – радиус-вектор точек тела относительно точки A. Проекции вектора KA на оси координат называются моментами
импульса относительно осей x, y, z:
|
|
A |
i |
, |
|
|
A |
j |
, |
|
|
|
|
|
|
Kx K |
Ky K |
Kz K |
Ak , |
(9.14) |
где i , j , k – орты осей координат.
Кинетический момент (момент импульса) относительно точки А дискретной системы вычисляется по формуле
n
KA AQ.
1
Моменты импульса относительно осей x, y, z вычисляются по
формулам (9.14).
При изменении центра A, относительно которого вычисляется кинетический момент, момент относительно нового центра B вычисляется по формуле
где вектор BA – расстояние между точками A и B.
В векторно-матричном виде кинетический момент твердого тела (системы), вращающегося вокруг точки 0, записывается в виде
K0 I , ωx , ωy , ωz p, q, r ,
а в развернутом виде получим
Kox Ix p Ixyq Ixz r;
Koy Ixy p I y q I yz r;
Koz Ixz p I yz q Iz r.
Закон об изменении импульса (количества движения)
ференциальной форме имеет вид
ddQt R e