Научная статья на тему 'О выборе параметров механизма Беннетта и моделировании его в Matlab'

О выборе параметров механизма Беннетта и моделировании его в Matlab Текст научной статьи по специальности «Математика»

CC BY
397
64
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
КИНЕМАТИКА МЕХАНИЗМОВ / ПРОСТРАНСТВЕННЫЙ ЧЕТЫРЕХЗВЕННЫЙ МЕХАНИЗМ / МЕХАНИЗМ БЕННЕТТА / ЧИСЛО СТЕПЕНЕЙ СВОБОДЫ / МОБИЛЬНОСТЬ / ОБОБЩЕННЫЕ КООРДИНАТЫ / УРАВНЕНИЯ СВЯЗИ / МАТРИЦА НАПРАВЛЯЮЩИХ КОСИНУСОВ / MATLAB / KINEMATICS OF MECHANISMS / SPATIAL FOUR-LINK MECHANISM / BENNETT MECHANISM / NUMBER OF DEGREES OF FREEDOM / GENERALIZED COORDINATES / COUPLING EQUATIONS / DIRECTION COSINE MATRIX

Аннотация научной статьи по математике, автор научной работы — Иоффе Марк Львович

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

On the Choice of Parameters of the Bennett Mechanism and Modelling It in MATLAB

The article deals with the analysis of the so-called Bennett mechanism representing a four-link mechanism consisting of four links connected by flat hinges. From the point of view of mathematical apparatus and practical application, the approach proposed by the author is considerably simpler than those described in the literature. When such a mechanism is implemented, the choice of parameters that ensure its mobility plays a significant role as the number of degrees of freedom is negative. A method for the mechanism analysis and the choice of the mechanism parameters is proposed. The method comprises two approaches based on the generalized coordinates and the generalized velocities respectively. Using the proposed method, the author has obtained parameters of the mechanism that are essentially different from the ones proposed by Bennett. The kinematic model of the mechanism is developed in the MATLAB environment.

Текст научной работы на тему «О выборе параметров механизма Беннетта и моделировании его в Matlab»

УДК 531.134 DOI 10.18698/0536-1044-2016-2-11-19

О выборе параметров механизма Беннетта и моделировании его в MATLAB

М.Л. Иоффе

США, Нью-Йорк, 2 Авеню, д. 444

On the Choice of Parameters of the Bennett Mechanism and Modelling It in MATLAB

M.L. Ioffe

444, 2nd Avenue, New York, USA e-mail: ioffe.mark@gmail.com

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

Ключевые слова: кинематика механизмов, пространственный четырехзвенный механизм, механизм Беннетта, число степеней свободы, мобильность, обобщенные координаты, уравнения связи, матрица направляющих косинусов, MATLAB.

The article deals with the analysis of the so-called Bennett mechanism representing a four-link mechanism consisting of four links connected by flat hinges. From the point of view of mathematical apparatus and practical application, the approach proposed by the author is considerably simpler than those described in the literature. When such a mechanism is implemented, the choice of parameters that ensure its mobility plays a significant role as the number of degrees of freedom is negative. A method for the mechanism analysis and the choice of the mechanism parameters is proposed. The method comprises two approaches based on the generalized coordinates and the generalized velocities respectively. Using the proposed method, the author has obtained parameters of the mechanism that are essentially different from the ones proposed by Bennett. The kinematic model of the mechanism is developed in the MATLAB environment.

Keywords: kinematics of mechanisms, spatial four-link mechanism, Bennett mechanism, number of degrees of freedom, generalized coordinates, coupling equations, direction cosine matrix, MATLAB.

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

В соответствии с известной формулой число степеней свободы такого механизма равно 3-6 - 4-5 = -2. Это объясняется тем, что четыре обобщенные координаты — углы поворота в шарнирах — должны удовлетворять шести

Рис. 1. Пространственный четырехзвенный механизм

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

