ft
Л ^v
ПОТОСКУЕВ Сергей Викторович
аспирант кафедры «Колесные машины» (МГТУ им. Н.Э. Баумана)
УДК 629.11.012
Повышение скорости выполнения маневров автомобиля с помощью принудительного задания соотношения скоростей вращения колес
С.В. Потоскуев
Установлена взаимосвязь угловых скоростей колес разных бортов, которая позволяет увеличить максимальную скорость автомобиля при криволинейном движении. Приведено описание математической модели криволинейного движения двухосного автомобиля с различными схемами трансмиссии. С использованием математического моделирования получены и приведены значения радиуса поворота в зависимости от скорости выполнения поворота.
Ключевые слова: поворот, управляемость, устойчивость, математическая модель, имитационное моделирование.
The interrelation of angular speeds of different sides wheels has been established, which allows to increase an automobile maximum speed in curvilinear motion. The description of a mathematical model of curvilinear motion of biaxial automobile with different transmission schemes is presented. The values of cornering radii against the speed of turning have been obtained by means of mathematical modelling and are presented.
Keywords: turning, manoeuvrability, stability, mathematical model, simulation.
Быстроходность автомобиля оценивается средней скоростью движения, на которую влияет множество факторов, в том числе кривизна пути. Для автомобиля максимальная скорость выполнения поворота на твердом основании определяется устойчивостью и управляемостью.
Наиболее распространенным способом поворота двухосного легкового автомобиля является поворот плоскостей вращения колес передней оси. В то же время с развитием электронных систем все более широко применяют системы управления продольными силами, улучшающие характеристики криволинейного движения автомобиля. Помимо известных систем динамической стабилизации, которые управляют тормозными механизмами и блокировкой дифференциалов, используют системы с межколесным дифференциалом, при полной блокировке которого колеса вращаются с различными скоростями (Mitsubishi, Honda, Magna, BMW, Audi). Это позволяет увеличить максимальные боковые реакции за счет перераспределения продоль-
ных реакций, тем самым увеличивается максимальная скорость выполнения маневра.
Для оценки эффективности подобных систем в статье рассмотрены следующие схемы трансмиссии:
— 4 х 4 с симметричными межколесными и межосевым дифференциалами;
— 4 х 4 с симметричным межосевым дифференциалом и симметричным межколесным дифференциалом на передней оси и принудительным заданием соотношения скоростей вращения колес задней оси;
— 4 х 2 с симметричным межколесным дифференциалом на задней оси;
— 4 х 2 с принудительным заданием соотношения скоростей вращения колес задней оси.
Исследования выполнены с помощью имитационного математического моделирования. Принято допущение, что движение происходит по твердой ровной горизонтальной поверхности. Поэтому движение автомобиля моделируется как плоское движение твердого тела. Для удобства составления уравнений вводится несколько систем координат: неподвижная система координат: О'х'у, подвижные системы ко-
ординат: Оху, связанная с корпусом автомобиля, и О"х"у" связанная с /-м колесом. Расчет-
7 1 I 17
ная схема движения автомобиля представлена на рис. 1.
Начало координат Оху связано с центром масс автомобиля, ось х направлена вдоль продольной оси автомобиля, ось у — вдоль поперечной. Начало координат О/х/у/ связано с проекцией центра /-го колеса на опорное основание, ось х" направлена вдоль проекции продольной плоскости колеса на основание, ось у" — вдоль проекции оси вращения колеса.
Движение автомобиля в системе координат Оху описывается следующей системой дифференциальных уравнений:
ах =
' 4 \
2** - р
\/=1
ау =
7
/ 7
2 *
\ / = 1
4
< Л =2 мпш+2 м(к /)
(1)
= 1
/ = 1
где ах, ау — проекции вектора ускорения центра масс автомобиля (абсолютной производной
Рис. 1. Схема поворота двухосного автомобиля с передними управляемыми колесами
4
от вектора скорости центра масс автомобиля); Т1 — момент инерции автомобиля относительно оси г; т — полная масса автомобиля; — угловая скорость автомобиля вокруг вертикальной оси; R( — вектор силы взаимодействия 1-го колеса с опорным основанием; Pw — силы сопротивления воздуха; Ыш — момент сопротивления повороту 1-го колеса.
Силу сопротивления воздуха Рк принято оценивать сосредоточенной силой, являющейся равнодействующей всех составляющих силы сопротивления воздуха, приложенной в точке, называемой центром парусности [1]. Поскольку скорость автомобиля в боковом направлении мала, а ветром целесообразно пренебречь, в модели учитывается только проекция силы сопротивления воздуха на ось х, Рух. Центр парусности расположен на высоте Н№ от опорной поверхности в продольной плоскости симметрии автомобиля. Величина силы сопротивления воздуха определяется по формуле Бура:
V2
Р^х = СхРРв~2~,
(2)
где сх — коэффициент аэродинамического сопротивления; ^ — площадь проекции на плоскость, перпендикулярную ее продольной оси; рв — плотность воздуха; Ух — проекция вектора скорости центра масс автомобиля на ось х системы координат Оху.
Для вычисления проекций вектора скорости центра масс V в системе Оху интегрированием необходимо определить относительную произ-~У
водную вектора-, выраженную с помощью
т
формулы Бура [2]:
йУ ~У ^ „
-Г = ~Т + ®У, ш ш
где ю — вектор угловой скорости поворота автомобиля, а
ЗУ йУх . йУу
Л Л
-I +-
сИ
1.
На основании допущения о плоском движении автомобиля значения координат вектора относительной производной от вектора скоро-
сти центра масс в системе Оху определяются по следующим формулам:
йК
сИ
йУ^
сИ
- ах + ;
■ = ау - югух ,
у проекции вектора скорости цент-
где Ух и У у
ра масс на оси системы координат Оху.
Сила взаимодействия колеса с опорным основанием определяется согласно подходу, изложенному в работе Ю.Л. Рождественского и К.Ю. Машкова [3]. Отсюда следует, что сила взаимодействия направлена в обратную сторону скорости скольжения. Скорость скольжения определяется как векторная сумма переносной Упер и относительной Уотн скоростей точки колеса, находящейся в контакте с опорным основанием и являющейся проекцией центра колеса на опорное основание:
У - У +У
(3)
Поскольку колеса передней оси автомобиля управляемые, параметры взаимодействия /-го колеса с опорным основанием удобней рассматривать в подвижной системе координат О"х"у", которая связана с /-м колесом. Рассмотрим 1-е колесо автомобиля, расчетная схема которого представлена на рис. 2.
Рис. 2. Расчетная схема сил и моментов, действующих на колесо со стороны опорной поверхности
Переносной скоростью точки /-го колеса Упер, является скорость центра /-го колеса. По-
скольку скорость какой-либо точки фигуры при ее плоском движении равна векторной сумме скорости полюса и относительной скорости этой точки от вращения фигуры вокруг полюса,
Vпер, = V + ®Рк, ,
где рк — радиус-вектор i-го колеса в системе координат Oxy. Проекции вектора переносной скорости на оси системы координат Oxy:
V =V - ra у ,
пер xi x z" к i '
Vuерy¡ =Vy + razXki ,
где xк,, ук. — координаты i-го колеса в системе координат Oxy.
Проекции вектора переносной скорости на оси системы координат O"x"y":
V^x-, ={VX - ra гук ( ) cos 0 к ( +(Vy + ra Л ( ) sin 0 ^,
Vпер у, =-{Vx - ra гук i ) sin 0 K i +(Vy + ra Л ( ) cos 0 ^ ,
где 0к — угол поворота i-го колеса относительно корпуса автомобиля.
Относительная скорость точки i-го колеса, находящейся в контакте с опорным основанием, Уотн, определяется произведением угловой скорости вращения колеса raR. и радиуса качения i-го колеса в свободном режиме rKC . Принято допущение, что радиус качения в свободном режиме постоянный и равный для всех колес. Проекции вектора относительной скорости на оси системы координат O"x"y":
V ,, =- ra r ,
отнх ¡ к ¡ ке;- '
V „ = 0.
отну ¿
Проекции вектора скорости скольжения VCK на оси системы координат O"x"y'í и его модуль определяют по следующим выражениям:
VCKX'r=(Vx -razyk¡ )cos0k¡ +
+(Vy + razXKi ) sin0Ki - rak rK0¡ ,
VcKy-i =-(Vx - razyKi) sin0Ki +
+(Vy + razX к i )c0s 0 к i ,
V = \v2 + V2
ск i у скх i ску ¡ '
Сила взаимодействия определяется нормальной реакцией и коэффициентом трения частичного скольжения:
* = Ц sЛi,
где ц^ — коэффициент трения частичного скольжения 1-го колеса; — нормальная реакция, которая действует на /-е колесо.
Коэффициент трения частичного скольжения зависит от многих параметров, но основной является зависимость от коэффициента скольжения 3к. Данная зависимость определяется экспериментально для колеса и опорного основания (рис. 3).
Рис. 3. Зависимость коэффициента трения частичного скольжения от коэффициента скольжения
Коэффициент скольжения: V
3 к
®к Гкс
Проекции вектора силы взаимодействия 1-го колеса с опорным основанием на оси системы координат О"х"у" имеет вид:
RX=- Ri cos a i; R,, =- R: sin a ¡,
где
V.
cos a:
Ve
V„
sin a =
с.ку "
Проекции вектора силы взаимодействия i-го колеса на оси системы координат Oxy:
RX, _ RX "к COS 0 к i - Ry "к SÍn 0 к i ,
Ry¡ _ Rx,,¡ sin 0 к¡ + Ry,, cos 0 к¡.
Нормальные реакции колес перераспределяются вследствие действия силы сопротивления воздуха, моментов сопротивления качению колес и ускорения центра масс:
4
2Rz¡ _ mg,
¡=1
4
4
^RíX, + 2 Mft cos 0 K¡ + PwxHw _
¡=i
¡=i
(4)
= -maxH7,
y¡ + 2 Mft sin 0к¡ =- mayH
¡=1 i _ i
где х/5 у1 — координаты /-го колеса в подвижной системе координат Оху; Н — высота центра масс автомобиля.
Предположив, что подвеска каждого колеса эквивалентна идеальной линейной пружине с жесткостью к, не сопротивляющейся поперечно, и пренебрегая динамикой колебаний корпуса, можно получить следующее соотношение:
= z + xi tg9+ y¡ tgy.
(5)
Совместное решение (4) и (5) позволяет определить значения нормальных реакций.
Связь колеса с опорным основанием не является удерживающей. Если при решении уравнений некоторые нормальные реакции получают отрицательные значения, это соответствует отрыву колес. В таком случае эти реакции принимаются равными нулю и система решается заново.
Момент сопротивления качению /-го колеса определяется выражением
Мл - , (6)
где — коэффициент сопротивления качению /-го колеса, принимается постоянным и рав-
ным для всех колес; гкЫ — радиус качения /-го колеса в ведомом режиме.
В разработанной модели приложенный к /-му колесу момент Мк. описывается выражением
(8)
Mt _ Rx rд ¡ + Mf,
где гд. — динамическим радиус /-го колеса.
Согласно [1], в среднем динамический радиус мало отличается от статического гст, поэтому принимается
r д. = Гст , = const.
Момент сопротивления повороту обусловлен скольжением отпечатка шины по опорной поверхности. При решении поставленной задачи принимается допущение, что отпечаток колеса имеет эллиптическую форму и ее изменением под действием нормальной реакции можно пренебречь. Величина момента сопротивления поворота /-го колеса определяется по следующей формуле[1]:
М„,
1 + 0,15
здесь Мт
^ пк max ■ — момент сопротивления повороту i-го колеса на месте,
Mnk max ¡ = 0,375Ц
R ЫА
s max Rz¡ V 4
где 1к — длина отпечатка шины; Rпк¡ — радиус траектории, по которой катится /-е колесо; Ьк — ширина отпечатка шины.
Координаты (ХП — смещение полюса поворота, УП = RП — радиус поворота [1, 4]) мгновенного центра поворота находятся решением системы алгебраических уравнений:
Y _ Vrn
Y П _ V
- (ХП - X К1)
' + y к1
пер y 1
V
Yn П.
Радиус траектории, по которой катится /-е колесо, рассчитывается по формуле
RnK ¡ _J(x к ¡ - X П )2 +( y к ¡ -Yп )2.
b
к
МАШИНОСТРОЕНИ
Для уменьшения бокового скольжения неуправляемых колес центр поворота должен лежать на задней оси. Если предположить, что колеса жесткие в боковом направлении, то минимальное боковое скольжение всех колес будет обеспечено следующей взаимосвязью углов и кинематического радиуса поворота:
В В
* = Ьс^ - — = Ьс^а 3 + —,
где *к — кинематический радиус поворота [1]; Ь — база автомобиля; а 1 — угол поворота /-го колеса; В — ширина колеи автомобиля.
При расчетах параметров установившегося движения принимаются следующие допущения:
1) двигатель поддерживает постоянную частоту вращения без ограничения по крутящему моменту;
2) маневр выполняется без переключения передачи в коробке передач.
Трансмиссия с колесной формулой 4 х 4 с симметричными межколесными и межосевым дифференциалами (рис. 4) описывается следующей системой уравнений:
* > ,1 = мпм + М1,
*к3® к3 =
*к2® к2 =
2
м„
2
М„.
-+м,
*к4 ® к4
2 м
2
со к1 + со
+ М,
+ М
к3
(6)
* тпС0 тп
* тз ® тз
® к2 + ® к4
2
Мс 2
Мс 2
+ со т
- М
- М.
сот
® д =■
2
* д ® д = М д - Мс
где 1Ш — момент инерции /-го колеса; сотп — угловое ускорение входного вала переднего
моста; сотз — угловое ускорение входного вала заднего моста; *д — приведенный момент инерции двигателя; *тп — приведенный к валу переднего моста момент инерции деталей, связанных с ним механически с постоянным передаточным отношением; *тз — приведенный к валу заднего моста момент инерции деталей, связанных с ним механически с постоянным передаточным отношением; Мд — момент на валу двигателя; Мпм — момент на входном валу переднего моста; Мзм — момент на входном валу заднего моста; сод — угловое ускорение двигателя.
Рис. 4. Трансмиссия с колесной формулой 4 х 4 с симметричными межколесными и межосевым дифференциалами
Трансмиссия с колесной формулой 4 х 4 с симметричным межосевым дифференциалом и симметричным межколесным дифференциалом на передней оси и принудительным заданием соотношения скоростей вращения колес задней оси (рис. 5) описывается следующей системой уравнений:
Рис. 5. Трансмиссия с колесной формулой 4 х 4 с симметричным межосевым дифференциалом и симметричным межколесным дифференциалом на передней оси и принудительным заданием соотношения скоростей вращения колес задней оси
Мпм
Т к ® к1 + М1,
М
Тк®к3 + М3,
Мзм - М2 (2 - ы) + МАи,
®4 - ® двU, ®2 - ® дв (2 - U),
® к1 + ® к3
® тп - 2 '
■ - М -
Ттп®тп 2 Мп.м ,
Т тз ®тз 2 Мз.м ,
■ ® тп + ® тз
ю д-—Г"
(7)
Тд®д - Мд - Мс,
где и — соотношение скоростей вращения левого и правого колес заднего моста.
Рис. 6. Трансмиссия с колесной формулой 4 х 2 с симметричным межколесным дифференциалом
Рис. 7. Трансмиссия с колесной формулой 4 х2 с принудительным заданием соотношения скоростей вращения колес задней оси
Трансмиссия заднеприводного автомобиля с колесной формулой 4 х 2 с симметричным межколесным дифференциалом (рис. 6) описывается следующей системой уравнений:
Тк2®к2 -
М
М
Т к4 ® к4 2
+ М2,
+ М4,
® к2 + ® к4
® д -
Т д ® д - Мд - Мз.
2
2
Г, 1 1 а) б) в) г) д) /
■......'......— .....' \ 4 4 //,.....*......~.. \ \ \ V /V: ..........,.,..............................................................................V......X/../\/...............................
: \ \ -у/// У/У/ ...... ............ ..............................*.....................................................- г/. / * ; ,/,/ у //£ / /.*
1
...................................../ />................................. //+ /А*
...................................................... /X/ X/* \ /¿^
" 1 ¡. 1
V, км/ч
Рис. 8. Функция минимального радиуса поворота Я(У) для автомобиля с различными видами трансмиссий:
а — 4 х 4 с симметричными межколесными и межосевым дифференциалами; б — 4 х 4 с симметричным межосевым дифференциалом и симметричным межколесным дифференциалом на передней оси и принудительным заданием соотношения скоростей вращения колес задней оси; в — 4 х 2 с симметричным межколесным дифференциалом на задней оси; г — 4 х 2 с принудительным заданием соотношения скоростей вращения колес задней оси; д — максимально возможная кривизна траектории, ограниченная максимальным значением коэффициента трения
скольжения (^кр = м )
Трансмиссия с колесной формулой 4 х 2 с принудительным заданием соотношения скоростей вращения колес задней оси (рис. 7) описывается системой уравнений вида
1 д ® д = Мд - Мз.м>
Мс = М2 ( 2 - и) + М4 и,
®4 = ® двU, ®2 = ® дв (2 - и).
(9)
Задаваемое соотношение определяется таким образом, чтобы обеспечить минимальное скольжение при качении колес по концентрическим окружностям:
1
и =
1-
в
2Я
(10)
Для построения функции минимального радиуса, ограниченного заносом, определяют фактический радиус поворота Я и скорость ус-
тановившегося движения Упри различных задаваемых кинематических радиусах Як = 4...110 м. По результатам обобщения данных заездов строится функция минимального радиуса поворота, ограниченного заносом, от скорости установившегося движения (рис. 8).
Основные технические характеристики автомобиля, используемые в имитационной модели, представлены в табл. 1.
Таблица 1
Технические характеристики исследуемого автомобиля
№ Показатель Значение
1 Полная масса, кг 1500
2 Момент инерции машины относительно оси г, кг-м2 120
3 База, мм 2 800
4 Ширина колеи, мм 1 630
5 Расстояние от оси переднего колеса до центра масс, мм 1 400
6 Высота центра масс, мм 400
7 Свободный радиус колеса, мм 280
8 Коэффициент аэродинамического сопротивления 0,8
Таблица 2
Максимальные скорости прохождения поворота радиусом 120 м для различных схем трансмиссии автомобиля 4 х4
Схема трансмиссии 4x4 с симметричными межколесным и межосевым дифференциалами 4х4 с принудительным заданием соотношения скоростей вращения колес задней оси
Vmaxf км/ч 93,5 97,3
Таблица 3 Максимальные скорости прохождения поворота радиусом 120 м для различных схем трансмиссии автомобиля 4 2
Схема трансмиссии 4x2 c симметричным межколесным дифференциалом на задней оси 4х2 с принудительным заданием соотношения скоростей вращения колес задней оси
Vmax, км/ч 94,7 99,3
Вывод
Из полученных зависимостей следует, что принудительное задание соотношения скоростей вращения колес задней оси для автомобилей со схемами трансмиссий 4x4 и 4x2 при движении по ровной твердой горизонтальной поверхности позволяет повысить максимальную скорость выполнения маневра в сравнении
с дифференциальной трансмиссией. При установившемся движении на высоких скоростях преимущество для полноприводного автомобиля составляет 4%, для заднеприводного — 5% (см. табл. 2, табл. 3). На низких скоростях движения (до 40 км/ч) установившиеся скорости движения одинаковы. Для более точной и полной оценки трансмиссии с принудительным заданием соотношения скоростей вращения колес необходимо провести исследование поведения автомобиля при входе на поворот и выходе из него.
Литература
1. Смирнов Г.А. Теория движения колесных машин: Учеб. для студентов машиностроительных специальностей вузов. 2-е изд., доп. и перераб. М.: Машиностроение, 1990. 352 с.
2. Добронравов В.В., Никитин Н.Н. Курс теоретической механики: Учеб. для машиностроительных специальностей. вузов. 4-е изд., перераб. и доп. М.: Высш. шк., 1983. 575 с.
3. Рождественский Ю.Л., Машков К.Ю. О формировании реакций при качении упругого колеса по недеформи-руемому основанию // Труды МВТУ. 1982. № 390. С. 56—64.
4. Фаробин Я.Е. Теория поворота транспортных машин. М.: Машиностроение, 1970. 176 с.
Статья поступила в редакцию 10.03.2011 г.