Модели механики в задаче управления силовым взаимодействием робота с поверхностью неопределенного профиля

DOI: 10.5862/JCSTCS.224.7 УДК 681.51

С.Ф. Бурдаков, О.Б. Шагниев

модели механики в задаче управления силовым взаимодействием робота с поверхностью неопределенного профиля

S.F. Burdakov, O.B. Shagniev

mechanics models in the control problem of the force interaction between a robot and a free-formed surface

На примере задачи управления силовым взаимодействием очувствленного робота с поверхностью неопределенного профиля показана решающая роль моделей механики при стремлении получить высокое качество отработки желаемых движений. Адекватные модели механики необходимы как при согласовании инерционных и жесткостных характеристик робота и силового датчика, так и при выборе структуры и настроек системы позиционно-силового управления, а также при выборе реализуемых желаемых движений. Рассмотрены режим позиционирования с выходом робота на контакт с поверхностью и режим управления силовым взаимодействием робота с поверхностью. Исследованы алгоритмы позиционно-силового управления с ПИД-регулятором и интегральным компенсатором. Показана возможность дополнительного повышения качества управления с помощью коррекции желаемых движений методом обучения.

РОБОТ; ПОВЕРХНОСТЬ НЕОПРЕДЕЛЕННОГО ПРОФИЛЯ; СИЛОВОЕ ОЧУВСТВЛЕНИЕ; МОДЕЛИ МЕХАНИКИ; ПОЗИЦИОННО-СИЛОВОЕ УПРАВЛЕНИЕ; ОБУЧЕНИЕ.

The example of the problem of controlling the force interaction of a robot with a free-formed surface shows the crucial role of mechanics models in obtaining high-quality implementations of the desired motions. Adequate mechanics models are needed both for coordinating the inertial and stiffness characteristics of the robot and the force sensor, and for choosing the structure and configuration of the position-force control, as well as when selecting the implemented desired motions. The article describes a positioning mode with the release of the robot to contact the surface and a force control mode of the robot interacting with surface. We investigated the algorithms of position-force control with a PID-controller and an integral compensator. The article shows the possibility of further enhancing the quality of control by correcting the required motion by learning.

ROBOT; FREE-FORMED SURFACE; FORCE SENSING; MECHANICS MODELS; POSITIONFORCE CONTROL; LEARNING.

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

еще остается достаточно широким. Прежде всего это относится к роботам для экстремальных условий, для целого ряда технических задач (обработка деталей, сборка и т. п.), а также к медицинским роботам [2—4].

Фундаментальные аспекты использования информации об усилиях при управлении роботами рассмотрены в работе [5]. Наиболее распространенный на практике метод связанного (гибридного) позиционно-

силового управления, когда управление перемещениями происходит в функции от силы взаимодействия, рассмотрен в [6]. На основе этого метода в [7] предложен алгоритм связанного позиционно-силового управления с сохранением базового контура управления по положению. Этот алгоритм позволил после некоторой модернизации (установки датчика сил и моментов) использовать серийные роботы для решения задач управления силомоментным взаимодействием роботов с внешними объектами. Это не только существенно расширяет области применения серийных роботов [3, 8], но и приводит к появлению новых областей, например, в медицине [9]. Однако при подобном очувствлении возникает проблема правильного согласования механических характеристик робота и силомо-ментного датчика с настройками системы управления и выбором желаемых движений. Эта проблема не может быть решена без адекватных моделей механики.

Вопрос о создании робота-массажера с использованием информации об усилиях взаимодействия впервые был поставлен российскими учеными на Втором симпозиуме по медицинской робототехнике в Гей-дельберге в 1997 г. В работах [4, 9] представлена реализация робота-массажера, в алгоритме управления которого важную роль играет механическая модель взаимодействия робота с мягкими тканями и суставами.

Результаты научных исследований и примеры практического применения робо-

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

Постановка задачи

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

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

