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

Мехатроника и динамика мини-роботов

.pdf
Скачиваний:
2
Добавлен:
24.11.2025
Размер:
4.17 Mб
Скачать

f1 f2 f3 f – коэффициенты трения колес; M N – момент, вращающий колесо 1z;

M5 – момент управляющий колесом 2z;

M0 – момент сопротивления в паре трения «колесо-дорога»,

возникающего при повороте колеса 2z;

l, l1, l2 – соответствующее расстояние, вытекающее из геомет-

рии системы;

r1 r2 r3 r – радиусы соответствующих колес.

Динамический анализ проведем на основе уравнения

 

 

 

T

(q) .

(19.5)

M (q)q

C(q,q)q B(q) J

 

Выступающие в (19.5) матрицы М, С, В следуют из уравнения (19.4). Дальше принимаем, что m первых столбцов матрицы J (q)

(19.3) образуют матрицу

J1(q) размером m m . Разложим вектор

координат q :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q [q , q ]T ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

 

 

 

 

 

 

где q Rm , q Rn m.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда уравнения связей запишем в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[J (q), J

2

(q)]

 

q1

 

0,

d e t J (q) 0,

 

1

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

q2

 

 

 

 

 

 

 

 

 

 

откуда следует

 

 

 

q1 J12

(q) q2 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

r cos ,

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r sin ,

0

 

J12 (q) J1 1(q)J2 (q) r

tg ,

0

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

0

 

 

 

 

 

 

 

 

 

 

cos

 

 

 

460

Тогда

q JI12 (q) q2 ,

n m

где In m – единичная матрица.

Векторы скорости и ускорения можно также записать в виде

q T (q) E2q2;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q T (q) E2 q2 T (q, q) E2 q2;

 

 

 

 

 

 

T (q)

I

m

J (q)

 

 

 

12

.

 

 

0

 

In m

 

Используя также зависимости, вытекающие из проведенных преобразований, запишем уравнения (19.5) в виде

 

 

 

 

T

(q) ;

(19.6)

M12 (q)q2

C12 (q, q2 )q2

B1(q) J1

M22 (q)q2 C22 (q, q2 )q2 B2 (q) ;

(19.7)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q1

J12 (q)q2 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

T

 

T

 

 

T

 

 

M

 

E

T

T

T

M T E ;

 

 

(T

 

T

CT )E2

;

12

 

 

 

C12 E1

 

MT

 

 

1

 

 

 

 

2

 

 

 

T

 

T

 

 

T

 

 

M

 

E

T

T

T

M T E

;

 

 

(T

 

T

CT )E2

;

22

2

 

C22 E2

 

MT

 

 

 

 

 

 

 

2

 

B2 E2T TT B;

 

 

 

 

B

ET TT B;

 

 

 

 

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

 

 

Im

 

Rn m ,

 

E

0

 

Rn n m) .

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

In m

 

 

 

 

 

 

Матрица Jm – единичная.

Рассмотренные выше методы позволяют решать прямые и обратные задачи динамики.

461

В дальнейшем анализе полагаем

q1 [xA, yA, , ]T , q2 [ , ]T ,

тогда уравнение (19.7) запишем в виде

 

2

 

1

 

 

 

(m1 m2 m3 m5 )r

 

Iz1 Iz2 Iz3

cos

2

 

Ix3

 

 

 

 

 

 

Ix1 Ix2 Ix3 2Iz1h12 Iz5 (m1 m2 )l12 m3l2 m5l22

rl tg

2

rl2 tg2 (19.8)

 

 

I

 

 

I

 

I

 

2I

h2

I

 

 

(m

m )l

2

m l2

m l2

r

2

I

 

 

 

 

 

x1

x2

x3

z5

 

 

z3

 

 

 

 

2

 

 

 

 

 

 

 

 

z1 1

 

 

 

1

2 1

 

3

5 2

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tg

 

M N N1 f1

N2 f2

N3 f3

 

1

 

