FSAE赛车悬架的设计与研究.pdf

第2 0 期AUTOMOBILE APPLIED TECHNOLOGY设i计研究10.16638/ki.1671-7988.2023.020.011FSAE赛车悬架的设计与研究安嘉晨，臧利国*，陈博文12，钱程1，柏宇星1（1.南京工程学院汽车与轨道交通学院，江苏南京2 1116 7；2.北汽重型汽车有限公司，江苏常州2 13133）摘要：基于大学生方程式汽车大赛（FSAE）比赛规则，分析了比赛对赛车悬架性能的要求，计算了相关的动力学参数，确定了悬架系统的形式并使用CATIA进行建模。在ADAMS中构建了悬架的运动学模型，并分析了悬架模型在轮跳分析时相关性能指标参数的变化范围。仿真结果表明，所设计的悬架系统满足规则要求，但车轮定位参数随车轮上下跳动时的变化幅度较大，不利于操纵稳定性。然后在ADAMS/Insight 中基于车轮的定位参数对模型的硬点坐标进行试验设计（DOE）分析优化，优化前后的结果对比表明，优化后车轮跳动时车轮定位参数的变化量显著减小，赛车的性能得到了很好的改善。关键词：FSAE；悬架系统；车辆运动学仿真中图分类号：U469.6+96AN Jjachen,ZANG Liguo,CHEN Bowen2,QIAN Cheng,BAI Yuxing(1.School of Automobile and Rail Transit,Nanjing Institute of Technology,Nanjing 211167,China;2.BAIC Heavy Truck Company Limited,Changzhou 213133,China)Abstract:Based on the formula society of automotive engineers(FSAE)competition rules,therequirements of the race on the suspension performance of the racing car are analyzed,the relevantdynamic parameters are calculated,the form of the suspension system is determined and modeled inCATIA.In ADAMS,the kinematics model of the suspension is constructed,and the variation rangeof the relevant performance index parameters of the suspension model during the wheel jumpanalysis is analyzed.The simulation results show that the designed suspension system meets therequirements of the rules,but the wheel alignment parameters change greatly with the verticaldisplacement of the wheel,which is not conducive to the handling stability.In ADAMS/Insight,thedesign of experimen(DOE)analysis and optimization of the hard point coordinates of the model arecarried out based on the positioning parameters of the wheels.The comparison results show that the作者简介：安嘉晨（2 0 0 0 一），男，研究方向为汽车动力学，E-mail:。通信作者：臧利国（198 6 一），男，博士，教授，研究方向为车辆系统动力学与控制，E-mail:。基金项目：江苏省大学生实践创新训练计划项目（2 0 2 1112 7 6 0 2 7 Z）；南京工程学院科研基金（CKJA202205）。文献标识码：A文章编号：16 7 1-7 98 8(2 0 2 3)2 0-51-0 8Design and Research of FSAE Racing Suspension52amount of change in wheel alignment parameters when the wheel runout after optimization issignificantly reduced,and the performance of the car is well improved.Keyword:FSAE;Suspension system;Vehicle kinematics simulation悬架是大学生方程式赛车的重要组成部分之一，国内外各大高校与汽车生产厂商均给出了许多不同的悬架设计布置方案。石璇等 在设计悬架弹簧布置方式时，选取了双横臂搭配第三弹簧的设计结构，并借助第三弹簧完成解耦。刘凯华等2 采用了互联双横臂独立悬架的结构设计，使得车轮定位参数在一定范围内得到有效控制，从而保证车辆良好的操纵稳定性。王乐等3使用了多连杆式悬架，该悬架的空间自由度更为灵活，能为赛车提供更大的横向及侧向刚度，给予轮胎更好的抓地力。王德朝等4在CATIA建立了不等长双横臂前悬架模型，并把模型拆解为各个运动副做刚性连接，以简化建模过程。邓健贤等5利用ADAMS/Car建立悬架系统模型，先利用MATLAB制作一个设计界面，得出所有悬架系统参数后在ADAMS中构建最终模型。陈刚等6 在建立前悬架运动学模型时做了若干简化和假设，将悬架简化为一个多连杆机构，假设其零部件连接不考虑间隙，各运动副视为刚性连接，减震器阻尼特性为线性，车轮上下跳动时不考虑车身运动。吴振昕等7 结合数值计算方法和空间机构运动学建立了双横臂独立悬架运动学模型，并对悬架的K&C特性进行了仿真计算。张建等8 在分析悬架的运动响应时，设计了一种基于车辆加速度的悬架行程解耦算法以及一种通过悬架运动行程判别路面不平度的算法，并进行了道路试验验证该方案的可行性。SINDHWANI等9借助SolidWorks自主设计摇臂和横臂，以改善人体工程学设计，通过早期整车参数估算整车重量，并推算每个车轮所分得的质量进而推算整车的质心位置。BALENA等10 借助多系统耦合仿真将悬架、转向、制动等输入全部设置在同一物理环境下进行模拟，并将悬架与底盘视作刚体，以减少仿真时的工作量。本文以国内外高校大学生方程式赛车悬架系统的研究为支点，确定一套悬架系统的设计和选型方案。完成设计后将模型导入ADAMS中进行相关的运动学仿真及优化，以改善设计的合理性及实用性。1悬架结构的设计与选型悬架系统在设计初期需要考虑大量因素，包汽车实用技术括但不限于整车基本参数、各系统质心分布、悬架布置形式、悬架硬点参数、悬架刚度计算、传递比和载荷计算、悬架零部件设计等。1.1悬架的基本设计思路悬架的基本设计思路如图1所示。整车参数确定系统质心分布悬架选型悬架几何设计零部件设计运动学仿真零部件优化总结评价图1悬架设计技术路线图1.1.1 确定基本参数根据赛车已知的整车基本参数，在确定好发动机与传动系的空间位置关系后，确定轮距、轴距等基本参数。取设计轴距为17 0 0 mm；前轮距为12 7 0 mm；后轮距为12 40 mm。1.1.2确定各子系统质量分布本文借助虚拟样机规划整车各系统质量分布和前后轴的载荷配比，最终确定各系统的空间位置和质量分布。已知整车的设计满载质量约为330kg，同时根据计算得到质心距前轴距为933.53mm，质心距后轴距b为7 6 6.47 mm。简单计算可知前后轴载荷比分别为45%和55%，符合大学生方程式赛车的设计思路。2023年单体壳转向系解耦悬架车轮定位参数刚度计算阻尼计算杆系设计轮边零部件软件建模运动学仿真目标参数优化摇臂优化轮毂优化立柱优化(8)第2 0 期1.2悬架的选型根据大学生方程式汽车大赛（Formula Societyof Automotive Engineers,FSAE）比赛规则，赛车必须装有可自由工作、带有减震器的悬架；且悬架行程不得小于50 mm，当车手坐进赛车时能保证轮边有上下2 5mm的轮跳范围。综合考虑空间布置、成本计划等因素，双横臂式悬架是最理想的结构方案。由于赛车将采用半单体壳的设计方案，且综合考虑到悬架系统的重心及弹性元件外部的空间布置等因素，前悬架拟采用双横臂直推式悬架，借助防倾杆来增加侧倾刚度。后悬拟采用双横臂双弹簧解耦悬架，用一根斜置弹簧控制侧倾刚度，另一根横置弹簧控制俯仰刚度。1.3悬架几何参数设计1.3.1车轮定位参数车轮定位参数包括车轮外倾角、主销后倾角、主销内倾角及车轮前束角。初步拟定车轮定位参数如表1所示，通过设定基本车轮定位参数来建立悬架硬点坐标，后期根据车轮定位参数变化范围进行悬架硬点参数优化。表1初步拟定车轮定位角参数参数主销内倾角主销后倾角车轮外倾角车轮前束角取0.2 5m，代入参数计算得侧倾增益为-1.11（)/g;前悬4后悬1.951.3.2 刚度计算在计算悬架刚度之前需要先确定悬架的偏频，通常大学生方程式赛车的偏频在2.44Hz之间，并使前悬架略大于后悬架。根据公式计算并暂取前后悬架的偏频f=3.2Hz、f=2.8 H z；计算前后轴单侧载荷（去除车手质量），得m=58.5kg、mt=71.5kg；计算前后轴悬架适乘刚度，得K=23 649.15 N/m、K RR=2 2 130.0 2 N/m。车轮中心刚度指的是车轮中心相对车架垂直位移所受的垂向力，赛车拟采用13英寸轮胎，取径向刚度Kr=100 700N/m，根据式（1）、式（2）求解前后悬架车轮中心刚度KwF、K w R:KwF=K-KRFK.KRRKwR=TKt-KRR根据计算结果，选取前后车轮中心刚度KwF、KwR分别为30 9 0 9 N/m、2 8 36 3N/m。为研究整车在转弯时悬架的运动场景，需要安嘉晨，等：FSAE赛车悬架的设计与研究KoF2180RF12元KKPR2180RR式中，t为前轮轮距，取1.2 7 m；t 为后轮轮距，取1.2 4m。计算得出前后悬架侧倾角刚度KoF、K o R分别为435.1 Nm/()、38 0.58 Nm/()。考虑到赛车在场地内转弯工况的模拟，假设赛车以v=40 km/h的速度过弯，转弯半径R为9m，根据式（5）计算转向时横向加速度A,为1.4g。A,=v/(Rg)(5)但在很多情况下，赛车转弯工况贴近以1g的横向加速度过弯，赛车其以1g的横向加速度过弯时车架和车身侧倾转角的大小被称之为侧倾增益，根据式（6）计算侧倾增益。P=-mgxHA,Kor+KoR单位：（）式中，为车身侧倾角；H为质心到侧倾轴线距离，1-12.450KKRE53计算车身在经过单位转角扭转时总的弹性恢复力矩，根据式（3）、式（4）计算前后悬架侧倾角刚度KoFKoR:(3)(4)(6)-1一般大学生方程式赛车的侧倾增益约在1.0 01.8()/g，计算得出的1.11()/g在目标参数合理置信区间内，考虑到前悬架采用直推式悬架，会需要更多的侧倾角刚度，故拟添加防倾杆。由于转弯产生了横向加速度，从而导致侧向发生了力的偏移，称为载荷转移。根据式（7）、式（8）计算前后轴因横向加速度所引起的载荷转移。AWKopHWFA.WKRHazRRWRKor+KoR式中，ZRF、ZRR分别为前后悬架静态侧倾中心高度，分别取-6.410 m、-4.4x 10 3m；计算载荷转(1)移WF、W R分别为46 8 N、436 N。至此，根据载荷转移反算悬架适乘刚度和悬架偏频以证明拟定(2)偏频设计合理。通过悬架几何结构确定前后悬架的传递比MRi=1.2、M R2=1.3；计算前后悬架的弹簧刚度得Ks=44 508.96 N/m、K s R=47 933.47 N/m；计算前后悬架弹簧的实际行程得 ZsF=ZsR=0.02 m。2元KbZRE)KFKR(7)1154根据上述行程、刚度等计算挑选减震器弹簧，拟选择型号RCP2S-190的山地车减震器，可调磅数550 1bs，磅数调节范围方便后悬架解耦策略，需实车落地后调节预载荷来实现完全解耦。1.3.3 阻尼计算在选择好目标减震器弹簧后，根据其设计生产参数计算减震器的阻尼系数：S=2ymsw=4ymsf式中，为相对阻尼系数；ms为簧上质量，kg；f为悬架偏频，Hz。减震器的相对阻尼系数决定了其减震效果，一般弹簧在设计时会把压缩时的相对阻尼系数b取得更小，回弹时的相对阻尼系数取得更大，即减震器在压缩过程中的减震效果不如回弹时的减震效果明显；在取b、的时候，取平均相对阻尼系数a讨论取值范围。对于无内摩擦的弹性元件悬架，取a=0.250.35，而在有内摩擦的弹性元件悬架，a会取得略大；综合考虑后，确定压缩和回弹阻尼相对系数a=0.3、=0.4。根据式（10）一式（13）计算前后悬压缩、回弹阻尼系数SBF、O RF、O BR、O RR:0r=4Vgm.JfORF=4元V,m.f0BR=4元V,m.J,RR=4元V,m.J.将之前计算的前后轴单侧载荷mf、mr，前后悬架偏频f、f 代入方程计算可得：前悬压缩回弹阻尼系数分别为7 0 5Ns/m、9 41Ns/m；后悬压缩回弹阻尼系数分别为6 17 Ns/m、8 2 3Ns/m。1.4悬架零部件设计1.4.1 杆系设计1.横臂、推杆设计在设计悬架横臂时，需考虑杆系的尺寸选择以及连接方式：从不同工况下悬架横臂的应力和变形情况以及加工的难易程度来考虑杆件的尺寸设计；综合考虑后，拟采用外径16 mm壁厚1.5mm的钢管作横臂，配合螺纹焊接套、杆端轴承及吊耳焊接。并用12.9级M8塞打螺栓配合M6防松螺母对横臂与立柱、车架安装点进行装配与连接。因前悬架为直推式，当赛车静置于地面时，推杆多受到轴向力，故可适当减小推杆尺寸来减少悬架系统的整体质量；采用外径14mm壁厚1.5mm钢管配合螺纹焊接套进行焊接。汽车实用技术2.防倾杆、摇臂设计拟通过设计U形防倾杆增加前悬架侧倾角刚度。防倾杆通过摇臂与连接杆、扭臂、扭杆把左右悬架连接，将防倾杆布置在车内底板，借助预埋件进行吊装，合理布置空间。防倾杆的装配如图2 所示。(9)根据空间硬点图设计摇臂外形，借助塞打螺栓和防松螺母连接各点，并根据有限元软件仿真分析，在保证摇臂基本结构性能的情况下适当减料，以达到悬架系统轻量化的目的。3.杆系装配(10)杆系设计完成后在CATIA中进行整体装配，(11)并在运动模拟分析模块（DMU）验证悬架杆系运(12)动范围以及合理性。将减震器模型、摇臂、防倾杆、扭杆、推杆、横臂、吊耳等依次连接，最终(13)得到杆系装配如图3所示。图3前悬架杆系装配模型1.4.2轮边零部件设计轮边的正向设计思路为优先考虑制动系，并计算发动机至差速器、传动轴极限转速来对轴承尺寸选型；同时考虑转向系，确定立柱的大致轮廓和关键节点，在确定好上下横臂的安装点后设计其安装方式；通过查阅资料确定轴承类型及型号，方便后期端盖、保护罩、缸套的选型和设计；2023年图2 防倾杆装配模型第2 0 期然后是转向吊耳位置的设计，转向吊耳的高度应优先考虑转向梯形的布置；最后再考虑轮速传感器的安装位。1.立柱及轮毂设计立柱在轮边零部件中最为重要，其具备为转向拉杆提供连接点、为制动卡钳提供安装点、设置轴承安装座、保证悬架横臂硬点的安装、设置传感器安装点等多项功能。立柱属于簧下质量部分，在设计中对其进行空间拓扑优化分析，可去掉部分结构以适当减小其质量。轮毂通过螺栓和轮辋连接，其运动规律与轮辋相互绑定；同时还具有安装制动盘的功用。其中后轮毂需借助三球销和花键轴传递半轴的力和力矩，同时可设计钢套减少三球销与轮毂的直接磨损。2.轮边装配完成设计后，用CATIA的DMU模块检查无机械干涉即可。将轮边所有零部件进行装配，包括但不限于：立柱、轮毂、端盖、刹车盘、轮辋、轮胎、轴承等，具体装配效果如图4所示。安嘉晨，等：FSAE赛车悬架的设计与研究的模型如图5所示。2.2悬架系统运动学仿真FSAE规则规定：在有车手乘坐的情况下，赛车轮胎的跳动行程至少为50 mml。所以在建模完成后切换至ADAMS/VieW模块，预设轮胎空载时的车轮自由半径为2 50 mm；车轮胎端面宽度为200mm；扁平率为45%；垂向刚度为2 0 0 kN/m；径向刚度为19 0 kN/m。分别对悬架进行静载、平行轮跳、单侧轮跳分析。H55Car模块中建立悬架模型并添加运动副，最终建立(a)静载仿真图4轮边总装模型(b)平行轮跳仿真2悬悬架系统的动力学建模与仿真2.1悬架系统运动学建模在确定好悬架系统的总体设计后，在ADAMS/1677.5(c)单侧轮跳仿真图6 悬架运动学仿真示意图进入静态载荷仿真模块，由静载仿真结果可知当空载时，轮边立柱受到的力约为57 0 N，方向垂直地面向上，符合预期，受力分析如图6(a)所示。图5解耦悬架模型为模拟赛车在起步或制动时的后仰和点头情56况，对赛车进行平行轮跳仿真；根据规则要求设计轮跳行程为2 5mm的仿真条件，观察减震器和摇臂的运动状况。受力以及力矩情况如图6(b)所示。当赛车以1.4g的加速度过弯时，需要验证悬架系统的可靠性，设置单侧轮跳工况分析来模拟过弯时的悬架运动过程。观察悬架运动规律并分析解耦悬架中弹簧的实际运动方式。单侧轮跳分析时，单侧轮处于+2 5mm位置时立柱的受力及力矩情况如图6(c)所示，结果证明悬架系统可靠性良好。2.3轮胎定位参数分析将硬点参数进行平行轮跳仿真试验的数据进行后处理并利用ADAMS/Insight模块作图分析，横坐标为轮跳行程，纵坐标为对应车轮定位参数角变化；因车轮上下跳动工况一致，为优化仿真设置并减少计算步骤，主要以正轮跳部分作为仿真行程设置，具体分析如图7 所示。如图7(a)所示，悬架设计硬点参数在平行轮跳时引起的车轮外倾角变化为0-1。处于静平衡位置时外倾角为0，变化率为0.9()/2 5mm。如图7(b）所示，悬架设计硬点参数在平行轮跳时引起的主销后倾角变化为2.4-2.95。处于静平衡位置时主销后倾角为2.45，变化率为0.5()/25 mm。0.10.0-0.1-0.23op/-0.3-0.4$-0.5-0.60.7-0.8-0.9-5.03.02.9Bop/2.82.72.62.52.4-5.0汽车实用技术2.82.72.62.52.42.32.22.12.01.91.8-5.00.200.15iop/串0.100.050.0-0.05-5.0图7平行轮跳下车轮定位参数变化示意图如图7(c)所示，悬架设计硬点参数在平行轮跳时引起的主销内倾角变化为1.9-2.7。处于静平衡位置时主销内倾角为1.945，变化率为0.83()/25 mm。如图7(d)所示，悬架设计硬点参数在平行轮跳时引起的前轮前束角变化为0 0.17。处于静平衡位置时车轮前束角为0，变化率为0.17(）25 mm。3悬架系统运动学优化在进行轮胎定位参数的分析后，需进一步利用ADAMS/Insight模块进行试验设计（Design OfExperimen,DOE）优化分析，通过设置优化参数、0.05.0轮跳行程/mm(a)车轮外倾角0.05.0轮跳行程/mm(b)主销后倾角2023年0.05.0轮跳行程/mm(c)主销内倾角0.05.0轮跳行程/mm(d)车轮前束角10.015.010.015.010.010.020.025.020.025.015.015.0目标参数以及优化方法和权重来实现对悬架的目标参数优化。3.1目标参数优化设置通过修改悬架硬点坐标来进行车轮定位参数的优化，考虑到悬架系统与车架的连接点已经确定，车轮轮距、轴距等参数基本无法修改。因此，在设置的过程中，不会对硬点的Y坐标进行修改。因为左右悬架对称设计，故选取单边作为研究对象，拟挑选上下横臂前点和后点的X坐标及Z坐标、导向杆内点X坐标及Z坐标作为参数优化的变量，四个车轮定位参数作为目标优化值。设置参数后即可进入Insight模块进行参数优20.020.025.025.0第2 0 期化，设定变量变化范围在3mm之内，计算方式采用线性计算，初步拟定4个目标参数优化权重均为1.0。优化前后四个定位角参数如表2 所示。表2 优化前后车轮定位角参数对比单位：()主销内主销后参数倾角优化前1.95优化后2.713.2优化后硬点参数DOE分析对比将优化后硬点参数进行单侧轮跳仿真试验的数据进行后处理，直接在原车轮定位参数的数据图上添加优化后的变化曲线，如图8 所示。如图8(a)所示，优化后的悬架设计硬点参数在单侧轮跳时引起的车轮外倾角变化为-0.7 0.7 5。处于静平衡位置时，车轮外倾角为0，变化率为1.45()/50mm。经DOE分析优化后，整体变化范围减小，实现参数优化目标。如图8(b)所示，优化后的悬架设计硬点参数在单侧轮跳时引起的主销后倾角变化为2.2 2.7。处于静平衡位置时，主销后倾角为2.45，变化率为0.5()/50 mm。经DOE分析优化后，整体变化范围减小，实现参数优化目标。如图8(c)所示，优化后的悬架设计硬点参数在单侧轮跳时引起的主销内倾角变化为1.2 52.75。处于静平衡位置时，主销内倾角为1.95，变化率为1.5()/50 mm。经DOE分析优化后，整体变化范围减小，实现参数优化目标。如图8(d)所示，优化后的悬架设计硬点参数在单侧轮跳时引起的车轮前束角变化为-0.150.15。处于静平衡位置时，车轮前束角为0，变化率0.3()/50 mm。因后悬架无转向机构，仅设置导向杆用于约束车轮随动转向，图中纵坐标表示的束角变化范围变化值小于0.5，属于合理变化范围。1.00.530p/0.0-0.5-1.0-25.0-20.0-15.0-10.0-5.0安嘉晨，等：FSAE赛车悬架的设计与研究车轮车轮倾角外倾角2.4502.720.770.05.010.015.020.0 25.0轮跳行程/mm(a)车轮外倾角573.02.0前束角00.051.525.020.015.0-10.05.03.02.5Bop/2.01.51.025.020.015.010.05.00.05.00.20.1Sop/0.0-0.1-0.2-25.0-20.0-15.0-10.05.00.05.0图8 优化后单侧轮跳车轮定位参数变化示意图4总结1）基于FSAE比赛规则并采用正向设计思路提出了一套大学生方程式赛车悬架系统的设计方案，根据整车的布置特点确定了悬架具体的几何参数，进行了悬架设计参数的计算及构件选型，并在CATIA中进行了建模。2)在ADAMS/Car和ADAMS/View中建模并分析了悬架的运动学性能，仿真的结果表明，所设计的悬架系统满足比赛规则要求。3）在ADAMS/Insight中基于车轮的定位参数对悬架的硬点坐标进行了DOE分析优化。优化前后仿真结果的对比表明，优化后车轮跳动时赛车的性能得到了很好的改善。0.05.010.015.020.025.0轮跳行程/mm(b)主销后倾角10.015.020.025.0轮跳行程/mm(c)主销内倾角10.015.0020.025.0轮跳行程/mm(d)车轮前束角584）在后期悬架的实际制造及调试过程中已对模型构建方法及仿真分析数据进行了相关验证，证明了本文所述内容的合理性。参考文献1石璇,叶胤欣,彭育辉.基于第三弹簧结构的前悬架俯仰运动分析J.时代汽车,2 0 2 0(8):7 9-8 2.2 刘凯华，王庆,李强.装有互联双横臂独立悬架的方程式赛车仿真分析J.浙江科技学院学报,2 0 19,31(3):239-246,252.3 王乐,朱建军,田宇.大学生方程式赛车多连杆悬架设计及优化J.机械设计与制造,2 0 18(2):9-12.4王德朝，丁华锋,孟凡净，等.FSEC赛车前悬架系统仿真与优化J.重庆理工大学学报(自然科学),2 0 18,32(4):21-26.5 邓健贤，查云飞.大学生方程式赛车悬架系统设计CJ/2021中国汽车工程学会年会论文集.北京：机械工业出版社,2 0 2 1:6.汽车实用技术6陈刚,曹永青.基于ADAMS的大学生方程式赛车前悬架设计与优化J.三明学院学报,2 0 17,34(6):33-38.7吴振昕,付振,程超,等.基于运动学特性的双横臂悬架硬点位置优化J.汽车技术,2 0 15(2):39-43.8张建,刘秋铮,王康,等.基于悬架行程解耦的 SUV 路面不平度识别算法研究J.汽车工程,2 0 16,38(10):1232-1236,1199.9SINDHWANI R,BHATNAGAR A,SONI A,et al.Design and Optimization of Suspension for FormulaSociety of Automotive Engineers(FSAE)VehicleJ.Materials Today,Proceedings,2021,38:229-233.10BALENA M,MANTRIOTA G,REINA G.DynamicHandling Characterizationand Set-up Optimization fora Formula SAE Race Car via Multi-body SimulationJJ.Machines,2021,9(6):126.11 李理光.中国大学生方程式汽车大赛规则Z.北京：中国汽车工程学会,2 0 2 2.2023年