一种基于线性时变的滑移率跟踪控制方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN108181811A(43)申请公布日2018.06.19(21)申请号201711363171.0(22)申请日2017.12.18(71)申请人长春工业大学地址130012吉林省长春市延安大街2055(72)发明人李绍松李政郭陆平韩玲崔高健卢晓晖于志新王国栋(51)Int.Cl.G05B13/04(2006.01)权利要求书3页说明书7页附图2页(54)发明名称一种基于线性时变的滑移率跟踪控制方法(57)摘要一种基于线性时变的滑移率跟踪控制方法，其特征在于，该方法包括期望滑移率模块、轮胎数据处理器、MPC控制器、Carsim汽车模型、滑移率计算模块。期望滑移率模块用于确定期望的滑移率；轮胎数据处理器用于确定轮胎纵向力与轮胎纵向刚度；MPC控制器根据当前时刻汽车纵向速度、滑移率，并结合期望的滑移率，优化求解出汽车轮胎的制动力矩，输出给CarSim汽车模型；CarSim汽车模型用于输出汽车的实际运动状态信息，包括汽车纵向速度、车轮转速；滑移率计算模块根据CarSim汽车模型输出的汽车纵向速度、车轮转速，计算得到滑移率。CN108181811ACN108181811A权利要求书1/3页1.一种基于线性时变的滑移率跟踪控制方法，其特征在于，该方法包括期望滑移率模块、轮胎数据处理器、MPC控制器、Carsim汽车模型、滑移率计算模块；期望滑移率模块用于确定期望的滑移率；轮胎数据处理器用于确定轮胎纵向力与轮胎纵向刚度；MPC控制器根据当前时刻汽车纵向速度、滑移率，并结合期望的滑移率，优化求解出汽车轮胎的制动力矩，输出给CarSim汽车模型；CarSim汽车模型用于输出汽车的实际运动状态信息，包括汽车纵向速度、车轮转速；滑移率计算模块根据CarSim汽车模型输出的汽车纵向速度、车轮转速计算得到滑移率；该方法包括以下步骤：步骤1、根据轮胎纵向力与滑移率的关系曲线，确定期望的滑移率：λref＝λp(1)其中：λp为轮胎纵向力最大值所对应的滑移率；步骤2、设计轮胎数据处理器，为了获得轮胎的非线性特性，基于Pacejka轮胎模型，获取不同路面附着系数下的轮胎纵向力与滑移率的关系曲线，得到轮胎纵向滑移特性三维图；获取不同路面附着系数下的轮胎纵向刚度与滑移率的关系曲线，得到轮胎纵向刚度三维图；轮胎数据处理器将实际的滑移率和路面附着系数分别输入到轮胎纵向滑移特性三维图和轮胎纵向刚度三维图，通过线性插值法分别获得轮胎纵向力和轮胎纵向刚度并输出给MPC控制器；在每个控制周期轮胎数据处理器更新一次轮胎纵向力和和轮胎纵向刚度数据；其中：Pacejka轮胎模型如下：其中：Fx是轮胎纵向力；λ是滑移率；Bx，Cx，Dx和Ex取决于轮胎垂直载荷Fz；b0＝1.57；b1＝35；b2＝1200；b3＝60；b4＝300；b5＝0.17；b6＝0；b7＝0；b8＝0.2；步骤3、设计MPC控制器，其过程包括如下子步骤：步骤3.1、建立车轮动力学模型与滑移率模型：其中：M为汽车质量；V为汽车纵向速度；R为轮胎半径；J为轮胎的转动惯量；Tb为轮胎的制动力矩；Fx为轮胎的纵向力；ω为轮胎的角速度；λ为滑移率；2CN108181811A权利要求书2/3页步骤3.2、建立预测模型，其过程包括如下子步骤：步骤3.2.1、预测模型的运动微分方程表达式为：考虑到汽车紧急制动工况下滑移率较大，轮胎纵向力随滑移率的增大而减小，两者呈现出非线性变化关系，为了表征轮胎纵向力与滑移率间的这种非线性变化特性，构建轮胎纵向力表达式如下：其中：其中：是轮胎的残余纵向力；是基于轮胎纵向滑移特性三维图，通过线性插值法获得的轮胎的纵向力；是基于轮胎纵向刚度三维图，通过线性插值法获得的轮胎纵向刚度；λ*为当前时刻的滑移率；最终得到预测模型的运动微分方程表达式为：将其写成状态空间方程，用于设计预测方程，具体如下：其中：状态变量x是滑移率；控制输入u是轮胎制动力矩；系统干扰输入d是轮胎的残余纵向力；式中状态矩阵A，控制输入矩阵Bu，干扰输入矩阵Bd，如下所示：步骤3.2.2、建立预测方程，用于预测系统未来输出；为了实现滑移率的跟踪控制，将连续时间的预测模型转换成离散时间系统的增量式模型：其中：取样时间k＝int(t/Ts)，t是仿真时间，Ts是仿真步长；C＝1；步骤3.3、设计优化目标及约束条件，其过程包括如下子步骤：步骤3.3.1、用期望的滑移率和实际的滑移率误差的二范数作为滑移率跟踪性能指标，体现滑移率跟踪特性，其表达式如下：其中：λref是期望的滑移率；λ是实际的滑移率；P是预测时域；k表示当前时刻；Q是加权因子；3CN108181811A权利要求书3/3页步骤3.3.2、用控制量变化率的二范数作为制动平滑指标，体现滑移率跟踪过