Моделирование основных режимов движения автомоби

УДК 629.113-519.711.3

И. С. Сазонов, д-р техн. наук, проф., Н. Н. Гурский, канд. техн. наук, Д. А. Дубовик

МОДЕЛИРОВАНИЕ ОСНОВНЫХ РЕЖИМОВ ДВИЖЕНИЯ АВТОМОБИЛЯ, ОСНАЩЕННОГО АКТИВНЫМИ СИСТЕМАМИ ПОДРЕССОРИВАНИЯ

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

Введение

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

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

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

В работе рассматриваются математическая модель, пространственная расчетная схема многоосной мобильной машины, принцип построения активной системы подрессоривания шасси с использованием электрогидравлического привода и результаты моделирования различных режимов движения данного объекта с помощью вычислительной среды АБМОБ [1, 2].

Математическая модель колебаний многоосной мобильной машины

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

Рис. 1. Расчетная схема пространственных колебаний многоосной мобильной машины

Математическая модель [3-5], соответствующая данной расчетной схеме, имеет вид:

Мо X-^ТСП Р - ~) + Рш + Р + Ра= 0;

1=1 І=1 гк

М2 +І І І Р]й = 0;

1=1 і=1 ] =1

JУФ +І І І1аР]а = 0;

1=1 і=1 j=l

ЛИ+І І І ЬаР]а = 0;

1=1 І=1 j=1

-ггв+І І 11аРл = 0;

1=11=1 j=6

ОТЙ 2Й - (Р Й + Р2 І1 + Р3 І1) + Р4 І1 + Р5 І1 = 0;

Іе(Ре + СГ (Ре -РР) = Ме,

Ір -Сг(Р-Ч>П) + си(р - —) = 0.

І = 1, 2; І = 1 п; j = 1, 2,3 с начальными условиями при і > і0

2 і=і0 20, 2 і=і0 20;

Р|і=і0 р0, рІі=і0 р00; и\і=і0 =и иіі=і0 = и;

в

і=і0 =в0, в

= в) •

і=і0 ^0 5

Ра

і=і0 2 0,1, 2,1

р0 а, <рп

і=і0 = 2)й ;І =1 п1 =1,2;

где М0 - общая масса; М - подрессоренная масса; Уу - момент инерции относительно оси тангажа; J— - момент инерции относительно оси крена; Уг -

момент инерции относительно оси рыскания; 1 - индекс борта; индекс номера опоры; j - индекс характеристики упругодиссипативного элемента ( j = 1 соответствует упругому элементу, j = 2 -амортизатору, j = 3 - условному эле-

менту сухого трения); Р}.и - нелинейная

характеристика восстанавливающей силы элемента ,1-подвески; 1а - расстояние по горизонтали от центра масс до оси /1-го колеса; Ъа - расстояние в поперечной плоскости от центра 1-го колеса до вертикальной плоскости, проходящей через ось крена; Ме - крутящий

момент двигателя; 1е - момент инерции двигателя; ре - угол поворота маховика двигателя; С ^ - тангенциальная жесткость ,1-й шины (на рис. 2 не показана); I 1 - момент инерции ,1-го колеса; г'а -

радиус ,1-го колеса; рп - угол поворота диска 1-го колеса;

Рш = кшр • (—)2 ,

Рш - сила сопротивления среды, кк -

коэффициент формы лобовой площади объекта; ^ - лобовая площадь объекта;

Ра = М0£ віпа

Ра - сила сопротивления подъема опорной поверхности; а - угол подъема; Рк - сила на крюке.

Относительные перемещения и скорости, являющиеся аргументами для вычисления упругодиссипативных сил в элементах ,1-й подвески Р ^ (у = 1,2,3),

имеют вид:

Д ,1 = 2р,1 2,1 ; 1 = 2р, 1 2і 1

где

2ргІ = 2 + 1г]Р + ЪгИ ; 2ргІ = 2 +1у1 + ЪгИ ,

а соответствующие аргументы для сил Р]И (у = 4,5) ,1-го колеса:

2і1 %1;

= 2п Чи.

Аргументы А.; зависят от типа и

свойств подвески ходовой части объекта. При моделировании динамической

—

системы с активной подвеской следует заменить вычисление аргументов Д^ на

Ф - ^

Д,1 = 2р,1 - 2,1 + 2,1; Ді1 = 2р,1 - 2,1 + 2а,

где 21 и 2*1 - перемещение и скорость

активного исполнительного устройства подвески.

Схема системы управления колебаниями локальной опоры мобильной машины показана на рис. 2.

Система включает контроллер 1, который может быть представлен в виде

РГО либо ИНС (искусственной нейронной сети), датчик 6, преобразователь-усилитель мощности 2, исполнительный механизм 3, соединительный механизм 4 и объект управления 5. Сигнал с датчика 6 подается в контроллер, где он сравнивается с желаемым уровнем иг перемещений подрессоренной массы и вырабатывается управляющий сигнал и, поступающий на вход электрогидравли-ческого усилителя 2 (рис. 3).