Впервые рассмотренный механизм, обладающий мобильностью [1], был предложен английским математиком Дж.Г. Беннеттом, который изложил его геометрическую теорию в работе [2]. В ней были рассмотрены три возможных случая существования в таком механизме мобильности, причем с одной степенью свободы. В первом случае все оси шарниров параллельны, во втором оси пересекаются в одной точке. Третий случай — наиболее интересный и открытый автором — относится к механизму, у которого на каждой из осей существует по точке, в которой соответствующая ось перпендикулярна отрезкам прямых, соединяющих эту точку с точками, лежащими на последующей и предыдущей осях. Указанные точки вместе с соединяющими их отрезками и образуют Skew Isogram — пространственный четырехугольник, у которого противоположные стороны равны. Теории механизма Беннетта посвящено значительное число работ [3-7]. В последнее время интерес к нему усилился вследствие возрастающей роли робототехники, где механизмы с параметрической мобильностью находят все большее распространение.

В известных работах анализируется лишь предложенный Беннеттом способ обеспечения

мобильности и соответствующий ему выбор параметров, т. е. указанный третий случай. В зарубежных работах для анализа и синтеза механизма применяется достаточно громоздкий математический аппарат, с помощью которого сложно получить практические результаты. Например, в работе [5] авторы используют дуальные кватернионы и подпространства 7-мерного проективного пространства. Кроме того, в научных трудах отсутствуют четко сформулированные условия, обеспечивающие его мобильность.

Проанализируем условия, которым должны удовлетворять параметры механизма для обеспечения его мобильности. Математически механизм определяется четырьмя векторами АВ, ВС, СБ, АБ, которые соединяют точки, принадлежащие подвижным телам, и расположены на осях вращения и четырьмя единичными векторами 11, 12,13,14, определяющими оси вращения. На рис. 1 показаны точки А, В, С, Б, а векторам 11,12,13,14 соответствуют оси аь аг, а3, а4. Все перечисленные векторы заданы своими проекциями на связанные с вращающимися телами оси, т. е. в процессе движения изменения значений углов поворота в шарнирах остаются постоянными. Полагаем, что при нулевых углах связанные с подвижными телами трехгранники совпадают друг с другом и с неподвижным трехгранником. Точки А и Б расположены на неподвижном теле. Таким образом, в качестве параметров механизма можно принять проекции указанных векторов, а в качестве обобщенных координат — углы поворота в шарнирах ф1,ф2,ф3,ф4. Условия, накладываемые на параметры механизма для обеспечения его мобильности — это условия, накладываемые на проекции указанных восьми векторов.

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

Рассматриваемый механизм (см. рис. 1) описывается следующей системой уравнений:

AD —

АВ +1ВС + CUTX2/x3 ) ГХ1 /Х2

X/х1 >

Í

dTx

X/Х1

di

d9i di

- Zz

dTx

X1/X2

