发明名称 一种基于粒子群算法的车辆悬架系统参数优化方法
摘要 本发明提供一种基于粒子群算法的车辆悬架系统参数优化方法,属于汽车底盘技术领域,同时涉及群体智能算法领域。根据不同车型或某一车型对悬架系统不同使用性能的要求,折中考虑悬架参数优化过程中的各种矛盾,采用搜索能力较强的粒子群算法对悬架系统的参数进行寻优,建立作为寻优判别标准的适应度函数,得到粒子种群的历史全局最优值,作为满足所需求解的悬架系统参数。该方法算法简单,收敛速度快,鲁棒性强,目标函数灵活可调,为悬架系统参数的优化设计提供了一种新的思路。
申请公布号 CN103646280B 申请公布日期 2016.03.02
申请号 CN201310615996.2 申请日期 2013.11.28
申请人 江苏大学 发明人 陈龙;施德华;江浩斌;汪若尘
分类号 G06N3/00(2006.01)I 主分类号 G06N3/00(2006.01)I
代理机构 南京知识律师事务所 32207 代理人 汪旭东
主权项 一种基于粒子群算法的车辆悬架系统参数优化方法,其特征在于包括以下步骤:步骤1、将汽车简化为四分之一单轮车辆模型,根据车身质量m<sub>b</sub>,车轮质量m<sub>t</sub>,车身垂直位移z<sub>b</sub>、车身垂直速度<img file="FDA0000852476540000011.GIF" wi="59" he="62" />和车身垂直加速度<img file="FDA0000852476540000012.GIF" wi="78" he="63" />车轮垂直位移z<sub>t</sub>、车轮垂直速度<img file="FDA0000852476540000013.GIF" wi="59" he="54" />和车轮垂直加速度<img file="FDA0000852476540000014.GIF" wi="69" he="55" />悬架刚度k<sub>s</sub>,轮胎刚度k<sub>t</sub>,减振器阻尼系数c以及路面不平度z<sub>0</sub>,建立二自由度车辆动力学微分方程;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>m</mi><mi>b</mi></msub><msub><mover><mi>z</mi><mo>&CenterDot;&CenterDot;</mo></mover><mi>b</mi></msub><mo>+</mo><msub><mi>k</mi><mi>s</mi></msub><mrow><mo>(</mo><msub><mi>z</mi><mi>b</mi></msub><mo>-</mo><msub><mi>z</mi><mi>t</mi></msub><mo>)</mo></mrow><mo>+</mo><mi>c</mi><mrow><mo>(</mo><msub><mover><mi>z</mi><mo>&CenterDot;</mo></mover><mi>b</mi></msub><mo>-</mo><msub><mover><mi>z</mi><mo>&CenterDot;</mo></mover><mi>t</mi></msub><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mrow>]]></math><img file="FDA0000852476540000015.GIF" wi="796" he="103" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>m</mi><mi>t</mi></msub><msub><mover><mi>z</mi><mo>&CenterDot;&CenterDot;</mo></mover><mi>t</mi></msub><mo>+</mo><msub><mi>k</mi><mi>s</mi></msub><mrow><mo>(</mo><msub><mi>z</mi><mi>t</mi></msub><mo>-</mo><msub><mi>z</mi><mi>b</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>k</mi><mi>t</mi></msub><mrow><mo>(</mo><msub><mi>z</mi><mi>t</mi></msub><mo>-</mo><msub><mi>z</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>+</mo><mi>c</mi><mrow><mo>(</mo><msub><mover><mi>z</mi><mo>&CenterDot;</mo></mover><mi>t</mi></msub><mo>-</mo><msub><mover><mi>z</mi><mo>&CenterDot;</mo></mover><mi>b</mi></msub><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mrow>]]></math><img file="FDA0000852476540000016.GIF" wi="1085" he="95" /></maths>步骤2、根据汽车对悬架系统偏频的要求,确定悬架刚度k<sub>s</sub>的范围,记k<sub>s</sub>的取值范围的最小值为k<sub>smin</sub>,k<sub>s</sub>的取值范围的最大值为k<sub>smax</sub>;步骤3、建立粒子群算法的适应度函数,优化目标为在满足约束条件的可行域内使得适应度值最小,当粒子不满足约束条件时,采用惩罚函数将该粒子的当前位置剔除;适应度函数用车身垂直加速度均方根值和轮胎动载荷均方根值的加权和表示,记为:J=q<sub>1</sub>a<sub>w</sub>+q<sub>2</sub>DTL<sub>rms</sub>其中,J为适应度值;a<sub>w</sub>为车身垂直加速度均方根值;DTL<sub>rms</sub>为轮胎动载荷均方根值;q<sub>1</sub>为车身加速度均方根值的加权系数,q<sub>2</sub>为轮胎动载荷均方根值的加权系数;步骤4、设定粒子群算法参数,初始化粒子群,其中搜索空间维度N表示悬架系统待优化的参数个数,取N=2,第一维表示减振器阻尼系数,第二维表示悬架刚度系数;群体规模为M,i表示粒子序数,i=1,2,…,M;第i个粒子的位置状态表示为lz<sub>i</sub>=(lz<sub>i1</sub>,lz<sub>i2</sub>),速度状态表示为lv<sub>i</sub>=(lv<sub>i1</sub>,lv<sub>i2</sub>),最大速度v<sub>max</sub>,每个粒子个体最优位置为per<sub>i</sub>=(per<sub>i1</sub>,per<sub>i2</sub>),全局最优位置为global=(global<sub>1</sub>,global<sub>2</sub>),设定最大迭代次数为K,随机产生粒子位置状态和速度状态的初始值,设定惯性权重w的初始值和最终值,学习因子c<sub>1</sub>、学习因子c<sub>2</sub>的初始值和最终值;步骤5、迭代次数为k,k=1,2,…,K;在第k次迭代中,根据适应度函数,求出每个粒子当前的适应度值<img file="FDA0000852476540000017.GIF" wi="78" he="66" />据此判断当前粒子位置的适应度值是否优于历史最优值J<sub>imin</sub>,如果是,则该粒子的历史最优位置per<sub>i</sub>被当前位置所取代,在当前迭代次数k中,找出M个粒子历史最优适应度值中的最优解Jbest=(J<sub>1min</sub>,J<sub>2min</sub>,…,J<sub>Mmin</sub>),并将其对应的位置记为全局最优global;步骤6、对每个粒子进行位置状态和速度状态更新,若更新后的速度超过设定的最大速度v<sub>max</sub>,将其约束为v<sub>max</sub>;步骤7、位置状态和速度状态更新后,使迭代次数k=k+1,若迭代次数不超过K,返回步骤(5)重复执行操作;当迭代次数大于K时,迭代终止,优化过程结束;步骤8、根据优化得到的结果确定悬架系统的性能指标,采用该性能指标对悬架系统参数进行优化设计;所述步骤3中,约束条件和惩罚函数的建立方法如下:约束条件:当以舒适性为导向时,q<sub>1</sub>&gt;0且q<sub>2</sub>=0,悬架撞击限位块的概率为0.3%,车轮跳离地面的概率为0.3%,阻尼比ζ范围为(0.2,0.45),悬架刚度系数范围为(k<sub>smin</sub>,k<sub>smax</sub>);当以安全性为导向时,q<sub>1</sub>=0且q<sub>2</sub>&gt;0,车身垂直加速度超过0.8g的概率为0.3%,悬架撞击限位块的概率0.3%,阻尼比ζ范围为(0.2,0.45),悬架刚度系数范围为(k<sub>smin</sub>,k<sub>smax</sub>);当兼顾平顺性和舒适性时,q<sub>1</sub>≠0且q<sub>2</sub>≠0时,设定不同的q<sub>1</sub>和q<sub>2</sub>可以使悬架系统达到不同的综合性能,悬架撞击限位块的概率为0.3%,阻尼比ζ范围满足(0.2,0.45),悬架刚度系数范围为(k<sub>smin</sub>,k<sub>smax</sub>);惩罚函数:当解不满足约束条件时,则令粒子当前位置的适应度值远远大于其历史最优适应度值。
地址 212013 江苏省镇江市京口区学府路301号