CC BY
183
130
Общетехнические задачи и пути их решения
УДК 621.01 Ф. А. Доронин
ПРИМЕНЕНИЕ ТЕОРИИ ВИНТОВ В СТАТИКЕ И ДИНАМИКЕ ПРОСТРАНСТВЕННЫХ МЕХАНИЗМОВ ПАРАЛЛЕЛЬНОЙ СТРУКТУРЫ
Дата поступления: 12.05.2015 Решение о публикации: 09.07.2015
Цель: Продемонстрировать возможности применения винтового исчисления при статических и динамических расчетах пространственных механизмов, содержащих замкнутые кинематические цепи и часто называемых. Отмечена значительная роль отечественных ученых в создании, развитии и применении винтового исчисления в различных задачах теории машин и механизмов. Методы: Решение некоторых задач механики механизмов параллельной структуры основано на использовании принципа возможных скоростей и свойств силовых и кинематических винтов, характеризующих нагруженное состояние механизма. Применение математической среды Mathcad позволяет значительно облегчить численные расчеты. Результаты: Продемонстрированы способы определения положения мгновенной оси вращения выходного звена (платформы) механизмов параллельной структуры. Показан простой способ построения матрицы жесткости системы, позволяющей определить не только матрицу инерции, но и основные характеристики малых свободных колебаний платформы механизмов параллельной структуры (собственные частоты и формы колебаний) с учетом упругости опорных звеньев. Практическая значимость: Продемонстрированы способы определения положения мгновенной оси вращения выходного звена механизма при наличии различных комбинаций связей, наложенных на систему. Предложенная методика позволяет существенно уменьшить трудоемкость расчетов, связанных с определением основных параметров, характеризующих свободные колебания механизмов параллельной структуры.
Винтовое исчисление, силовой винт, кинематический винт, механизм параллельной структуры, мгновенная винтовая ось, относительный момент, плюккеровы координаты.
Feliks A. Doronin, Cand. Sci. (Eng.), associate professor, [email protected] (Petersburg State Transport University) APPLICATION OF THE THEORY OF SCREWS IN THE STATICS AND DYNAMICS OF SPATIAL MECHANISMS OF PARALLEL STRUCTURE
Objective: To demonstrate possibilities of using screw calculus in static and dynamic calculations for spatial mechanisms that contain closed kinematic chains and are frequently designated. The significant role of Russia’s researchers in creating, development and application of screw calculus to different tasks of theory of machines and mechanisms is noted. Methods: Solution of some tasks of mechanics of parallel structure mechanisms is based on deployment of the principle of virtual velocities and qualities of screws of forces and kinematic screws that characterise a mechanism’s loaded condition. Using Mathcad mathematical environment allows for significant easing of numerical calculations. Results: Methods for determination of position of instantaneous axis of rotation of output component (platform) of parallel structure mechanisms are shown. A simple method for construction of matrix of system rigidity which allows to establish not only inertia matrix but also primary characteristics of small-size autonomous oscillations of parallel structure mechanisms’ structure (proper frequencies and oscillation modes) that takes into account elasticity of support links. Practical importance: The study demonstrated methods of determination of position of instantaneous axis of rotation of a mechanism’s output component when different combinations of connections applied to the system are provided. The method proposed allows for significant decrease in complexity of calculations related to determining primary parameters that characterise autonomous oscillations of parallel structure mechanisms.
Screw calculus, screws of forces, kinematic screw, mechanisms of parallel structure, instantaneous screw axis, relative moments, Pluckers coordinates.
2015/3
Proceedings of Petersburg Transport University
Общетехнические задачи и пути их решения
131
Ирландский астроном сэр Роберт Стоелл Болл (1840-1913 гг.) в 1867 г. стал профессором прикладной математики в Королевском Колледже науки в Дублине. Там он читал лекции и писал работы по механике, в которых излагал теорию винтов. За трактат «Теория винтов» [7] Болл в 1879 г. был награжден Королевской ирландской Академией медалью Каннингема, вручаемой раз в три года за выдающиеся результаты в науке. В этом трактате он, опираясь на работы предшественников (в том числе французского математика и механика Луи Пуансо (1777-1859 гг.), французского геометра Мишеля Шаля (1793-1880 гг.), немецких математиков Августа Фердинанда Мёбиуса (1790-1868 гг.), Юлиуса Плюкке-ра (1801-1868 гг.) и Феликса Клейна (18491925 гг.)), ввел в механику твердого тела понятие о силовом винте и описал основные операции с силовыми винтами.
Экстраординарный профессор кафедры механики физико-математического факультета Новороссийского университета (с 1905 г. -демократически избранный ректор этого университета) Иван Михайлович Занчевский (1861-1928 гг.) опубликовал работу [4], в которой изложил теорию векторов, линейных комплексов, винтов и их групп, а также приложения теории винтов к механике.
В 1895 г. появился замечательный труд русского математика и механика Александра Петровича Котельникова (1865-1944 гг.) [6], в котором было последовательно изложено винтовое исчисление и впервые сформулирован «принцип перенесения» для векторной алгебры. Принцип перенесения позволяет изучать геометрию одного пространства, используя геометрию другого, более изученного (например, в кинематике этот принцип позволяет установить соответствие между сферическим и свободным движениями твердого тела). Котельникову удалось истолковать все формулы теории кватернионов как «неразвернутые» формулы теории бикватернионов, т. е. установить полную аналогию тех и других формул.
Профессор Боннского университета Эдуард Штуди (1862-1930 гг.) в работе [13] ука-
зывал на связь параболического бикватерниона с группой движений эвклидова пространства, а в 1901 г. в фундаментальном труде [12] сформулировал принцип перенесения, постулирующий аналогию между векторными и винтовыми операциями.
Идеи теории винтов были развиты учеником Штуди немецким математиком и механиком Ричардом Эдлером фон Мизесом (18831953 гг.). В 1924 г. вышли в свет его работы [9, 10], в которых изложено несколько иное направление теории винтов - моторное исчисление. Достоинством моторного исчисления является то, что оно приспособлено к решению задач динамики, правда, за счет утраты полной аналогии с векторным исчислением, которая характерна для винтового исчисления. Принцип перенесения в этом случае не работает. Мизес рассмотрел ряд приложений, в том числе в механике стержневых систем.
На основе созданного в 1895 г. А.П. Котельниковым винтового исчисления советский математик Д. Н. Зейлигер развил теорию линейчатых поверхностей и конгруэнций [5].
Винтовое исчисление в теории механизмов применил Ф. М. Диментберг. В трудах [1, 2] помимо подробного обзора литературы обстоятельно изложены основы теории винтов и моторного исчисления и приведено множество приложений этих теорий к решению различных задач анализа и синтеза пространственных механизмов.
В 1956 г. появилась статья В. Е. Гью (В. Е. Гауфа) [8], а в 1965 г. - статья Д. Стюарта [11], в которых были описаны принципиально похожие устройства, представляющие собой неподвижное основание, к которому шестью стержневыми опорами прикреплена подвижная платформа. При изменении длины опорных стержней платформа изменяет свое положение в пространстве. Впоследствии эти устройства в научной литературе стали называться платформами Гауфа - Стюарта (гексаподами), а также механизмами параллельной структуры (параллельными механизмами). Механизмы параллельной структуры благодаря своим свойствам (относительно малому
ISSN 1815-588Х. Известия ПГУПС
2015/3
132
Общетехнические задачи и пути их решения
весу, повышенной грузоподъемности, жесткости и высокой точности) применяются в различных областях современного машиностроения. Одним из важных преимуществ механизмов параллельной структуры является то, что они в отличие от традиционных манипуляторов содержат замкнутые кинематические цепи, поэтому воспринимают нагрузку подобно пространственным фермам.
Впоследствии появилось значительное число работ, в которых было показано, что во многих случаях винтовое исчисление оказывается весьма эффективным при исследовании статики, кинематики и динамики механизмов параллельной структуры.
Применение теории групп винтов к определению усилий в абсолютно твердых опорных стержнях платформы
Однородная прямоугольная горизонтальная платформа весом G = 2500 Н, находящаяся в покое под действием заданных сил P = 1000 Н и Q = 600 Н, прикреплена к неподвижному основанию шестью стержневыми опорами 1-6 (рис. 1). Поставим задачу определить усилия Sl-S6 в опорных стержнях 1-6 платформы, если известны размеры a = 2 м, b = 3 м, с = 4 м.
Применим принцип возможных скоростей. Согласно принципу освобождаемости от связей, мысленно удалим опорный стержень 1 (рис. 2). Платформа получит одну степень свободы, совершая в общем случае
Рис. 2. Схема платформы с удаленным стержнем
винтовое движение, описываемое пока неизвестным кинематическим винтом V. Сохранившиеся стержневые опоры накладывают на винт V следующее условие: скорости шарниров O, A и B должны лежать в плоскостях, перпендикулярных осям опорных стержней. Это означает, что единичные винты нулевого параметра e2-e6, оси которых совпадают с осями оставшихся стержней 2-6, после приведения к центру О являются взаимными с винтом V1 и образуют линейный комплекс, определяемый этим винтом.
Будем считать, что все опорные стержни растянуты и единичные векторы e2-e6 направлены от шарниров O, A и B, соответственно. Запишем плюккеровы координаты единичных винтов e2-e6, приведенных к началу О выбранной системы координат Oxyz (см. рис. 2):
e2 =
= (- cos у 0 - sin у b sin у 0 -b cos y);
e3 _ (0 0 -1 b о о
e4 = (0 - cos P - sin P о о 0);
e5 = (0 0 -1 0 о о
e6 = (cosу 0 - sin у 0 -c sin Y 0)
Перемещения шарниров O, A и B происходят в полярных плоскостях этих стержней. Иными словами, относительные моменты винтов, соответствующих стержням 2-6, с винтом V равны нулю.
2015/3
Proceedings of Petersburg Transport University
Общетехнические задачи и пути их решения
133
Это позволяет составить шесть уравнений с шестью неизвестными, которыми являются плюккеровы координаты единичного винта V (x, y, z, Mx, My, Mz). Пять из этих уравнений представляют собой равенства нулю относительных моментов V1*e2 = 0; V1*e3 = 0; V1*e4 = = 0; V1*e5 = 0; V1*e6 = 0, а шестое уравнение
нормирует модуль единичного винта V: 2 2 2 „ 1 x + y + z = 1. В рассматриваемом примере
эти уравнения принимают вид:
Vj * e2 = x • b sin у + y • 0 + z • (-b cos у) +
+Mx • (- cos Y) + My • 0 + Mz • (- sin Y) = 0;
Vj * e3 = x • b + y • 0 + z • 0 +
+Mx • 0 + My • 0 + Mz • (-1) = 0;
Vj* e4 = x • 0 + y • 0 + z • 0 + Mx • 0 +
+My • (- cos в) + Mz • (- sin в) = 0;
Vj* e5 = x • 0 + y • 0 + z • 0 +
+Mx • 0 + My • 0 + Mz • (-1) = 0;
V * e6 = x • 0 + y • (-c sin Y) +
+ z • 0 + Mx • cos Y + My • 0 +
+Mz • (- sin Y) = 0;
2,2,2 3 x + y + z = 1 .
В данных уравнениях значком «*» обозначена операция вычисления относительного момента двух винтов.
Решая полученную систему уравнений, находим плюккеровы координаты искомого кинематического винта:
V1 =(0 - cos в sin в -a cos в 0 0 ) .
Для определения усилия S1 по принципу возможных скоростей составим выражение для суммы возможных мощностей сил, действующих на систему, на кинематическом винте V1, т. е. найдем относительный момент
силового винта заданных сил Р вместе с неизвестным усилием S1:
P
P sin а
-P cos а ^
-Q - G - S1 G- + Qb + S1b
-GC - Q— - S1c \ 2 2 1
P cos a c
и кинематического винта V1:
-a cos в - P sin а + cos в x
f c c Л < G — + Q— + S1c + sinвх
V 22
xP cos а^ c = 0;
^ S1 = -1550 Н.
Заметим, что ось кинематического винта V1 направлена вдоль пунктирной прямой шх (рис. 2), которая совпадает с осью Риттера для стержня 1 (сумма моментов всех усилий S2-S6 относительно этой оси тождественно равна нулю). Нужно отметить, что в данном случае положение прямой шх даже без расчетов легко определяется из геометрических соображений.
Пользуясь аналогичным алгоритмом, находим плюккеровы координаты кинематических винтов V2 и V3, характеризующих возможные движения платформы при удалении опорных стержней 2 и 3:
V2 = (0 0 1 0 0 0); (1)
V3 =
f - sin а - cos а V-a cos а 0
0 ^ 0 ,
(2)
Положение мгновенной винтовой оси m2 платформы в момент удаления стержня 2, определяемое координатами (1), представлено на рис. 3. Как и в первом случае, его легко можно определить геометрическим способом. Для этого необходимо найти прямую линию, которая параллельна стержням 1, 3 и 5 и одновременно пересекает стержни 4 и 6.
ISSN 1815-588Х. Известия ПГУПС
2015/3
134
Общетехнические задачи и пути их решения
Рис. 3. Положение мгновенной винтовой оси m1
Рис. 5. Схема платформы, опирающейся на упругие стержни
Плюккеровы координаты (2) определяют мотор мгновенной оси вращения платформы в момент удаления стержня 3 (расстояние между мгновенной осью вращения и мгновенной винтовой осью m3 определяется по формуле d = ab2/(b2 + c2) [3]). При этом возможность найти геометрическим способом положение мгновенной винтовой оси m3 платформы (рис. 4) представляется не столь очевидной, как в двух предыдущих случаях.
Рис. 4. Положение мгновенной винтовой оси m2
в предыдущем примере, опирается на упругие в продольном направлении стержни Sx-S6 с коэффициентами жесткости cx = 105 Н/м, С2 = 2-10 4 Н/м, c3 = 3104 Н/м, c4 = 4 104 Н/м, c = 5 104 Н/м и c = 6 104 Н/м, соответственно. На рис. 5 система изображена в положении покоя. В первом приближении пренебрежем сопротивлением движению звеньев механизма.
Выберем систему координат Oxyz, начало которой совместим с центром масс O платформы, а оси x, у и z - с ее главными центральными осями инерции. За обобщенные координаты системы примем смещения x, y и z центра масс платформы и углы фх, фу и фz ее поворотов вокруг осей х, у и z, соответственно.
Пренебрежем кинетической энергией стержневых опор 1-6, сообщим системе обобщенные скорости X, у, z, фх, ф и фz и запишем выражение кинетической энергии системы:
Пользуясь изложенным алгоритмом, последовательно находим усилия во всех оставшихся опорных стержнях: S2 = 894,43 Н; S3 =
= -700 Н; S4 = -721,11 Н; S5 = -850 Н; S6 = 0 Н. Г х >
Применение винтового исчисления у z
к определению частот малых свободных колебаний платформы где q = ф X
ф у
Рассмотрим случай, когда однородная пря- 1ф z )
моугольная платформа (рис. 5), рассмотренная скоростей;
1 т
т=2 qт Aq,
вектор-столбец обобщенных
2015/3
Proceedings of Petersburg Transport University
Общетехнические задачи и пути их решения
135
( m 0 0 0 0 0 Л
0 m 0 0 0 0
A = 0 0 m 0 0 0
0 0 0 Jx 0 0
0 0 0 0 Jy 0
10 0 0 0 0 Jz у
инерции; m - масса платформы
матрица
J =
m(f2 + b2) 12
Jy =
m( f2 + с2) 12
r m(c2 + b2) ^
Jz =---—------моменты инерции платфор-
мы относительно ее главных центральных осей инерции.
Потенциальную энергию системы можно представить в матричной форме:
Из этих строк сформируем вспомогательную матрицу Pluck:
Pluck =
0 0 -1
-cosy 0 - ■sin у
0 0 -1
0 —cos р - -sinp
0 0 -1
cosy 0 - ■sin у
0,56
0,56 sin у
0,56
1 о СЛ о sin Р - 0,5/ cosp
-0,56
-0,56 sin у
1 T
п=2 q Cq,
где C - матрица инерции системы.
Запишем плюккеровы координаты единичных винтов et-e6, приведенных к началу O выбранной системы координат Oxyz (рис. 5):
е1 = (0 0 -1 0,5b -0,5с 0);
Є2 _
- cos у 0,5b sin у
0,5с sin y + 0,5 f cos у
0 Л
- sin у 0,5b cos уу
-0,5c
0,5c siny + 0,5 f cosy 0,5c
0,5c sin P 0,5c
-0,5csiny-0,5/ cosy 0
-0,56 cos у 0
-0,5c cos P 0
\ -0,56 cos у
Є4 =
е3 = (0 0 -1 0,5b 0,5с 0);
( 0 - cos P ^
- sin p -0,5b sin p-0,5 f cos P
Введем в рассмотрение диагональную матрицу С коэффициентов жесткости упругих стержней
е6 =
v 0,5с sin р -0,5с cos Р у ( с 1 0 0 0 0 0
е5 = (0 0 -1 -0,5b 0,5с 0); 0 с2 0 0 0 0
0 0 с3 0 0 0
( cos у 0 Л C1 = 0 0 0 с4 0 0
- sin у -0,5b sin у 0 0 0 0 с5 0
v-0,5с sin у - 0,5 f cos у -0,5b cos уу 10 0 0 0 0 сбу
ISSN 1815-588Х. Известия ПГУПС
2015/3
136
Общетехнические задачи и пути их решения
Искомая матрица жесткости системы определяется формулой
C = (C1Pluck)T Pluck .
Подстановка выражений для T и П в уравнения Лагранжа II рода приводит к известному дифференциальному уравнению свободных колебаний системы с n степенями свободы в матричной форме:
q + A-1Cq = 0.
Для численного решения задачи воспользуемся стандартным пакетом математических программ Mathcad. Вычисления проведем при следующих исходных данных: m = 254,84 кг, J = 191,34 кг-м2, J = 340 кг-м2, J = 530,92 кг-м2, a = 2 м, b = 3 м, с = 4 м, f= 0,1 м.
Собственные частоты системы можно найти с привлечением стандартной процедуры определения собственных чисел матрицы A-1C:
k =sjeigenvals(A 'C).
Получаем следующие значения частот свободных колебаний платформы:
kT= (6,54; 11,35; 16,49; 31,06;
35,84; 62,44) с-1.
Матрица собственных форм колебаний системы может быть записана в виде:
M =
1 1 1
75,734 0,001 -0,856
-174,373 0,898 -3,986
1297,895 -2,497 -0,479
-944,475 -2,084 -0,091
73,128 -0,370 -1,853
1 1 1
-3,222 0,296 2,095
0,065 0, 207 -3,255
0,157 0,318 -1,800
0,026 0, 228 -0,401
0,891 -0,057 3,336
Заключение
Приведенные в статье примеры применения винтового исчисления к решению задач статики и динамики пространственного механизма параллельной структуры (МПС) иллюстрируют некоторые особенности использования элементов теории винтов при решении подобных задач, а также показывают, что подобная методика позволяет существенно уменьшить трудоемкость расчетов, связанных с определением основных параметров, характеризующих свободные колебания МПС.
Библиографический список
1. Диментберг Ф. М. Теория винтов и ее приложения / Ф. М. Диментберг. - М. : Наука, 1978. -327 с.
2. Диментберг Ф. М. Теория пространственных шарнирных механизмов / Ф. М. Диментберг. - М. : Наука, 1982. - 336 с.
3. Доронин Ф. А. К вопросу об определении положения мгновенной оси вращения / Ф. А. Доронин // Бюл. результатов научных исследований. -2013. - Вып. 3 (8). - С. 23-35.
4. Занчевский И. М. Теория винтов и приложения ее к механике / И. М. Занчевский. - Одесса, 1889. - 157 с.
5. Зейлигер Д. Н. Комплексная линейчатая геометрия / Д. Н. Зейлигер. - М. : Гостехиздат, 1934. -196 с.
6. Котельников А. П. Винтовое счисление и некоторые приложения его к геометрии и механике / А. П. Котельников. - Казань, 1895. - 216 с.
7. Ball R. S. The Theory of Screws : A study in the dynamics of a rigid body / R. S. Ball. - Ireland, Dublin : Hodges, Foster and Co., 1876. - 194 p.
8. Gough V. E. Contribution to discussion of papers on research in Automobile Stability, Control and Tyre performance / V. E. Gough. - Proc. Auto Div. Inst. Mech. Eng., 1956-1957. - P. 392-394.
9. Mises R. Anwendungen der Motorrechnung / R. Mises // Zeitschr. fur angew. Math. und Mech. -1924. - B. 4, H. 3. - S. 193-213.
2015/3
Proceedings of Petersburg Transport University
Общетехнические задачи и пути их решения
137
10. Mises R. Motorrechnung, ein neues Hilfmittel der Mechanik / R. Mises // Zeitschr. fur angew. Math. und Mech. - 1942. - B. 4, H. 2. - S. 155-181.
11. Stewart D. A Platform with Six Degrees of Freedom / D. Stewart. - Proc. Inst. Mech. Eng., 1965. -P. 1-15.
12. Study Е. Von Bewegungen und Umlegun-gen / E. Study. - M. A. XXXIX Band., 1891.
13. Study Е. Geometrie der Dynamen / E. Study. -Lpz., 1901.
References
1. Dimentberg F. M. Teoriya vintov i yeye prilo-zheniya [Screw Theory and Its Applications]. Moscow, Nauka, 1978. 327 p.
2. Dimentberg F. M. Teoriya prostranstvennykh sharnirnykh mekhanizmov [Theory of Spatial Jointed Mechanisms]. Moscow, Nauka, 1982. 336 p.
3. Doronin F. A. Byulleten rezultatov nauch-nykh issledovaniy - Bull. Res. Results, 2013, no. 3 (8), pp. 23-35.
4. Zanchevskiy I. M. Teoriya vintov i prilozheniya yeye k mekhanike [Screw Theory and Its Applications in Mechanics]. Odessa, 1889. 157 p.
5. Zeylinger D. N. Kompleksnaya lineychataya geometriya [Complex Ruled Geometry]. Moscow, Gostekhizdat, 1934. 196 p.
6. Kotelnikov A. P. Vintovoye schisleniye i neko-toryye prilozheniya yego k geometrii i mekhanike [Screw Calculus and Some of Its Applications in Geometry and Mechanics]. Kazan, 1895. 216 p.
7. Ball R. S. The Theory of Screws: A study in the dynamics of a rigid body. Ireland, Dublin, Hodges, Foster and Co., 1876. 194 p.
8. Gough V. E. Contribution to discussion of papers on research in Automobile Stability, Control and Tyre performance. Proc. Auto Div. Inst. Mech. Eng., 1956-1957, pp. 392-394.
9. Mises R. Anwendungen der Motorrechnung. Zeitschr. fur angew. Math. und Mech., 1924, B. 4, H. 3, ss. 193-213.
10. Mises R. Motorrechnung, ein neues Hilfmittel der Mechanik. Zeitschr. fur angew. Math. und Mech., 1942, B. 4, H. 2, ss. 155-181.
11. Stewart D. A Platform with Six Degrees of Freedom. Proc. Inst. Mech. Eng., 1965, pp. 1-15.
12. Study Е. Von Bewegungen und Umlegungen. M. A. XXXIX Band., 1891.
13. Study Е. Geometrie der Dynamen. Lpz., 1901.
ДОРОНИН Феликс Александрович - канд. техн. наук, доцент, [email protected] (Петербургский государственный университет путей сообщения Императора Александра I).
ISSN 1815-588Х. Известия ПГУПС
2015/3