d(?2

di di

dTX2 /X3 d<$3

di di

d TX3/X4 d<$4

di

Tx-i /x Л TX

di

X3/X4 TX2/X3 TX1 /x2TX/X1

11

i-y

0

- Z 2 Z 2

i-y 0

- T 3

z

Z 3

i-y 0

-Z4

z

Z 4

i-y

Zz1 0

-/1

X

Z2

z 0 -/X

Z 3

z 0 -/X

Z 4

z 0

-l 4

X

- Z

1

7

l 1

X

T:

X / X1 >

0

y

- T2 Л

iy

Z 2

X 0

/

- 73 Л

i-y

Z 3

X 0

- Z

T

X

X1/X2 >

(1)

X2/X3 >

У

4 Л

7

Z 4

X

T

X

X3/X4 '

Tx/X —

1 0 0

л

где Тх/хх, ТХ1/Х2, ТХ2/Хз, ТХз/Х4 — матрицы перехода соответствующих трехгранников.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

При постоянной угловой скорости дифференциальное уравнение

dT;

X/х

di

— m

0

-Zz

Zz Z

\

Zy Zx

TX/X —(OLT,

X / x

с начальным условием

dTx/х (0) di

— E

имеет решение

TX/х — E + L sin mi + L2(1 - cos mi).

(2)

Рассмотрим второй подход, позволяющий выбрать параметры для обеспечения мобильности.

Скорость точки Б согласно системе (1) определяется выражением

dAD — d91

di

f

х

х

0

- V

Z1

t-y 0

-Z2 Z 2

i-y

V 7 f

di L

Z1 -71 Л y

0 Z1

-Zx 0 y

Z 2 - 72 Л y

0 Z 2

-Zx 0 y

T;

X/X1

d92 / di '

х

^^ + CDTx2/x3 )x

(3)

TX1 / X2 TX/X1

i

+CD di

0

- Z 3 Z 3

i-y

Z 3 0 Z3

- Z

3

y

Z 3

X 0

TX2/X3TX1/X2 TX/x1 — 0.

Уравнение (3) представим в виде следующей однородной системы линейных уравнений:

£А5 = ** н 1 + ^ н 2 + ^ Н 3 = 0;

Л Л Л Л

H —

f

х

АВ + (ВС + CDTx2/X3 )

Z1 -Z1Л

í-z

-I1 0

ZV Zx 0

X1/ X2

х

0

ly

Z1

X

T;

X/X1 j

у

V ly

H2 — (С + CDTx2/X3 )x 0

(4)

х

2

- z Z 2

i-y

H — CD

В первом подходе при заданных параметрах механизма, обеспечивающих его мобильность, необходимо, задавшись значением угла фь определить углы ф2, ф3, ф4, удовлетворяющие первому и последнему уравнениям системы (1). Первое векторное уравнение является системой трех уравнений с тремя неизвестными, которое благодаря выбранным параметрам имеет решение. Определив из него углы фх, ф2, ф3, из последнего уравнения системы (1) найдем угол ф 4.

Z 2 z - T2 ^ y

0 Z 2 1-х TX1¡X2TX/х1 ;

-Zx 0

f 0 Z 3 z - T3 ^ y

- Zz3 0 Z 3 TX2 /X3

Z 3 V ly -Z 3 X 0 y

5 TX / XO TX

Эта система имеет ненулевые решения только в том случае, когда векторы Н1, Н2, Н3 линейно зависимы, т. е. справедливо равенство

(Н1 X Н 2) Н 3 = 0. Если равенство (5) выполняется, то

Н3 = Н1 + ^2 Н2 .

(5)

(6)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Параметры Л.1 и Л.2 определяются выражениями

А-1 =

Я 2 =

Н1Н 2)-( Н1Н1 Н1Н1)-(Н з Н1

Н3Н1)(Н2Н2)-(Н3Н2)(Н1Н2 Н1Н 2 Н 3 Н 2

Н1Н2)(Н1Н2)-(Н1Н1)(Н2Н2

Н 2 Н 2 Н1Н 2

(7)

Общее решение линейной системы имеет

вид

¿ф1 ¿ф2 ¿ф3 Л ' '

■) = ^ (, Я2,-1),

(8)

йг йг

где ^ — произвольная постоянная.

Таким образом, если выбрать в качестве независимой скорости отношение йф^йг, то две остальные скорости описываются соотношениями

йф2 Х2 йф1 йг йфъ

Х1 йг 1 Аф1

(9)

Н1 = (АВ + ВС + СБ

Н 2 =

(ВС

Н 3 = СБ

+ СБ

-/3 У3

»у

-/ 2 У 2

ч 1У

У 3 0

-/3

»"Х

0 /1 -11 ^ у

У1 0 у 1

Т1 у у -/Х 0) )

У 2 - У2 ^ у

0 У 2 >

-/Х 0 )

(10)

- У3 >

»•у

У 3 0

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

¿ф1 - т ¿ф2 т т ¿ф3 7 т ¿ф4 7 = п -11 т--12 т--13 т--14 = 0;

йг йг йг йг

Аф1 у X АБ+72 х ВБ+У3 X СБ = 0,

йг йг йг

(11)

где X — знак векторного произведения.

С помощью формул (11) можно показать, что когда оси шарниров пересекаются в одной точке, механизм обладает мобильностью. В этом случае второе векторное уравнение (11) следует из первого и справедливы следующие выражения:

АБ = Х1 /1 - х 4 /4; ВБ = Х2 /2 - Х4 /4; СБ = Х3 /3 - Х4 /4.

(12)

йг Х1 йг

