Method for determining the stiffness of arbitrary-shaped anti-roll bar

Cover Page

Cite item

Full Text

Abstract

The article studies methods for stiffness calculation of wheeled vehicle suspension systems spring elements. The KAMAZ Innovation Center has developed a new method for antiroll-bar stiffness calculation is developed, taking into account its geometric shape. This method allows to achieve a high accuracy of calculation with an error within 5 %. The design scheme of the antiroll-bar is presented. The comparative assessment is made on the basis of the results of a series of anti-roll bar models calculations using a numerical method in a finite element analysis program.

The proposed method allows to clarify the calculations of vehicle dynamics by obtaining the stiffness value of the transverse stability stabilizer close to the real one for correct determination of the roll angle of the wheeled vehicle's sprung mass, which is an important criterion of vehicle motion safety.

Full Text

Введение

Система подрессоривания – одна из самых значимых систем транспортного средства (ТС), для оценки работоспособности которой необходимо иметь точную динамическую модель автомобиля, в рамках которой будут учитываться упругие элементы. Существуют различные способы задания упругих тел в твердотельные динамические модели транспортных средств, наиболее известными являются:

  1. внедрение редуцированной конечно-элементной модели модальным методом Крейга-Бэмптона [1-3];
  2. разделение единого твердого тела на несколько частей и соединение их шарниром с заданной жесткостью [4, 5].

Стабилизатор поперечной устойчивости является одним из податливых элементов, присутствующих в системе подрессоривания. Он предназначен для уменьшения углов поперечного крена подрессоренной массы ТС. При прямолинейном движении автомобиля все участки устройства поворачиваются на один и тот же угол, и стабилизатор не скручивается, а в момент поворота появляется крен ТС и две половины стабилизатора поворачиваются на разные углы, за счет чего появляется упругий момент, создающий сопротивление крену автомобиля.

При создании динамической модели ТС возникает проблема корректного определения жесткости стабилизатора поперечной устойчивости. Современные методики предполагают рассмотрение данного узла как упругого стержня, работающего на кручение [6]. При этом если геометрия стабилизатора представляет собой сложную форму возникают трудности с применением данных методик ввиду большого расхождения результатов аналитических расчетов с экспериментальными данными. Исходя из этого возникает необходимость в разработке нового метода определения жесткости стабилизатора поперечной устойчивости произвольной формы.

Целью работы является разработка метода расчета жесткости стабилизатора поперечной устойчивости произвольной формы.

Задачи исследования:

  1.  анализ геометрии стабилизатора поперечной устойчивости;
  2.  оценка типа нагружения участков стабилизатора;
  3.  приведение жесткостей участков к общей жесткости устройства;
  4.  проверка корректности подхода методом конечных элементов (МКЭ).

Метод определения жесткости стабилизатора поперечной устойчивости 

Стабилизатор представляет из себя профиль круглого сечения с направляющей произвольной формы, одни части которого работают на кручение, а другие – на изгиб (рис. 1) [7].

 На основе расчетной схемы (рис. 1) построен метод определения его жесткости, состоящий из следующих этапов.

  1. Стабилизатор разбивается на несколько простых участков.
  2. Определяются участки работающие на изгиб (s1, s3, s5, s7) и на кручение (l2, l4, l6).
  3. Рассчитываются жесткости участков стабилизатора в соответствии с типами их нагружения.

 

Рис. 1. Геометрические параметры стабилизатора поперечной устойчивости сложной формы\

Fig. 1. Geometrical parameters of an anti-roll bar of complex shape

 

Жесткости участков стабилизатора, работающих на изгиб, определяются по формуле (1):

ci=EIxsi, (1)

где i = 1, 3, 5, 7 – номера участков стабилизатора, работающих на изгиб; E– модуль упругости первого рода для материала стабилизатора, Па; Ix – центральный момент инерции сечения при изгибе, m4; si – длины участков стабилизатора, работающих на изгиб, м.

Жесткости участков стабилизатора, работающих на кручение, определяются по формуле (2):

Cj=GIplj, (2)

где j = 2, 4, 6– номера участков стабилизатора, работающих на кручение; G – модуль упругости второго рода для материала стабилизатора, Па; Ip – полярный момент инерции сечения, m4; lj; – длина участка стабилизатора, работающего на кручение, м.

  1. Определяется суммарная жесткость стабилизатора поперечной устойчивости.

Так как все участки стабилизатора между собой соединены последовательно, суммарная жестокость определяется по формуле (3):

1cсб=1c1+1c2+1c3+1c4+1c5+1c6+1c7, (3)

где c1, c3, c5, c7 – жесткости участков стабилизатора, работающих на изгиб, Нм/°; c2, c4, c6 – жесткости участков стабилизатора, работающих на кручение, Нм/°.

Проверка работоспособности аналитического метода

Объектом исследования являются стабилизаторы поперечной устойчивости простой и сложной геометрических форм с диаметрами прутка dсб=36 мм и 45 мм соответственно.

