Материал: 2426

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

[qn7]; [qn8В; qn8Н]; [qn9В; qn9Н]; [qn10В; qn10Н] – для начальной точки;

[qk7]; [qk8В; qk8Н]; [qk9В; qk9Н]; [qk10В; qk10Н] – для конечной точки. (4.162)

4. Построение полидистантной поверхности вокруг реальной по-

верхности препятствий по исходным данным ([YПР]; lзап_г и lзап_в) по методике, изложенной в разделе 2.3 [102, 94]. В результате построе-

ния формируется дискретная матрица высот [YЭ] того же размера, что

иисходная матрица препятствий [YПР].

5.Определение уточненных минимальных диапазонных значений

управляемых координат qn8Н, qn9Н, qn10Н, qk8Н, qk9Н, qk10Н с учетом непересечения подвижных звеньев ГПК с препятствиями.

Для начальной точки положения груза в цикле с шагом u8 значение управляемой координаты q8 меняется от qn8Н до qn8В. Изменению координаты q8 ставится в соответствие индекс iq8 [1; iq8max], где

iq8max = ë(qn8B qn8N ) u8 û.

(4.163)

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

Формируются матрицы переходов из локальных систем координат (базового шасси, поворотной колонки, стрелы, телескопического звена A1, A2, A3, A4 соответственно) в локальные системы координат более высокого иерархического уровня (находящиеся ближе к инерциальной системе координат). A1=T1 имеет вид (4.157), A2, A3, A4 определяются по зависимостям, аналогичным (4.156), с подстановкой значений управляемых координат q7, q8 и q9:

где

écosq7

0

-sinq7

 

 

ê

0

1

0

A

=ê

w2

ê sinq7

0

cosq7

 

ê

0

0

0

 

ë

A2=Ax2Ay2Az2Aγ2Aν2 Aω2;

A3=Ax3Ay3Az3Aγ3Aν3 Aω3; (4.164) A4=Ax4Ay4Az4Aγ4Aν4 Aω4,

0ù

 

écosq8

sinq8

0

0ù

 

é1 0 0 q9ù

ú

 

ê

 

 

 

ú

 

ê

ú

0ú; A

-sinq8

cosq8

0

0ú; A =ê0 1 0

0ú. (4.165)

ú

v3

ê

0

0

1

ú

x4

ê

ú

0ú

 

ê

0ú

 

ê0 0 1

0ú

1û

 

ë

0

0

0

1û

 

ë0 0 0 1û

Формируются матрицы перехода из локальных систем координат стрелы и телескопического звена в инерциальную систему координат (T3 и T4 соответственно):

240

T3=A1A2A3; T4=T3A4.

(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:

yio3YЭ(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

 

 

 

yio4YЭ(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. Определяется диапазон [q7; q7] допустимых значений угла поворота поворотной колонки q7 для генерации вершин графа дорожной карты, который в общем случае должен превышать диапазон из начального и конечного значений данной координаты [qn7; qk7], чтобы обеспечить обход возможных препятствий, однако будет меньше угла в 360°.

242

Для этого определяется приращение диапазона координаты q7:

dq7=δq7–(|qk7qn7|/2),

(4.178)

где δq7 – заданная постоянная величина.

Диапазон допустимых значений угла поворота поворотной колонки будет равен

[q7 ; q7

]= íì[(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} графа дорожной карты, представляющих собой точки в пространстве конфигураций ГПК, т.е. возможные положения груза в

пределах диапазонов координат [q7; q7] [q8min; q8max] [q9min; q9max]; [q10min; q10max], в которых он не пересекается с препятствиями.

Для создания дорожной карты при помощи генератора случайных чисел создается ng точек в пространстве конфигураций с координатами

sp=(q7p, q8p, q9p, q10p), p [2; ng–1],

(4.183)

где

q7p=q7+Rand∙(q7q7); q8p=q8min+Rand∙(q8maxq8min); (4.184)

q9p=q9min+Rand∙(q9maxq9min); q10p=q10min+Rand∙(q10maxq10min).

Значения 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