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

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

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

Из структурной схемы получим

 

 

 

p2

2d0 p 02 02 з p 2d0 1,

(18.33)

где 2d

0

k ;

2 k

0

.

 

 

 

1

 

0

 

 

 

Передаточные числа k0 и k1

выбираются из условия получения

заданных значений частоты Ω0

и коэффициента затухания d.

Системе (18.32) свойственна статическая ошибка, определяемая главным образом членом 2d0 1.

При формировании закона управления (18.31) используется измеренное значение дистанции , а производная p может быть

получена дифференцированием сглаженного сигнала . Из структурной схемы (рис. 18.12) имеем уравнение

p2 2d0 p 02 02 з p 1.

Рис. 18.12. Структурная схема системы управления позиционной дистанцией

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

18.3. Динамика систем автоматического управления группой

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

450

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

Рассмотрим некоторые случаи автоматического управления дистанцией между четырьмя дронами в строю при условии, что тяга двигателя устанавливается в соответствии с законом управления без запаздывания. Пусть дроны движутся строем «кильватер» и последующий дрон ориентируется в строю по предыдущему. Тогда уравнения строя в соответствии с системой (18.6)–(18.9) будут

p 12 WV p p1 p2

f12 ;

 

 

 

 

p 23

WV p p2

p3

 

f23

 

 

(18.34)

;

p

 

W

p

 

 

 

 

 

 

f

 

 

 

 

34

p3

p4

34

,

 

 

V

 

 

 

 

 

 

 

где f12 FV12 и т. д.

Уравнение автомата тяги (автомата дистанции) возьмем в виде

 

δp W p λ k0, p λ3 ,

 

 

 

 

 

 

 

(18.35)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W p k

 

 

 

k

0

k p k

2

p2

 

 

, k

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

, k

 

 

k

 

(18.36)

д p 1

 

 

 

p 1

 

 

 

 

 

 

2

0

 

 

 

0

 

д

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

для статического регулятора и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W p

k

д

 

 

k

0

k p k

2

p2

 

k ,

 

k

д

 

 

 

 

 

 

 

 

 

1

 

 

 

 

,

 

 

 

k

 

(18.37)

p

 

 

p 1

 

p 1

p

 

 

 

2

 

0

 

 

0

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

для астатического регулятора; kд – коэффициент усиления;

k0 , k1 , k2 – периодические числа;

1 и 2 – постоянные времени, характеризующие несовершен-

ство измерительных и исполнительных устройств автомата дистанции, т. е. передачу полезных сигналов через инерционные звенья.

451

Будем полагать, что дистанция λ измеряется, а производные pλ

и p2λ получаются дифференцированием сигнала λ.

Если в структуре закона управления взять сигналы, пропорциональные и p , тогда, учитывая наличие погрешностей τ1 и 2

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

Решая совместно уравнения (18.34) и (18.35), найдем

 

 

 

 

12 f12 ;

 

 

 

 

 

 

2

23

W W f f

23

;

 

 

 

V

12

 

 

3

34

W 2W 2 f

W W

23

 

 

V

12

V

 

 

где p WVW;

WV ( p1 k0kд з).

 

 

 

 

(18.38)

2 f34 ,

Из уравнений (18.38) следует, что возмущения, вносимые в систему управления изменением силы тяги ведущего дрона p1 или перена-

стройкой регулятора λ3 , передаются через систему как через фильтр, состоящий из одинаковых ячеек. Заметим, что заданная дистанция λ3

может бытьустановленаиизмененатолько наведущем дроне. Уравнения (18.38) являются чрезвычайно сложными и их анали-

тическое решение затруднено. На рис. 18.13–18.17 приводятся графики переходных процессов, полученные моделированием этих уравнений.

а

б

Рис. 18.13. Графики переходных процессов в статической системе управления при возмущении перенастройкой (а) и тягой (б)

452

а

б

Рис. 18.14. Графики переходных процессов в статической системе управления при возмущении перенастройкой (а) и тягой (б)

а

б

Рис. 18.15. Графики переходных процессов в статической системе управления при возмущении перенастройкой (а) и тягой (б)

а

б

Рис. 18.16. Графики переходных процессов в статической системе управления при возмущении перенастройкой (а) и тягой (б)

453

а б

Рис. 18.17. Графики переходных процессов в статической системе управления при возмущении перенастройкой (а) и тягой (б)

При моделировании уравнений менялись передаточные числа k0 , k1 и k2 закона управления с целью выбора приемлемых зна-

чений. Переходные процессы построены для дронов строя, имею-

щих Tc 50 с,

kc 0,5.

Постоянные времени 1 и 2

приняты

1 2 0,2 с.

Графики

построены для статической

системы

(см. рис. 18.13–18.15) с законом управления (18.36) и астатической системы (см. рис. 18.16 и 18.17) с законом управления (18.37). Зна-

чения передаточных чисел k0 , k1 и k2 , для которых построены

графики, даны в табл. 18.1 (статическая система) и табл. 18.2 (астатическая система).

Таблица 18.1 Передаточные числа для статической системы

Рисунок

k0

k1, c

k1, c2

Рис. 18.13

2

8

8

Рис. 18.14

1

2

0

Рис. 18.15

1

4

0

 

 

 

Таблица 18.2

Передаточные числа для статической системы

 

 

 

 

Рисунок

k0, с–1

k1

k2, c2

Рис. 18.16

1

4

4

Рис. 18.18

0,5

2

2

454

