CC BY
8
DOI: 10.53078/20778481_2021_3_108 УДК 629.114.2 Э. И. Ясюкович
ИМИТАЦИОННОЕ МОДЕЛИРОВАНИЕ УПРАВЛЯЕМОГО КУРСОВОГО ДВИЖЕНИЯ МНОГООСНЫХ КОЛЕСНЫХ МАШИН
E. I. Yasyukovich
SIMULATION OF CONTROLLED DIRECTIONAL MOTION OF MULTI-AXLE WHEELED MACHINES
Аннотация
Рассматриваются разработанные математическая модель и программное обеспечение для имитационного моделирования управляемого курсового движения и вертикальной динамики многоосных колесных машин. Обсуждаются полученные результаты расчетных исследований трех- и четырехосных колесных машин по моделируемой недеформируемой опорной поверхности со случайным макро- и микропрофилем.
Ключевые слова:
многоосная колесная машина, математическая модель, имитационное моделирование, микро-и макропрофиль опорной поверхности, вертикальная динамика, управляемое курсовое движение, курсовая устойчивость, управляемость.
Для цитирования:
Ясюкович, Э. И. Имитационное моделирование управляемого курсового движения многоосных колесных машин / Э. И. Ясюкович // Вестник Белорусско-Российского университета. - 2021. - № 3(72). -С. 108-116.
Abstract
The developed mathematical model and software for simulation of controlled directional motion and vertical dynamics of multi-axle wheeled vehicles are considered. The results of computational studies of three- and four-axle wheeled vehicles on a simulated non-deformable supporting surface with a random macro- and microprofile are discussed.
Keywords:
multi-axle wheeled vehicle, mathematical model, simulation modeling, micro-and macro-profile of the bearing surface, vertical dynamics, controlled directional movement, directional stability, controllability.
For citation:
Yasyukovich, E. I. Simulation of controlled directional motion of multi-axle wheeled machines / E. I. Yasyukovich // The Belarusian-Russian university herald. - 2021. - № 3(72). - P. 108-116.
Введение
Повышение эффективности многоосных грузовых колесных машин, обладающих высокой маневренностью, курсовой устойчивостью и управляемостью, требует рационального выбора
©.Ясюкович!::). И., 2021
их массогеометрических параметров, а также упругих характеристик подвески и шин. В связи с этим была поставлена задача разработки имитационной математической модели управляемого курсового движения с учетом вертикальной динамики назван-
ных машин при движении их по дорогам с микро- и макропрофилем.
Представлены основные характеристики разработанных математического и программного обеспечений виртуальных испытаний управляемого курсового движения многоосных грузовых колесных машин по недефор-мируемой опорной поверхности. Программное обеспечение (ПО) позволяет оценить влияние основных массогео-метрических параметров, упруго-диссипативных характеристик подвески и шин колесной машины на ее курсовую устойчивость и управляемость. При моделировании маневров разработанное ПО дает возможность выбирать рациональные скорости поступательного движения для обеспечения движения колесной машины без отрыва шин ее колес от опорной поверхности и без их скольжения.
Расчетные схемы и математическая модель управляемого курсового движения с учетом вертикальной динамики многоосных колесных машин
Для построения имитационной математической модели исследуемой механической системы были разработаны расчетные схемы управляемого курсового движения многоосной колесной машины (рис. 1) и ее вертикальной динамики (рис. 2).
В обозначениях параметров и независимых координат моделируемой механической системы используются двухзначные индексы, первые цифры которых относятся к ее левому борту, а вторые - к правому.
Вывод уравнений управляемого курсового движения и вертикальной динамики многоосной колесной машины выполнен на основе математической схемы Лагранжа второго рода.
Разработанная имитационная математическая модель управляемого курсового движения исследуемой меха-
нической системы содержит динамические уравнения и уравнения кинематических связей эластичных шин колесных движителей с опорной поверхностью.
Для моделирования управляемого курсового движения колесной машины на каждом шаге интегрирования уравнений ее движения необходимо вычислять значения боковых реакций Yij опорной поверхности на шину каждого управляемого колеса. Эти реакции можно вычислить с использованием различных теорий качения деформируемой шины по недеформируемой опорной поверхности, например таких, как теория Келдыша, Грейдануса, Рокара.
Вычисление реакций Yij производилось на базе теории Рокара, согласно которой боковая реакция шины определяется на основе ее угла увода 5у и коэффициента сопротивления боковому уводу kuij [1].
Поскольку коэффициенты kuij зависят от вертикальной нагруженности шин, то математическая модель курсового движения колесной машины разработана с учетом вертикальной динамики ее подрессоренной и непод-рессоренных масс - движителей.
Приведенная на рис. 1 расчетная схема курсового движения исследуемой механической системы содержит следующие геометрические параметры: li, lmc - расстояния от центра масс колесной машины до центров ее i-й оси и до проекции точки О на ее продольную ось; 0ii, 0i2...0ni, 0n2 - углы поворота управляемых колес i-й оси; dkii, dki2 - левая и правая половины ширины i-й колеи; Vii, Vi2...Vni, Vn2 - продольные скорости перемещения центров левого и правого колес i-й оси.
Для вывода динамических уравнений вертикальной динамики колесной машины разработана расчетная схема (см. рис. 2), которая содержит следующие упругие и диссипативные характеристики: cii, Ci2... Cni, Cn2 - коэффициенты нормальной жесткости шин;
от, сй2...сш, от — коэффициенты тангенциальной жесткости шин; Ор/1, Ор/2...Орп1, Орп2 - коэффициенты жесткости элементов подвески; к/1, ка...кп1, кп2 — коэффициенты нормального демп-
В качестве независимых координат моделируемой системы были выбраны следующие: Хо, уо, го — перемещения центра масс колесной машины по продольной 0х, поперечной 0у и вертикальной 0г осям; г/1, га...гп1, гп2 —
фирования шин; кт, ки2...кш, кш — коэффициенты тангенциального демпфирования шин; крп, кр/2...крп1, крп2 — коэффициенты демпфирования элементов подвески.
X
вертикальные перемещения центров неподрессоренных масс (колес); ф, у, Ф — курсовой угол, углы бокового крена и тангажа остова; 5/1, 5/2...5п1, 5п2 — углы увода шин /-й оси.
У
Рис. 1. Расчетная схема курсового движения многоосной колесной машины
^1,2
Рис. 2. Расчетная схема вертикальной динамики многоосной колесной машины
Для вычисления углов 5//, определяемых как разность между углами поворота ободьев колес вокруг их осей шкворней и плоскостью, проходящей через среднюю линию шин, были составлены уравнения неголономных связей шин с недеформируемой опорной поверхностью [5].
Таким образом, разработанная математическая модель многоосной колесной машины представлена в виде системы обыкновенных дифферен-
циальных уравнений (1), которая объединяет динамические уравнения курсового движения - три первых уравнения; вертикальной динамики - четвертое и пятое уравнения; продольно- и поперечно-угловых колебаний подрессоренной массы - шестое и седьмое уравнения, а также последнее восьмое уравнение, описывающее неголономные связи /-го колеса /-й оси с недефор-мируемой опорной поверхностью.
п 2
X = -
81П (ф ± 0//+5/ )(( - Р ) (Ф ±0/ +5/ )} т;
¿=1 / =1 п2
Ус = И7 ссз(ф±9// + 5//)) + (( -Рч)п(ф±0// + 5//)}/т;
/=1 /=1 п2
ф = У(ЦС05( -5//) + 51п( - 5/))(( -Р/) х
/=1 /=1
х (Ц 81П (0// - 5//) (Ц 008 (9// - 5/) ± ^ 008 (9// + 5//)}/ /;
((/ )
п2
¿с = ЦР/ т;
/=1 /=1 п2
^ = Р^ - РА)//;
/=1 /=1
п2
т
-, /= 1,4, / = 1,2;
9 = Ш( - РК)//;
/=1 /=1
5// К 008(Ф + 0/ - 5// ) + ус 81П (ф +0/ - 5/ ) + Ф(1 81П (0/ -5/ ) +
Л/ 008(0/ - 5/ ))} = . 81П (ф + 0/ - 5/ ) Ус 008 (ф + 0/ - 5/ )
ф (( 008(0// - 5//) + 81П (05/ ))-ф2(181П (0/ - 5// ) +
008(0/ - 5/ )) + (Ф + 0) ( 008(ф + 0/ - 5/ ) + Ус 81П(Ф + 0/ -5/ ) + Ф (( 81П (05/ ) + 008(0/. - 5/ )), 7= l, n, / = l,2, / = l, п.
+
+
(1)
В системе уравнений (1) приняты следующие параметры и обозначения: т, т/ - подрессоренная масса колесной машины и неподрессоренные массы ее /-го колеса /-й оси; /у, /г - моменты инерции остова колесной машины относительно ее продольной, поперечной и вертикальной осей; Ь - расстояние от центра масс до /-й оси колесной машины; Лку - половина ширины /-й колеи /-й оси; кс - расстояние от поверхности дороги до центра масс; Р/ - вертикальное усилие в /-м элементе подвески /-й оси; Рк/ - вертикальное усилие в /-й шине /-й оси; Ра/ - боковое усилие в /-й шине -й оси; Рс - центробежная
сила, действующая относительно продольной оси, проходящей через центр масс колесной машины, Рс = ту2 / гг.
Технология разработки математи-
ческих моделей курсового движения колесных машин подробно рассмотрена в [4, 5].
Боковые реакции опорной поверхности, действующие в пятнах контактов колес движителей с опорной поверхностью и определяющие характер курсового движения колесной машины, вычислялись по формуле
У = км.. 5..,
/ / V'
/ = 1,2, / = 1...п.
(2)
Коэффициенты кщ изменяются в зависимости от приложенных к колесам нормальных и тангенциальных сил, давления воздуха в шинах, характеристик дорожной поверхности по условиям сцепления и т. д. [1, 5]. Поэтому при моделировании курсового движения исследуемой механической системы на протяжении длительного маршрута
значение коэффициентов кт необходимо уточнять на каждом шаге интегрирования уравнений движения.
Математическая модель (1) разработана для случая отсутствия бокового проскальзывания шин при движении по криволинейной траектории, который требует, чтобы нормали к проекциям
средних линий шин каждого колеса на опорную поверхность пересекались в одной точке О (см. рис. 1). Это обусловливается следующими уравнениями, определяющими углы поворота каждого управляемого колеса относительно его оси шкворня:
0/1 = агс tan
ГI
(/+1)/2 1аП0/-1,1 ^
0 / 2 =
-111(/+1)/20/-1,1
/ 2 ,2___2
11 0/-1,1 +1( /+1)/2^П 0/-1,1
/ = 3, 5, ..., п — 1;
(3)
0Я = агс1ап
Г 1 г/2 ^0,-1,1 Л 1г/2 + ^ ^^-,1
/2 А
'//2°/-1,2
и 2
( 1/2 + ^ tan0/-1,2 ) С°§2 0,2 +1,2 Sin2 0/-1,2
/ = 2, 4,
п.
(4)
Алгоритм решения задачи исследования курсового движения многоосной колесной машины сводится к численному интегрированию уравнений движения (1) и содержит следующие блоки: ввод исходных данных и начальных условий интегрирования; считывание и привязка к реальным условиям движения ординат неровностей макро-и микропрофиля опорной поверхности; вычисление скоростей изменения этих неровностей; считывание параметров управления курсовым движением автомобиля и вычисление углов поворота каждогоу-го колеса /-й оси.
В процессе интегрирования уравнений движения разработанное ПО формирует файл результатов моделирования, содержащий параметры: моменты времени; ординаты неровностей микропрофиля опорной поверхности и их скорости; значения обобщенных координат моделируемой системы и их скорости; значения углов увода каждого колеса; нормальные и боковые реакции опор-
ной поверхности на шины колес и др.
Алгоритм предусматривает также управление продольной скоростью движения многоосной колесной машины при выполнении маневров и фиксацию момента начала бокового скольжения ее шин.
В качестве исходных данных при проведении расчетных экспериментов, кроме вышеперечисленных, использовались параметры закона изменения скорости угла поворота его переднего левого управляемого колеса, который задается с помощью специальной таблицы [5]. В первой строке таблицы указываются моменты времени начала действия заданных во второй строке скоростей изменения угла поворота переднего левого управляемого колеса, в третьей — его предельно допустимое значение, а в четвертой — предельно допустимая скорость движения колесной машины при совершении маневра, используемая при моделировании курсового движения по заданной траектории.
На основе этих данных вычислялись углы поворота каждого управляемого колеса в соответствии с геометрией поворота многоосной колесной машины, что позволило организовать моделирование различных режимов движения, таких как: прямолинейное движение, движение по круговой траектории, смена полосы движения, обгон и др.
Расчетные эксперименты управляемого курсового движения проводились на интервале времени от 0 до 250 с и более по дорогам с различными параметрами макро- и микропрофиля опорной поверхности.
Результаты моделирования выводились в файл на жесткий диск компьютера в виде численных значе-
ний: моментов времени; обобщенных координат и их скоростей; боковых реакций опорной поверхности на движители колесной машины; углов увода шин и др. По результатам имитационного моделирования формировались графические зависимости изменения во времени различных параметров курсового движения исследуемой многоосной колесной машины.
Один из вариантов моделирования курсового движения с начальной скоростью движения 25 м/с по горизонтальной дороге со случайным микропрофилем по траектории «Смена полосы движения» на интервале времени от 0 до 10 с приведен на рис. 3.
Рис. 3. Результаты имитационного моделирования курсового движения четырехосной колесной машины по траектории «Смена полосы движения» на начальной скорости 25 м/с: 1 — скорость движения; 2 — курсовой угол; 3 — перемещение центра масс по оси 0у; 4-7 — углы поворота правых управляемых колес первой, второй, третьей и четвертой осей колесной машины; 8-11 — углы увода правых управляемых колес первой, второй, третьей и четвертой осей
Из рис. 3 видно, что в начале регистрации движения в момент времени I = 0 с колесная машина двигалась со
скоростью 25 м/с (линия 1), а через 0,5 с она начала снижение скорости, т. к. через 1 с должен начинаться ее маневр —
поворот против часовой стрелки. Снижение скорости движения до 12 м/с прекратилось через 1,7 с, а в момент времени 5,4 с после завершения маневра началось повышение ее скорости до 25 м/с, и далее продолжалось прямолинейное движение по новой полосе движения (линия 3). При этом курсовой угол (линия 2) через 1 с после начала движения начал увеличиваться и достиг значения 1 рад. После 3,5 с курсовой угол начал уменьшаться и примерно при I = 9 с после начала движения достиг нулевого значения.
На рис. 3 изменение углов поворота управляемых колес отображено линиями 4-7 - левого борта, а 8-11 -
правого.
На рис. 4 представлена фазовая траектория перемещения центра масс колесной машины при совершении маневра «Смена полосы движения» на интервале движения от 0 до 460 м.
Разработанное программное обеспечение позволяет моделировать движение по криволинейной дороге, а также по траектории разворота (рис. 5).
Результаты имитационного моделирования динамики движения четырехосной колесной машины на скорости 25 м/с по дороге со случайным микропрофилем по круговой траектории представлены на рис. 6.
Рис. 4. Фазовая траектория перемещения центра масс при совершении маневра «Смена полосы движения»
NN NN. \\
\\ \\ || ¡1
2\ К 1/ // у/ уу V
-40 -20 0 20 М 40 Хс-►
Рис. 5. Фазовая траектория перемещения центра масс при совершении маневра «Разворот»:
1 - траектория перемещения центра масс; 2, 3 - левого и правого колес третьей оси
Рис. 6. Фазовая траектория перемещения центра масс колесной машины при совершении маневра «Движение по кругу»
Можно отметить, что разработанное программное обеспечение позволяет моделировать курсовое движение многоосной колесной машины по криволинейным дорогам различных категорий, что подтверждают вышеприведенные результаты.
Таким образом, разработаны математическая модель и программное обеспечение для моделирования управ-
ляемого курсового движения многоосной колесной машины со всеми управляемыми колесами с учетом вертикальной динамики, позволяющие моделировать курсовое движение по различным траекториям и категориям дорог.
Приведенные результаты моделирования подтверждают работоспособность программного обеспечения.
СПИСОК ИСПОЛЬЗОВАННОЙ ЛИТЕРАТУРЫ
1. Антонов, А. А. Теория устойчивости движения многоосных автомобилей / А. А. Антонов. -Москва: Машиностроение, 1979. - 216 с.: ил.
2. Динамика колесных машин : монография / И. С. Сазонов [и др.]. - Могилев: Белорус.-Рос. ун-т, 2006. - 462 с.: ил.
3. Литвинов, А. С. Управляемость и устойчивость автомобиля / А. С. Литвинов. - Москва: Машиностроение, 1971. - 416 с.: ил.
4. Ясюкович, Э. И. Виртуальные испытания курсовой устойчивости и управляемости четырехосных колесных машин / Э. И. Ясюкович // Вестн. Белорус.-Рос. ун-та. - 2017. - № 1. - С. 110-118.
5. Ясюкович, Э. И. Имитационное моделирование курсового движения трехосного автомобиля со всеми управляемыми колесами / Э. И. Ясюкович // Вестн. Белорус.-Рос. ун-та. - 2009. - № 4. - С. 60-67.
Статья сдана в редакцию 20 мая 2021 года
Эдвард Игнатьевич Ясюкович, канд. техн. наук, доц., Белорусско-Российский университет. E-mail: ed50@tut.by.
Edvard Ignatyevich Yasyukovich, PhD (Engineering), Associate Prof., Belarusian-Russian University. E-mail: ed50@tut.by.
