摘要
针对高速公路自主换道场景,开展考虑驾驶员不满度换道决策和轨迹跟踪控制的辅助驾驶功能研究. 首先,在换道决策中考虑驾驶员心理因素,建立基于驾驶员不满度和最小安全距离的换道决策模型. 然后,基于跟踪误差模型建立横向线性二次调节器(Linear Quadratic Regulator,LQR)控制器,得到最优反馈前轮转角. 由于存在权重矩阵Q和R需要反复调试的问题,提出基于遗传算法(Genetic Algorithm, GA)的横向LQR控制器,获得最优的权重矩阵Q和R,提升算法的鲁棒性. 同时建立纵向LQR控制器保持期望车速. 利用PreScan、MATLAB/Simulink和CarSim联合仿真平台对自主换道策略进行验证,在不同固定车速换道工况下验证横向LQR控制器的鲁棒性和可靠性. 在换道安全性不满足工况下,期望车速能够被纵向LQR控制器稳定跟随,最大车速误差绝对值为0.69 km/h. 比较了LQR控制器和GA+LQR控制器的横向控制性能,采用GA+LQR控制器使横向和航向跟踪精度分别提升了66.7%和27%. 结果证实纵向LQR控制器和优化后的横向LQR控制器在不同换道工况下均具有良好的跟踪精度.
在高速公路交通流量日益增大的背景下,高速公路的交通安全面临更大压力.车辆换道是高速公路中较为普遍的驾驶行为,驾驶员对周围交通环境信息的错误判断而做出的换道操作往往容易引起交通事故.据统计,换道行为导致的交通事故占交通事故总数的4%~10%,因驾驶员对周围道路信息的误判而引起的交通事故占交通事故总数的75
车辆换道决策模型主要分为2种,基于人工智能方法和基于规则方法. 早期的车辆换道决策模型是Gipps模
轨迹跟踪控制是实现车辆安全行驶的重要环节,不仅要保证车辆安全、舒适地行驶,还要使车辆平稳地跟踪目标轨迹. 目前,控制方法多数基于车辆运动学或动力学制定相关控制策略从而得到控制量. 常用的控制算法有:比例-积分-微分(Proportional Integral Differentiation,PID)控制、模糊控制、滑模控制、模型预测控制(Model Predictive Control,MPC)、线性二次调节器(Linear Quadratic Regulator,LQR)等. Kusuma
本文针对上述问题,基于高速公路环境,提出一种基于驾驶员不满度的车辆自主换道决策模型,考虑驾驶员在实际驾驶中的心理因素,将决策模型类人化. 根据现有的驾驶员不满度模型和最小安全距离模型判断车辆换道的可行性,保证车辆在换道过程中的安全性. 基于二自由度车辆动力学模型设计基于跟踪误差最优反馈的LQR控制器,设计前馈控制消除曲率导致的稳态误差. 采用遗传算法对Q和R矩阵进行离线优化,得到Q和R矩阵的最优值,来提高传统LQR控制器的鲁棒性. 基于位置误差和速度误差设计纵向LQR控制器,最终通过纵横向LQR控制器实现车辆的纵横向协同控制. 在PreScan、MALTAB/Simulink和CarSim联合仿真平台上对所提出的换道决策模型和纵横向轨迹跟踪算法进行了验证.
1 换道决策模型
1.1 驾驶员不满度
车辆在行驶过程中,驾驶员需要根据前方车辆来调整自车车速. 当前方车辆持续低速,致使自车无法达到预期速度,驾驶员产生厌烦心理,容易激发其换道意图. 将上述驾驶员的厌烦心理称为驾驶员不满度,为了量化驾驶员不满度,采用速度不满的积累作为评价指
考虑驾驶员不满度在车辆行驶过程中的积累过程,当超过预定阈值时,驾驶员产生换道意图,驾驶员不满度表达式为:
(1) |
式中:为当前时刻驾驶员不满度;为上一时刻驾驶员不满度;为目标车速;为前车车速;T为采样时间;为增益系数,与前方车辆类型相关.在实际驾驶过程中,当前方车辆为大型车辆时,驾驶员更容易产生换道意图,驾驶不满度积累越快,取值更大;当前方车辆为小型轿车时,则反之.
最小安全距离指的是前方车辆突然刹车或前方有障碍物时,自车在刹车后,自车能够停在距前车一定的距离,保证不会发生追尾行为. 一般在停车时,自车会和前方障碍物保持一定距离,本文选取5 m.文献[
(2) |
式中:Db=0.012 2 V+0.058 5V
当车辆和前车的车距小于最小安全距离,并且跟车距离不断减小时,驾驶员不满度开始积累,当驾驶员不满度超过阈值Hthr时,驾驶员产生换道意图, 车辆执行换道操作,车辆开始执行换道后,驾驶员不满度清零.
(3) |
1.2 换道可行性分析
为了简化换道场景,只考虑交通流密度较低的高速公路环境,同时道路曲率连续且比较小,以确保车辆高速行驶.当驾驶员产生换道意图后,需要判断周围的车辆情况是否可以进行换道操作.在换道过程中,要考虑当前车道前车、目标车道前车以及目标车道后车对自车的影响.一般情况下,由于车辆速度不满而向左侧换道,因此,在考虑换道可行性时将道路简化为同向两车道.由于右换道和左换道相似,因此,本文方法也适用于多车道.
简化后车辆换道场景如

图1 简化后车辆换道场景
Fig.1 Simplified vehicle lane change scenario
在换道可行性分析中,通常用换道最小安全距离(Minimum Safety Spacing,MSS)模型,MSS表示在换道时刻需要和周围车辆保持的最小距离. MSS通过自车和周围车辆的距离、相对速度、加速度以及换道时间来判断换道是否安
在现有的研究中,车辆模型被简化为矩形模型、大圆模型、动态圆模型以及椭圆模型.本文采用

图2 椭圆车辆模型
Fig.2 Elliptical vehicle model
由于车辆的纵向长度对换道行为影响较大,因此采用椭圆长半轴Lx建立车辆模型,即
(4) |
式中:W为车辆宽度;L为车身长度;为后车的车速;为驾驶风格的系数. 值越小,椭圆的长半轴越长,说明所包含的车辆安全行驶范围越大,驾驶员越谨慎,需要更大的安全距离. 将驾驶员分为谨慎型、普通型和激进型,根据文献[
为了保证车辆换道安全,保证自车与自车道前车不发生斜向碰撞,因此需要满足:
(5) |
式中:为最小安全距离;为车辆M在换道过程中行驶的距离;为车辆L0在换道过程中行驶的距离;为车辆驾驶方向和车道线之间的夹角.
当车辆L0的车速大于车辆M的车速时,随着换道时间的增加,两车的间距会越来越大,换道行为不会引起两车碰撞, 此时车辆L0和车辆M的最小安全距离可设置为d,d一般取2~5
(6) |
式中:为换道时间.
加速换道的其他两种情形和上述匀速换道情形相似,因此,本文仅给出车辆M和车辆L0匀速换道过程中的最小安全距离模型的计算过程.
1.3 换道决策流程
换道决策流程如

图3 换道决策流程
Fig.3 Lane change decision process
2 换道跟踪控制
2.1 车辆动力学模型
针对车辆的横向运动控制,建立

图4 二自由度车辆模型
Fig.4 Two-degree-of-freedom vehicle model
考虑车辆侧向和横摆动力学特性,并假设轮胎侧偏特性处于线性范围,建立如下车辆动力学模型:
(7) |
式中:lf、lr分别为前、后轮到质心的距离;cf、cr分别为前、后轮侧偏刚度;分别为前、后轮的侧偏角;为横摆角速度;为前轮转角;IZ为绕Z轴的转动惯量,Z轴表示垂直方向,与XY平面垂直.
2.2 车辆横向误差模型建立
将建立的车辆动力学模型用误差变量重新定义,如

图5 道路信息
Fig.5 Road information
由于质心侧偏角较小,方向误差可近似为横摆角的误差,因此,横向误差和方向误差满足:
(8) |
式中:为道路曲率. 由于方向误差足够小,则横向误差变化率为:
(9) |
由于足够小,因此忽略此项. 以横向误差、横向误差变化率、方向误差、方向误差变化率为状态变量,以前轮转角为控制变量,为后续建立控制器方便描述,将状态空间写成如下形式:
(10) |
式中:
2.3 横向LQR反馈控制器设计
将连续系统离散化,由积分中值定理可得:
(11) |
式中:.
通常采用欧拉法将连续问题离散化,欧拉法一般可分为3种形式:向前欧拉法、中点欧拉法、向后欧拉法. 设定u(t)为当前时刻控制量,u(t+1)为下一时刻控制量,则u(t)是LQR控制算法需要求解的控制量,而u(t+1)为未知量,因此对中采用向前欧拉法,对中采用中点欧拉法,可得:
(12) |
式中:dt为采样时间,取0.01 s;I为单位矩阵. 得到离散后的微分方程为:
(13) |
式中:;.
根据所建立的车辆动力学模型,以横向速度、横摆角速度、横向误差以及方向误差为输入变量,设置控制变量u=,则LQR性能函数为:
(14) |
式中:Q和R分别为状态变量和控制变量的权重矩阵.Q矩阵中变量对应数值越大,表示该变量在中越重要,R矩阵中数值越大,表示控制变量的约束越大.
Q和R矩阵必须是正定或半正定的Hermit矩阵来确保优化的凸性,Q和R矩阵分别为:
, |
LQR控制是能找到一个控制输入,从而让LQR性能函数值最小,最终获得最优的前轮转角. 设存在唯一最优控制变量使LQR性能函数达到最值. 则表达式为:
. | (15) |
反馈增益矩阵K为:
. | (16) |
P为黎卡提方程的唯一正定解.黎卡提方程为:
(17) |
将状态空间和控制量代入
(18) |
2.4 横向前馈控制器设计
将控制律代入系统状态方程,可得:
(19) |
由
(20) |
因此,需要通过前馈补偿LQR反馈控制器所求出的最优控制量,使系统在趋于稳定的同时,系统的横向稳态误差也为0. 前轮转角的控制量为:
(21) |
式中:为前馈控制器输出的控制量.
当系统趋于稳定,为0,则稳态跟随误差为:
(22) |
可以得到和的表达式为:
(23) |
当为0时,前馈控制器的前轮转角控制量为:
(24) |
2.5 基于遗传算法的横向LQR控制器优化
遗传算法在车辆控制方面有诸多应用,Pizarro-Lerma
设定种群的数量为60,遗传进化代数为100,交叉概率为0.8,采用约束自适应变异,设置变异概率为0.09,适应度值函数偏差为1×1

(a) 适应度值

(b) 最佳个体值
图6 遗传算法迭代过程
Fig.6 Genetic algorithm iterative process
横向LQR算法优化前的Q、R矩阵分别为:
Q=, R=1 000
采用遗传算法优化横向LQR算法后,Q、R矩阵分别为:
QGA=
RGA=19 025.150 0
在仿真验证时,横向LQR算法均采用基于遗传算法优化后的Q、R矩阵.
需要解释的是,在设计遗传算法时,将横向误差和方向误差作为性能指标函数,并进一步将该性能指标函数转换为适应度函数进行约束,优化得到的Q、R矩阵使横向误差和方向误差往变小的方向进化. 虽然人为调整也能够起到改善控制精度的作用,但是不同权重之间影响关系较难确认,会导致调整比较困难,因此,采用遗传算法最终达到提升控制性能的目的.
2.6 纵向LQR控制器设计
本文采用LQR作为纵向控制器的控制算法,将实际轨迹和期望轨迹的车速误差和位置误差作为状态空间变量. 因此,建立状态空间并对其进行离散化处理,得到:
(25) |
式中:和分别为车辆在t时刻实际轨迹和期望轨迹的位置误差和车速误差;为加速度反馈控制量;为采样时间,取0.01 s.
定义,根据横向误差和方向误差最小化,可得到性能指标的表达式为:
(26) |
式中:为状态误差的权重矩阵;为控制量的权重矩阵; 为当前时刻加速度反馈控制量,采用2.4节所述方法可求得.当前时刻前馈控制量的表达式如下:
(27) |
式中:为车辆纵向加速度;为车辆横向速度;为横摆角速度. 最终加速度控制量
3 仿真验证与结果分析
3.1 联合仿真平台
本文设计了基于驾驶员不满意度的换道决策模块、基于五次多项式的轨迹规划模块、基于LQR的横、纵向控制模块. 基于PreScan、MATLAB/Simulink、CarSim搭建联合仿真平台,其中PreScan模拟高速公路场景,CarSim提供车辆动力学模型,利用MATLAB/Simulink搭建车辆自主换道决策、控制算法. 在PreScan中,车辆M周围的环境车辆设置1.2节所述的换道场景,如

图7 PreScan高速公路换道场景
Fig.7 Lane-change scenario of highway in PreScan
在PreScan中选用Audi_A8作为控制车辆,即M车,其整车参数如
参数 | 数值 |
---|---|
簧上质量m/kg | 1 820 |
前轮到质心距离/mm | 1 265 |
后轮到质心距离/mm | 1 682 |
轴距L/mm | 2 947 |
前轮侧偏刚度/(N·ra | -175 016 |
后轮侧偏刚度/(N·ra | -130 634 |
质心绕Z轴转动惯量/(kg· | 4 095 |
车宽W/mm | 2 040 |
设计的换道决策模块需要通过雷达探测车辆M与周围车辆的距离来判断是否满足换道的安全性, 因此,需要给车辆M配置相应传感器. 选用雷达传感器作为距离检测传感器.
首选设置长距离雷达,最大探测距离为150 m,水平探测角度为10°,根据当前场景,设置最大检测目标为2. 长距离雷达的探测角度相对较小,当两侧的车辆距离较近时,则无法探测到该车辆,因此需要配置短距离雷达来探测两侧车辆,设置最大探测距离为30 m,水平探测角度为80°,设置最大检测目标为2. 2种雷达的波形均为角锥形. 由于车辆在制动过程中会出现前倾的现象,因此垂直探测角不宜过小,在此处设置为10°. 由于换道过程中需要考虑与前、后车辆的距离,因此在车辆前方和后方均设置长距离和短距离雷达,共设置4个探测雷达,采样频率均为20 Hz.

(a) 车辆M和车辆L0位置关系

(b) 车辆M和车辆Ld位置关系

(c) 车辆M和车辆Fd位置关系
图8 换道过程中车辆间位置关系图
Fig.8 Relationship of each vehicle during lane-change process
3.2 自主换道仿真与结果分析
3.2.1 不同车速换道工况
本节基于前文所搭建的横纵向控制器验证自主换道的安全性及控制精度. 换道轨迹采用五次多项式换道轨迹,换道时间为4 s. 驾驶员不满意度0.02 s积累一次,即T为0.02 s.
本节暂不考虑驾驶员不满意度的影响,当车辆M和车辆L0之间的距离小于最小安全距离时,便执行换道.为了验证匀速换道情况下不同车速对横纵控制器的影响,设定车辆M分别以90 km/h、100 km/h、110 km/h的速度执行换道. 为了更好地表示换道安全性,设定车辆Fd和车辆Ld的车速与车辆M的车速相同,并保持匀速行驶. 车辆Fd初始时刻和车辆M相距60 m,车辆Ld初始时刻和车辆M相距30 m,车辆L0初始车速为80 km/h,初始时刻和车辆M相 距80 m.

(a) 横向误差
(b) 质心侧偏角

(c) 横摆角速度
(d) 横向加速度

(e) 车速误差
图9 不同车速自主换道仿真结果
Fig.9 Simulation results of autonomous lane-change at different vehicle speeds
3.2.2 换道安全性不满足工况
本工况主要考虑在产生换道意图时,安全性不满足的情况下驾驶员不满意度的变化, 同时关注产生换道意图时是否会出现消失或震荡,验证在跟踪换道轨迹时,遗传算法优化的有效性. 设置车辆M的初始车速为100 km/h,车辆Fd车速为110 km/h,并保持匀速行驶,初始时刻和车辆M相距50 m;车辆L0车速为80 km/h,并保持匀速行驶,初始时刻和车辆M相距100 m,目标车道前车Ld车速为100 km/h,初始时刻和车辆M相距30 m. 不满意度阈值为55.2. 车辆纵向距离与车速仿真结果如

(a) 车辆M与车辆L0纵向距离

(b) 车辆M与车辆L0车速
图10 车辆纵向距离与车速仿真结果
Fig.10 Vehicle longitudinal distance and speed simulation results
由
车辆换道决策仿真结果如

(a) 驾驶员不满意度与换道意图

(b) 换道可行性
图11 车辆换道决策仿真结果
Fig.11 Vehicle lane change decision simulation results
车辆轨迹跟踪仿真结果如

(a) 横向误差

(b) 方向误差

(c) 车速误差
图12 车辆轨迹跟踪仿真结果
Fig.12 Vehicle track tracking simulation results
不同控制器误差对比如
控制器 | 最大横向误差绝对值/m | 最大方向误差绝对值/rad | 平均横向误差绝对值/m | 平均方向误差绝对值/rad |
---|---|---|---|---|
LQR | 0.072 | 0.010 0 | 0.006 4 | 0.000 93 |
GA+LQR | 0.024 | 0.007 3 | 0.001 8 | 0.000 64 |
4 结 论
1)本文针对自主车辆换道,建立了基于驾驶员不满意度的换道决策模型.建立了基于车辆的道路误差模型的横向LQR控制器,提出了用遗传算法离线优化LQR控制器的方法.由于本文采用的遗传算法优化LQR控制器为离线优化,虽然获得最优权重矩阵时需要等待遗传算法优化结果,但对控制器的实时性并无影响. 同时建立纵向LQR控制器控制车速.
2)仿真结果表明,车辆在车速分别为90 km/h、100 km/h、110 km/h工况下,最大横向误差绝对值分别为0.028 m、0.034 m、0.054 m,验证了横向LQR控制器的鲁棒性和可靠性. 3种车速下,最大车速误差绝对值分别为0.24 km/h、0.35 km/h、0.35 km/h,车辆能够按照预期的车速行驶. 在换道安全性不满足工况下,LQR控制器和GA+LQR控制器最大横向误差绝对值分别为0.072 m和0.024 m,优化了66.7%;平均横向误差绝对值优化了71.9%. LQR控制器和GA+LQR控制器最大方向误差绝对值分别为0.010 0 rad和0.007 3 rad,优化了27%;平均方向误差绝对值优化了31.2%. 优化后的横向LQR控制器在不同车速和工况下的换道均具有较小的跟踪误差,纵向LQR控制器可以稳定跟随期望车速,具有较小的速度误差.
3)本文建立的自主换道决策和控制模型可以较好地实现换道行为,未来拟在建立考虑不同驾驶人员风格的换道决策模型基础上,进一步提升换道控制精度.
参考文献
AMMOUN S,NASHASHIBI F,LAURGEAU C.An analysis of the lane changing manoeuvre on roads:the contribution of inter-vehicle cooperation via communication[C]//2007 IEEE Intelligent Vehicles Symposium, June 13-15, 2007.Istanbul,Turkey:IEEE,2007:1095-1100. [百度学术]
GIPPS P G. A model for the structure of lane-changing decisions[J]. Transportation Research Part B:Methodological, 1986, 20(5): 403-414. [百度学术]
RATROUT N T, Rahman S M. A comparative analysis of currently used microscopic and macroscopic traffic simulation software[J]. Arabian Journal for Science and Engineering, 2009, 34: 121-133. [百度学术]
SEO J, YI K. Robust mode predictive control for lane change of automated driving vehicles[R]. Warrendale: SAE Technical Paper, 2015. [百度学术]
CHEN Y Y, LI G F, LI S, et al.Exploring behavioral patterns of lane change maneuvers for human-like autonomous driving[J].IEEE Transactions on Intelligent Transportation Systems, 2022, 23(9): 14322-14335. [百度学术]
TALEBPOUR A,MAHMASSANI H S,HAMDAR S H.Modeling lane-changing behavior in a connected environment:a game theory approach[J]. Transportation Research Procedia,2015,7:420-440. [百度学术]
WANG X,FAN J H,LIU N.A novel decision-making algorithm of autonomous vehicle based on improved SVM[C]//Proceedings of the 2020 International Conference on Aviation Safety and Information Technology,October 14-16,2020. Weihai, China:ACM,2020:643-648. [百度学术]
李斌. 汽车快速行驶换道场景分析与决策方法研究[D]. 长春:吉林大学, 2020. [百度学术]
LI B. Research on the scenario analysis and decision-making method of the lane change during fast driving[D]. Changchun: Jilin University, 2020. (in Chinese) [百度学术]
陈慧,王洁新.基于驾驶人不满度的高速公路自动驾驶换道决策[J]. 中国公路学报,2019, 32(12): 1-9. [百度学术]
CHEN H, WANG J X. A decision-making method for lane changes of automated vehicles on freeways based on drivers’ dissatisfaction[J].China Journal of Highway and Transport,2019,32(12):1-9.(in Chinese) [百度学术]
KUSUMA D H,ALI M,SUTANTRA N.The comparison of optimization for active steering control on vehicle using PID controller based on artificial intelligence techniques[C]//2016 International Seminar on Application for Technology of Information and Communication (ISemantic),August 5-6,2016. Semarang,Indonesia: IEEE, 2016: 18-22. [百度学术]
SAEEDI M A. A new robust combined control system for improving manoeuvrability, lateral stability and rollover prevention of a vehicle[J]. Proceedings of the Institution of Mechanical Engineers, Part K: Journal of Multi-body Dynamics, 2020, 234(1): 198-213. [百度学术]
李玉治, 李刚, 张志华. 极限工况下无人驾驶四轮转向汽车横向跟踪控制策略[J].重庆理工大学学报(自然科学), 2023,37(1): 66-74. [百度学术]
LI Y Z,LI G, ZHANG Z H. Lateral tracking control of a driverless four-wheel steering vehicle under extreme conditions[J]. Journal of Chongqing University of Technology (Natural Science),2023,37(1):66-74. (in Chinese) [百度学术]
YANG T,BAI Z W,LI Z Q,et al. Intelligent vehicle lateral control method based on feedforward+ predictive LQR algorithm[J].Actuators,2021,10(9): 228. [百度学术]
段敏,孙小松,张博涵.基于模型预测控制与离散线性二次型调节器的智能车横纵解耦跟踪控制[J].汽车技术,2022(8):38-46. [百度学术]
DUAN M,SUN X S,ZHANG B H.Horizontal and vertical decoupling tracking control based on MPC and DLQR[J].Automobile Technology,2022(8):38-46.(in Chinese) [百度学术]
高琳琳, 唐风敏, 郭蓬, 等. 自动驾驶横向运动控制的改进LQR方法研究[J].机械科学与技术, 2021, 40(3): 435-441. [百度学术]
GAO L L, TANG F M, GUO P,et al. Research on improved LQR control for self-driving vehicle lateral motion[J]. Mechanical Science and Technology for Aerospace Engineering, 2021, 40(3): 435-441.(in Chinese) [百度学术]
SEILER P, SONG B, HEDRICK J K. Development of a collision avoidance system[J]. Neurosurgery, 1998, 46(2): 492. [百度学术]
PHAN T T, LE L B, NGODUY D. A cooperative space distribution method for autonomous vehicles at a lane-drop bottleneck on multi-lane freeways[J]. IEEE Transactions on Intelligent Transportation Systems, 2022, 23(4): 3710-3723. [百度学术]
SUN B H,DENG W W,WU J,et al.Research on the classification and identification of driver’s driving style[C]//2017 10th International Symposium on Computational Intelligence and Design (ISCID), December 9-10, 2017. Hangzhou,China: IEEE, 2017: 28-32. [百度学术]
刘玉雯. 高速公路换道风险判别与个性化预警研究[D]. 西安: 长安大学, 2020: 10-26. [百度学术]
LIU Y W. Research on lane-changing risk discrimination and personalized forewarning of expressway[D]. Xi’an:Chang’an University, 2020:10-26. (in Chinese) [百度学术]
PIZARRO-LERMA A O,GARCÍA-HERNÁNDEZ R,SANTIBÁÑEZ V A. Fine-tuning of a fuzzy computed-torque control for a 2-DOF robot via genetic algorithms[J].IFAC-PapersOnLine, 2018, 51(13): 326-331. [百度学术]
郭景华, 胡平, 李琳辉, 等. 基于遗传优化的无人车横向模糊控制[J]. 机械工程学报, 2012, 48(6): 76-82. [百度学术]
GUO J H,HU P,LI L H,et al.Study on lateral fuzzy control of unmanned vehicles via genetic algorithms[J]. Journal of Mechanical Engineering, 2012, 48(6): 76-82.(in Chinese) [百度学术]
MARADA T,MATOUSEK R,ZUTH D.Design of linear quadratic regulator (LQR) based on genetic algorithm for inverted pendulum[J].MENDEL Soft Computing Journal,2017,23(1): 149-156. [百度学术]
SHIRAHATTI A,PRASAD P S S,PANZADE P,et al.Optimal design of passenger car suspension for ride and road holding[J].Journal of the Brazilian Society of Mechanical Sciences and Engineering,2008,30(1):66-76. [百度学术]