УДК 629.113
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ДЛЯ ИССЛЕДОВАНИЯ УСТОЙЧИВОСТИ ДВУХОСНОГО АВТОМОБИЛЯ С ИНДИВИДУАЛЬНЫМ ПРИВОДОМ ПРИ РАЗГОНЕ НА ПРЯМОЛИНЕЙНОМ УЧАСТКЕ
С.Н. Марченко, А.В. Келлер, Г.Д. Драгунов
Разработана математическая модель для исследования устойчивости автомобиля с индивидуальным приводом при разгоне максимальной интенсивности и обоснован рациональный характер распределения мощности между ведущими колесами двухосного автомобиля при разгоне на прямолинейном участке.
В настоящее время появились конструкции различных систем распределения мощности, которые позволяют регулировать крутящие моменты на каждом ведущем колесе в соответствии с условиями движения, независимо от условий движения других колес. Такие конструкции создали базу, благодаря которой стал возможен дальнейший качественно новый этап в развитии теории и разработке более совершенных оптимальных схем привода.
Новый этап может быть сформулирован как создание управления эксплуатационными свойствами путем разработки систем распределения мощности с задаваемыми управляемыми характеристиками. Развить данное направление теории простым суммированием результатов исследований, уже выполненных отдельными авторами, невозможно в силу следующих причин:
- положения существующей теории ориентированы на использование традиционных для крупносерийного машиностроения технических решений: постоянно включенный дифференциальный или подключаемый блокированный привод, использование дифференциальных механизмов, часто имеющих различную схему блокировки и т. д.;
- исследования по отдельным разделам теории мобильных колесных машин носят причинно-следственный характер. Это означает, что в методологическом плане такие исследования относятся к работам анализа, поскольку содержащиеся в них рекомендации по выбору характеристик и параметров машин получены путем сопоставительного анализа.
Одним из важнейших вопросов при разработке новой конструкции системы распределения мощности является обоснование правильного распределения моментов, которое можно охарактеризовать коэффициентом распределения крутящих моментов км. Не трудно убедиться, что большее значение км обеспечивает повышение тяговых качеств автомобиля за счет рационального использования сцепного веса ведущего моста. А это, в свою очередь, обеспечивает лучшую динамичность и проходимость автомобиля в условиях различных сил сцепления ведущих колес с дорожной поверхностью. Однако значительная величина км может оказать отрицательное воздействие на устойчивость, управляемость и топливную экономичность автомобиля.
Следовательно, степень распределения моментов нужно ограничивать, исходя из условий, при которых указанное отрицательное воздействие минимально или отсутствует.
Очевидно, что влияние степени распределения моментов на поведение автомобиля наиболее существенно при движении на повороте, а также при прямолинейном движении в случае трога-ния с места и разгоне, или движении автомобиля на подъем в условиях, когда силы сцепления ведущих колес с поверхностью различны.
Чтобы получить рациональный характер распределения крутящих моментов между колесами двухосного автомобиля, на основе уравнений Лагранжа второго рода были составлены дифференциальные уравнениями движения материальной системы (рис. 1) по трем степеням свободы X, У, Д
™ ■ Jx' = +2^1 8‘п01 ’
т • У/ ~тн(Ь~Ьн)/3 = Куг +2008^ “ 2вшвх, (1)
= -Гх2)0,5В2+(Я'х1 -*",)(),5^0*4 ~ХЯУ2Ь-
~Хкуіасо^і ~2Хіа8іп0і мсш> где т - полная масса автомобиля; тн - масса неподрессоренных частей; Jz - момент инерции автомобиля вокруг вертикальной оси, проходящей через центр масс О'; Я‘х/ - сила тяги на колесах автомобиля; К'У] - боковые реакции на колесах автомобиля; В} - ширина колеи соответствующего моста; ЕМСШ - суммарный стабилизирующий момент на колесах; а, Ъ - координаты центра масс в продольном направлении; Ън - расстояние центра неподрессоренных масс до задней оси; в - среднее значение угла поворота управляемых колес.
Ускорения центра масс по соответствующим осям определяют из кинематических соотношений:
Л' = ох. - оуР, уу = иу. - охф, где и] - скорость центра масс по соответствующей координате; /3' - производная по углу разворота автомобиля.
Для определения величин действующих вертикальных реакций Е!2] в месте контакта колес с дорогой составим основные уравнения движения центра подрессоренных масс по двум степеням свободы а, ф. Пренебрегая в этих уравнениях членами второго порядка, получим:
Jnкyа + 2(гр1а2к + гр2Ъ2к )а + 2(С\а2к + С2Ь2к )а + М%+ Мпуу (Ик —гк) = 0,
(Зпюс + МгЯр № + куЧ + (Су- Спккр )\у - Jnhв$ - МД^у'у +М%= О,
где У1щу - момент инерции подрессоренных масс относительно поперечной оси, проходящей через центр масс Ок; М„ - подрессоренная масса автомобиля; гр— коэффициент сопротивления амортизатора; Ькр - плечо крена кузова; Л// - момент сухого трения в подвеске в продольном направлении; - момент инерции подрессоренных масс относительно продольной оси, проходящей через центр масс Ок; - центробежный момент инерции подрессоренных масс; гк - радиус колеса; аК, Ьк, Ик - координаты центра тяжести подрессоренных масс; С, - приведенная жесткость рессоры в вертикальном направлении; - приведенная жесткость подвески в поперечном направлении; к,:, - приведенный коэффициент демпфирования в поперечном направлении
к¥ - 0,5^ гр1^1/ ; ~ момент сухого трения подвески в поперечном направлении
В случае пробоя передней подвески дифференциальное уравнение движения (2) по координате а запишется в виде:
зпкуа + Р2^1& + 2{СшХа\ + С2Ъ\)а - 2а к (Сш1 - С\)1{ + 2^р2^к + Мп.1х'0гк ~гк)~ (3)
где Г\- ход передней части кузова до пробоя рессоры по оси передних колес; Сш - вертикальная жесткость тины переднего колеса.
Уравнения (2) и (3) позволяют определить усилия, передаваемые через элементы подвески автомобиля на подрессоренные части. Учитывая, что угловые колебания кузова записаны от положения статического равновесия, определим вертикальные составляющие на колесах автомобиля в процессе разгона:
+0,5 К%±м$,
где $щст - величина статической нагрузки на колесе; Е.\- - динамическая нагрузка, передаваемая через рессору при дифференте кузова; АВУ% -динамическая нагрузка, передаваемая через рессору при крене кузова.
После этого получим:
11
= 2ак (Сха + г ^й), при а < —.
ак
К
Д*1 = 2ак\Сш1а -{Сш1 -С^^-З, при а >
ак ак
= 2Ък(С2а + гр2а),
ДЛ^=0,5(Л'.-Л^) =
С,„у/ + КщуУ + Мру ы§а\(г
Траектория движения автомобиля во время разгона определяется наличием углов увода переднего и заднего мостов, которые возникают под действием боковых сил К у , приложенных в пятне контакта колес с дорогой. Используя кинематические соотношения, определим величины боковых реакций на каждом колесе автомобиля для малых углов:
К-у\ -кбУїіУхіів + Оу/
иу> + оза
соЬ-оу'
Ку2 =к5у'г2у1Х2
где д-ф - угол разворота колес управляемого моста вследствие крена кузова; кд - значение коэффициента сопротивления увода колеса при оптимальной нагрузке; - корректирующие ко-
эффициенты.
В зависимости от величины силы тяги и нагрузки на колесе коэффициент сопротивления боковому уводу изменяется по известному закону [Л]. Значения корректирующих коэффициентов в расчетах определяли по следующим зависимостям:
Л2
У* =2,4
К
П
-1,8
Яі-
+ 0,4
К,
^3
я
%}0]0
V %
где Яф„т - оптимальная нагрузка на колесо; ф - коэффициент сцепления под колесами соответствующего борта автомобиля.
Для определения законов изменения основных силовых факторов, действующих на автомобиль в процессе разгона с индивидуальным приводом ведущих колес, найдем основные парамет-
ры траектории движения центра масс, которые получим при решении системы дифференциальных уравнений (1). Для удобства систему приведем к виду со = А ; ох. =£ -Vуб); 0у> = В - охгсо,
Ед.,асоз9, 2Х», (д;,-д;,)д. , (д;,-д;,)в,
и
где А
Б = В =
Д =
J J 2J
Z Z Z
IX +2Хcos^ -IXsin^,.
т
+IX,cos^ -IX.sin^..
т
м},Ф-К)
J
м
со - j3.
Тогда координаты траектории движения центра масс автомобиля можно определить как t t t /З - Jcodt; X = j(ux' cos /3 - uy> sin /3)dt; у = j[oy cos [3 + vx' sin [3)dt,
0 0 0 где x, у - координаты центра масс механической системы на плоскости; (3 - угол разворота продольной оси.
В процессе разгона рост сил тяги Kxj (рис. 2) на ведущих колесах сопровождается дифферентом кузова а и увеличением вертикальных нагрузок на колеса задней оси.
16000 4,0 5,0
R, Н . У м/с " jx, м/с2
/ -4,0
12000 - / 3,0-
R ~ ■3,5
10000 - у/ 2 5 і -3,0
8000- Кху' ~~~~ 2,0- -2,5
6000- 1,5- -2,0
4000- 1,0 - -1,5
/ .їх -1,0
2000- 0,5- -0,5
0 V— 1 1 1 1 0 1 1 1 1
0 1
t, с
0
1
г, с
Рис. 2. Изменение параметров автомобиля при разгоне
Вследствие различных сил тяги на ведущих колесах в горизонтальной плоскости возникает поворачивающий момент, который обусловливает появление боковых реакций на колесах и углов увода, а кроме того, приводит к развороту передних управляемых колес за счет податливости деталей рулевого привода, что способствует появлению дополнительного угла разворота управляемых колес вру, который в данном случае выступает как дестабилизирующий фактор.
Отмеченные выше явления приводят к началу движения автомобиля по криволинейной траектории (рис. 3). Движение по криволинейной траектории сопровождается действием на подрессоренные массы поперечной СИЛЫ инерции за счет возникновения бокового ускорения /(,. Это приводит к крену кузова и соответствующему перераспределению вертикальной нагрузки (см. рис. 2) в сторону колес с большими действующими силами тяги, что в свою очередь при работе системы управления приведет к еще большему увеличению на них сил тяги. Следовательно, при наличии неравномерности коэффициентов сцепления ведущих колес имеем возмущенный характер движения автомобиля. При этом в процессе разгона автомобиля возможно достижение предельной силы сцепления на колесах переднего моста, что может привести к потере ими способности воспринимать боковые усилия без скольжения и потери вследствие этого управляемости. Начавшееся регулирование крутящих моментов на ведущих колесах уже в переходной ста-
дии может обеспечить автомобилю требуемый уровень управляемости при высокой динамике разгона.
Снижение величины начальной неравномерности сил тяги на ведущих колесах способствует уменьшению отклонений центра масс от направления движения. Однако при больших значениях неравномерности коэффициентов сцепления ведущих колес выравнивание характеристик процесса заканчивается при некотором положительном угле /3 разворота продольной оси, под которым впоследствии автомобиль продолжает движение. В случае отсутствия корректирующего воздействия со стороны водителя в данной ситуации возможен выход автомобиля из своего ряда движения.
Рис. 3. Изменение параметров, характеризующих траекторию автомобиля при разгоне
Чтобы определить крайние точки автомобиля при теоретическом анализе, воспользуемся координатой смещения центра масс в боковом направлении у и углом разворота продольной оси /3. С учетом габаритных размеров запишем отклонения крайних точек для передней и задней частей кузова:
ух - у + аг sin /? + 0,5Вг cos j.3 sign /3, у2 = у - Ъг sin /3 - 0,55г cos /? sign /?,
где аг, Ьг - расстояние от центра масс соответственно до передней и задней части кузова; Вг - ширина кузова автомобиля.
Тогда значение коэффициента устойчивости определим из отношения разности величин допускаемой полуширины полосы дороги (назовем ее шириной коридора безопасности) и максимального отклонения крайних точек кузова автомобиля в процессе разгона от осевой линии к разности этих же величин в исходном положении при нахождении центра масс на осевой линии полосы:
_ 0>5Да - W
Vi~ 0,5{Вд-В г)' где Вд - ширина коридора безопасности.
При этом потеря устойчивости наблюдается при т\,<0, причем степень неустойчивости характеризуют отрицательные величины коэффициента 7]j.
Таким образом, разработанная математическая модель позволяет определить рациональный характер распределения крутящего момента по ведущим колесам и оценить его влияние на устойчивость движения автомобиля.
Литература
Антонов Д.А. Теория устойчивости движения многоосных автомобилей. — М.: Машиностроение, 1978. - 216 с.