Проанализируем графики, характеризующие переходные процессы в системе управления строем дронов.

В статической системе управления (см. рис. 18.13–18.15) имеется статическая ошибка, которую можно уменьшить путем увеличения

передаточного числа k0 . Процессы в системе при перенастройке автомата на новую дистанцию λ3 и при возмущении тягой p1 ве-

дущего дрона могут быть неколебательными (см. рис. 18.13), если передаточное число по позиционному сигналу k0 2 при значи-

тельных передаточных числах по скорости k1 8 и ускорению k2 8. Реализовать такие большие передаточные числа по скорости

и ускорению затруднительно из-за помех.

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

T k1 / k0 должна быть T 10 c.

Астатическая система управления при прочих равных условиях более склонна к колебаниям (см. рис. 18.16 и 18.17), чем статическая, что в общем случае объяснимо. При увеличении передаточно-

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

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

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

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

455

Глава 19. ДИНАМИКА КОЛЕСНЫХ МИНИАТЮРНЫХ РОБОТОВ

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

Рассмотрим модель трехколесной платформы робота, изображенную на рис. 19.1. Движение рассматривается в плоскости 0XY. Робот со смонтированным на платформе оборудованием моделируется однородной прямоугольной пластиной, главные оси инерции которой совпадают с геометрическими главными центральными осями.

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

Колеса 1 и 2 действительной конструкции заменим эффективным колесом 1z, расположенным в точке А, середине оси. Колесо 3 обозначим колесом 2z.

Обобщенными координатами, описывающими движение платформы являются xA , yA , , , , , где xA , yA – координаты точ-

ки А, расположенной на оси посередине между двумя колесами, и проходят точках В, С.

456

Поворотное колесо с центром в точке F характеризуется углами

и углом , вращение рулевого колеса характеризуется углом , ,

аведомых колес – углом .

Кинематические соотношения для неголономных связей имеютвид

xA r1 cos 0;

yA r1 sin 0;

(19.1)

xA l sin r3 cos( ) 0;

yA l cos r3 sin( ) 0.

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

Уравнения Лагранжа II рода динамики неголономной системы запишем в виде

d

E T

 

E T

 

 

 

 

 

 

Q JT (q) ,

 

 

 

 

 

 

(19.2)

 

dt

q

 

q

 

 

 

 

 

 

где – q [xA , yA , , , , ] вектор обобщенных координат; Q – вектор обобщенных сил;

J (q) – реакции связей;

– вектор множителей Лагранжа, которые являются силами трения, лежащими в плоскостях контакта колес с плоскостью.

Кинематическая энергия Т имеет вид

E 6 12Tr(TioBiTioT ),

i 1

457

где В – матрица инерции i-й части;

Tio – матрица преобразования локальной i-й системы координат

к основной Oxyz.

Уравнения неголономных связей запишем в виде

J (q)q 0,

где J (q) имеет вид

 

 

 

 

 

 

1

0

0

0

r cos 0

 

 

 

 

 

 

 

 

J (q) 0

1

0

0

r sin 0 .

(19.3)

 

0

l sin

r cos( )

0

 

 

1

0

 

0

1

l cos

r sin( )

0

0

 

Вектор обобщенных сил Q определяется соотношениями

 

 

 

6

ri

 

 

 

Q

q

 

 

 

 

,

j 1, ..., 6,

j

FT

 

j

 

i

 

 

 

 

 

 

 

 

i 1

 

 

 

 

 

 

 

 

 

i

 

 

где Fi – вектор сил, действующих на i-й член, редуцированных

кточке начала i-й системы;

ri – вектор виртуального перемещения, являющегося нача-

лом i-й системы;

i – вектор виртуального вращения i-й части.

Зависимость между векторами ri , а также виртуального вращения i соответствующими обобщенными перемещениями, вы-

текающим из разложения скоростей точек А, В, С, которые определяют зависимость между соответствующими угловыми скоростями:

1 h1 ;

2 h1 .

458

По выполнению необходимых математических операций, вытекающих из уравнения Лагранжа II рода, движение анализируемой модели описывает систему дифференцируемых уравнений

(m1 m2 m3 m5 )xA (m1 m2 )l1( cos 2 sin )

(m3l m5l2 )( sin 2 cos ) 1 3;

(m1 m2 m3 m5 )yA (m1 m2 )l1( sin 2 cos )

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

(m1

 

 

 

2

m3l

2

2

 

m2 )(xA cos yA sin )l1 (m1 m2 )l1

 

m5l2

 

 

 

 

 

 

 

 

 

 

 

 

(m3l m5l2 )( xA sin yA cos )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

Ix1 Ix2 Ix3 Iz5 (Iz1 Iz2 )h1

(Iz1

Iz2 )h1 Ix3

 

 

Ms M0 s g n 3l sin 4l cos ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Iz3 N3 f3

3r cos( ) 4r sin( );

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(Iz1 Iz2 ) (Iz1 Iz2 )h1 MN N1 f1 N2 f2

 

1r cos 2r sin ;

(19.4)

Ix3 ( ) Ms M0 sgn ,

где m1, m2, m3, m5 – массы соответствующих частей; m3 – масса колеса 3, актуатора 6;

Ix1, Ix2 , Ix3, Iz1, Iz2 , Iz3, Iz5 – моменты инерции соответствующих частей осей;

Ix3 – учитывает момент инерции актуатора 6, определенным относительно оси za ;

N1, N2 , N3 – силы движения колес;

459