Научная статья на тему 'Программное управление движением робота с дифференциальным приводом при разных моделях трения'

Программное управление движением робота с дифференциальным приводом при разных моделях трения Текст научной статьи по специальности «Строительство и архитектура»

CC BY
15
0
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
программное управление роботом / управление в неголономной механике / модель трения Кулона / поликомпонентное трение / ошибка угловой скорости вращения / software control of a robot / control in nonholonomic mechanics / Coulomb friction model / multicomponent friction / error in angular velocity of rotation

Аннотация научной статьи по строительству и архитектуре, автор научной работы — Салимов Максим Сергеевич, Меркурьев Игорь Владимирович

Введение. Проектирование систем управления движением мобильных роботов требует построения математических моделей. Исследователи неоднократно обращались к этой теме. В частности, опубликованы работы, посвященные расчетам мультифизических процессов, моделированию движения различных видов колес в тех или иных условиях. Кроме того, рассматривались динамика деформируемых контактирующих тел при скольжении, качении и вращении, вопросы автономности и управляемости мобильных роботов. Отметим, что динамика и точность позиционирования колесных роботов во многом определяется трением. В литературе не представлены исследования динамики робота с дифференциальным приводом с учетом взаимосвязей эффектов трения скольжения, верчения и качения на основе теории поликомпонентного трения. Изыскания в этой области способны выявить новые динамические эффекты. Основываясь на полученных таким образом данных, можно заняться повышением точности позиционирования при построении математических моделей управления. Цель представленной работы — исследовать движение автоматического устройства с дифференциальным приводом с учетом трех моделей контакта: неголономного, кулоновского трения и поликомпонентного. Материалы и методы. В качестве базовой приняли схему двухколесного робота с дифференциальным приводом и безотрывным движением по опорной поверхности. Перемещение устройства обеспечивалось программным управлением. Динамику описали в виде уравнений Аппеля. Для вычислений использовали математические модели, по-разному учитывающие трение. Координация действий механизма формировалась на динамическом уровне. Управляющие воздействия — моменты двигателей колес. При визуализации исследуемых моделей задействовали встроенные численные методы системы Wolfram Mathematica («Вольфрам математика») с минимальной точностью 10 –6. Результаты исследования. При построении математической модели определили равенства для угловых скоростей колес. Учли наличие площадки контакта и вывели уравнения динамики робота с дифференциальным приводом. Элементами системы были проекции силы и момент, показатели верчения платформы, массы, угловые ускорения и инерция колес. Показали, как формируются управляющие воздействия в рамках неголономной механики. Описали модель двигателей, которые создают момент управления на ведущих колесах. Вывели решение как взаимосвязь индуктивности проводников электродвигателей и работы электропитания. Детально рассмотрели три модели, описывающие динамику робота с дифференциальным приводом. Первая — неголономная. Вторая и третья включали систему уравнений динамики робота с дифференциальным приводом для общего случая, при наличии площадки контакта. При этом во второй модели игнорируется время переключений в двигателе и задействуется трение Кулона. В третьей ввели параметр, определяющий скорость переходных процессов в двигателе, и задействовали разложения Паде. Это модель споликомпонентным трением. Итоги расчетов показали в виде графиков. На них исследуемые модели визуализировали в виде кривых разного цвета. Сопоставление графиков показало, в каких случаях после завершения переходных процессов управление обеспечит требуемую точность. Это модели 1 и 2. В модели 3 программное управление генерирует ошибку в угловой скорости вращения платформы. Данную ошибку нельзя спрогнозировать в рамках 1-й и 2-й модели. Во всех рассмотренных системахскорость скольжения колес в поперечном направлении падает до нуля. Получено и обосновано условие безотрывного движения опорного колеса. Обсуждение и заключение. Программное управление допустимо в моделях, не учитывающих трение колес при одновременном скольжении, верчении и качении (общий случай пространственного движения). Однако важно учитывать взаимосвязь данных процессов и поликомпонентного трения. Это необходимо для более точного выполнения роботом программных движений. Установлено, что программное управление в модели, учитывающей трение верчения и качения, приводит к отклонениям от программных значений угловой скорости платформы. Полученные результаты можно использовать при построении системы управления с прогнозирующими моделями.

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

Похожие темы научных работ по строительству и архитектуре , автор научной работы — Салимов Максим Сергеевич, Меркурьев Игорь Владимирович

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

Software Control of the Movement of a Differential Drive Robot for Different Friction Models

