RESEARCH ON TRAILING ARM SUSPENSION SYSTEM OF FULL X-BY-WIRE CONTROL CHASSIS BASED ON DATA DRIVE
-
摘要: 汽车在越野类极限路况下行驶, 对车身高度有一定范围的调节需求, 传统悬架方案与全线控底盘进行技术融合时, 存在机构运动干涉、底盘升降过程中车轮外倾程度过大、车轮发生侧向位移等现象, 易导致轮胎过度磨损, 致使行驶失稳. 将车身高度变化对轮胎侧向参数的影响转化为车轮纵向滚动, 进而实现稳定的大行程车身高度调节, 是解决上述问题的关键. 本研究建立整车七自由度动力学模型, 对悬架系统导向机构展开力学分析, 集两者作为系统研究的输入信息; 通过正弦波激振台对弹性元件、减振器进行相关特性参数获取, 基于数据驱动开展一体化电动轮的运动学仿真测试, 包括对悬架系统关键铰接位置进行力学性能分析、对电动轮整体结构进行运动学特性研究, 以此定义系统关键性能指标, 结合理论研究与仿真测试, 确定双纵臂式主动悬架系统方案. 仿真结果与实车验证综合表明, 搭载本研究系统方案的全线控平台, 进行大行程高度调节过程中, 车轮外倾问题得到有效解决, 一体化电动轮具备良好的独立运行能力, 本研究对提高车辆在极限路况下的通过性具有重要意义.Abstract: The car is driven under extreme road conditions such as off-road, and there is a certain range of adjustment requirements for the body height. When the traditional suspension scheme is technically integrated with the full-wire chassis, there are mechanism motion interference, excessive wheel camber during the chassis lifting process, and wheel generation. Lateral displacement and other phenomena can easily lead to excessive tire wear, resulting in unstable driving. The key to solving the above problems is to convert the influence of vehicle height changes on the lateral parameters of tires into longitudinal rolling of the wheels, thereby achieving stable and long-travel vehicle height adjustment. . In this study, a seven-degree-of-freedom dynamic model of the whole vehicle is established, and the mechanical analysis of the suspension system guiding mechanism is carried out, and the two are collected as the input information of the system research; the relevant characteristic parameters of the elastic elements and shock absorbers are carried out through the sine wave vibration table. Acquire, carry out the kinematics simulation test of the integrated electric wheel based on the data drive, including the mechanical performance analysis of the key hinge position of the suspension system and the kinematic characteristics study of the overall structure of the electric wheel, so as to define the key performance indicators of the system, and combine the theory Research and simulation tests are carried out to determine the scheme of the double trailing arm active suspension system. The simulation results and the real vehicle verification show that the wheel camber problem can be effectively solved in the process of adjusting the height of the large stroke with the full-wire control platform equipped with the system scheme of this study. , the integrated electric wheel has good independent operation ability, and this research is of great significance to improve the vehicle's passability under extreme road conditions.
-
引 言
智能网联汽车逐渐成为汽车工业的研究热点, 底盘线控执行控制成为其发展的关键技术[1]. 底盘动力学域控技术通过整合全车传感器信息, 充分发挥出底盘性能[2]. 集成传统车辆驱动-制动-转向和主动悬架系统的全矢量底盘技术具有重要的发展潜力, 车辆是一类典型的多系统耦合复杂非线性动力学系统[3], 凭借其拥有更多的可控自由度从而更加有利于实现车辆的智能化控制[4]. 对于多自由度模型, 当自由度取得越多时, 理论上越逼近真实系统[5]. 从目前发展趋势以及各种驱动技术的特点来看, 电动汽车是实现交通领域节能减排和污染防治的重要举措[6], 全线控电动汽车成为解决复杂交通环境问题的重要工具, 代表了未来电动汽车的发展趋势[7-8].
目前, 车辆行驶平顺性正从随机路面、匀速直线行驶工况延展到凸台、凹坑、变速等越野极限工况[9], 这对悬架系统性能提出了重要挑战, 通过吸收和衰减来自路面不平而产生的振动和车轮波动, 悬架系统可以让车辆平顺地行驶[10-11], 其对车辆横向动力学密切相关的振动进行了阻尼作用[12]. 国内从20世纪末开始对悬架运动学、弹性运动学, 即 K&C特性进行系统研究[13].
李韶华等[14]搭建刚柔耦合的重型汽车整车模型, 以车身垂向加速度和各轮垂向轮胎力为评价指标, 提升了减速带连续通过的行驶平顺性. 李小彭等[15]基于振动功率流理论对减振元件衬套进行优化选型, 有效改善车辆悬架系统的综合减振性能. 王震等[16]建立完整的机械-气体耦合多自由度动力学联合仿真模型, 为侧翻事故发生率极高的全地形车设计了一种采用囊式空气弹簧的交叉型双气室空气互联悬架. 纪秀业等[17]对悬架系统刚度与阻尼的匹配计算及验证, 为后续悬架系统设计及操纵舒适性提升形成理论基础. 张智等[18]对悬架的几何参数进行分析, 以变化量处于更为合理范围内的车轮定位参数为目标进行优化设计.
王军年等[19]对传统汽车底盘进行电动轮改进, 利用全局非归一化的多目标遗传优化算法对优化解集实现验证. 李凡杰等[20]基于机电相似理论, 对汽车悬架结构及参数进行优化设计, 进一步提升了汽车悬架的减振性能. 袁宝峰等[21]在摇臂式悬架的基础上, 基于火星车地面力学等特点对主动悬架进行优化设计. 张农等[22]对液压互联悬架的主要参数进行了全局灵敏度分析, 通过NSGA-Ⅱ算法进行多目标优化, 显著提升了综合性能.
陈龙等[23]建立新型悬架系统前轮定位参数的理论计算模型和多体动力学模型, 解决传统悬架不适合轮毂驱动的难题. 陆建辉等[24]将遗传优化算法与多体运动学分析方法相结合, 对麦弗逊前悬架的设计参数进行优化. 李芳[25]设计出改进型遗传算法对悬架参数进行了优化, 提高了被动状态下的汽车性能.
本文从电动汽车线控技术的发展趋势出发, 基于作者团队提出多自由度集成的全矢量概念, 针对传统悬架方案存在的机械结构配合限位等相关问题, 开展承载全矢量线控平台的双纵臂独立悬架结构方案设计, 基于数据驱动背景, 通过悬架子系统参数获取以及参数化设计变量的动态调整, 进而确定最优方案.
1. 一体化电动轮概述
传统汽车可控输入量较少, 故对其动力学控制难度较大, 路面的激励作用在单个车轮上可体现成3个相对独立的力, 从理论角度出发, 传统四轮车辆系统的独立输入集合可扩展至12个, 基于此, 作者团队提出如下定义: 如果单个车轮的横向、纵向和垂向作用力都可进行独立控制, 那么就称为全矢量控制(full vector control, FVC)汽车[1]. 作为自动驾驶的关键执行系统, 线控技术可以实现代替驾驶员的物理操纵而对车辆发出行驶指令. 传统的线控底盘未对车辆关键工作系统进行耦合协作, 在此, 作者团队设计驱动-制动-转向-悬架一体化的多功能电动轮系统, 如图1, 面向线控底盘技术, 进一步提高整车灵活性和可控性.
一体化电动轮系统作为全矢量线控底盘的重要搭建板块, 其对多子系统进行高度集成于一体, 突破传统车辆单个系统间的耦合难题, 进而实现全矢量线控技术给智能移动平台带来的技术空缺. 对于多功能电动轮, 其驱动系统采用高功率密度轮毂电机方案; 制动系统采用再生制动与线控液压制动复合制动方式; 转向系统采用转向电机经过减速器直接带动转向臂旋转方案, 可实现车轮360°独立转向功能.
主动悬架对平顺性和操纵稳定性之间的冲突实现较好的平衡, 显著提升乘坐舒适性和车辆行驶安全性[26-27], 因此, 本系统选取主动悬架为设计对象. 悬架系统装配空气悬架可实现主动车身高度调节, 装配磁流变阻尼器可实现阻尼线控调节, 进而实现车轮垂向力的主动可控, 悬架系统设计是本文重点开展研究的内容.
弹性元件采用空气弹簧形式, 空气弹簧内部充满压缩空气, 以压缩空气作为工作介质, 满载工作压力通常为0.8 MPa左右, 当空气弹簧压缩时, 气囊沿活塞座轮廓面向下挠曲变形[28].
磁流变阻尼器是一种阻尼可控的执行器, 阻尼力可控范围大, 响应时间可精确至毫秒级[29], 因此, 采用磁流变阻尼器作为悬架系统的减振器模块, 凭借其特性, 可实现对外部激励的有效吸能与缓冲减振.
簧下空间直接影响悬架系统的结构尺寸, 一体化电动轮通过高密度轮毂电机实现驱制动, 轮毂电机具有快速响应、转矩精确可控等特点, 通过驱动转矩协调分配[30], 簧下质量大导致轮胎动载荷增大, 电机电磁力直接影响车辆垂向振动[31-32], 使车轮动载荷变大, 直接影响平顺性评价指标三个指标中的车轮相对动载, 对电动汽车平顺性造成了一定的负面效应[33]. 为了解决上述问题, 考虑簧下所占空间、搭载质量受一体化电动轮中轮毂电机的影响而发生量化, 将转向中心线与转向电机轴线重合设计, 即主销倾角设置为0, 在不发生运动干涉的同时实现全方位转向.
2. 悬架系统动力学与力学分析
2.1 车辆系统建模
本文搭建如图2 的车辆系统动力学模型, 用以研究车辆垂向动力学的相关问题. 采用7 DOF 动力学模型, 包括车身垂向、俯仰和侧倾三个自由度; 四个车轮独立垂向四个自由度, 车辆坐标系为 OXYZ. 在下文中, FL, FR, RL 和 RR 或者数字 1 ~ 4 分别表示左前轮、右前轮、左后轮和右后轮. B 为轮距, a和b分别为质心到前后轴距离.
车辆与车轮垂向动力学模型
$$ {m_{\text{s}}}\ddot z = {F_1} + {F_2} + {F_3} + {F_4} - {m_{\text{s}}}g $$ (1) $$ {m_{{\text{u}}i}}{\ddot z_{{\text{u}}i}} = {k_{{\text{t}}i}}\left( {{z_{{\text{u}}i}} - {z_{{\text{r}}i}}} \right) - {F_i} $$ (2) 式中,
$ {{z}} $ 为车辆质心垂向位移;$ {m_{\text{s}}} $ 与$ {m_{{\text{u}}i}} $ 为簧上质量与车轮质量;$ {z_{{\text{u}}i}} $ 与$ {z_{{\text{r}}i}} $ 为车轮垂向位移与路面输入;$ {F_i} $ 为悬架力.各悬架力可表示为
$$ {F_i} = {k_{{\text{s}}i}}\left( {{z_{{\text{s}}i}} - {z_{{\text{u}}i}}} \right) + {c_i}\left( {{{\dot z}_{{\text{s}}i}} - {{\dot z}_{{\text{u}}i}}} \right) + {F_{{\text{a}}i}} $$ (3) 式中,
$ {k_{{\text{s}}i}} $ 为悬架刚度系数,$ {c_i} $ 为悬架阻尼系数,$ {F_{{\text{a}}i}} $ 为空气悬架主动作动力,$ {z_{{\text{s}}i}} $ 悬架与车身安装点得垂向位移.车身俯仰运动和侧倾运动微分方程为
$$ {I_y}\ddot \theta = \left( {{F_3} + {F_4}} \right)gb - \left( {{F_1} + {F_2}} \right)ga + {m_{\text{s}}}{a_x}{h_{\text{p}}} $$ (4) $$ {I_x}\ddot \phi = \left( {{F_1} + {F_3}} \right)\frac{B}{2} - \left( {{F_2} + {F_3}} \right)\frac{B}{2} + {m_{\text{s}}}{a_y}{h_{\text{r}}} $$ (5) 式中,
$ {I_x} $ 和$ {I_y} $ 为车身绕x和y轴的转动惯量,$ \ddot \theta $ 和$ \ddot \phi $ 为俯仰、侧倾角加速度,$ {a_x} $ 和$ {a_y} $ 为纵、侧向加速度,$ {h_{\text{p}}} $ 和$ {h_{\text{r}}} $ 为车身质心到俯仰、侧倾中心的垂向距离.当车身俯仰角较小时, 悬架与车身连接点位的垂向位移与车身垂向位移的关系如下
$$ {z_{{\text{s}}1}} = z - a\theta + \frac{B}{2}\phi $$ (6) $$ {z_{{\text{s}}2}} = z - a\theta - \frac{B}{2}\phi $$ (7) $$ {z_{{\text{s}}3}} = z + a\theta + \frac{B}{2}\phi $$ (8) $$ {z_{{\text{s}}4}} = z + a\theta - \frac{B}{2}\phi $$ (9) 式中,
$ a $ 为车辆质心到前轴距离,$ B $ 为轮距,$ \theta $ 为车身俯仰角,$ \phi $ 为车身侧倾角.2.2 悬架系统力学分析
根据总体设计方案, 考虑结构紧凑、整车布置空间等建立出面向一体化电动轮系统的悬架结构如图3所示. 其中转向臂上端与整车转向机构相连, 受转向电机作用力带动一体化电动轮做周转运动; 双纵臂一端与转向臂铰接, 另一端与转向节铰接, 悬架工作过程中, 分别绕铰接点做往复摆动, 实现车轮的垂向行程; 减振器上端与转向臂铰接相连, 下端固定在下纵臂, 装配轮毂电机、制动卡钳以实现驱、制动功能. 其中, 分别对6处关键铰接位置进行命名, 上下纵臂与转向节铰接处分别为A和B, 与转向臂铰接处分别为D和C; 减震器上端铰接点为
$O_1$ , 下端为$O_2$ .减振器作为悬架系统中重要工作部件, 通过其往复运动实现对路面不平带来的冲击的缓冲, 其安装位置应作首要考虑, 通过对其进行静力学分析, 可以判断出关键连接点的位置. 根据减振器的结构尺寸, 将减振器上端与转向臂铰接, 下端与下纵臂进行铰接. 保持上下纵臂静态位置时水平, 则上纵臂在车轮纵向平面内并不受力, 因此在对结构进行平面分析时, 可将上纵臂略去. 不考虑部件自身质量, 当空气弹簧气囊压力一定时, 在静平衡状态下可将减振器简化为受力连杆, 结构简化如图4所示.
其中B为转向节与下纵臂铰接点,
$ O_2$ 为减振器下铰接点, C为下纵臂与转向臂的铰接点,$ O_1 $ 为减振器上铰接点.$ L_1$ 为B和$ O_2$ 两铰接点的距离,$ L_2$ 为B和C两铰接点的距离,$ L_3 $ 为减振器上下两铰接点的垂向距离,$ L_4 $ 为减振器上铰接点到转向臂右侧距离, mg为单个电动轮载重, β为减振器倾斜角度, 其中$$F L_{3}+F_{O 1} L_{4} \cos \beta +F_{O 1} L_{3} \sin \beta -{mg} {L}_{2}=0$$ (10) $$F_{O 1}=\frac{{mg} {L}_{2}-F L_{3}}{L_{4} \cos \beta + L_{3} \sin \beta } $$ (11) $$ F_{O 1}=\frac{{mg} {L}_{2}-F L_{3}}{\left(L_{4}+ L_{3} \tan \beta \right) \cos \beta} $$ (12) $$ F_{O 1}=\frac{{mg} {L}_{2}-F L_{3}}{\left(L_{4}+\dfrac{L_{2}-L_{1}-L_{4}}{L_{3}} L_{3}\right) \cos \beta} $$ (13) 对于整个系统
$$ \sum M_{B}=F L_{3}=0 $$ (14) 故
$$ \left. \begin{array}{l} F=0\\ F_{O 1}=\dfrac{{L}_{2}}{\left(L_{2}-L_{1}\right)} \dfrac{{mg}}{\cos \beta} \end{array} \right\} $$ (15) 同理
$$ F_{O 2}=\frac{{L}_{2}}{\left(L_{2}-L_{1}\right)} \frac{{mg}}{\cos \beta} $$ (16) 根据导向机构与减振器关键铰接位置的受力情况分析可以发现, 影响减振器受力的主要因素是上下铰接位置, 其之间的差异也就是安装角度. 对减振器受力随其下端铰接位置的影响, 可知下纵臂与转向臂铰接位置到其与转向节铰接位置间距在下纵臂总长中的占比越小, 减振器受力越小, 但结合本系统实现垂向行程可控范围最大化的设计目标, 同时考虑过长的下纵臂导致转弯半径增加, 降低灵活性.
3. 悬架系统参数获取
3.1 磁流变阻尼器参数获取
磁流变阻尼器特性测试主要包括阻尼力特性测试及阻尼特性测试两方面, 分别为: 阻尼力Fd随位移X的变化关系和随速度v的变化关系.
假设振动试验台以100次/分的振动频率给所测减振器施加正弦激励, 测试行程X为100 mm (即振幅 ± 50 mm), 则阻尼器测试行程位移为
$$ X = {X_{\max}} \sin(wt) $$ (17) 活塞与缸筒的相对振动速度为
$$ \begin{split} v =& \frac{{{\rm{d}}X}}{{{\rm{d}}t}} = w {X_{\max}} \cos(wt) =\\ &2{\text{π}} f {X_{\max}} \cos(2{\text{π}} ft) \end{split} $$ (18) $$ {v_{\max}} = 2{\text{π}} f {X_{\max}}\approx 0.52\;{\text{m}}/{\rm{s}} $$ (19) 在某一时刻, 阻尼系数乘以该时刻的活塞与缸筒的相对运动速度可以得到该时刻的阻尼力大小
$$ \begin{split} \sin{^2}(2{\text{π}} ft) +& \cos{^2}(2{\text{π}} ft) =\\ &{\left( {\frac{X}{{{X_{\max}}}}} \right)^2} + {\left( {\frac{{{F_d}}}{{2{\text{π}} fc{X_{\max}}}}} \right)^2} = 1 \end{split} $$ (20) $$ {F_d} = cv = 2{\text{π}} fc{X_{\max}}\cos(2{\text{π}} ft)\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\; $$ (21) 因此, 理想的阻尼器阻力特性曲线近似为关于行程位移X和阻尼力Fd的椭圆方程, 该椭圆方程描述了给定位移的阻尼力大小.
对于磁流变阻尼器而言, 改变电流可以获得阻尼力变化, 因此阻尼系数并非恒定常数, 但在一个循环周期内, 其消耗能量的大小和被动阻尼器相等, 磁流变阻尼器的等效阻尼系数为
$$ {W_d} = \oint {{F_d}{\rm{d}}X} = \int_0^{2{\text{π}} /w} {{F_d}} \frac{{{\rm{d}}X}}{{{\rm{d}}t}}{\rm{d}}t $$ (22) 同时, Wd也可写作如下形式
$$ \qquad\begin{split} {W_d} =& \oint {{c_{eq}}\frac{{{\rm{d}}X}}{{{\rm{d}}t}}{\rm{d}}X} = {c_{eq}}\int_0^{2{\text{π}} /w} {{{\left(\frac{{{\rm{d}}X}}{{{\rm{d}}t}}\right)}^2}} {\rm{d}}t{\text{ = }} \\ & {c_{eq}}X_{\max}^2{w^2}\int_0^{2{\text{π}} /w} {\cos{{(wt)}^2}} {\rm{d}}t{\text{ = }} \\ & {c_{eq}}X_{\max}^2{w^2} \frac{{\text{π}} }{w} = {c_{eq}}{\text{π}} wX_{\max}^2 \end{split} $$ (23) 故磁流变阻尼器等效阻尼系数为
$$ {c_{eq}} = \frac{{{W_d}}}{{2{{\text{π}} ^2}fX_{\max}^2}} $$ (24) 其中, Wd是一个周期内Fd-X曲线所围成的面积, 表示阻尼器运动一个周期所消耗的能量, f为振动频率, Xmax为激励振幅. 由于磁流变阻尼器是通过改变输入控制电流的大小来改变磁场, 从而改变磁流变液的黏度, 使阻尼系数获得大范围改变, 从而提高磁流变阻尼器的耗能能力, 可以获得磁流变阻尼器的阻尼力大小为
$$ \begin{split} {F_d} = {c_{eq}}v =& \frac{{{W_d}}}{{2{{\text{π}} ^2}fX_{\max}^2}}2{\text{π}} f{X_{\max}}\cos(2{\text{π}} ft){\text{ = }} \\ & \frac{{{W_d}}}{{{\text{π}} {X_{\max}}}}\cos(2{\text{π}} ft) \end{split} $$ (25) 磁流变阻尼器的阻尼力大小和耗散能量Wd、振动频率f、振幅Xmax三个因素相关, 由于耗散能量和控制电流输入大小相关, 因此影响磁流变阻尼器阻尼力的直接因素主要包括控制电流i、振动频率f和振幅Xmax.
簧上质量与簧下质量相互远离时, 需要较大的阻尼力消除振荡, 在两者相互靠近时, 需要较小的阻尼力减小对车架的冲击. 在一个循环周期内, 压缩阻尼力小于复原阻尼力, 通过示功机可同时测出压缩阻尼力Fy和复原阻尼力Ff (图5).
取压缩阻尼力为
$$ {F_y} = \frac{{2{F_d}}}{{1 + \beta }} $$ (26) 复原阻尼力为
$$ {F_f} = \frac{{2\beta {F_d}}}{{1 + \beta }} $$ (27) 其中β为伸缩系数, 常取为10~15.
由于阻尼器内含有少量气体, 需要测量出气体反拔力并在测试结果中去除, 得到真实的油液阻尼力. 在0~2.1 A内以0.3 A间隔改变施加于阻尼器的控制电流, 固定振动频率为1 Hz, 固定振幅为50 mm, 获得的阻力特性曲线和阻尼特性曲线如图6 ~图7所示.
随电流以间隔0.3 A增加, 磁流变效应显著提升, 屈服应力不断增加, 示功曲线由内到外所围面积逐渐增大, 表明阻尼器耗能能力增强.
从阻尼特性曲线看出, 随电流增加, 曲线斜率不断增加, 表明电流的增加使得阻尼器阻尼不断变大. 当电流增大到2.1 A时, 阻力相较于1.8 A时变化已经不明显, 表明磁流变阻尼器阻尼力变化已达到饱和, 饱和电流为1.8 A.
因此, 在对磁流变阻尼器进行控制时, 可将电流调节范围限制在1.8 A以下. 同时可以看出在控制电流为0时, 减振器作为被动阻尼器的阻尼系数计算如下
$$ {C_d} = \frac{{\Delta {F_d}}}{{\Delta v}} = 633.8 \; \rm N \cdot s/m $$ (28) 由于气囊在筒壁上卷曲时具有一定的摩擦阻力和材料卷曲阻尼, 根据经验适当增大阻尼系数, 获得最终的被动阻尼器阻尼系数
$$ {C_s} = 635\;{\rm{N}} \cdot {\text{s}}/{\rm{m}} $$ 在0.5~2.5 Hz内以0.5 Hz间隔改变施加于阻尼器的激励频率, 固定控制电流为1 A, 固定振幅为50 mm, 获得的阻力特性曲线和阻尼特性曲线如图8~图9所示.
通过改变振动频率, 示功曲线所围面积增加, 最大阻尼力也逐渐增大, 表明磁流变阻尼器对外界振动频率较为敏感, 而外界振动来源于路面起伏, 使得悬架动挠度发生改变, 造成簧上与簧下质量相对运动速度变化, 相对运动速度变化越快, 则需要更多的能量来衰减振动, 从而使得最大阻尼力不断增大.
在10~50 mm内以20 mm间隔改变施加于阻尼器的振幅, 固定控制电流为0.3 A, 固定振动频率为1 Hz, 获得的阻力特性曲线和阻尼特性曲线如图10和图11所示.
综上所述, 通过三组实验测试结果分析了控制电流i, 振动频率f和振动幅度Xmax对阻尼力大小的影响, 由于控制电流i, 振动频率f的变化对阻尼力的影响远大于振动幅度Xmax的变化, 因此, 本文仅将控制电流i和振动频率f作为磁流变阻尼器模型的两输入, 对这两种因素进行适当控制, 最终获得期望的阻尼力, 故将磁流变阻尼器视作双输入单输出非线性模型, 其中建立正向模型表达式如下
$$ {F_d} = \sum\limits_{i = 0}^n {({a_i} + {b_i}I)} {v^i} $$ (29) 其中, Fd为磁流变阻尼力大小, I为电流幅值, v为阻尼器活塞速度, ai和bi为模型系数, n为模型阶次,为同时保证拟合精度并降低拟合难度, 取n=7. 建立逆向模型表达式为
$$ {F_d} - \sum\limits_{i = 0}^n {{a_i}} {v^i} = \sum\limits_{i = 0}^n {{b_i}I} {v^i} $$ (30) $$ I{\text{ = }}\frac{{{F_d} - \displaystyle\sum\limits_{i = 0}^n {{a_i}} {v^i}}}{{\displaystyle\sum\limits_{i = 0}^n {{b_i}} {v^i}}} $$ (31) 3.2 空气弹簧参数获取
设置正弦波激振台振动频率为0.02 Hz, 同时测试过程中保证气路的密闭性, 充入最大气压设置为0.6 MPa且采用最大伸张状态下单向激励的方式进行实验,弹簧初始高度为140 mm, 记录不同充气压力下空气弹簧作用力和囊内气压随弹簧高度变化间的关系, 计算出此时空气弹簧有效面积如表1.
表 1 空气弹簧测试结果Table 1. Air spring test resultsInitial
pressure/
MPaDownward pressure
displacement of the
exciter head/mmAir spring
force/
NIntracapsular
air pressure/
MPaEffective
area/
mm20.3 0 859 0.30 2863.3 10 1064 0.31 3432.3 20 1136 0.32 3550 30 1196 0.33 3624.2 40 1263 0.34 3714.7 50 1338 0.36 3716.7 60 1412 0.38 3715.8 70 1496 0.40 3740 80 1624 0.42 3866.7 90 1726 0.44 3922.7 100 1817 0.465 3907.5 0.6 0 1997 0.60 3328.3 10 2347 0.62 3785.5 20 2476 0.64 3868.8 30 2594 0.665 3900.8 40 2780 0.70 3971.4 50 2852 0.73 3906.8 60 2966 0.75 3954.7 70 3083 0.77 4003.9 80 3286 0.82 4007.3 90 3463 0.84 4122.6 100 3672 0.90 4080.0 通过三次样条曲线拟合获得空气弹簧作用力和囊内气压随空气弹簧高度变化的关系, 如图12和图13所示.
三次样条拟合模型表达式为
$$ P(h) = A + Bh(t) + C{h^2}(t) + D{h^3}(t) $$ (32) 根据数据拟合得到初始压力为0.3和0.6 MPa气囊内压力随高度的变化关系式如下
$$ \begin{split} P(h) =& \;0.560\;9 - 0.002\;41h(t) + \\ & 1.456\;88 \times {10^{ - 6}}{h^2}(t) + 3.885 \times {10^{ - 8}}{h^3}(t) \end{split}\;\;\; $$ (33) $$ \begin{split} P(h) = &\;1.126\;99 - 0.007\;3h(t) + \\ & 4.152\;1 \times {10^{ - 5}}{h^2}(t) - 1.175\;21 \times {10^{ - 7}}{h^3}(t) \end{split} $$ (34) 三次样条拟合模型表达式为
$$ {f_p}(h) = A + Bh(t) + C{h^2}(t) + D{h^3}(t) $$ (35) 根据数据拟合得到初始压力分别为0.3和0.6 MPa时空气弹簧作用力随高度的变化关系式如下
$$ \qquad\begin{split} {f_p}(h) = &2746.970\;3 - 31.923\;66h(t) + \\ & 0.270\;63{h^2}(t) - 9.801\;86 \times {10^{ - 4}}{h^3}(t) \end{split} $$ (36) $$\qquad \begin{split} {f_p}(h) =& 5509.067\;6 - 66.579\;49h(t) + \\ & 0.621\;33{h^2}(t) - 0.002\;31{h^3}(t) \end{split} $$ (37) 可以发现, 充入0.6 MPa时拟合得到的曲线形状和0.3 MPa的很类似, 两者相差2倍. 得到充入初始压力为x MPa时, 空气弹簧作用力随高度的变化关系式如下
$$ \begin{split}{f_p}(h) \approx & \frac{{{P_0}}}{{0.3}}\Big[ 2746.970\;3 - 31.923\;66h(t) + \Big.\hfill \\ & \Big. 0.270\;63{h^2}(t) - 9.801\;86 \times {10^{ - 4}}{h^3}(t) \Big] \end{split} $$ (38) 充入初始压力后, 假设充入0.6 MPa且达到静态平衡位置后, 任一时刻空气弹簧作用力和弹簧高度直接相关. 则充入初始压力P0后的弹簧刚度随高度的变化如下
$$\begin{split} k = &\frac{{{\rm{d}}{f_p}(h)}}{{{\rm{d}}h}} = {f_p}(h) \approx\frac{{{P_0}}}{{0.3}}\Big[ - 31.923\;66 + \Big.\\ & \Big. {\text{0}}{\text{.541\;26}}h(t) - {\text{2}}{\text{.940\;558}} \times {10^{ - 3}}{h^2}(t) \Big] \end{split}$$ (39) 由于高度为140 mm时的数据点相比于其余测试数据点偏差较大, 在拟合时将其屏蔽, 得到初始压力分别为0.3和0.6 MPa时空气弹簧有效面积随高度的变化关系式及有效面积变化率如下
$$ A(h) = 4140.103\;64 - 4.953\;45h(t) $$ (40) $$ A(h) = 4223.164\;85 - 3.093\;94h(t) $$ (41) 可以发现有效面积基本上都是关于弹簧高度线性变化的, 且和囊内初始气压大小无关. 综合两种初始气压实验结果, 因此可以获得空气弹簧有效面积变化率
$$ \frac{{{\rm{d}}A(h)}}{{{\rm{d}}h}} = \frac{{ - 4.953\;45 - 3.093\;94}}{2} \approx - 4 $$ (42) 综上所述, 可将空气弹簧也抽象称为双输入单输出的非线性模型, 即空气弹簧刚度大小和囊内初始气压及气囊高度有关, 囊内初始气压、气囊高度可根据悬架气压传感器和高度传感器进行实时测量得到, 最终得到弹簧刚度值.
4. 仿真验证与实物验证
4.1 参数化模型搭建
基于上文建立的悬架三维模型进行运动学仿真分析, 引入数据驱动下的参数化变量设计.三维物理模型输入运动学仿真环境后, 首先对悬架系统的关键硬点坐标进行确定, 这里包括: 上、下纵臂与转向节铰接位置、上下纵臂与转向臂铰接位置、减振器上下两端铰接位置共6处关键硬点. 通过对上述6个关键铰接位置的参数化变量设计, 可以得到多组不同数据对车轮垂向运动特性的影响对比, 较高的可调整性有利于最优位置参数的确定.双臂间距通过调整A和D两点的z坐标值; 臂长度通过调整A和B两点y坐标值; 减振器上铰接
$ O_{1} $ 处通过调整y坐标值; 减振器下铰接$ O_{2} $ 处通过调整x坐标值. 基础悬架空间位置坐标如表2.运动学问题仿真环境下, 进行对象融合时, 将零部件定义为刚体; 忽略运动副间摩擦. 由于设计初期将主销倾角定义为0, 故重点开展结构参数与车轮外倾角之间的关系, 同时对悬架工作过程中车轮动能以及纵向位移进行综合考虑.
表 2 基础悬架空间位置坐标值(X, Y, Z)Table 2. The coordinate value of the space position of the base suspension (X, Y, Z)Critical hard points X Y Z upper longitudinal arm front articulation point A −102.96 −8 115 lower longitudinal arm front articulation point B −102.96 −8 −65 lower longitudinal arm posterior articulation point C −54.21 −200 −65 upper longitudinal arm rear articulation point D −54.21 −200 115 hinge point up the shock absorber O1 −55.46 −108.5 488 hinge point under the shock absorber O1 −55.46 −70.5 −25 4.2 参数化增量分析
基于上文对减振器安装位置的初步确定, 将基础铰接位置参数进行调整, 其中减振器下铰接位置沿纵向调整3.3 mm. 首先展开对减振器安装角度对上下铰接位置的力学影响. 通过对上铰接位置做正、负增量调整来实现多角度状态对比, 得到结果如图14 .
根据减振器角度的增量调整后的结果可以得知, 做负调整时, 下铰接位置受力明显增加; 正调整及倾角为0处, 上下位置处受力相较于基础结果均有一定幅度的增加. 对于
$ O_{1} $ 处的受力变化曲线, 增量调整后虽在0 ~ 1.5 s时段内的受力低峰值有所降低, 但在2.0 s左右的高峰值却均有上升, 长期使用易导致减振器上下位置摆动的流畅性受到影响.随后, 基于上文对减振器安装角度的初步确定, 开展对A~D四处铰接位置的受力影响分析, 通过对双臂间距及长度增量调整, 得到四位置处的受力变化结果, 如图15 .
图15(a)~图15(d)分别对应A~D四处铰接位置的受力随增量调整的变化情况. 从中可以发现, 增加双臂的长度可以有效降低四个位置的受力; 减少臂间距虽可降低B, C两位置的受力, 但对A, D两位置的受力起到副作用, 并且降低的力远小于增加的力, 故不考虑调整臂间距. 增长双臂会导致车轮转弯时的灵活性受到影响, 并且影响减振器的上下端受力情况, 图16 为针对调整臂长, 给减振器受力带来的影响结果.
根据结果可以发现, 增加臂长虽然能一定幅度的降低A~D四处位置受力, 但会对减振器产生明显的附加垂向作用力, 这不利于减振器的正常工作, 长期使用会降低垂向调节行程. 同时, 在保证减振器倾斜角度的前提下, 增加臂长意味着臂间距随之增加, 故综合考虑, 造成四处铰接位置受力受到影响.
纵臂式悬架的合理设计可以在提升垂向调节范围的同时, 将轮胎侧向运动转变为纵向运动. 即搭载此悬架的车辆, 在底盘升降过程中, 前后轴车轮产生纵向相对运动, 从而缩短轴距, 进一步减小转弯半径, 提升机动灵活性. 因此, 悬架结构对车轮纵向移动量的影响作为重点研究对象, 经过前期开展对减振器铰接位置, 双纵臂间距、长度对车轮纵向位移的灵敏度分析, 确定出主要影响因素为臂间距和臂长度两者, 图17(a)为初步判断主要影响因素结果.
首先分别对臂间距以及臂长度进行增量调整, 为了确保准确性, 采用正负两方向对位置坐标进行改变. 可以发现, 0 ~ 1 s减振器压缩过程中, 车轮向整车质心方向发生位移, 对臂间距做正增量及臂长度做负增量调整时, 车轮位移均可达到50 mm以上; 1.5 ~ 2.5 s底盘举升过程中, 对臂长度及臂间距分别做负增量和正增量时, 车轮位移在基础位置结果的上下5 ~ 10 mm波动, 基于此结果, 结合考虑多变量的均一化调整影响整体准确定性, 故对参数化变量进行交叉融合分析, 得到图17(b)结果.
根据均一调整增量以及交叉融合两项迭代结果可以清晰发现, 对纵臂间距和长度作出调整时, 底盘降低和举升两个极限位置处, 车轮发生的纵向位移均有所偏颇, 易导致车轮做出不同程度的移动量、变形量, 相比于添加量调整增量的结果, 基础位置参数表现较为稳定, 峰值相对均衡.
悬架作为缓解路面冲击的关键系统, 通过导向机构来实现对外界力的传递、缓解, 因此, 在底盘升降过程中, 车轮动能变化可以准确反应出悬架性能表现. 首先分别通过对双臂间距及双臂长度进行正负增量调整, 得到图18 所示结果.
通过动能变化曲线反应的结果, 可以推断, 臂间距及臂长度无论从均一调整增量或是进行交叉融合分析, 都会导致车轮的动能变化出现波动, 其中, 调整双臂长度可以实现在0 ~ 1 s和2 ~ 3 s时段动能变化的相对平缓, 但1 ~ 2 s内, 相较于基础位置结果, 变化幅度略大, 在实际应用过程中, 车轮缓解冲击对整车的行驶平顺性起到决定性作用.
基于数据驱动模式, 采用参数化设计, 通过对增量的均一化、交叉化方式分别对车轮纵向位移以及车轮动能变化的影响因素进行综合比较, 最终确定出悬架系统关键铰接的位置坐标. 对于悬架系统性能优劣的判断, 应结合悬架动挠度、悬架动载荷等进行综合判断. 本文设计的悬架系统可有效解决车轮外倾的现象, 同时对侧向位移进行消除, 车轮在(−125 mm, 125 mm) 行程变化中的相关特性曲线如图19 .
为了确保垂向行程大幅提升的同时, 轮胎上胎面与转向臂不会发生运动干涉, 对此间距进行分析, 可以发现, 距离最近处为车轮动行程到最高点, 此处的间距保证了撞击悬架限位概率为0.
分别进行轮胎外倾角、悬架动挠度以及车轮纵-垂-侧三向速度变化分析可以发现, 悬架系统工作过程中, 轮胎外倾角虽出现微量波动, 但整体处于0位, 最大程度杜绝外倾现象; 悬架动挠度应充足, 避免颠簸路面带来的机构冲击; 通过路面对车轮垂-纵-侧作用力随轮跳行程的变化情况可以发现, 路面对轮胎侧向作用力为0, 保证了车轮在工作过程中始终保持垂直于地面, 轮胎接地面积实现最大化, 车轮不存在横向运动激励. 根据仿真结果可以得知, 实现垂向行程的同时, 降低车轮外倾的程度; 将车轮纵向位移保持在相对合理的空间, 提升车辆行驶性能, 强化了悬架系统对缓解路面不平度带来的冲击能力, 平缓的动能释放对行驶平顺性提供重要保障.
4.3 实物验证
为验证本研究悬架系统性能情况, 基于上文开展的理论基础工作, 搭建如图20的一体化电动轮, 为确保本研究可靠性, 利用如图21的全矢量线控平台进行实车验证. 通过控制执行拨杆实现对空气弹簧内气体的充、放, 以此达到底盘大行程升降性能要求, 得到电动轮安装位置处车架高度变化情况. 本研究主要针对车辆垂向动力学特性, 忽略车辆左右两侧的差异性, 故只针对单侧悬架系统进行动态测试, 根据高度传感器测试信号反馈, 得到如图22 ~ 图24实验结果.
通过对车辆左、右前悬架进行充、放气实验, 得车架高度变化, 结果可以发现: ①举升过程中, 车架高度变化平缓稳定, 证明本研究对悬架开展工作的有效性; ②前轴举升与下降过程中, 左后轮位置处的高度未发生变化, 表现出电动轮良好的独立性能.
作为高性能悬架系统, 需具备单侧车身高度可调节功能, 故对左侧前后两悬架系统进行充放气实验, 可以发现: ①左后轮在举升时首先介入对车架高度的变化影响, 下降时对其影响迟后于左前轮对应位置处; ②在同侧调节车身高度时, 变化过程中相对可靠, 高度稳定后产生微量波动, 处在允许范围内, 这是由于轮胎接地时地面产生法向作用力导致.
通过对所有悬架系统的升降触发, 实现整车通过性能的进一步提升, 得到如下实验结果, 结合前面两组结论, 可以发现: 车架高度调节范围在220 mm附近, 与初期测算理论范围(−125 mm ~ 125 mm)相比较, 较为合理, 同时发现, 左前轮在车身下降过程中的反应更为迅速, 这是由于车身质心到前轴范围内重量偏大, 同时也为后续开展车身姿态调节研究提供有效参考.
5. 结论
本文面向全线控平台, 结合多功能一体化电动轮结构特性, 考虑全线控平台的应用场景, 对主动悬架特性进行研究. 基于数据驱动背景, 开展悬架系统关键参数获取与辨识, 提出一种双纵臂式主动悬架系统方案, 可扩大垂向行程调节范围, 为搭载此悬架系统的全线控底盘实现大行程升降奠定基础.
针对 7 DOF 整车动力学模型进行搭建, 分析轮胎垂向运动与车身运动间的关系, 结合导向机构力学分析, 为悬架系统方案设计提供了理论依据. 利用数据驱动方式, 对悬架各项子系统进行特性测试与运动学分析, 并据此精确推算出三项系统的关键数据. 仿真测试与实车验证结果综合表明, 本研究提出的悬架系统在车身高度调节过程中的作用呈现正向影响, 满足高性能底盘的功能需求. 本项研究同时为后续开展车身姿态调节稳定性控制工作奠定基础.
-
表 1 空气弹簧测试结果
Table 1 Air spring test results
Initial
pressure/
MPaDownward pressure
displacement of the
exciter head/mmAir spring
force/
NIntracapsular
air pressure/
MPaEffective
area/
mm20.3 0 859 0.30 2863.3 10 1064 0.31 3432.3 20 1136 0.32 3550 30 1196 0.33 3624.2 40 1263 0.34 3714.7 50 1338 0.36 3716.7 60 1412 0.38 3715.8 70 1496 0.40 3740 80 1624 0.42 3866.7 90 1726 0.44 3922.7 100 1817 0.465 3907.5 0.6 0 1997 0.60 3328.3 10 2347 0.62 3785.5 20 2476 0.64 3868.8 30 2594 0.665 3900.8 40 2780 0.70 3971.4 50 2852 0.73 3906.8 60 2966 0.75 3954.7 70 3083 0.77 4003.9 80 3286 0.82 4007.3 90 3463 0.84 4122.6 100 3672 0.90 4080.0 表 2 基础悬架空间位置坐标值(X, Y, Z)
Table 2 The coordinate value of the space position of the base suspension (X, Y, Z)
Critical hard points X Y Z upper longitudinal arm front articulation point A −102.96 −8 115 lower longitudinal arm front articulation point B −102.96 −8 −65 lower longitudinal arm posterior articulation point C −54.21 −200 −65 upper longitudinal arm rear articulation point D −54.21 −200 115 hinge point up the shock absorber O1 −55.46 −108.5 488 hinge point under the shock absorber O1 −55.46 −70.5 −25 -
[1] 李亮, 王翔宇, 程硕等. 汽车底盘线控与动力学域控制技术. 汽车安全与节能学报, 2020, 11(2): 143-160 (Li Liang, Wang Xiangyu, Cheng Shuo, et al. Technologies of control-by-wire and dynamic domain control for automotive chassis. Journal of Automotive Safety and Energy Saving, 2020, 11(2): 143-160 (in Chinese) doi: 10.3969/j.issn.1674-8484.2020.02.001 Li Liang, Wang Xiangyu, Cheng Shuo, et al. Technologies of control-by-wire and dynamic domain control for auto-motive chassis. Journal of Automotive Safety and EnergySaving, 2020, 11(02): 143-160 (in Chinese) doi: 10.3969/j.issn.1674-8484.2020.02.001
[2] 陈国迎, 郑宏宇. 分布式全线控电动汽车的底盘集成控制. 华南理工大学学报(自然科学版), 2015, 43(11): 87-95 (Chen Guoying, Zheng Hongyu. Integrated chassis control of distributed drive-by-wire electric vehicles. Journal of South China University of Technology (Natural Science Edition) , 2015, 43(11): 87-95 (in Chinese) C-hen Guoy-ing, Zheng Hongyu. Integrated chassis control of distributed drive-by-wire electric vehicles. Journal of South China Univer-sity of Technology (Natural Science Edition), 2015, 43(11): 87-95 (in Chinese)
[3] 李韶华, 王伟达. 车辆动力学与控制研究进展. 动力学与控制学报, 2021, 19(3): 1-4 (Li Shaohua, Wang Weida. Research advance in vehicle dynamics and control. Journal of Dynamics and Control, 2021, 19(3): 1-4 (in Chinese) Li Shaohua, Wang Weida. Research adv-ance in vehicle dynamics and control. Jour-nal of Dynamics and Control, 2021, 19(3): 1-4. (in Chinese)
[4] 张友朋. 全线控四轮独立电动汽车状态估计与路径跟踪控制研究. [博士论文]. 长春: 吉林大学, 2019 Zhang Youpeng. Res-earch on state estimation and path following control of full X-by-wire electric vehicles. [PhD Thesis]. Changchun: Jilin University, 2019 (in Chinese)
[5] 姜鹏. 汽车悬架系统的仿真分析与参数优化设计. [博士论文]. 杭州: 浙江大学, 2006 Jiang Peng. The simulation study and parameter optimization on vehicle suspension system. [PhD Thesis]. Hangzhou: Zhejiang University, 2006 (in Chinese)
[6] 丁晓林, 王震坡, 张雷. 四轮轮毂电机驱动电动汽车驱动系统参数多目标优化匹配. 机械工程学报, 2021, 57(8): 195-204 (Ding Xiaolin, Wang Zhenpo, Zhang Lei. Powertrain sizing for four-wheel-independent-actuated electric vehicles based on multi-objective optimization. Journal of Mechanical Engineering, 2021, 57(8): 195-204 (in Chinese) Ding Xiaolin, Wang Zhenpo, Zhang Lei. Powe-rtrain sizing f-or four-wheel-independent-actuated electric vehicles based on multi-objective optimization. Journal of Mechanical Engineerin-g, 2021, 57(08): 195-204 (in Chinese
[7] 张泽星. 基于驾驶风格的个性化全线控电动汽车控制策略研究. [博士论文]. 长春: 吉林大学, 2020 Zhang Zexing. Study on personalized control strategy of full X-by-wireelectric vehicle based on driving style. [PhD Thesis]. Changchun: Jilin University, 2020 (in Chinese)
[8] 张冰. 全线控电动汽车不同行驶模式的轨迹跟踪控制研究. [博士论文]. 长春: 吉林大学, 2020 Zhang Bing Trajectory tracking control research on full X-by-Wire electric vehicle based on different maneuvering modes. [PhD Thesis]. Changchun: Jilin University, 2020 (in Chinese)
[9] 刘欢, 李韶华, 张培强. 刚柔耦合特种车辆越障行驶动力学分析及悬架优化. 动力学与控制学报, 2021, 19(3): 74-82 (Liu Huan, Li Shaohua, Zhang Peiqiang. Rigid flexible coupling special vehicle dynamic analysis and suspension optimization of obstacle crossing. Journal of Dynamics and Control, 2021, 19(3): 74-82 (in Chinese) Liu Huan, Li Shaohua, Zhang Peiqiang. Rigid flexi-ble coupling s-pecial vehicle dynamic analysis and suspe-nsionopt-imization o-f obstacle crossing. Journal of Dyna-mics and Control, 2021, 19(03): 74-82. (in Chinese)
[10] Mahmoodabadi MJ, Farhadi F, Sampour S. Firefly algorithm based optimum design of vehicle suspension systems. International Journal of Dynamics and Control, 20-19, 7(1): 134-146
[11] 周创辉. 一种新型液电馈能式汽车悬架系统的设计与研究. [博士论文]. 长沙: 湖南大学, 2018 Zhou Chuanghui. Design and study on a novel hydraulic-electrical regenerative suspension system for vehicles. [PhD Thesis]. Changsha: Hunan University, 2018 (in Chinese)
[12] Zhang L, Zhang Z, Wang Z, et al. Chassis coordinated control for full X-by-wire vehicles-A review. Chinese Journal of Mechanical Engineering, 2021, 34(1): 1-25 doi: 10.1186/s10033-020-00524-5
[13] 郭孔辉. 汽车操纵动力学. 长春: 吉林科学技术出版社, 1991 Guo Konghui. Vehicle Handling Dynamics. Changchun: Jilin Science and Technology Press, 1991 (in Chinese)
[14] 李韶华, 张兵, 黄玉亭. 刚柔耦合重型汽车建模及通过连续减速带的平顺性分析. 动力学与控制学报, 2018, 16(5): 397-402 (Li Shaohua, Zhang Bing, Huang Yuting. Modeling of rigid-flexible coupled heavy-duty vehicle and ride comfort analysis when passing through continuous speed bumps. Journal of Dynamics and Control, 2018, 16(5): 397-402 (in Chinese) Li Shaohua, Zhang Bing, Huang Yuting. Modeling of rigid-flexible coupled heavy-duty vehicle a-nd ride comfort analysis when passing through continuo-us speed bumps. Journal of Dynamics and Control, 201-8, 16(05): 397-402. (in Chinese))
[15] 李小彭, 李凡杰, 杨舲雪等. 车辆悬架系统的优化设计与动力学特性分析. 东北大学学报(自然科学版), 2020, 41(8): 1097-1102 (Li Xiaopeng, Li Fanjie, Yang Lingxue, et al. Optimization design and dynamics analysis of vehicle suspension system. Journal of Northeastern University (Natural Science) , 2020, 41(8): 1097-1102 (in Chinese) Li Xiaopeng, Li Fanjie, Yang Lingxue, et al. Optimizatio-n design and dynamics analysis of veh-icle suspension system. Journal of Northeastern Univers-ity (Natural Science), 2020, 41(08): 1097-1102 (in Chinese)
[16] 王震, 祝恒佳, 陈晓宇等. 基于交叉型双气室空气互联悬架的全地形车侧倾特性研究. 力学学报, 2020, 52(4): 996-1006 (Wang Zhen, Zhu Hengjia, Chen Xiaoyu, et al. Roll characteristics of all-terrain vehicle based on cross-type double-chamber air interconnection suspension. Chinese Journal of Theoretical and Applied Mechanics, 2020, 52(4): 996-1006 (in Chinese) Wang Zhen, Zhu Hengjia, Chen Xiaoyu, et al. Roll characteristics of all-terrain vehicle based on cross-type double-chamber air interconnection suspension. Chinese Journal of Theoretical and Applied Mechanics, 2020, 52(04): 996-1006 (in Chinese)
[17] 纪秀业, 侯献晓, 袁晓红. 基于ADAMS的某乘用车悬架系统匹配计算. 机械设计, 2021, 38(7): 90-95 (Ji Xiuye, Hou Xianxiao, Yuan Xiaohong. Matching calculation of the passenger vehicle’s suspension system based on AD-AMS. Journal of Machine Design, 2021, 38(7): 90-95 (in Chinese) Ji Xiuye, Hou Xianxia-o, Yuan Xiaohong. Matching calculation of the passenger ve-hicle’s suspension system based on AD-AMS. Journal of Mac-hine Design, 2021, 38(07): 90-95(in Chinese)
[18] 张智, 施晓芬, 李俊文. 汽车悬架系统的运动仿真及优化设计. 机械设计, 2015, 32(9): 30-33 (Zhang Zhi, Shi Xiaofen, Li Junwen. Kinematic simulation and structuraloptimization of vehicle suspension. Journal of Machine Design, 2015, 32(9): 30-33 (in Chinese) Zhang Zhi, ShiXiaofen, Li Junwen. Kinematic simulation and structuraloptimization of v-ehicle suspension. Journal of MachineDesign, 2015, 32(09): 30-33 (in Chinese)
[19] 王军年, 刘鹏, 杨钫等. 轮毂电机驱动电动汽车双横臂前悬架运动学优化. 汽车工程, 2021, 43(3): 305-312 (Wang Junnian, Liu Peng, Yang Fang, et al. Kinematics optimization of double wishbone front suspension for in-wheel motor-driven electric vehicles. Automotive Engineering, 2021, 43(3): 305-312 (in Chinese) W-ang Junnian, Liu Peng, Yang Fang, et al. Kinematics opt-imization of double wishbone front suspension for In-wheel Motor⁃Driven electric vehicles. Automotive Engineerin-g, 2021, 43(03): 305-312 (in Chinese)
[20] 李凡杰, 李小彭, 沃旭等. 混联汽车悬架系统减振性能分析与优化. 东北大学学报(自然科学版), 2021, 42(8): 1098-1104 (Li Fanjie, Li Xiaopeng, Wo Xu, et al. Analysisand optimization of damping performance of suspensionsystem of hybrid connected vehicle. Journal of Northeastern University (Natural Science) , 2021, 42(8): 1098-1104 (in Chinese) Li Fanjie, Li Xiaopeng, Wo Xu, et al. Analysisand optimization of damping performance of suspensionsystem of hybrid connec-ted vehicle. Journal of Northea-stern University (Natural Science), 2021, 42(08): 1098-1-104 (in Chinese)
[21] 袁宝峰, 王成恩, 邹猛等. 火星车主动悬架设计及蠕动脱困策略. 吉林大学学报(工学版), 2021, 51(1): 154-162 (Yuan Baofeng, Wang Chengen, Zou Meng, et al. Design of the active suspension of the mars rover and the strategy of creeping relief. Journal of Jilin University Engineering and Technology Edition) , 2021, 51(1): 154-162 (in Chinese) Yuan Baofeng, Wang Chengen, Zou Meng, et al. Design of the active suspension of the mars rover and the strategy of creeping relief. Journal of Jilin University Engineering and Technology Edition), 2021, 51(1): 154-1-62 (in Chinese)
[22] 张农, 王少华, 张邦基等. 液压互联悬架参数全局灵敏度分析与多目标优化. 湖南大学学报(自然科学版), 2020, 47(10): 1-9 (Zhang Nong, Wang Shaohua, Zhang Bangji, et al. Global sens-itivity analysis and multi-objective optimization of hydraulically interconnected suspension parameters. Journal of Hunan University (Natural Sciences) , 2020, 47(10): 1-9 (in Chinese) Zhang Nong, Wang Shaohua, Zhang Bangji, et al. Global sens-itivity analysis and multi-objective opti-mization of hydraulica-lly interconnected suspension para-meters. Journal of Hunan University (Natural Sciences), 2020, 47(10): 1-9 (in Chinese)
[23] 陈龙, 董红亮, 李利明. 适合轮毂电机驱动的新型悬架系统设计. 振动与冲击, 2015, 34(8): 174-180 (Chen Long, Dong Hongliang, Li Liming. A new type suspension design suitable for an in-wheel motor driving system. Journal of Vibration and Shock, 2015, 34(8): 174-180 (in Chinese) Chen Long, Dong Hon-gliang, Li Liming. A new type suspension de-sign suitable for an in-wheel motor driving system. Jou-rnal of Vibration and Shock, 2015, 34(08): 174-180 (in Chinese)
[24] 陆建辉, 周孔亢, 郭立娜等. 电动汽车麦弗逊前悬架设计及参数优化. 机械工程学报, 2012, 48(8): 98-103 (Lu Jianhui, Zhou Kongkang, Guo Lina, et al. Design and parametric optimization of mcPherson front suspension of electric vehicle. Journal of Mechanical Engineering, 2012, 48(8): 98-103 (in Chinese) LuJianhui, Zho-u Kongkang, Guo Lina, et al. Design and p-arametric optimiz-ation of mcPherson front suspension of electric vehicle. Journa-l of Mechanical Engineering, 201-2, 48(8): 98-103 (in Chinese)
[25] 李芳. 悬架结构参数优化及基于卡尔曼滤波的主动悬架控制研究. [博士论文]. 长春: 吉林大学, 2017 Li Fang. Optimization of suspension structure parameters and control strategy research of active suspension based on kalman filte. [PhD Thesis]. Changchun: Jilin University, 2017 (in Chinese)
[26] 李杰, 贾长旺, 成林海等. 轮毂电机电动汽车主动悬架约束状态反馈H(∞)控制. 吉林大学学报(工学版), 2022, 出版中, doi: 10.13229/j.cnki.jdxbgxb20210528 Li Jie, Jia Changwang, Cheng Linhai, et al. Constrained state feedback H(∞) control for active suspension of in-wheel motor electric vehicle. Journal of Jilin University (Engineering and Technology Edition), 2022, in press, doi: 10.13229/j.cnki.jdxbgxb20210528 (in Chinese)
[27] 潘惠惠. 汽车主动悬架系统的非线性控制研究. [博士论文]. 哈尔滨: 哈尔滨工业大学, 2017 Pan Huihui. Nonlinear control for vehicle active suspension systems. [PhD Thesis]. Harbin: Harbin Institute of Technology, 2017 (in Chinese)
[28] 方瑞华, 解跃青, 雷雨成. 空气悬架理论及其关键技术. 同济大学学报(自然科学版), 2003, 9: 1072-1076 (FangRuihua, Xie Yueqing, Lei Yucheng. Air suspension technology and developling trend. Journal of Tongji University, 2003, 9: 1072-1076 (in Chinese) doi: 10.3321/j.issn:0253-374X.2003.09.014 FangRuihua, Xi-e Yueqing, Lei Yucheng. Air suspension tech-nology and dev-elopling trend. Journal of Tongji Univers-ity, 2003(09): 1072-1076 (in Chinese)) doi: 10.3321/j.issn:0253-374X.2003.09.014
[29] 白先旭, 杨森. 磁流变半主动落锤冲击缓冲系统的“软着陆”控制试验与分析. 机械工程学报, 2021, 57(1): 121-127 Bai Xianxu, Yang Sen. Experimental test and analysis of “soft-landing” control for drop-induced shock systems using magnetorh-eological energy absorber. Journalof Mechanical Engineering, 2021, 57(1): 121-127 (in Chinese)
[30] 王震坡, 丁晓林, 张雷. 四轮轮毂电机驱动电动汽车驱动防滑控制关键技术综述. 机械工程学报, 2019, 55(12): 99-120 (Wang Zhenpo, Ding Xiaolin, Zhang Lei. Overview on key technologies of acceleration slip regulation for four-wheel-independently-actuated electric vehicles. Journal of Mechanical Engineering, 2019, 55(12): 99-120 (in Chinese) Wa-ng Zhenpo, Ding Xiaolin, Zhang Lei. Overvie-w on key tech-n-ologies of acceleration slip regulation forfour-wheel-independ-ently-actuated electric vehicles. Journ-al of Mechanical Engine-ering, 2019, 55(12): 99-120 (in Chinese)
[31] 李韶华, 罗海涵, 冯桂珍等. 电动汽车机-电-路耦合系统建模及动力学分析. 机械工程学报, 2021, 57(12): 51-61 (Li Shaohua, Luo Haihan, Feng Guizhen, et al. Modeling and dynamic analysis of mechanic-electro-road coupling system of electric vehicles. Journal of Mechanical Engineering, 2021, 57(12): 51-61 (in Chinese) Li Shaohu-a, Luo Haihan, Feng Guizhen, et al. Modelin-g and dynamic analysis of mechanic-electro-road coupli-ng system of electric vehicles. Journal of Mechanical E-ngineering, 2021, 57(12): 51-61 (in Chinese))
[32] 李韶华, 冯桂珍, 丁虎. 考虑胎路多点接触的电动汽车路面-耦合系统振动分析. 力学学报, 2021, 53(9): 2554-2568 (Li Shaohua, Feng Guizhen, Ding Hu. Vibration analysis of electric vehicle-road coupling system considering tire-road multipoint contact. Chinese Journal of Theoretical and Applied Mechanics, 2021, 53(9): 2554-2568 (in Chinese) Li Sha-ohua, Feng Guizhen, Dinghu. Vibration analys-is of electric ve-hicle-road coupling system considering t-ire-road multi-point c-ontact. Chinese Journal of Theoreti-cal and Applied Mechanic-s, 2021, 53(09): 2554-2568 (in Chinese))
[33] Öksüztepe E. In-wheel switched reluctance motor design for electric vehicles by using a pareto-based multiobjective differential evolution algorithm. IEEE Transactions on Vehicular Technology, 2016, 66(6): 4706-4715