УДК 629.113.073
В. Г. МИХАИЛОВ
О НЕКОТОРЫХ ПОДХОДАХ МОДЕЛИРОВАНИЯ АВТОМОБИЛЯ
НА СИМУЛЯТОРАХ
Целью данной работы является рассмотрение и отработка подходов, методики, схем реализации моделирования автомобиля на симуляторах на основе выбора более совершенных моделей движения, колебаний, управляемости автомобиля в пакете Matlab/Simulink для получения и оценки параметров автомобиля, его на-груженности. Рассмотрены особенности имитационного моделирования автомобиля на симуляторе, требования к нему, сушествующие методы такого моделирования на основе компьютерных моделей, гидропульсаторов и движущейся платформы с мониторами/проекторами, имитирующей дорожную обстановку, где учитываются реакции на нее водителя, существенно влияющие на режимы движения и нагруженность автомобиля с дальнейшим использованием полученных данных при стендовых испытаний. Разработаны математические модели шины, движения, колебаний, управляемости автомобиля, обеспечивающие более корректное моделирование в реальном масштабе времени.
Предложены новые подходы и методика в реализации имитационного моделирования движения, управляемости и нагруженности автомобиля, базирующиеся на комбинированном использовании мощного компьютера и мини-компьютера Raspberry для считывания параметров органов управления, управления визуализацией дорожной обстановки, модуля S-Function Builder, созданного ПО на С/0++ и объединенного массива макро и микропрофиля реальной дороги, организации обмена информацией через три общих файла на SD диске Raspberry с определенной периодичностью (5 Гц), разработана блок-схема в Matlab/Simulink и ПО. Проведенным моделированием автомобиля 6х6 с ma = 21200 кГустановлено, что он имеет расход топлива 63,56 л/100 км, среднюю скорость 66,1 км/ч и число переключений КПП 7 (при t>10 c). Причем в основном переключение КПП происходит на 3-х последних высших передачах.
Ключевые слова: автомобиль, математическая модель, имитационное моделирование, симулятор, блок-схема реализации, движение, нагруженность, управляемость, продольный профиль, микропрофиль дороги, Matlab/Simulink, S-Function Builder.
Введение
Сейчас за рубежом крупные автомобильные фирмы перешли от методов компьютерного моделирования на методы имитационного моделирования на симуляторах на основе компьютерных моделей, гидропульсаторов и движущейся платформы с мониторами/проекторами, имитирующей дорожную обстановку, где учитываются реакции на нее водителя (рис. 1), существенно влияющие на режимы движения и нагруженность автомобиля с дальнейшим ис-
пользованием полученных данных при стендовых испытаний [1, 2].
Это позволяет проработать на моделях большое количество дорожных ситуаций и вариантов конструкции и испытать их в стендовых условиях, отработать эргономику, значительно сократить время создания новых моделей машин, не подвергая риску водителей и испытателей.
Особенностью таких испытаний является проведение компьютерного моделирования в ре-
Рис. 1. Имитационные испытания автомобилей за рубежом
альном масштабе времени. Визуализация дорожной обстановки обычно осуществляется на основе компьютерной анимации, требующей больших затрат на ее создание и доступно только крупным фирмам. Подобные испытания больше проводятся для отработки систем безопасности: рулевого управления, тормозных систем в экстремальных условиях и решения задач эргономики, а также используются для отладки систем управления беспилотными автомобилями и для моделирования движения и нагруженности автомобиля.
В СНГ подобных систем нет из-за больших затрат на их создание.
Для имитационного моделирования в основном используется пакет Matlab/Simulink, обладающий большими возможностями решения нестандартных ситуаций и задания внешнего возмущения.
В работах ведущей шведской расчетной компании NIRA Dynamics AB выполняющих заказы крупных автопроизводителей и являющейся крупным центром обучения [2], а также преподавателей, аспирантов, выпускников ВУЗов, прошедших стажировку на зарубежных фирмах используются известные линейные модели колебаний и управляемости автомобиля [2, 5] аналогичные применяемым в СНГ. Последние приводят к значительному расхождению по колебаниям с экспериментом (до 50— 100%).
При имитационном моделировании необходимо использование более совершенных моделей, обеспечивающих погрешность 10-15% с учетом спектра на реальной дороге.
Целью данной работы является рассмотрение и отработка подходов, методики, схем реализации моделирования автомобиля на симу-ляторе на основе выбора более совершенных моделей движения, колебаний, управляемости автомобиля в пакете Matlab/Simulink для получения и оценки параметров автомобиля, его нагруженности.
1. Задание внешнего возмущения
Для имитационного моделирования использовался полученный экспериментально массив данных продольного и микропрофиля дороги (96000 точек) [6]. Реализация возмущения с шагом 0,125 м осуществлена с использованием модуля S-Function Builder, позволяю-
щего создать пользовательскую программу на языке С/С++.
2. Реализация моделирования колебаний автомобиля
Для моделирования колебаний автомобиля использовалась плоская модель, описанная в работе [7] с применением комбинированной модели трения с упругими элементами и «вязкого» трения, которая более точнее описывает динамическую характеристику подвесок и обеспечивает хорошую сходимость в октав-ных полосах частот (погрешность 10-15%).
3. Управляемость и движение автомобиля
Исходя из анализа работ [4, 5] использовалась модель шины в виде точечного контакта, как сочетание параллельного включения упругого и демпфирующего элемента в вертикальной плоскости. Боковой увод рассматривался как параллельное включения нескольких упругих элементов с ограничением их максимальных величин [4]. Его можно описать следующими уравнениями,
F =
А увш
^уш.дин11 прк ^па ) dt Flpl > 0
t
> ^уш.дин J (^прк — ^апа ) dt > — -^гр.усил 0
t
-^гр.усил ^сли ^шу.дин1 { прк — ^па ) dt > -^гр.усил 0
t
—-^гр.усил^шу.дин! {прк — ^па ) dt < — Fo
огр.усил-
При моделировании использована следующая схема сил, действующих на автомобиль 6х6 при движении и управлении ТС (рис. 2).
На основании второго закона Ньютона движение и управляемость ТС может быть описана следующими тремя дифференциальными уравнениями:
в направлении продольной оси х
тах=соэ (5 ^)- ^ Sin (5 м)+ ^ С08 (5 №2 )-^2 sin (5^2 ) + *Хэ + ^4 + + *Хб - ^ - F»;
ах = ^+^ у;
т ах =
2^х1е°8 (5м)-2Fylsin (5^1)+4^з- Fв- ^;
а б
Рис. 2. Схема сил, действующих на автомобиль при движении (а) и управлении ТС (б)
где Fв - сила сопротивления воздуха Fв=KвАаu2; Кв - коэффициент обтекаемости автомобиля, Н с2/м4; Аа - лобовая площадь автомобиля, м2; va - текущая скорость движения автопоезда, м/с, ух - производная продольной скорости, Уу -боковая скорость; у - угловая скорость вращения автомобиля вокруг ЦТ; - сила дорожного сопротивления = 9,810та; та - полная масса автомобиля, кг; 0 - коэффициент дорожного сопротивления,
0 = ^соБа +sinа,
где а - угол наклона текущего элементарного участка маршрута, ^ - коэффициент сопротивления качению шин; Fj - инерционная сила сопротивления разгону,
р = mЛag, где а - ускорение автомобиля м/с2.
- Jeng Пт (икп ивм ) X Jk 1 =-;- + '
m„r,2 mlrrl
где © - угол поперечного наклона дороги; уу-производная продольной скорости, ух - боковая скорость; ...р^. усилия на колесах.
Вращение автомобиля относительно оси 2 вокруг ЦТ
•аЧ^- рх1 С0Б(5) -ру1 sin( 5М)-Fx2 С0Б(5^2 )-
в
ру2 (5^2) - Fx3 + Рх4 - Рх5 + ^сбЬ.
При равенстве сил Fx1 = Fx
x 2'
Fx3 = Fx 4
Fx5 = Fx6
Fy1 = Fy 2'
.5
'Л '"/Л
где l - коэффициент учета вращающихся масс, Jeng - момент инерции вращающихся масс двигателя и сцепления, кгм2, J/ - суммарный момент инерции колес автопоезда и приводных валов, кгм2.
в направлении поперечной оси y
may =(Fx1 sin (5 w1)+Fx2 sin (5w2 )+Fy1 COS (5 w1) Fy1 COS (5w2 ) + Fx2 sin (5 w2 ) + Fx2 sin (5w2 ) + Fy 2cOs(5
w2
) + Fy2 cOs (5 w2 ) + Fgr;
Fgr = -m g sin (©)
та y = (2Fxisin (5wi) + 2 Fyi cos (5wi) + Fgr;
ay = vy+vxV;
Ja Y = [-2Fx1 cos (5w1) - 2 Fy1 sin (5w1)] - .
Значение для sin углов 5w1,5w2 целесообразно заменить на 5wi (в радианах), чтобы учесть изменение знака величины. И принять углы 5w3 ... 5w6 = 0. Значение cos (5w1 2) меняются в приделах 1-0,94 при 0-10°. Их также можно заменить на 1, включая cos (5w3 6 )= 1.
Как видно из уравнений для их решения требуется определить усилия на колесах с учетом углов поворота передних колес. Они могут быть определены исходя из тяговой динамики автомобиля, описываемой известным дифференциальным уравнением второго порядка [5]
dVr,
dvx ЕЛ - F - F»
dt dt та 5
где Рк - сила тяги с учетом угла поворота управляемых колес и распределения массы на них; Fxi = К^к 5; для 1-го, 2-го колеса, Fyj = К2Fk - для 3-6-го колеса
р = Микпивм Пт .
рк =-;
г
'к
где К1 - коэффициент распределения массы на переднюю ось; К2 - коэффициент распределе-
ния массы на ведомые оси; М - текущий крутящих моментов двигателя внутреннего сгорания (ДВС), Н м, пт - КПД трансмиссии на соответствующей передаче КП; гк - радиус качения ведущих колес автомобиля, м; икп - текущее значение передаточного числа КПП; ивм - передаточное число главной передачи ведущего моста; - текущая скорость движения автомобиля, м/с,
п пг
3° UKnUBM
где n - текущая частота вращения коленчатого вала ДВС, мин-1.
Логика переключения передач КПП реализована на основе инверсированных зависимостей загрузки двигателя, предложенных в работе [5] с помощью модуля S-Function Builder. Переход на новую передачу осуществляется с задержкой t = 2 c через нейтральную, при которой F}i = 0.
5. Реализации имитационного моделирования автомобиля на симуляторе
Наиболее идеальным вариантом было бы использование схемы на рис. 1. Но она требует больших затрат. Поэтому для реализации моделирования движения и управляемости автомобиля на симуляторе разработана и использована более простая схема, представленная на рис. 3, которую в дальнейшем предполагается поэтапно усложнить.
Разработанная методика заключалось в одновременном проведении имитационного моделирования на двух компьютерах: с помощью мини-компьютера Raspberry Pi 3 B+ с платой АЦП_ЦАП ADS1256_DAC8532 и 2-х ядерно-
го основного компьютера с разделением задач. Миникомпьютер Raspberry Pi 3 использовался для визуализации дорожной обстановки (воспроизведения видеозаписи дороги) пропорциональной скорости движения и получения управляющих сигналов от педали газа, КПП и рулевого колеса. Визуализация дорожной обстановки в Raspberry Pi 3 осуществлялась с помощью разработанной программы на языке C/C++ (компилятор GCC 8.2) и библиотеки OpenCv версии 4.0.1. Обмен информации между компьютерами осуществлялся по сети с использованием файлового сервера Samba через три бинарных файла на общем (share) SD-диске Raspberry с определенной периодичностью (5 Гц) на основе счетчиков циклов, чтобы не снижать быстродействие системы. Новизной предложенного подхода является использование разработанной программы на C/C++ в модуле S-Function Builder для задания и воспроизведения массива продольного профиля и микропрофиля дороги (96000 точек), поворотов, дорожных условий, реализации многопараметрических зависимостей двигателя, учета передаточных чисел трансмиссии, нелинейных характеристик двигателя, логики переключения передач и обмена информацией, а также использование созданного ПО в Raspberry для просмотра реальной видеозаписи дорожной ситуации с переменной частотой просмотра, а не анимации [1, 2], требующей больших затрат при ее создании и соответствующего программного обеспечения..
Для макетного образца симулятора использовались руль и педаль газа от игровой приставки типа xbox 360 PS2, имеющие потенцио-метрические датчики (5,6 кОм), которые запи-
а б
Рис. 3. Общая схема (а) и блок-схема на Simulink (б) реализации имитационного моделирования автомобиля
таны напряжением 5 В. Сигналы напряжений с датчиков подаются на Raspberry Pi 3 через АЦП платы расширения High-Precision AD/ DA Board.
6. Полученные предварительные результаты
Исследование проводилось с использованием данных реальной местной дороге под Минском: Белробот-Приморье-Семков горо-док-Чучаны-МКАД2 общей протяженностью 10 км (4 км асфальтобетонного разбитого шоссе, 6 км асфальтированного шоссе хорошего качества) [8], характеристик двигателя ЯМЗ-7511 (240 кВт) для автомобиля 6х6 (ma = 21200 кГ) и КПП ЯМЗ-202 (8-ми ступенчатой) i = 8,73; 6,13; 4,51; 3,50; 2,78; 1,75; 1,00; 0,795 (i^. = 6,59), радиус качения ведущего колеса rk = 0,522 м, время переключения передачи t = 2,0 с. На данном этапе траектория движения была задана в модулях yvd_S1, yvd_S2.
Результаты моделирования приведены на рис. 4. На данном участке дороги получены следующие результаты: расход топлива 63,56 л/ 100 км, средняя скорость 66,1 км/ч и число переключений КПП (при t > 10 c) 7. Причем в основном переключение КПП происходит на 3-х последних высших передачах. Несколько повышенный расход топлива обуславливается завышенной мощностью (240 кВт) для данного автомобиля и его массы, работой двигателя
при неполной загрузке и более высоких оборотах из-за задаваемого режима высоких скоростей движения.
Проведенные исследования при скорости иа = 60-75 км/ч показали, что в этом случае можно обойтись линейной зависимостью характеристики увода. И аналогичное наблюдается при повороте при снижение иа до 25 км/ч.
Анализ показывает, что в дальнейшем целесообразно использовать проектор с экраном вместо монитора, поворачивающего изображение на экране и имитирующего траекторию дороги как это делается за рубежом.
Использование предложенной системы моделирования позволяет на стадии проектирования определить режимы движения и нагру-женности автомобиля, отработать его эргономику, оценить его параметры, влияние параметров трансмиссии на динамику, экономичность и использовать полученные данные при испытаниях на стенде.
Создание целесообразно проводить поэтапно. На первом этапе создание макетного настольного образца с рулем, педалями, управлением КПП, визуализацией, где осуществляется отработка программного обеспечения. Вторым этапом создание простого имитационного стенда на основе гидропульсаторов. А на третьем этапе движущейся платформы с кабиной. Но это по плечу только ведущим отечественным автомобилестроительным фирмам (МАЗ,
h, м О -20 -40
»1
-5
h'iS о
-20 1 а км/ч 80 60 40 20
У, ч б 10 20 30 40
V - „--
!.Ч
^--\
/ \
60 701, сек
\
\
\
ч
40 20 0 -20 -40
а*"/4 80
60
40
20
0
а)
А
А л л
\ Г \ / \r J v/ V \
V
в) Л л
\ "" V V
/—' \
/
/
в)
I—1 1—1 S 1
и.. 2000 1500 1000 500
г) \_ А ------ —
' N4 у
-V
450 t, сек
Рис. 4. Траектория перемещения ЦТ автомобиля по оси Y при выполнении маневров «переставка» (а), «поворот» (б) в зависимости от скорости и времени, профиль дороги по высоте (в), скорость движения (г), переключение передач КПП (д), обороты двигателя (е)
МЗКТ, ГАЗ, ВАЗ), имеющие гидропульсаторы Schenk или подобные им.
Заключение
1. Рассмотрены особенности имитационного моделирования автомобиля на симулято-ре, требования к нему, предложены математические модели шины, движения, колебаний, управляемости автомобиля, обеспечивающие более корректное моделирование в реальном масштабе времени.
2. Предложены новые подходы и методика в реализации имитационного моделирования движения, управляемости и нагруженности автомобиля, базирующиеся на комбинированном использовании мощного компьютера и мини-компьютера Raspberry для считывания параме-
тров органов управления, управления визуализацией и модуля S-Function Builder, созданного ПО на C/C++ и объединенного массива макро и микропрофиля реальной дороги, разработана блок-схема в Matlab/Simulink и ПО.
3. Установлено, что разработанный вариант обмена информацией через три общих файла на SD диске Raspberry с определенной периодичностью (5 Гц) позволяет не снижать быстродействие и осуществлять моделирование в реальном масштабе времени.
4. Моделирование автомобиля 6х6 с ma = 21200 кГ показало, что он имеет расход топлива 63,56 л/100 км, среднюю скорость 66,1 км/ч и число переключений КПП 7 (при t > 10 c). Причем в основном переключение КПП происходит на 3-х последних высших передачах.
ЛИТЕРАТУРА
1. Mercedes-Benz Innovation Vehicle Developing. https://www.mercedes-benz.com/en/mercedes-benz/next/advanced-engineering / [Электронный ресурс /Electronic resource] / - Режим доступа / Access mode: 22.07.2018.
2. Emanuele Obialero, A Refined Vehicle Dynamics Model for Driving Simula... rs // Charhalmers University of Technology / Göteborg, Sweden 2013. Master's thesis, P. 120.
3. Выгонный, А. Г. Закон переключения передач, оптимальный по тягово-скоростным и топливо-экономическим показателям магистрального автопоезда / А. Г. Выгонный, О. С. Руктешель, Р. Э. Шейбак// Актуальные вопросы машиностроения, 2015, Т. 4, С. 182-188.
4. Хачатуров А. А. Динамика системы дорога-шина-автомобиль-водитель / под ред. Хачатурова А. А. - М.: Машиностроение, 1976. - 535 с.
5. Hans B. Pacejka. Tyre and Vehicle Dynamics / The Netherlands, Second edition 2006. P. 642.
6. Михайлов, В. Г. Получение и использование единого массива продольного профиля и микропрофиля дороги для моделирования ТС// журнал автомобильных инженеров № 2, 2018, с. 4-7.
7. Михайлов В. Г. О колебательной модели грузового автомобиля / Д. В. Мишута, // Автомобильная промышленность. - 2016, № 7.
REFERENCE
1. Mercedes-Benz Innovation Vehicle Developing. https://www.mercedes-benz.com/en/mercedes-benz/next/advanced-engineering / [Электронный ресурс /Electronic resource] / - Режим доступа/Access mode: 22.07.2018.
2. Emanuele Obialero, A Refined Vehicle Dynamics Model for Driving Simula... rs // Charhalmers University of Technology / Göteborg, Sweden 2013. Master's thesis, P. 120.
3. Vygonnyj A. G. Zakon pereklyucheniya peredach, optimal'nyj po tyagovo-skorostnym i toplivo-ekonomicheskim po-kazatelyam magistral'nogo avtopoezda / A. G. Vygonnyj, O. S. Rukteshel', R. E. Shejbak // Aktual'nye voprosy mashinostro-eniya 2015, T. 4, S. 182-188.
4. Hachaturov A. A. Dinamika sistemy doroga-shina-avtomobil'-voditel' / pod red. Hachaturova A. A. - M. Mashinos-troenie, 1976. - 535 s.
5. Hans B. Pacejka. Tyre and Vehicle Dynamics / The Netherlands, Second edition 2006. P. 642.
6. Mikhailov V. G. Poluchenie i ispol'zovanie edinogo massiva prodol'nogo profilya i mikroprofilya dorogi dlya mod-elirovaniya TS// zhurnal avtomobil'nyh inzhenerov № 2, 2018, s. 4-7.
Поступила После доработки Принята к печати
18.02.2019 11.07.2019 01.10.2019
MIKHAILOV V. G.
ABOUT SOME APPROACHES MODELING OF THE VEHICLE ON SIMULATOR
The purpose of the given work is consideration and working off of approaches, techniques, schemes of realisation of modelling of the vehicle on simulators on the basis of a choice of more perfect models of movement, fluctuations, roadability of the vehicle in package Matlab/Simulink for reception and an estimation ofparametres of the vehicle, it loading. Features of imi-
tating modelling of the vehicle on a simulator, requirements to it, using methods of such modelling on the basis of computer models, hydropulsators and a moving platform with the monitors/projectors, simulating road conditions where reactions to it of the driver, movements essentially influencing modes and loading the vehicle with the further use of the received data at bench tests are considered are considered. Mathematical models of the tyre, movement, fluctuations, the roadabilities of the vehicle providing more correct modelling in real time are developed.
New approaches and technique in realisation of imitating modelling of movement, roadability and loading the vehicle, based on combined use of the powerful computer and minicomputer Raspberry for reading ofparametres of controls, steering of visualisation of road conditions, module S-Function Builder created ON on C/C++ both an incorporatedfile longinal profile and a microprofile of real road, the information interchange organisation through three general files on SD disk Raspberry with certain periodicity (5 Hz) are offered, the block diagramme in Matlab/Simulink and is developed software. By the spent modelling of the vehicle 6x6 with ma = 21200 kg it is established that it has fuel consumption of km of 63,56 l/100, average speed of 66,1 km/h and number of switchings of a check point 7 (at t> 10 c). And basically check point switching occurs on last 3 higher drives.
Keywords: the Vehicle, a mathematical model, simulation modeling, the unit the circuit of implementation of its driving, loading and controllability, a longitudinal profile, a road microprofile, Matlab/Simulink, S-Function Builder.
Михайлов Владимир Георгиевич. Специалист в области разработки систем CALS/PLM (PDM, ERP) PDM, ERP на Oracle на C++, PL/SQL, Java, программирования (32-х летний опыт) на C/C++, Windows, Rasbian, MC, Ar-duino/Raspberry, автомобилестроения, большой опыт моделирования динамических систем в пакете SIMULIK с использованием S-Function Builder, оценки напряженно-деформированного состояния в пакете ANSYS, испытаний подвесок, рам ТС, пневматики, гидравлики.
Mikhailov Vladimir. The expert in the field of system engineering CALS/ PLM (PDM, ERP), programming application Oracle on C++, PL/SQL, Java, programming MC, motor industry, modeling of dynamic systems in package MAT-LAB\SIMULIK (with S-Function Builder), modeling dynamic systems, vibration. estimations of the intense-deformed state in packet ANSYS, to tests of suspension, frames of the vechicle, a pneumatics, hydraulics.