发明名称 硬件在回路仿真的车辆动力学模型的建模方法
摘要 本发明提供了一种硬件在回路仿真的车辆动力学模型的建模方法,包括驾驶员子模型的建模,在驾驶员子模型的建模中,对第一类道路,采用比例积分微分PID控制进行仿真,对第二类道路,采用基于与时间有关的常微分方程组的最优控制进行仿真,其中,所述最优控制的控制量是通过下述步骤得到的:采用复变函数的留数定理计算与所述常微分方程组对应的矩阵指数;根据所述矩阵指数计算所述常微分方程组的解;将得到的所述常微分方程组的解作为所述最优控制的控制量。本发明采用复变函数法求矩阵指数积分,提高了车辆模型在硬件在回路仿真中的精度和实时性。
申请公布号 CN101789038B 申请公布日期 2012.06.06
申请号 CN201010100875.0 申请日期 2010.01.25
申请人 北京经纬恒润科技有限公司 发明人 张明;薛俊亮;万亮;范志卿
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 北京康信知识产权代理有限责任公司 11240 代理人 余刚
主权项 1.一种硬件在回路仿真的车辆动力学模型的建模方法,包括驾驶员子模型的建模,其特征在于,在驾驶员子模型的建模中,对第一类道路,采用比例积分微分PID控制进行仿真,对第二类道路,采用基于与时间有关的常微分方程组的最优控制进行仿真;其中,所述最优控制的控制量是通过下述步骤得到的:采用复变函数的留数定理计算与所述常微分方程组对应的矩阵指数;根据所述矩阵指数计算所述常微分方程组的解;将得到的所述常微分方程组的解作为所述最优控制的控制量,其中,所述步骤包括:根据牛顿第二定律确定的整车侧向和横摆运动的方程如下:<maths num="0001"><![CDATA[<math><mrow><msub><mi>m</mi><mi>Vehicle</mi></msub><msub><mover><mi>v</mi><mo>&CenterDot;</mo></mover><mrow><mi>y</mi><mo>,</mo><mi>Vehicle</mi></mrow></msub><mo>=</mo><msub><mi>F</mi><mrow><mi>y</mi><mo>,</mo><mi>Wheel</mi><mo>,</mo><mi>Front</mi></mrow></msub><mi>cos</mi><msub><mi>&delta;</mi><mi>Steering</mi></msub><mo>+</mo><msub><mi>F</mi><mrow><mi>y</mi><mo>,</mo><mi>Wheel</mi><mo>,</mo><mi>Rear</mi></mrow></msub><mo>-</mo><msub><mi>v</mi><mrow><mi>x</mi><mo>,</mo><mi>Vehicle</mi></mrow></msub><msub><mover><mi>&gamma;</mi><mo>&CenterDot;</mo></mover><mi>Vehicle</mi></msub><msub><mi>m</mi><mi>Vehicle</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msub><mi>J</mi><mi>Vehicle</mi></msub><msub><mover><mi>&gamma;</mi><mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mrow></mover><mi>Vehicle</mi></msub><mo>=</mo><msub><mi>aF</mi><mrow><mi>y</mi><mo>,</mo><mi>Wheel</mi><mo>,</mo><mi>Front</mi></mrow></msub><mi>cos</mi><msub><mi>&delta;</mi><mi>Steering</mi></msub><mo>-</mo><msub><mi>bF</mi><mrow><mi>y</mi><mo>,</mo><mi>Wheel</mi><mo>,</mo><mi>Rear</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,m<sub>Vehicle</sub>-整车质量,<img file="FDA0000083002590000013.GIF" wi="133" he="54" />-质心侧向加速度<img file="FDA0000083002590000014.GIF" wi="107" he="49" />-车辆的横摆角速度,J<sub>Vehicle</sub>-车辆的横摆转动惯量,<img file="FDA0000083002590000015.GIF" wi="108" he="49" />-车辆的横摆角加速度,δ<sub>steering</sub>-前轮转角,v<sub>vehicle</sub>-车辆质心的速度,v<sub>x</sub>-质心速度沿车身纵向的分量,F<sub>y,Wheel,Front</sub>-前轮所受地面的侧向力,F<sub>y,Wheel,rear</sub>-后轮所受地面的侧向力,a-质心到前轴的距离,b-质心到后轴的距离;将二自由度运动方程线性化,所述运动方程为标准的状态空间方程组:<maths num="0003"><![CDATA[<math><mrow><mover><mi>x</mi><mo>&CenterDot;</mo></mover><mo>=</mo><mi>Ax</mi><mo>+</mo><mi>Bu</mi></mrow></math>]]></maths>y=Cx+Du    (3)其中,x-状态向量,该向量是[车辆质心的侧向位移,车辆质心的侧向速度,车辆的横摆角速度,车辆的横摆角],y-输出向量,该向量是[车辆质心的侧向位移],u-控制输入向量,该向量是[前轮转角],A、B、C、D-整理成标准状态空间方程过程中生成的系数距阵;所述式(3)的离散时间解是:y(kT)=CΦ(kT)x(0)+CΓ(kT)Bu    (4)其中:Φ(kT)=e<sup>AkT</sup>,<img file="FDA0000083002590000022.GIF" wi="386" he="133" />k-第k个采样点,T-采样周期,x(0)-初始时刻的状态向量;采用复变函数的留数定理计算与常微分方程组对应的矩阵指数Φ(kT)=e<sup>AkT</sup>,其中,所述常微分方程组为所述式(3);根据所述矩阵指数计算所述常微分方程组的解y(kT)=CΦ(kT)x(0)+CΓ(kT)Bu;根据目标函数的极值求控制量u的过程如下:所述目标函数:<img file="FDA0000083002590000023.GIF" wi="342" he="114" />其中m为计算选取的参考采样点数;目标函数极值处:<maths num="0004"><![CDATA[<math><mrow><mfrac><mi>dJ</mi><mi>du</mi></mfrac><mo>=</mo><mn>0</mn><mo>,</mo></mrow></math>]]></maths>从而求出最优的控制输入u,使得目标函数取极值;将得到的所述最优的控制输入u作为最优控制的控制量。
地址 100101 北京市朝阳区安翔里23号商业楼1-2层
您可能感兴趣的专利