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

2426

.pdf
Скачиваний:
7
Добавлен:
07.01.2021
Размер:
8.54 Mб
Скачать

Вывод коэффициентов полинома четвертой степени для системы уравнений (4.15)–(4.18).

В уравнении тора (4.18), согласно выражению (4.16), в принятой системе координат x2+z2=ρ2. В то же время, согласно (4.30), ρ2 может быть выражено через y. Подставив (4.30) в (4.18), получим уравнение с одним неизвестным:

(h12 + (y × к + y0 )2 + y2 + R2 - r2 )2 -

 

- 4 × R2 × (h12 + (y × к + y0 )2 )= 0.

 

(4.74)

Дальнейшие преобразования приводятся без пояснений:

(h 2 + y2

× к2 + 2 × y × к × y

0

+ y 2 + y2 + R2

- r2 )2 -

1

(h 2 + y2 × к2

 

0

 

+ y 2 )= 0;

 

- 4 × R2 ×

+ 2 × y × к × y

0

 

 

1

 

 

 

0

 

 

(y2 ×(к2 +1)+y×(2×к×y )+(h2

+y 2 +R2

-r2))2 -

 

 

 

0

1

 

0

 

 

 

 

-4×R2 ×h

2 -4×R2 ×y2 ×k2

-4×R2 ×y×2×k×y -4×R2 ×y 2

=0;

1

 

 

 

 

0

0

 

y4 ×(к2 +1)2 + y3 ×(2×(к2 +1)× 2×к × y0 )+

+ y2 ×(2×(h 2 + y 2 + R2 - r 2 ))+ y2 ×(2×к × y )2 +

+y ×(2× 2×к × y0 ×(h12 + y02 + R2 - r 2 ))+

+(h12 + y02 + R2 - r 2 )2 - 4× R2 × h12 - 4× R2 × y2 ×к2 -

- 4× R2 × y × 2× к × y0 - 4 × R2 × y02 = 0.1 0 0

Если обозначить

(к2 +1)2 = s1;

(2 × (к2 +1)× 2 × к × y0 )= s2 ;

2× (h12 + y02 + R2 - r2 )+ 22 × к2 × y02 - 4× R2 × к2 = s3;

4×к × y0 ×(h12 + y02 + R2 - r2 )- 4× R2 ×2×к × y0 = s4 ;

(h12 + y02 + R2 - r2 )2 - 4× R2 ×h12 - 4× R2 × y02 = s5 ,

то может быть получена классическая форма полинома четвертой степени:

y4·s1+y3·s2+y2·s3+y·s4+s5=0,

(4.75)

200

 

уравнение которого может быть решено численными методами. Моделирование описанной методики с использованием ПК под-

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

Упрощенная методика определения управляемых координат ГПК по известным координатам груза при нулевых углах наклона шасси.

В случае обеспечения горизонтального положения опорной платформы при помощи систем автоматического горизонтирования

[58, 79, 84, 98, 118, 125, 126, 156, 158, 162] может быть использована упрощенная методика определения управляемых координат ГПК по известным координатам груза, которая заключается в описанной ниже последовательности шагов [103]. Предполагаются известными декартовы координаты точки подвеса груза [xгр; yгр; zгр], заданные в неподвижной системе координат O0Х0Y0Z0, а также координаты точки O1 начала системы координат O1Х1Y1Z1, связанной с базовым шасси ГПК:

[q1; q2; q3; 0; 0; q6].

1. Определяются координаты груза в полярной системе координат, связанной с базовым шасси (значение координаты угла поворота q7 и вылета ρ). Используются схема, приведенная на рис. 4.18, а, и зависимости для перевода декартовых координат в полярные на плоскости [31]:

ρ =

 

.

 

(xгр - q1 )2 + (zгр - q3 )2

(4.76)

ìarctg((zгр - q3 ) (xгр - q1)) ïïarctg((zгр - q3 ) (xгр - q1))+π

q6,7 = ïíarctg((zгр - q3 )(xгр - q1))-π

ïïπ2 ï-π 2

î

при (xгр - q1)> 0;

при (xгр - q1)< 0; (zгр - q3 )³ 0;

при (xгр - q1)< 0; (zгр - q3 )< 0; (4.77)

при (xгр - q1)= 0; (zгр - q3 )> 0; при (xгр - q1)= 0; (zгр - q3 )< 0.

