УДК 531.8; 621.01 doi: 10.18698/0536-1044-2020-8-21-28
Влияние особых положений
на жесткость механизмов параллельной структуры
К.Г. Эрастова
МГТУ им. Н.Э. Баумана
The Influence of Special Positions on the Rigidity of Parallel Structure Mechanisms
K.G. Erastova
Bauman Moscow State Technical University
Применение роботов на основе механизмов с замкнутыми кинематическими цепями становится перспективным вследствие резкого повышения точности изготовления и контроля изделий машиностроения. Одним из преимуществ механизмов параллельной структуры над таковыми с последовательной структурой является повышенная жесткость конструкции. Отсутствие закономерности изменения жесткости по рабочей зоне и влияние так называемых особых положений обусловливают необходимость разработки новых методов оптимизационного синтеза механизма, основанных на этом параметре. Предложен подход, базирующийся на использовании матрицы Якоби механизма, позволяющий выявлять жесткость в наихудшем направлении приложения внешней силы в каждой точке рабочего пространства. Такой подход позволяет определять форму и размеры эффективной рабочей зоны. На примере плоского пятизвенного механизма продемонстрировано применение этого подхода: получены матрицы жесткости звеньев, построен параллелограмм сил и определено распределение жесткости по рабочей зоне. Показано, что при приближении к особому положению жесткость механизма уменьшается в десятки раз, а непосредственно в нем равна нулю. Предложенный подход может быть использован при оптимизации геометрических параметров механизма из условия максимизации эффективной рабочей зоны.
Ключевые слова: механизмы параллельной структуры, жесткость механизма, матрица Якоби, оптимизационный синтез, робототехнические системы, рабочая зона
Due to a sharp increase in the accuracy of manufacturing and quality control of mechanical engineering products, robots based on mechanisms with closed-loop kinematic chains have become a promising area of engineering. One of the advantages of parallel structure mechanisms over serial mechanisms is their increased rigidity. The absence of patterns in the change of rigidity in the working zone and the influence of the so-called special positions necessitate the development of new methods of optimization synthesis of the mechanism. This paper presents an approach based on the Jacobian matrix using which it is possible to determine rigidity in the worst direction of the external force application in each point of the working zone. The same approach can also be used for determining the shape and size of the effective working zone. An example of a flat five-chain mechanism is used to demonstrate the application of the method and obtain a matrix of chain rigidity, construct the parallelogram of forces and determine rigidity distribution in the working zone. It is shown that when approaching a special position, the rigidity of the mechanism decreases ten-fold, and in the special position it equals a zero. The proposed method can be used for optimizing the geometrical parameters of the mechanism on the conditions of maximising the effective working zone.
Keywords: mechanisms of parallel structure, mechanism rigidity, Jacobian matrix, optimization synthesis, robotic systems, working zone
Механизмы параллельной структуры (МПС) нашли широкое применение в различных областях техники. Их используют как в любительских 3D-принтерах, так и в высокоскоростных манипуляционных комплексах, медицинских роботах, и в установках, работающих в агрессивных средах.
Повышенный интерес к механизмам такого типа обусловлен их характеристиками, пониженной металлоемкостью и повышенными грузоподъемностью и жесткостью по сравнению с механизмами последовательной структуры. МПС воспринимают нагрузку подобно пространственным фермам, что позволяет им манипулировать более массивными объектами [1], а жесткость замкнутой кинематической цепи гораздо выше незамкнутой [2].
Несмотря на перечисленные достоинства, МПС обладают существенным недостатком: в некоторых положениях, называемых особыми (или сингулярными), кинематическая структура механизма вырождается, и выходное звено (или отдельные звенья) манипулятора могут потерять управляемость или подвижность [2, 3].
Негативные явления (снижение точности позиционирования и жесткости, повышение нагрузки на приводы) начинают проявляться уже при приближении к особым положениям [4], что обусловливает необходимость тщательного синтеза и анализа робота-манипулятора и подбора геометрических параметров, минимизирующих размеры данных областей.
Цель работы — исследование влияния особых положений на жесткость МПС, а также создание алгоритмов вычисления формы и размера эффективной рабочей зоны манипулятора по жесткости. Это позволит сформировать практические рекомендации для конструкторов робототехнических систем.
Эффективная рабочая зона манипулятора определяется как множество точек, в которых может находиться выходное звено МПС согласно конструктивным ограничениям реального механизма и требованиям по нужной ориентации выходного звена при условии сохранения в этих точках необходимых конструктивных параметров [5].
В качестве критического параметра будем рассматривать жесткость, которая является
важнейшей характеристикой для различных типов машин и роботов. Однако отсутствие методики оптимизационного синтеза механизмов из условия достижения наилучших жесткостных характеристик при заданных конструктивных и эксплуатационных условиях обусловливает актуальность данной работы. Дополнительно осложняет задачу синтеза отсутствие закономерности изменения жесткости — изоморфности жесткости по рабочей зоне [6].
Рассмотрим некую обобщенную модель МПС, имеющего п степеней свободы и п кинематических цепей, каждая из которых состоит из звеньев, соединенных одноподвижными шарнирами. Каждая кинематическая цепь содержит только одну приводную пару. Известно [7], что при исследовании только кинематики механизма можно связать усилия в приводных парах и силу на выходном звене, генерируемую механизмом, через следующее соотношение:
J-T T = -P,
(1)
где | т — обратная транспонированная матрица Якоби механизма; Т = (Ть..., Тп )т — массив усилий в приводных парах каждой цепи; Р = (Рх1,..., Рхп )т — компонента внешней нагрузки, соответствующая выходной координате Хк (к = 1, ..., п).
Обратная транспонированная матрица Якоби механизма определяется как
(
r-T -
dF1/dx1 dFJdqi
dF1/dx„ dFJdqi
dF„ /dx1 Л dF„ /dq„
dF„ /dx„
dF„ /dq
„
где Рг (х1,..., хп, qi) = 0 — уравнения связи механизма; qi — входные координаты.
Максимальный момент в приводе МПС достигается в случае, когда г-я строка }т матрицы ]т и вектор внешней силы Р являются кол-линеарными:
Т шах | = }ТР. (2)
Известно, что угловые и линейные перемещения выходного звена связаны с вектором
действующей на выходное звено силои с помощью матрицы жесткости К:
М = Р,
(3)
где d = (йА,..., йх„)т — упругая деформация цепи, соответствующая координате Хк.
Объединяя выражения (1) и (3), получаем
-I -т Т = Kd. (4)
Матрицы I и К зависят только от кинематической структуры, геометрических параметров и текущего положения выходного звена МПС. Таким образом, из выражений (2) и (4) видно, что направление вектора внешней силы, вызывающее максимальный момент в г-м приводе, совпадает с направлением внешней силы, обусловливающим максимальную упругую деформацию г-й цепи:
хк г г
= к-1 ^ ,к1 Р = кк ~Р, К
(5)
где к к — к-я строка матрицы жесткости К.
Матрицу жесткости механизма можно составить, зная матрицы жесткостей отдельных звеньев в глобальной системе координат (ГСК) [8]: К =
( К1 1Ч-11
К 21
К 22 + К 21
о
К 22
о ^
о
о
0 0 0 .. 0 К -н К12 0
0 0 0 К >-1 • • 21 К 2-1 + К11 К12
0 0 0 .. 0 К 21 к 22
(6)
Здесь К' — в ГСК,
К > =
ч к 21
к 12 ^
к 22,
(7)
К > = Б тК ЦЬ,
Ь =
V0 0 ■■• Иу
г 2С
с
В\ ^^ /С |
X \ 1 \ \ ФИ
\ \>-—-О ^ >. \ \>_ '1- /\Ф2
^■у//)//. \ у/////.
А | | Л2
матрица жесткости ]-го звена
хс х
Рис. 1. Кинематическая схема пятизвенного механизма: А{ — приводные пары (г — номер цепи); — пассивные (неприводные) пары; С — выходное звено; ф,- — углы поворота входных звеньев (входные координаты); хс, хс — координаты точки С
Таким образом, подход позволяет при заданных геометрических параметрах с помощью итерационного метода [5, 9-11] оценить смещение выходного звена (и промежуточных шарниров) под действием внешней нагрузки в наихудшем направлении. Возможно решение обратной задачи: поиск оптимальных геометрических параметров манипулятора из условия непревышения максимально допустимого смещения.
Рассмотрим предложенный подход на примере плоского пятизвенного механизма. Схема механизма показана на рис. 1, а его 3D-визуа-лизация — на рис. 2. Физические и геометрические параметры механизма (см. таблицу) подобраны заведомо неоптимальным образом для наилучшей визуализации описанного подхода. Уравнения связи и анализ рабочей зоны такого механизма приведены в статьях [12, 13]. Для расчета в первом приближении вращательные шарниры примем абсолютно жесткими.
Геометрические и физические параметры механизма для расчета
где Ь — матрица перехода от локальной системы координат (ЛСК) звена в ГСК; К \оссй — матрица жесткости _/-го звена в ЛСК.
Матрица перехода от ЛСК звена в ГСК имеет вид
(И 0 ■.. 0^ 0 И ■.. 0
Параметр Значение
Расстояние между входными 100 парами /А[А2, мм
Длина звена /, мм 175
Момент инерции сечения относительно оси у 1у, мм4 3,41-103
Момент инерции сечения относительно оси х 1х, мм4 1,33-103
Полярный момент инерции сечения 1р, мм4 3,26-103
Модуль Юнга Е, МПа 68,9-103
Модуль сдвига С, МПа 25,9-103
где И — матрица поворота.
Коэффициент Пуассона ц
0,33
Рис. 2. 3Б-визуализация манипулятора на основе пятизвенного механизма: хгугхг — локальные системы координаты звеньев
Для описания каждого звена будем использовать пространственный балочный элемент. Матрица жесткости в соответствующей ЛСК имеет вид [8]
E
K local = J3 M,
где М — матрица размером 12x12, m11 = = Ш7.7 = Al2, mi,7 = Ш7,1 = -Al2, m2,2 = m.8,8 = = 12IZ, m2,6 = m2,i2 = m6,2 = mi2,2 = 6lIz, m2,8 = = m8,2 = -12IZ, тз,з = m9,9 = 12Ir, m3,5 = m3,n = = m5,3 = mn,3 =-6lIy, тз,9 = т9,з = -12Iy,
m4,4 = тю,ю = GIpl2 IE, mi>w = mw,4 =-GIpl 2IE, m5,5 = mn,n = 4l2 Iy, m5,9 = m8,n = m9,5 = mn,8 = = 6lIy, m5,n = mn,5 = 2l2 Iy, m6,6 = m12,12 = 4l2 Iz, m6,8 = m8,6 = m8,12 = m12,8 = -6lIz, m6,12 = mn,6 = = 2l2 Iz, остальные элементы матрицы М — нулевые; А — площадь сечения балочного элемента.
Эта матрица жесткости записана для балочного элемента с заделкой на обоих концах. В рассматриваемой задаче необходимо составить матрицы жесткости для двух типов балок: заделка — одноподвижный шарнир (первое звено в каждой цепи, привод условно принят абсолютно жестким) и одноподвижный шарнир — одноподвижный шарнир (второе звено в каждой цепи). Формирование новых матриц жесткости проводят следующим образом:
• матричное уравнение (3) записывают в виде системы линейных уравнений;
• силовой фактор Рг, соответствующий подвижности шарнира, приравнивают к нулю;
• из ¿-го уравнения выражают ¿-е перемещение, которое подставляют в остальные уравнения системы;
• из новой системы линейных уравнений формируют уравнение в матричной форме и, соответственно, матрицу жесткости.
В соответствии с рис. 1 и 2 матрицы жесткости для первого и второго звена каждой цепи в ЛСК записываются следующим образом:
K 1ocal =— N,
E l3
где N — матрица размером 12x12, nu =
= П7,7 = Al2, П1,7 = «7,1 = - Al2 , П2,2 = П8,8 = 12IZ ,
«2,6 = «2,12 = «6,2 = «12,2 = 6lIz, «2,8 = «8,2 =-12Iz,
«3,3 = «9,9 = 3Iy, «3,5 = «5,3 = -3lIy, «3,9 = «9,3 = = -3Iy, «4,4 = «10,10 = GIpl2 IE, «4,10 = «10,4 =
= -GIp l2 IE,
«6,6 = 4l2 Iz,
«5,5 = 3l2 Iy , «5,9 = «9,5 = 3lIy , «6,8 = «8,6 = «8,12 = «12,8 =-6lIz ,
«6,12 = «12,6 = 2l2 Iz, «12,12 = 4l2 Iy, остальные элементы матрицы N — нулевые;
E l3
Klocal = H,
где H — матрица размером 12x12, h1,1 = = h7,7 = Al2, Ни = h7,1 = - Al2 , h2,2 = h8,8 = 12Iz, h2,6 = h2,12 = h6,2 = h12,2 = 6lIz, h2,8 = Н8,2 = -12Iz, h4,4 = Ню,ю = GIpl2 IE, h4,10 = hw,4 =-GIpl 2IE,
^6,6 = 412!х , Нб,8 = ^8,6 = ^8,12 = ^12,8 = 6 ,
йб,12 = Й12,6 = 212/г, ^12,12 = 4/21у, остальные элементы матрицы Н — нулевые.
Блок-схема алгоритма расчета жесткости механизма по рабочей зоне приведена на рис. 3. Согласно алгоритму, в каждой точке рабочей зоны матрицы жесткости записываются в ГСК в соответствии с выражениями (7), где матрица поворота имеет вид
(
И
V
со8 а 0
- в1п а
81и а 0
со8 а
\
/
Здесь для первого и четвертого звеньев а = ф, (см. рис. 1, 2), для второго и третьего звеньев
а = агс^2(гс-г,ц, хс -хв,),
где гв,, хв, — координаты точек Вг.
После проведения ансамблирования согласно выражению (6) из матрицы жесткости удаляются пустые строки и столбцы, а также первые и последние шесть строк и столбцов, соответствующие заделкам. Таким образом, формируется окончательная матрица жесткости системы размером 15x15.
Внешняя сила прикладывается в плоскости хОх в наихудшем направлении согласно выражению (5). При этом разработанная трехмерная модель позволит в дальнейшем проводить математическое моделирование, включающее в себя всевозможные силовые факторы.
Для визуализации этого этапа расчета построим в некоторой точке рабочего пространства параллелограмм сил (рис. 4) [14], показывающий, какую максимальную силу может сгенерировать выходное звено манипулятора в каждом из направлений при известном значении максимально допустимого момента в приводе Тшах (в данной задаче Тшах = 2000 Н-мм).
Наихудшее направление внешней силы определяется меньшей диагональю параллелограмма и соответствует направлению к приближению к особому положению. На рис. 5 показаны наихудшие направления внешней силы для обоих приводов. Синим выделена зона повышенной нагрузки на первый привод, красным — на второй. Синим пунктиром обозначено примерное расположение зоны особых положений типа ПУ [3, 15], в которых выходное звено механизма теряет управляемость.
В каждой точке рабочей зоны вычислено перемещение в направлении действующей внеш-
Вывод массивов
С
Конец
3
Рис. 3. Блок-схема алгоритма расчета жесткости
механизма по рабочей зоне: хшах, ушах и хшп, уш1п — максимальные и минимальные значения хи у
4,15 Н
28,05 Н
Рис. 4. Параллелограмм сил пятизвенного манипулятора в точке [50, 150]
z, м
-300 -200 -100
Рис. 5. Наихудшие направления внешней силы, приложенной к выходному звену, для каждой точки рабочей зоны
ней силы и соответствующая жесткость. Для наглядности результат визуализирован в виде графика, приведенного на рис. 6, демонстрирующего распределение жесткости по рабочей зоне в плоскости механизма.
Как видно из рис. 6, средняя жесткость по рабочей зоне составляет около 100 Н/мм, в то время как при приближении к особым положениям она резко падает до нуля. Таким образом, работа выходного органа вблизи особых положений невозможна, и эти области не могут быть включены в эффективную рабочую зону механизма.
Выводы
1. Предложен подход, основанный на использовании матрицы Якоби механизма, позволяющий определить наихудшее значение
0 0,5 1,0 1,5 2,0 2,5 3,0 Жесткость (абсолютное значение, Н/мм), 1о§10
Рис. 6. Распределение жесткости по рабочей зоне в плоскости механизма
жесткости механизма в каждой точке рабочего пространства и таким образом определить размер и форму эффективной рабочей зоны манипулятора.
2. Разработанный подход решает актуальную задачу исследования неизоморфных жесткост-ных характеристик манипуляторов параллельной структуры. Применение подхода подробно продемонстрировано на примере плоского пя-тизвенного механизма.
3. В дальнейших исследованиях планируется изучение влияния жесткости приводов и шарниров, проведение натурного эксперимента и разработка автоматизированных процессов оптимизации.
Литература
[1] Крайнев А.Ф. Механика машин. Фундаментальный словарь. Москва, Машинострое-
ние, 2000. 904 с.
[2] Merlet J.-P. Parallel Robots. Springer, 2006. 394 p.
[3] Ларюшкин П.А. Оценка близости к особым положениям механизмов параллельной
структуры путем дифференцирования уравнений связи. Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение, 2019, № 1, с. 71-83, doi: 10.18698/0236-3941-2019-1-71-83
[4] Arakelian V., Briot S., Glazunov V. Increase of singularity-free zones in the workspace of
parallel manipulators using mechanisms of variable structure. Mechanism and Machine Theory, 2008, vol. 43, no. 9, pp. 1129-1140, doi: 10.1016/j.mechmachtheory.2007.09.005
[5] Эрастова К.Г., Ларюшкин П.А. Рабочие зоны механизмов параллельной структуры и
способы определения их формы и размеров. Известия высших учебных заведений. Машиностроение, 2017, № 8, с. 78-87, doi: 10.18698/0536-1044-2017-8-78-87
[6] Глазунов В.А., Есина М.Г., Быков Р.Э. Управление механизмами параллельной струк-
туры при переходе через особые положения. Проблемы машиностроения и надежности машин, 2004, № 2, с. 78-84.
[7] Ларюшкин П.А., Глазунов В.А., Эрастова К.Г. Определение максимальных усилий в
приводах манипуляторов параллельной структуры по заданной величине внешней нагрузки. Машиностроение и инженерное образование, 2016, № 2(47), с. 40-46.
[8] Kassimali A. Matrix A«alysis of Structures. Cengage Learning, 2012. 657 p.
[9] Ларюшкин П.А., Палочкин С.В. Рабочая зона манипулятора параллельной структуры с
тремя степенями свободы. Известия высших учебных заведений. Технология текстильной промышленности, 2012, № 3(339), с. 92-96.
[10] Пащенко В.Н. Построение рабочей зоны шестистепенного манипулятора параллельной структуры на базе кривошипно-шатунного механизма. Интернет-журнал «Науковедение», 2016, т. 8, № 3, doi: 10.15862H42TVN316. URL: http://naukovedenie.ru/PDF/ 142TVN316.pdf (дата обращения 28 января 2020).
[11] Merlet J.-P. Workspace-oriented methodology for designing a parallel manipulator. Pro-ceedi«gs — IEEE I«ter«atio«al Co«fere«ce o« Robotics a«d Automatio«, 1996, vol. 4, pp. 3726-3731.
[12] Захаров М.Н., Ларюшкин П.А., Эрастова К.Г. Определение геометрических параметров плоского пятизвенного механизма из условия обеспечения его устойчивой работы в заданной области. Вестник машиностроения, 2017, № 11, с. 10-13.
[13] Эрастова К.Г., Ларюшкин П.А. Оптимизация геометрических параметров механизмов параллельной структуры методом SQP из условия максимизации эффективной рабочей зоны на примере пятизвенного манипулятора. Машиноведение и инновации. Тр. XXIX междунар. конф, Москва, 06-08 декабря 2017, Москва, ИМАШ РАН, 2018, с. 356-359.
[14] Lynch K.M., Park F.C. Moder« Robotics: Mecha«ics, Pla««i«g, a«d Co«trol. Cambridge University Press, 2017. 624 p.
[15] Глазунов В.А., Нгуен Нгок Хуэ, Нгуен Минь Тхань. К анализу особых положений механизмов параллельной структуры. Машиностроение и инженерное образование, 2009, № 4, с. 11-16.
References
[1] Kraynev A.F. Mekha«ika mashi«. Fu«dame«tal'«yy slovar' [Mechanics of cars. Fundamental
Dictionary]. Moscow, Mashinostroyeniye publ., 2000. 904 s.
[2] Merlet J.-P. Parallel Robots. Springer, 2006. 394 p.
[3] Laryushkin P.A. Estimating How dose Parallel Mechanisms are to their Singularities by Dif-
ferentiating Constraint Equations. Herald of the Bauma« Moscow State Tech«ical U«iversi-ty. Series Mecha«ical E«gi«eeri«g, 2019, no. 1, pp. 71-83 (in Russ.), doi: 10.18698Ю236-3941-2019-1-71-83
[4] Arakelian V., Briot S., Glazunov V. Increase of singularity-free zones in the workspace of
parallel manipulators using mechanisms of variable structure. Mecha«ism a«d Machi«e Theory, 2008, vol. 43, no. 9, pp. 1129-1140, doi: 10.1016Ij.mechmachtheory.2007.09.005
[5] Erastova K.G., Laryushkin P.A. Workspaces of Parallel Mechanisms and Methods of Deter-
mining Their Shape and Size. Proceedi«gs of Higher Educatio«al I«stitutio«s. МаЛ« Buildi«g, 2017, no. 8, pp. 78-87 (in Russ.), doi: 10.18698Ю536-1044-2017-8-78-87
[6] Glazunov V.A., Esina M.G., Bykov R.E. Controlling mechanisms of a parallel structure when
passing through special positions. Journal of Machi«ery Ma«ufacture a«d Reliability, 2004, no. 2, pp. 78-84 (in Russ.).
[7] Laryushkin P.A., Glazunov V.A., Erastova K.G. Computation of the maximal actuation ef-
forts in parallel manipulators at a specified value of external load. Mashi«ostroye«iye i i«-zhe«er«oye obrazova«iye, 2016, no. 2(47), pp. 40-46 (in Russ.).
[8] Kassimali A. Matrix A«alysis of Structures. Cengage Learning, 2012. 657 p.
[9] Laryushkin P.A., Palochkin S.V. Workspace of the parallel structure manipulator with three
degrees of freedom. Izvestiya vysshikh ucheb«ykh zavede«iy. Tekh«ologiya tekstil'«oy promyshle««osti, 2012, no. 3(339), pp. 92-96 (in Russ.).
[10] Pashchenko V.N. The construction working area shestistrunnoi manipulator of the parallel structure on the basis of the crank mechanism. Internet-zhurnal "Naukovedeniye", 2016, vol. 8, no. 3, doi: 10.15862/142TVN316. Available at: http://naukovedenie.ru/PDF/ 142TVN316.pdf (accessed 28 January 2020).
[11] Merlet J.-P. Workspace-oriented methodology for designing a parallel manipulator. Proceedings — IEEE International Conference on Robotics and Automation, 1996, vol. 4, pp. 3726-3731.
[12] Zakharov M.N., Laryushkin P.A., Erastova K.G. Determination of geometric parameters of plain five-link mechanism from condition of assurance of its operation stability in specified area. Vestnik mashinostroyeniya, 2017, no. 11, pp. 10-13 (in Russ.).
[13] Erastova K.G., Laryushkin P.A. Optimization of the geometric parameters of parallel structure mechanisms using the SQP method from the condition of maximizing the effective working area using the example of a five-link manipulator. Mashinovedeniye i innovatsii. Tr. XXIX mezhdunar. konf. [Machine science and innovation. Proceedings of the XXIX International Conference]. Moscow, 2018, pp. 356-359.
[14] Lynch K.M., Park F.C. Modern Robotics: Mechanics, Planning, and Control. Cambridge University Press, 2017. 624 p.
[15] Glazunov V.A., Nguyen Ngok Khue, Nguyen Min' Tkhan'. Singular configuration analysis of the parallel mechanisms. Mashinostroyeniye i inzhenernoye obrazovaniye, 2009, no. 4, pp. 11-16 (in Russ.).
Информация об авторе
ЭРАСТОВА Ксения Георгиевна — ассистент кафедры «Основы конструирования машин». МГТУ им. Н.Э. Баумана (105005, Москва, Российская Федерация, 2-я Бауманская ул., д. 5, стр. 1, e-mail: erastovakg@gmail.com).
Статья поступила в редакцию 29.02.2020 Information about the author
ERASTOVA Ksenia Georgievna — Research Assistant, Scientific Research Institute of Automation of Production Processes. Bauman Moscow State Technical University (105005, Russian Federation, 2nd Baumanskaya St., Bldg. 5, Block 1, e-mail: erastovakg@gmail.com).
Просьба ссылаться на эту статью следующим образом:
Эрастова К.Г. Влияние особых положений на жесткость механизмов параллельной структуры. Известия высших учебных заведений. Машиностроение, 2020, № 8, с. 21-28, doi: 10.18698/0536-1044-2020-8-21-28
Please cite this article in English as: Erastova K.G. The Influence of Special Positions on the Rigidity of Parallel Structure Mechanisms. BMSTU Journal of Mechanical Engineering, 2020, no. 8, pp. 21-28, doi: 10.18698/0536-1044-2020-8-21-28