Введение. В последние годы автономные мобильные роботы значительно усовершенствовались, и особое внимание уделяется автоматизированным управляемым транспортным средствам (АУТС) [1–5]. Одной из главных проблем, с которыми сталкиваются такие системы, является улучшение маневренности. Одним из эффективных решений этой проблемы является использование колес Илона – конструктивного элемента, представляющего собой ступицу с роликами, расположенными под углом на окружности [6, 7].Колеса Илона значительно расширяют возможности движения мобильных роботов, позволяя выполнять сложные траектории и адаптироваться к различным условиям окружающей среды [8–12]. Кинематические и динамические модели роботов с такими колесами подробно исследовались в работах, где анализировалось влияние параметров робота на его поведение [6, 8, 9]. Применение этих колес позволяет реализовать высокоточные алгоритмы управления движением, включая избегание препятствий в динамической среде [10] и эффективное передвижение по сложным поверхностям, таким как лестницы и пересеченная местность [12].Несмотря на обширное количество публикаций, посвященных мобильным роботам с тремя или четырьмя колесами Илона, работы, исследующие модели с шестью колесами, остаются крайне редкими [13–16]. Этот аспект открывает перспективы для разработки новых моделей и анализа их поведения в различных условиях.В большинстве исследований применяются платформы прямоугольной формы, которые обеспечивают равномерное распределение нагрузки и обладают преимуществами в плане грузоподъёмности. Однако в некоторых сценариях, таких как сложные маневры или управление на высоких скоростях, круглая платформа может демонстрировать лучшие характеристики. Одним из ключевых преимуществ круглой платформы является её способность обеспечивать более стабильное распределение центробежных сил при выполнении резких поворотов. В прямоугольных конструкциях возможны неравномерные нагрузки на отдельные колёса во время маневров, что может привести к увеличенному износу или ухудшению динамических характеристик. В случае круглой платформы распределение нагрузки остаётся более симметричным, что способствует улучшенной устойчивости при вращении. Кроме того, круглая форма платформы может упростить алгоритмы навигации и управления, особенно в условиях, когда требуется равномерное поведение робота при смене направления. Такая конструкция снижает влияние ориентации корпуса на траекторию движения, что может быть полезным в задачах, требующих высокой точности позиционирования.Целью данной работы является создание динамической модели мобильного робота с шестью колесами Илона. Кроме того, работа направлена на проверку модели через моделирование движения по прямолинейной и криволинейной траекториям.Остальная часть исследования включает вывод уравнений динамической модели, моделирование управления движением, а также обсуждение полученных результатов.Методы, оборудование, материалы. Исследование выполнено с использованием методов теоретической механики и численного моделирования. Для разработки динамической модели мобильного робота применялись уравнения Лагранжа второго рода, что позволило учесть основные параметры движения и взаимодействия робота с окружающей средой.Моделирование проводилось в программной среде MATLAB с использованием инструмента Simulink и специализированного набора библиотек Robotics Toolbox. Эти инструменты обеспечили возможность анализа движения робота по различным траекториям, включая прямолинейные и криволинейные маршруты.В процессе разработки модели движение робота рассматривалось на основе предположения о чистом качении колёс без проскальзывания, что упрощает выведение уравнений системы. Особое внимание было уделено точности формулировки кинетической энергии системы и выведению уравнений движения, описывающих поведение робота в горизонтальной плоскости.Описание динамики робота. Графическое изображение рассматриваемого робота представлено на рис. 1. Рис. 1. Аналитическая модель мобильного роботас колесами Илона [15] В данной работе использованы уравнения Лагранжа второго рода. Для их формулировки предполагается следующий вектор обобщенных координат            q=&φ1&φ2&φ3&φ4&φ5&φ6                               (1)В свою очередь, кинематические уравнения робота выводятся с учетом отсутствия скольжения между землей и колесами Илона. Для получения динамических уравнений движения характеристической точки робота был использован уравнение Лагранжа. Уравнение Лагранжа второго рода можно записать в следующей матричной форме:ddt∂L∂q-∂L∂q=Q                   (2)где L=L(q,q) – кинетический потенциал системы; q=q1,q2,q3,q4,q5,q6T – вектор обобщенных координат; Q=Q1,Q2,Q3,Q4,Q5,Q6T – вектор обобщенных сил.Кинетический потенциал L системы удовлетворяет следующему соотношению:            L=E-U                               (3)где E – полная кинетическая энергия системы, U – потенциальная энергия системы.Предполагается, что робот движется по плоской поверхности, поэтому уравнение (3) упрощается до следующего вида:            L=E                                   (4)Затем была определена полная кинетическая энергия системы в предположении, что ролики являются пассивными элементами в анализируемой системе. Следовательно, полная кинетическая энергия системы E составляет:            E=Ep+Ew                           (5)где Ep – кинетическая энергия платформы, Ew – кинетическая энергия колес Илона.Платформа робота движется в плоскости, поэтому кинетическая энергия платформы Ep может быть записана в следующем виде:            Ep=12mpvx2+vy2+12Ipβ2             (6)где mp – масса платформы, Ip – массовый момент инерции платформы робота, определенный относительно оси zL, проходящей через точку o.колеса Илона перемещаются в сложном движении, состоящем из плоского подъемного движения платформы робота и вращательного движения колес, определяемого вокруг оси собственного вращения колеса Илона. Следовательно, кинетическая энергия колес Ew удовлетворяет следующему соотношению: Ew=612mwvx2+vy2+12IwzLβ2+12Iwφ12+φ22+φ32+φ42+φ52+φ62                       (7) где mw – масса колеса Илона, Iw – массовый момент инерции колеса Илона, определяемый относительно оси, собственное вращение колеса, IwzL - массовый момент инерции колеса Илона, определяемый относительно оси zL, проходящей через точку o. Кинетический потенциал системы можно записать в следующем виде:             L=12mp+3mwvx2+vy2+12Ip+3IwzLβ2+12Iwφ12+φ22+φ32+φ42+φ52+φ62  (8)Будет предполагаться, что:            mpс=mp+6mw                                                             (9)            Ipс=Ip+6IwzL                                                            (10)Подставляя уравнения (9) и (10) в уравнение (8), получаем:            L=12mpcvx2+vy2+12Ipcβ2+12Iwφ12+φ22+φ32+φ42+φ52+φ62              (11)Подставляя уравнения кинематической модели [15] в уравнение (11), получаем:            L=12mpcrφ14+φ28+φ38+φ44+φ58+φ682++r-φ14+φ28+φ38-φ44+φ58+φ682++12Ipcr4R3+3-3+1φ1+3+1φ2+-3+1φ3+3+1φ4-2φ5+2φ62++12Iwφ12+φ22+φ32+φ42+φ52+φ62                 (12)Будем считать, что:            A=mpcr2128                                                                   (13)            B=Ipcr232R23+63+9=Ipcr232R212+63=Ipcr2192R22+3                                    (14)            C=Iw2                                                                      (15)Подставляя уравнения (13), (14) и (15) в уравнение (12), получаем:            L=A2φ1+φ2+φ3+2φ4+φ5+φ62+-2φ1+φ2+φ3-2φ4+φ5+φ62++B-3+1φ1+3+1φ2-3+1φ3+3+1φ4-2φ5+2φ62++Cφ12+φ22+φ32+φ42+φ52+φ62 (16)⇒L=8A+4+23B+Cφ12+2A+4+23B+Cφ22++2A+4+23B+Cφ32+8A+4+23B+Cφ42+2A+4B+Cφ52++2A+4B+Cφ62-24+23Bφ1φ2+24+23Bφ1φ3++16A-24+23Bφ1φ4+43+1Bφ1φ5-43+1Bφ1φ6++4A-24+23Bφ2φ3+24+23Bφ2φ4+4A-43+1Bφ2φ5++4A+43+1Bφ2φ6-24+23Bφ3φ4+4A+43+1Bφ3φ5++4A-43+1Bφ3φ6-43+1Bφ4φ5+43+1Bφ4φ6+4A-8Bφ5φ6  (17) Чтобы установить конкретные обобщенные силы, действующие на отдельные колеса, силы и моменты сил, действующие на каждое колесо, графически показаны на рисунке 2а и 2б.    (a) (б)Рис. 2. (а) Силы и моменты сил, действующих на колеса № 1, 3, 5,(б) Силы и моменты сил, действующих на колеса № 2, 4, 6  i-е колесо Илона приводится в движение за счет тягового момента τi (i = 1,2,3,4,5,6), исходящего от i-й системы передачи энергии. Кроме того, предполагалось, что вес платформы колесного мобильного робота равен Gp, тогда как вес i-го колеса Илона равен Gi, и оно движется по плоской, шероховатой (коэффициент сухого трения μi) и деформируемой (коэффициент трения качения fi) поверхность без заноса. При анализе также учитывается сила сухого трения Ti и сила давления Ni для i-го колеса.Работа системы сил составляет:             D=Qiqi,(i=1,2,..,6)=τ1-T1rφ1+τ2-T2rφ2+τ3-T3rφ3++τ4-T4rφ4+τ5-T5rφ5+τ6-T6rφ6                  (18) Величину сухого трения Ti можно вывести из рисунка 3, который показывает соотношение сил, действующих на колесо:            T=RFsinθ                       (19)где RF – сила реакции. Поскольку θ - малый угол, можно найти следующее:             T=RFsinθ≈RFtanθ≈RFθ  (20)Теперь, умножив и разделив r1 на r1:             T=RFθr1r1   но    f=r1θ⇒T=RFfr1⇒Tr1=RFf                                   (21)Возвращаясь к рисунку 4, можно найти:RF2=G2+T2=(-N)2+T2=N2+T2⇒RF=N2+T2                         (22)Подставляя уравнение (22) в уравнение (21): Tr1=N2+T2f⇒T2r12=N2+T2f2⇒T2r12-f2=N2f2⇒T=Nfr12-f2=Nfr   (23) Рис. 3. Силы, действующие на колесоУравнение (23) компенсируется уравнением (18):D=τ1-N1f1rrφ1+τ2-N2f2rrφ2+τ3-N3f3rrφ3++τ4-N4f4rrφ4+τ5-N5f5rrφ5+τ6-N6f6rrφ6             ⇒D=τ1-N1f1φ1+τ2-N2f2φ2+τ3-N3f3φ3++τ4-N4f4φ4+τ5-N5f5φ5+τ6-N6f6φ6                             (24)Таким образом, мощность обобщенных сил составляет:            P=Qiqi,(i=1,2,..,6)=τ1-N1f1φ1+τ2-N2f2φ2+τ3-N3f3φ3++τ4-N4f4φ4+τ5-N5f5φ5+τ6-N6f6φ6            (25)Но при расчете величины сухого трения необходимо учитывать направление угловой скорости (знак угловой скорости), поэтому уравнение (25) принимает следующий вид:            P=Qiqi,(i=1,2,..,6)=τ1-N1f1sgn(φ1)φ1+τ2-N2f2sgn(φ2)φ2++τ3-N3f3sgn(φ3)φ3+τ4-N4f4sgn(φ4)φ4++τ5-N5f5sgn(φ5)φ5+τ6-N6f6sgn(φ6)φ6         (26)Уравнение (26) показывает, что вектор обобщенных сил можно записать в следующем виде:Q=Q1,Q2,Q3,Q4,Q5,Q6T=&τ1-N1f1sgn(φ1),τ2-N2f2sgn(φ2),τ3-N3f3sgn(φ3),&τ4-N4f4sgn(φ4),τ5-N5f5sgn(φ5),τ6-N6f6sgn(φ6)T  (27)Наконец, используя уравнения Лагранжа второго рода, получаются следующие системные уравнения:28A+4+23B+Cφ1-24+23Bφ2+24+23Bφ3++16A-24+23Bφ4+43+1Bφ5-43+1Bφ6=τ1-N1f1sgn(φ1)  (28)22A+4+23B+Cφ2-24+23Bφ1++4A-24+23Bφ3+24+23Bφ4++4A-43+1Bφ5+4A+43+1Bφ6=τ2-N2f2sgn(φ2)                (29)            22A+4+23B+Cφ3+24+23Bφ1++4A-24+23Bφ2-24+23Bφ4++4A+43+1Bφ5+4A-43+1Bφ6=τ3-N3f3sgn(φ3)                (30)            28A+4+23B+Cφ4+16A-24+23Bφ1++24+23Bφ2-24+23Bφ3+-43+1Bφ5+43+1Bφ6=τ4-N4f4sgn(φ4)        (31)      22A+4B+Cφ5+43+1Bφ1+4A-43+1Bφ2++4A+43+1Bφ3-43+1Bφ4+4A-8Bφ6=τ5-N5f5sgn(φ5)     (32)22A+4B+Cφ6-43+1Bφ1+4A+43+1Bφ2++4A-43+1Bφ3+43+1Bφ4+4A-8Bφ5=τ6-N6f6sgn(φ6)      (33) Итак, мы можем написать:            &Q1&Q2&Q3&Q4&Q5&Q6=M&φ1&φ2&φ3&φ4&φ5&φ6                      (34)Учитывая, что vx,vy — центральные скорости робота в локальных координатах, их можно выразить в обобщенных координатах следующим образом:            vx=xgcosβ+ygsinβ            (35)            vy=-xgsinβ+ygcosβ           (36)Здесь уравнения (35) и (36) должны быть выведены по времени, а затем приравнены к прямой кинематической модели после ее вывода по времени. Таким образом, мы получаем:             xgcosβ-xgβsinβ+ygsinβ+ygβcosβ=rφ14+φ28+φ38+φ44+φ58+φ68                  (37)            -xgsinβ-xgβcosβ+ygcosβ-ygβsinβ=r-φ14+φ28+φ38-φ44+φ58+φ68                (38)            β=r4R3+3&-3+1φ1+3+1φ2-3+1φ3+&+3+1φ4-2φ5+2φ6                                        (39) Здесь мы должны найти уравнения для φ1,φ2,φ3,φ4,φ5,φ6  из уравнения (34), затем подставить полученные уравнения в уравнения (37), (38), (39), чтобы найти требуемую динамическую модель.Из уравнения (34) получаем: φ1=63+20BC2+323+128ABC+C3+8AC2283+24BC3+1283+384ABC2+C4+16AC3(τ1-N1f1sgn(φ1))++2+3BC83+24BC2+C3(τ2-N2f2sgn(φ2))-2+3BC83+24BC2+C3(τ3-N3f3sgn(φ3))++2+3BC2-163+64ABC-4AC283+24BC3+1283+384ABC2+C4+16AC3(τ4-N4f4sgn(φ4))++-(1+3)B(83+24)BC+C2(τ5-N5f5sgn(φ5))+(1+3)B(83+24)BC+C2(τ6-N6f6sgn(φ6))         (40) φ2=2+3BC83+24BC2+C3(τ1-N1f1sgn(φ1))++63+20BC2+323+112ABC+C3+6AC2283+24BC3+643+192ABC2+C4+8AC3(τ2-N2f2sgn(φ2))++2+3BC-8AB-AC83+24BC2+643+192ABC+C3+8AC2(τ3-N3f3sgn(φ3))++-2+3BC83+24BC2+C3(τ4-N4f4sgn(φ4))++1+3BC2-16ABC-AC283+24BC3+643+192ABC2+C4+8AC3(τ5-N5f5sgn(φ5))++-(1+3)BC2-(163+32)ABC-AC2(83+24)BC3+(643+192)ABC2+C4+8AC3(τ6-N6f6sgn(φ6))                    (41)φ3=&-(2+3)BC(83+24)BC2+C3(τ1-N1f1sgn(φ1))+&+(2+3)BC-8AB-AC(83+24)BC2+(643+192)ABC+C3+8AC2(τ2-N2f2sgn(φ2))+&+(63+20)BC2+(323+112)ABC+C3+6AC22((83+24)BC3+(643+192)ABC2+C4+8AC3)(τ3-N3f3sgn(φ3))+&+(2+3)BC(83+24)BC2+C3(τ4-N4f4sgn(φ4))+&+-(1+3)BC2-(163+32)ABC-AC2(83+24)BC3+(643+192)ABC2+C4+8AC3(τ5-N5f5sgn(φ5))+&+(1+3)BC2-16ABC-AC2(83+24)BC3+(643+192)ABC2+C4+8AC3(τ6-N6f6sgn(φ6))                  (42)            φ4=2+3BC2-163+64ABC-4AC283+24BC3+1283+384ABC2+C4+16AC3(τ1-N1f1sgn(φ1))++-(2+3)BC(83+24)BC2+C3(τ2-N2f2sgn(φ2))+(2+3)BC(83+24)BC2+C3(τ3-N3f3sgn(φ3))++63+20BC2+323+128ABC+C3+8AC2283+24BC3+1283+384ABC2+C4+16AC3(τ4-N4f4sgn(φ4))++(1+3)B(83+24)BC+C2(τ5-N5f5sgn(φ5))+-(1+3)B(83+24)BC+C2(τ6-N6f6sgn(φ6))       (43)            φ5=&-(1+3)B(83+24)BC+C2(τ1-N1f1sgn(φ1))+&+(1+3)BC2-16ABC-AC2(83+24)BC3+(643+192)ABC2+C4+8AC3(τ2-N2f2sgn(φ2))+&+-(1+3)BC2-(163+32)ABC-AC2(83+24)BC3+(643+192)ABC2+C4+8AC3(τ3-N3f3sgn(φ3))+&+(1+3)B(83+24)BC+C2(τ4-N4f4sgn(φ4))+&+(83+20)BC2+(483+112)ABC+C3+6AC22((83+24)BC3+(643+192)ABC2+C4+8AC3)(τ5-N5f5sgn(φ5))+&+2BC2-(83+8)ABC-AC2(83+24)BC3+(643+192)ABC2+C4+8AC3(τ6-N6f6sgn(φ6))                 (44)            φ6=&(1+3)B(83+24)BC+C2(τ1-N1f1sgn(φ1))+&+-(1+3)BC2-(163+32)ABC-AC2(83+24)BC3+(643+192)ABC2+C4+8AC3(τ2-N2f2sgn(φ2))+&+(1+3)BC2-16ABC-AC2(83+24)BC3+(643+192)ABC2+C4+8AC3(τ3-N3f3sgn(φ3))+&+-(1+3)B(83+24)BC+C2(τ4-N4f4sgn(φ4))+&+2BC2-(83+8)ABC-AC2(83+24)BC3+(643+192)ABC2+C4+8AC3(τ5-N5f5sgn(φ5))+&+(83+20)BC2+(483+112)ABC+C3+6AC22((83+24)BC3+(643+192)ABC2+C4+8AC3)(τ6-N6f6sgn(φ6))                  (45)Теперь, подставляя уравнения (40, …, 45) в уравнения (37), (38), (39), получаем:            xg=-ygβ+r&cosβ-sinβ&τ2+τ3+τ5+τ6-N2f2sgnφ2-N3f3sgnφ3-&-N5f5sgnφ5-N6f6sgnφ6*&*83+24B+C1683+24BC+643+192AB+C2+8AC+&+cosβ+sinβτ1+τ4-N1f1sgnφ1-N4f4sgnφ4*&*83+24B+C883+24BC+1283+384AB+C2+16AC  (46)             yg=xgβ+r&sinβ+cosβ&τ2+τ3+τ5+τ6-N2f2sgnφ2-N3f3sgnφ3-&-N5f5sgnφ5-N6f6sgnφ6*&*83+24B+C1683+24BC+643+192AB+C2+8AC+&+sinβ-cosβτ1+τ4-N1f1sgnφ1-N4f4sgnφ4*&*83+24B+C883+24BC+1283+384AB+C2+16AC  (47)             β=r4R3+3&&-3+143+16BC+643+256AB+C2+16AC283+24BC2+1283+384ABC+C3+16AC2+&+8+233+1B83+24BC+C2*&*τ1-τ4-N1f1sgnφ1+N4f4sgnφ4+&+&-22+33+1B83+24BC+C2+&+43+83+1BC+963+160AB+3+1C2+83+1AC283+24BC2+643+192ABC+C3+8AC2*&*τ2-τ3-N2f2sgnφ2+N3f3sgnφ3+&+&23+12B83+24BC+C2+&+-8-43BC+-323-64AB-C2-8AC83+24BC2+643+192ABC+C3+8AC2*&*τ5-τ6-N5f5sgnφ5+N6f6sgnφ6                           (48)Из уравнений (46), (47), (48) можно записать следующую общую форму:            q=Uq,qq+Hqτ+WqSNF                                          (49)где:q=&xg&yg&β,q=&xg&yg&β,Uq,q=0-β0β00000,τ=&τ1&τ2&τ3&τ4&τ5&τ6,f=&f1&f2&f3&f4&f5&f6, N=N1000000N2000000N3000000N4000000N5000000N6,S=sgnφ1000000sgnφ2000000sgnφ3000000sgnφ4000000sgnφ5000000sgnφ6,HTq=cosβ+sinβh1sinβ-cosβh1h3cosβ-sinβh2cosβ+sinβh2h4cosβ-sinβh2cosβ+sinβh2-h4cosβ+sinβh1sinβ-cosβh1-h3cosβ-sinβh2cosβ+sinβh2h5cosβ-sinβh2cosβ+sinβh2-h5,Wq=-Hq гдеh1=r83+24B+C883+24BC+1283+384AB+C2+16AC,h2=r83+24B+C1683+24BC+643+192AB+C2+8AC, h3=r4R3+3-3+143+16BC+643+256AB+C2+16AC283+24BC2+1283+384ABC+C3+16AC2+8+233+1B83+24BC+C2, h4=r4R3+3-22+33+1B83+24BC+C2+43+83+1BC+963+160AB+3+1C2+83+1AC283+24BC2+643+192ABC+C3+8AC2, h5=r4R3+323+12B83+24BC+C2+-8-43BC+-323-64AB-C2-8AC83+24BC2+643+192ABC+C3+8AC2 . Из уравнения (49) можно написать: τ=H-1qq-H-1qUq,qq+SNF     (50)Полученная математическая модель была использована для проведения численных моделирований движения выбранной характерной точки робота.Численное моделирование. В ходе численных исследований была проведена симуляция для отображения движения характерной точки робота по двум траекториям: прямолинейной и криволинейной. Моделирование уравнения (50) выполнялось в среде Matlab Simulink. Для проведения численного моделирования движения выбранной точки робота были приняты следующие геометрические параметры: δ = π/4 [рад], R = 0.5 [м], r = 0.052 [м]. Другие константы, использованные в симуляции: f1 = f2 = f3 = f4 = f5 = f6 = 0.002 [м], N1 = N2 = N3 = N4 = N5 = N6 = 1455.1 [Н], Iw = 0.4688 [кг·м²], mp = 20 [кг], mw = 0.5 [кг].Были заданы следующие начальные условия: x00=0,  y0(0)=0 , vx0=0, vy0=0,  β0=0,        β(0)=0  в симуляции. В свою очередь, каркас робота движется с постоянной конфигурацией в течение симуляции: β(t)=β(0) .1. Движение по прямолинейной траектории:На рисунке 4 показана смоделированная прямолинейная траектория.  Рис. 4. Смоделированная прямолинейная траектория На рисунке 5 показано, как ведет себя набор скоростей для выбранной точки робота: vx(t),vy(t),β(t) .На рисунке 6 показано, как ведет себя набор ускорений выбранной точки робота: vx(t),vy(t),β(t) .Рисунки 7 и 8 иллюстрируют решение задачи обратной кинематики: поведение угловых скоростей и угловых ускорений отдельных колес.На рисунке 9 показано решение задачи обратной динамики: поведение приводных моментов отдельных колес.Рисунок 10 демонстрирует величину ошибки между эталонной скоростью и скоростью, полученной в модели прямой динамики.  Рис. 5. Набор скоростей для выбранной точки роботаРис. 6. Набор ускорений для выбранной точки роботаРис. 7. Поведение угловых скоростей колесРис. 8. Поведение угловых ускорений колесРис. 9. Поведение приводных моментов колес Рис. 10. Ошибка скорости   2. Движение по криволинейной траекторииРисунок 11 показывает смоделированную криволинейную траекторию.На рисунке 13 показано, как ведет себя набор ускорений выбранной точки робота: vx(t),vy(t),β(t) . На рисунке 12 показано, как ведет себя набор скоростей для выбранной точки робота: vx(t),vy(t),β(t) .Рисунки 14 и 15 иллюстрируют решение задачи обратной кинематики: поведение угловых скоростей и угловых ускорений отдельных колес.На рисунке 16 показано решение задачи обратной динамики: поведение приводных моментов отдельных колес.Рисунок 17 демонстрирует величину ошибки между эталонной скоростью и скоростью, полученной в модели прямой динамики.Из рисунков 10 и 17 видно, что ошибки очень малы. Они проявляются только в переходных процессах, когда скорость изменяется. Почти их нет, когда скорость постоянна, что подтверждает правильность полученных результатов. Рис. 11. Сымитированная криволинейная траектория Рис. 12. Набор скоростей для выбранной точкироботаРис. 13. Набор ускорений для выбранной точкиробота Рис. 14. Поведение угловых скоростей колес Рис. 15. Поведение угловых ускорений колесРис. 16. Поведение приводных моментов колес   Рис. 17. Ошибка скоростиВывод.В рамках данного исследования была разработана динамическая модель мобильного робота с шестью колесами Илона. Модель была построена на основе уравнений Лагранжа второго рода с использованием множителей. Для проверки корректности модели было проведено численное моделирование в среде Matlab Simulink, включающее движение робота по двум траекториям: прямолинейной и криволинейной.Результаты моделирования показали:На прямолинейной траектории наблюдались разные знаки приводных моментов колес после стабилизации скорости: моменты τ2, τ3, τ5, τ6 были положительными, а τ1, τ4 – отрицательными.На криволинейной траектории все знаки моментов были положительными после стабилизации скорости.Данный результат демонстрирует, что прямолинейное движение способствует снижению отрицательных эффектов трения на колеса. Дополнительно, малые значения ошибки скорости подтвердили точность и корректность предложенной модели.Разработанная модель может быть использована в будущем для проектирования и управления мобильными роботами в условиях, требующих высокой маневренности.В перспективе дальнейшие исследования могут быть направлены на разработку обобщённой модели мобильных платформ с шестью колесами Илона, включающей различные геометрические конфигурации, такие как круг, овал и прямоугольник. Такой подход позволит избежать необходимости проведения отдельных исследований для каждой возможной конструкции и обеспечит возможность их сопоставления с точки зрения динамических характеристик. Кроме того, перспективным направлением является оптимизация геометрических параметров платформы и численный анализ её поведения в сравнении с традиционными прямоугольными платформами, представленными в работах [13]-[14]. Это позволит глубже понять преимущества различных конфигураций и повысить эффективность проектирования мобильных роботов.Таблица 1 Список обозначенийАУТСАвтоматически управляемые транспортные средстваO,xl,yl,zlсистема координатδугол между осью колеса и осью качения [рад]αугол между осью xl и осью колеса a [рад] φ угловая скорость колеса Илона [рад/с] rрасстояние между осью колеса и землей [м] vx скорость центра робота (O) в направлении x [м/с] vy скорость центра робота (O) в направлении y [м/с] β угловая скорость робота [рад/с] ax,ay координаты центра колеса [м]Oцентр массы роботаRрадиус платформы робота [м]qвектор обобщенных координатQвектор обобщенных силLкинетический потенциал системы [кг·м²·с⁻²]Eполная кинетическая энергия системы [кг·м²·с⁻²]Uпотенциальная энергия системы [кг·м²·с⁻²]Epкинетическая энергия платформы [кг·м²·с⁻²]Ewкинетическая энергия колес Илона [кг·м²·с⁻²]mpмасса платформы [кг]Ipмассовый момент инерции платформы робота [кг·м²]mwмасса колеса Илона [кг]Iwмассовый момент инерции колес Илона относительно оси вращения колеса [кг·м²]IwzLмассовый момент инерции колеса Илона относительно оси ez, проходящей через точку o [кг·м²]mpсмасса робота [кг]Ipсмассовый момент инерции робота [кг·м²]τкрутящий момент двигателя колеса [кг·м²·с⁻²]Gpвес платформы мобильного робота [кг·м·с⁻²]Giвес i-го колеса Илона [кг·м·с⁻²]μкоэффициент сухого тренияfкоэффициент трения качения [м]Tсила сухого трения [кг·м·с⁻²]Nсила давления [кг·м·с⁻²]Dработа системы сил [кг·м²·с⁻²]RFреакционная сила [кг·м·с⁻²] Pмощность обобщенных сил [кг·м²·с⁻³]