T3=A1∙A2∙A3; T4=T3∙A4. |
(4.166) |
Каждый из векторов множества { Rio3 } характерных точек стрелы вида
Rio3 = [xio3 yio3 zio3 1]T ; io3 [1; co3], |
(4.167) |
где xio3, yio3, zio3 – координаты точки io3 стрелы в собственной локальной декартовой системе координат (№ 3), переносится в инерциальную систему координат:
Rio3,0 =T3 × Rio3 ; io3 [1; co3], |
(4.168) |
здесь Rio3,0 – вектор положения точки io3 в инерциальной системе координат вида
R |
= [x |
y |
z |
1]T . |
(4.169) |
io3,0 |
io3,0 |
io3,0 |
io3,0 |
|
|
Для каждой точки io3 [1; co3] выполняется проверка условия превышения ее вертикальной координаты над соответствующей вертикальной координатой полидистантной поверхности препятствий с теми же координатами xio3, zio3:
yio3≥YЭ(xio3, zio3). (4.170)
Аналогичная проверка выполняется для множества { Rio4 } харак-
терных точек телескопического звена стрелы с использованием матрицы T4:
Rio4,0 |
=T3 × Rio4 ; |
io4 [1; co4]; |
(4.171) |
||
R |
= [x |
y |
z |
1]T ; |
(4.172) |
io4,0 |
io4,0 |
io4,0 io4,0 |
|
|
|
|
yio4≥YЭ(xio4, zio4). |
|
(4.173) |
||
В случае выполнения условия (4.170) для всех точек стрелы множества { Rio3 }, где io3 [1; co3], а также выполнения условия (4.173) для всех точек телескопического звена множества { Rio4 }, где io4 [1;
co4] значение отдельной компоненты с индексом iq8 вектора Cross индикатора пересечения звеньев ГПК с препятствиями принимается равным 0, что соответствует отсутствию пересечений. В противном случае значение отдельной компоненты вектора Cross индикатора пересечения звеньев ГПК с препятствиями принимается равным 1, что соответствует пересечению с препятствиями:
241
ìCross(iq8)=0 при
ïí("(io3Î [1,co3]); yio3 ³YЭ(xio3, zio3))Ù("(io4Î [1,co4]); yio4 ³YЭ(xio4, zio4));(4.174)
ïîCross(iq8)=1 в остальных случаях.
Учитывается специфика задачи, которая выражается в допущении, что подъем стрелового оборудования (увеличение значений координат q8, q9, q10 в пределах их возможных диапазонов при сохранении постоянных декартовых координат груза) в случае отсутствия пересечений с препятствиями не может приводить к пересечению с препятствиями. Тогда условие нахождения минимальных диапазонных значений управляемых координат qn8Н, qn9Н, qn10Н в начальной точке положения груза с учетом непересечения подвижных звеньев ГПК с препятствиями будет на отдельной итерации цикла iq8 [2; iq8max] выглядеть следующим образом:
qn8Н=q8(iq8); qn9Н=q9(iq8); qn10Н=q10(iq8) |
(4.175) |
при Cross(iq8)=0 и Cross(iq8–1)=1. |
|
По зависимостям, аналогичным (4.163)–(4.175), в цикле iq8 [1; iq8max] определяются уточненные минимальные диапазонные значения управляемых координат в конечной точке положения груза qk8Н, qk9Н, qk10Н. Максимальное значение индекса iq8max в этом случае определится следующим образом:
iq8max = ë(qk8B − qk8N ) u8 û. |
(4.176) |
Условие нахождения минимальных диапазонных значений управляемых координат qn8Н, qn9Н, qn10Н в конечной точке положения груза с учетом непересечения подвижных звеньев ГПК с препятствиями будет на отдельной итерации цикла iq8 [2; iq8max] выглядеть аналогично:
qk8Н=q8(iq8); qk9Н=q9(iq8); qk10Н=q10(iq8) |
(4.177) |
при Cross(iq8)=0 и Cross(iq8–1)=1. |
|
6. Определяется диапазон [qnд7; qkд7] допустимых значений угла поворота поворотной колонки q7 для генерации вершин графа дорожной карты, который в общем случае должен превышать диапазон из начального и конечного значений данной координаты [qn7; qk7], чтобы обеспечить обход возможных препятствий, однако будет меньше угла в 360°.
242
Для этого определяется приращение диапазона координаты q7:
dq7=δq7–(|qk7–qn7|/2), |
(4.178) |
где δq7 – заданная постоянная величина.
Диапазон допустимых значений угла поворота поворотной колонки будет равен
[qnд7 ; qkд7 |
]= íì[(qn7 − dq7 );(qk7 + dq7 )] |
при qn7 ≤ qk 7 ; |
(4.179) |
|||||||||||||
|
î |
[(q |
n7 |
+ dq |
7 |
);(q |
k 7 |
- dq |
7 |
)] |
при q |
n7 |
> q |
k7 |
. |
|
|
|
|
|
|
|
|
|
|
|
|||||||
7. Используя вложенные циклы по индексам iqn8 [1; iqn8max] и iqk8 [1; iqk8max], варьируются значения управляемой координаты q8 в начальной и конечной точках положения груза qn8 и qk8 соответственно в пределах уточненных диапазонных значений управляемых координат вида (4.162) с шагом u8. При этом значения управляемых координат q9 и q10, определяемые в начальной и конечной точках положения груза по методике раздела 4.3, также будут варьироваться в пределах собственных уточненных диапазонных значений.
Для каждого сочетания значений [qn8; qn9; qn10] и [qk8; qk9; qk10] выполняется поиск оптимальной по значению целевой функции траек-
тории перемещения груза.
7.1. Окончательное значение целевой функции L* (T*, Ae* либо C*) как общий результат работы алгоритма принимается равным бесконечно большой величине:
Ae*=∞; T*=∞; C*=∞. |
(4.180) |
7.2. Определяются максимальные значения индексов iqn8 и iqk8 по уточненным диапазонным значениям:
iqn8max = ë(qn8B − qn8N ) u8 û; iqk8max = ë(qk8B − qk8N ) u8 û. |
(4.181) |
7.3. На каждой итерации индексов iqn8 и iqk8 по их значениям определяются текущие значения координаты q8 в начальной и конечной точках (qn8 и qk8):
qn8=qn8N+(iqn8–1)∙ u8; qk8=qk8N+(iqk8–1)∙ u8. |
(4.182) |
7.4. Определяются текущие значения координат qn9, qk9 и qn10, qk10, соответствующие текущим значениям координат qn8 и qk8 по методике раздела 4.3 [103, 107, 109, 112, 113, 115, 116].
243
7.5. Генерируется случайным образом множество вершин Sr={s2,…, sng–1} графа дорожной карты, представляющих собой точки в пространстве конфигураций ГПК, т.е. возможные положения груза в
пределах диапазонов координат [qnд7; qkд7] [q8min; q8max] [q9min; q9max]; [q10min; q10max], в которых он не пересекается с препятствиями.
Для создания дорожной карты при помощи генератора случайных чисел создается ng точек в пространстве конфигураций с координатами
sp=(q7p, q8p, q9p, q10p), p [2; ng–1], |
(4.183) |
где
q7p=qnд7+Rand∙(qkд7–qnд7); q8p=q8min+Rand∙(q8max–q8min); (4.184)
q9p=q9min+Rand∙(q9max–q9min); q10p=q10min+Rand∙(q10max–q10min).
Значения q7p, q8p, q9p, q10p, полученные для каждого значения индекса p, должны удовлетворять условиям непересечения подвижных звеньев ГПК и груза с эквидистантной (полидистантной) поверхностью [YЭ], что выражается проверкой по условиям (4.170) и (4.173) для подвижных звеньев ГПК, а также приведенным ниже аналогичным зависимостям для груза.
Формируется матрица перехода A5 из локальной системы координат груза в локальную систему координат телескопического звена:
где |
|
|
|
|
A5=Ax5 Aν5 Ay5Aω5, |
|
|
|
|
|
|
(4.185) |
|||||
é cos(− q8 ) |
sin(− q8 ) |
0 |
0ù |
|
|
|
é1 |
0 |
0 |
0 |
ù |
||||||
|
|
|
|
||||||||||||||
A |
ê- sin(- q |
|
) |
cos(- q |
|
) |
0 |
0ú |
; A |
= |
ê0 |
1 |
0 |
(-q )ú |
|||
= ê |
8 |
|
|
8 |
|
|
ú |
ê |
|
|
10 |
ú; |
|||||
v5 |
ê |
0 |
|
|
0 |
|
|
1 |
ú |
|
y5 |
|
ê |
|
1 0 |
ú |
|
|
ê |
|
|
|
|
0ú |
|
|
|
ê0 0 |
ú |
||||||
|
ë |
0 |
|
|
0 |
|
|
0 |
1û |
|
|
|
ë0 0 |
0 1 |
û |
||
|
|
|
|
é–cosq11 |
|
0 |
−sinq11 |
0ù |
|
|
|
|
|
||||
|
|
|
|
ê |
|
|
|
|
|
|
|
ú |
|
|
|
|
|
|
|
A |
= ê |
0 |
|
1 |
|
0 |
|
0ú |
, |
|
|
|
(4.186) |
||
|
|
w5 |
|
ê |
sinq11 |
|
0 |
|
cosq11 |
|
ú |
|
|
|
|
|
|
|
|
|
|
ê |
|
|
0 |
|
|
|
|
|
|||||
|
|
|
|
0 |
|
0 |
|
0 |
|
|
ú |
|
|
|
|
|
|
|
|
|
|
ë |
|
|
|
1û |
|
|
|
|
|
||||
здесь q11 – поворот груза вокруг вертикальной оси грузового каната относительно стрелы крана.
Формируется матрица перехода из локальной системы координат груза в инерциальную систему координат:
244