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

УДК 621.864.8

РАЗРАБОТКА АЛГОРИТМА КОМПЬЮТЕРНОГО МОДЕЛИРОВАНИЯ ДВИЖЕНИЯ МОБИЛЬНОГО МИНИРОБОТА, ПЕРЕМЕЩАЮЩЕГОСЯ С ОТРЫВОМ ОТ ОПОРНОЙ ПОВЕРХНОСТИ

© 2011 А Н. Рукавицын, И.В. Лупехина

Юго-Западный государственный университет, г. Курск

Поступила в редакцию 11.11.2011

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

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

Мобильные автономные устройства, движущиеся с отрывом от опорной поверхности (так называемые прыгающие роботы), являются объектами исследования ученых многих стран, поскольку скачкообразный способ перемещения при движении по сильно пересеченной местности удобнее скольжения или качения. Большинство подобных устройств создаются на основе принципов бионики - копирования движения некоторых животных, способных накапливать энергию в мышцах и затем во время прыжка быстро её высвобождать [1, 2]. Другие методы проектирования основываются на применении подвижных внутренних масс, встраиваемых в корпус робота [3, 4]. В [5] было представлено мобильное устройство, двигающиеся с отрывом от поверхности, за счет периодического движения внутренней вращающейся массы (см. рис. 1). Это позволяет осуществлять управление движением за счет изменения одного параметра - угловой скорости вращения дебаланса.

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

Для установления области параметров системы, при которых корпус робота отрывается от поверхности и начинает вращение вокруг точки приложения реакции, составим уравнение моментов относительно центра масс робота -точки С с относительными координатами (^с,Пс) (рис. 2).

Рукавицын Александр Николаевич, кандидат технических наук, доцент кафедры теоретической механики и мехатроники. E-mail: alruk75@mail.ru Лупехина Ирина Владимировна, аспирантка. E-mail: teormeh@inbox.ru

Рис. 1. Схема движения прыгающего миниробота

Y

п L,' N

г Пс • 1

О' 1 1

хЛ

¿.с

Рис. 2. Схема приложения сил и моментов

Плечо А нормальной реакции N должно

находиться в интервале > ^з ), где Lх -длина корпуса. Уравнение моментов и в случае покоя, и в случае скольжения имеет вид:

N-A-Ffrlc = 0

(1)

Принимая во внимание определение центра масс робота, для его координат имеем:

п

Ща = +Е тг^г,

1=1 п

М^а = тПк тгПг ■

г=1

(2)

Здесь М = т + ^ т - общая масса корпуса и

г=1

масс; пК - координаты центра тяжести корпуса. Тогда, учитывая ограничения для плеча, придем к неравенству

0 < ^ + Щс < Жх

(3)

с < тт{^с, 4 -£с }

(5)

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

Очевидно, что в зависимости от суммарных проекций сил инерции корпус робота может оставаться неподвижным, может перемещаться по предоставленной поверхности без отрыва или двигаться с отрывом от нее. При этом в общем случае движение робота в пространстве является свободным. Физические условия контакта корпуса робота с поверхностью, по которой он перемещается, таковы, что поверхность сопротивляется проникновению корпуса внутрь, но не сопротивляется отрыву. Поэтому сила нормального давления Щ, вызванная контактом корпуса робота с плоскостью опоры, может быть направлена вертикально вверх или равняться нулю, т.е., удовлетворять условию N<0. Если в дальнейшем движение корпуса происходит без отрыва от плоскости опоры, то, в соответствии с (2), получим ограничение на относительное ускорение внутреннего тела, колеблющегося вдоль вертикали

т + Х т• 18 +Х - 0

1=1 ) 1=1

(6)

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

Если внутренние массы движутся таким образом, что в некоторые промежутки времени неравенство (6) нарушается, то корпус миниро-бота будет отрываться от поверхности, а движение робота приобретет прыжкообразный характер [6]. Тем не менее, выполнение неравенства (6) не гарантирует прижатия всей грани корпуса, так как нормальная реакция поверхности является силой, распределенной по площади соприкосновения. Ее интенсивность п(х) может принимать только неотрицательные значения, поскольку отсутствует прилипание, и условно картина распределения имеет вид, аналогичный показанному на рис. 3.

Рис. 3. Схема распределения нормальной реакции поверхности

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

как интеграл М0 (Щ) = | п(х) хйх, равен нулю, что

о

при х>0 выполнится лишь при п(5)=0 во всех точках касания, кроме нулевой. Если точка приложения является крайней правой точкой нагрузки, то разность между моментами равнодействующей и нагрузки равна:

N ■ Ьх + | п(х)хdх = | п(х)(£х - х)аХ

(6)

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

п

п

о

о

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

N ®(t)

« \

7//////// ///////// 7 mg

О (t) = О0 sin(^t + )

(7)