Отсюда при известном значении q6

 

q7= q6,7 q6.

(4.78)

2. Определяются максимальные и минимальные значения диапазонов управляемых координат ГПК [q8В q8Н]; [q9В q9Н].

201

 

q1

 

xгр

X0

O0

 

 

X2, X3, X4

zгр

 

 

 

 

4

 

0, 5, 6, 7

3

 

 

Z4

 

 

 

 

 

 

q6,7=q6+q7

 

а)

 

q7

 

 

 

 

q3

1

X1

 

q6

 

 

X0ʹ

 

Z1

2

Z2

 

 

 

Z3

 

Z0

Y4

 

q9

X3, X4

4

y4,43

 

 

 

 

 

 

3

 

5

 

 

 

 

Y1,

Y2

x3,33

 

R2,5

 

 

q10 y5,7

 

 

 

 

 

 

 

 

 

 

б)

Y3

 

 

 

 

q8

 

 

 

 

 

 

 

 

q8,1

 

 

 

6

 

 

 

 

 

 

O

q8,2

 

 

 

 

y6,7

 

X2

 

 

 

 

 

 

 

 

 

 

2

3

 

 

 

 

 

 

 

 

 

 

y1,2

O1,O2

Y0

X2

y0,6=yгр

7

 

 

y0,7=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

q2

 

 

q7

 

0

=q2+y1,2

 

 

Z1

 

X1

 

 

 

 

 

 

 

 

Z0 O0

 

 

 

 

 

 

 

 

 

x1,2

X0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ

 

 

 

 

 

Рис. 4.18. Схема для определения управляемых координат ГПК при нулевых углах наклона шасси

2.1. Для этого определяются значения вылета ρ, соответствующие четырем возможным сочетаниям максимальных и минимальных допустимых конструкцией ГПК значений координат q8 и q9 (точки 1, 2, 3, 4 положения оголовка стрелы):

202

Точка 1 (ρ=ρ1):

q8=q8max; q9=q9min;

(4.79)

Точка 2 (ρ=ρ2):

q8=q8max; q9=q9max;

(4.80)

Точка 3 (ρ=ρ3):

q8=q8min; q9=q9min;

(4.81)

Точка 4 (ρ=ρ4):

q8=q8min; q9=q9max.

(4.82)

Для каждого сочетания зна- Y1 2 чений q8 и q9 используется после-

довательность вычислений:

 

 

 

 

q8,1=arctg(y4,43/(x3,33+q9));

 

1

 

 

R2,5 = (x3,33 + q9 )2 + y4,432 ;

 

 

 

 

 

 

q8,2=q8q8,1;

(4.83)

 

 

 

4

ρ=R2,5∙cos q8,2 – x1,2.

(4.84)

 

 

 

По (4.83)–(4.84) определяют-

 

 

3

ρ

 

 

ся последовательно 4 постоянных

ρ1

ρ2

ρ3

ρ4

значения: ρ1, ρ2, ρ3, ρ4.

 

 

Рис. 4.19. К определению диапазонов

2.2. По заданному текущему

значению

ρ

определяются

гра-

 

 

вылета ρ

 

ничные

значения

диапазонов

 

 

 

 

 

 

 

 

управляемых координат [q8В

q8Н];

При ρ1ρ<ρ2

 

[q9В q9Н].

 

 

 

 

 

 

 

 

 

 

q8В=q8max;

q8Н=arctg(y4,43/(x3,33+q9min))+arccos((x1,2+ρ)/ y4,432 +(x3,33 + q9min )2 ); (4.85)

q9В=(x1,2+ρ–y4,43∙sin q8max)/cos q8max x3,33;

q9Н=q9min.

При ρ2ρ<ρ3

q8В=arctg(y4,43/(x3,33+q9max))+arccos((x1,2+ρ)/ y4,432 + (x3,33 + q9max )2 );

q8Н=arctg(y4,43/(x3,33+q9min))+arccos((x1,2+ρ)/

 

 

);

 

 

y4,432 +(x3,33 +q9min )2

(4.86)

q9В=q9max; q9Н=q9min.

 

 

При ρ3ρρ4

 

 

q8В=arctg(y4,43/(x3,33+q9max))+arccos((x1,2+ρ)/

y4,432 + (x3,33 + q9max )2

);

q8Н=q8min; q9В=q9max;

(4.87)