Рис. 1. Расчетные схемы робота в режиме без контакта

является гладким, но заранее не определен. Поверхность не деформируема.

На рис. 1 введены следующие обозначения: т1 и т2 — приведенная масса двигателя и масса руки робота; с — эквивалентная жесткость робота; у1 и у' — координаты двигателя и руки робота; т, и с — масса и жесткость датчика силы; у' — координата датчика силы; уур (х) — координата поверхности.

В режиме позиционирования (рис. 1 а) математическая модель очувствленного робота как объекта управления имеет вид:

т^ + ¿1 у + Ь(у - У') + + с(У1 - У') = Р - т^, т2У2 + Ь2у2 + Ь(У2 - Ю + С(У2 - У1) + (1) + ¿(у2 - Ю + С(У2 - У,) = -т2§,

ту , +Ь, (У - :у2) +с (у,- У2) =

где Р — сила, развиваемая двигателем робота; т^, m2g, т^ — силы тяжести; Ь, Ь1, Ь2, Ь — коэффициенты вязкого трения.

В отличие от координат у' и у' (рис. 1 а), координаты у2 и у, в уравнениях (1) учитывают конструктивные размеры робота и датчика силы, а также статические удлинения упругих элементов с и с,.

Переход к позиционно-силовому режиму управления происходит при безотрывном контакте датчика силы с поверхностью:

У, = у', Р ^ р,

где Р = Ь(у, - у2) + с,(у* - у22 — показание

датчика силы (Р, > 0 соответствует сжатию упругого элемента датчика); Р,, — величина, характеризующая зону переходных процессов по силе Р,, при переходе от режима позиционирования к позиционно-силовому режиму; у '()) = у^ (х )| х=х" ()) — не определенная заранее координата поверхности при условии, что движение по координате х происходит по желаемому закону х"()) (рис.1 б).

Математическая модель робота в позиционно-силовом режиме имеет вид:

щух + ¿1 у + Ь(у - у2) + + с( у1 - у2) = р - т g, (3) т2у2 + Ь2у2 + Ь(у2 - у1) + с(У2 - У1) +

+ Ь> (у2 - у, ) + с, (у2 - У, ) = -m2g, у> = у'()).

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

Будем также считать, что в режиме позиционирования обратная связь формируется только по измерениям датчика положения, а датчик силы используется в индикаторном режиме для определения момента времени перехода к позиционно-силовому режиму. В позиционно-силовом режиме обратная связь формируется по измерениям обоих датчиков.

Задающим воздействием для режима позиционирования является у" (?), а для позиционно-силового режима — Р ,' " (?). Возмущающим воздействием для позиционно-силового режима является у '(I). Дополнительно вводится требование о недопустимости нарушения контакта робота с поверхностью после его установления.

В дальнейшем при компьютерном моделировании будут использоваться следующие значения параметров робота:

т1 = т2 = 1 кг, т,, = 0,2 кг, с = 2,2  105Н/м,

с, = 104 Н/м, Ь = 10 Нс/ м,

Ь1 = Ъ2 = 20 Нс/м, = 5 Нс/м.

По этим параметрам могут быть получены оценки механических характеристик очувствленного робота:

• для режима позиционирования

т1 + т2 + т,

;—2-- - (0,05 + 0,1) с — меха-

Т -

1ш —

Ч -Г и2

ническая постоянная времени,

л>2

т1 + т2

т.

- (30 35) Гц

* /

первая собственная частота,

т т + т, собственная частота;

- 100 Гц — вторая

для позиционно-силового режима

т1 + т2

- (10 ^ 12) Гц — первая

собственная частота,

(

^ 12

= 100 Гц — вторая соб-

— + —

V т т у ственная частота.

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

Режим выхода робота на контакт с поверхностью

В рассматриваемом режиме объектом управления является многомассовая система (см. рис. 1 а). Требуется перевести ее из начального положения у10 в сторону поверхности у*( )) до контакта с ней силового датчика. Для решения этой задачи при одном датчике обратной связи, измеряющем координату у1, можно использовать ПИД-регулятор, т. е. сформировать управляющую силу Е по закону:

Г

Е =

К + к, - + ко М , Р 1 + N1

\

где кр, к,, к0 — коэффициенты обратных связей; N — параметр дифференцирующего фильтра.

Замкнутая система (1), (4) обладает астатизмом второго порядка, поэтому задание для режима позиционирования целесообразно выбрать в следующем виде:

Уа

где у0 = у10 — начальное значение задания, согласованное с начальным положением робота; ¥у — варьируемая скорость убывания задания.

Настройка в МаНаЪ ПИД-регулятора на максимум полосы рабочих частот замкнутой системы при принятых выше значениях параметров объекта управления дает кр = 9,58-103 Н/м; к, = 1,63 • 104 Н/мс; к0 = 4,71 • 102 Нс/м; ^ 3,69 • 102 1/с.

Выбранная настройка ПИД-регулятора в режиме позиционирования согласована с механическими характеристиками очувствленного робота так, что ю с « 10ТМ, <Ю <юс < ю2, где юс — частота среза частотной характеристики замкнутой системы от входа у0 до выхода у1.

Начало отсчета для координаты поверхности 0 < у) < у0 соответствует рис. 1 а. Учитывая, что координата поверхности у*()) заранее не определена, выход робота на контакт с поверхностью может произойти в любой момент времени ) *

