CN105137758B - 电动助力转向***的多学科优化设计方法 - Google Patents
电动助力转向***的多学科优化设计方法 Download PDFInfo
- Publication number
- CN105137758B CN105137758B CN201510560843.1A CN201510560843A CN105137758B CN 105137758 B CN105137758 B CN 105137758B CN 201510560843 A CN201510560843 A CN 201510560843A CN 105137758 B CN105137758 B CN 105137758B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- optimization
- parameter
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Steering Control In Accordance With Driving Conditions (AREA)
Abstract
本发明涉及一种电动助力转向***的多学科优化设计方法,具体步骤是:1.车辆***建模,建立二自由度车辆模型,并对EPS模型的各部件进行受力分析,得到各自的运动方程;2.优化目标的建立,采用横摆角速度与转向角之比来表示转向灵敏度;3.EPS参数优化:(1)目标函数的确定,(2)优化算法:采用遗传算法和改进的粒子群算法,(3)优化步骤。本发明结合多学科协同优化对EPS***机械***和控制***的参数进行集成优化设计,在此基础上,以操纵路感和转向灵敏度为优化性能指标,分别运用遗传算法和粒子群算法两种优化方法进行优化,最终使EPS***得以全局优化,达到提高整车操纵稳定性的目的,实现了多科学多目标优化设计。
Description
技术领域
本发明涉及一种电动助力转向***该***,尤其是一种用于提高转向的轻便性和灵敏性,为驾驶者提供合适的路感的动助力转向***该***的优化设计方法。
背景技术
车辆操纵稳定性是汽车研究的一个重要方面。其中,转向***对驾驶员的力反馈对车辆操纵稳定性影响明显。由于电动助力转向***是复杂的机电耦合子***,其设计、参数优化以及与整车的匹配受到众多机械与控制参数的影响和制约,在EPS***的初始设计时,通常将机械结构参数和控制器性能参数单独设计。这样优化出的设计参数往往不能使转向***性能达到全局最优,从而影响车辆操稳性和行驶安全性。由于近年来,众多品牌的乘用车都出现了各种各样的EPS***问题,若能从***工程分析的角度对EPS进行全局优化设计,势必可以极大的降低故障发生的概率,而针对EPS***参数的优化设计选取必须通过综合考虑机械、控制参数的多学科全局优化方法得到,并通过整车操稳性的主客观评价方法进行匹配验证。
在传统的电动助力转向***设计中,先是通过理论来设计EPS的机械结构参数,然后采用相应的控制策略来设计控制器。这种设计方法把一个机电***的设计一分为二,虽然在前后两步设计过程中都使用了优化设计思想,但机械结构和电机控制器之间存在错综复杂的相互关系,使得机械结构与控制结构之间存在耦合,因此在设计之初就需要考虑这两者之间的关系,以得到全局最优参数。另一方面,由于被控对象与控制器之间又是相互联系的。因此对两个***同时优化设计,比单独设计更合理,也期望得到更佳的结果。EPS是一个高度协调的机电一体***,为了能够更加提高转向***性能并且减少能源消耗。
发明内容
本发明提出一种电动助力转向***的多学科优化设计方法,结合多学科协同优化对EPS***机械***和控制***的参数进行集成优化设计,在此基础上,以操纵路感和转向灵敏度为优化性能指标,分别运用遗传算法和粒子群算法两种优化方法进行优化,最终使EPS***得以全局优化,达到提高整车操纵稳定性的目的,实现了多科学多目标优化设计。这种通过机械参数和控制参数的多目标集成优化,为EPS***的优化提供一种新的解决方案,这是提出该策略的诱因。
本发明的技术方案是:一种电动助力转向***的多学科优化设计方法,包括以下具体步骤:
1.车辆***建模
建立二自由度车辆模型,并对EPS模型的各部件进行受力分析,得到各自的运动方程:转向柱:助力电机:输出轴:齿条:Tc=Ks(θh-θe) T1=K1δ Tm=Kaia Ta=Km(θm-N1δ) (5)
其中
由扭矩传感器测量的车速和转向柱力矩为输入量,通过试验得出助力脉谱特性,由助力电流计算出助力扭矩;
2.优化目标的建立
采用横摆角速度与转向角之比r(s)/θh(s)来表示转向灵敏度,即车辆方向盘转角到车辆横摆角速度响应的传递函数表示:
由公式1和EPS模型公式得转向灵敏度评价传递函数:
Uo(s)=Ks(Ls+R+KvN1Kt) (10)
其中:Xo2(s)=Cf(amVs+lCr)
式中:m为整车质量;V为车速;Cr,f分别为前后轮胎侧偏刚度;β为质心侧偏角;γ为横摆角速度;a为前轮到质心的距离;b为后轮到质心的距离;l=a+b;Iz为车身横摆转动惯量;N2为齿轮齿条转向器传动比;Kv为助力增益;Kt为电动机的电磁转矩常数。
求出操纵路感评价函数:
其中:
Th≈Ts=Ks(θh-θe)=-Ksθe (12)
根据劳斯判据判断稳定性,获得***的特征方程和各项系数,并将***特征方程的系数排成Routh表,并得出***稳定的必要条件;
3.EPS参数优化
(1)目标函数的确定:多目标优化问题的基本模型
min F(x)=(f1(x),f2(x),.....,f3(x))T (13)
L≤x≤U x=(x1,x2,...xn)
式中F(x)为目标函数,g(x)为不等式约束函数,h(x)为等式约束函数,x为决策矢量,L、U为x的上下界。
建立以操纵路感和转向灵敏度为目标函数的优化模型,操纵路感的频域能量:
转向灵敏度频域能量:
由于优化算法中进行优化时是使目标函数最小化,所以取优化目标函数为:
目标函数1:min-f1(Ks,Im,N1,Kp,Ki,Kd)
目标函数2:min-f2(Ks,Im,N1,Kp,Ki,Kd)
约束条件的确定,开始优化前根据分析确定参数的上下界:40≤Ks≤180 5≤N1≤30 0.0001≤Im≤0.001 0≤Kp≤20 0≤Ki≤10 0≤Kd≤10
(2)优化算法:
采用遗传算法和改进的粒子群算法,来寻找全局最优值,更新各粒子速度和位置的公式为:
vid t+1=wvid t+c1r1[pid t-xid t]+c2r2[average(pgd t)-xgd t] (16)
xid t+1=xid t+vid t+1 (17)
其中:w表示惯性权重;c1、c2为学习因子;r1、r2为介于0到1之间均匀分布的随机函数自动生成;t表示第t次迭代,d代表决策变量的维数,average(pgd t)代表对“记忆体”中当前最优Pareto解两两相互比较所得的较优粒子;
(3)优化步骤:
首先通过Simulink模型,求得一组名义***参数,即Ks,N1,Im,Kp,Ki,Kd的初始值,设定***在2秒内稳定,横摆角速度峰值在0.3左右,且增益幅度不超过20%的条件下,确定一组名义模型参数:
1)建立以操纵路感和转向灵敏度为目标函数的模型,将控制器参数取名义定值,利用两优化算法分别进行优化机械***参数Ks、N1、Im;
2)根据EPS控制和优化的需要,设计PID控制器控制策略,将机械结构参数取名义定值,优化控制器参数Kp、Ki、Kd;
3)将机械结构参数和控制***参数同时设为自变量,利用两种算法进行协同优化。代入simulink模型进行验证并比较最终优化效果。
本发明的有益效果是:本发明结合多学科协同优化对EPS***机械***和控制***的参数进行集成优化设计,在此基础上,以操纵路感和转向灵敏度为优化性能指标,分别运用遗传算法和粒子群算法两种优化方法进行优化,最终使EPS***得以全局优化,达到提高整车操纵稳定性的目的,实现了多科学多目标优化设计。
附图说明
图1是本发明集成优化框图;
图2a是EPS折线型助力特性图;
图2b是EPS助力特性试验Map图;
图3是机械参数多目标优化图;
图4是控制参数多目标优化图;
图5是多目标参数集成优化图;
图6是机械参数多目标粒子群优化结果图;
图7是控制参数多目标粒子群优化结果图;
图8是多目标集成多目标粒子群优化结果图;
图9是采用遗传算法的横摆角速度优化结果图;
图10是采用改进的粒子群算法的横摆角速度优化结果图;
图11是采用遗传算法的质心侧偏角优化结果图;
图12是采用改进的粒子群算法的质心侧偏角优化结果图;
图13是优化前后的转向灵敏度对比曲线图;
图14是优化前后的转向路感对比曲线图。
具体实施方式
下面结合附图与实施例对本发明作进一步说明。
如图1所示,本发明的电动助力转向***的多学科优化设计方法,包括以下具体步骤:
1.车辆***建模
建立二自由度车辆模型,并对EPS模型的各部件进行受力分析,得到各自的运动方程:转向柱:助力电机:输出轴:齿条:Tc=Ks(θh-θe) T1=K1δ Tm=Kaia Ta=Km(θm-N1δ) (5)
其中
合理的转向助力特性曲线不仅可保持汽车低速行驶时转向轻便灵活,而且可保持中高速行驶时的路感和操纵稳定性。该车初始设计采用的是折线型助力特性曲线图(如图2a),车速和转向柱力矩(由扭矩传感器测量)为输入量,通过试验得出助力脉谱特性,由助力电流计算出助力扭矩,实际助力Map特性如图2b所示。
2.优化目标的建立
在本发明中,为了更加直接客观的反映转向***和整车***的综合性能,拟采用横摆角速度与转向角之比r(s)/θh(s)来表示转向灵敏度,即车辆方向盘转角到车辆横摆角速度响应的传递函数表示。结合上面公式可知:
由公式1和EPS模型公式得转向灵敏度评价传递函数:
Uo(s)=Ks(Ls+R+KvN1Kt) (10)
其中:Xo2(s)=Cf(amVs+lCr)
式中:m为整车质量;V为车速;Cr,f分别为前后轮胎侧偏刚度;β为质心侧偏角;γ为横摆角速度;a为前轮到质心的距离;b为后轮到质心的距离;l=a+b;Iz为车身横摆转动惯量;N2为齿轮齿条转向器传动比;Kv为助力增益;Kt为电动机的电磁转矩常数。
求出操纵路感评价函数:
其中:
Th≈Ts=Ks(θh-θe)=-Ksθe (12)
根据劳斯判据判断稳定性,获得***的特征方程和各项系数,并将***特征方程的系数按下列形式排成Routh表,并得出***稳定的必要条件。
3.EPS参数优化
目标函数的确定:多目标优化问题的基本模型
min F(x)=(f1(x),f2(x),.....,f3(x))T (13)
L≤x≤U x=(x1,x2,...xn)
式中F(x)为目标函数,g(x)为不等式约束函数,h(x)为等式约束函数,x为决策矢量,L、U为x的上下界。
建立以操纵路感和转向灵敏度为目标函数的优化模型,操纵路感的频域能量:
转向灵敏度频域能量:
由于优化算法中进行优化时是使目标函数最小化,所以取优化目标函数为:
目标函数1:min-f1(Ks,Im,N1,Kp,Ki,Kd)
目标函数2:min-f2(Ks,Im,N1,Kp,Ki,Kd)
约束条件的确定,由于EPS优化的参数均为有实际意义的物理参数,故在开始优化前需确定其上下限。根据分析确定参数的上下界:40≤Ks≤180 5≤N1≤30 0.0001≤Im≤0.001 0≤Kp≤20 0≤Ki≤10 0≤Kd≤10
(1)优化算法:
遗传算法(genetic algorithm,GA)是一种进化优化算法,是将获取问题的最优解模拟成一个生物的进化过程来进行求解。其主要特点是:群体搜索策略和群体中个体之间的信息交换,搜索不依赖于信息梯度,具有良好的全局优化性能,较好地解决了传统算法易于先入局部最优的问题,可以求解目标函数比较复杂的问题。
粒子群算法(PSO)是一种基于群体的随机优化算法。它是通过粒子追寻搜索到的当前最优值来寻找全局最优值。单纯的粒子群优化很难解决多目标优化问题,本发明利用“记忆体”概念,通过不断记忆最优粒子,来寻找全局最优值,更新各粒子速度和位置的公式为:
vid t+1=wvid t+c1r1[pid t-xid t]+c2r2[average(pgd t)-xgd t] (16)
xid t+1=xid t+vid t+1 (17)
其中:w表示惯性权重;c1、c2为学习因子;r1、r2为介于0到1之间均匀分布的随机函数自动生成;t表示第t次迭代,d代表决策变量的维数,average(pgd t)代表对“记忆体”中当前最优Pareto解两两相互比较所得的较优粒子。
(2)优化步骤:
首先通过Simulink模型,求得一组名义***参数,即Ks,N1,Im,Kp,Ki,Kd的初始值。设定***在2秒内稳定,横摆角速度峰值在0.3左右,且增益幅度不超过20%的条件下,确定一组名义模型参数。
1)建立以操纵路感和转向灵敏度为目标函数的模型,将控制器参数取名义定值,利用两优化算法分别进行优化机械***参数Ks、N1、Im。
2)根据EPS控制和优化的需要,本文设计了PID控制器控制策略,将机械结构参数取名义定值,优化控制器参数Kp、Ki、Kd。
3)将机械结构参数和控制***参数同时设为自变量,利用两种算法进行协同优化。代入simulink模型进行验证并比较最终优化效果。
遗传算法优化:
如图3、4、5分别为机械结构参数、控制参数和多目标集成优化图,横轴为转向路感频域能量值,纵轴为转向灵敏度频域能量值。两能量值在取最优值时是相互冲突的,优化时就需要在两值之间进行适当的均衡处理,使两值均达到适度最优。MATLAB遗传算法工具箱在进行多目标优化时,优化结果常用Pareto图表示,Pareto前沿即Pareto最优化解的集合,当所得解越接近Pareto前沿,说明该解越趋于最优解。
表1遗传优化结果
根据图形的结果得到表格1,可以看出:集成优化得到的目标函数分别是机械参数优化的3.484倍、3.64倍;控制参数优化的5.554倍、3.602倍。
改进的粒子群算法优化:
选取种群个数为100,w取值在[0.5,1.5]之间,c1、c2取值为2,r1、r2为介于0到1之间均匀分布的随机函数自动生成。根据前面已经计算出的目标函数和优化参数,结合改进的粒子群程序,运行得到最终优化结果,如图6至图8所示。
表2改进的粒子群优化算法结果
根据图形的结果得到表格2,可以看出:集成优化得到的目标函数分别是机械参数优化的3.74倍、4.55倍;控制参数优化的5.39倍、4.40倍。
应用例:
关于转向操纵稳定性的实验评价,通常情况下在时域主要评价指标为质心侧偏角和横摆角速度。由于在实际工况中很难模拟方向盘阶跃响应,本发明以接近90度的斜率开始J-turn操作,获得如图9,10所示结果。
图9~12为遗传算法和改进的粒子群算法优化结果,与单独优化相比,集成优化后两者的超调量明显降低,***过渡时间缩短,更快地由瞬态达到稳态。
图13优化前后的转向灵敏度对比曲线图14优化前后的转向路感对比曲线
为了进一步从频域客观评价优化前后EPS对整车操纵稳定性的影响,文中采用了国际通用的操稳性四参数分析方法(Mimuro T,Ohsaki M,Yasunaga H,et al.Fourparameter evaluation method of lateral transient response[R].SAE TechnicalPaper,1990.),得到了图13和图14转向灵敏度幅频与相频特性曲线。
表3优化结果对比
优化指标 | 优化前 | 优化后(GA) | 优化后(OPS) |
横摆角速度增益(dB) | 0.189 | 0.257 | 0.274 |
共振峰频率(Hz) | 1.18 | 1.31 | 1.37 |
共振时增幅比(b/a) | 2.21 | 1.39 | 1.53 |
横摆角速度峰值/稳态横摆角速度=1/(2*ζ) | 0.226 | 0.359 | 0.327 |
f=0.1Hz时相位滞后角(°) | 0.023 | 0.003 | 0.010 |
f=0.6Hz时相位滞后角(°) | -22.4 | -14.34 | -14.30 |
f=1Hz时相位滞后角(°) | -30.1 | -25.3 | -24.1 |
由表3可知:优化后的汽车稳态横摆角速度增益适度增大、共振峰频率变高、增幅比降低、转向灵敏度响应带宽增大,***响应速度和车辆稳定性提高。同时,优化后***阻尼比增大,超调量减小,幅频特性曲线较平缓,通频带宽增加,相位滞后角绝对值适当减小,汽车的响应速度和***稳定性得到进一步提高。
此外,由图14转向路感幅频与相频特性曲线可知:
优化后的伯德曲线中频段增益平坦,相位和幅值裕度充足。相位裕度和带宽均增加,相位滞后减少,汽车操纵稳定性得到提高。
总之,通过机械与控制多学科多目标协同优化方法对EPS***进行集成优化,找到了最优的EPS***参数组合,最大限度降低了二者在运行情况下的相互作用和制约,使电动助力转向***实现全局最优性能。
Claims (1)
1.一种电动助力转向***的多学科优化设计方法,其特征在于,包括以下具体步骤:
一、车辆***建模
建立二自由度车辆模型,并对EPS模型的各部件进行受力分析,得到各自的运动方程:转向柱:助力电机:输出轴:齿条:Tc=Ks(θh-θe);T1=K1δ;Tm=Kaia;Ta=Km(θm-N1δ) (5)
其中
由扭矩传感器测量的车速和转向柱力矩为输入量,通过试验得出助力脉谱特性,由助力电流计算出助力扭矩;
二、优化目标的建立
采用横摆角速度与转向角之比γ(s)/θ(s)来表示转向灵敏度,即方向盘转角到车辆横摆角速度响应的传递函数表示:
<mrow>
<mfrac>
<mrow>
<mi>&gamma;</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&theta;</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&gamma;</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&delta;</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mi>&delta;</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&theta;</mi>
<mi>e</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mfrac>
<mrow>
<msub>
<mi>&theta;</mi>
<mi>e</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&theta;</mi>
<mi>h</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
由公式(1)和EPS模型公式得转向灵敏度评价传递函数:
<mrow>
<mfrac>
<mrow>
<mi>&gamma;</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&theta;</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>X</mi>
<mrow>
<mi>o</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>X</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mfrac>
<mn>1</mn>
<msub>
<mi>N</mi>
<mn>2</mn>
</msub>
</mfrac>
<mfrac>
<mrow>
<msub>
<mi>U</mi>
<mi>o</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>X</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>mI</mi>
<mi>z</mi>
</msub>
<msup>
<mi>Vs</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<mo>&lsqb;</mo>
<mi>m</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>a</mi>
<mn>2</mn>
</msup>
<msub>
<mi>C</mi>
<mi>f</mi>
</msub>
<mo>+</mo>
<msup>
<mi>b</mi>
<mn>2</mn>
</msup>
<msub>
<mi>C</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>I</mi>
<mi>z</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mi>f</mi>
</msub>
<mo>+</mo>
<msub>
<mi>C</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mi>s</mi>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>C</mi>
<mi>f</mi>
</msub>
<msub>
<mi>C</mi>
<mi>r</mi>
</msub>
<msup>
<mi>l</mi>
<mn>2</mn>
</msup>
</mrow>
<mi>V</mi>
</mfrac>
<mo>-</mo>
<mi>m</mi>
<mi>V</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>aC</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<msub>
<mi>bC</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>N</mi>
<mn>1</mn>
</msub>
<mn>2</mn>
</msup>
<msub>
<mi>I</mi>
<mi>m</mi>
</msub>
<mo>+</mo>
<msub>
<mi>I</mi>
<mi>e</mi>
</msub>
<mo>)</mo>
</mrow>
<msup>
<mi>s</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>N</mi>
<mn>1</mn>
</msub>
<mn>2</mn>
</msup>
<msub>
<mi>C</mi>
<mi>m</mi>
</msub>
<mo>+</mo>
<msub>
<mi>C</mi>
<mi>e</mi>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<msup>
<msub>
<mi>N</mi>
<mn>1</mn>
</msub>
<mn>2</mn>
</msup>
<msub>
<mi>K</mi>
<mi>b</mi>
</msub>
<msub>
<mi>K</mi>
<mi>a</mi>
</msub>
</mrow>
<mrow>
<mi>L</mi>
<mi>s</mi>
<mo>+</mo>
<mi>R</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mi>s</mi>
<mo>+</mo>
<msub>
<mi>K</mi>
<mi>s</mi>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>K</mi>
<mi>v</mi>
</msub>
<msub>
<mi>N</mi>
<mn>1</mn>
</msub>
<msub>
<mi>K</mi>
<mi>a</mi>
</msub>
<msub>
<mi>K</mi>
<mi>b</mi>
</msub>
</mrow>
<mrow>
<mi>L</mi>
<mi>s</mi>
<mo>+</mo>
<mi>R</mi>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
<mrow>
<mo>(</mo>
<mi>L</mi>
<mi>s</mi>
<mo>+</mo>
<mi>R</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mrow>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>Q</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mi>L</mi>
<mi>s</mi>
<mo>+</mo>
<mi>R</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
Uo(s)=Ks(Ls+R+KvN1Kt) (10)
其中:Xo2(s)=Cf(amVs+lCr);
式中:是Ta拉氏变换的结果;m为整车质量;V为车速;Cr,Cf分别为前后轮胎侧偏刚度;β为质心侧偏角;γ为横摆角速度;a为前轮到质心的距离;b为后轮到质心的距离;l=a+b;Iz为车身横摆转动惯量;N2为齿轮齿条转向器传动比;Kv为助力增益;Kt为电动机的电磁转矩常数;s为拉氏算子;
求出操纵路感评价函数:
其中:
Th≈Ts=Ks(θh-θe)=-Ksθe (12)
根据劳斯判据判断稳定性,获得***的特征方程和各项系数,并将***特征方程的系数排成Routh表,并得出***稳定的必要条件;
三、EPS参数优化
(1)目标函数的确定:多目标优化问题的基本模型
minF(k)=(f1(k),f2(k),.....,f3(k))T (13)
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>g</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo><</mo>
<mn>0</mn>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>.</mo>
<mo>,</mo>
<mi>m</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>h</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>.</mo>
<mo>,</mo>
<mi>n</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
</mtable>
</mfenced>
L≤k≤U;k=(k1,k2,...kn)
式中F(k)为目标函数,gi(k)为不等式约束函数,hj(k)为等式约束函数,k为决策矢量,L、U为k的上下界;
建立以操纵路感和转向灵敏度为目标函数的优化模型,操纵路感的频域能量:
<mrow>
<mi>max</mi>
<mo>:</mo>
<msub>
<mi>f</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>K</mi>
<mi>s</mi>
</msub>
<mo>,</mo>
<msub>
<mi>I</mi>
<mi>m</mi>
</msub>
<mo>,</mo>
<msub>
<mi>N</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>K</mi>
<mi>p</mi>
</msub>
<mo>,</mo>
<msub>
<mi>K</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msub>
<mi>K</mi>
<mi>d</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
<mi>&omega;</mi>
</mrow>
</mfrac>
<munderover>
<mo>&Integral;</mo>
<mn>0</mn>
<msub>
<mi>&omega;</mi>
<mn>0</mn>
</msub>
</munderover>
<mo>|</mo>
<mfrac>
<mrow>
<msub>
<mi>T</mi>
<mi>h</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<msubsup>
<mo>|</mo>
<mrow>
<mi>s</mi>
<mo>=</mo>
<mi>j</mi>
<mi>&omega;</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mi>d</mi>
<mi>&omega;</mi>
<mo>,</mo>
<mi>&omega;</mi>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<mn>0</mn>
<mo>,</mo>
<msub>
<mi>&omega;</mi>
<mn>0</mn>
</msub>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
</mrow>
转向灵敏度频域能量:
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mo>:</mo>
<msub>
<mi>f</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>K</mi>
<mi>s</mi>
</msub>
<mo>,</mo>
<msub>
<mi>I</mi>
<mi>m</mi>
</msub>
<mo>,</mo>
<msub>
<mi>N</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>K</mi>
<mi>p</mi>
</msub>
<mo>,</mo>
<msub>
<mi>K</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<msub>
<mi>K</mi>
<mi>d</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
<mi>&omega;</mi>
</mrow>
</mfrac>
<munderover>
<mo>&Integral;</mo>
<mn>0</mn>
<msub>
<mi>&omega;</mi>
<mn>0</mn>
</msub>
</munderover>
<mo>|</mo>
<mfrac>
<mrow>
<mi>&gamma;</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&theta;</mi>
<mi>h</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<msubsup>
<mo>|</mo>
<mrow>
<mi>s</mi>
<mo>=</mo>
<mi>j</mi>
<mi>&omega;</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mi>d</mi>
<mi>&omega;</mi>
<mo>,</mo>
<mi>&omega;</mi>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<mn>0</mn>
<mo>,</mo>
<msub>
<mi>&omega;</mi>
<mn>0</mn>
</msub>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>15</mn>
<mo>)</mo>
</mrow>
</mrow>
其中:Kp、Ki、Kd分别为比例增益参数、积分增益参数和微分增益参数;ω为频率;
由于优化算法中进行优化时是使目标函数最小化,所以取优化目标函数为:
目标函数1:min-f1(Ks,Im,N1,Kp,Ki,Kd)
目标函数2:min-f2(Ks,Im,N1,Kp,Ki,Kd)
约束条件的确定,开始优化前根据分析确定参数的上下界:40≤Ks≤180;5≤N1≤30;0.0001≤Im≤0.001;0≤Kp≤20;0≤Ki≤10;0≤Kd≤10
(2)优化算法:
采用遗传算法和改进的粒子群算法,来寻找全局最优值,更新各粒子速度和位移公式为:
vid t+1=wvid t+c1r1[pid t-xid t]+c2r2[average(pgd t)-xgd t] (16)
xid t+1=xid t+vid t+1 (17)
其中:w表示惯性权重;c1、c2为学习因子;r1、r2为介于0到1之间均匀分布的随机函数自动生成;t表示第t次迭代,d代表决策变量的维数,average(pgd t)代表对“记忆体”中当前最优Pareto解两两相互比较所得的较优粒子;vid,pid,xid,xgd分别表示第i个粒子在d维空间的速度,最优位置,位置和全局最优位置;
(3)优化步骤:
首先通过Simulink模型,求得一组名义***参数,即Ks,N1,Im,Kp,Ki,Kd的初始值,设定***在2秒内稳定,横摆角速度峰值在0.3左右,且增益幅度不超过20%的条件下,确定一组名义模型参数:
1)建立以操纵路感和转向灵敏度为目标函数的模型,将控制器参数取名义定值,利用两优化算法分别进行优化机械***参数Ks、N1、Im;
2)根据EPS控制和优化的需要,设计PID控制器控制策略,将机械结构参数取名义定值,优化PID控制器参数Kp、Ki、Kd;
3)将机械结构参数和控制***参数同时设为自变量,利用两种算法进行协同优化,代入simulink模型进行验证并比较最终优化效果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510560843.1A CN105137758B (zh) | 2015-09-06 | 2015-09-06 | 电动助力转向***的多学科优化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510560843.1A CN105137758B (zh) | 2015-09-06 | 2015-09-06 | 电动助力转向***的多学科优化设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105137758A CN105137758A (zh) | 2015-12-09 |
CN105137758B true CN105137758B (zh) | 2017-12-26 |
Family
ID=54723143
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510560843.1A Expired - Fee Related CN105137758B (zh) | 2015-09-06 | 2015-09-06 | 电动助力转向***的多学科优化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105137758B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106891992A (zh) * | 2017-01-16 | 2017-06-27 | 南京航空航天大学 | 一种复合转向***及其多目标优化方法 |
CN108595853B (zh) * | 2018-04-28 | 2021-06-04 | 吉林大学 | 一种基于遗传算法的并联混合动力汽车参数优化设计方法 |
CN108984860B (zh) * | 2018-06-27 | 2023-04-18 | 江苏大学 | 一种复合电源eps的电源参数优化方法 |
CN108942936A (zh) * | 2018-07-25 | 2018-12-07 | 重庆交通大学 | 一种机器人自适应可靠性优化方法 |
CN108828952B (zh) * | 2018-08-02 | 2021-03-05 | 中国民航大学 | 基于小脑神经网络的舵机电动加载***智能控制方法 |
CN109614703B (zh) * | 2018-12-11 | 2020-07-28 | 南京天航智能装备研究院有限公司 | 一种汽车电液复合转向***的多学科集成建模及优化方法 |
CN110615032B (zh) * | 2019-09-16 | 2021-10-01 | 南京航空航天大学 | 一种基于信息物理***的主动前轮转向参数优化方法 |
CN113276942B (zh) * | 2020-02-19 | 2022-07-26 | 广州汽车集团股份有限公司 | 一种后轮主动转向控制方法及其***、控制设备 |
CN113408156B (zh) * | 2021-08-18 | 2021-12-31 | 天津所托瑞安汽车科技有限公司 | 转向助力***结构参数的优化方法、电子设备 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101883708A (zh) * | 2007-07-27 | 2010-11-10 | 通用汽车环球科技运作公司 | 电动助力转向控制 |
CN102616269A (zh) * | 2011-01-25 | 2012-08-01 | 福特全球技术公司 | 齿条转向装置和装有齿条转向装置的伺服转向*** |
CN103562049A (zh) * | 2011-05-25 | 2014-02-05 | 三菱电机株式会社 | 电动助力转向的控制装置 |
CN104044631A (zh) * | 2013-03-15 | 2014-09-17 | 通用汽车环球科技运作有限责任公司 | 减小电动助力转向***中的方向盘振动的***和方法 |
CN104520168A (zh) * | 2012-08-20 | 2015-04-15 | Zf操作***有限公司 | 用于运行用于机动车的电动助力转向***的方法 |
CN104709348A (zh) * | 2013-12-11 | 2015-06-17 | 通用汽车环球科技运作有限责任公司 | 与eps控制器结合的碰撞避免控制 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9150246B2 (en) * | 2013-12-06 | 2015-10-06 | GM Global Technology Operations LLC | Algorithm for steering angle command to torque command conversion |
-
2015
- 2015-09-06 CN CN201510560843.1A patent/CN105137758B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101883708A (zh) * | 2007-07-27 | 2010-11-10 | 通用汽车环球科技运作公司 | 电动助力转向控制 |
CN102616269A (zh) * | 2011-01-25 | 2012-08-01 | 福特全球技术公司 | 齿条转向装置和装有齿条转向装置的伺服转向*** |
CN103562049A (zh) * | 2011-05-25 | 2014-02-05 | 三菱电机株式会社 | 电动助力转向的控制装置 |
CN104520168A (zh) * | 2012-08-20 | 2015-04-15 | Zf操作***有限公司 | 用于运行用于机动车的电动助力转向***的方法 |
CN104044631A (zh) * | 2013-03-15 | 2014-09-17 | 通用汽车环球科技运作有限责任公司 | 减小电动助力转向***中的方向盘振动的***和方法 |
CN104709348A (zh) * | 2013-12-11 | 2015-06-17 | 通用汽车环球科技运作有限责任公司 | 与eps控制器结合的碰撞避免控制 |
Also Published As
Publication number | Publication date |
---|---|
CN105137758A (zh) | 2015-12-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105137758B (zh) | 电动助力转向***的多学科优化设计方法 | |
Zhao et al. | Yaw and lateral stability control for four-wheel steer-by-wire system | |
Ren et al. | MPC-based yaw stability control in in-wheel-motored EV via active front steering and motor torque distribution | |
Wang et al. | Coordination control of differential drive assist steering and vehicle stability control for four-wheel-independent-drive EV | |
Xiong et al. | Vehicle dynamics control of four in-wheel motor drive electric vehicle using gain scheduling based on tyre cornering stiffness estimation | |
CN106372758B (zh) | 一种辅助泊车***的路径跟随方法及装置 | |
Zhang et al. | An analytical approach to improve vehicle maneuverability via torque vectoring control: theoretical study and experimental validation | |
CN107719372A (zh) | 基于动态控制分配的四驱电动汽车动力学多目标控制*** | |
CN108482363A (zh) | 车辆横摆稳定性预测模型控制方法 | |
Kang et al. | Coordinated vehicle traction control based on engine torque and brake pressure under complicated road conditions | |
Chen et al. | Dynamic drifting control for general path tracking of autonomous vehicles | |
CN103895704B (zh) | 基于后轮主动转向的变传动比控制方法 | |
CN107972667A (zh) | 一种车道偏离辅助***的人机协调控制方法及其控制*** | |
CN109911004A (zh) | 一种电动助力转向装置的回正控制方法及装置 | |
Cui et al. | Design optimization of vehicle EHPS system based on multi-objective genetic algorithm | |
Fahami et al. | Development of estimation force feedback torque control algorithm for driver steering feel in vehicle steer by wire system: Hardware in the loop | |
Zhang et al. | Human-centered torque vectoring control for distributed drive electric vehicle considering driving characteristics | |
Katsuyama et al. | A state-of-the-art review: toward a novel vehicle dynamics control concept taking the driveline of electric vehicles into account as promising control actuators | |
CN104908814A (zh) | 一种汽车线控转向***的分数阶pid控制方法 | |
Liu et al. | Research on Curved Path Tracking Control for Four‐Wheel Steering Vehicle considering Road Adhesion Coefficient | |
Ni et al. | Handling performance control for hybrid 8-wheel-drive vehicle and simulation verification | |
Chen et al. | Dynamics integrated control for four-wheel independent control electric vehicle | |
CN103439883A (zh) | 汽车底盘集成***神经网络广义逆解耦控制器及构造方法 | |
Zhang et al. | Cooperative game-based driver assistance control for vehicles suffering actuator faults | |
CN110968911B (zh) | 一种基于新型趋近律的汽车abs滑模控制器设计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171226 Termination date: 20200906 |
|
CF01 | Termination of patent right due to non-payment of annual fee |