Дьяков И.Ф., д-р техн. наук, проф.
Ульяновский государственный технический университет
ОПТИМАЛЬНЫЙ ВЫБОР ПОДВЕСКИ ТРАНСПОРТНЫХ СРЕДСТВ
В статье произведен анализ системы подрессоривания транспортных средств. Основными характеристиками подрессоренной системы являются частоты собственных колебаний и коэффициенты их относительного затухания.
Цель данной работы - выбор подвески по критериям оптимальности на стадии проектирования и испытания в дорожных условиях. При анализе подвески для упрощения расчетов транспортное средство рассматривается как плоская модель. Степень защиты подрессоренной массы от вибраций оценивается по критериям оптимальности и коэффициентом виброзащиты.
Ключевые слова: подвеска, спектр частот, передаточная функция, момент инерции балансира, коэффициент демпфирования, виброзащита, торсион._
Состояние вопроса. Современные транспортные средства оснащены различными системами подрессоривания [1, 3, 5], для которых характерны динамические взаимодействия различного рода упругих и диссипативных сил. При движении все массы автомобиля совершают более или менее интенсивные линейные и угловые колебания в различных плоскостях. Возбуждение колебаний происходит в результате взаимодействия автомобиля с дорогой, а также из-за неравномерности работы двигателя и агрегатов трансмиссий [2]. Для гигиенической оценки вибрации, действующей на человека используют интегральный показатель и раздельно-частотный [4].
Исследованиями установлено, что ускорения элементов автомобиля при его движении имеют широкий спектр частот (свыше 500 Гц). Весь спектр частот ускорений можно разделить на два диапазона: низкочастотный ( до 25 Гц) и высокочастотный (свыше 25 Гц). Низкочастотные колебания в основном вызываются взаимодействием автомобиля с дорогой. Их уровень зависит от условий движения, геометрических параметров транспортного средства и его подвески. Высокочастотные колебания обусловлены неравномерностью работы двигателя и агрегатов трансмиссии. Они зависят от скорости движения (частоты вращения трансмиссии и нагрузки двигателя), степени удаленности рассматриваемого элемента от источника возмущения и конструкции элементов, передающих вибрации. Колебания в высокочастотном диапазоне практически не зависят от микропрофиля дороги и параметров подвески.
Сопоставление расчетных и экспериментальных спектров ускорений показывает, что они достаточно близки только в низкочастотном диапазоне. Любая расчетная схема включает инерционные элементы (массы), соединенные упругими элементами, а также элементы, обес-
печивающие рассеивание энергии при колебаниях масс.
Массы элементов, связывающих подрессоренные и неподрессоренные части (упругие элементы, рычаги направляющего устройства, амортизаторы, тяги рулевого привода и карданный вал), относят частично к подрессоренным и неподрессоренным массам. Исследованиями установлено, что при расчетах колебаний в низкочастотном диапазоне все подрессоренные массы могут быть объединены в одну массу Ш с моментом инерции относительно поперечной оси, проходящей через центр масс, и — относительно продольной. Неподрессоренная масса каждого моста Ш рассматривается отдельно и считается сосредоточенной.
Основными упругими элементами в расчетных схемах являются различные подвески и шины. При расчетах плавности хода принимается, что все эти элементы расположены в плоскостях колеса, а их упругие свойства оцениваются приведенными характеристиками.
Учитывая, что в наибольшей степени плавность хода определяется колебаниями в вертикальной плоскости, в ряде случаев при ее анализе для упрощения расчетов рассматривают только эти колебания. При этом пространственная модель автомобиля заменяется плоской, в которой совмещаются правые и левые подвески и колеса мостов, а высоты неровностей дороги q считают равными полусумме высот неровностей, находящихся в определенный момент времени под правым и левым колесами моста
Ч = 0,5 (Чпр + Чш )' где Чпр;, — высота неровностей дороги соответственно под правым и левым колесами \ -го моста.
Расчетные схемы подвески, применяемые при проектировании транспортных средств в
продольной плоскости приведены на рис. 1.
Рис. 1. Расчетные модели, применяемые при анализе плавности хода транспортных средств
На схеме показаны: подрессоренная масса то с моментом инерции относительно поперечной оси; неподрессоренные массы щ и ш2 упругие элементы подвесок с жесткостями 2ср1 и 2с 2 ; амортизаторы с сопротивлениями 2&р1 и 2£ 2; силы трения и ; упругие элементы с жесткостями 2сш1 и 2сш2 и условные амортизаторы с сопротивлениями 2£ш1 и 2кш2, моделирующие шины; Ь — база автомобиля. На колеса действуют неровности дороги высотой ^ и
д2, которые возбуждают вибрации масс. На стадии проектирования по выбранной схеме подрессоривания транспортного средства в рамках обобщенного подхода предлагается оценка межпарциальных связей рычажной системы подвески. Примем, что транспортное средство
движется с постоянной скоростью, поэтому массы автомобиля в продольном направлении не имеют ускорений. В продольно-вертикальной плоскости подрессоренная масса совершает перемещение по координатам Z и ф , а неподрес-
соренная - по координатам ^ и . Перемещение подрессоренной массы можно рассматривать и в других координатах, например в координатах Zl и Z2 (перемещения точек (А{ )
подрессоренной массы, которые лежат над передним и задним мостом соответственно).
Решение. На первом этапе проектирования подвески плавность хода транспортного средства можно рассматривать при условии, когда обеспечивается минимум передаточной функции, представляющей собой отношение потенциальной энергии неподрессоренной массы к подрессоренной в виде
* ( ^ )
( ¿нУ
= е
2п
л тг
с
V у
V, 1 (115 + 112 ) — ^2 (122 + 125 )
<
2 (Zп + V Л — |>п +(V л — ф,зАб )])
[* ],
(1)
где 2 — отношение неподрессоренных масс к
подрессоренным; , — скорости вертикальных перемещений соответственно не-подрессоренной массы и подрессоренной; /тр
— коэффициент трения межлистовой рессоры; Сш, Ср — вертикальные жесткости соответ-
ственно шины и рессоры; V, V — угловые перемещения соответственно переднего и заднего
и—
моста; V —угловые перемещения рамы;
расстояние от центра тяжести до точки опоры соответственно /-й массы (/ =1...3, 1- передний, 2-задний мосты и 3 - рама транспортного средства).
Согласно выражению (1) по геометрическим параметрам подвеску можно подобрать так, чтобы не возникал резонанс колебаний собственной частоты и частоты колебаний подрес-
соренной массы от неровности дорожного полотна. На рис. 2 приведена амплитудно-фазо-частотная характеристика угловых колебаний рамы транспортного средства.
Рис. 2. Амплитудно-фазо-частотная характеристика угловых колебаний рамы
Как видно, рама при движении по неровностям дороги имеет четыре колебательных контура при переходе из одной частоты в другую.
При движении со скоростью V неровности дороги изменяются по закону х = vt, а возмущение, передаваемое от неровностей дороги на транспортное средство по закону косинусоиды
9 = 9 [1 - со§ {уг)] , где V — частота возмущения, связанная со скоростью движения автомобиля и длиной неровности; I — длина неровности дороги.
Циклическая частота связана с путевой соотношением X = 2п /1. Путевую и циклическую частоты можно представить как частоты возмущений, передаваемых на транспортное средство от неровности дороги при его движении со скоростью 1 м/с: путевую — в радианах в секунду, циклическую - в герцах. Циклическая (путевая) частота X соответствует частоте возмущений при скорости движения 1м/с. Если транспортное средство движется со скоростью V, частота возмущения V при тех же длинах неровностей увеличивается в V раз. В общем виде v = XV, что равносильно удлинению оси абсцисс графи-
ка и увеличению площади под кривой спектральной плотности, равной дисперсии высот неровностей в п раз. Компоновки современных транспортных средств таковы, что для большинства конструкций горизонтальной оси координат а и Ь относительно центра тяжести (см.рис. 1) транспортного средства принимают условие аЬ = р 2 и соблюдается с погрешностью не более 10 %. Это позволяет рассматривать колебания передней и задней частей транспортного средства независимо друг от друга. Поэтому в ряде случаев, особенно когда необходимо проследить влияние характеристик отдельных элементов подвески на интенсивность колебаний, можно использовать расчетную схему, приведенную на рис. 1. Расчетные схемы многоосных транспортных средств отличаются от рассмотренных только наличием большего количества мостов. Если на транспортном средстве применяется балансирная рессорная подвеска, в расчетных схемах необходимо учитывать момент инерции
балансира. Возможность выбора балансир-
ной подвески из трех практически равнозначными схемами показана на рис. 3.
Рис. 3. Моделирование разновидной балансирной подвески
Как видно на первой схеме (рис.3, а упругий элемент установлен между подрессоренной массой и жестким балансиром, момент инерции которого относительно оси качения равен моменту инерции реального балансира. На второй схеме (рис.3, б) жесткий балансир соединен с подрессоренной массой шарнирно, а с неподрес-соренной — упругими элементами с жесткостью, составляющей половину жесткости рессоры балансира. На третьей схеме (рис.3, в) колеса связаны жестким безынерционным балансиром, а масса балансира распределена между подрессоренной и неподрессоренными массами так, чтобы моменты инерции реальных и расчетных неподрессоренных масс и балансира относительно оси качания были одинаковы. В некоторых случаях при исследовании низкочастотных колебаний подрессоренных масс в диапазоне ниже 5 Гц можно не учитывать инерционные воздействия неподрессоренных масс. В этом случае расчетная схема двухосного автомобиля может быть представлена в виде балки, установленной на двух пружинах с амортизаторами. Жесткость пружины определяется суммарной
жесткостью рессор С и шин С
Основными
(^ + 1) Д (^ + 2 ) Д (<
характеристиками любой колебательной системы являются частоты собственных колебаний и коэффициенты их относительного затухания. Они связаны с частотным диапазоном возмущений, при которых систему можно рассматривать как динамическую модель.
Вычисление частот собственных колебаний проводим в определенной последовательности: составляем уравнения движения системы (по схемам без затухания и возмущений): разделяем каждое уравнение на инерционный коэффициент; представляем уравнения в операторной форме; составляем определитель системы, заменив в нем оператор дифференцирования s на мнимую частоту Ю, и приравниванием его нулю; раскрытием определителя и решение его относительно О . Вычисление определителя для системы с п степенями свободы сводится к решению уравнения, имеющего п корней (значений квадратов частот собственных колебаний).
Уравнения без членов, характеризующих затухание и возмущения в системе представляем в операторной форме:
ш
51
+пД - = 0; :+ пД - п?252 = 0;
) 51 - у ,1 ^ = 0;
(2)
(^2+ш52) 52 - у2 = 0;
где шz1 =
2
Ш22 =
(Ь2+р2) т;
2 (ср2 + сш2 ),
т0
•2 •П!' =(Ь! + Р2) т„ •
После замены оператора (2) дифференцирования на мнимую частоту получаем следующий определитель системы, вычислением которого
ш2 -О2
Ш й =
2c.lL1
2Ср 2 £
(а2 + р2) то'
аЬ - р\ 2 2(
п, L; ш21 =—
ср1 + сш1
Ь + р2
2cp2L2
2с
. 2 _ . _ ~~р1 .
' П £2 = ? ~ ' ^ ,1 ="
(Ь2 + Р2) то'
т,
т
2с
У ,1 = — т2
могут быть получены частоты собственных колебаний системы
п,
п,
,1
0
ш2 2 -О2
,2
После аналогичных преобразований определитель системы уравнений (3)
ш2 -О2 -п, ш 2 -О2
-п5 0
ш^ -О2 0
п52
ш -О2
= 0.
(3)
Решая уравнение (4), получим частоты собственных колебаний системы:
п5
= 0
(4)
О2 = 0,5
(ш2 + ш^ )±^(ш2 - ш^ )2 + 4п,п
где ш2 = 2ср/тЛ; =2(ср + сш)т; п, = 2ср/то1;
щ
2ор / т
Если принять перед корнем знак плюс будем иметь более высокую частоту
0
0
0
а =. 0,5
(«^+«2 )+^( «2—«2 )2+
, (5)
при знаке минус - низкую
а=. 0,5
(«^ + «2 ) — «2 — «2 )2 +
.(6)
При этом связанную двухмассовую систему можно рассматривать две одномассовые парциальные. Парциальную систему можно представить как многомассовую, в которой все массы закреплены, кроме одной, которая может совершать движение по одной координате. Из выражений (5), (6) видно, что при наличии связи меняется частота собственных колебаний - увеличивается более высокая и уменьшается низкая.
На втором этапе проектирования проводится исследование плавности хода транспортного средства в дорожных условиях. Критерием оптимальности подвески принимаем время успокоения £ , за которое амплитуда колебаний
уменьшится до заданной минимальной величины:
£усп
1
е У «о
-1п
1
X,(1 — еУ)
^ Ш1П ,
где еу = ¿усп / (2тп«о ) = ¿усп / () —
степень успокоения; ¿усп — коэффициент успокоения, численно равное силе сопротивления подвески при скорости подрессоренной массы, равной единице (2 = 1); ю0 = /тп — коэффициент, зависящий от упругих свойств подвески; Сп — жёсткость подвески; тп — подрессоренная масса; Хк — коэффициент точности установки подрессоренных масс относительно центра тяжести автомобиля.
Результаты.
Результаты показали, что нагрузочные характеристики подвески имеют слабо выраженную нелинейность и нельзя получить достаточную точность оценки плавности хода транспортного средства при различных нагрузках на ось, так как полезная нагрузка меняется в широких пределах. В такой подвеске для подрессори-вания транспортного средства, применяют тор-сион рабочее тело, которого выполнено в виде конической формы, а соотношение наибольшего и наименьшего диаметров выбраны таким, чтобы обеспечить нагрузочную характеристику подвески с прогрессивной характеристикой, при
которой частоты малых собственных колебаний подрессоренных масс могли остаться неизменными, независимо от загрузки транспортного средства (изохронная подвеска). Расчеты параметров торсиона и определение нагрузочной характеристики задней подвески транспортного средства подтверждают возможность создания изохронной подвески. Для получения наивысшей плавности хода требуется обеспечить равенство частот собственных колебаний подвески при любой загрузке транспортного средства
300
V = У2 или, как следует из формулы V =- ,
11
где /1 = /2— прогиб подвески при частотах
V = V собственных колебаний. Это достигается подбором диаметров правого ^ и левого й2 торсиона. Степень защиты подрессоренной массы от вибраций часто оценивается коэффициентом виброзащиты у, равным отношению амплитуды колебаний подрессоренной массы к амплитуде неподрессоренной. Виброзащита будет удовлетворительной при у =1. Это условие выполняется, если отношение частоты колебания рамы автомобиля «р к частоте колебания
неподрессоренной массы «н будет выражено
неравенством. Обычно это отношение находится в пределах от 2,5 до 5. Если отношение частот соответствует неравенству « / «н < 1,4,то
виброзащита ухудшает условия эксплуатации, так как в этом случае у >1 и амплитуда
юр /шн > 1,4, а колебания подрессоренной массы больше амплитуды колебания неподрессо-ренной массы. При низких частотах амортизаторы плохо предохраняют от ударов подрессоренную массу, поэтому обычно выбирают амортизаторы такими, чтобы коэффициент демпфирования находился в пределах 0,25...0,5. В этих пределах может быть торсионная система подрессоривания. Для исключения влияния случайных факторов на результаты испытаний проводилось по три заезда на каждой скорости движения автомобиля.
На рис. 4 приведена запись осциллограммы вертикальных колебаний подрессоренных масс автомобиля УАЗ-3303.
При обработке данных вычислялось среднее арифметическое и среднеквадратическое значения вертикальных виброускорений. Измерение уровней виброускорений относительно нулевого порога в вертикальном направлении проводилось при работе двигателя без нагрузки с частотой вращения коленчатого вала 800 и 3600 мин-1.
Рис.4. Осциллограмма вертикальных колебаний автомобиля при У=30км/ч
Наиболее эффективна система защиты агрегатов от вибрационного воздействия со стороны рамы, на которую установлен двигатель является опора крепления, которая должна иметь частоту f в 1,5 ...2 раза ниже частоты возбуждения, чем все собственные частоты f виброизолированного объекта. Исследования показали, что не всегда удается добиться того, чтобы в одном или двух направлениях при при
n = 800 ...1000 мин 1 получить f /f > 2 . В
этих случаях необходимо стремиться к тому, чтобы резонанс не возникал в тех направлениях, в которых действуют основные возмущающие факторы. При рассмотрении вынужденных колебаний с одной степенью свободы принимают условие, что внешнее усилие, вызывающее колебания, изменяется гармонически по закону синуса F sin rot. Кроме внешней возмущающей силы в системе действует восстанавливающая сила c z, стремящаяся вернуть систему в положение равновесия. При наличии трения в опорах крепления массы свободные колебания быстро затухают. В системах с вынужденными колебаниями они будут продолжаться с амплитудой
z = F sin rot / (c — mro2).
Максимальную амплитуду этих колебаний находили из выражения
zm* = F / (c _ mro2 ) = F / c f 1 _ m
m 7
--ro
c
Отношение F / c представляет собой максимальное перемещение гст, вызванное статическим действием силы F . Обеспечение данного условия при создании конкретных конструкций виброизоляторов значительно осложняется рядом особенностей работы автомобильного транспорта, а зачастую - практически невоз-
можно. Во-первых, двигатель - один из основных источников шума и вибрации, как правило, имеет переменный режим работы. Отношение максимальных рабочих скоростей вращения коленчатого вала двигателя к минимальным составляет 6,4 и боле. Во-вторых, виброизолированный двигатель является динамической системой со многими степенями свободы и совершает сложные пространственные взаимосвязанные колебания. Кроме того масса двигателя с коробкой передач в два и более раза больше массы несущей системы, а податливость последней приводит к статическому прогибу ее до нескольких миллиметров. В связи с этим упругое основание существенно влияет на характер колебаний кузова, и это влияние необходимо учитывать. Если принять во внимание еще конструктивные и технологические ограничения, то создание приемлемой конструкции виброизолятора является технически сложной задачей. В силу отмеченных особенностей транспортных средств практически не представляется возможным избежать резонансных режимов. При этом задача сводится к решению вопроса о выборе оптимального демпфирования в упругой опоре или применению других методов (например, с переменной жесткостью, с ограничителем хода виброизолированного агрегата) для ограничения амплитуд вынужденных колебаний на резонансе.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Дьяков И.Ф. Основы оптимизации в автомобилестроении. М.: Машиностроение, 2012. 380 с.
2. Жеглов Л.Ф. Виброакустика колесных машин. Изд. МГТУ им. Н. Э. Баумана, 2013. 170 с.
3. Жеглов Л.Ф. Спектральный метод расчета систем подрессоривания колесных машин. Изд-во МГТУ им. Н. Э. Баумана, 2009. С. 150.
4. Ружичка Дж. Дерби Т. Виброизоляция с нелинейным демпфированием /Пер. с анг. К. А. Выгодского, ОНТИ № 609, 1973. 42 с.
5. Раймпель Й. Шасси автомобиля: конструкция подвески / Пер. с нем. В. П. Агапова. М.: Машиностроение, 1989. 328 с.
Dyakov I.F.
OPTIMAL SELECTION OF VEHICLE SUSPENSION
The article analyzes springing system of vehicles. The main characteristics of the springing system are frequency of natural oscillations and their relative settlement coefficients.
The purpose of this work is choosing the suspension according to the criteria of optimality at the stage of designing and testing on road. When analyzing the suspension to facilitate the calculation the vehicle is considered as a plane model. Protection of the sprung mass vibration is evaluated according to the criteria of optimality and vibration protection factor.
Key words: suspension, frequency spectrum, transfer function, moment of inertia of the balance, the damping coefficient, vibration protection, torsion bar.
Дьяков Иван Федорович, доктор технических наук, профессор кафедры основы проектирования машин и инженерная графика.
Ульяновский государственный технический университет. Адрес: Россия, 432027, Ульяновск, ул. Северный Венец, д. 32. E-mail: [email protected].