0 <) * < у0 V.

Задание (5) с ограниченной скоростью убывания V позволяет установить при-

а)

у, м [19 0.3 07 ОБ 0.5 0.4

о.э 0.2 01

-Координата робота ---Поверхность контакта

б)

), с

7

), с

Рис. 2. Графики процессов у ()) (а) и Е ()) (б) в режиме позиционирования

Е, Н

емлемыи уровень

Fs = sup Fs показаний

ie[0,t"]

датчика силы, который можно принять за индикатор установления безотрывного контакта робота с поверхностью. Заметим, что в режиме позиционирования датчик силы работает только как индикатор перехода к позиционно-силовому режиму управления. Уровень F > 0 устанавливается путем компьютерного моделирования или экспериментально.

На рис. 2 приведены графики процессов ys(t) и Fs (t) в режиме позиционирования

пРи Уо = У10 = У20 = Ло = 1 м и vy = 0>2 м/с-Поверхность постоянного профиля находилась на уровне y* = 0,1 м. Видно, что контакт робота с поверхностью происходит через t* = 4,5 с. По переходным процессам при t < t находим F = sup Fs = 5 Н.

te[0,t']

При t > t упругая сила Fs = cs (y s - y2) быстро нарастает, т. к. при контакте ys = y * = 0,1 = const, а координата y2 про-

должает уменьшается в соответствии с заданием (5).

Чтобы предотвратить поломку робота или повреждение поверхности необходимо перейти при ^ > )? к управлению силовым взаимодействием робота с поверхностью.

Режим позиционно-силового управления

В соответствии с концепцией связанного позиционно-силового управления [6, 7] контур обратной связи по положению (4) сохраняется без изменения. Однако вместо задания у" при выполнении условия (2) автоматически вводится новое задание у", которое операторным образом зависит от измерений датчика силы Р :

=W(p)(F/ - Fs),

где Р " — задание по силе взаимодействия с поверхностью, которое должен отрабатывать робот; Т(р) — передаточная функция регулятора силы.

Рис. 3. Частотные характеристики базового контура от задания у" до показаний датчика силы Р

Таким образом, в режиме позиционно-силового управления вводится еще одна обратная связь. При этом базовый контур с обратной связью по положению (3), (4) является для нее объектом управления. Базовый контур положения из-за контакта робота с поверхностью существенно меняет свои динамические свойства по сравнению с тем, что было в режиме позиционирования. На рис. 3 приведены частотные характеристики базового контура от задания у0 до показаний датчика силы Е5 . Видно, что в диапазоне рабочих частот 0 < ю < юс, который был достигнут при позиционировании, из-за контакта робота с поверхностью происходят изменения частотных характеристик, требующие существенной коррекции замкнутой системы в сторону расширения диапазона рабочих частот. Если пренебречь высокочастотной динамикой, обусловленной упругой податливостью робота, то соответствующая передаточная функция базового контура будет иметь вид:

5 (Р) =

с5

тръ + (Ь + ка)р2 + (С5 + к )р + к

.

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

р = —V = -0,83 1/с.

С5 + Кр

Структурная схема замкнутой системы (6), (7) приведена на рис. 4.

Рассмотрим некоторые способы реализации обратной связи по силе [5, 7, 11].

Управление при помощи второго ПИД-регулятора

л 1 N

ж(Р) = + - + ^-

р 1 + N

1

Настройка в МаНаЪ регулятора (8) при принятых выше значениях параметров на максимальный диапазон рабочих частот замкнутой системы (рис. 4) дает:

wp = 3,93 • 10-5 м/Н, = 1,24 • 10-2 м/сН, =-1,11 • 10-7мс/Н, N = 3, 55 -102 1/с.

Управление при помощи интегрального

компенсатора

ж (р) = К^-1( р),

Р

где к — коэффициент усиления; 5* (р) — физически реализуемое приближение оператора _1( р), обратного к 5(р).

Идея этого метода состоит в том, что при полной компенсации динамики объекта управления, т. е. при 5 (р) = 5_1(р), передаточная функция замкнутой системы (рис. 4) приобретает вид:

Н (р) =

Ж (р)Б (р) = к 1 + Ж (р)Б (р) р + к'

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

(С5 + Кр ) р + К

5 (р) =

с (Кр2 + к»р + К У

Приближенная компенсация не допускает неограниченного роста коэффициента усиления к из-за потери устойчивости замкнутой системы (6), (7). Диапазон приемлемых значений при принятых выше значениях параметров составляет к = (30 ^ 100) 1/с.

Замкнутая система (6), (7) обладает свойством астатизма, поэтому целесообразно выбрать задание Е/ по силе Е5 с варьируемой скоростью ¥е :

Е/ = Е + V,а - /*), / е [/*, /"],

где ) =) +

Е - Е

V,

время выхода зада-

ния Е 0 на требуемый уровень Е*.

На рис. 5 приведены графики переходных процессов выхода робота на контакт с поверхностью у * = 0,1 м и установления силового взаимодействия робота с поверхностью на уровне Е* = 40 Н. Рисунок 5 а соответствует варианту со вторым ПИД-

Рис. 4. Структурная схема замкнутой системы в позиционно-силовом режиме управления

Рис. 5. Графики переходных процессов при использовании второго ПИД-регулятора

и интегрального компенсатора

регулятором. Рисунок 5 б соответствует варианту с интегральным компенсатором.

Оба варианта дают вполне приемлемые результаты. Сплошные линии на графиках Р \ (^) соответствуют ступенча-

тому заданию. Пунктирные линии соответствуют заданию с конечной скоростью ¥Р = 40 Н/с. Плавное задание почти полностью снимает перерегулирование по силе Р5.

а)

ур м 0.3

Е, > Н

ДО

Е, Н ео

02 0 3 0 4 05 ОБ 07 О.Э 03 X, м

б) Л, м 03

0 25

0.15

0 05

Е, Н

л'

10

Э5 30 25 20 15 10 5 О

02 0.3 04 05 06 0.7 0.8 О.Э

X, м

01 02 0.3 0.4 0.5 0.6 0 7 0.8 О.Э

X, м

Е, Н 60 40 20 О -20 -40 60

-т-1-1-1-р-Г"

V

0.1 0.2 О.Э 04 05 0.6 0.7 О.Э 0.9

X, м

Рис. 6. Графики координаты робота, показаний датчика силы и управляющей силы

без обучения

Обучение роботов перед выполнением циклических операций

Для очувствленных роботов, выполняющих циклические операции, имеется дополнительный ресурс повышения точности поддержания требуемого уровня Р * силы Р с помощью процедуры обучения на пробных циклах движения [12, 13].

На (рис. 6 а) приведены графики вертикальной координаты робота у1{)?, показаний датчика силы Р\ (^) и управляющей силы Р()) в случае гармонического профиля у*()) = 0,1 + 0, 05 8т10х| = " (Хорошо

Iх=х {?)

видно, что ошибка отработки задания по силе взаимодействия еР = Р5 - Р* тоже меняется по гармоническому закону с амплитудой порядка 2 Н. Такая точность может оказаться недостаточной. В целях повышения точности будем использовать итеративную процедуру обучения с низкочастотным фильтром [13]:

Г+1) = 1С * к =1,2,...,

где у < 1 — настраиваемый весовой коэффициент.

В соответствии с (9) для каждого (к+1)-го цикла движения коррекция задания производится после завершения к-го цикла по сохраненной в памяти системы управления ошибке вР.

На рис. 6 б приведены графики про-

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

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

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

1. Юревич Е.И. Интеллектуальные роботы. М.: Машиностроение, 2007. 360 с.

2. Попов А.В., Юревич Е.И. Позиционно-силовое управление манипуляторами: состояние и перспективы // Мехатроника, автоматизация и управление. 2008. № 5. Приложение. С. 2—6.

3. Егоров И.Н. Позиционно-силовое управление робототехническими и мехатронными устройствами. Владимир: Изд-во Владимир. гос. ун-та, 2010. 191 с.

4. Головин А.ф., Архипов М.В., Журавлев В.В. Проблемы развития робототехники в восстановительной медицине // Мехатроника, автоматизация, управление. 2009. № 9. С. 49—53.

5. Гориневский Д.М., формальский А.М., Шнейдер А.Ю. Управление манипуляционными системами на основании информации об усилиях. М.: Физматлит, 1994. 368 с.

6. Hogan N. Impedance control: an approach to manipulation // ASME J. Dynamic Systems, Measr. and Control. 1985. Vol. 107. Pp. 1-24.

7. Lawrence D., Stronghton R.M. Position-based impedance control: Achieving stability in practice // Proc. AIAA Cont. on Guidance, Navigation and Control. 1997. Pp. 221-226.

8. Обработка турбинной лопатки [электронный ресурс]/ URL: www.youtube.com.

9. Головин А.ф., Архипов М.В., Журавлев В.В. Метод силового обучения при планировании траекторий робота для восстановительной медицины // Мехатроника, автоматизация, управление. 2009. № 10. С. 29-30.

10. An C., hollerbach J. The role of dynamic models in cartesian force control of manipulators // Int. J. of Robotics Res. 1989. Vol. 8(4). Pp. 51-72.

11. Rocco P., Ferretti G., Magnani G. Implicit force control for industrial robots in contact with

stiff surfaces // Proc. of the 13th Triennial World Congress. San Francisco, 1996. Pp. 355—360.

12. Первозванский А.А. Обучаемое управление и его приложения // Автоматика и телемеханика. 1995. № 11, 12.

13. Бурдаков С.ф., Смирнова Н.А. Обучаемое управление упругим манипулятором при выполнении циклических операций // Изв. РАН. Теория и системы управления. 2002. № 4. С. 142-149.

references

1. Yurevich Ye.I. Intellektualnyye roboty [Intelligent robots]. Moscow: Mashinostroyeniye Publ., 2007, 360 p.

2. Popov A.V., Yurevich Ye.I. Pozitsionno-silovoye upravleniye manipulyatorami: sostoyaniye i perspektivy [Position-force control of manipulators: Status and Prospects]. Mekhatronika, avtomatizatsiya i upravleniye [Mechatronics, automation and control], 2008, No. 5, Pp. 2-6.

3. Yegorov I.N. Pozitsionno-silovoye upravleniye robototekhnicheskimi i mekhatronnymi ustroystvami. [Position-force control of robotic and mechatronic devices], Vladimir: VGU Publ., 2010, 191 p.

4. Golovin A.F., Arkhipov M.V., Zhuravlev V.V. Problemy razvitiya robototekhniki v vosstanovitelnoy meditsine [Problems of development of robotics in rehabilitation medicine]. Mekhatronika, avtomatizatsiya, upravleniye [Mechatronics, automation and control], 2009, No. 9, Pp. 49-53.

5. Gorinevskiy D.M., formalskiy A.M., Shneyder A.Yu. Upravleniye manipulyatsionnymi sistemami na osnovanii informatsii ob usiliyakh [Office handling systems on the basis of information on efforts]. Moscow: Fizmatlit Publ., 1994, 368 p.

6. Hogan N. Impedance control: an approach to manipulation. ASME J. Dynamic Systems, Measr. and Control, 1985, Vol. 107, Pp. 1-24.

7. Lawrence D., Stronghton R.M. Position-based impedance control: Achieving stability in

practice. Proc. AIAA Cont. on Guidance, Navigation and Control, 1997, Pp. 221-226.

8. Obrabotka turbinnoy lopatki [Processing turbine blade]. Available: www.youtube.com

9. Golovin A.F., Arkhipov M.V., Zhuravlev V.V. Metod silovogo obucheniya pri planirovanii trayektoriy robota dlya vosstanovitelnoy meditsiny [The method of power training in the planning of robot trajectories for regenerative medicine]. Mekhatronika, avtomatizatsiya, upravleniye [Mechatronics, automation and control], 2009, No. 10, Pp. 29-30.

10. An C., Hollerbach J. The role of dynamic models in cartesian force control of manipulators, Int. J. of Robotics Res., 1989, Vol. 8(4), Pp. 51-72.

11. Rocco P., ferretti G., Magnani G. Implicit force control for industrial robots in contact with stiff surfaces. Proc. of the 13th Triennial World Congress, San Francisco, 1996, Pp. 355-360.

12. Pervozvanskiy A.A. Obuchayemoye upravleniye i yego prilozheniya [Educational management and its application]. Automatica i telemechanica [Automatics and telemechanics], 1995, No. 11, 12.

13. Burdakov S.F., Smirnova N.A. Obuchayemoye upravleniye uprugim manipulyatorom pri vypolnenii tsiklicheskikh operatsiy [Educational Management resilient arm when the cyclic operations]. Izv. RAN. Teoriya i sistemy upravleniya [Izvestia RAS. Theory and control systems], 2002, No. 4, Pp. 142-149.

БУРДАКОВ Сергей федорович — профессор кафедры механики и процессов управления Института прикладной математики и механики Санкт-Петербургского политехнического университета Петра Великого, доктор технических наук.

195251, Россия, Санкт-Петербург, ул. Политехническая, д. 29.

E-mail: burdakov.s@yandex.ru

BURDAKOV Sergey F. Peter the Great St. Petersburg Polytechnic University.

195251, Politekhnicheskaya Str. 29, St. Petersburg, Russia.

E-mail: burdakov.s@yandex.ru

ШАГНИЕВ Олег Булатович — инженер 2 категории кафедры механики и процессов управления Института прикладной математики и механики Санкт-Петербургского политехнического университета Петра Великого.

195251, Россия, Санкт-Петербург, ул. Политехническая, д. 29.

E-mail: shagnoleg@yandex.ru

SHAGNIEV Oleg B. Peter the Great St. PetersburgPolytechnic University. 195251, Politekhnicheskaya Str. 29, St. Petersburg, Russia. E-mail: shagnoleg@yandex.ru

© Санкт-Петербургский политехнический университет Петра Великого, 2015