Из _системы_уравнений (4) следует, что векторы Н1, Н2, Н± зависят от параметров системы векторов АВ, ВС, СБ, /1, /2, /3 и углов поворота в шарнирах ф1, ф2, ф3. Как уже указывалось, в начальный момент времени углы ф1,ф2,ф3 равняются нулю, а оси связанных трехгранников параллельны и одинаково направлены. В соответствии с системой уравнений (4) имеем выражения для векторов:

Подставив выражения (12) во второе уравнение системы (11), получим

^ У X аБ+^ у X ВБ + у X СБ = йг йг йг

(йф1 7 ,#27 , ¿ф3 7 ^ X (л %\

= -Х4 I —— /1 +—— /2 — /3 IX /4. (13) ^ йг йг йг )

В соответствии с первым уравнением системы (11)

йф1 7 йф^ 7 йф 3 7 _ йф 4 7 • /1 т--/2 т--/3 =--/4 .

йг

йг

йг

йг

(14)

Из выражений (13) и (14) непосредственно следует второе уравнение (11).

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

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

Два равенства системы (11) в совокупности составляют однородную систему шести линей-

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

Предположим, что три вектора

m1 = AD, m2 = BD, m3 = CD

(15)

неколлинеарны. Полагаем:

11 = xm + yim3;

12 = x2mi + у2тъ;

13 = x3mi + узт2.

(16)

dt

-x1m1 x m2 - y1m1 x m3j +

(17)

йф2 , _ _ _ _ \

+-(Х2т1 X т2 - у2т2 Xт3) +

йг

йф3 , _ _ _ _ \ _

+-(Х3т1 Xт3 + у3т2 Xт3) = 0.

йг

Необходимым и достаточным условием существования ненулевых производных углов, удовлетворяющих выражению (17), является линейная зависимость векторов — коэффициентов при этих производных, т. е.

(- Х1т1 X т2 - у1т1 X т3) X (Х2т1 X т2 - у2т2 X т3) *

(18)

