主权项: |
1.一种轮毂液压马达辅助驱动系统的模型预测控制方法,其特征在于,包括下列步骤: 第一步,计算变量泵和液压马达效率 液压泵组件包含主泵和补油泵,主泵采用斜盘式变量柱塞泵;变量泵效率受油液粘度、工作压力、斜盘开度及变量泵转速的影响,液压马达效率受油液粘度、工作压力及马达转速的影响,二者均为时变参量;假设液压油为牛顿流体,油液流动均为层流状态,忽略工作时部件间隙的变化和液体压缩性的影响,变量泵的容积效率和机械效率由下式(1)、(2)计算 式中:ηpv——变量泵容积效率 ηpm——变量泵机械效率 Cps——变量泵层流泄露系数 Cpv——变量泵层流阻力系数 Cpf——变量泵机械阻力系数 Tc——扭矩损失 μ——油液动力粘度 Dp——变量泵斜盘开度 Δp——变量泵进出口压差 np——变量泵转速 Vpmax——变量泵最大排量 液压马达容积效率和机械效率由下式(3)、(4)计算 式中:ηmv——液压马达容积效率 ηmm——液压马达机械效率 Cms——液压马达层流泄露系数 Cmf——液压马达机械阻力系数 Cmv——液压马达层流阻力系数 nm——马达转速 公式(1)、(2)、(3)、(4)中的油液动力粘度μ由下式(5)计算 μt=μt0e-λ(t-t0)··························(5) 式中:μt——液压油t温度下的动力粘度 μt0——油液在t=t0时的动力粘度 λ——液压油粘温系数 第二步,计算液压路径输入输出转矩 根据公式(6)计算得变量泵的理论输入转矩 式中:Tp_th——变量泵理论转矩 根据公式(7)计算得补油泵的理论输入转矩 式中:Tpc_th——补油泵理论转矩 Vpc——补油泵定排量 Δppc——补油泵进出口压差 利用公式(6)、(7),得到发动机传递至液压路径的扭矩 式中:TeH——发动机传递至液压路径的扭矩 ip——取力器速比 考虑液压系统容积效率的流量连续性原理,可得 DpnpVpmaxηpvηmv=2nmVm···················(9) 式中:Vm——马达排量 发动机转速ne和泵转速np的关系 np=ne/ip········(10) 忽略管道内的流量损失,利用公式(9)、(10)可得前轮转速与液压马达转速 式中:nf——前轮转速 系统中两个液压马达的总输出转矩为 利用公式(11)、(12),可得两个马达的总转矩为 根据现有公开技术中的前后轮速跟随思想,即保证搭载轮毂液压马达辅助驱动系统的车辆的前轮轮速与后轮轮速相等的控制思想,可得 nf=nr···············(14) 式中:nr——后轮转速 后轮轮速与发动机转速关系 nr=ne/igi0··························(15) 式中:ig——变速箱速比 i0——主减速器速比 利用式(13)、(14)、(15),液压马达总转矩可推导为 第三步,搭建整车传动系统动力学模型 发动机输出扭矩Te分别传递给前轴液压路径TeH和后轴机械路径TeM,即 Te=TeH+TeM······················(17) 利用式(8)可得发动机传递至机械路径的扭矩 后轮机械路径动力学方程为 前轮液压路径动力学方程为 式中:Tvr——后轮阻力矩 Jωr——后轮转动惯量 Cvr——后轮阻尼 ωr——后轮转速 Tvf——前轮阻力矩 Jωf——前轮转动惯量 Cvf——前轮阻尼 ωf——前轮转速 结合式(16)和式(18),得到传动系统动力学微分方程为 式中: 系统动力学模型的微分方程表达式为 式中,ηpv、ηpm、ηmv和ηmm分别通过理论公式(1)、(2)、(3)、(4)计算获取,其余各系数为定常系数; 第四步,建立预测模型 选取前后轮轮速ωf、ωr为状态变量,DpΔp、Te-VpcΔppc/(20πip)为控制变量,前后轮轮速差ωf-ωr为输出量,并对控制量u1、u2进行归一化处理,建立面向控制的系统状态空间方程 式中:yc=Δω=ωf-ωr 这里阻力矩Tvf、Tvr是可测扰动,可以通过现有技术获取;选择采样周期Ts,将系统模型式(23)离散化为 其中: 将模型改为获取控制增量的增量形式,如式(25)所示 第五步,推导预测输出方程 推导预测输出方程之前,需做以下假设 (1)控制时域之外,控制量不变 Δu(k+i|k)=0,i=u,u+1,......,p-1············(26) (2)可测干扰在k时刻之后不变 Δd(k+i|k)=0,i=1,2,......,p-1·······(27) 式中:p——预测时域 u——控制时域,且u≤p 根据k时刻检测到的系统状态x(k),计算预测系统未来动态的起点Δx(k),如式(28)所示,并根据Δx(k)来预测模型输出 Δx(k)=x(k)-x(k-1)······(28) 首先,计算模型状态的未来变化 接着,根据系统未来状态预测系统的未来输出 Y(k+1|k)=SxΔx(k)-SuΔU(k)+SdΔd(k)-Icy(k)······(30) 其中系数矩阵为 第六步,求解反馈控制律 定义参考序列 Rc(k+1)=[αΔω(k) α2Δω(k) α3Δω(k) ... ... αiΔω(k) ... αpΔω(k)]T(p)×1····(32) 根据轮毂液压马达辅助驱动系统的多目标需求制定优化问题 式中:Q、R为权重系数矩阵,性能指标J由以下两部分组成 (1)J1=||Q(Y(k+1|k)-Rc(k+1))||2将预测输出和参考值的差值最小作为优化控制的目标,该项使系统朝前后轮速差尽量小的方向优化,Q越大,轮速跟随响应越快; (2)J2=||RΔU(k)||2是包含控制量动作变化量的部分,用来惩罚控制动作的变化量,将控制量的动作变化率控制在较低范围作为优化目标,R越大,系统冲击越小; 求解上述如式(33)所示优化问题,实际上就是寻找最优的ΔU(k)使得性能指标J最小; 首先定义辅助量 然后将预测输出方程式(30)带入式(34),得到 其中, 所以上述性能指标J的表达式(33)可推导为 为了寻找使J最小的ΔU(k),对上式(37)求解一阶导数和二阶导数 由上式(38)可知,性能指标J的二阶导数大于0,则J的一阶导数为零的点即为使其取得最小值的解;因此上述优化问题的解,即k时刻的最优控制序列为 其中,Ep(k+1|k)由式(36)计算得到; 将上述最优控制序列的第一个分量作用于系统,得到系统闭环控制律为 定义预测控制增益 Kmpc=[I 0 ··· 0](SuTQTQSu+RTR)-1SuTQTQ··············(41) 则控制增量计算式为 ΔU(k)=KmpcEp(k+1|k)············(42) 其中,Ep(k+1|k)可由式(36)在线计算,Kmpc采用在线求解方式; 在下一采样时刻,获取新的观测值,又由预测方程重新计算未来输出,从而得到新的最优控制序列,如此重复上述优化计算过程。 |