Introduction. Designing motion control systems for mobile robots requires the construction of mathematical models. Researchers have repeatedly addressed this topic. In particular, works have been published on the calculations of multiphysical processes, modeling the movement of various types of wheels under certain conditions. In addition, the dynamics of deformable contacting bodies during sliding, rolling and rotation, issues of autonomy and controllability of mobile robots were considered. Note, however, that the dynamics and positioning accuracy of wheeled robots is largely determined by friction. The literature does not present studies on the dynamics of a robot with a differential drive taking into account the interrelationships of sliding, spinning and rolling friction effects based on the theory of multicomponent friction. Research in this area can reveal new dynamic effects. Based on the data obtained in this way, it is possible to improve the accuracy of positioning in building mathematical control models. The presented work aims at investigating the movement of an automatic device with a differential drive taking into account three contact models: nonholonomic, Coulomb friction, and multicomponent models. Materials and Methods. The scheme of a two-wheeled robot with differential drive and continuous movement on the support surface was adopted as the basic one. The movement of the device was provided through software control. The dynamics was described in the form of Appel equations. Mathematical models were used for calculations, taking into account friction in different ways. Coordination of the actions of the mechanism was formed at a dynamic level. The control actions were the moments of the wheel motors. When visualizing the models under study, the built-in numerical methods of the Wolfram Mathematica system were used with a minimum accuracy of 10 -6 . Results. When building a mathematical model, the equations for the angular velocities of the wheels were determined. The authors took into account the presence of a contact site and derived the equations of dynamics of a differential drive robot. The elements of the system were force and moment projections, indicators of platform spin, masses, angular accelerations, and inertia of the wheels. It was shown how control actions were formed within the framework of nonholonomic mechanics. The model of engines that created a moment of control on the driving wheels was described. The solution was derived as the relationship between the inductance of the conductors of electric motors and the operation of the power supply. Three models describing the dynamics of a differential drive robot were examined in detail. The first model was nonholonomic. The second and third included a system of equations for the dynamics of a differential drive robot for a general case with a contact platform. At the same time, in the second model, the switching time in the engine was ignored and the Coulomb friction was involved. In the third model, a parameter to determine the speed of transients in the engine was introduced, and Pade decomposition was involved. This was a model with multicomponent friction. The calculation results were shown in the form of graphs. On them, the studied models were visualized in the form of curves of different colors. Comparison of the graphs showed in which cases, after the completion of transients, the control provided the required accuracy. These were models 1 and 2. In model 3, the software control generated an error in the angular velocity of rotation of the platform. This error could not be predicted within the framework of the 1st and 2nd models. In all the systems considered, the sliding speed of the wheels in the transverse direction dropped to zero. The condition of continuous motion of the support wheel was obtained and validated. Discussion and Conclusion. Software control is acceptable in models that do not take into account wheel friction during simultaneous sliding, spinning and rolling (general case of spatial motion). However, it is important to consider the relationship between these processes and multicomponent friction. This is required for the robot to perform program movements more accurately. It was established that software control in a model that takes into account the friction of spinning and rolling caused deviations from the program values of the angular velocity of the platform. The results obtained can be used in the building of a control system with predictive models.

Текст научной работы на тему «Программное управление движением робота с дифференциальным приводом при разных моделях трения»

МЕХАНИКА MECHANICS

0

H) Check for updates

В

a

о

T3

<D f

УДК 531.36 Научная статья

https://doi.org/10.23947/2687-1653-2023-23-4-346-355

Программное управление движением робота с дифференциальным приводом при разных моделях трения

М.С. Салимов И, И.В. Меркурьев

Национальный исследовательский университет «Московский энергетический институт», г. Москва, Российская Федерация И SalimovMS@mpei.ru

Аннотация

Введение. Проектирование систем управления движением мобильных роботов требует построения математических моделей. Исследователи неоднократно обращались к этой теме. В частности, опубликованы работы, посвященные расчетам мультифизических процессов, моделированию движения различных видов колес в тех или иных условиях. Кроме того, рассматривались динамика деформируемых контактирующих тел при скольжении, качении и вращении, вопросы автономности и управляемости мобильных роботов. Отметим, что динамика и точность позиционирования колесных роботов во многом определяется трением. В литературе не представлены исследования динамики робота с дифференциальным приводом с учетом взаимосвязей эффектов трения скольжения, верчения и качения на основе теории поликомпонентного трения. Изыскания в этой области способны выявить новые динамические эффекты. Основываясь на полученных таким образом данных, можно заняться повышением точности позиционирования при построении математических моделей управления. Цель представленной работы — исследовать движение автоматического устройства с дифференциальным приводом с учетом трех моделей контакта: неголономного, кулоновского трения и поликомпонентного. Материалы и методы. В качестве базовой приняли схему двухколесного робота с дифференциальным приводом и безотрывным движением по опорной поверхности. Перемещение устройства обеспечивалось программным управлением. Динамику описали в виде уравнений Аппеля. Для вычислений использовали математические модели, по-разному учитывающие трение. Координация действий механизма формировалась на динамическом уровне. Управляющие воздействия — моменты двигателей колес. При визуализации исследуемых моделей задействовали встроенные численные методы системы Wolfram Mathematica («Вольфрам математика») с минимальной точностью 10-6.

