Транспортное машиностроение. 2022. № 12(12). С. 36-47. ISSN 2782-5957 (print) Transport Engineering. 2022. no. 12(12). P. 36-47. ISSN 2782-5957 (print)
Научная статья Статья в открытом доступе УДК 621.01
doi: 10.30987/2782-5957-2022-12-36-47
АНАЛИЗ СОБСТВЕННЫХ ЧАСТОТ И ФОРМ КОЛЕБАНИЙ ПЛОСКОГО МАНИПУЛЯТОРА ПАРАЛЛЕЛЬНОЙ СТРУКТУРЫ
Андрей Константинович Толстошееви
Брянский государственный технический университет, Брянск, Россия andrei.tolstosheev@yandex.ru, https://orcid.org/0000-0002-9417-775X
Аннотация
Цель исследования заключается в анализе собственных частот и форм малых колебаний плоского манипулятора параллельной структуры и определении влияния упругих и инерционных свойств промежуточных звеньев и привода на собственные частоты и формы. Статья посвящена решению задачи обеспечения виброустойчивости манипуляторов. Использованы аналитические методы теории колебаний и вычислительные эксперименты в программном комплексе «Универсальный механизм». Новизна работы состоит в определении зависимости собственных частот и форм колебаний плоского манипулятора параллельной структуры от соотношения упругих и инерционных свойств звеньев и привода. В результате исследования разработаны аналитическая и несколько компьютерных моделей плоского манипулятора с тремя степенями свободы, определены собственные частоты и формы колебаний для разработанных моделей. Выводы: замена абсолютно жёстких промежуточных звеньев упругими стержнями привела
Ссылка для цитирования:
Толстошеее А.К. Анализ собственных частот и форм колебаний плоского манипулятора параллельной структуры / А.К. Толстошеее // Транспортное машиностроение. - 2022. - № 12. - С. 36 - 47. doi: 10.30987/27825957-2022-12-36-47.
Original article Open Access Article
ANALYSIS OF OWNFREQUENCIES AND OSCILLATION FORMS OF A PLANE PARALLEL MANIPULATOR
Andrey Konstantinovich Tolstosheev^
Bryansk State Technical University, Bryansk, Russia andrei.tolstosheev@yandex.ru, https://orcid.org/0000-0002-9417-775X
Abstract
The study objective is to analyze the own frequencies and small-amplitude oscillations of a plane manipulator having parallel structure and to determine the influence of elastic and inertial properties of intermediate links and the drive on own frequencies and forms. The paper is devoted to solving the problem of ensur-
к появлению собственных частот и соответствующих им собственных форм, на которых выходное звено колеблется в вертикальной плоскости, перпендикулярной плоскости движения звеньев манипулятора; колебания выходного звена в вертикальной и горизонтальной плоскостях на собственных формах разделены; значения собственных частот, на которых выходное звено колеблется в вертикальной плоскости, не зависят от коэффициента жёсткости привода, а определяются упругими и инерционными свойствами стержней; значения собственных частот, на которых выходное звено колеблется в горизонтальной плоскости, зависят от соотношения податливостей привода и стержней; влияние масс упругих звеньев на значения собственных частот определяется видом собственных форм и соотношением инерционных параметров выходного звена и стержней.
Ключевые слова: манипулятор, структура, колебания, частоты и формы, моделирование.
ing vibration stability of manipulators. Analytical methods of the theory of oscillations and computational experiments in the software package "Universal Mechanism" are used. The novelty of the work is in defining the dependence of the own frequencies and oscillation forms of a plane parallel manipulator on the ratio of
36
© Толстошеев А.К., 2022
elastic and inertial properties of the links and the drive. As a result of the research, analytical and several computer models of a plane manipulator with three degrees of freedom are developed, own frequencies and oscillation forms for the developed models are determined. Conclusions: the replacement of absolutely rigid intermediate links with elastic rods led to the appearance of own frequencies and own modes corresponding to them , on which the output link oscillates in the vertical plane perpendicular to the plane of the manipulator links movement; the oscillations of the output link in the vertical and horizontal planes on their own forms are separated; the values of own frequencies at which
the output link oscillates in the vertical plane do not depend on the coefficient of rigidity of the drive, but are determined by the elastic and inertial properties of the rods; the values of own frequencies at which the output link oscillates in the horizontal plane depend on the ratio of the drive and rods compliance; the influence of the elastic link masses on the values of their own frequencies is determined by the type of their own forms and the ratio of inertial parameters of the output link and rods.
Keywords: manipulator, structure, oscillations, frequencies and forms, modeling.
Reference for citing:
Tolstosheev AK. Analysis of own frequencies and oscillation forms of a plane parallel manipulator. Transport Engineering. 2022; 12:36-47. doi: 10.30987/2782-5957-2022-12-36-47.
Введение
Манипуляторы нового типа с параллельной структурой благодаря высоким показателям по быстродействию, грузоподъёмности и точности находят всё большее применение в различных транспортных и технологических машинах [1, 2]. Одной из задач, возникающих при проектировании манипуляторов, является обеспечение виброустойчивости их конструкции. Решение этой задачи связано с определением собственных частот и форм малых колебаний манипулятора. На движение выходного звена манипулятора как параллельной, так и последовательной структуры, существенное влияние оказывают колебания, обусловленные упругостью промежуточных звеньев кинематической цепи [3, 4]. Однако при определении собственных частот колебаний манипуляторов параллельной структуры используют, как правило, аналитические методы и упрощённые расчётные схемы, в которых упругие и инерционные свойства промежуточных звеньев не учитывают. Так, в работе [5] на примере плоского механизма разработана методика расчёта аналитическим методом собственных частот колебаний манипулятора с параллельной структурой, в которой учитывают инерционные свойства выходного звена и упругие свой-
ства приводов, а промежуточные звенья считают невесомыми и абсолютно жёсткими. Предложенная методика не предусматривает определение собственных форм колебаний. В работе [6] с использованием данной методики выполнен расчёт собственных частот пространственного поступательно направляющего манипулятора с тремя степенями свободы, у которого выходное звено описывается как материальная точка. Аналогичные допущения приняты при расчёте свободных и вынужденных колебаний поступательно направляющих манипуляторов иной параллельной структуры [7-9].
Настоящая работа посвящена анализу собственных частот и форм малых колебаний плоского манипулятора параллельной структуры и определению влияния упругих и инерционных свойств промежуточных звеньев, а также податливости привода на собственные частоты и формы. Для вычислительных экспериментов использован программный комплекс «Универсальный механизм» [10, 11]. С целью подтверждения достоверности численных результатов в работе выполнен сравнительный расчёт собственных частот колебаний манипулятора аналитическим и численным методами.
Материалы, модели, эксперименты и мет Аналитическое моделирование.
Аналитический расчёт собственных частот колебаний манипулятора с абсолютно жёсткими звеньями выполняем в соответ-
ствии с методикой [5]. Рассмотрим плоский манипулятор с тремя степенями свободы в фиксированном положении в гори-
зонтальнои плоскости xoy и заторможенными приводами (рис. 1).
Рис. 1. Кинематическая схема манипулятора Fig. 1. Kinematic diagram of the manipulator
Начало неподвижной декартовой системы координат расположим в центре масс выходного звена, изображённого в виде треугольника A1A2A3, массой m и моментом инерции I. Массой промежуточных звеньев, выполненных в виде стержней одинаковой длины L, пренебрегаем. В качестве обобщённых координат qi, q2, q3 выбираем углы поворота входных звеньев B1C1, B2C2, B3C3. Упругость приводов учитываем заданием в приводных кинематических парах восстанавливающих пар сил с коэффициентом жёсткости с. Для составления дифференциальных уравнений сво-
бодных колебаний манипулятора используем уравнение Лагранжа второго рода
где П = с(ц1 + + ц^)/2 , Т = (тх2 + ту2 + 1ф2)/2 - соответственно потенциальная и кинетическая энергия системы; х, у ,ф - проекции скорости центра масс выходного звена на оси координат и его угловая скорость.
Вектор скоростей выходного звена и и вектор обобщённых скоростей ^ связаны между собой соотношением:
й = Уц, (2)
где и = (х, у, ф)Т; ц = ^, ц3)т; У -матрица Якоби, У = А-1 В, где
A =
/dF1 dF1 dF^x дх ду д<р dF2 дх dF3
dF2 ду
dF2 д<р
\ дх ду dF.
dF, dF,
B = -
/0*1 ' aqi
0
0 IT
dq2
dtp/ 0
0
\
0
0 — I
дцз/
^ - функции, связывающие в неявной форме обобщённые координаты и перемещения выходного звена,
Pi = (*Ai - xBi-L cos qi)2 + (yAi — yBi-L sin qt)2, i = 1, 2, 3, где, в частности, при i = 1 имеем хА1 = х + Rt sin ф ; уА1 = у — R\ cos ф.
Принимаем радиусы Я = 0,1 м и Я2 = 0,2 м, длины стержней Ь = 0,15 м, координаты точек В1, В2, Вз соответственно (0; -0,2), (0,173; 0,1), (- 0,173; 0,1), точек А1, А2,
A3 соответственно (0; - 0,1), (0,0866; 0,05), (- 0,0866; 0,05). Тогда матрицы A, B и Y имеют вид
0,2828 0,100
A =
0,0548 0,2282
0,2828n
0,02828
0
-0,2950 -0,2829); B = 0,1949 -0,2829
0,0592 -0,0115
0 -0,02829 00 0,0478n
Y = (-0,021 0,0618 -0,0408 ,0,3334 0,3333 0,3333
0 0
0,02829
Подставив значения величин с учётом выражения (2) в уравнение (1), получим в матричной форме систему дифференциальных уравнений свободных колебаний манипулятора с абсолютно жёсткими звеньями в горизонтальной плоскости:
Mq + Cq = 0,
(3)
где ц = Ц2, ЧъУ, Ч = 42, Чъ)Т, M и C - соответственно матрицы инерционных и упругих коэффициентов.
При массе платформы m = 1 кг, моменте инерции I = 0,01 кгм2 и коэффициентах жёсткости приводов c = 100 Н-м матрицы M и C имеют вид
M = 10-
/ 5,058 -0,8676 -0,8616\ /100 (-0,8674 5,062 -0,8606); C = ( 0 \-0,8618 -0,8608 5,060 / V 0
0 100 0
0 0 100
Частные решения системы дифференциальных уравнений (3) ищем в виде
q¿ = Ai sin wt, í = 1, 2 3, (4)
Подставив решения (4) в уравнения (3), получим систему однородных алгебраических уравнений относительно амплитуд А1, A2, Aз:
(100000 - 5,058^2)Л1 + 0,8676^2Л2 + 0,8616^2Л3 = 0^ 0,8674^2Л1 + (100000 - 5,062^2)Л2 + 0,8606^2Л3 = 0 0,8618^2Л1 + 0,8608^2Л2 + (100000 - 5,060^2)Л3 = 0
(5)
Приравняв к нулю определитель, составленный из коэффициентов при Al, A2, Aз, получим частотное уравнение
- 638 • 102^4 + 130 • 107^2 - 855 • 1010 = 0. (6)
Корнями частотного уравнения являются круговые частоты свободных колебаний.
Компьютерное моделирование. При построении компьютерных моделей для вычислительных экспериментов манипулятор представляется в трёхмерном пространстве системой твёрдых тел, образующих между собой подвижные соединения. В соответствии с требованиями программного комплекса «Универсальный механизм» принимаем, что все твёрдые тела в кинематической цепи механизма обладают ненулевыми инерционными параметрами, которые или задаются в виде массы и тензора инерции или вычисляются в программе по созданным графическим образам.
Каждое промежуточное упругое звено (стержень) в кинематической цепи манипулятора рассматривается как совокупность п абсолютно твёрдых тел, образующих между собой четырёхподвижные соединения. Каждое такое соединение накладывает на относительное движение тел две геометрические связи, запрещая линейные перемещения в плоскости, перпендикулярной продольной оси стержня. В
четырёхподвижное соединение по каждой обобщённой координате i введён линейный упругий элемент с коэффициентом жёсткости Ci. Коэффициентам жёсткости присваиваем следующие значения:
С1 = с2 = си, с3 = ск, с4 = сп
Коэффициенты крутильной Ск, продольной Сп и изгибной Си жёсткости рассчитываются по известным формулам
_ С/р _ £5 _
Ск = г ,Сп = г ,Си = г , (7)
где E и G - соответственно модуль упругости первого и второго рода материала стержня; £ - площадь поперечного сечения стержня; I - длина абсолютно жёсткого элемента стержня, I = Ып; 3 - полярный и осевой моменты инерции поперечного сечения стержня.
В работе принято, что стержни могут быть выполнены из стали или из алюми-
ния, имеют круглое сплошное или кольцевое поперечное сечение.
Для круглого сечения диаметром Б
Масса mi и главные центральные моменты инерции каждого абсолютно жёст-
где р - плотность материала стержня; х -ось местной системы координат, проходящей вдоль стержня. Для стержня с круглым поперечным сечением момент инерции Лз = miD1l8. Принимаем для стали Е = 210 ГПа, G = 80 ГПа, р = 7800 кг/м3; для алюминия - Е = 71 ГПа, G = 26 ГПа, р = 2700 кг/м3.
В тестовом примере программного комплекса «Универсальный механизм» при расчёте собственных частот изгибных колебаний консольная балка разделялась на десять подвижных частей, и первые пять собственных частот определялись с относительной погрешностью менее 2,5 % по сравнению с аналитическим решением [12]. Поэтому принимаем в компьютерной модели каждого упругого стержня в кинематической цепи механизма п = 10.
Анализ собственной формы на каждой собственной частоте выполняем с использованием анимационного окна программы, в котором изображается движение графического образа механизма.
В работе все расчёты с использованием программного комплекса «Универсальный механизм» выполняются для двух геометрически подобных кинематических схем манипулятора (см. рис. 1). Одна схема выполнена с короткими стержнями (Ь = 0,15 м), а другая - с длинными (Ь = 1,5 м). В кинематических схемах для устранения избыточных связей согласно [13, 14] шарниры В и С1 описаны как вращательные кинематические пары, а шарниры Ai - как трёхподвижные сферические пары. На основе этих кинематических схем построены восемь компьютерных моделей манипулятора, отличающихся геометрическими,
5=^£! ,;=££ ,/р = 2/. (8)
Для трубки с внешним диаметром Б и внутренним диаметром ё
(9)
кого ьго элемента стержня в виде трубки определяются по известным формулам
(10)
упругими и инерционными свойствами звеньев. Модели № 1 и № 2 применяются для расчёта манипуляторов с абсолютно жёсткими невесомыми стержнями и сравнения численного решения с аналитическим, модели № 3 и № 4 используются для анализа спектра собственных частот и оценки влияния на собственные частоты манипулятора коэффициентов жёсткости привода, модели № 5-8 применяются для анализа влияния упругих и инерционных свойств стержней. Во всех компьютерных моделях для выходного звена задаётся масса т и тензор инерции в виде диагональной матрицы Т, для которой элементы
¿11 = ¿22 = tзз = I.
Модель № 1 соответствует плоской расчётной схеме с короткими стержнями, принятой в аналитическом решении: для выходного звена т = 1 кг, I = 0,01 кг-м2; коэффициенты жёсткости приводов одинаковы с = 100 Н-м; стержни абсолютно жёсткие малоинерционные с пренебрежимо малыми массой (10-7 кг) и тензором инерции (10-9 кг-м2).
Модель № 2 соответствует модели № 1, но кинематическая схема выполнена с длинными стержнями.
Модель № 3 выполнена на основе модели № 1, но имеет следующие отличия: стержни манипулятора выполнены в виде малоинерционных алюминиевых трубок диаметрами Б = 26 мм и ё = 12 мм, коэффициенты жёсткости приводов варьируются от 1 до 1-1010 Н-м.
Модель № 4 выполнена на основе модели № 2, но имеет следующие отличия: стержни выполнены в виде алюминиевых трубок с диаметрами Б = 50 мм и ё =
^ср = — .-Ь = 4
■'=£ (1- £)■*■ = 2/ ■
т1*2 I _ т1°ср
= Р^, = = 12 , ^ =
4
30 мм. Для каждого абсолютно жёсткого элемента стержня длиной и = 0,15 м задаются масса mi = 0,509 кг и ненулевые диагональные элементы тензора инерции Т: ¿11 = 1х = 2,04-10-4 кг-м2, ¿22 = 1у = 9,54-10-4 кг-м2, ¿33 = Ь = 9,54-10-4 кг-м2, для выходного звена т = 2 кг и I = 0,01 кг-м2, коэффициенты жёсткости приводов варьируются от 1 до 1^1010 Н-м.
Модель № 5 выполнена на основе модели № 1, имеет следующие отличия: с = 1000 Н-м, малоинерционные алюминиевые стержни манипулятора выполнены со сплошным круглым сечением
Результаты
Сравнение собственных частот и форм для разных моделей. Для аналитической модели манипулятора, решая частотное уравнение (6), получим круговые частоты: Ю1 = 124,9 рад/с, Ю2 = 135,6 рад/с, Ю3 = 172,7 рад/с. Собственные частоты колебаний /1А = 19,89 Гц, /2А = 21,59 Гц, /за = 27,50 Гц.
Для компьютерной модели № 1 собственные частоты манипулятора с малоинерционными стержнями /1 = /х = 20,59 Гц, /2 = /у = 20,74 Гц, /3 = /и = 27,53 Гц. Здесь и далее /х, /у, /щ - собственные частоты, ко-
диаметром Б, который варьируется от 8 до 44 мм.
Модель № 6 аналогична модели № 5, но массы и тензоры инерции стержней вычисляются в программе автоматически.
Модель № 7 выполнена на основе модели № 2, но имеет следующие отличия: стальные малоинерционные стержни имеют круглое поперечное сечение диаметром Б, который варьируется от 20 до 100 мм, для выходного звена т = 2 кг и I = 0,1 кгм2, с = 10000 Н-м.
Модель № 8 аналогична модели № 7, но инерционные параметры стержней вычисляются в программе автоматически.
торым соответствуют собственные формы х, у, фг. На собственной форме х выходное звено совершает возвратно-поступательное движение в горизонтальной плоскости с частотой /х с максимальной амплитудой вдоль оси х неподвижной системы координат. Аналогично, на собственной форме у выходное звено колеблется с частотой /у с максимальной амплитудой вдоль оси у. На собственной форме фг выходное звено колеблется вокруг вертикальной оси г. Максимальная относительная погрешность численного решения
Д/= (/2л - /2)100 %//"2Л = (21,59 - 20,74)100 %/21,59 = 3,9 %.
Для компьютерной модели № 2, имеющей геометрически подобную кинематическую схему с увеличенными в 10 раз длинами стержней, собственные частоты /х = 2,059 Гц и /у = 2,074 Гц уменьшились в 10 раз, а частота /фЪ = 27,53 Гц не изменилась.
Увеличение в 10 раз момента инерции I, при сохранении без изменения остальных параметров компьютерной модели № 2, привело к уменьшению в 3,16 раза частоты /фЪ = 8,71 Гц. Остальные частоты /х = 2,059 Гц и /у = 2,074 Гц остались без изменения, что подтверждает независимость поступательного и вращательного движения выходного звена на собственных формах для рассмотренного положения механизма с малоинерционными жёсткими звеньями.
Для компьютерной модели № 3 с упругими стержнями при с = 100 Н-м общее число собственных частот равно 219. Из них первые шесть низших собственных частот имеют значения: /х = 20,58 Гц, /у = 20,85 Гц, /щ = 27,51 Гц, /фх = 137,7 Гц, /Ру = 138,0 Гц, /г = 195,1 Гц, которые существенно отличаются от ближайшей большей частоты /7 = 35565 Гц. Здесь и далее / фх и /фу - собственные частоты, которым соответствуют собственные формы фх и фу, на которых выходное звено имеет максимальные амплитуды колебаний соответственно вокруг осей х и у. На собственной форме г выходное звено совершает возвратно-поступательное движение с частотой/2 вдоль оси г. Замена абсолютно жёстких стержней упругими привела к появлению собственных частот и форм, на кото-
рых выходное звено колеблется в вертикальной плоскости. Колебания выходного звена в вертикальной и горизонтальной плоскостях на собственных формах разделены. Разделены также вращательные и поступательные движения выходного звена в каждой плоскости, а колебания на собственных формах с близкими частотами/ и/у, а такжеи/<ру взаимосвязаны.
Для компьютерной модели № 4 с упругими инерционными стержнями при с = 10000 Н-м спектр первых девяти собственных частот более плотный чем в модели № 3: / = 3,97 Гц, /у = 3,98 Гц, = 4,56 Гц, / = 8,56 Гц, /фХ = 10,03 Гц, /9у = 10,04 Гц, /?2 = 13,9 Гц, /фх2 = 14,569 Гц, /фу2 = 14,572 Гц. Далее следуют три частоты 52,8 Гц, 53,7 Гц и 53,8 Гц, на собственных формах которых выходное звено практически неподвижно, а стержни совершают изгибные колебания в горизонтальной плоскости. На следующих трёх собственных формах с частотами /Щ2 = 64,8 Гц, /2 = 66,2 Гц, /2 = 66,3 Гц выходное звено колеблется в горизонтальной плоскости вследствие изгиба стержней и с меньшими амплитудами, чем на частотах /фг, /х, /у. В конце спектра, начиная с частоты /175 = 5997 Гц и заканчивая /219 = 7126 Гц, собственные формы определяются продольными колебаниями стержней.
Численные эксперименты с варьированием коэффициента жёсткости привода. В процессе вычислительных экспериментов на модели № 3 коэффициенты жёсткости привода изменялись от 1 до 1-1010 Н •м. Установлено, что собственные частоты / / фх, /у, соответствующие собственным формам с колебаниями выходного звена в вертикальной плоскости, не зависят от величины коэффициента жёсткости привода.
Собственные частоты /х, /у, /фг, на которых выходное звено колеблется в горизонтальной плоскости, зависят от величины коэффициента жёсткости привода аналогичным образом (табл. 1 ). Полулогарифмический график /Х(1ё с) (рис. 2) имеет три участка с различной кривизной, что можно объяснить свойствами последовательного соединения двух упругих элементов: привода и стержней, связывающих
выходное звено и неподвижное основание в каждой из трёх упругих связей. Общая податливость последовательного соединения элементов есть сумма податливостей привода и стержней. На первом участке, для которого значения коэффициента жёсткости привода малы с < 1000 Н-м (1§ с < 3), суммарная податливость системы определяется податливостью привода, а податливость стержней пренебрежимо мала. Это объясняет, почему при с =100 Н-м в модели № 3 с упругими звеньями получены значения собственных частот /х = 20,58 Гц, /у = 20,85 Гц, /фг = 27,51 Гц, мало отличающиеся от значений для модели № 1 с жёсткими звеньями. На третьем участке при 1§ с > 6 податливость привода пренебрежимо мала по сравнению с податливостью стержней, поэтому дальнейшее увеличение жёсткости и, соответственно, уменьшение податливости привода не изменяет собственные частоты колебаний. На промежуточном втором участке податливости стержней и привода сопоставимы. Для сравнения, коэффициенты жёсткости в соединении двух абсолютно жёстких элементов упругого стержня в соответствии с формулами (7), (9) си = 1,01-105 Н-м, ск = 7,42-104 Н-м, сп = 1,98-109 Н/м.
Рис. 2. Зависимость /i(lg с)
Fig. 2. Dependence of fx(lg с)
Для манипулятора с абсолютно жёсткими звеньями собственные частоты постоянно растут с ростом коэффициента жёсткости привода. Например, при с = 1-109 Н •м для модели № 1 fx = 65107 Гц, fy = 65582 Гц, fpz = 87073 Гц. Второй и третий участки графика fx(lg с) в этом случае отсутствуют.
Аналогичная картина сохраняется и для модели № 4 при других значениях геометрических, инерционных и упругих параметров промежуточных звеньев. Низшие собственные частоты с собственными формами, соответствующими колебаниям выходного звена в вертикальной плоско-
сти, равны fz = 8,56 Гц, /щ = 10,03 Гц, /фу = 10,04 Гц, не зависят от величины коэффициента жёсткости привода. Значения собственных частот для двух собственных форм с колебаниями выходного звена в горизонтальной плоскости приведены в табл. 1.
f Собственная форма Номер модели № lg c, Н>м
1 2 3 4 5 6 7 8 9 10
fx X 3 6,91 20,6 64,2 181 330 376 382 383 383 383
4 0,15 0,434 1,36 3,97 7,89 9,36 9,56 9,58 9,58 9,58
9z 3 8,22 27,5 86,0 243 442 504 512 513 513 513
4 0,16 0,50 1,60 4,60 9,10 10,9 11,2 11,2 11,2 11,2
Таблица 1
Значения низших собственных частот f для моделей № 3 и № 4, Гц
Table 1
Values of the lowest natural frequencies f for models No. 3 and No. 4, Hz
Численные эксперименты с варьированием диаметра сечения стержней.
На первом этапе для оценки влияния только упругих свойств стержней выполняем расчёты на модели № 5 при с = 1000 Н-м для разных значений диаметра Б коротких алюминиевых стержней круглого поперечного сечения с пренебрежимо малыми массой и тензором инерции. В соответствии с формулами (7), (8) при увеличении диаметра увеличиваются коэффициенты изгибной, крутильной и продольной жёст-
кости. Массы и тензоры инерции звеньев, при этом, остаются без изменения. На втором этапе для оценки влияния инерционных свойств стержней выполняем аналогичные расчёты на модели № 6, для которой масса и тензор инерции каждого жёсткого элемента стержня автоматически рассчитываются в программе по формулам (10) и зависят от диаметра. Результаты вычислительного эксперимента для моделей № 5 и № 6 приведены в табл. 2.
Таблица 2
Значения низших собственных частот f для моделей № 5 и № 6, Гц Values of the lowest natural frequencies f for models No. 5 and No. 6, Hz
Table 2
f Собственная форма Номер модели № Диаметр стержня D, мм
8 12 16 20 24 28 32 36 40 44
fx x 5 32,5 51,6 59,7 62,7 63,9 64,5 64,7 64,9 65,0 65,0
6 31,4 47,6 52,1 51,2 48,7 45,8 42,8 40,0 37,4 35,0
фг 5 43,6 69,1 80,0 84,0 85,6 86,3 86,7 86,8 87,0 87,0
6 41,4 61,6 66,0 63,6 59,4 54,9 50,6 46,8 43,3 40,2
fz z 5 18,3 42,5 75,6 118 170 232 302 383 473 572
6 18,7 41,8 73,2 112 150 211 268 330 394 462
f ФХ фХ 5 13,4 30,0 53,4 83,4 120 164 214 270 334 404
6 13,3 29,8 52,5 81,3 116 156 200 250 303 359
Значения частот /, /у и /фа зависят от диаметра Б аналогичным образом. Значения частоты /х при невесомых стержнях
(рис. 3, график 1) асимптотически приближаются к максимальному значению 65,1 Гц, полученному при с = 1000 Н-м для
модели № 1 с абсолютно жёсткими звеньями. Вид кривой /х(Ф) для невесомых стержней также, как и рис. 2, объясняется сложением податливостей привода и стержней.
Для собственных частот / /х, /фу с собственными формами, на которых выходное звено колеблется в вертикальной плоскости, зависимости от диаметра Ф при невесомых стержнях существенно отличаются, так как эти колебания не зависят от жёсткости привода (рис. 4, график 1). При увеличении диаметра стержня в два раза, с 12 мм до 24 мм, коэффициенты изгибной и крутильной жёсткости стержней увеличились по формулам (7) и (8) в 16 раз, а собственная частота / увеличилась в четыре раза с 42,5 Гц до 170 Гц, что хорошо согласуется с теорией колебаний. Кривая 1 на рис. 4 близка к параболе.
Введение в расчётную схему масс звеньев приводит к снижению значений собственных частот. Вид зависимостей /х(Ф) и /Ф) с учётом масс звеньев (график 2, рис. 3, 4) можно объяснить тем, что при малых значениях диаметра Ф масса выходного звена для модели № 6 существенно превышает суммарную массу стержней и определяет массу динамической системы, как в модели № 5 (график 1, рис. 3, 4). Например, при Ф = 20 мм суммарная масса шести стержней составляет 0,76 кг, а масса выходного звена - 1 кг. При больших значениях диаметра Ф масса системы определяется уже суммарной массой стержней, которая пропорциональна квадрату диаметра стержня, поэтому графики 1 и 2 на рис. 3, 4 расходятся с ростом Ф. Начиная с Ф = 24 мм, график 2 на рис. 3 близок к гиперболе, а график 2 на рис. 4 - к прямой.
Рис. 3. Зависимость fx(D) для моделей № 5 и № 6:
1 - без учёта масс стержней; 2 - с учётом масс Fig. 3. The dependence of fx(D) for models No. 5 and No. 6: 1 - without taking into account the masses of the rods; 2 - taking into account the masses
fi, Гц
500
400
300
200
100
1 ^
^ 2
12
16
20
24
28
32
36 40 D. мм
Рис. 4. Зависимость fz(D) для моделей № 5 и № 6:
1 - без учёта масс стержней; 2 - с учётом масс Fig. 4. Dependence fZ(D) for models No. 5 and No. 6: 1 - without taking into account the masses of the rods; 2 - taking into account the masses
Для моделей № 7 и № 8 с длинными стальными стержнями результаты вычислений приведены в табл. 3.
Таблица 3
Значения низших собственных частот для моделей № 7 и № 8, Гц
Table 3
Values of the lowest natural frequencies for models No. 7 and No. 8, Hz
f Собственная форма Номер модели № Диаметр стержня D, мм
20 30 40 50 60 70 80 90 100
f* x 7 7,68 11,8 13,5 14,1 14,4 14,5 14,5 14,6 14,6
8 2,76 2,90 2,49 2,08 1,76 1,53 1,34 1,19 1,08
ftyz 9z 7 45,9 70,8 80,7 84,3 85,7 86,3 86,6 86,8 86,9
8 3,26 3,28 2,77 2,30 1,95 1,68 1,47 1,31 1,18
fz z 7 4,59 10,3 18,4 28,7 41,3 56,2 73,4 92,9 115
8 2,86 4,8 6,70 8,56 10,4 12,2 14,0 15,8 17,6
fФ* ф* 7 14,5 32,6 58,0 90,6 131 178 232 294 362
8 3,46 5,27 7,06 8,85 10,6 12,4 14,2 16,0 17,7
Результаты расчёта хорошо согласуются с данными для моделей № 5 и № 6. Например, для модели № 7 при увеличении диаметра Б стержня в 2 раза (с 20 мм до 40 мм) собственная частота / с собственной формой г увеличилась с 4,59 Гц до 18,4 Гц - в 4 раза, а коэффициенты из-гибной и крутильной жёсткости, согласно формулам (7), (8), увеличились в 16 раз. Так как масса системы не меняется, то собственная частота на собственной форме
20 30 40 50 60 70 80 90 А мм
Рис. 5. Зависимость fX(D) для моделей № 7 и № 8: 1 - без учёта масс стержней; 2 - с учётом масс Fig. 5. The dependence of fX(D) for models No. 7 and No. 8: 1 - without taking into account the masses of the rods; 2 - taking into account the masses
г должна увеличиться примерно в 4 раза, что и наблюдается.
В модели № 8 в отличие от модели № 6 суммарная масса стержней превышает массу выходного звена уже при небольших значениях диаметра, так как стержни стальные и длинные. Поэтому графики 1 и 2 (рис. 5, 6) с ростом диаметра расходятся сразу, начиная с Б = 20 мм. Кривая 2 на рис. 6 - прямая, а кривая 2 на рис. 5, начиная с Б = 50 мм, близка к гиперболе.
20 30 40 50 60 70 80 90 О, мм
Рис. 6. Зависимость fz(D) для моделей № 7 и № 8:
1 - без учёта масс стержней; 2 - с учётом масс Fig. 6. Dependence fZ(D) for models No. 7 and No. 8: 1 - without taking into account the masses of the rods; 2 - taking into account the masses
Заключение
Построена аналитическая модель плоского манипулятора параллельной структуры с тремя степенями свободы, а также его компьютерные модели в программном комплексе «Универсальный механизм» для расчёта собственных частот и форм малых колебаний.
Определены три собственные частоты /х, /у, /щ колебаний манипулятора с абсолютно жёсткими невесомыми промежуточными звеньями в горизонтальной плоскости движения механизма аналитическим и численным методом. Для рассмотренной конфигурации плоского манипулятора поступательное и вращательное движения выходного звена на собственных формах разделены и не зависят друг от друга. Допустимая разность аналитического и численного решений подтверждает адекватность принятых расчётных схем и алгоритмов.
Замена абсолютно жёстких промежуточных звеньев упругими стержнями привела к появлению собственных частот /, /фх, /фу и соответствующих им форм, на которых выходное звено колеблется в вертикальной плоскости, перпендикулярной плоскости движения манипулятора. Поэтому, для анализа собственных колебаний плоского манипулятора недостаточно использовать плоскую кинематическую схему с абсолютно жёсткими звеньями. Установлено, что колебания выходного звена в вертикальной и горизонтальной плоскостях на собственных формах разделены. Учёт инерционных свойств упругих промежуточных звеньев привёл к появлению собственных частот и форм в нижней части спектра, на которых выходное звено совершает те же движения, что и на более низких частотах, но с меньшей амплитудой.
В результате вычислительных экспериментов с использованием разработанных компьютерных моделей определены зависимости значений собственных частот, соответствующих шести первым неизменным формам колебаний, от коэффициента жёсткости привода и от диаметра D стержней круглого поперечного сечения. Установлено, что для манипулятора с упругими стержнями значения частот fx, f, fz зависят от соотношения податливостей привода и стержней, образующих последовательное соединение, и с ростом коэффициента жёсткости привода не увеличиваются до бесконечности, как при абсолютно жёстких звеньях, а асимптотически приближаются к постоянным значениям парциальных частот упругих колебаний стержней. Поэтому введение упругих муфт в привод может существенно снизить собственные частоты колебаний манипулятора в горизонтальной плоскости. Значения собственных частот f, fx,fqy не зависят от коэффициента жёсткости
СПИСОК ИСТОЧНИКОВ
1. Merlet J.P. Parallel Robots. 2nd edition. Dordrecht: Springer, 2006. 417 p.
2. Ганиев Р.Ф., Глазунов В.А. Манипуляционные механизмы параллельной структуры и их приложения в современной технике. Доклады Академии наук. 2014;459(4):428. DOI 10.7868/S086956521434009X.
3. Маслов А.Н. Позиционирование нежесткого звена робота-манипулятора с учетом ограничений на управление. Вестник Московского энергетического института. Вестник МЭИ. 2011;2:5-9.
4. Антонов А.В., Глазунов В.А. Влияние упругих сил на точность движения манипулятора параллельной структуры. Экстремальная робототехника. 202;1(1):47-55.
5. Хейло С.В., Ширинкин М.А., Глазунов В.А. Определение собственных частот колебаний манипулятора параллельной структуры. Известия высших учебных заведений. Технология текстильной промышленности. 2011;4(333):120-124.
6. Носова Н.Ю. Разработка и исследование пространственных механизмов параллельной структуры с шарнирными параллелограммами с различным числом степеней свободы: специальность 05.02.18 «Теория механизмов и машин»: дис. на соискание учёной степени канд. техн. наук / Носова Наталья Юрьевна; Институт машиноведения им. А.А. Благонравова Российской академии наук. Москва, 2021. 152 с. Библиогр.: с. 119-130.
7. Демидов С.М., Артеменко Ю.Н., Глазунов В.А. и др. Анализ динамических свойств механизмов па-
привода, определяются упругими и инерционными свойствами стержней и увеличиваются с ростом диаметра невесомых стержней по квадратичному закону.
Влияние масс упругих звеньев на значения собственных частот определяется соотношением инерционных параметров выходного звена и стержней и видом собственных форм. В случае, когда масса системы определяется массой стержней, то значения частот/X, /у, обратно пропорциональны диаметру Ф, а значения частот /, _/*, /фу - прямо пропорциональны.
В результате исследования спектра собственных частот, сравнения значений коэффициентов изгибной, крутильной и продольной жёсткости в подвижных соединениях, анализа зависимости значений собственных частот от диаметра стержней установлено, что продольная жёсткость стержней не оказывает влияния на значения собственных частот в нижней части спектра.
раллельной структуры. Машиностроение и инженерное образование. 2012;1(30):36-41.
8. Скворцов С.А., Лысогорский А.Е., Глазунов В.А. Динамический анализ механизма параллельной структуры, выполняющего поступательные перемещения. Известия Юго-Западного государственного университета. Серия: Техника и технологии. 2015;2(15):70-79.
9. Глазунов В.А., Хейло С.В., Костюков А.М. Исследование колебаний механизма параллельной структуры. Вибрационные технологии, мехатро-ника и управляемые машины: сб. тр. XII междунар. конф.: в 2 частях. Курск, 18-20 мая 2016 года. Курск: Юго-Западный государственный университет. 2016. С. 23-28.
10. Pogorelov D. Y. On numerical methods of modeling large multibody systems. Mechanism and Machine Theory. 1999;34(5):791-800. DOI 10.1016/S0094-114X(98)00055-X.
11. Погорелов Д.Ю., Толстошеев А.К., Ковалёв Р.В. Динамический анализ и синтез механизмов с использованием программы UM. Брянск: Изд-во Брянского государственного технического университета, 1997. 46 с.
12. Универсальный механизм 9. Руководство пользователя: начинаем работать. 2021. 75 с. URL: http://www.universalmechanism.com/download/90/ru s/gs_um.pdf (дата обращения: 30.05.2022).
13. Толстошеев А.К., Татаринцев В.А. Структурный анализ механизмов роботов-станков с параллельной кинематикой. Вестник Брянского государ-
ственного технического университета. 2017;1(54):33-43. DOI 10.12737/24889. 14. Ширинкин М.А. Структурный анализ пространственных механизмов параллельной структуры с
REFERENCES
четырьмя и шестью степенями свободы. Машиностроение и инженерное образование. 2011;2:17-21.
1. Merlet JP. Parallel Robots. 2nd ed. Dordrecht: Springer; 2006.
2. Ganiev RF, Glazunov VA. Manipulation mechanisms of parallel structure and their applications in modern technology. Dokladi Akademii Nauk. 2014;459(4):428. DOI 10.7868/S086956521434009X.
3. Maslov AN. Positioning of non-rigid link of the robot manipulator taking into account control restrictions. Vestnik Moscovskogo Energeticheskogo Instituta. Vestnik MEI. 2011;2:5-9.
4. Antonov AV, Glazunov VA. Influence of elastic forces on the movement accuracy of the parallel manipulator. Extreme Robotics. 202;1(1):47-55.
5. Heylo SV, Shirinkin MA, Glazunov VA. Determination of natural oscillation frequencies of a parallel manipulator. Proceedings of Higher Educational Institutions. Textile Industry Technology. 2011;4(333):120-124.
6. Nosova NYu. Development and study of spatial mechanisms of a parallel structure with articulated parallelograms with various number of degrees of freedom [dissertation]. [Moscow (RF)] Mechanical Engineering Research Institute of the Russian Academy of Sciences; 2021.
7. Demidov SM, Artemenko YuN, Glazunov VA. Analysis of dynamic properties of parallel mechanisms. Mashinostroenie I Inzhenernoe Obrazovanie. 2012;1(30):36-41.
8. Skvortsov SA, Lysogorsky AE, Glazunov VA. Dynamical analysis of parallel structure mechanism
Информация об авторе:
Толстошеев Андрей Константинович - кандидат
технических наук, доцент, доцент кафедры «Трубопроводные транспортные системы» Брянского
executing translational motions. Proceedings of the Southwestern State University. Series: Engineering and Technologies. 2015;2(15):70-79.
9. Glazunov VA, Heylo SV, Kostyukov AM. Study of oscillations of parallel structure mechanism. Collection of Scientific Papers of the XII International Conference, May 18-20, 2016: Vibration Technologies, Mechatronics and Controlled Machines; Southwest State University. Kursk; 2016. p. 23-28.
10. Pogorelov DY. On numerical methods of modeling large multibody systems. Mechanism and Machine Theory. 1999;34(5):791-800. DOI 10.1016/S0094-114X(98)00055-X.
11. Pogorelov DYu, Tolstosheev AK, Kovalev RV. Dynamic analysis and synthesis of mechanisms using UM programs. Bryansk: Publishing House of Bryansk State Technical University; 1997.
12. Universal mechanism 9. User Manual: Getting started [Internet]; 2021 [cited 2022 May 30]. Available from: http://www.universalmechanism.com/download/90/ rus/gs_um.pdf
13. Tolstosheev AK, Tatarintsev VA. Structural analysis of robot machine mechanisms with parallel kinematics. Bulletin of Bryansk State Technical University. 2017;1(54):33-43. DOI 10.12737/24889.
14. Shirinkin MA. Structural analysis of spatial mechanisms of a parallel structure with four and six degrees of freedom. Mashinostroenie I Inzhenernoe Obrazovanie. 2011;2:17-21.
государственного технического университета, тел. В-953-293-04-51,
Tolstosheev Andrey Konstantinovich - Candidate of Technical Sciences, Associate Professor of the De-
Статья опубликована в режиме Open Access. Article published in Open Access mode.
partment of Pipeline Transport Systems at Bryansk State Technical University, phone: 8-953-293-04-51.
Статья поступила в редакцию 06.07.2022; одобрена после рецензирования 23.09.2022; принята к публикации 21.11.2022. Рецензент - Куц В.В., доктор технических наук, профессор кафедры «Машиностроительные технологии и оборудование» Юго-Западного государственного университета, член редсовета журнала «Транспортное машиностроение»
The article was submitted to the editorial office on 06.07.2022; approved after review on 23.09.2022; accepted for publication on 21.11.2022. The reviewer is Kuts V.V, Doctor of Technical Sciences, Professor, Head of the Department of Engineering Technologies and Equipment at the Southwest State University, member of the Editorial Board of the journal Transport Engineering.