(Ms M0 sgn

 

r

 

 

 

cos2

 

cos

)

l

tg ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

r

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ix3

 

 

 

tg Ix3 Ix3

 

 

 

 

 

Ms M0 sgn ,

 

(19.9)

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

l

cos2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

уравнения (19.8), а также (19.9) позволяют провести анализ прямых и обратных задач динамики без знания множителей Лагранжа. Уравнения же (19.6) запишем в следующей форме:

 

 

 

 

 

 

 

 

 

 

 

 

 

r 2

 

 

 

 

 

 

 

 

 

(m1 m2 m3

m5 )r( cos

 

 

 

 

 

tgφsin )

(m3l m5l 2 )

 

l

 

 

 

r

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

sinβ

 

 

(19.10)

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

( tgφcos

l

 

 

 

tg

 

φsin φ

cos2

φ

) = λ1

+ λ3;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

2

 

 

 

 

 

 

 

 

 

 

(m1 m2 m3

m5 )r( sin

l

 

tg cos )

(m3l m5l 2 )

 

r

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

cos

 

 

(19.11)

 

 

 

 

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( tg cos

 

 

 

tg

sin

 

 

 

) 2

4 ;

 

l

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

462

Ix1 Ix2 Ix3 2Iz1h12 Iz5 (m1 m2 )l12 m3l2 m5l22

 

 

r

 

 

 

1

 

) Ix3 (m3l m5l2 )

r2

2

tg