На рис. 2 представлен стабилизатор простой геометрической формы, состоящий из трех простых участков.

 

Рис. 2. Геометрические параметры стабилизатора поперечной устойчивости простой формы

Fig. 2. Geometrical parameters of an anti-roll bar of simple shape

 

Участки с длинами s1, s3 работают на изгиб, а участок с длиной l2 на кручение.

Момент инерции круга при изгибе вычисляется по формуле (4) [8]:

Ix=πdсб464 , (4)

где dсб диметр стабилизатора.

Полярный момент инерции круга вычисляется по формуле (5) [8]:

Ip=πdсб432.(5)

Жесткости этих участков рассчитываются по формулам (1) и (2) соответственно.

Учитывая симметричную форму стабилизатора и подставляя (1), (2), (4) и (5) в основную формулу (3), получим следующее выражение (6):

1cсб=264s1Eπdсб4+32l2Gπdсб4 . (6)

Далее необходимо рассмотреть стабилизатор поперечной устойчивости сложной геометрической формы (рис. 1). Стабилизатор состоит из семи простых участков. Участки с длинами s1, s3, s5, s7 работают на изгиб, а l2, l4, l6 на кручение. Жесткости этих участков рассчитывается по формулам (1) и (2) с учетом формул (4) и (5). Учитывая симметрию стабилизатора и подставляя (1), (2), (4) и (5) в основную формулу (3), получим выражение (7) для определения жесткости стабилизатора сложной геометрической формы:

1cсб=264s1Eπdсб4+232l2Gπdсб4+264s3Eπdсб4+32l4Gπdсб4 (7)

Для подтверждения корректности метода аналитического расчета жесткости стабилизатора поперечной устойчивости, необходимо воспользоваться расчетами МКЭ [9-11].

Граничные условия для стабилизатора поперечной устойчивости во всех четырех расчетных случаях выбраны в соответствии с закреплениями в конструкции подвески натурного транспортного средства (рис. 3).

 

Рис. 3. Граничные условия: а) резинометаллический шарнир; б) ограничение по всем степеням свободы, кроме вращения вокруг оси Х; в) вращательное кинематическое воздействие относительно оси Х на 1º

Fig. 3. Boundary conditions: a) rubber-metal hinge; b) restriction on all degrees of freedom except rotation around the X-axis; c) rotational kinematic action relative to the X-axis by 1º

 

Результаты расчетов жесткостей стабилизатора поперечной устойчивости при помощи МКЭ для всех расчетных случаев представлены на рис. 4 и 5.

 

 Рис. 4. Распределение моментов реакций для стабилизатора поперечной устойчивости простой геометрической формы: а) стабилизатор dсб =36 мм; б) стабилизатор dсб = 45 мм

Fig. 4. Distribution of reaction moments for the anti-roll bar of simple geometric shape: a) anti-roll bar dсб = 36 mm; b) anti-roll bar dсб = 45 mm

  

Рис. 5. Распределение моментов реакций для стабилизатора поперечной устойчивости сложной геометрической формы: а) стабилизатор dсб =36 мм; б) стабилизатор dсб = 45 мм

Fig. 5. Distribution of reaction moments for the anti-roll bar of complex geometric shape:  a) anti-roll bar dсб = 36 mm; b) anti-roll bar dсб = 45 mm

 

Численные значения результатов расчетов жесткостей стабилизаторов поперечной устойчивости аналитическим методом и МКЭ представлены в табл. 1.

 

Таблица 1.Результаты расчетов жесткостей стабилизаторов поперечной устойчивости

Table 1. Results of stiffness calculations of anti-roll bars

 

Аналитический расчет

МКЭ

Диаметр стабилизатора dсб, мм

36

45

36

45

Жесткость стабилизатора простой геометрии Cсб, Нм/°

203

496

211

515

Жесткость стабилизатора сложной геометрии Cсб, Нм/°

220

538

215

525

 

Из табл. 1 видно, что отклонения результатов аналитических расчетов от результатов расчета численным методом составляют менее 5 %, что свидетельствует о корректности проведенных расчетов и работоспособности разработанного метода определения жесткости стабилизатора поперечной устойчивости.

Заключение

Разработан аналитический метод определения жесткости стабилизатора поперечной устойчивости произвольной геометрической формы. Проведен анализ геометрии и оценка типа нагружения участков стабилизатора. Определена общая жесткость устройства. Отклонения результатов, полученных при помощи нового метода, от результатов численного расчета составляют не более 5 %, что подтверждает корректность метода определения жесткости стабилизатора поперечной устойчивости произвольной формы.

×

About the authors

M. K. Мавлонов

KAMAZ Innovation Center LLC

Author for correspondence.
Email: mirzoodil.mavlonov@mail.ru
ORCID iD: 0009-0005-9193-4996
Russian Federation, Moscow

M. V. Chetverikov

KAMAZ Innovation Center LLC

Email: mihchet@gmail.com
ORCID iD: 0000-0003-3723-1171
Russian Federation, Moscow