q9Н=(x1,2+ρ–y4,43∙sin q8min)/cos q8min x3,33.

 

 

В случае, если значение вылета ρ находится вне интервала конструктивно допустимых значений

203

ρ [ρ1; ρ4 ],

(4.88)

делается вывод о невозможности обеспечить требуемые координаты груза при текущем положении базового шасси ГПК, алгоритм завершает свою работу.

2.3. По текущему значению ρ и граничным значениям координат [q8В q8Н]; [q9В q9Н] определяются граничные значения диапазона управляемой координаты [q10В q10Н].

ρ [ρ14 ], используя (4.83), получим

y5,7=tg(q8–arctg(y4,43/(x3,33+q9)))∙(ρ+x1,2).

(4.89)

Согласно схеме на рис. 4.18, б,

 

q10=y5,7 y6,7 = y5,7 – (yгрy0,7)= y5,7 – (yгр– (q2+y1,2))=

 

= tg(q8–arctg(y4,43/(x3,33+q9)))∙(ρ+x1,2) – (yгр– (q2+y1,2)),

(4.90)

где y4,43 x3,33 x1,2 y1,2 – постоянные конструктивные значения.

 

Соответственно значения [q10В q10Н] определятся как

 

q10Н = tg(q8Н –arctg(y4,43/(x3,33+q9Н)))∙(ρ+x1,2) – (yгр– (q2+y1,2));

 

q10В = tg(q8В –arctg(y4,43/(x3,33+q9В)))∙(ρ+x1,2) – (yгр– (q2+y1,2)).

(4.91)

2.4. В случае, если полученное по (4.91) значение q10В меньше минимальной конструктивно возможной длины грузовой лебедки от оголовка стрелы

q10В<q10min,

(4.92)

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

Если выполняется условие

q10Вq10min q10Н<q10min,

(4.93)

где – знак логического умножения (конъюнкции), т.е. точка подвеса груза расположена внутри пространства возможных положений оголовка стрелы (позиция 4 на рис. 4.7), q10Н корректируется:

q10Н=q10min. (4.94)

После этого также последовательно корректируются значения q9Н и q8Н по (4.97) и (4.96) с подстановкой значения q10Н.

204

Пуск

1

3

Задание исходных данных: координат точки подвеса груза xгр, yгр, zгр,

координат шасси q1q6, постоянных конструктивных параметров кра-

на x1,2, y1,2, x3,33, y4,43, q8max, q8min, q9max, q9min, q10max, q10min

 

 

 

 

 

 

 

 

5

 

ρ =

(xгр q1 )2 + (zгр q3 )2

 

 

 

 

 

 

 

 

 

7

 

 

 

 

 

 

 

 

 

Определение q6,7 по (4.77)

 

9

 

 

 

 

 

 

 

 

 

 

 

 

q7= q6,7 q6

 

 

 

 

 

 

 

 

 

 

 

 

11

 

 

 

 

 

 

 

q8=q8max; q9=q9min

13

 

 

 

 

 

 

 

 

 

 

 

 

 

q8,1=arctg(y4,43/(x3,33+q9));

 

 

 

 

R2,5 =

(x3,33 + q9 )2 + y4,432

;

 

 

 

q8,2=q8q8,1; ρ1=R2,5∙cos q8,2 – x1,2

q8=q8max; q9=q9max 15 17

q8,1=arctg(y4,43/(x3,33+q9));

R2,5 = (x3,33 + q9 )2 + y4,432 ; q8,2=q8q8,1; ρ2=R2,5∙cos q8,2 – x1,2

19

q8=q8min; q9=q9max 21

q8,1=arctg(y4,43/(x3,33+q9));

R2,5 = (x3,33 + q9 )2 + y4,432 ; q8,2=q8q8,1; ρ3=R2,5∙cos q8,2 – x1,2

23

q8=q8min; q9=q9min 25

q8,1=arctg(y4,43/(x3,33+q9));

R2,5 = (x3,33 + q9 )2 + y4,432 ; q8,2=q8q8,1; ρ4=R2,5∙cos q8,2 – x1,2

Нет

59

Да

Задается q9

 

 

 

 

 

 

27

 

 

 

q8=q8max; q9=q9min

 

 

 

Нет

 

29

 

 

 

 

 

ρ1ρ<ρ2

 

 

31

 

 

 

 

Да

 

 

 

Определение [q8В

q8Н]; [q9В

q9Н]по (4.85)

 

