CC BY
55
УДК 621.015
МОДЕЛІ КОЛИВАНЬ КУЗОВА АВТОМОБІЛЯ У ПРОЦЕСІ ГАЛЬМУВАННЯ В.О. Богомолов, професор, д.т.н., В.О. Гелло, аспірант, ХНАДУ
Анотація. Запропоновано моделі коливань чотириколісної машини під дією моментів гальмових сил. Моделі дозволяють визначити вертикальні та кутові переміщення непідресореної маси в одній та двох площинах. На основі розроблених моделей виконано розрахунки у середовищі MATLAB/SIMULINK.
Ключові слова: модель, коливання, гальмові сили, комп ’ютерне моделювання.
МОДЕЛИ КОЛЕБАНИЙ КУЗОВА АВТОМОБИЛЯ В ПРОЦЕССЕ
ТОРМОЖЕНИЯ
В.А. Богомолов, профессор, д.т.н., В.А. Гелло, аспирант, ХНАДУ
Аннотация. Предложены модели колебаний четырехколесной машины под действием моментов тормозных сил. Модели позволяют определить вертикальные и угловые перемещения не-подрессоренной массы в одной и двух плоскостях. На основе разработанных моделей выполнены расчеты в среде MATLAB/SIMULINK.
Ключевые слова: модель, колебания, тормозные силы, компьютерное моделирование.
MODELS OF VEHICLE BODY OSCILLATIONS DURING BRAKING
V. Bogomolov, Professor, Doctor of Technical Science, V. Gello, post-graduate,
KhNAHU
Abstract. In the article the models considering the vehicle body oscillations under the impact of braking force torques for a four-wheel vehicle are proposed. The models can determine the vertical and rotational movements of the vehicle body in one or two planes. The simulation of the developed models in MATLAB/SIMULINK is performed.
Key words: model, oscillations, braking forces, simulation.
Вступ
Виникнення коливальних процесів, викликаних нерівномірною дією гальмових сил на колеса автомобіля, є достатньо цікавою проблемою. Якщо розглядати причини та наслідки таких процесів, то з теорії робочих процесів автомобіля стає зрозумілим, що відносно великі значення амплітуди та частоти коливань можуть вказувати на недосконалість конструкції або зношеність вузлів підвіски та гальмової системи. На автомобілях із пневматичною гальмовою системою, обладнаних регуляторами гальмових сил, такі коливання ще можуть вказувати на не-
стабільність роботи регулятора гальмових сил. У статті розглядається саме така модель, що дозволяє аналізувати коливання кузова автомобіля під дією гальмових сил.
Аналіз публікацій
Існує багато моделей коливань кузова автомобіля, що дозволяють дослідити реакцію підвіски на дорожні нерівності. Такі моделі розглянуто у [1-3]. У [1] розглянуто основні умови виникнення коливань та запропоновано одномасову модель оцінки вільних коливань автомобіля. У [2] розглянуто метод Лагранжа та наведено приклад його використан-
ня для аналізу коливань спрощеної однома-сової коливальної системи. У [3] проведено ґрунтовний аналіз поведінки автомобіля у повороті та наведено рекомендації щодо створення просторових динамічних моделей автомобіля.
Мета та постановка задачі
Метою дослідження є визначення характеру коливань кузова автомобіля під дією моментів гальмових сил та встановлення основних факторів, що сприяють виникненню таких коливань. Виникнення коливань є одночасно наслідком та причиною перехідних процесів у підвісці та гальмовому керуванні автомобіля; при цьому властивість гасити такі перехідні процеси є одним із проявів стабільності гальмового керування.
Об’єктом дослідження є автомобіль, що складається з кузова, підвіски та коліс, при цьому до коліс прикладені гальмові сили.
Моделювання коливальної системи
Розглянемо двоосний чотириколісний автомобіль як коливальну систему, що складається з підресореної маси (кузова автомобіля), непідресореної маси, пружин та демпферів.
Введемо наступні позначення (як показано нарис. 1).
Нехай 2 = 2 ^), ф = ф(ґ), Р = Р^) - відповідно вертикальне переміщення центру мас, поздовжнє кутове переміщення підресореної маси відносно центру мас та поперечне кутове переміщення підресореної маси відносно центру мас; т, Іх, Іу - відповідно підресорена маса та моменти інерції підресореної маси у поздовжньому та поперечному напрямку; к1,^,п1,п3 - коефіцієнти жорсткості та демпфірування передніх коліс; ^, ^, п2, п4 - коефіцієнти жорсткості та демпфірування задніх коліс; а, Ь - відповідно відстані від передньої та задньої осі до центру мас автомобіля; с, d - відповідно відстані від центрів правих та лівих коліс до центру мас; h - висота центру мас відносно шляху; Ть Т2, Т3, Т4 - гальмові сили, прикладені до коліс відповідно передньої та задньої осі; Мх, Му -
відповідно поздовжній та поперечний крутні моменти гальмових сил.
Величини а, Ь та с, d є динамічними параметрами, що характеризують розподіл ваги між осями. З конструктивних міркувань бажано, щоб автомобіль зберігав рівне положення за будь-якого рівня завантаженості, тому на величини коефіцієнтів жорсткості по бортах автомобіля накладаються обмеження (2) та
(3) [1].
Рис. 1. Модель коливальної системи
Припустімо, що величини коефіцієнтів ^, п, і = 1,4 по колесах підпорядковуються наступним обмеженням
ь
а'
к1 = к3, к2 = к4,
(1)
(2)
Для визначення поздовжнього крутного моменту Му припустімо, що напрямом відліку
для нього є напрям руху автомобіля. Тоді крутний момент Му дорівнюватиме сумарному моменту сил Fi на плечі h
му =1 іь
І=1
(6)
с = і
(3)
а вертикальні переміщення кожного з коліс автомобіля пов’язані з узагальненими координатами наступними співвідношеннями
21 = 2 + аф + ф,
22 = 2 + аф- dp,
23 = 2 - Ьф + ф,
24 = 2 - Ьф- dp.
(4)
Для визначення поперечного крутного моменту Мх припустімо, що момент виникає у повороті за рахунок дії бокових проекцій гальмових сил на передніх колесах, тоді для кутів повороту передніх коліс а1, а2 і неке-рованих задніх коліс момент дорівнюватиме
Мх = (Т1 sin а1 + Т3 sin а2)h.
Маємо остаточний вираз
(7)
При цьому умова (1) виконується для порожнього автомобіля, а умова (2) визначає рівність жорсткості лівого та правого бортів.
Побудуємо систему рівнянь Лагранжа для розглянутої коливальної системи. Відомо, що рівняння Лагранжа для консервативної механічної системи має вигляд
d
( дт Л
дТ
дд
дП дФ
дд д&.
(5)
.^(дОдТ дП дФ=0
_ ~Шч~~д2 +~32 + ~2~ ,
d (дТ Л дТ дП дФ
--------1-------+
дф дф
-(Т + Т2 + Т3 + Т4) h = 0,
d (дТ Л дТ дП дФ
-------1------+
др др
-(Т sinа1 + Т3 sinа2)h = 0.
(8)
Тут П - повна потенційна енергія коливальної системи; Т - її повна кінетична енергія; Ф - повна енергія, що розсіюється в демпферах; д, & - узагальнене переміщення та його перша похідна.
Для розглянутої коливальної системи узагальненими переміщеннями є 2, ф, р . Тому система складається з трьох рівнянь відносно цих переміщень, а для побудови цих рівнянь необхідно визначити енергії П, Т, Ф та їх часткові похідні по 2, ф, р. Таким чином формулюються ліві частини рівнянь системи (5).
Оскільки гальмові сили Т1,Т2,Т3,Т4 спрямовані у площині, паралельній напрямку руху автомобіля, то вертикальна збурююча сила, що впливає на переміщення кузова 2 , відсутня. Тому права частина першого рівняння (8) дорівнює нулю.
Вираз для повної потенційної енергії коливань системи П матиме вигляд
П = 2 [ ^ ^ ^ ] 2 +
+-[^а + + kзb + ]ф +
+1 [V2 + ^2 + ^с2 + k4d2 ] .р2 + (9)
+2 ф[а (^ + k2) — Ь (k3 + k4 )] +
+2 р[с (^ + kз) — d ( k2 + k4 )].
Повна кінетична енергія коливань системи Т складається з поступального переміщення підресореної маси 2 та її кутових переміщень у двох площинах з кутами ф, р. Тому маємо
2
2
(10)
Аналогічно виразимо повну енергію дисипації Ф
Ф = п1 (2+ а&+ с&) / 2 (2+ а&- d&)2/2 +
(2- Ь&+ с&)2/2 +
(2-Ь&- d|&)2 /2.
+
+п
+п
(11)
Підставивши похідні від (9-11) у (8), отримаємо шукану систему рівнянь (12-14)
т.
п1 + п2 + п3 + п4 ) 2+
+(п1а + п2а - п3Ь - п +(п1с - п2d + п3с - п
+(к1 + к2 + к3 + к4 ) 2 +
+[а (к1 + к2) — Ь ( к3 + к4 )] ф + +[с (к1 + к3) - d (к2 + к4 )]р = 0,
Іу ф+( п1а + п2а + п3Ь + п4Ь2 +(п1а + п2а - п3Ь - п4Ь ) 2+
+(п1ас - n2аd - п3Ьс + n4Ьd ) &+ +(ка + к2а + к3Ь + к4Ь )ф + +[а (к + к2) - Ь ( к3 + к4 )] 2 -
-(Т + Т2 + Т3 + Т4 ) h = 0,
(12)
(13)
спрощення розрахунків у цьому випадку можна використовувати спрощену модель коливань, що оперує тільки двома переміщеннями.
Розглянемо коливання в одній площині, що описується тільки параметрами (2, ф). Тоді система рівнянь Лагранжа матиме наступний вигляд
d (дТЛ дТ дП дФ = 0
_ ~дт~д2+~д2+ а&= ,
d (дТЛ-дТ ш 5Ф дф& дф дф дф& -(Т1 + Т2 + Т3 + Т4 ) h = 0.
(15)
Оскільки для спрощеної моделі можна прийняти, що і передні, і задні колеса жорстко з’єднані мостами, то у виразах для потенційної енергії кінетичної енергії ті енергії дисипації доцільно використовувати приведені до мостів значення жорсткості та в’язкості. Ці значення дорівнюють сумарній жорсткості та в’язкості усіх коліс одного моста (у випадку чотириколісного автомобіля - по два на міст). Відповідно вирази матимуть наступний вигляд
П — — (к + к2 + к3 + к4 ) 2 +
+ 2 [(к1 + к3 ) а 2 +(к2 + к4 ) Ь2 ]ф2 + (16)
+ [(к1 + к3 ) а-( к2 + к4 ) Ь ]( 2 + У )ф
Іхр*+( п1с2 + п2^ + п3с2 + п^2 +(п1с - п2d + п3с - п4d) 2+
+(п1ас - n2ad - п3Ьс + n4Ьd) &+ +(к^с + к2d + к3с + k4d )в + +[с (к + к3) — d ( к2 + к4 )] 2 — -(Т sinа1 + Т3 sinа2) = 0.
(14)
Під час прямолінійного руху автомобіля передні колеса направлені уздовж вектора швидкості автомобіля, тому бокова проекція гальмових сил із передніх коліс дорівнює нулю. У цьому випадку Му = 0 і коливання
у поперечній площині не виникають. Для
Ф = 0,5 (п1 + п3)(2+ аф& +0,5(п2 + п4)(2-Ь&) .
+
(17)
(18)
Взявши часткові похідні від (16-18) та підставивши їх у (15), отримаємо спрощену систему рівнянь відносно двох головних координат 2 та ф без урахування поперекових переміщень кузова Р, але з урахуванням впливу статичного прогину кузова f на вертикальні та поздовжні кутові коливання:
2
2
mZ& (n1 г n2 г n3 г n4
г ( k + k2 + k3 + k4 ) Z +
+[(n + Пз ) a _(n2 + n4 ) b ]&+ +[(k1 + k3) a _( k2 + k4 ) b ]ф = 0,
I(&-l (n + n3)a2 +(n2 + n4)b
(19)
+[(k1 + k3) a2 +(k2 + k4 ) b2 ]ф + +[(n1 + n3) a _(n2 + n4 ) b ] Z+ +[(k1 + k3) a — ( k2 + k4 ) b ] Z +
+[(k1 + k3) a — ( k2 + k4 ) b ] f —
_ (T1 + T2 + T3 + T4 ) h = 0.
На основі спрощеної моделі у середовищі MATLAB / SIMULINK було розроблено імітаційну модель та проведено розрахунки для вантажного автомобіля КАЗ-4540.
Імітаційна модель та чисельний розв’язок
Розроблена модель складається з блоків введення початкових даних, підсистем розв’язку рівнянь та блоків виводу. Розв’язок рівнянь виконується за допомогою подвійного чисельного інтегрування інтегратором ode45. Блоки введення даних реалізовані у вигляді констант, значення яких можна ввести перед стартом розрахунку. Блоками виведення є стандартний осцилограф (Scope) та вивід у простір моделі (ScopeData).
Розглянемо чисельний розрахунок для автомобіля КАЗ-4540. Згідно з [4] беремо дані для розрахунків
- для завантаженого автомобіля
а = 1.80(т),
Ь = 1.80(т), h = 1.44(т), т = 12000(к£),
/ст = 0.09(т),
V = 1.5( Нг);
- для порожнього автомобіля
а = 1.168(т),
Ь = 2.432(т), h = 1.1(т), т = 6500(к£),
/сп = 0.09(т),
V = 1.5( Нг).
Приймемо, що співвідношення виконується для порожнього автомобіля, а також приймемо масу мостів рівною 600 кг (передній міст) та 700 кг (задній міст). Тоді для порожнього автомобіля підресорена маса буде дорівнювати 4000 кг, звідки
к3 = к4 = 70729Н / т, к1 = к2 = 147271Н / т.
Екранну форму з моделлю у середовищі SIMULГNK наведено на рис. 2.
Екранні форми з отриманими графіками коливань передньої та задньої частин автомобіля, а також коливання центру мас наведено на рис. 3.
Рис. 2. Модель у середовищі SIMULINK
0.08 -0.06 -0 04 —
0.02 -
-О 02 -
01 23456789 10
Рис. 3. Графіки коливань
З останнього рисунка видно, що під час гальмування виникають коливальні процеси, пов’язані з характером взаємодії двох систем «пружина-демпфер», і характер цих процесів (амплітуда та частота коливань) залежить від співвідношень між вихідними значеннями коефіцієнтів жорсткості та демпфірування. На рисунку видно, що розраховані значення амплітуди коливань 10-12 мм для передньої осі та 15-18 мм для задньої з частотою 1,50 Гц, що є допустимим для даного вантажного автомобіля.
Висновки
У статті розглянуто моделі коливань кузова автомобіля під дією моментів, що виникають під час гальмування. Описано запропоновану модель коливань у двох площинах та спрощену модель коливань в одній площині. Для побудови моделей використовується метод рівнянь Лагранжа для консервативних механічних систем, а для комп’ютерного моделювання процесу гальмування - середовище SIMULINK.
Література
1. Литвинов А.С. Автомобиль. Теория эксплуатационных свойств / А.С. Литвинов. - М. : Машиностроение, 1989. -400 с.
2. Основы прикладной теории колебаний: учеб. пособие / А.Н. Туренко, В.И. Клименко, В.А. Богомолов и др. - Харьков: ХНАДУ, 2002. - 130 с.
3. Thomas Gillespie. Fundamentals of
Vehicle Dynamics / T. Gillespie. -Michigan: SAE Press,1997. - 640 p.
4. Грузовик КАЗ-4540. - Режим доступа : www/URL: http: //www .gruz-inform. interdalnoboy.com/description/kaz-4540-bortovoy-4x2.html - 05.05.2011р. - Загл.
з екрану.
Рецензент: В.П. Волков, професор, д.т.н.,
ХНАДУ.
Стаття надійшла до редакції 17 червня 2011 р.