References

  1. Рубанов, П.С. Оценка влияния учета податливости рамы фронтального погрузчика на возникающие нагрузки в системе динамики твердых тел / П.С. Рубанов, Р.Б. Гончаров, Г. И. Скотников и др. // Известия МГТУ МАМИ. 2023. Т. 17, № 4. С. 401-409. – doi: 10.17816/2074-0530-472077. – EDN ZYYKEM.
  2. Vdovin, D. Prediction of fatigue life of suspension parts of the semi-trailer in the early stages of design / D. Vdovin, Y. Levenkov, I. Chichekin // IOP Conference Series: Materials Science and Engineering: Design Technologies for Wheeled and Tracked Vehicles, MMBC 2019, Moscow, 01–02 октября 2019 г. Vol. 820. – M.: Institute of Physics Publishing, 2020. P. 012002. – doi: 10.1088/1757-899X/820/1/012002. – EDN HNLDCE.
  3. Горелов, В.А. Анализ нагрузок рамы грузового автомобиля методом динамики систем тел с использованием конечно-элементной модели / В.А. Горелов, А.И. Комиссаров, Д.С. Вдовин, О. И. Чудаков // Транспортные системы. 2020. № 4(18). С. 4-14. – EDN GLXUZD.
  4. Максимов, Р.О. Виртуальный стенд для определения нагрузок в пневматической подвеске задней тележки грузового автомобиля на ранних стадиях проектирования / Р.О. Максимов, И.В. Чичекин // Известия МГТУ «МАМИ». 2021. №3 (49). С. 76 – 86. doi: 10.31992/2074-0530-2021-49-3-76-86.
  5. Левенков, Я.Ю. Определение параметров модели рессоры для анализа нагрузок и оценки прочности элементов подвески в системе расчета динамики твердых тел / Я.Ю. Левенков, И.В. Чичекин // Инженерный вестник. 2016. № 12. С. 4. – EDN XICZAL.
  6. Афанасьев, Б.А. Проектирование полноприводных колесных машин: учебник для вузов. Т. 3 / Б.А. Афанасьев, Б.Н. Белоусов, Л.Ф. Жеглов и др. – М.: МГТУ им. Н.Э. Баумана, 2008. – 432 с.
  7. Тягунов, А.Е. Расчет и обоснование конструктивных параметров автомобильных механических стабилизаторов поперечной устойчивости с регулируемой жесткостью / А.Е. Тягунов, А.Б. Карташов // Известия МГТУ «МАМИ». 2017. Т. 11. №2. C. 64-71. – doi: 10.17816/2074-0530-66918.
  8. Феодосьев, В.И. Сопротивление материалов: учебник для вузов. – М.: МГТУ им. Н.Э. Баумана, 1999. – 592 с.
  9. Рубанов, П.С. Совершенствование конструкции несущей системы полуприцепа на стадии проектирования путем параметрической оптимизации / П.С. Рубанов, Р.Б. Гончаров, А.О. Поляков, И.О. Хренов // Труды НГТУ им. Р.Е. Алексеева. – 2024. – № 3(146). – С. 108-119.
  10. Petyukov A.V., Goncharov R.B., Gonsales Astua A.V. Numerical simulation of airbag module operating features. AIP Conf. Proc. 6 October 2023; 2833 (1): 020043. https://doi.org/10.1063/5.0151985
  11. Рубанов, П.С. Методика синтеза геометрии продольного профиля и конструктивных параметров листовой рессоры с применением метода конечных элементов / П.С. Рубанов, Р.О. Максимов, М.В. Четвериков // Тракторы и сельхозмашины. 2024. Т. 91. № 3. C. 331-340.

Supplementary files

Supplementary Files
Action
1. JATS XML
2. Fig. 1. Geometrical parameters of an anti-roll bar of complex shape

Download (90KB)
3. Fig. 2. Geometrical parameters of an anti-roll bar of simple shape

Download (68KB)
4. Fig. 3. Boundary conditions: a) rubber-metal hinge; b) restriction on all degrees of freedom except rotation around the X-axis; c) rotational kinematic action relative to the X-axis by 1º

Download (182KB)
5. Fig. 4. Distribution of reaction moments for the anti-roll bar of simple geometric shape: a) anti-roll bar dсб = 36 mm; b) anti-roll bar dсб = 45 mm

Download (190KB)
6. Fig. 5. Distribution of reaction moments for the anti-roll bar of complex geometric shape: a) anti-roll bar dсб = 36 mm; b) anti-roll bar dсб = 45 mm

Download (205KB)

Copyright (c) 2025 Мавлонов M.K., Chetverikov M.V.

Creative Commons License
This work is licensed under a Creative Commons Attribution 4.0 International License.

СМИ зарегистрировано Федеральной службой по надзору в сфере связи, информационных технологий и массовых коммуникаций (Роскомнадзор).
Регистрационный номер и дата принятия решения о регистрации СМИ: серия ПИ № ФС 77 - 56417 от 11 декабря 2013.