F, = -О ~ (t) = -О0 sin(^t + (р0) cos а

(8)

Fr = N

J ' fYlQV "

(9)

значение силы трения покоя, т.е., перестанет выполняться неравенство

О (t)|< jJN

(10)

В этом случае тело приобретает ускорение, определяемое разностью проекций вынуждающей силы и силы трения

mx = Ох (t) - f sgn Ох (t) • N

(11)

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

Рис. 4. Схема сил, приложенных к роботу

Пусть модель перемещается по горизонтальной шероховатой плоскости. Будем рассматривать движение под действием силы тяжести mg, нормальной реакции N, силы сухого трения Ffr , вязкого сопротивления R и периодической силы Ф(), направленной к поверхности под углом а (рис.2). Модуль вынуждающей силы изменяется по закону:

mx = Ох (t) - f sgn о • N

(12)

где ш - частота, ф0 - начальная фаза.

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

Направление силы трения в покое противоположно проекции силы Ф(?). Пока тело остается неподвижным значение силы трения не превышает максимального, определяемого равенством

где f - коэффициент трения, N = mg - О (t)sina - значение нормальной реакции, которое предполагаем всегда положительным.

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

Здесь и в дальнейшем будем рассматривать кусочно-непрерывную функцию проекций Х скоростей или внешних сил (кроме трения) на ось х: f(X) = sgn X, которая позволяет определить коэффициент и направление трения одновременно. Скорость и координату, как в момент начала скольжения, так и в процессе движения будем определять численно при помощи разложения в степенной ряд.

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

Алгоритм вычисления основывается на следующих допущениях. Прежде всего, в движении выделяем две фазы: I - полет, II - движение по поверхности. Переход из одной фазы в другую происходит или при выполнении условий отрыва (II—^1), или при выполнении условий приземления (I—II). Условием отрыва является равенство нулю нормальной реакции поверхности, условие приземления - невыполнение неравенства у1>у1ть для значений вертикальной координаты центра проекции корпуса. Величина минимального допускаемого значения координаты центра проекции у 1шП зависит

от выбора начала отсчета. Например, если горизонтальная координатная плоскость проходит через центр, то у1т/п=0. Программа организована в виде цикла, выход из которого осуществляется,

как только расчетное значение времени превысит заданную величину. Для численного расчета параметров движения использовался метод Эйлера II порядка.

Рис. 5. Блок-схема алгоритма

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

у(х) = у0 + у'(х0)(х - х0) + 1 у"(х0)(х - х0)2 +...

2 .(13)

Таким образом, ускорение каждой из обобщенных координат qj, ]=1,2,3, будем находить из правой части уравнений системы (2) в случае полета, или из уравнений (1) в случае

движения по поверхности. Тогда, если (Я , Х- ускорение обобщенной координаты в начале /-го шага, т.е., определенное для предшествующего шагу момента времени 1м=А (/ - 1), то скорость и координата, согласно разложению, определятся выражениями:

(Ч,) = (Ч,) г-1 + (9,) г-1 'А

1

(14)

) г = ) г-1 + (Ч,) г-1 + ^(9}) г-1 ■ А2

' ' ' 2 ' . (15)

Здесь А - выбранный шаг по времени.

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

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

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

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

Работа выполнена в рамках реализации Федеральной целевой программы «Научные и научно-педагогические кадры инновационной России» на 2009-2013 годы (гос. регистр. № П699, шифр НК-617П-4).

СПИСОК ЛИТЕРАТУРЫ:

1. Miyazaki, M. Jumping via robot body deformation -Mechanics and mechanism for higher jumping / M. Miyazaki, S. Hirai // Advances in Mobile robotics. Proceedings of the 11 International Conference on Climbing and Walking Robots and the Support Technologies for Mobile Machines. Coimbra. Portugal, 2008. P. 373-380.

2. Larin, V.B. Concerning the designing of the hopping apparatus / V.B. Larin, V.M. Matiyasevich // Proceedings of the Fifth International Conference on Climbing and Walking Robots and their Supporting Technologies. CLAWAR 2002. P. 365-372.

3. Черноусько, Ф.Л. О движении тела, содержащего подвижную внутреннюю массу / Ф.Л. Черноусько // ДАН. 2005. Т. 405, № 1. С. 1-5.

4. Болотник, Н.Н. Динамика управляемых движений вибрационных систем / Н.Н. Болотник, И.М. Зейдис, К. Циммерманн и др. // Изв. РАН. ТиСУ. 2006. №5. С. 157-167.

5. Яцун, С.Ф. Исследование движения прыгающего миниробота / С.Ф. Яцун, И.В. Лупехина, А.Н. Ру-кавицын // Известия высших учебных заведений. Северо-кавказский регион. Серия: Технические науки. 2011. №2. С. 3-9.

6. Лупехина, И.В. Моделирование движения прыгающего вибрационного микроробота / И.В. Лу-пехина, К.А. Сапронов, С.Ф. Яцун // Известия Курского государственного технического университета. - Курск: КурскГТУ. 2009. №5. С. 2531.

7. Журавлев, В.Ф. Основы теоретической механик -М.: Физматлит, 2001. 320 с.

DEVELOPMENT OF COMPUTER SIMULATION ALGORITHM OF MOBILE MINIROBOT MOTION, MOVING WITH LIFT-OFF FROM

THE BASIC SURFACE

© 2011 A.N. Rukavitsyn, I.V. Lupehina

SouthWest State University, Kursk

In article questions of formulation the algorithm of the decision the differential equations describing the motion dynamics of mobile two-mass mechanical system, moving with a lift-off from basic surface are considered. The presented model of the springing minirobot consists of exterior body (case) and internal (debalance), change of which angular speed allows to carry out traffic control.

Key words: springing robot; mathematical model, algorithm, dynamic parameters, motion equation

Alexander Rukavitsyn, Candidate of Technical Sciences, Associate Professor at the Theoretical Mechanics and Mechatronics. E-mail: alruk75@mail. ru

Irina Lupehina, Post-graduate Student. E-mail: teormeh@inbox.ru