Нет

 

33

 

 

 

 

 

 

ρ2ρ<ρ3

 

 

35

 

 

 

 

Да

 

 

 

Определение [q8В

q8Н]; [q9В

q9Н]по (4.86)

 

Нет

 

37

 

 

 

 

 

ρ3ρρ4

 

 

39

 

 

 

 

Да

 

 

 

Определение [q8В

q8Н]; [q9В

q9Н]по (4.87)

 

Вывод сообщения о невоз-

41

43

 

 

 

Останов

можности обеспечить тре-

 

буемые координаты груза

45 q10Н=tg(q8Н–arctg(y4,43/(x3,33+q9Н)))∙(ρ+x1,2)–(yгр

– (q2+y1,2)); q10В=tg(q8В–arctg(y4,43/(x3,33+q9В)))∙(ρ+x1,2)–(yгр

– (q2+y1,2))

Да

47

 

Нет

 

 

 

 

49

 

 

51

q10В<q10min

 

 

Да

(q10Вq10min)

 

 

 

 

 

 

 

 

 

 

 

 

 

(q10Н<q10min)

 

 

 

 

Коррекция q10Н, q9Н и q8Н

 

 

 

 

 

 

 

 

 

 

 

Нет

 

 

 

 

по (4.94), (4.97) и (4.96)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

53

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Нет

Задается q8

55

 

 

 

 

 

 

 

Да

 

 

 

 

q9=(x1,2+ρ–y4,43∙sin q8)/cos q8 x3,33

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

57

 

 

 

 

 

 

Определение q10 по (4.90)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q8=arctg(y4,43/(x3,33+q9))+arccos((x1,2+ρ)/

71

 

 

 

 

/

y4,432 + (x3,33 + q9 )2

)

 

 

61

Определение q10 по (4.90)

 

q9 = y4,432 + ( (q10 + yгр q2 y1,2 )2 + (x1,2 + ρ)2 )2 x3,33

63

q8=arctg(y4,43/(x3,33+q9))+arccos((x1,2+ρ)/ 65

/ y4,432 + (x3,33 + q9 )2 )

67

Вывод q8В, q8Н, q9В, q9Н,

q10В, q10Н, q7, q8, q9, q10

Останов 69

Рис. 4.20. Блок-схема алгоритма определения управляемых координат ГПК по упрощенной методике при нулевых углах наклона шасси

205

3. По значению одной из управляемых координат (q8, q9, q10), заданному внутри соответствующего диапазона [q8В q8Н], [q9В q9Н], [q10В q10Н], определяются значения двух оставшихся управляемых координат.

 

 

 

Задана координата q8.

 

 

 

q9=(x1,2+ρ–y4,43∙sin q8)/cos q8 x3,33;

(4.95)

q10 определяется по (4.90).

 

 

 

 

 

 

 

 

Задана координата q9.

 

 

 

 

 

 

 

 

);

 

