Кукрск, Россия
Курск, Курская область, Россия
Курск, Курская область, Россия
Цель исследования состоит в разработке конструкции и исследовании системы управления колёсного подводного робота с регулируемой плавучестью для подлёдного мониторинга подводных объектов. В рамках исследования ставится задача разработки математической модели и её компьютерной реализации для управляемого движения колёсного робота по внутренней поверхности льда, а также постановка вычислительных экспериментов, демонстрирующих работоспособность предлагаемой системы автоматического управления антипробуксовочного устройства. В работе используются как классические методы построения и исследования математических моделей и систем управления мобильных роботов, так и оригинальные методики численного интегрирования систем нелинейных дифференциальных уравнений, описывающих эффекты проскальзывания колёс робота относительно опорной поверхности льда. Новизна работы заключается в использовании оригинального способа перемещения колёсного робота на границе сред с регулированием нормальной реакции за счёт управления объёмом колёсных воздушных баллонов и контролем проскальзывания с помощью системы автоматического управления, обеспечивающей повышение энергоэффективности работы приводов в процессе разгона. В ходе исследования были разработаны: оригинальная конструкция и способ перемещения робота для мониторинга подводных объектов, расчётная схема устройства и структурная схема системы управления приводной системы робота. Создана математическая модель, описывающая управляемое движение объекта и выполнен комплекс вычислительных экспериментов, доказывающих применимость предлагаемых методов реализации движения без пробуксовки ведущих колёс мобильного робота. Предлагаемый в работе способ перемещения имеет ряд ключевых преимуществ при осуществлении мониторинга подводных объектов, однако реализация преимуществ на практике сопряжена с необходимостью использования системы контроля проскальзывания колёс робота с опорной поверхностью. В работе предложена структура и алгоритмы системы управления, позволяющие обеспечить движение с минимальным проскальзыванием, а вычислительные эксперименты на математической модели устройства доказывают её применимость и эффективность.
колёсный подводный робот, подводный мониторинг, контроль проскальзывания, система управления, математическое моделирование
Введение
Известно, что немалая часть объектов хозяйственной инфраструктуры сегодня находится под водой. Это кабели связи и интернета, трубопроводы, опоры мостов и причалов, подводные части гидросооружений, таких как причалы, плотины, дамбы и т.д. [1 – 2] Поддержание эксплуатационных свойств таких объектов, предотвращение возможных аварий требует выполнения регулярных осмотровых работ и мониторинга за их состоянием [3 – 4]. Однако в Арктической зоне, а также в зимние месяцы многие водоёмы покрыты льдом, что делает невозможным выполнение данных работ без погружения под воду, а привлечение водолазов сопряжено с повышенными рисками для их жизни и дополнительными расходами [5]. Для обеспечения мониторинга подводных объектов весьма рациональным является использование мобильных роботов, оснащенных необходимым оборудованием и способных выполнять требуемые операции [6 – 12]. Наличие льда на поверхности водоёма, открывает возможности использование границы раздела сред для скоростного и энергоэффективного перемещения за счёт гусеничных или колесных движителей [13 – 15]. Перемещение робота по внутренней поверхности льда имеет преимущество при выполнении видеосъёмки или других типов диагностики за счёт надежного позиционирования в неподвижном состоянии без потребления дополнительной энергии на компенсацию внешних воздействий – течений, волнения и т.д. Кроме того, подводный робот, способный эффективно перемещаться по внутренней поверхности льда, может найти применение для поисковых задач и задач сопровождения, например, ледоколов при прохождении участков с возможной подлёдной опасностью [14].
Предлагаемое в работе роботизированное устройство способно изменять собственную плавучесть для обеспечения управления нормальными реакциями в точках контакта колёс с поверхностью льда. Это открывает дополнительные возможности по преодолению препятствий – трещин, выступов, полыньи и т.д. Вместе с тем перемещение по поверхности льда требует контроля процесса взаимодействия колёс со льдом, предотвращения проскальзывания для обеспечения энергоэффективного перемещения.
В работе ставится задача создания математической модели и системы управления мобильного подлёдного колёсного робота, позволяющих установить закономерности влияния различных параметров устройства на кинематические и динамические характеристики движения подо льдом, определить их рациональные значения и настройки системы управления колесным приводом.
Материалы, модели, эксперименты и методы
Рассмотрим более подробно конструкцию и принцип управляемого движения подводного робота для подлёдного перемещения. На рис. 1 показаны следующие основные конструктивные элементы: 1 – корпус робота; 2 – колесный диск с грунтозацепами; 3 – заднее ведомое колесо; 4 – резервуар со сжатым воздухом; 5 – модули электроники; 6 – аккумуляторные батареи; 7 – управляемый пневматический распределитель; 8 – колёсные приводы постоянного тока; 9 – навесное диагностическое оборудование; 10 – колёсные баллоны с изменяемым объемом для управления нормальной реакцией в точках касания колёс с опорной поверхностью льда.
Рис. 1. Трехмерная модель прототипа подлёдного робота
Fig. 1. Three-dimensional model of the under-ice robot prototype
В начальный момент времени баллоны внутри колёс спущены, робот при этом имеет слабо положительную плавучесть, что позволяет оператору без труда запустить его под лёд через прорубь. Оказавшись под водой, робот задействует распределители, закачивая в колесные баллоны воздух из резервуара высокого давления. Баллоны под действием давления воздуха расширяются, увеличивая свой объём и заполняя собой внутреннюю полость колёсного диска, тем самым увеличивая плавучесть. После обеспечения требуемых нормальных реакций и соответствующих им сил сцепления ведущих колёс робота с внутренней поверхностью льда, задействуется приводная система колёс. Создаваемые крутящие моменты через силы трения сцепления обеспечивают движение робота, что фиксируется набором датчиков, в том числе датчиком угла поворота ведущих колёс и ведомого заднего колеса. Задание параметров движения может осуществляться как заранее оператором, так и осуществляться в режиме реального времени по проводному интерфейсу робота.
В отличие от управляемого движения колёсного робота по поверхности на суше, при движении в вязкой водной среде, где поджатие к поверхности льда обусловлено только силой Архимеда, возникающей за счет увеличения объема колесных баллонов, обеспечение минимизации проскальзывания является сложной задачей, требующей как развитой системы грунтозацепов и протектора, реализующих хорошее сцепление колёс со льдом, так и специального алгоритма системы управления, учитывающего особенности динамики процесса.
Будем рассматривать задачу управления как обеспечение требуемых скоростных параметров движения: Q = (V*, Ω*), где V* – требуемая линейная скорость движения робота; Ω* – скорость поворота робота относительно неподвижной системы координат.
Обобщенная структурная схема системы управления робота и его расчетная схема представлены на рис. 2.
На данной схеме приняты следующие обозначения: RV, RΩ, RωL, RωR – регуляторы, соответственно линейной скорости, скорости разворота корпуса, частоты вращения левого и правого колёс; EncL, EncR, EncС – инкрементальные энкодеры, определяющие фактическую частоту вращения колес. На схеме также показаны основные силы, действующие на робота: FAС, mg – приведенные силы Архимеда и силы тяжести корпуса робота с колёсами со спущенными колёсными баллонами; N1L, N1R, N2 – нормальные реакции в точках контакта, соответственно левого, правого и ведомого колёс; Ffr1L, Ffr1R – силы трения сцепления колес; FAWL и FAWR – силы Архимеда, обусловленные объёмом левого и правого колёсных баллонов; Fμ – сила вязкого сопротивления движению, зависящая от площади лобового сечения и обтекаемости робота; Md1L, Md1R и Mr1L, Mr1R – соответственно моменты электроприводов левого и правого колёс и моменты сил трения качения левого и правого ведущих колёс; Mr2 – момент трения качения ведомого колеса (далее будем считать этот момент пренебрежимо малым и не будем его учитывать при моделировании и составлении уравнений).
Рис. 2. Схема системы управления и расчетная динамическая схема подлёдного робота:
1 – лед; 2 – ведущие колёса; 3 – корпус робота; 4 – ведомое измерительное колесо
Fig. 2. Scheme of the control system and design of the under-ice robot:
1 – ice; 2 – driving wheels; 3 – robot body; 4 – passive measuring wheel
Согласно принятой расчетной схеме, центр масс устройства располагается в точке С и смещен относительно центра (т. А) объёма конструкции. За счет этого и обеспечивается наклонное (колесами вниз) нейтральное положение под водой. Однако после накачки колесных баллонов и создания дополнительных сил Архимеда, робот стремится развернуться, принимая горизонтальное положение и обеспечивая требуемую силу поджатия к поверхности льда, необходимую для возможности движения.
Как видно из схемы системы управления, движение вдоль желаемой траектории может быть осуществлено заданием желаемых законов изменения скорости поступательного движения и поворота корпуса робота. При этом фактическая линейная скорость робота определяется пересчетом угловой скорости ведомого заднего колеса, а угловая скорость разворота может быть получена дифференцированием сигнала с магнитометра, расположенного на корпусе робота.
После расчета требуемых угловых скоростей вращения колес ωL*, ωR*, контуры системы управления приводов рассчитывают напряжения питания для двигателей с учётом показаний датчиков-энкодеров, расположенных непосредственно на ведущих колесах.
Предварительные результаты численного моделирования показали, что в случае превышения предельной силы сцепления между колесом и поверхностью льда происходит проскальзывание, что вызывает нарушение работы системы управления. Очевидно, что для достижения наиболее эффективного энергопотребления и точности движения по траектории необходимо недопущение проскальзывания колёс робота относительно поверхности льда. Ситуация затрудняется тем, что под водой определение фактической линейной скорости робота затруднено. Здесь возможно использование глобальной навигационной спутниковой системы для определения положения устройства и слежения за скоростью, но такая система имеет низкую точность, кроме того, сигнал от спутников может не проходить через толстый слой непрозрачного льда. Использование различных оптических систем слежения за перемещением также может быть затруднено из-за отражающих свойств льда. Наиболее простым и точным способом измерения фактической скорости является оценка угла поворота ведомого колеса, расположенного в задней части робота (см. рис. 2). Использование бесконтактного инкрементального энкодера позволяет обеспечить минимальный момент сопротивления вращению оси колеса.
Далее рассмотрим более подробно эффект проскальзывания ведущих колёс и его влияние на динамику робота. Так как система управления левого и правого колёс идентичны, далее будем рассматривать плоскую задачу, расчетная схема которой с соответствующей структурной схемой системы управления представлена на рис. 3.
Рис. 3. Схема одноканальной системы управления и плоская расчётная схема подлёдного робота:
1 – лёд; 2 – ведущее колесо; 3 – корпус робота; 4 – пассивное колесо
Fig. 3. Planar scheme of the control system and the under-ice robot:
1 – ice; 2 – drive wheel; 3 – robot body; 4 – passive wheel
Обозначения, принятые для данной схемы, соответствуют приведенным ранее (для сил и моментов, действующих на ведущее колесо, индекс не указан). Дополнительно обозначены: Δω = (ω – ωс) – величина проскальзывания колес относительно опорной поверхности, оцениваемая как разность вращений ведущих и ведомого колес; RΔω – регулятор проскальзывания; lB – расстояние между точками контакта колёс робота с поверхностью льда K1 и K2; lС – расстояние от оси вращения ведущих колес до центра масс С; lA – расстояние от оси вращения ведущих колёс до центра (т. А) приложения силы Архимеда; rK – радиус ведущих колёс робота.
В общем случае в качестве регуляторов RV, RΩ, RωL, RωR, RΔω – могут быть использованы как классические ПИД-регуляторы, так и более сложные, а параметры могут быть получены, например, методом численного зондирования параметрического пространства в ходе математического моделирования [16, 17]. В рамках настоящей работы для регулятора величины проскальзывания (RΔω) предлагается использовать П-регулятор, а в качестве регулятора линейной скорости (RV) – инкрементный цифровой ПД-регулятор [18, 19]. При записи выражения для напряжения примем следующие обозначения ei = (Vi*-Vi) – ошибка линейной скорости, на i-том такте работы регулятора. В данном выражении Vi* – заданная оператором линейная скорость движения робота; Vi – текущее значение линейной скорости (на i-том такте работы регулятора), определяемое по скорости вращения ведомого колеса: Vi = ωcirk (ωci – скорость вращения ведомого колеса); Δei = (Vi*-Vi)- (Vi-1*-Vi-1) – изменение ошибки линейной скорости на текущем i-том и предыдущем такте работы регулятора.
Тогда, согласно предлагаемой стратегии управления, напряжение питания двигателя можно представить в виде:
где kp, kd – пропорциональный и дифференциальный коэффициенты регулятора линейной скорости; kpω – коэффициент П-регулятора проскальзывания.
Для получения рациональных параметров регулятора, а также анализа влияния параметров системы на динамику движения, в частности на энергоэффективность движения, в работе выполнено математическое моделирование динамики движения рассматриваемого робота с учётом возможного проскальзывания ведущих колёс.
Согласно схеме, показанной на рис. 3, можем записать выражения для нормальных реакций, возникающих в точках контакта колёс с внутренней поверхностью льда, следующим образом:
Анализ данных формул позволяет понять, что изменение объема баллонов в колёсах приводит к увеличению дополнительной архимедовой силы FAW и реакции N1 ведущих колес, что необходимо для создания сил сцепления и движения объекта по поверхности льда.
Для состояния покоя и движения без проскальзывания, т.е. когда выполняется условие Ffr1 < f1N1, динамика робота определяется системой:
Для движения с проскальзыванием, когда силы сцепления достигают своих предельных значений (Ffr1 = f1N1), динамика определяется системой:
В данных уравнениях: Fμ = μ
Выражение для силы сцепления представим в виде:
При моделировании динамики электропривода воспользуемся уравнением для крутящего момента в следующей форме:
где km – моментный коэффициент двигателя привода колеса; ke – коэффициент противо-ЭДС, ω – частота вращения колеса; R – сопротивление якорной обмотки двигателя; iG, η – передаточное число и коэффициент полезного действия редуктора; U – напряжение питания двигателя, формируемое системой управления в соответствии с выбранной стратегией.
Численные значения основных параметров математической модели робота представлены в табл. 1.
Таблица 1
Основные параметры модели колёсного подлёдного робота
Table 1
Main parameters of the wheeled under-ice robot model
Параметр |
Значение |
Масса корпуса робота m, кг |
7,7 |
Масса колеса робота mk, кг |
0,15 |
Водоизмещение корпуса робота, л |
8 |
Водоизмещение колесных баллонов, л |
0...1,5 |
База опорных элементов по длине корпуса lB, м |
0,5 |
Радиус ведущих колес rK, м |
0,2 |
Расстояние от оси ведущих колёс до центра объема корпуса lA, м |
0,25 |
Расстояние от оси ведущих колёс до центра масс lC, м |
0,2 |
Коэффициент сцепления ведущих колёс f1 |
0,75 |
Коэффициент трения качения ведущих колёс δ, м |
0,005 |
Коэффициент вязкого сопротивления среды μ, Нс/м |
1,2 |
Моментный коэффициент km, Нм/А |
0,0141 |
Коэффициент противо-ЭДС ke, В.с/рад |
0,0139 |
Сопротивление обмотки двигателя R, Ом |
0,196 |
Передаточное число редуктора iG |
4 |
КПД механической передачи η |
0,85 |
Одним из способов обеспечения отсутствия проскальзывания ведущих колёс при разгоне является плавное нарастание требуемой скорости на ведущих колёсах, которое может быть реализовано гладким законом изменения напряжения для приводов. При этом максимальное значение напряжения, а также время разгона, не допускающего проскальзывания, будет определяться предельной силой сцепления, которая, в свою очередь, зависит от физических свойств контактирующих элементов и величины поджатия колёс робота к поверхности льда.
Результаты
Покажем примеры работы системы без контроля проскальзывания с гладким законом изменения заданной линейной скорости V* (рис. 4). Для удобства анализа результатов моделирования, величину проскальзывания представим в виде разности линейной скорости точки на ободе и скорости корпуса робота. При формировании закона изменения задающей скорости используем кусочно-полиномиальную функцию [20].
Рис. 4. Результаты моделирования динамики колесного робота с программно-заданным законом
изменения линейной скорости V*:
1 – угловая скорость вращения ведущего колеса (рад/с); 2 – линейная скорость корпуса робота (м/с); 3 – сила сцепления ведущих колес (Н); 4 – крутящий момент на ведущих колесах (Нм); 5 – величина проскальзывания (м/с); 6 – механическая мощность привода робота
Fig. 4. Results of modeling the dynamics of a wheeled robot with a software-specified law
of change in the linear velocity V*:
1 – velocity of the drive wheel (rad/s); 2 – linear speed of the robot body (m/s); 3 – friction force (N); 4 – actuator torque (Nm); 5 – slippage value (m/s); 6 – actuator mechanical power
Как следует из анализа результатов моделирования, при заданной максимальной желаемой скорости движения (5 м/с) и указанных параметров привода, даже плавный разгон вызывает появление проскальзывание на 7 секунде. Далее робот продолжает разгон, но уже с пробуксовывающими колёсами. При этом максимальный крутящий момент достигает значения 1,9 Нм, а максимальное проскальзывание составляет 1,4 м/с. Так как рост частоты вращения двигателя не ограничен системой управления, то через 45 с приводы выходят на установившуюся частоту вращения около 55 рад/с, что соответствует механической мощности 57 Вт.
Для наглядной демонстрации эффективности работы антипробуксовочной системы, покажем результаты моделирования при неоптимальном коэффициенте регулятора проскальзывания
Рис. 5. Результаты моделирования динамики колесного робота с П-регулятором проскальзывания (
Fig. 5. Results of simulating the dynamics of a wheeled robot with a slip-controller (
Как очевидно из результатов моделирования, введение ограничения подводимой мощности в соответствии с формулой (1), приводит к значительному снижению угловой скорости колеса, хотя и не устраняет проскальзывание полностью. Однако использование даже неоптимального параметра регулятора, позволило сократить максимальное значение проскальзывания до 0,35 м/с, а требуемую мощность до 34 Вт. При этом темп нарастания скорости корпуса практически не изменился, так как динамика разгона определяется исключительно свойствами среды и сцеплением робота с поверхностью при обеспечении достаточного крутящего момента со стороны приводной системы.
На рис. 6. показаны результаты численного моделирования системы с оптимизированным значением коэффициента пропорциональности в контуре регулирования проскальзывания. Дальнейшее увеличение коэффициента приводит к появлению колебательности и потери устойчивости системы.
Рис. 6. Результаты моделирования динамики колесного робота с оптимизированным П-регулятором проскальзывания (
Fig. 6. Simulation result s of the dynamics of a wheeled robot with an optimized slip-controller (
Результаты данного эксперимента показали возможность практически полного исключения проскальзывания (максимальное значение 0,08 м/с) с сохранением динамики разгона робота. Максимальная потребляемая мощность при этом не превышает 27 Вт, что практически в 2 раза ниже, чем в первом эксперименте, где система контроля проскальзывания не была задействована.
Для оценки энергоэффективности предложенной антипробуксовочной системы управления приводами ведущих колёс мобильного робота на рис. 7 приведены временные диаграммы:
а) б)
Рис. 7. Эффективность преобразования энергий при работе с регулятором (1)
и без регулятора проскальзывания (2)
Fig. 7. Energy conversion efficiency with regulator (1) and without slip-regulator (2)
Согласно результатам моделирования, использование системы контроля проскальзывания в 2 раза увеличивает к.п.д механического преобразования движителя, значение которого достигает уровня около 94 %. Общая энергоэффективность, выражаемая отношением механической энергии движения к затрачиваемой электрической, при использовании регулятора повышается с 19 % до 26 %.
Заключение
Роботизированный мониторинг объектов подводной инфраструктуры водоёмов, покрытых льдом, может быть осуществлён посредством колёсного робота, перемещающегося по внутренней поверхности льда. Перемещение на границе раздела сред дает ряд преимуществ, среди которых высокое быстродействие и эффективность, точность перемещений, возможность сохранять своё положение под действием внешних сил, не используя приводную систему.
Однако указанные преимущества требуют контроля проскальзывания ведущих колёс относительно поверхности льда. Предложенная в работе структура системы управления, позволяющая ограничивать подводимую мощность в соответствии с величиной проскальзывания, позволяет минимизировать энергопотери за счет пробуксовки. В работе предложена математическая модель, представленная системой нелинейных дифференциальных уравнений, позволяющих моделировать процесс управляемого движения как в режиме движения без проскальзывания, так и с его учётом. Результаты численных экспериментов на математической модели устройства доказали эффективность предлагаемой методики. Для используемых в модели параметров с помощью контроля проскальзывания удалось сократить энергопотребление в 2 раза с сохранением динамики разгона робота.
1. Харин И.О., Князев Н.П., Скуднев Д.М. Межконтинентальные сети коммуникации // Информационные технологии в процессе подготовки современного специалиста. - 2017. - С. 113-118.
2. Бабкин Ю.А., Федоров А.В. Автоматизация процессов измерений и контроля линейного тракта подводных волоконно-оптических линий связи для Арктической зоны Российской Федерации // Арктика: история и современность. - 2017. - С. 200-209.
3. Юха С. Арктические подводные коммуникационные кабели и региональное развитие северных территорий // Арктика и Север. - 2018. - №. 32. - С. 63-81.
4. Миддлтон А. Сила коммуникационного взаимодействия в Арктике: участие граждан в арктических институтах // Арктика и Север. - 2021. - №. 42. - С. 170-185.
5. Piazza P., Cummings V. J., Lohrer D. M. et al. Divers-operated underwater photogrammetry: Applications in the study of Antarctic benthos // International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences. - 2018. - Т. 42. - С. 885-892.
6. Матвиенко Ю.В., Ремезков А.В. Концепция создания роботизированного комплекса обследования и мониторинга технического состояния объектов подводной добычи // Газовая промышленность. - 2020. - №. S2. - С. 172-179.
7. Шматков А.А., Шматков А.Г., Сергейчиков О.А. Современное состояние производства и перспективы развития серийных российских телеуправляемых необитаемых подводных аппаратов // Газовая промышленность. - 2020. - №. S2. - С. 156-160.
8. Яцун С.Ф., Лушников Б.В., Казарян К.Г. Конструктивные особенности бионического робота-рыбы // Известия Юго-Западного государственного университета. Серия: Техника и технологии. - 2017. - Т. 7. - №. 2. - С. 94-102.
9. Jatsun S., Lushnikov B., Politov E. et al. Underwater floating robot-fish: a comparative analysis of the results of mathematical modelling and full-scale tests of the prototype // MATEC Web of Conferences - 2017. С. 02014.
10. Wu Y., Ta X., Xiao R. et al. Survey of underwater robot positioning navigation // Applied Ocean Research. - 2019. - Т. 90. - С. 101845.2
11. Cui R., Chen L., Yang C.et al. Extended state observer-based integral sliding mode control for an underwater robot with unknown disturbances and uncertain nonlinearities // IEEE Transactions on Industrial Electronics. - 2017. - Т. 64. - №. 8. - С. 6785-6795.
12. Petillot Y.R., Antonelli G., Casalino G. et al. Underwater robots: From remotely operated vehicles to intervention-autonomous underwater vehicles // IEEE Robotics & Automation Magazine. - 2019. - Т. 26. - №. 2. - С. 94-101.
13. Цыфанский С.Л. Новый способ движения на границе природных сред подледный транспортный аппарат // Наука, техника и образование. - 2015. -№4 (10). - С. 43-49.
14. Jakuba M.V., German C.R., Bowen A.D. et al. Teleoperation and robotics under ice: Implications for planetary exploration // 2018 IEEE Aerospace Conference. - IEEE, 2018. - С. 1-14.
15. Berisford D.F., Leichty J.M., Klesh A.T. et al. Remote under-ice roving in Alaska with the buoyant rover for under-ice exploration // AGUF all Meeting Abstracts. - 2013. - Т. 2013. - С. C13C-0684.
16. Jatsun S., Malchikov A., Yatsun A. Adaptive Control System for DC Electric Drive under Uncertainty //2020 International Conference on Industrial Engineering, Applications and Manufacturing (ICIEAM). - IEEE, 2020. - С. 1-5.
17. Jatsun S., Malchikov A., Yatsun A. Comparative analysis of the industrial exoskeleton control systems // Proceedings of 14th International Conference on Electromechanics and Robotics “Zavalishin's Readings” ER (ZR). - Springer Singapore, 2020. - С. 63-74.
18. Денисенко В.В. ПИД-регуляторы: принципы построения и модификации // Современные технологии автоматизации. - 2006. - № 4. - С. 66-74; 2007. - № 1. - С. 78-88.
19. Borase R.P., Maghade D.K., Sondkar S.Y. et al. A review of PID control, tuning methods and applications // International Journal of Dynamics and Control. - 2021. - Т. 9. - С. 818-827.
20. Jatsun S., Malchikov A., Yatsun A. et al. Simulation of a walking robot-exoskeleton movement on a movable base // Journal of Artificial Intelligence and Technology. - 2021. - Т. 1. - №. 4. - С. 207-213.