(

l

 

cos2

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ms M0 sgn 3l sin 4l cos ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

sin

 

 

 

 

 

Iz3

 

 

 

 

Iz3

 

N3 f3

r 3 cos( )

 

 

cos

 

 

 

cos2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r 4 sin( ).

(19.12)

(19.13)

Из системы уравнения (19.10)–(19.13) можно определить множители Лагранжа, т. е. определить, как в течение анализируемого состояния изменяются составляющие сил сухого трения, возникающих в плоскостях взаимодействия колес и поверхности.

Рис. 19.2. Модель трехколесного робота: вид сверху и сбоку

463

19.2. Динамика модели четырехколесного миниатюрного робота

Рассмотрим модель четырехколесного мини-робота, изображенного на рис. 19.3. Два передних колеса 3 и 4 и задних 1, 2 заменяем эффективными колесами 2z и 1z соответственно, расположенными посередине осей в точках А и F. В качестве обобщенных координат берем, как и в случае трехколесного робота: xA, yA, , , , .

Уравнения неголономных связей, наложенных на робот, имеют вид (19.1). Уравнения Лагранжа II рода записываются в виде (19.2).

Определение обобщенных сил Q , кинетической и потенциальной энергии также проводится аналогично случаю трехколесного робо-

та. Предположим,

что m1 m2 ,

m3 m4 ,

Ix1 Ix2 ,

Ix3 Ix4 ,

Iz1 Iz2 Iz3 Iz4 ,

M03 M04 – моменты сопротивления при по-

вороте рулевых колес 3 и 4.

 

 

 

Рис. 19.3. Модель четырехколесного робота

464

Уравнения динамики четырехколесного робота имеют вид

2(m1 m3 12 m5 )xA (2m3l m5l2 )( sin 2 cos ) 1 3;

(19.14) 2(m1 m3 12 m5 ) yA (2m3l m5l2 )( sin 2 cos ) 2 4.

В реальной конструкции колеса 3 и 4 снабжены отдельными приводами (актуаторами) так, что моменты M3, M4 , прикладывае-

мые к ним, различные. Предположим, что мощности двигателей актуаторов равны и имеет место равенство

 

 

(19.15)

M3 3

M4 4.

Учитывая зависимости, следующие из кинематики робота:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r1

 

r2

 

1

h1 1;

 

 

 

2

h1 1;

h1

l1

 

l2

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

l1(l1 sin 2l cos )sin ;

 

 

 

 

3

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

l1(l1 sin 2l cos )sin ;

 

 

(19.16)

 

4

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

l2

 

l1(l1 sin 2l cos )sin ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

2

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l2

l1(l1 sin 2l cos )sin

 

 

 

 

 

и подставляем 3, 4 из (19.16) в (19.15), получим

M

 

M

W3 .

(19.17)

 

4

 

3 W

 

 

 

 

4

 

465

В реальной конструкции колёса также управляются независимыми двигателями актуаторов. Моменты, управляющие этими колесами обозначаются M35, M45 . Полагаем, что мощности двигате-

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

M35 3

M45 4.

(19.18)

 

 

 

С учетом уравнений (19.16) из (19.18) получим соотношение

M

45

M

W1 .

(19.19)

 

 

35 W

 

 

 

 

2

 

Предполагаем, что давление на колеса, взаимодействующие с опорной поверхностью, имеет статический характер и определено при предположении, что давление на передние колеса равно давлению на задние (давления на оси). Решая системы уравнений (19.14)– (19.19), можем проводить анализ прямых и обратных задач динамики мини-робота.

Характеризуя обратные задачи динамики, т. е. зная параметры движения, интерпретируем значениями моментов вращающих и управляющих. Как видно из уравнений (19.14), их решение представляет трудности, поэтому, как и в случае трехколесного робота, запишем (19.14) в векторно-матричном виде:

 

 

T

(q) .

(19.20)

M (q) q

C(q, q) q B(q) J

 

Вид матриц M, C, B, J вытекает из уравнений (19.14). Проводя аналогично тому, как это было сделано для трехколесного робота, преобразования, отделяющие множители Лагранжа от моментов, запишем уравнение (19.20) в виде

 

 

 

 

T

(q) ;

(19.21)

M12 (q)q2

C12 (q, q2 )q2

B1(q) J1

M22 (q)q2

C12 (q, q2 )q2 B2 (q) ;

 

 

 

 

 

 

 

 

466

q1 J12 (q)q2;

q1 [xA, yA, , ]T ; q2 [ , ]T .

Расписывая уравнения (19.20), (19.21) в координатах аналогично тому, как это было сделано для трехколесного робота, получим для (19.20)

 

 

 

 

 

(2m1 2m3

m5)r

2

2Iz1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cos

 

2Iz1

 

 

 

 

2Ix1 2Ix3 2Iz1h12 Iz5 2m1l12 2m3(l2 l12 ) m5l22

r2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

l

2

 

 

2

I

 

lr(w w )sin (2m

 

 

 

 

 

2

 

2I

 

1

 

sin

 

2m m ) r

 

 

z3

l

2

 

 

 

 

 

x3

 

1

 

2

 

 

 

 

 

1

 

3

 

5

 

 

 

 

 

 

 

 

 

 

cos 2Ix1 2Ix3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

Iz5

 

2

 

 

 

 

 

2Iz1 sin

2Iz1h1

2m1l1

 

 

(19.22)

 

 

 

 

2

2

 

 

2

 

r2

 

 

 

 

 

 

l

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

2m3(l

 

l1 )

m5l2

 

2

sin cos 2Iz3

 

2

 

sin cos

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

2

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w

 

 

 

 

 

 

 

2Ix3 l l1 r

(w1 w3 w2 w4 )sin (M3 N3 f3) l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

(M4 N4 f4 ) l 1w2 (N1 f1 N2 f2 )cos

(M3s M30 sgn )l2w1 (M4s M40 sgn )l2w2 rl sin ;

 

 

 

 

 

2

2

 

Ix3lr (w1 w2 )sin Ix3l4 (w1

w2 )

 

 

l2

 

2

 

 

 

 

1

 

 

 

 

 

2Iz3

 

 

 

sin cos Ix3lr(w1

w2 ) cos

(19.23)

l2

 

 

 

 

 

 

 

 

 

2Ix3l4l1 2 (w13w3 w23w4 )

(M3s M30 sgn )l2w1 (M4s M40 sgn )l2w2.

467

Для уравнения (19.37) получим соответственно

 

(2m1

2m3 m5 )rcos β cos φ (2m3l m5l2 )

r

 

 

 

sin β sin φ

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(19.24)

(2m1 2m3 m5 )r( β sin β cos φ φ cos β sin φ)

 

 

 

 

 

 

 

 

 

 

 

(2m3l m5l2 )

r

 

 

 

 

 

 

 

 

(β cos β sin φ + φ sin β cos φ) 1 3;

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

(2m1 2m3

m5 )r sin cos (2m3l m5l2 )

 

 

l

cos sin

 

 

 

 

 

 

 

 

 

 

 

(2m1

 

 

 

 

 

 

 

 

(19.25)

 

2m3 m5 )r( cos cos sin sin )

 

 

 

 

 

 

 

 

 

 

 

(2m3l m5l2 ) rl ( cos sin sin sin ) 2 4;

2Ix1 2Ix3 2Iz1h12 Iz5 2m1l12 2m3 (l2 l12 ) m5l22

rl sin Ix3l2 (w1 w2 ) (2m3l m5l2 )r cos

2Ix1 2Ix3 2Iz1h12 m1l12 2m3 (l2 l12 ) m5l22 rl cos

2

2

2

2

 

 

 

2

w1

 

(19.26)

2Ix3l

l1

(w1 w3

w2 w4 )

(M3s M30 sgn )l

 

 

 

 

 

2

w2

3lsin 4l cos ;

 

 

(M4s M40 sgn )l

 

 

 

2Iz1 cos 2Iz1 sin N1 f1 N2 f2 1r cos 2r sin .

 

 

 

 

 

 

 

 

 

 

 

 

Структура уравнений (19.22) и (19.23) с учетом (19.18), (19.19)

позволяет анализировать прямые и обратные задачи динамики. Из уравнений (19.24), (19.25), (19.26) определяем значения множителей Лагранжа, т. е. составляющие сил сухого трения, действующих в плоскостях касания колес 1z и 2z с покрытием, чтобы реализовать программное движение мобильного колесного мини-робота:

468

2m1l12 2m3 (l2 l12 ) m5l22 2Ix1 2Iz1h12 2Ix3 Iz5 (2m3l m5l2 )

( xA sin y cos ) Ix3l2 (w1 w2 ) 2Ix3l2l1 2 (w12w3 w22w4 )

(M3s M30 sgn )w1l2 (M4s M40 sgn )w2l2 3l sin 4l cos ;

 

 

 

 

 

 

 

 

 

 

 

 

 

N2 f2 1r cos 2r sin ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2Iz1 N1 f1

 

 

 

 

 

 

 

 

 

 

 

 

l2

l2 sin2

 

 

 

 

 

 

 

l2

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

1

 

 

 

4Iz

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2Iz3

 

 

 

l2

 

 

 

3 l2 sin cos (M3 N3 f3 ) l

w

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

(M4 N4 f4 )

 

1

 

 

3r cos( ) 4r sin( );

 

 

 

 

 

 

 

 

 

l w

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

w1 2l

2

 

2

 

2

 

 

 

2

 

 

2

w2 2l

2

2

2

 

 

2

 

 

Ix3 ( l

 

 

l1

w1 w3 )w1l

 

 

Ix3 ( l

 

 

l1

w2 w4 )w2l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

2

 

1

 

 

 

sin cos (M3s M30 sgn )w1l

(M4s M40 sgn )w2l

 

2Iz3

 

 

 

 

 

 

;

l2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w1 l2

 

l1(l1 sin 2l cos )sin 1 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w2 l2

l1(l1 sin 2l cos )sin 1 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w3 l (l1 sin 2l cos )cos ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w4 l (l1 sin 2l cos )cos .

 

 

 

 

 

 

 

 

 

 

19.3. Динамика модели двухколесного миниатюрного робота

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

Рассмотрим модель, когда два колеса 1 и 2 расположены на одной оси, а третьеопорное 3, как в схеме трехколесного робота, рис. 19.4.

469