q8=arctg(y4,43/(x3,33+q9))+arccos((x1,2+ρ)/ y4,43

2 + (x3,33 + q9 )2

(4.96)

q10 определяется по (4.90).

 

 

 

 

 

 

 

 

Задана координата q10.

 

 

 

 

 

 

x3,33 ;

 

 

 

2 + (

 

 

 

)2

 

q9 = y4,43

(q10 + yгр q2 y1,2 )2 + (x1,2 + ρ)2

(4.97)

q8 определяется по (4.96).

Блок-схема алгоритма определения управляемых координат ГПК при нулевых углах наклона шасси приведена на рис. 4.20.

Рассмотренный алгоритм позволяет моделировать управляемые координаты ГПК для заданного положения груза, при нулевых углах наклона шасси; позволяет в автоматизированном режиме решать обратную задачу кинематики ГПК [103]. Алгоритм может использоваться при планировании траекторий перемещения груза, при решении задач кинематического анализа и синтеза конструктивных параметров ГПК и технологических параметров его рабочего процесса.

4.4.Методика проверки положения автомобильного крана

впространстве конфигураций по ограничению на устойчивость

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

206

использована как составляющая часть общей методики поиска оптимальной траектории в пространстве конфигураций ГПК.

Математическая модель динамической системы ГПК [104, 105], сведенная в систему из 11-ти дифференциальных уравнений второго порядка с постоянными коэффициентами с 11 неизвестными функциями – обобщенными координатами, позволяет решать задачи статики и динамики ГПК.

Статическая составляющая разработанной модели использовалась для оценки устойчивости ГПК и проверки положения ГПК в пространстве конфигураций по ограничению на устойчивость [86, 87, 93, 96, 111, 117, 119, 123, 127, 128].

Практическую значимость представляют при этом зависимости обобщенных координат qi, i [1;6], которые использовались для определения значений показателей устойчивости ГПК. Полученные обобщенные координаты также использовались для моделирования значений информационных

параметров для вычисления крите-

Рис. 4.21. Прямоугольный опорный рия устойчивости [93, 96]. контур ГПК

Для проверки устойчивости положения ГПК, имеющего прямоугольный опорный контур (рис. 4.21), в пространстве его конфигура-

ций предлагается использовать следующий критерий оценки устойчивости:

ξ = min

ìk ;

k

 

;

1

;

1

ü

,

(4.98)

2

 

 

ý

 

í

1

 

 

k1

k2

 

 

 

î

 

 

 

 

þ

 

 

где k1, k2, 1/k1, 1/k2 – показатели устойчивости для четырех осей опрокидывания, входящих в опорный контур, вычисляемые на основе нормальных реакций в опорных элементах:

k =

R1

+ R2

;

1

= R3 + R4 ;

k

2

=

R2 + R3

;

1

=

R1 + R4

,

(4.99)

 

 

 

 

 

 

1

R3

+ R4

k1

R1 + R2

 

 

R1

+ R4

 

k2

 

R2 + R3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

207

 

 

 

 

 

 

 

где R1, R2, R3, R4 – нормальные реакции на опорных элементах ГПК

(см. рис. 4.21).

Из этих показателей только два (k1, k2) являются независимыми. Аналитические выражения для вычисления критерия устойчиво-

сти (sgn – функция знака числа):

ξ1=k2∙sgn(1+ sgn(k1k2))+k1∙sgn(1+ sgn(k2k1));

 

ξ2=(1/k2)∙sgn(1+ sgn(1/k1–1/k2))+ (1/k1)∙sgn(1+ sgn(1/k2–1/k1)).

(4.100)

Тогда

 

ξ=ξ2∙sgn(1+ sgn(ξ1ξ2))+ξ1∙sgn(1+ sgn(ξ2ξ1)).

(4.101)

После подстановки (4.100) в (4.101) полное аналитическое выражение критерия будет иметь вид:

ξ=((1/k2)∙sgn(1+sgn(1/k1–1/k2))+(1/k1)∙sgn(1+sgn(1/k2– –1/k1)))∙sgn(1+sgn((k2∙sgn(1+sgn(k1k2))+k1∙sgn(1+sgn(k2k1)))–

– ((1/k2)∙sgn(1+sgn(1/k1–1/k2))+(1/k1)∙sgn(1+sgn(1/k2– –1/k1)))))+(k2∙sgn(1+sgn(k1k2))+k1∙sgn(1+sgn(k2

k1)))∙sgn(1+sgn(((1/k2)∙sgn(1+sgn(1/k1–1/k2))+(1/k1)∙sgn(1+sgn(1/k2– –1/k1)))–(k2∙sgn(1+sgn(k1k2))+k1∙sgn(1+sgn(k2k1))))). (4.102)

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

Для получения нормальных реакций R1, R2, R3, R4 на опорных элементах ГПК как первичных информационных параметров использовались отклонения обобщенных координат q2, q4, q5, q6, полученные в результате решения задачи статики при заданных больших значениях обобщенных координат q1q11.

Вертикальное смещение δyio,0 характерной точки выносной гидравлической опоры io [1;4] ГПК вдоль оси O0Y0 инерциальной системы координат является второй компонентой вектора Rio,0 , задающе-

го смещения характерной точки опоры вдоль координатных осей инерциальной системы координат:

v

= [δxio,0

δyio,0 δzio,0 1]T ; i [1;4],

(4.103)

Rio,0

Вектор Rio,0 определяется по формуле

 

 

 

R

= R

 

-T * × R

,

(4.104)

 

io,0

io,1

1 io,1

 

 

 

 

 

208

 

 

где Rio,1 – вектор положения характерной точки опоры io в локальной системе координат базового шасси, к которому относится масса каж-

дой выдвижной опоры, R

= [x

y

z

1]T ; T1* – матрица перехо-

io,1

io,1

io,1

io,1

 

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

ординат (qj+qj), j=2;4;5;6.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пуск

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ввод исходных данных:

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q1q11; m1m5; { Ris }

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Расчет статических отклонений обобщенных

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

координат q1q11 по модели ГПК

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

é cos(q6 + q6 )×cos(q5 + q5 )

 

sin(q5 + q5 )

 

 

 

-sin(q6 + q6 )×cos(q5 +q5 )

0

ù

 

 

 

 

 

 

 

 

 

ê

sin(q +q )×sin(q +q

)-

 

 

cos(q5 + q5 )´

 

 

cos(q + q

)×sin(q + q

)+

 

 

 

 

 

ú

 

 

 

 

 

 

 

 

 

 

 

6

 

6

 

4

 

4

 

 

 

 

 

 

 

 

6

 

 

 

 

6

 

4

4

 

 

 

 

 

(q + q )ú

 

 

 

 

 

*

ê-cos(q + q

)×cos(q + q )´

 

 

 

+sin(q + q )×cos(q +q )´

 

 

 

 

 

T1 = ê´sin(q + q )

 

 

 

4

 

4

 

 

´cos(q4 +q4 )

 

 

´sin(q + q

)

 

 

 

4

 

4

 

 

 

2

 

ú .

 

 

 

 

 

 

 

 

ê

 

6

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

6

 

 

 

 

 

 

 

 

 

2 ú

 

 

 

 

 

 

 

 

 

 

 

5

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ê

sin(q + q )×cos(q + q

 

)+

 

-cos(q5 +q5 )´

 

cos(q

 

+ q

)×cos(q

+ q

 

)-

 

 

 

 

ú

 

 

 

 

 

 

 

 

 

ê

 

6

 

6

 

4

4

 

 

 

 

 

 

6

 

 

 

 

6

 

4

4

 

 

 

 

 

 

ú

 

 

 

 

 

 

 

 

 

ê+cos(q + q

)×sin(q +q )´

 

 

-sin(q +q

)×sin(q +q )´

0

ú

 

 

 

 

 

 

 

 

 

ê

 

6

6

 

 

 

4

 

4

 

 

´sin(q4 + q4 )

 

 

 

 

 

 

 

6

 

 

6

 

 

 

4

 

4

 

 

 

 

 

ú

 

 

 

 

 

 

 

 

 

 

 

 

+ q5 )

 

 

 

 

 

 

 

 

 

 

 

 

 

