Алгоритмы управления движениями точки и робота-манипулятора. Учебное пособие
.pdfДелим одно уравнение на другое:
1 |
d r |
cos |
VM |
1 |
. |
|
r |
d |
sin |
||||
sin |
V |
|
||||
|
|
|
D |
|
|
Обозначим VM K , тогда
VD
d r |
cos |
|
1 |
|
|
r |
|
K |
|
|
d . |
sin |
|
||||
sin |
|
|
|||
Интегрируем последнее уравнение, используя начальные условия, и получаем
|
r r |
sin |
0 |
tgK 0 |
tg K |
1 |
. |
|||||
|
|
|
||||||||||
|
0 |
|
|
|
2 |
|
2 sin |
|||||
|
|
|
|
|
|
|
||||||
Назовем величину D0 sin 0 |
H относительной высотой. Тогда |
|||||||||||
|
|
|
|
|
0 |
|
cosK 1 |
|||||
|
r H tgK |
|
|
|
2 |
. |
|
|||||
|
|
|
|
|
|
|
||||||
|
|
|
|
2 |
|
2sin |
K 1 |
|||||
|
|
|
|
|
|
|
2 |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
||
Если K 1 , то |
r 0 . Значит, |
чтобы обеспечить попадание в цель, |
||||||||||
необходимо обеспечить K 1, |
т.е.VM VD . Попадание в цель проис- |
|||||||||||
ходит при развороте вектора скорости точки М на угол , т.е. точка М догоняет цель и попадает в хвост.
Пример 4. Жесткое касание
Пусть цель движется с постоянной скоростью VD 1238 м/c на по-
стоянной высоте H = 10 000 м. Начальная абцисса обнаружения цели xD = 5000 м. Точка М выпускается из начала системы координат с
постоянной по модулю скоростью VM 1500 м/c.
21
|
1 104 |
|
VD |
|
VD |
|
y |
|
|
|
|
D |
|
(t) 8 103 |
|
|
VМ |
|||
|
|
|
|
|||
|
6 103 |
|
|
|
||
|
|
|
|
|
|
|
|
4 103 |
|
|
|
|
|
|
2 103 |
|
M |
|
V |
|
|
0 |
|
М |
|
||
|
4 103 |
|
|
4 103 |
|
|
|
8 103 |
0 |
|
x(t) |
||
|
|
Рис. 1.12 |
|
|
||
На рис. 1.12 показана траектория точки М. Аналогично, как и в предыдущем методе, силу Р, обеспечивающую полученную зависимость скорости точки М, найдем из уравнения движения
m |
dVM |
P mg sin . |
(1.20) |
|
dt |
||||
|
|
|
1.3. МЕТОД НАВЕДЕНИЯ ПО ЛУЧУ
D |
Это самый распространенный метод |
наведения. Он состоит в том, что оператор |
|
|
О, точка М и цель D все время располага- |
ются на одной линии. Рассмотрим случай,
Окогда точка О неподвижна, а точки M и D
движутся. Получим кинематические урав-M
нения метода.
Уравнения движения
Поместим начало координат в точку О и считаем ее неподвижной. В этом методе вектор r | | rM , направлен вдоль rM . Угол, опреде-
ляющий положение радиуса-вектора |
rM , обозначим |
|||
наведения угол (t) считается известным. |
|
|||
|
|
|
|
|
Так как rM |
rM , то дифференцируя rМ , получим |
|||
|
drM |
drM r |
d . |
|
|
dt |
dt |
M dt |
|
. В процессе
(1.21)
22
Здесь
d , dt
|
|
|
|
|
|
|
|
|
|
|
где – угловая скорость луча локатора. |
|
|
||||||||
|
y |
|
|
|
|
|
|
VD |
|
D |
|
|
|
|
|
|
|
|
ц |
|
|
|
|
|
|
|
VM |
|
|
r |
rD |
D |
|
|
|
rM |
|
|
|
|
|
H |
|
|
|
|
|
|
|
|
|
|||
n |
|
|
M |
|
|
|
|
x |
||
|
|
|
|
|
|
|
||||
|
O |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 1.13 |
|
|
|||
Перепишем (1.21): |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
drM |
|
|
|
|
(1.22) |
|
|
VM |
dt |
|
rM n . |
|
||||
|
|
|
|
|
|
|
|
|
|
|
Проецируем (1.22) на касательную и нормаль: |
||||||||||
|
|
drM V |
cos( ), |
|
|
|||||
|
|
|
dt |
|
M |
|
|
|
|
(1.23) |
|
|
|
|
|
|
|
|
|
||
|
|
|
d |
|
|
|
|
|
||
|
|
r |
V |
|
sin( ). |
|
|
|||
|
|
M dt |
M |
|
|
|
|
|||
Обозначая через |
rD |
радиус-вектор цели, получим |
||||||||
|
|
|
r |
d V |
sin . |
|
(1.24) |
|||
|
|
|
D dt |
|
D |
|
|
|
||
Так как rD H / sin , то в трех уравнениях (1.23) и (1.24) будет три |
||||||||||
неизвесгных: rM (t), |
VM , |
|
(t) . |
|
|
|
|
|
|
|
23
Свяжем угол со временем. Для этого подставим выражение для rD в (1.24). Получим дифференциальное уравнение для :
|
|
d |
VD sin2 , |
|
|
|
|
(1.25) |
|||||||
|
|
dt |
H |
|
|
|
|
|
|
|
|
|
|
||
интегрируя которое, находим |
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
ctg 0 |
ctg |
VD |
t . |
|
|
|
|
(1.26) |
|||||
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
H |
|
|
|
|
|
|
||
Необходимо найти связи r ( ) |
и ( ) . |
Подставляя |
d |
во вто- |
|||||||||||
|
|||||||||||||||
|
|
|
|
M |
|
|
|
|
|
|
|
|
dt |
||
рое уравнение (1.23), получим |
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
sin( ) r |
VD |
sin2 |
|
. |
|
(1.27) |
||||||||
|
|
|
|
||||||||||||
|
|
|
|
M HV |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
M |
|
|
|
|
|
|
|||
Дифференцируем уравнение (1.27) по времени: |
|
|
|
||||||||||||
cos( ) d |
d |
drM |
|
VD |
|
sin2 r |
|
|
VD |
2sin cos d . |
|||||
|
|
|
|
||||||||||||
|
|
dt |
HVM |
|
|
|
M |
|
HVM |
|
|
dt |
|||
dt |
dt |
|
|
|
|
|
|
|
|
||||||
Разделим обе части последнего уравнения на ddt , с учетом (1.25) получим
d |
|
cos( ) 2sin( )ctg . |
cos( ) |
1 |
|
d |
|
|
Получили уравнение, связывающее и . Разрешим это уравнение относительно dd :
d |
2(1 ctg tg ( )) 2 sin cos( ) cos sin( ) |
, |
d |
sin cos( ) |
|
24
или
d |
|
1 sin cos cos sin sin |
|
1 sin2 |
(1 ctg ctg ) . |
||
d |
|||||||
|
2 |
sin |
|
2 |
|
||
Последнее уравнение линейно относительно ctg :
|
|
d ctg |
1 1 ctg ctg . |
|
|
|
||||||||
|
|
d |
2 |
2 |
|
|
|
|
|
|
|
|
||
Интегрируя уравнение (1.28), получим |
|
|
|
|
|
|
|
|
||||||
|
|
1 ctg d |
|
|
1 |
|
|
|
1 |
ctg d |
|
|
||
ctg e |
|
2 0 |
|
ctg 0 |
|
2 |
|
e2 0 |
|
d |
, |
|||
|
|
|
|
|
|
|
0 |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
sin 0 |
|
|
|
|
1 |
|
|
|
|
|
|
|
|
ctg |
ctg 0 |
|
|
|
|
|
sin d . |
|||||||
|
|
|
|
|
|
|||||||||
|
|
sin |
|
|
2 |
sin 0 |
|
0 |
|
|
|
|||
|
|
|
|
|
|
|
|
|
||||||
(1.28)
(1.29)
Уравнение (1.29) можно непосредственно решить численно: задавая , вычислить .
Таким образом, задавая 1, 2 ,... , по формуле (1.29) находим соответствующее значение . Далее по формуле
rM VM H sin( ) VD sin2
находим модуль радиуса-вектора, определяющего положение точки М, а по формуле
rD sinH
вычисляем модуль радиуса-вектора, определяющего положение цели.
25
Пример 5. Жесткое касание
Пусть цель движется с постоянной скоростью VD 1238 м/с на по-
стоянной высоте H = 10 000 м. Начальная абсцисса обнаружения цели xD = 5000 м. Точка М выпускается из начала системы координат с по-
стоянной скоростью VM 1500 м/c. На рис. 1.14 показана траектория
движения точки М. |
|
|
|
|
|
|
|
|
|
|
y(t) |
VD |
|
|
VD |
|
VD |
|
|
|
1 104 |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
5 103 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
5 103 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4 |
|
3 |
3 |
|
3 |
x(t) |
|
|
1 10 |
6 10 |
2 10 |
||||||
|
|
|
2 10 |
|
|
||||
|
|
|
|
|
|
Рис. 1.14 |
|
|
|
Сравнивая методы наведения по времени движения точки М до встречи с целью, отметим, что меньше времени надо по методу параллельного сближения – 7 с, при наведении по лучу – 9 с, по методу погони – 10 с.
Для вычисления силы Р, необходимой для полученной зависимости скорости точки М, используем уравнение движения (1.20).
26
2. УПРАВЛЕНИЕ ДВУХЗВЕННЫМ МАНИПУЛЯТОРОМ
СДВУМЯ СТЕПЕНЯМИ СВОБОДЫ
2.1.КИНЕМАТИКА
2.1.1. Постановка задачи кинематики
Робот-манипулятор состоит из двух звеньев (рис. 2.1) с известными размерами:
|
|
|
l1 OA, |
l2 AМ. |
|
|
Y |
|
|
|
|
VD |
|
|
|
|
|
|
|
|
|
|
|
|
|
VМ |
D |
|
|
|
|
|
y |
|
|
|
|
M |
|
|
|
|
|
|
|
|
|
|
|
|
|
C2 |
|
l2 y |
|
|
|
|
2 |
2 |
|
|
|
|
|
M A2 |
lC 2 y |
|
|
|
|
M A1 |
|
D0 |
||
|
|
2 |
|
|||
|
|
|
|
|
||
1 |
1 |
|
A |
|
l1y |
|
|
|
|
|
X |
||
O |
|
|
||||
|
|
|
|
|||
|
|
lC 2 x |
|
|
|
|
M 01 |
|
|
|
|
|
|
l1x |
l2 x |
|
x |
|
||
|
|
|
||||
|
|
|
|
Рис. 2.1 |
|
|
Управление манипулятором, осуществляемое двумя двигателями в шарнирах О и А, должно обеспечить за время τ сближение схвата М
27
с движущейся деталью D. Деталь движется по произвольной траектории, задаваемой координатным способом. Начальное неподвижное положение манипулятора задано углами поворота звеньев φ10, φ20, положение детали D0 – координатами xD(0), yD(0). К моменту времени t = τ требуется достичь совмещения точек М и D с относительной погрешностью δ, которая определяется [5] через рассогласование их перемещений:
x( ) y( ) .
x(0) y(0)
Для решения поставленной задачи необходимо:
1)составить уравнение движения точки М, а также уравнения угловых движений звеньев механизма;
2)выбрать параметр управления, обеспечивающий за время τ
сближение точек М и D с заданной погрешностью δ;
3) составить алгоритм численного интегрирования методом Эйлера уравнений движения на интервале времени τ.
2.1.2. Уравнения движений
Уравнения движения детали D
Деталь движется по произвольной траектории, задаваемой координатным способом:
xD (t) fDx (t), |
yD (t) fDy (t) . |
(2.1) |
Тогда скорости находятся через производные от координат:
V |
(t) |
dxD |
f |
(t), V |
(t) |
dyD |
f |
(t) . |
(2.2) |
|
|
||||||||
Dx |
|
dt |
VDx |
Dy |
|
dt |
VDy |
|
|
|
|
|
|
|
|
|
|
||
Например, если деталь движется по прямой с постоянной скоростью под углом α к оси x, то ее скорости и координаты определяются явно:
|
|
|
|
VD const, |
|
||
VDx VD cos , |
VDy VD sin , |
||
xD (t) xD (0) VDxt, |
yD (t) yD (0) VDyt. |
||
28
Уравнения движения схвата М
Введем радиусы-векторы точек М и D (рис. 2.2):
|
rM ixM jyM , |
|
|
|
rD ixD jyD |
|
||||
и вектор рассогласования координат |
r: |
|
|
|||||||
rD rM r , |
|
r rD rM , |
|
|||||||
r i (xD xM ) j ( yD yM ) i x j y. |
|
|||||||||
Здесь |
|
|
|
|
|
|
|
|
|
|
x x |
D |
x , |
y y |
D |
y |
M |
, |
r2 x2 y2 . |
(2.3) |
|
|
|
M |
|
|
|
|
|
|||
Скорость точки |
М представим в виде векторной суммы перенос- |
|||||||||
ной скорости, равной скорости точки D, |
и относительной |
V: |
||||||||
|
|
|
VM VM пер VM отн VD V , |
|
||||||
отсюда |
|
|
|
|
|
|
|
|
|
|
|
V VM VD , |
|
|
|
|
|
(2.4) |
|||
|
Vx VMx VDx , |
Vy VMy VDy . |
||||||||
|
|
|||||||||
|
|
|
|
|
|
VM |
|
|
||
|
|
|
|
VD |
|
|
|
|
VD |
|
|
|
|
|
M |
|
r |
|
D |
|
|
|
|
|
|
V |
|
|
|
|||
|
|
|
|
|
|
|
|
|
||
Y rM
rD
O 
X
Рис. 2.2
29
В качестве стратегии управления движением точки М примем два положения:
– во-первых, вектор относительной скорости должен быть направлен по вектору рассогласования r на точку D;
– во-вторых, по модулю вектор относительной скорости должен уменьшаться по мере приближения к точке D, чтобы в момент захвата обеспечить мягкое касание с нулевой относительной скоростью.
Для этого можно ввести линейную связь между относительной скоростью и рассогласованием, что равносильно линейной комбинации рассогласований и их производных:
|
|
k r , |
|
|
|
V |
|
(2.5) |
|||
Vx k x, |
Vy k y. |
||||
|
|||||
Здесь положительная константа k размерностью (1/с) является параметром управления и пока не определена. В более общем случае можно было бы ввести два параметра управления – отдельно по каждому направлению kx и ky, но тогда вектор относительной скорости будет направлен не по вектору рассогласования.
Из (2.4) и (2.5) можем определить скорости точки М:
VMx VDx k x |
VMx VDx k x, |
(2.6) |
|||
VMy VDy k y |
|
VMy VDy k y. |
|||
|
|||||
Отсюда видно, что чем больше параметр |
k, тем значительнее ско- |
||||
рость схвата М превышает скорость детали |
D на начальном участке |
||||
погони, когда рассогласования |
x и |
y велики. Параметр управле- |
|||
ния k определен далее, там же даны его количественные оценки в зависимости от требуемой погрешности δ.
Добавим к скоростям их связи с координатами и получим уравне-
ния движения точки М в форме уравнений Коши: |
|
|
||||||||
dxM V |
, |
V |
V |
k(x |
D |
x |
|
), |
||
dt |
Mx |
Mx |
Dx |
|
M |
|
||||
|
|
|
|
|
|
|
|
|
|
|
dyM V |
, |
V V |
k( y |
D |
y |
M |
). |
|||
dt |
My |
My |
Dy |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
30