Рис. 3. Схема электрогидравлической опоры

Для моделирования пространственных колебаний мобильной машины, оснащенной активными опорами, систему уравнений (1) необходимо дополнить дифференциальными уравнениями (2), описывающими работу блоков 2, 3, 6.

ТХЯ + 41 - Ки = 0; тиЯ.2 + 42-%1 =0; т^з + 4з - киХ = 0

где

ТХ,кХ,

г I г г

ТX, к *, -

г I г I

постоянные

времени и коэффициенты усиления соответственно блоков 2, 3, 6.

Результаты моделирования

На рис. 4 приведена компьютерная модель мобильной машины с активной системой подрессоривания шасси и указаны ее основные массово-инерционные и геометрические параметры, принятые при проведении вычислительного эксперимента.

Рис. 4. Основные параметры колесной машины с активной подвеской

Для сравнительной оценки поведения машины с пассивной и активной системами подрессоривания были проанализированы следующие режимы движения исследуемого объекта:

1) переезд единичной неровности одновременно колесами левого и правого бортов;

2) переезд единичной неровности, сдвинутой во времени по отношению к колесам левого и правого бортов;

3) режим разгона и торможения с

заданным ускорением/замедлением;

4) движение на повороте;

5) движение по грунтовой дороге с последующим торможением.

Профиль единичной неровности задавался выражением

Я =

1 - соь^ж^-

Ь V.

V.

если ? >--------,

V.

где Ут - курсовая скорость машины;

Ь - длина неровности. В вычислениях принималось: Ут = 10 м/с; А = 0,1 м;

Ь = 2 м .

Режим 1. Как видно из рис. 5, активная система подрессоривания позволяет более чем в 2,5 уменьшить откло-

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

Время, с

Рис. 5. Графики продольных угловых колебаний шасси машины

Режим 2. Подвеска с активной системой позволяет снизить амплитуду поперечных угловых колебаний примерно в 5 раз (рис. 6).

Режим 3. Поведение корпуса машины при торможении с замедлением 3,7 м/с2 приведено на рис. 7. Видно, что при пассивном подрессоривании наблюдается наклон корпуса и этот угол сохраняется до останова машины. Активная подвеска удерживает горизонтальное положение корпуса в течение всего времени торможения, и лишь в начальный и конечный моменты отмечается незначительный короткий клевок корпуса.

Режим 4. Вхождение машины в поворот и движение по траектории с радиусом 12 м показано на рис. 8. Как видно, активная система подрессорива-ния удерживает корпус машины без крена и только на этапе вхождения в

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

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

Угол, град Угол, град

Угол, град

0.5 1 1.5 2 2.5 3 3.5 3.999 4.499 4.999

Время, с

Рис. 6. Графики поперечных угловых колебаний шасси машины

Время, с

Рис. 7. Графики продольных угловых колебаний шасси машины при торможении

Время, с

Рис. 8. Графики крена шасси машины при движении по криволинейной траектории

о

к

д

о

I?

о

о

л

о

С

8 10 Время, с

Рис. 9. Графики вертикальных колебаний центра масс шасси машины

15 999 17 999 19999

Л

Н

О

о

д

н

о

ц

д

д

Л

к

сз

&

<и

а

о

Частота, Гц

Рис. 10. Спектральные плотности вертикальных перемещений центра масс машины

Заключение

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

СПИСОК ЛИТЕРАТУРЫ

1. Гурский, Н. Н. Моделирование и оптимизация колебаний многоопорных машин: монография / Н. Н. Гурский, Р. И. Фурунжиев. -Минск : БНТУ, 2008. - 296 с.

2. Гурский, Н. Н. Виртуальное проектирование ходовой части мобильных машин: монография / Н. Н. Гурский, А. М. Захарик, А. М. Захарик. - Минск : БНТУ, 2010. - 174 с.

4. Смирнов, Г. А. Теория движения колесных машин / Г. А. Смирнов. - М. : Машиностроение, 1990. - 352 с.

5. Фурунжиев, Р. И. Автоматизированное проектирование колебательных систем / Р. И. Фурунжиев. - Минск : Выш. шк., 1977. -458 с.

Белорусско-Российский университет Белорусский национальный технический университет

ОИМ НАН Беларуси Материал поступил 21.04.2011

I. S. Sazonov, N. N. Gursky,

D. A. Dubovik

Modeling of basic modes of motion of vehicles equipped with active cushioning systems

This paper considers a mathematical model of the car, the chassis of which is presented in the form of concentrated masses united by elastic-dissipative elements, and the functional diagram of the system of active control of chassis vibrations. The comparative results of the simulation of the main modes of the vehicle propulsion with passive and active chassis cushioning are presented.

3. Ким, В. А. Методология создания адаптивных систем активной безопасности автотранспортных средств на основе силового анализа / В. А. Ким ; под ред. д-ра техн. наук, проф. Р. И. Фурунжиева. - Могилев : Белорус.-Рос. ун-т, 2003. - 344 с.