´sin(q5 +q5 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ê´sin(q5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ú

 

 

 

 

 

 

 

 

 

ê

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

1

 

ú

 

 

 

 

 

 

 

 

 

ë

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

û

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

R1=co∙δy1o,0; R2=co∙δy2o,0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

= R

 

 

 

 

- T * × R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R3=co∙δy3o,0; R4=co∙δy4o,0

 

 

 

 

 

 

io ,0

io ,1

io ,1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

11

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

15

 

 

 

 

 

 

 

 

 

k1 =

R1 + R2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Да

 

 

 

 

 

 

 

Нет

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R2 + R3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξ1<ξ2

 

 

 

 

 

R3 + R4

 

 

 

 

 

 

 

 

 

 

k2 =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R1 + R4

 

 

 

 

 

 

 

 

 

 

 

 

 

16

 

 

 

 

 

 

17

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Да

 

k1<k2

 

Нет

 

 

 

 

 

 

 

 

 

1

1

12

 

 

 

 

 

 

 

 

 

 

ξ=ξ1

 

 

 

 

 

 

 

ξ=ξ2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Да

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Нет

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10

 

 

 

 

 

 

 

9

 

 

 

 

 

 

 

 

 

k

<

k

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

18

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξ1=k1

 

 

 

 

 

ξ1=k2

 

 

 

 

 

 

 

13

 

1

 

 

 

 

 

 

14

 

 

 

 

 

 

 

Да

 

 

 

 

 

 

 

 

 

Нет

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξ < ξкрит

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξ2=1/k1

 

 

 

 

 

 

 

ξ2=1/k2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

19

 

20

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

21

 

Opr =1

 

 

Opr =0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

22

 

 

 

 

 

Вывод

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Останов

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

результатов: Opr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 4.22. Блок-схема алгоритма определения критерия устойчивости ГПК на основе нормальных реакций

209

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]