CC BY
81
УДК 62-503.5
ЭТАПЫ ДВИЖЕНИЯ ЧЕТЫРЕХЗВЕННОГО РОБОТА, ПЕРЕМЕЩАЮЩЕГОСЯ С ОТРЫВОМ ОТ ПОВЕРХНОСТИ
Локтионова Оксана Геннадьевна д.т.н., доцент
Юго-Западный государственный университет,
Курск, Россия
e-mail: te ormeh®inbox. ru
Яцун Сергей Федорович д.т.н., профессор
Юго-Западный государственный университет,
Курск, Россия,
e-mail: teormehuhinbox. ru
Серебровский Вадим Владимирович д.т.н., профессор
Юго-Западный государственный университет, Курск, Россия
Волкова Людмила Юрьевна аспирантка
Юго-Западный государственный университет, Курск, Россия
e-mail: milal80888(a)yandex.ru
Приведена расчетная схема четырехзвенного робота, перемещающегося с отрывом от поверхности, прыжок устройства описан в виде последовательности этапов, на каждом из которых звенья робота совершают определенные виды движения, разработана математическая модель в виде дифференциальных уравнений, описывающих движение робота на всех этапах, представлены результаты численного моделирования прыжка в соответствии с указанной последовательностью этапов
Ключевые слова: ПРЫГАЮЩИЙ МНОГОЗВЕННЫЙ РОБОТ, ЭТАПЫ ПРЫЖКА, ПОСЛЕДОВАТЕЛЬНОСТЬ ЭТАПОВ, ХАРАКТЕРИСТИКИ ПРЫЖКА, ВЫСОТА, ДЛИНА
UDC 62-503.5
STAGES OF MOVEMENT OF THE FOURLINK ROBOT MOVING WITH THE JUMP FROM THE SURFACE
Foktionova Oksana Gennadievna Dr.Sci.Tech., associate professor
South-West State University, Kursk, Russia e-mail: teormehuhinbox. ru
Jatsun Sergey Fedorovich Dr.Sci.Tech., professor
South-West State University, Kursk, Russia e-mail: teormeh®inbox. ru
Serebrovsky Vadim Vladimirovich Dr.Sci.Tech., professor
South-West State University, Kursk, Russia
Volkova Fyudmila Yurievna postgraduate student
South-West State University, Kursk, Russia e-mail: milal80888(a>vandex. ru
The settlement scheme of the four-link robot moving with a jump from a surface is provided, the jump of the device is described in the form of sequence of stages, on each of which links of the robot make certain types of movement, the mathematical model in the form of the differential equations describing movement of the robot at all stages is developed, results of numerical modeling of a jump according to the specified sequence of stages are presented
Keywords: MUFTIFINK JUMPING ROBOT, STAGES OF JUMP, SEQUENCE OF STAGES, CHARACTERISTICS OF THE JUMP, HEIGHT, FENGHT
Введение
Интерес к прыгающим роботам в последние годы все время возрастает, проектируются и исследуются различные по конструкции и способу прыжка роботы, параметры прыжка (высота и длина) которые также сильно отличаются. Общим для
большинства ученых, занимающихся прыгающими роботами, является представление прыжка робота в виде некоторой последовательности этапов, каждый из которых характеризуется особым видом движения звеньев. Так в [1] рассматриваются четыре этапа: подготовка к прыжку при нахождении объекта на поверхности, разгон устройства до отрыва от поверхности, движение при отсутствии взаимодействия с поверхностью, приземление. Аналогичные этапы прыжка описываются в [2]. Один прыжок робота представлен в виде последовательности восьми этапов в [3]: подготовка к прыжку, разгон, переход из состояния, когда объект взаимодействует с поверхностью, к отрыву от поверхности, отрыв от поверхности, полет вверх, полет вниз, приземление. Миниатюрный робот-кузнечик, описанный в [4], совершает один прыжок за три этапа: подготовка к прыжку, разгон, движение в воздухе. В [5] прыжок робота характеризуется двумя этапами: нахождение на поверхности и движение с отрывом от нее, которые связывают два перехода: разгон и приземление.
Прыжок описанного в данной работе многозвенного прыгающего робота представлен в виде последовательности восьми этапов, каждому из которых соответствуют определенные виды движения звеньев, описанные в виде математической модели с учетом наложенных на этапах связей.
Расчетная схема робота
Исследуемый прыгающий робот состоит из четырех звеньев /=1^-4, связанных между собой приводами вращательного 5, 6 и поступательного 7 движения, звено 1 соответствует стопе, звенья 2 и 3 - ноге, звено 4 - корпусу (рис. 1). Звено 3 установлено в центре симметрии звена 4 - точке С4 [6-9].
%
О
Рис. 1 Расчетная схема многозвенного прыгающего робота: 1 - 4 - звенья,
5 - 7 - приводы
Движение робота происходит в вертикальной плоскости Оху, в которой звенья 1 =1-КЗ представляют собой стержни длинами ^, а звено 4 имеет вид прямоугольника с размерами 2я х 26 . Массы звеньев т> сосредоточены в центрах их симметрии - точках Сг. Углы наклона звеньев к оси Ох равны ф1, Фг, ‘Рз =Фг Ф4 5 относительное перемещение звеньев 2 и 3 - ^ .
Взаимодействие робота с опорной поверхностью
При нахождении робота на поверхности, взаимодействие с последней может происходить звеном 1 или 4. При этом возникают распределенные силы трения и нормальной реакции. Будем считать, что если на поверхности находится звено 1, то
нормальные реакции будут приводиться к крайним точкам 0\ и О2 ( ^°], ^02) звена, а
сила трения & - к точке О2, что обусловлено конструктивными особенностями опорного элемента в этой точке, имеющего форму иглы.
Когда на поверхности находится звено 4, взаимодействие с последней может
происходить в одной из точек А->, ■! 1 - 4, или в двух точках, если одна из сторон
корпуса параллельна оси Ох. В первом случае точка контакта представляет собой
центр приведения сил трения и нормальной реакции, во втором случае приведенные нормальные реакции будут сосредоточены в крайних точках стороны корпуса, находящейся на поверхности, а приведенная сила трения только в одной из данных точек - в последней точке отрыва корпуса от поверхности или в первой точке
приземления. Это объясняется тем, что точки А-> звена 4 оснащены опорными элементами (рис. 2), силы трения в которых с поверхностью являются управляемыми величинами, что достигается применением специальных приводов, изменяющих свойства этих элементов. То есть, в одной из точек контакта звена 4 с поверхностью сила трения может иметь некоторое ненулевое значение, а в другой - быть равной нулю.
Рис. 2 Опорные элементы в точках А] корпуса: 1 - электромагнитный привод, 2 -опорный элемент при ненулевой силе трения, 3 - опорный элемент при нулевой силе
трения
Для удобства точку, в которой возникает сила трения, будем обозначать К. Взаимодействие этой точки с поверхностью может описываться различными моделями силы трения: сухое кулоново трение, трение с уменьшенной динамической составляющей, трение с непрерывно падающей характеристикой (некулоново трение) [10,11]. Положим, что возникающая сила трения описывается моделью Кулона [12]:
*> =
-А*$да.{хк), если хкФ 0;
— если *К = I ^0х I— ,
— J/¡Vsgn(_Fo), если *К = 0’ 11*0X |>
где ^0х - горизонтальная проекция равнодействующей всех сил, кроме силы сухого трения, приложенных к звену, взаимодействующему с поверхностью; ? -
коэффициент трения скольжения; N - нормальная реакция в точке К; Хк - скорость
точки К вдоль оси Ох. Рассмотрим случай, когда значение силы трения всегда меньше предельного
К\<Р^ (2)
то есть, скольжения звена робота, контактирующего с поверхностью, не происходит.
Математическая модель прыжка робота
Для разработки математической модели запишем вектор координат q, которым будет определяться положение звеньев робота
?=(*! У\ Х4 У 4 Ф1 Ф2 Ф4 къТ а (3)
где х1> У1’Х4>У4 - координаты центров масс звеньев 1 и 4, ф1, ф2 и ф4 - углы поворота
звеньев 1, 2 и 4, /г’ - длина ноги робота. Углы ф1 и ф2 будем определять как углы
между соответствующими звеньями и положительным направлением оси Ох, угол ф4
- как угол, под которым сторона корпуса А\А4 наклонена к горизонтальной оси.
Система дифференциальных уравнений движения робота в матричной форме имеет вид
А(д)д + В(д)д2 + С(^)ф2/23 = ^ (4)
где , в, с- матрицы коэффициентов, Р - матрица обобщенных сил.
Этапы прыжка робота
Прыжок устройства начинается из положения, в котором звено 4 робота находится на поверхности, звенья 1, 2 и 3 полностью втянуты в корпус (звено 4) -
/2з=/2зтш - и расположены под некоторыми произвольными углами к горизонту: ф1, ф2,
Фз =Фг и представляет собой последовательность семи этапов (к= 1^7).
1 этап. Звенья 1 - 3, жестко связанные между собой, поворачиваются относительно точки С4 за счет момента М43, создаваемого приводом 5, до выполнения
_ о о
условия ф2 _ф2, где ф2 - угол наклона звеньев 2 и 3 к оси Ох, под которым действует
сила, отрывающая робота от поверхности. Обобщенной координатой выступает угол Ф2 ;
а6бФ2 + ^6бФ2 = ^61 _ (5)
2 этап. Звено 1 поворачивается относительно точки Oj под действием момента
_ о _
формируемого приводом 6, до тех пор, пока не выполнится условие ф1 ф1 71 . Дифференциальное уравнение по ф1 имеет вид
а55Ф1 = ^51 (6)
3 этап. Звенья 1 и 2, жестко связанные между собой, перемещаются относительно
звена 3 при помощи силы ^32, генерируемой приводом 7, до взаимодействия звена 1 с
поверхностью, что описывается дифференциальным уравнением
^88^23 = ^81 _ (7)
4 этап. Жестко связанные между собой звенья 3 и 4 перемещаются относительно
звена 2 за счет создаваемой приводом 7 силы ^23 до отрыва звена 1 от поверхности. К
/°
концу этого этапа все звенья робота приобретают скорость 23, которая является начальной скоростью этапа полета. В качестве обобщенной координаты выступает ^23
а88^23 = ^81 _ (8)
5 этап. Полет робота, во время которого все звенья движутся с отрывом от поверхности. При этом звенья 1 и 2, связанные между собой, втягиваются в корпус за
J7 л 1 — ]т*п
счет силы 32, формируемой приводом 7, до выполнения условия 23 23 . Движение
робота на данном этапе описывается четырьмя обобщенными координатами
Я33Х4 + Я37Ф4 + Яз8/2з + ¿37Ф4 + сзз4зФ4 - *31
О44У4 + Я47Ф4 + «48^23 + ^47Ф4 + ^44^23Ф4 “ ^41
77 ф 4 “t" (173¿4 “Ь 743^4 + ^78^23 ^77^4 ^77^2зФ4 — ^ 71
а88^23 + а83*4 + а84^4 + а87Ф4 + ^87Ф4 “ ^81
(9)
(10)
(П)
(12)
6 этап. Приземление робота на одну или две точки корпуса. В первом случае
скорость точки контакта становится равной нулю, а все звенья робота приобретают
*
некоторую угловую скорость ф4 , после чего прыжок переходит в следующий этап. Во втором случае при приземлении на две точки корпуса их скорости обнуляются, и объект становится неподвижным, прыжок завершается.
7 этап. Все звенья робота поворачиваются относительно точки контакта до тех
пор, пока вторая точка корпуса не начнет взаимодействовать с поверхностью. Дифференциальное уравнение по обобщенной координате щ записывается
следующим образом
Я77Ф4 + 677Ф4 ^ (13)
В формулах (5) - (13) использованы обозначения
к
тп 5 т=1^-8, п=1^-8 - коэффициенты при вторых производных обобщенных координат,
тк
тп 5 /й=1ч-8, п=1-^8 - коэффициенты при квадратах первых производных обобщенных координат,
трк
т„, /и=1-^8, п=1-^8 - обобщенные силы по обобщенным координатам.
Управляющие моменты и силы
Для исследуемого робота управление параметрами прыжка происходит на этапах позиционирования объекта и его разгона до отрыва от поверхности (этапы 1 - 4), а также непосредственно во время полета (этап 5). Этапы приземления 6 и 7 неуправляемые.
Управляющими воздействиями являются моменты М43, М21, создаваемые приводами 5 и 6, посредством которых соответствующие звенья робота
поворачиваются на углы ф2 и ф1 , и сила 32, формируемая приводом 7, за счет которой звенья 2 и 3 поступательно перемещаются друг относительно друга.
Будем рассматривать указанные величины как кусочно-постоянные функции,
численные значения которых на каждом из этапов прыжка являются постоянными.
Диаграмма логических уровней управляющих параметров представлена на рис. 3,0-
отсутствие сигнала, ±1 - наличие сигнала.
Зтап | прыжка I
щ 1
о
М2!\
1
Б
-1
6.?!
1
О
-1
к ! I
к I и
Рис. 3 Циклограмма управляющих моментов и сил
Высота и длина прыжка варьируются за счет изменения временных интервалов ^
о
- t5 этапов 1-5, что позволяет формировать различные значения угла ф2, под
7°
которым происходит разгон устройства, скорости 23 отрыва объекта от поверхности, а также скорости и ускорения относительного перемещения звеньев 2 и 3 во время полета.
Моделирование движения робота
В соответствии с последовательностью этапов движения прыгающего робота на геометрические размеры звеньев наложены следующие ограничения:
/£ш +1Л<Ь /23Ш + 1л<а и<Ь 12 < а 13<Ь 13< а 23 эшср® ^
При численном моделировании прыжка устройства использовался переход к безразмерным параметрам при помощи масштабных коэффициентов М=0.05 кг, 7=0.1 с, /.=0.1 м. Объектом
моделирования является прыгающий робот, обладающий следующими массогабаритными параметрами: т\=\, т2=\, гпу=\, т4=1, <а?=1, Ь=1, /1=0.5,4=0.9, /3=0.9, /23тш=0.4, /2з(Н.8, = 4;Ф2
=7г/4. Моделирование проводится при начальных условиях: {= 0 , х4 = а ? х4 = 0 5 У4 = ь 5 = 0 ?
ср4= 0 ср4 = 0 Фх =п ф1 = 0 ф2=7г/2 ср2 = 0
? ? ? ? ?
На рис. 4-7 приведены временные зависимости координат, посредством которых описывается движение прыгающего робота, а также показаны этапы одного прыжка. По данным графикам видно, что моделируемое движение робота соответствует разработанной и описанной последовательности этапов.
Рис. 4 Графики зависимостей: а б - у\(0
*4 У4
Рис. 5 Графики зависимостей: а -х4(/% б -^(0
На первых трех этапах, во время которых происходит позиционирование робота, центр масс корпуса неподвижен (х4=соп81:, ^4=соп81:, ф4 =соп81:). Звено 1 на первом
этапе поворачивается вместе со звеном 2 (Лф1 = 4ф2) до тех пор, пока не станет
_ о о
справедливым условие ф2 _ф2 (в рассматриваемом случае ф2 =л/4), что также
отражается на временных зависимостях координат центра масс стопы.
Во время второго этапа происходит поворот звена 1 против часовой стрелки до
достижения значения ф| =тг, траектория у\(х\) на этом этапе имеет форму дуги окружности.
т
Фг
5 Г
Рис. 6 Графики зависимостей: а - ф1 (0, б - ф2 (?)
Ср,,
1
0.5
-0.5
-1
0 2 4 6
а
1 2 ! 3 1 5 6
1
1 !
! 1
1.5
0.5
1\
1 3 1 2 1 1 5 6
Рис. 7 Графики зависимостей: а - ф4 (0, б - /2з(0
На третьем этапе наблюдается выдвижение звеньев 1 и 2 до взаимодействия звена
1 с поверхностью, что находит отражение на графиках рис. 4, 6а, 76, причем
изменение координат точки С\ равно ^ = А/гз C0S(P2, A^i = А4з sin(P2 ^ углы наклона звеньев к оси Ох не меняются. Завершению третьего этапа соответствует выполнение условия >>1=0, это означает, что звено 1 лежит на поверхности. Изменение величины /2з под действием постоянной силы F32 происходит по ветви параболы, что видно по графику рис. 76, характер изменения Х\ и у{ аналогичен.
Четвертый этап происходит при взаимодействии звена 1 с поверхностью, что подтверждается графиками рис. 4 (xi=const, >'i=0) и характеризуется перемещением звеньев 3 и 4 относительно звена 2 в положительных направлениях осей Ох и Оу, т.е. увеличением расстояния /2з. На этом этапе наблюдается возрастание координат
центра масс корпуса, причем справедливы следующие соотношения Ах'4 = А/23 C0S(P2 5
Ду4 = Д/23 sincp2 увеличение ДЛИНЫ /23 происходит ПОД Действием ПОСТОЯННОЙ СИЛЫ /*23, значит, характер изменения этой обобщенной координаты так же, как и значений х4 и у4, параболический. Поворота звеньев на данном этапе не происходит, что видно по графикам рис. 6, 7а.
Пятый этап, во время которого робот движется с отрывом от поверхности, характеризуется изменением координат центров масс всех звеньев робота, а также втягиванием звеньев 1 и 2 в корпус. В моделируемом прыжке линии действия разгоняющей силы и силы, посредством которой звенья 1 и 2 втягиваются в корпус, проходят через центр масс последнего, поэтому вращения корпуса относительно его центра масс в полете не происходит, что находит отражение на графике рис. 76. Звенья 1 и 2 перемещаются относительно корпуса до тех пор, пока значение /2з не
7 _ /min
станет равным начальному ( 23 “ 23 ), это наблюдается до момента времени, отмеченного на всех графиках пунктирной линией. На этом участке пятого этапа (до пунктирной линии) изменения координат центров масс звеньев 1 и 4 различны и Axi>Ax4, Ау1>Ау4, т.к. скорость звена 1 больше скорости звена 4 и в проекциях на оси
Ох и Оу равна xi = х4 + къ C0S(P2 5 У\ = Уа + къ sin(P2 После того, как втягивание ноги завершится, все звенья робота будут двигаться с одинаковыми скоростями и
ускорениями, поэтому изменения координат точек С\ и С4 будут равными: Лх1=Лх4, Д>'1=>'4. На графиках у\{1) в момент времени, когда прекращается втягивание ноги, виден излом, до которого подъем графиков более крутой, а после - более пологий, что подтверждает правильность моделирования.
Движение робота переходит в шестой этап, когда одна или две точки корпуса начнут взаимодействовать с опорной поверхностью. В связи с тем, что в рассматриваемом случае поворота корпуса робота в полете не происходит, приземление будет осуществляться на те же точки корпуса, с каких начинался первый этап прыжка, т.е. на точки А\ и А4. На графике у4(1) видно, что в момент начала этапа 6 значение ординаты центра масс корпуса равно значению у4 в момент начала прыжка, это означает что точки А\ и А4 лежат на поверхности. Изменения координат центров масс звеньев 1 и 4 во время этого этапа не наблюдается, т.к. в соответствии с принятой моделью трения скорости точек контакта с поверхностью равны нулю. Углы поворота звеньев также остаются постоянными.
Анализ представленных в работе графиков свидетельствует о правильности разработанных алгоритма и программы моделирования прыжка робота, что подтверждается соответствием характера изменения координат вектора ц описанной последовательности этапов движения устройства.
Заключение
В статье рассматривается прыгающий робот, один прыжок которого представлен в виде последовательности восьми этапов. Разработана математическая модель поэтапного движения робота, условия перехода из одного этапа в другой, а также связи, наложенные на движения звеньев робота на каждом из этапов.
В соответствии с данной математической моделью проведено численное моделирование прыжка робота, результаты которого позволяют подтвердить правильность разработанной программы и проследить последовательность этапов движения на временных зависимостях характеристик робота.
Список литературы
1. Zufferey J.-C. First Jumps of the 3D Bow Leg Hopper // Micro-engineering, EPFL. 2001. 65 p.
2. Ishikawa A. I. М., Neki М., Нага J. Energy preserving control of a hopping robot based on hybrid port-controlled hamiltonian modeling // ШЕЕ Conference on Control Applications, Yokohama, Japan. 2003. Vol. 2. P. 1136-1141.
3. Zeglin G. Uniroo: A One-Legged Dynamic Hopping Robot // Massachusetts Institute of Technology. 1991. 65 p.
4. Sato A. A Planar Hopping Robot with One Actuator: Design, Simulation, and Experimental Results // Department of Mechanical Engineering, Montreal, Canada. 2004. 106 p.
5. Kovac М., Fuchs М., Guignard A. et al. A miniature 7g jumping robot // Proc. of the IEEE Intern. Conf. on Robotics and Automation. Pasadena. 2008. P. 373-378.
6. Волкова Л.Ю., Серебровский В.В., Яцун С.Ф. Исследование движения многозвенного робота, перемещающегося с отрывом от поверхности // Известия Юго-Западного государственного университета. 2012. № 1 (40). Ч. 2. С. 12-18.
7. Яцун С.Ф., Рукавицын А.Н., Волкова Л.Ю. Исследование динамики управляемого прыжка робота // Известия Самарского научного центра Российской академии наук. 2012. Т. 14. №4(5). С. 1355 - 1358.
8. Яцун С.Ф., Волкова Л.Ю. Управление высотой и длиной прыжка робота путем его позиционирования и разгона // Известия ЮЗГУ. Серия Техника и технологии. 2012. № 2. Ч. 1 С. 210-213.
9. Яцун С.Ф., Волкова Л.Ю. Исследование влияния закона втягивания ноги в полете на высоту и длину прыжка многозвенного робота // Materialy IX mezinarodni vedecko - prakticka conference «Moderni vymozenosti vedy - 2013». Praha, Str. 77-81.
10. Лушников Б.В. Компьютерное моделирование динамики элемента сухого некулонова трения // Известия Самарского научного центра Российской академии наук. 2010. Т. 12. №1(2). С.439-444.
11. Лушников Б.В. Испытательный стенд для идентификации динамических параметров сухого некулонова трения// Фундаментальные исследования. 2012. №3. Ч.З. С. 238-241.
12. Болотник Н. Н, Зейдис И.М., Циммерман К. и др. Динамика управляемых движений вибрационных систем//Изв. РАН.ТиСУ. 2006. №5. С. 157-167.
References
1. Zufferey J.-C. First Jumps of the 3D Bow Leg Hopper // Micro-engineering, EPFL. 2001. 65 r.
2. Ishikawa A. I. М., Neki М., Нага J. Energy preserving control of a hopping robot based on hybrid port-controlled hamiltonian modeling // IEEE Conference on Control Applications, Yokohama, Japan. 2003. Vol. 2. P. 1136-1141.
3. Zeglin G. Uniroo: A One-Legged Dynamic Hopping Robot // Massachusetts Institute of Technology. 1991. 65 r.
4. Sato A. A Planar Hopping Robot with One Actuator: Design, Simulation, and Experimental Results // Department of Mechanical Engineering, Montreal, Canada. 2004. 106 r.
5. Kovac М., Fuchs М., Guignard A. et al. A miniature 7g jumping robot // Proc. of the IEEE Intern. Conf. on Robotics and Automation. Pasadena. 2008. P. 373-378.
6. Volkova L.Ju., Serebrovskij V.V., Jacun S.F. Issledovanie dvizhenija mnogozvennogo robota, peremeshhajushhegosja s otryvom ot poverhnosti // Izvestija Jugo-Zapadnogo gosudarstvennogo universiteta. 2012. № 1 (40). Ch. 2. S. 12-18.
7. Jacun S.F., Rukavicyn A.N., Volkova L.Ju. Issledovanie dinamiki upravljaemogo pryzhka robota // Izvestija Samarskogo nauchnogo centra Rossijskoj akademii nauk. 2012. T. 14. №4(5). S. 1355 - 1358.
8. Jacun S.F., Volkova L.Ju. Upravlenie vysotoj i dlinoj pryzhka robota putem ego pozicionirovanija i
razgona // Izvestija JuZGU. Serija Tehnika i tehnologii. 2012. № 2. Ch. 1 S. 210-213.
9. Jacun S.F., Volkova L.Ju. Issledovanie vlijanija zakona vtjagivanija nogi v polete na vysotu i dlinu pryzhka mnogozvennogo robota // Materialy IX mezinarodni vedecko - prakticka conference «Moderni vymozenosti vedy - 2013». Praha, Str. 77-81.
10. Lushnikov B.V. Komp'juternoe modelirovanie dinamiki jelementa suhogo nekulonova trenija // Izvestija Samarskogo nauchnogo centra Rossijskoj akademii nauk. 2010. T.12. №1(2). S.439-444.
11. Lushnikov B.V. Ispytatel'nyj stend dlja identifikacii dinamicheskih parametrov suhogo nekulonova trenija//Fundamental'nye issledovanija. 2012. №3. Ch.3. S. 238-241.
12. Bolotnik N. N., Zejdis I.M., Cimmerman K. i dr. Dinamika upravljaemyh dvizhenij vibracionnyh sistem // Izv. RAN.TiSU. 2006. №5. S. 157-167.