Результаты исследования. При построении математической модели определили равенства для угловых скоростей колес. Учли наличие площадки контакта и вывели уравнения динамики робота с дифференциальным приводом. Элементами системы были проекции силы и момент, показатели верчения платформы, массы, угловые ускорения и инерция колес. Показали, как формируются управляющие воздействия в рамках неголономной механики. Описали модель двигателей, которые создают момент управления на ведущих колесах. Вывели решение как взаимосвязь индуктивности проводников электродвигателей и работы электропитания. Детально рассмотрели три модели, описывающие динамику робота с дифференциальным приводом. Первая — неголономная. Вторая и третья включали систему уравнений динамики робота с дифференциальным приводом для общего случая, при наличии площадки контакта. При этом во второй модели игнорируется время переключений в двигателе и задействуется трение Кулона. В третьей ввели параметр, определяющий скорость переходных процессов в двигателе, и задействовали разложения Паде. Это модель с поликомпонентным трением. Итоги расчетов показали в виде графиков. На них исследуемые модели визуализировали в виде кривых разного цвета. Сопоставление графиков показало, в каких случаях после завершения переходных процессов управление обеспечит требуемую точность.

Это модели 1 и 2. В модели 3 программное управление генерирует ошибку в угловой скорости вращения платформы. Данную ошибку нельзя спрогнозировать в рамках 1-й и 2-й модели. Во всех рассмотренных системах

© Салимов М.С., Меркурьев И.В. 2023

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

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

Ключевые слова: программное управление роботом, управление в неголономной механике, модель трения Кулона, поликомпонентное трение, ошибка угловой скорости вращения

Благодарности. Авторы выражают признательность кафедре «Робототехника, мехатроника, динамика и прочность машин» ФГБОУ ВО «НИУ «МЭИ» за помощь в подготовке статьи. Благодарим рецензентов за конструктивную критику и предложения, которые способствовали значительному повышению качества статьи.

Для цитирования. Салимов М.С., Меркурьев И.В. Программное управление движением робота с дифференциальным приводом при разных моделях трения. Advanced Engineering Research (Rostov-on-Don). 2023;23(4):346-355. https://doi.org/10.23947/2687-1653-2023-23-4-346-355

Original article

Software Control of the Movement of a Differential Drive Robot for Different Friction Models

Maksim S. Salimov M, Igor V. Merkuryev

National Research University "Moscow Power Engineering Institute", Moscow, Russian Federation И SalimovMS@mpei.ru

Abstract

Introduction. Designing motion control systems for mobile robots requires the construction of mathematical models. Researchers have repeatedly addressed this topic. In particular, works have been published on the calculations of multiphysical processes, modeling the movement of various types of wheels under certain conditions. In addition, the dynamics of deformable contacting bodies during sliding, rolling and rotation, issues of autonomy and controllability of mobile robots were considered. Note, however, that the dynamics and positioning accuracy of wheeled robots is largely determined by friction. The literature does not present studies on the dynamics of a robot with a differential drive taking into account the interrelationships of sliding, spinning and rolling friction effects based on the theory of multicomponent friction. Research in this area can reveal new dynamic effects. Based on the data obtained in this way, it is possible to improve the accuracy of positioning in building mathematical control models. The presented work aims at investigating the movement of an automatic device with a differential drive taking into account three contact models: nonholonomic, Coulomb friction, and multicomponent models.

Materials and Methods. The scheme of a two-wheeled robot with differential drive and continuous movement on the support surface was adopted as the basic one. The movement of the device was provided through software control. The dynamics was described in the form of Appel equations. Mathematical models were used for calculations, taking into account friction in different ways. Coordination of the actions of the mechanism was formed at a dynamic level. The control actions were the moments of the wheel motors. When visualizing the models under study, the built-in numerical methods of the Wolfram Mathematica system were used with a minimum accuracy of 10-6.

Results. When building a mathematical model, the equations for the angular velocities of the wheels were determined. The authors took into account the presence of a contact site and derived the equations of dynamics of a differential drive robot. The elements of the system were force and moment projections, indicators of platform spin, masses, angular accelerations, and inertia of the wheels. It was shown how control actions were formed within the framework of nonholonomic mechanics. The model of engines that created a moment of control on the driving wheels was described. The solution was derived as the relationship between the inductance of the conductors of electric motors and the operation of the power supply. Three models describing the dynamics of a differential drive robot were examined in detail. The first model was nonholonomic. The second and third included a system of equations for the dynamics of a differential drive robot for a general case with a contact platform. At the same time, in the second model, the switching time in the engine was ignored and the Coulomb friction was involved. In s the third model, a parameter to determine the speed of transients in the engine was introduced, and Pade decomposition was ах

