The method of synthesis of the geometry of the longitudinal profile and the design parameters of the leaf spring using the finite element method
- 作者: Rubanov P.S.1, Maksimov R.O.1,2, Chetverikov M.V.1,2
-
隶属关系:
- KAMAZ Innovation Center
- Moscow Polytechnic University
- 期: 卷 91, 编号 3 (2024)
- 页面: 331-340
- 栏目: Theory, designing, testing
- ##submission.dateSubmitted##: 17.01.2024
- ##submission.dateAccepted##: 27.06.2024
- ##submission.datePublished##: 06.09.2024
- URL: https://journals.eco-vector.com/0321-4443/article/view/625745
- DOI: https://doi.org/10.17816/0321-4443-625745
- ID: 625745
如何引用文章
详细
BACKGROUND: In most cargo vehicles, leaf springs are used as an elastic element in the suspension, therefore, improving approaches to the calculation and synthesis of vehicle suspensions with leaf springs to reduce vibration load and to increase ride comfort is a relevant issue. Using the synthesis of the longitudinal profile of the spring leaves of complex shape, it is possible to achieve high strength properties of a leaf spring with a sufficiently low stiffness by applying calculations and optimizations using the finite element method (FEM), which helps to create a more perfect spring shape in terms of ride smoothness of the vehicle.
AIM: Development of a new method for the synthesis of a leaf spring with a variable profile of its longitudinal section and obtaining its characteristics using modern engineering methods based on the use of the FEM.
METHODS: The solution of the task is carried out in the NX software package including the Simcenter 3D CAE software. Topological optimization is applied to obtain the geometry of the longitudinal profile of the leaf spring, and then a strength test calculation is performed using the FEM to obtain the characteristics of the leaf spring.
RESULTS: During the work carried out at the KAMAZ Innovation Center, a method for forming the longitudinal profile of a leaf spring (regardless of the number of leaves) was developed and the dependences of the stiffness of the leaf spring on its parameters were constructed. According to the obtained dependencies, the optimal geometry of the longitudinal profile of the spring was synthesized, in which the stiffness was reduced by 33% compared to the prototype of the spring, while maintaining the bearing capacity of the vehicle.
CONCLUSION: This method of synthesis of the geometry of the longitudinal profile and the design parameters of the leaf spring can be used in the design processes of vehicle suspensions and to be further applied during research work.
全文:
ВВЕДЕНИЕ
В современных методах проектирования широко распространяется использование метода конечных элементов [1–5]. Наибольшее применение этот метод находит в проведении прочностных расчётов, с помощью которых проводится проверка работоспособности и безопасности конструктивных решений, применяемых в транспортных средствах [6–10]. Также МКЭ применяется для уточнения динамических моделей транспортных средств путём внедрения в систему динамики твёрдых тел податливых элементов [11–15] и проведения оптимизации конструкции для повышения прочностных характеристик [16–17].
Листовые рессоры являются важнейшей частью систем подрессоривания транспортных средств (ТС), которые эксплуатируются в любых климатических условиях [18]. Для исследования влияния свойств рессорных подвесок на показатели вибронагруженности и комфорта применяются динамические модели, в которых важно точное описание модели рессоры для получения достоверных результатов [19]. Известно много работ, в которых МКЭ применяется для расчёта листовых рессор на прочность. С помощью проводимых расчётов появляется возможность точно определить максимальные прогибы, возникающие в листах рессоры, распределение напряжений и спрогнозировать срок службы (ресурс) рессоры [20–21].
Целью работы является создание инновационной методики синтеза геометрии продольного профиля рессоры при помощи методов топологической оптимизации и дальнейшего получения ее напряженно-деформируемого состояния для определения зависимости характеристик рессоры от ее геометрических параметров.
МЕТОДЫ
Формирование продольного профиля рессоры
Использование топологической оптимизации в данной работе позволяет получить наиболее оптимальное распределение материала в направлении листа рессоры при условии обеспечения равной прочности по всей длине рессоры.
При синтезе геометрии продольного профиля рессоры выдвигаются следующие требования:
- достижение наибольшей мягкости рессоры для улучшения показателей вибронагруженности и комфорта водителя;
- способность рессоры выдерживать динамическую нагрузку (силу, действующую на рессору во время движения);
- сила, необходимая для удара в отбойник, должна быть меньше статической нагрузки на рессору;
- места крепления рессоры соответствуют реальным расположениям креплений в конструкции автомобиля;
- при воздействии динамических нагрузок не должен превышаться предел текучести материала рессоры;
- напряжения возникающие в листе рессоры распределены равномерно по всей длине листа.
Для получения геометрии продольного профиля рессоры необходимо провести ряд топологических оптимизаций упрощённой модели-заготовки для дальнейшего синтеза приблизительной формы листовой рессоры. В качестве исходной модели принят параллелепипед с точками крепления, расположенными в местах, аналогичных местам крепления рессоры в конструкции системы подрессоривания натурного грузового транспортного средства. Ширина выбрана аналогично ширине рессор, устанавливаемых на грузовые транспортные средства с нагрузкой на ось 9–10 тонн, а высота задана произвольно. Исходная модель для оптимизации формы рессоры представлена на рис. 1.
Рис. 1. Исходная модель для оптимизации листовых рессор.
Fig. 1. The initial model for optimization leaf springs.
Топологическая оптимизация данной модели проведена в программном комплексе NX в приложении для конечно-элементных расчётов Simcenter 3D. В первой итерации топологической оптимизации (см. рис. 2 а) видны места, где необходимо оставлять материал в модели для обеспечения меньшей массы и сохранения жёсткости. Большая часть материала, которую рекомендуется оставить, находится в нижней части модели. Материал, который находится сверху, возникает вследствие того, что происходит изгиб и возникают напряжения при растяжении в верхней части исходной модели. Следовательно, положение верхней части модели можно изменить и сместить вниз. В нижней части следует убрать материал по краям и уменьшить расстояние до нижней грани от точек закрепления, чтобы напряжения распределялись более равномерно по нижней области. Далее итерационно изменялась геометрия изначальной модели в соответствии с получаемыми результатами, как показано на рис. 2.
Рис. 2. Итерации топологической оптимизации: a) первая итерация топологической оптимизации; b) вторая итерация топологической оптимизации; c) третья итерация топологической оптимизации.
Fig. 2. Iterations of topological optimization: a) the first iteration of topological optimization; b) the second iteration of topological optimization; c) the third iteration of topological optimization.
Чтобы определить влияние количества листов рессоры на рассчитываемые параметры её жёсткости и прочности, на основе полученной топологической оптимизацией формы продольного профиля рессоры были построены упрощённые модели двухлистовой и трехлистовой рессор. На рис. 3 показаны полученные формы профилей однолистовой рессоры и упрощённые 3D-модели двухлистовой и трехлистовой рессор.
Рис. 3. Построенные с учетом топологической оптимизации 3D-модели рессор.
Fig. 3. The 3D models of leaf springs built after the topological optimization.
Расчёт параметров и характеристик рессор
Конечно-элементные модели рессор, показанные на рис. 4, сформированы из конечных элементов типа HEXA(8). В местах крепления рессор созданы точки, которые соединены с опорной площадкой рессоры RBE2 элементами (показаны синим цветом).
Рис. 4. КЭ модели листовых рессор.
Fig. 4. The FE models of leaf springs.
В проушинах рессор заданы ограничения, приложенные к узлам в серединах отверстий. Левое ограничение допускает только вращение вокруг оси Z, а правое ограничение — перемещение вдоль оси X и вращение вокруг оси Z. Между листами двухлистовой и трехлистовой рессоры задано контактное взаимодействие. Трение пренебрежимо мало, так как расчёт проводится на упрощённой модели, используемой только для получения зависимостей жёсткости и прочности от геометрических параметров рессоры. В современных конструкциях применяется ряд конструктивных решений, приводящих к снижению межлистового трения. В расчётах в качестве нагружения задаётся прогиб рессоры, для всех случаев одинаковый и соответствующий максимальному ходу (до полного сжатия отбойника) в составе передней подвески автомобиля.
Для представленных моделей проведены расчёты МКЭ для определения влияния геометрических параметров рессор на их характеристики жёсткости и прочности. Для построения зависимостей в моделях варьировались и задавались значения высоты и ширины поперечного сечения в центральной части рессоры (пакета рессор). Результаты расчётов представлены на рис. 5.
Рис. 5. Результаты расчётов рессор с помощью МКЭ.
Fig. 5. Results of leaf spring calculations using the FEM.
Для рассчитанных типов рессор построена зависимость жёсткости от количества листов рессоры, показанная на рис. 6.
Рис. 6. Зависимость жёсткости рессоры от количества листов.
Fig. 6. Dependence of spring stiffness on the number of leaves.
Для определения зависимостей характеристик рессор от их геометрических параметров проведены аналогичные расчёты, в которых варьировались значения высоты H и ширины B рессор.
В результате расчётов определялись: необходимая сила F для создания прогиба X, соответствующего полному ходу рессоры подвески, жёсткость рессоры и напряжения, возникающие в рессоре. Полученные зависимости показаны на рис. 7–10.
Рис. 7. Зависимости необходимой силы для обеспечения полного хода рессоры от высоты рессоры в поперечном сечении её центральной части.
Fig. 7. The dependence of the required force to ensure the full stroke of the leaf spring on the leaf spring height in the cross section of its central part.
Рис. 8. Зависимости жёсткости рессоры от высоты рессоры в поперечном сечении её центральной части.
Fig. 8. Dependences of the leaf spring stiffness on the leaf spring height in the cross section of its central part.
Рис. 9. Зависимости возникающих в рессоре максимальных напряжений от высоты рессоры в поперечном сечении её центральной части.
Fig. 9. The dependence of the maximal stresses arising in the leaf spring on the leaf spring height in the cross section of its central part.
Рис. 10. Зависимости жёсткости рессоры от ширины рессоры в поперечном сечении её центральной части.
Fig. 10. Dependences of the leaf spring stiffness on the leaf spring width in the cross section of its central part.
Зависимости напряжений, возникающих в рессоре, от ширины рессоры в поперечном сечении ее центральной части не строились, так как влияние ширины на напряжения достаточно мало.
На рис. 11 представлена методика синтеза геометрии продольного профиля рессоры в виде блок-схемы. Благодаря полученным зависимостям, пользуясь составленным алгоритмом синтеза формы рессоры, мы можем выбрать наиболее подходящую для транспортного средства конфигурацию, исходя из предъявляемых требований: по массе транспортного средства, на котором будет использоваться рессора; по требуемому сроку службы; по предъявляемому коэффициенту динамичности в зависимости от условий эксплуатации.
Рис. 11. Методика синтеза геометрии продольного профиля рессоры.
Fig. 11. The method of synthesis of the geometry of the longitudinal profile of a leaf spring.
По проведённым расчётам можно выделить 3 конфигурации рессоры, которые обладают наиболее выгодными характеристиками для применения на грузовых транспортных средствах с нагрузкой на ось 9–10 тонн:
- однолистовая рессора высотой 43 мм, шириной 90 мм;
- двухлистовая рессора высотой 60 мм, шириной 90 мм;
- трехлистовая рессора высотой 70 мм, шириной 90 мм.
При использовании однолистовой рессоры удаётся избежать возникновения межлистового трения, однако уменьшается коэффициент динамичности и возрастают напряжения по Мизесу. Для однолистовой рессоры достигнуто снижение жёсткости на 33% по сравнению с прототипом. При применении двухлистовой и трехлистовой рессор присутствует межлистовое трение в листах рессоры, но за счет описанных геометрических параметров рессоры удаётся уменьшить жёсткость рессоры.
РЕЗУЛЬТАТЫ
В результате проведённых расчётов получены:
- Зависимости сил полного сжатия, жёсткостей и напряжений рессоры от ее геометрических параметров.
- Характеристики жёсткости рессоры в зависимости от количества листов.
- Предложен новый метод синтеза геометрии продольного профиля и подбора требуемых характеристик рессор подвесок, применяемых на транспортных средствах.
ЗАКЛЮЧЕНИЕ
С увеличением требований, предъявляемых к системам подрессоривания транспортных средств необходимо внедрение новых методов проектирования. В работе представлен новый подход к синтезу рессор, являющихся важнейшей частью системы подрессоривания. Разработанный подход к проектированию позволяет подобрать оптимальную геометрию продольного профиля рессоры и рассчитать характеристики рессоры, а также использовать их для создания динамических моделей систем подрессоривания, применяемых в транспортных средствах для дальнейших расчётов динамики транспортных средств.
ДОПОЛНИТЕЛЬНАЯ ИНФОРМАЦИЯ
Вклад авторов: П.С. Рубанов — проведение прочностных расчетов и оптимизации, написание текста и редактирование статьи; Р.О. Максимов — обзор литературы, сбор и анализ литературных источников, редактирование статьи; М.В Четвериков — редактирование текста рукописи, экспертная оценка. Авторы подтверждают соответствие своего авторства международным критериям ICMJE (все авторы внесли существенный вклад в разработку концепции, проведение исследования и подготовку статьи, прочли и одобрили финальную версию перед публикацией).
Конфликт интересов. Авторы декларируют отсутствие явных и потенциальных конфликтов интересов, связанных с проведённым исследованием и публикацией настоящей статьи.
Источник финансирования. Авторы заявляют об отсутствии внешнего финансирования при проведении исследования и подготовке публикации.
ADDITIONAL INFORMATION
Authors’ contribution. P.S. Rubanov — carrying out strength calculations and optimization, writing and editing the text of the manuscript; R.O. Maksimov — literature review, collection and analysis of literary sources, editing the text of the manuscript; M.V. Chetverikov — editing the text of the manuscript, expert opinion. All authors made a substantial contribution to the conception of the work, acquisition, analysis, interpretation of data for the work, drafting and revising the work, final approval of the version to be published and agree to be accountable for all aspects of the work.
Competing interests. The authors declare that they have no competing interests.
Funding source. This study was not supported by any external sources of funding.
作者简介
Pavel Rubanov
KAMAZ Innovation Center
编辑信件的主要联系方式.
Email: rubanov_ps@bk.ru
ORCID iD: 0009-0000-2055-2046
SPIN 代码: 6955-1901
Design Engineer of the Engineering Calculations and Modeling Service
俄罗斯联邦, MoscowRoman Maksimov
KAMAZ Innovation Center; Moscow Polytechnic University
Email: romychmaximov@gmail.com
ORCID iD: 0009-0003-4947-790X
SPIN 代码: 7384-6758
Postgraduate of the Ground Vehicles Department; Design Engineer of the Engineering Calculations and Modeling Service
俄罗斯联邦, Moscow; MoscowMikhail Chetverikov
KAMAZ Innovation Center; Moscow Polytechnic University
Email: mihchet@gmail.com
ORCID iD: 0000-0003-3723-1171
SPIN 代码: 7949-0814
Postgraduate of the Ground Vehicles Department; Design Engineer of the Engineering Calculations and Modeling Service
俄罗斯联邦, Moscow; Moscow参考
- Khlepitko AS, Mavleev IR. Features of the use of rational methods of designing the load-bearing system of a truck. In: XII Kama readings : collection of reports of the All-Russian scientific and practical conference of students, postgraduates and young scientists, Naberezhnye Chelny, November 20, 2020 / Kazan Federal University, Naberezhnye Chelny. Naberezhnye Chelny; 2020:361–367. (In Russ). EDN: EAAWRV
- Taupek IM, Polozhentsev KA. Analysis of finite element modeling of the forging process at the RKM. In: Modern problems of the mining and metallurgical complex. Science and Production : Proceedings of the Nineteenth All-Russian Scientific and Practical Conference with international participation, Stary Oskol, December 07, 2022. Stary Oskol: NITU “MISIS”; 2023:309–314. (In Russ). EDN: DUJGZJ
- Fischev AV. Investigation of stress concentrations in nodes of metal structures. Fundamental principles of mechanics. 2023;11:64–67. (In Russ). EDN: OCNHND doi: 10.26160/2542-0127-2023-11-64-67
- Dergachev DA. Improving product quality through design automation. In: SNK-2022 : Proceedings of the LXXII open International Student Scientific Conference of the Moscow Polytechnic University, Moscow, April 04–22, 2022. Moscow: MosPolytech; 2022:90–94. (In Russ). EDN: PUFAVD
- Sutyagin AN, Kolesova VI. On the issue of specialized software performing calculations based on the finite element method. Bulletin of the Russian State University named after P. A. Solovyov. 2022;4(63):107–112. (In Russ). EDN: XKGSLE
- Gonzalez AA, Goncharov RB, Petyukov AV. Physico-mathematical modeling of the process of interaction of a passenger car airbag with an anthropomorphic mannequin. Bulletin Bauman MSTU. The Natural Sciences series. 2022;4(103):4–21. (In Russ). EDN: NJQZLI doi: 10.18698/1812-3368-2022-4-4-21
- Chetverikov MV, Goncharov RB, Butarovich DO. Investigation of the residual stress-strain state of the mini-loader load-bearing system under repeated loading according to the requirements of the ROPS safety standard. Proceedings of NAMI. 2023;1(292):46–55. (In Russ). EDN: DETBGE doi: 10.51187/0135-3152-2023-1-46-55
- Goncharov RB, Zuzov VN. Improving the structures of truck cabins at the design stage to ensure the requirements of passive safety in case of impact and minimizing mass. Proceedings of NAMI. 2019;4(279):28–37. (In Russ). EDN: XXVGQA
- Goncharov RB, Zuzov VN. Features of the search for optimal parameters of amplifiers of the rear part of the cabin of a truck based on parametric and topological optimization in order to ensure the requirements for passive safety according to international rules and obtain its minimum mass. Proceedings of the NSTU named after R.E. Alekseev. 2019;2(125):163–170. (In Russ). EDN: ZTSJEL doi: 10.46960/1816-210X_2019_2_163
- Levenkov YaYu, Vdovin DS, Alexandrov DA. Development of ROPS from aluminum alloys for front loaders. Machines and installations: design, development and operation. 2023;3:1–15. (In Russ). EDN: ALZYPP
- Vdovin DS, Chichekin IV, Levenkov YaYu. Forecasting the fatigue life of semi-trailer suspension elements at the early stages of design. Proceedings of NAMI. 2019;2(277):14–23. EDN: QTMBXK
- Gorelov VA, Komissarov AI, Vdovin DS, Chudakov OI. Analysis of the loads of a truck frame by the method of dynamics of body systems using a finite element model. Transport systems. 2020;4(18):4–14. (In Russ). EDN: GLXUZD doi: 10.46960/62045_2020_4_4
- Zhu SH, Xiao ZJ, Li XY. Vehicle frame fatigue life prediction based on finite element and multi-body dynamic. Applied Mechanics and Materials. 2012;141:578–585.
- Yudakov AA. Principles of constructing general equations of dynamics of elastic bodies based on the Craig–Bampton model and their practically significant approximations. Vestn. Udmurtsk. un-ta. Matem. Fur. Computer science. 2012;3:126–140. (In Russ).
- Goncharov RB, Ryabov DM. Methodology for calculating the loads acting in the guiding elements of the car suspension when overcoming obstacles. Izvestiya MSTU MAMI. 2015;1(3(25)):129–135. (In Russ). EDN: UXKHEZ
- Goncharov RB, Zuzov VN. Topological optimization of the car bumper design under impact from the standpoint of passive safety. Izvestiya MGTU MAMI. 2018;2(36):2–9. (In Russ). EDN: XUWXVB
- Shabolin ML, Vdovin DS. Reducing the requirements for the strength of the material of the truck subframe with independent suspension by topological optimization of the structural and power circuit. Izvestiya MGTU MAMI. 2016;4(30):90–96. (In Russ). EDN: XDEHED
- Yakovleva SP, Buslaeva II, Makharova SN, Levin AI. The influence of structural changes on the resistance to brittle fracture of the metal of the KAMAZ car springs during operation in the conditions of the North. Problems of mechanical engineering and machine reliability. 2019;3:65–73. (In Russ). EDN: WBWKBG doi: 10.1134/S0235711919030155
- Rykova OA. Issues of friction modeling in leaf springs. Proceedings of the Fraternal State University. Series: Natural and Engineering Sciences. 2022;1:207–213. (In Russ). EDN: UIEXYO
- Artyomov II, Kelasyev VV, Generalova AA. Experimental studies of destruction of leaf springs of vehicles. News of higher educational institutions. The Volga region. Technical sciences. 2009;2(10):145–155. (In Russ). EDN: KWKSHH
- Chekurda VV, Nozdrin MA. Durability of a multi-leaf car spring. In: Reliability and durability of machines and mechanisms : Collection of materials of the XIII All-Russian scientific and practical conference, Ivanovo, April 14, 2022. Ivanovo: IPSA GPS MChS Rossii; 2022;426–431. (In Russ). EDN: COUUUW
补充文件