* (x3m1 x m3 + y3m2 x m3 j =

= (у1Х2у3 - Х1 у2Х3 )[( Xт2 )* т3 ] ,

где * — знак скалярного произведения.

Рассмотрим частный случай. Выберем параметры Х1, у1, Х2, у2, Х3, у3, удовлетворяющие следующим выражениям:

Х2 = Х1;

Х3 = у1; (19)

у3 = у2.

Эти значения удовлетворяют уравнению (18).

Подставив выражения (19) во второе уравнение (17), получим:

¿ф1 _ йф 2- , йф 3 _

-п1 +-п2 +-п3 = 0;

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

йг йг йг

п1 = -Х1т1 хm2 -у1т1 Xт3; п2 = Х1т1 X т2 - у2т2 X т3; (20)

п3 = у1т1 X т3 + у2т2 X т3; п1 + п2 + п3 = 0.

Таким образом, второе уравнение (17) справедливо при произвольной величине йф1/йг и при значениях двух других производных углов, равных

¿ф 2 йф1

dt dt dф 3 dф1

(21)

Подставив выражения (16) и (15) в систему (11), получим

dф1 / _ _ \ _ dф2 / _ _ \

-(x1m2 + y1m3 jx m1 +-(x2m1 + y2m3 jx

dt dt

_ dфз , _ _ \ _ _

x m2 +-(x3m1 + y3m2 jx m3 = 0;

dt

dt dt

Свободные параметры^, y1, y2 определим из условия, что векторы l1, l2, i3 являются единичными. В соответствии с выражениями (16) получим следующую систему уравнений:

xm + y2m32 + 2x1 y1m2 * m3 = 1; x2m12 + y2m32 + 2x1 y2m1 * m3 = 1; (22) y12m12 + y22m22 + 2 y1 y2m1 * m2 = 1.

Простейшее решение системы (22) получим, приняв следующие соотношения:

x1 = y1 = y2 = x;

|m11 = |m2| = |m3| = m; (23)

m1 * m2 = m1 * m3 = m2 * m3 = m2 cos 9.

Система (22) в этом случае будет справедлива при произвольных величинах m и cos 0 и при значении x, равном

1 1

x =-■ =-. (24)

my] 2(1 + cos 9) 2m cos 9/2

В соответствии с выражениями (16) получим

11 = x (m2 + m3 j;

12 = x (m1 + m3 j; (25)

13 = x (m1 + m2 j.

Векторы m1, m2, m3 удовлетворяют уравнениям (23), а параметр x — уравнению (24). При этих условиях значение dф1/dt может быть произвольным, а значения двух других произ-

водных углов равняются этому значению. Значения йф]^/^ и 14 определим из первого уравнения (17):

d91/_ ч d94 -

2х-(m1 + m2 + m3) —--Z4.

di ' di

Тогда

Z4 —

1

m1 + m2 + m3 j;

y¡3myj 1 + 2cos 6

^ — 2x7^/172^9 ^ di di

(26)

Рассмотрим пример механизма, параметры которого определяются формулами (25), (26) и следующими выражениями:

m1 — m (0,0,1);

m2 — m (sin 6, 0, cos б);

m3 — m (sin 6 cos P, sin 6 sin p, cos б); (27)

„ cos 6-cos2 6 cos P —-rT6-.

sin2 6

В этом случае имеем:

АВ — m1 - m2;

вс — m2 - m3;

CD — m3; AD — m1; Z1 — x (m2 + m3); Z2 — x (m1 + m3); Z3 — x (mm + m2);

1

x —-

2w cos 9/2 Очевидно, что

H1 + H 2 + H 3 = 0.

Тогда

Ai = ^2 = -1

и в соответствии с выражениями (9) получим:

dt dt (28) dф3 _ dф1

dt dt

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

Для решения этого вопроса используем модель механизма, созданную с помощью программного обеспечения МЛТЬЛВ. В начальный момент времени положение механизма определяется формулами (27). Таким образом, механизм задается двумя параметрами: т и 0. По

ним_находят векторы т1, т2, т3, 11, 12, 13,

АВ, ВС, СБ. Начальные значения углов ф равняются нулю, а соответствующие матрицы перехода являются единичными. Полученные векторы определены своими проекциями относительно связанных со звеньями трехгранников и не меняются при движении механизма. Параметры механизма выбраны таким образом, что выполняется_ условие мобильности, т. е. векторы Н1, Н 2, Н 3 линейно зависимы, система имеет одну степень свободы и мы можем задать угловую скорость первого звена произвольно. Две другие скорости будут при этом равняться первой. Определив параметры механизма, начальные условия и начальные угловые скорости, можно найти положение механизма через некоторый промежуток времени и проверить условие подвижности для этого нового положения. В случае его выполнения расчет можно продолжить, приняв новое положение в качестве начального. Таким образом, расчет проводится реккуррентно.

Реккурентный расчет на каждом шаге заключается в следующем. Определяются матрицы Тх/Х!, Тх1/Х2, Тх2/Х3 по следующим формулам:

(

L1 —

L —

L3 —

0

- Tz1

Z1

t-y

- Z 2

z

Z1

v y

f 0

- V

Z 3

y

Zz1 0

- Zx

Z 2

0

- Z1 Z 3

lz

0

- T 3

«-X

- Z1

y

Л

- Z

2

y

Z 2

X

0

- Z

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

j

3 Л

y

Z 3

X

Txц = \E + 1\ +

йг

+1?

1 | dфl 1 - С081 -1

йг

Тх/^(0);

Tx1/x2 = ] E + L2sin| йф2г | +

+ L22

1 , йф2 . 1 - ^ | — г

йг

(29)

>Тх1/х2 (0);

Т

1 X

х2/х3 = ^ Е + Lз sin| йфзг | +

Л

+ 4

1 I йфз .

1 - ^ | — г йл

ТХ2/Х3 (0).

_Вычисляются проекции векторов АВ, ВС,

СБ на оси неподвижного трехгранника ОХХХ:

АВ' = АВТх/Х1;

ВС = ВСТХ1/Х2 ТХ/Х1; (30)

СБ = СБТХ2 / Х3 ТХ1/Х2 ТХ/ Х1 .

Определяются векторы Н1, Н 2, Н 3:

Н1 =

L{T:

X/х1 ;

Проверяется условие коллинеарности векторов Н1, Н2, Н3, т. е. вычисляется тройное векторно-скалярное произведение

р = (Н1 X Н2) * Н 3. (32)

Если р = 0, то условие коллинеарности выполняется, т. е. механизм обладает подвижностью. В этом случае по формулам (7) и (8) вычисляются новые значения угловых скоростей. Параметры Л.1 и Л.2 можно определить также по формуле

(ь^2 ) =

\(Н1 *Н1) ((*Н2)) 1 \(Н1 *Н3]

(Н1 * Н2 ) ((2 * Н2

Н 2 * Н 3

.(33)

АВ + (ВС + СБТх2/х3 )ТХ1/х2

Н2 =(ВС + СБТх2/х3 )ТХ1/Х2ТХ/х1 ; (31)

Н 3 = СБ-£<3ТХ2/Х3ТХ1/Х2 ТХ/Х1.

Рассмотренная реккурентная модель была реализована в функции Ьеппей.ш программы МАТЬЛВ для механизма с параметрами ш = 1, 9 = л/4. Критичным для работы модели оказался выбор значение приращения угла на шаге вычисления, т. е. значение (йф:/йг)г. Проверка условия подвижности на каждом шаге состояла в проверке справедливости неравенства abs (р) > 1Е-8. При таком значении порога максимальная величина приращения угла составила 0,000003. При большем значении условие мобильности благодаря погрешностям расчета не выполнялось уже после первого шага. Процесс вычисления оказался достаточно медленным и для получения конечного значения пер-

0,50

0,45

0,40

0,35

0,30

& 0,25 к

Ц 0,20

0,15 0,10 0,05

0

*

* е

+

"Г"

-5 " у ф

* + ^ ж*

*ф * <ь

г

10 15 Шаг

а

20

25

0,50

0,45

0,40

0,35

0,30

& 0,25 н

Ц 0,20

0,15 0,10 0,05

0

*

-э *

+

+

ф г

+ * ^ |>ф

ж®

^ ф >ф

г

Ф

10 15 Шаг

б

20

25

Рис. 2. Углы поворота ф1 (+), ф2 (*), ф3 (О) механизма Беннетта, полученные при реккурентном (я)

и нереккурентном (б) расчетах

вого угла, равного 0,3, понадобилось 100 000 шагов.

Результаты расчета трех углов фь ф2, ф3 механизма Беннетта приведены на рис. 2, а. По оси абсцисс отложены прореженные номера шагов, т. е. из общего числа 100 000 взят равномерно 21 шаг. По оси ординат отложены значения углов поворота (рад) в первых трех шарнирах механизма.

Рассмотренный реккурентный алгоритм при известных параметрах механизма не является единственным методом вычисления углов. Для этой же цели можно использовать конечные формулы, т. е. выражения, в которых угол поворота относительно первой оси является известным, а углы поворота относительно второй и третьей осей определяются решением соответствующей нелинейной системы уравнений:

АВ = АВТх/х1; ВС = ВСТх1 / х2 Тх/х1;

СБ = СБТх2 /х3 Тх1 /х2 Тх/х1 ;

АВ + ВС + СБ' - т1 = 0; (34)

Тх / х1 = Е +118ш(ф1) + Ь21 [1 - С08(ф1)];

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Тх1/х2 = Е +12 8т(ф2) + 4 [1 - С08(ф2)];

Тх2/х3 = Е + 13 8ш(ф3) + Ь23 [1 - С0в(ф3 )].

Литература

В среде МЛТЬЛВ нелинейные системы решаются с помощью функции ^попНп.ш, которая использовалась в разработанной функции ЬеппеШ.ш для реализации описанных в формуле (34) вычислений углов поворота механизма с теми же параметрами, что и при рекку-рентном расчете. Результаты расчета углов приведены на рис. 2, б. По оси абсцисс отложены номера каждого из 21 шага приращения первого угла. Приращение угла принято равным 0,3/20. По оси ординат отложены значения углов поворота (рад) в первых трех шарнирах механизма.

Выводы

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

2. Разработана модель кинематики механизма в среде МЛТЬЛВ.

[1] Bennett G.T. A new mechanism. Engineering, 1903, vol. 76, pp. 777-778.

[2] Bennett G.T. The scew isogram-mechanism. Proceedings of the London Mathematical Society,

1913-1914, vol. 13, pp. 151-173.

[3] Hunt K.H. Kinematic Geometry of Mechanisms. Oxford engineering science series, University

Press, Oxford. 1978. 484 p.

[4] Perez A., McCarthy J.M. Dimensional synthesis of Bennett linkages. Proceedings of the ASME

Design Engineering Technical Conferences, 2000, Baltimore, Maryland, USA, paper DETC2000/MECH-14069.

[5] Brunnthaler K., Schrocker H.-P., Husty M. A New Method for the Synthesis of Bennett

Mechanisms, Proceedings of CK2005, International Workshop on Computational Kinematics, Cassino, May 4-6, 2005, pp. 1-8.

[6] Яруллин М.Г., Мингазов М.Р. Синтез структурных модификаций механизма Беннет-

та. URL: http://www.mmf.spbstu.ru/mese/2014/129.pdf (дата обращения 2 сентября 2015).

[7] Яруллин М.Г., Мингазов М.Р. Краткий анализ модификации механизмов Беннетта.

Проблемы механики современных машин. Мат. У междунар. конф., Улан-Удэ, Изд-во ВСГУТУ, 2012, т. 1, с. 177-181.

[8] Яруллин М.Г., Мингазов М.Р. Структурный синтез двухподвижного пространственно-

го 5R механизма и элементы следящего управления. Известия Самарского научного центра Российской академии наук, 2014, № 6, с. 214-220.

References

[1] Bennett G.T. A new mechanism. Engineering, 1903, vol. 76, pp. 777-778.

[2] Bennett G.T. The scew isogram-mechanism. Proceedings of the London Mathematical Society,

1913-1914, vol. 13, pp. 151-173.

[3] Hunt K.H. Kinematic Geometry of Mechanisms. Oxford engineering science series, University

Press, Oxford. 1978. 484 p.

[4] Perez A., McCarthy J.M. Dimensional synthesis of Bennett linkages. Proceedings of the ASME

Design Engineering Technical Conferences, 2000, Baltimore, Maryland, USA, paper DETC2000/MECH-14069.

[5] Brunnthaler K., Schrocker H.-P., Husty M. A New Method for the Synthesis of Bennett

Mechanisms, Proceedings of CK2005, International Workshop on Computational Kinematics, Cassino, May 4-6, 2005, pp. 1-8.

[6] Iarullin M.G., Mingazov M.R. Sintez strukturnykh modifikatsii mekhanizma Bennetta [Syn-

thesis of structural modifications mechanism Bennett]. Available at: http:// www.mmf.spbstu.ru/mese/2014/129.pdf (accessed 2 September 2015).

[7] Iarullin M.G., Mingazov M.R. Kratkii analiz modifikatsii mekhanizmov Bennetta [A brief

analysis of the modification arrangements Bennett]. Problemy mekhaniki sovremennykh mashin. Materialy 5 mezhdunarodnoi konferentsii [Problems of modern machines. Proceedings of the 5 International Conference]. Ulan-Ude, VSGUTU publ., 2012, vol. 1, pp. 177-181.

[8] Iarullin M.G., Mingazov M.R. Strukturnyi sintez dvukhpodvizhnogo prostranstvennogo 5R

mekhanizma i elementy slediashchego upravleniia [Structure Synthesis of Differential 5R Linkage and Tracking Control Elements]. Izvestiia Samarskogo nauchnogo tsentra Rossiis-koi akademii nauk [Proceedings of the Samara Scientific Center of the Russian Academy of Sciences]. 2014, no. 6, pp. 214-220.

Информация об авторе

ИОФФЕ Марк Львович (Нью-Йорк) — кандидат технических наук. США, Нью-Йорк, 2 Авеню, д. 444, e-mail: ioffe.mark@gmail.com).

Статья поступила в редакцию 14.12.2015 Information about the author

IOFFE Mark Lvovich (New York) — Candidate of Science (Eng.) (444, 2nd Avenue, New York, USA, e-mail: ioffe.mark@gmail.com).

i Надоели баннеры? Вы всегда можете отключить рекламу.