е

involved. This was a model with multicomponent friction. The calculation results were shown in the form of graphs. On them, g the studied models were visualized in the form of curves of different colors. Comparison of the graphs showed in which cases, after the completion of transients, the control provided the required accuracy. These were models 1 and 2. In model 3, the

В

a

о

T3

<D f

software control generated an error in the angular velocity of rotation of the platform. This error could not be predicted within the framework of the 1 st and 2nd models. In all the systems considered, the sliding speed of the wheels in the transverse direction dropped to zero. The condition of continuous motion of the support wheel was obtained and validated. Discussion and Conclusion. Software control is acceptable in models that do not take into account wheel friction during simultaneous sliding, spinning and rolling (general case of spatial motion). However, it is important to consider the relationship between these processes and multicomponent friction. This is required for the robot to perform program movements more accurately. It was established that software control in a model that takes into account the friction of spinning and rolling caused deviations from the program values of the angular velocity of the platform. The results obtained can be used in the building of a control system with predictive models.

Keywords: software control of a robot, control in nonholonomic mechanics, Coulomb friction model, multicomponent friction, error in angular velocity of rotation

Acknowledgements. The authors appreciate the Department of Robotics, Mechatronics, Dynamics and Strength of Machines, National Research University MPEI, for assistance in preparing the article. We would like to thank the reviewers for their criticism and suggestions, which contributed to a significant improvement in the quality of the article.

For citation. Salimov MS, Merkuryev IV. Software Control of the Movement of a Differential Drive Robot for Different Friction Models. Advanced Engineering Research (Rostov-on-Don). 2023;23(4):346-355. https://doi.org/10.23947/2687-1653-2023-23-4-346-355

Введение. В современном мире растет востребованность автоматических мобильных колесных аппаратов, активно исследуются возможности координации их работы. Значительный научный и прикладной интерес имеют вопросы автономности и управляемости мобильных роботов. В частности, рассматриваются подходы к расчетам мультифизических процессов, моделированию движения конкретных видов колес в разных условиях. Некоторые изыскания фокусируются на динамике деформируемых контактирующих тел при скольжении, качении и вращении. Результаты этих исследований находят применение в робототехнике [1-3]. Управление мобильными колесными аппаратами рассматривается на примере электрических самокатов «Сигвей» (Segway) [4].

Мобильные механизмы с простыми и эффективными системами управления движением колес нужны в первую очередь для работы в узких пространствах [5], для проведения исследований в сложных промышленных условиях [6].

В [7-11] показаны возможности точного моделирования динамики колесных аппаратов с учетом скольжения, верчения и качения колес.

При этом в литературе не описаны явления, возникающие при комбинации эффектов трения. Цель представленной работы — исследовать влияние поликомпонентного трения на динамику робота с дифференциальным приводом. Задачи: построить математическую модель движения робота с поликомпонентным трением и смоделировать программные движения.

Материалы и методы. Рассматривается схема двухколесного робота с дифференциальным приводом (рис. 1). Предполагается безотрывное движение колес радиуса R по опорной поверхности.

Л"

Z

а) б)

Рис. 1. Кинематическая схема робота: а — платформа; б — одно из ведущих колес

3

Центр тяжести платформы С смещен относительно геометрического центра платформы А по оси X на расстояние ^ , по оси У — на расстояние h . Точки с1 и с2 — центры масс колес (индекс 1 соответствует

левому колесу, индекс 2 — правому).

Вводятся следующие обозначения: т1 — масса платформы робота; ]Сг — момент инерции платформы

относительно центра масс С; 01, х1, у1 — неподвижная система координат; хА, уА определяют положение точки

А в 01, х1, у1; у — угол курса платформы; ф. — углы поворота колес (индекс 1 соответствует левому колесу,

индекс 2 — правому); АС1 = АС2 = Ь — расстояние между колесами; М1,М2 — моменты двигателей,

приложенные к колесам.

Предполагается, что на опорное колесо 4 (рис. 1) действует только нормальная реакция опоры. Во время движения опорное колесо располагается на оси X на расстоянии R0 от центра платформы (Я0 < L). Результаты исследования

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

и х = 0, и,2х = 0, и ,1 у = 0, и у = 0, а)

где и, , и, , иКу, и — проекции на оси !и У скоростей точек К1 и К2, являющихся центрами площадок

контакта колес с опорной поверхностью.

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

и„ = о-О£-Rю,,

К1х 1'

иКх = и + О£ -Rю2, (2)

ик,у = и,2 у = иСу = И.

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

и-О Ь и + ОЬ

(3)

1 R " R

В общем случае, когда предполагается наличие площадки контакта, получим уравнения динамики робота с дифференциальным приводом:

m1 [и-hyQ-Q(u +hxQ)] + 2тк (U-Qu) = F^ cosa1 + F cosa2 -F^sin a1 -F^sin a2, m1 [U + hxQ + Q(u-hyQ)] + 2mK (u + Qu) = F sina1 + F sina2 -F±icos a1 -F±2cos a2, Q (JCz + 2JKz + m1 (h2 + h2y ) + 2m KL ) + m1 [(U + Qu)hx - hy (U - Qu)] = Mс + MC^ +(F cosa2 -F±2sin a2 -F cosa1 + F^sin a1)L, JKy( = M1 - My -(F cosa1 - F±i sina1) R,

JKy(2 = m2 - My2 -(Ft cosa2 - fi2 sina2)R. Здесь F,F — проекции силы трения на направления скоростей скольжения центров пятен контакта колес; F± , F± — проекции силы трения на направления, перпендикулярные скоростям центров пятен контакта колес; Mc,Mc — слагаемые момента трения при верчении для двух колес; Q — угловое ускорение верчения платформы; (1,(2 — угловые ускорения колес, проецируемые на их ось вращения; mK — масса каждого колеса; JKy и JKz — моменты инерции колеса относительно соответствующих осей.

Для формирования управляющих воздействий будем использовать модель динамики робота, описанного в рамках неголономной механики [12]:

m1 [и -hy Q -Qhx Q J + 2mK U

2JKy . M1 + M 2 s

> +—и = —1--, Я

R2 R S

a

о

T3

<D

m1 [(+//,6 + О (и - hyО))hx -hy (и - hyО - О/О)] +

■ I I 2 Л ь2 ■ (м2 - м) ь (5)

О(/сг + 2JKI + т $ + /у2) + 2ткИ) + 2JKy— О ^ 2 я '' .

Ее также можно представить в виде:

, мО = ( м 2- м 1) ь, (6)

и я О я

где Fu и мО — левые части уравнения (5).

Модель двигателей на ведущих колесах. Опишем модель двигателей, способных генерировать момент управления на ведущих колесах м,, м 2.

Уравнения движения колесных двигателей постоянного тока можно выразить формулами [13]:

^ -ёг+я'к=и^ (7)

М к = ^, (k = 1,2).

Здесь Ье — индуктивность проводников электродвигателей; Я — сопротивление электрическому току в цепи двигателя; Пк — электропитание для двигателя к; iк — ток, проходящий через цепь якоря двигателя; с1, с2 — конструкционные постоянные двигателей.

Предполагается, что двигатели ведущих колес имеют одинаковую характеристику. Из (7) видно, что

к = ^ (8)

С2

Учитывая (8), в первом уравнении из системы (7) получим дифференциальное уравнение относительно момента мк :

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

-Мк + ЯМк = ик - ^. (9)

с2 с2 я

Уравнение (9) разделим на множитель — и получим:

с2

ЬеС2Мк + Мк = С2 ик - ^ Шк. (10)

с2 Я Я Я к

Введем обозначения:

t = k = k2 = (11)

1 c2 R 1 R 2 R

В таком случае управляющие моменты можно описать:

тМк + Мк = киик - КШк, (12)

где К1, К2 — константы; Т1 — параметр, определяющий скорость переходных процессов в электрической части двигателей.

При игнорировании времени переключений в электромеханической системе двигателя:

Мк = К ик - Кш. (13)

Авторы [14] получили похожие формулы.

Исследуемые модели. Рассмотрим модели, описывающие динамику робота с дифференциальным приводом. Модель 1 описывается уравнениями (6) и (12), когда длительность переходных процессов полагается исчезающе малой T1 = 0 :

г. , ^ „п „ . 2JK К и - К.ю. + К 2и 2 - К.ш7 т. I и-А„О-О/О +2тки +—^и = -LJ-—-—,

1 [ у х J к я2 Я

3 ш1 [(+hО + О (u- hyО))hx -hy (и - hyО - OhxО)J + 2Л,, — О +

L •

Ky R2

,( T + 2 T + (h2 + h2 ) + 2 ,2 ) ((K 2U 2 - K1^2 )-(K 2U1 - K1^1 )) L (JCz + 2JKz + (К + hj, ) + 2^L )=-R-■

(14)

Hi +О (.......... . .

R

> Здесь ,ш2 вычисляются по формулам (3).

^ Для описания модели 2 применяются уравнения (4), (13), а также модель трения Кулона:

^ Я=- Ж-, F1 = 0, Мс = 0. (15)

Модель 3 представлена уравнениями (4), (12). Кроме того, учитывается поликомпонентное трение, представленное в [12] с помощью разложений Паде:

F = - fN

F, =- fNr

__2I1 - krI2 cos I — a | ,

V^KX+^Kr+kr/112Q-^-- cos I - -a

I 2

i

2I1 - krI2 cos I — a

U * + UK, , +-

I

к n Q

(16)

(17)

2/1 - kr12 cos I — a kr I112-:—--- cos | — -a +ui„ + |Q|

/„

2-a

Ky

kr 13 cos b-a |/

Mr =-2 fnl 2

2I

y2 -krcos --a- >/°K, x +uK, , +Q

-krcos-a| Uu2kx + UK,y + |Q|

M, =-kr nI3 T^i,. ' Ю,.

(18)

(19)

Здесь 11,12, /3 — постоянные величины, которые находятся по формулам в [12]; а — угол, характеризующий направление линейных скоростей скольжения относительно точек контакта колес К1,К2; кг — коэффициент, знак которого зависит от направления качения.

Для определения программного управления следует указать законы движения робота (/), ОЛ (/). По (5)

можно выразить величины = (/), = /О (/), МО = МО (() в зависимости от времени. Используем известные значения /,Ма для вычисления управляющих напряжений Ц и и2, подаваемых на приводы колес:

K2 - 2

№ =-L W F, + Mq| + K>21.^1 i4i F„- Mq| + a>,

L - 12 - 1 K2-2 - L Движение системы с заданными начальными условиями и( 0 ) = 0,5—, u (0 ) = 0,2—, Q (0 ) = 0,2

(20) 1

■ М ^ „ м

—, u ( 0 ) = 0,2—

с с

смоделировали в математическом пакете «Вольфрам математика» (Wolfram Mathematica).

Рассмотрим и (t) = 0,5 —, Od (t) = 0,3-. Для интеграции системы автоматически выбирали шаг

c c

интегрирования и задействовали встроенные численные методы Wolfram Mathematica с точностью результатов не менее 10-6.

v(t), м/c

0,5

0,4

0,3 0,2

0,1

t, c

5

10

15

20

Рис. 2. Зависимость продольной скорости и. Зеленая кривая — неголономная модель 1, красная кривая — модель 2 с трением Кулона, синяя кривая — модель 3 с поликомпонентным трением.

Синяя и красная кривые совпадают

ак и

и

е

0

u(t), м/с

t, с

Рис. 3. Зависимость поперечной скорости и. Зеленая кривая — неголономная модель 1, красная кривая — модель 2 с трением Кулона, синяя кривая — модель 3 с поликомпонентным трением.

Зеленая, синяя и красная кривые совпадают

^(0, 1/с

t, с

5

10

15

20

Рис. 4. Зависимость угловой скорости робота. Зеленая кривая — неголономная модель 1, красная кривая — модель 2 с трением Кулона, синяя кривая — модель 3 с поликомпонентным трением

В

а

о

43

(D

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

I

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

Для проверки выполнения условий безотрывного движения опорного колеса запишем теорему изменения кинетического момента системы в проекциях на ось У:

Jку(А + ю2) =М1 + м2 -МУ1 -мУг -N3К + С0^ -^ япа + ^ С^а2 -^ япа2)К

Согласно (4):

М1 + М 2 = - N3 К0 + т1 ghx.

Из последнего уравнения получаем условие безотрывного движения опорного колеса:

т1 ghx - (М + М2) > 0.

С учетом (5):

2J.

m1 ghx -R(m1 + 2mk +—+ Rm1 [hyQ + hxQ2 ] > 0. R

Ускорения робота не превышали 1 м/c2, поэтому данное неравенство выполняется в каждый момент времени.

Обсуждение и заключение. Итоги моделирования позволили сделать ряд выводов. Программное управление, основанное на неголономной модели 1, приемлемо для роботов, конструкция которых не учитывает контактное трение колес, возникающее при комбинации скольжения, верчения и качения. Однако при формировании управляющих воздействий для достижения более точного следования программным законам движения необходимо учитывать модель поликомпонентного трения. Полученное решение можно использовать при построении системы управления с прогнозирующей моделью (model predictive control). Такая разработка будет предметом дальнейшего исследования.

Список литературы

1. Gerasimov K.V., Zobova A.A., Kosenko I.I. Omni-Vehicle Dynamical Models Mutual Matching for Different Roller-Floor Contact Models. In: Proc. European Congress on Computational Methods in Applied Sciences and Engineering "MultibodyDynamics". Cham: Springer; 2019. P. 511-517. http://doi.org/10.1007/978-3-030-23132-3 61

2. Zobova A.A., Habra T., Van der Noot N., Dallali H., Tsagarakis N.G., Fisette P., et al. Multi-Physics Modelling of a Compliant Humanoid Robot. Multibody System Dynamics. 2017;39:95-114. https://doi.org/10.1007/s11044-016-9545-4

3. Gerasimov K.V., Zobova A.A. On the Motion of a Symmetrical Vehicle with Omniwheels with Massive Rollers. Mechanics of Solids. 2018;53:32-42. https://doi.org/10.3103/S0025654418050060

4. Corke P. Robot Arm Kinematics. In book: Robotics, Vision and Control, 2nd ed. Cham: Springer; 2017. P. 193-228. https://doi.org/10.1007/978-3-319-54413-7 7

5. Junjie Shen, Dennis Hong. OmBURo: A Novel Unicycle Robot with Active Omnidirectional Wheel. In: Proc. IEEE Int. Conf. on Robotics and Automation (ICRA). New York City: IEEE; 2020. P. 11-19. https://doi.org/10.1109/ICRA40945.2020.9196927

6. Munitsyn L.V. Vibrations of a Rigid Body with Cylindrical Surface on a Vibrating Foundation. Mechanics of Solids. 2017;52:675-685. https://doi.org/10.3103/S0025654417060073

7. Goryacheva I.G., Zobova A.A. Dynamics of Deformable Contacting Bodies with Sliding, Rolling, and Spinning. International Journal of Mechanical Sciences. 2022;216:106981. https://doi.org/10.1016/Uimecsci.2021.106981

8. Zobova A.A., Goryacheva I.G. Effect of Contacting Bodies' Mechanical Properties on the Dynamics of a Rolling Cylinder. ActaMechanica. 2020;232:1971-1982. http://doi.org/10.1007/s00707-020-02800-w

9. Киреенков А.А., Жаворонок С.И., Нуштаев Д.В. О моделях шины, учитывающих как деформированное состояние, так и эффекты сухого трения в области контакта. Компьютерные исследования и моделирование. 2021;13(1): 163-173. https://doi.org/10.20537/2076-7633-2021-13-1-163-173

10. Kireenkov A.A., Zhavoronok S.I. Coupled Dry Friction Models in Problems of Aviation Pneumatics' Dynamics. International Journal of Mechanical Sciences. 2017;127:198-203. http://doi.org/10.1016/uimecsci.2017.02.004

11. Kireenkov A.A., Fedotenkov G.V., Shiriaev A., Zhavoronok S.I. Application of the Theory of the Multicomponent Dry Friction in Some of Control Robot Systems. International Journal for Computational Civil and Structural Engineering. 2022;18(1):14-23. https://doi.org/10.22337/2587-9618-2022-18-1-14-23

12. Salimov M.S., Saypulaev G.R., Merkuriev I.V. Four-Dimensional Integral Model of Dry Friction on the Example of Wheel Movement. Journal of Physics: Conference Series. 2021;2096:012043. https://doi.org/10.1088/1742-6596/2096/1/012043

13. Alaci S., Ciornei F.C., Pentiuc R.D., Ciornei M.C., Romanu I.C. Indetermination versus Incompatibility in Dynamic Systems with Dry Friction. Journal of Physics: Conference Series. 2020;1426:012011. https://doi.org/10.1088/1742-6596/1426/1/012011

14. Long Chen, Jun Liu, Hai Wang, Youhao Hu, Xuefeng Zheng, Mao Ye, et al. Robust Control of Reaction Wheel Bicycle Robot via Adaptive Integral Terminal Sliding Mode. Nonlinear Dynamics. 2021;104:2291-2302. https://doi.org/10.1007/s11071-021-06380-9

References

1. Gerasimov KV, Zobova AA, Kosenko II. Omni-Vehicle Dynamical Models Mutual Matching for Different Roller-Floor Contact Models. In: Proc. European Congress on Computational Methods in Applied Sciences and Engineering "MultibodyDynamics". Cham: Springer; 2019. P. 511-517. http://doi.org/10.1007/978-3-030-23132-3 61

2. Zobova AA, Habra T, Van der Noot N, Dallali H, Tsagarakis NG, Fisette P, et al. Multi-Physics Modelling of a Compliant Humanoid Robot. Multibody System Dynamics. 2017;39:95-114. https://doi.org/10.1007/s11044-016-9545-4 §

К

3. Gerasimov KV, Zobova AA. On the Motion of a Symmetrical Vehicle with Omniwheels with Massive Rollers. |§ Mechanics of Solids. 2018;53:32-42. https://doi.org/10.3103/S0025654418050060 |

4. Corke P. Robot Arm Kinematics. In book: Robotics, Vision and Control, 2nd ed. Cham: Springer; 2017. P. 193-228. https://doi.org/10.1007/978-3-319-54413-7 7

5. Junjie Shen, Dennis Hong. OmBURo: A Novel Unicycle Robot with Active Omnidirectional Wheel. In: Proc. IEEE Int. Conf. on Robotics and Automation (ICRA). New York City: IEEE; 2020. P. 11-19. https://doi.org/10.1109/ICRA40945.2020.9196927

6. Munitsyn LV. Vibrations of a Rigid Body with Cylindrical Surface on a Vibrating Foundation. Mechanics of Solids. 2017;52:675-685. https://doi.org/10.3103/S0025654417060073

7. Goryacheva IG, Zobova AA. Dynamics of Deformable Contacting Bodies with Sliding, Rolling, and Spinning. International Journal of Mechanical Sciences. 2022;216:106981. https://doi.org/10.1016/uimecsci.2021.106981

8. Zobova AA, Goryacheva IG. Effect of Contacting Bodies' Mechanical Properties on the Dynamics of a Rolling Cylinder. ActaMechanica. 2020;232(5): 1971-1982. http://doi.org/10.1007/s00707-020-02800-w

9. Kireenkov AA, Zhavoronok SI, Nushtaev DV. On Tire Models Accounting for Both Deformed State and Coupled Dry Friction in a Contact Spot. Computer Research and Modeling. 2021;13(1): 163-173. https://doi.org/10.20537/2076-7633-2021-13-1-163-173

10. Kireenkov AA, Zhavoronok SI. Coupled Dry Friction Models in Problems of Aviation Pneumatics' Dynamics. International Journal of Mechanical Sciences. 2017;127:198-203. http://doi.org/10.1016/uimecsci.2017.02.004

11. Kireenkov AA, Fedotenkov GV, Shiriaev A, Zhavoronok SI. Application of the Theory of the Multicomponent Dry Friction in Some of Control Robot Systems. International Journal for Computational Civil and Structural Engineering. 2022;18(1):14-23. https://doi.org/10.22337/2587-9618-2022-18-1-14-23

12. Salimov MS, Saypulaev GR, Merkuriev IV. Four-Dimensional Integral Model of Dry Friction on the Example of Wheel Movement. Journal of Physics: Conference Series. 2021;2096:012043. https://doi.org/10.1088/1742-6596/2096/1/012043

13. Alaci S, Ciornei FC, Pentiuc RD, Ciornei MC, Romanu IC. Indetermination versus Incompatibility in Dynamic Systems with Dry Friction. Journal of Physics: Conference Series. 2020;1426:012011. https://doi.org/10.1088/1742-6596/1426/1/012011

14. Long Chen, Jun Liu, Hai Wang, Youhao Hu, Xuefeng Zheng, Mao Ye, et al. Robust Control of Reaction Wheel Bicycle Robot via Adaptive Integral Terminal Sliding Mode. Nonlinear Dynamics. 2021;104:2291-2302. https://doi.org/10.1007/s11071-021-06380-9

Поступила в редакцию 20.09.2023 Поступила после рецензирования 18.10.2023 Запланирована в номер 07.11.2023

Об авторах

Максим Сергеевич Салимов, ассистент кафедры робототехники, мехатроники, динамики и прочности машин Национального исследовательского университета «МЭИ» (111250, РФ, г. Москва, ул. Красноказарменная, 14), SPIN-код: 7721-7434, ORCID, SalimovMS@mpei.ru

Игорь Владимирович Меркурьев, доктор технических наук, доцент, заведующий кафедрой робототехники, мехатроники, динамики и прочности машин Национального исследовательского университета «МЭИ» (111250, РФ, г. Москва, ул. Красноказарменная, 14), SPIN-код: 1608-9638, ScopusID, ORCID, MerkuryevIV@mpei. ru

Заявленный вклад соавторов

М.С. Салимов — формирование цели и задачи исследования, подготовка текста, расчеты.

И.В. Меркурьев — научное руководство, анализ результатов исследований, доработка текста, корректировка выводов.

Конфликт интересов: авторы заявляют об отсутствии конфликта интересов.

Все авторы прочитали и одобрили окончательный вариант рукописи.

В

3 Received 20.09.2023

! Revised 20.10.2023

2 Accepted 07.11.2023

•s

s e

v/

I

About the Authors:

Maksim S. Salimov, Teaching assistant of Department of Robotics, Mechatronics, Dynamics and Strength of Machines, National Research University MPEI (14, Krasnokazarmennaya St., Moscow, 111250, RF), SPIN-code: 7721-7434, ORCID, SalimovMS@mpei.ru

Igor V. Merkuryev, Dr.Sci. (Eng.), Associate Professor, Head of the Department of Robotics, Mechatronics, Dynamics and Strength of Machines, National Research University MPEI (14, Krasnokazarmennaya St., Moscow, 111250, RF), SPIN-code: 1608-9638, ScopusID, ORCID, MerkuryevIV@mpei.ru

Claimed contributorship:

MS Salimov: research objectives and tasks formulation, text preparation, computational analysis.

IV Merkuryev: academic advising, analysis of the research results, revision of the text, correction of the conclusions.

Conflict of interest statement: the authors do not have any conflict of interest.

All authors have read and approved the final manuscript.

3

s

и

<D

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