主权项: |
1.一种基于综合效率最优的蠕行模式泵排量控制方法,其特征在于,包括以下步骤: 步骤一、泵斜盘开度搜索区间确定,根据当前车辆蠕行车速,确定车轮处需求功率,根据流量一致性原理求解得到轮毂液压马达转速,同时结合当前闭式回路液压压力差,确定变量泵斜盘开度的搜索区间; 根据流量一致性原理,蠕行模式下两个轮毂液压马达的流量等于变量泵输出流量,如式(1)所示: ωpβVpmaxηpvηvv=2ωmVm/ηmv (1) 式中,ωp、ωm分别表示液压变量泵转速以及轮毂液压马达转速,ηvv表式液压控制阀组及管道的效率损失; 此时,轮毂液压马达转速ωm与液压变量泵转速ωp之间满足关系: 同时,计算两个前轮轮毂液压马达的输出转矩和,如式(3)所示: 进一步,得到车轮处轮毂液压马达的输出功率,如式(4)所示: 此外,轮毂液压混合动力系统中发动机与变量泵之间通过PTO进行连接,因此发动机转速和变量泵转速之间满足以下关系: ωp=ωe/ip (5) 式中,ωe表示发动机转速,ip表示PTO速比; 因此,车轮处液压马达的输出功率与油路压力差ΔP以及当前发动机转速ωe的关系,如式(6)所示: 可见,蠕行模式下轮毂液压马达的功率输出主要取决于变量泵斜盘开度β、液压油路压力差ΔP以及发动机转速ωe,其中,油路压力主要取决于前轮的负载转矩,当车辆在蠕行模式下的目标行驶车速确定的情况下,此时车辆行驶需求功率以及需求转矩可以相应的确定,进而可以确定当前需求的液压油路压力差ΔP,那么根据式(6)所示的关系,通过调整变量泵的排量(即斜盘开度β)即可调节发动机转速ωe的工作区间; 基于任意允许的蠕行行驶车速vcreep,在确定了当前车轮处的需求功率Pcreep,轮毂液压马达转速ωm以及闭式回路液压压力差ΔP后,可以进一步确定变量泵斜盘开度的搜索区间根据发动机的最优转速区间[ωe,opt,min,ωe,opt,max],结合式(2),可以确定蠕行模式下对应当前蠕行车速的变量泵的最小斜盘开度βcreep,min与最大斜盘开度βcreep,max,如式(7)所示: 进而可以确定当前蠕行车速下对应的泵斜盘开度搜索区间,[βcreep,min,…,βcreep,n,…,βcreep,max]; 步骤二、发动机需求功率确定,对应任意的泵斜盘开度以及当前闭式回路液压压力差,结合变量泵与轮毂液压马达的效率计算公式,得到当前泵斜盘开度下对应的发动机需求功率; 利用变量泵与轮毂液压马达的效率计算公式,计算得到当前泵斜盘开度对应的液压变量泵的容积效率ηcreep,pv,n,机械效率ηcreep,pm,n,以及轮毂液压马达的容积效率ηcreep,mv,n,机械效率ηcreep,mm,n;进而得到当前泵斜盘开度下对应的发动机需求功率Pe,req,n,如式(8)所示: Pe,req,n=Pcreep/ηcreep,pv,n/ηcreep,pm,n/ηcreep,mv,n/ηcreep,mm,n (8) 步骤三、发动机工作点确定,基于当前轮毂液压马达转速和泵斜盘开度信息,求解得到不同泵斜盘开度对应的发动机转速,利用步骤二中的需求功率,求解得到当前泵斜盘开度对应的发动机转速转矩工作点; 步骤四、设置目标寻优函数,基于综合效率最优原则,选取当前泵斜盘开度对应的发动机油耗作为目标优化函数,目标优化函数最小的点即为系统综合效率最优的点; 基于综合效率最优的思想,迭代搜索得到当前蠕行车速下对应的最优泵斜盘开度以及最优发动机工作点,选取当前泵斜盘开度对应的发动机油耗Be,n作为目标函数,如式(9)所示,显然,目标函数最小的点即为系统综合效率最优的点,即最优泵排量控制目标以及发动机最优转速、转矩控制目标; Fn(βcreep,n)=Be,n=be,nPe,req,n (9) 式中,be,n表示当前发动机工作点下对应的燃油消耗率。 2.按照权利要求1所述的一种基于综合效率最优的蠕行模式泵排量控制方法,其特征在于,所述步骤五迭代寻优求解具体包括以下内容: 根据式(9)所示的目标寻优函数利用黄金分割算法迭代寻优,在泵斜盘开度搜索区间内快速求解当前蠕行车速对应的最优泵斜盘开度,具体如下所示: 1)选择已确定好的初始泵斜盘开度搜索区间,[a(1),b(1)]=[βcreep,min,1],设置搜索精度要求tol,黄金分割系数T=0.618,循环计数器初始值k=1; 2)令c(k)=a(k)+(1-T)(b(k)-a(k)),d(k)=b(k)-(1-T)(b(k)-a(k)),计算Fc=F(c(k)),Fd=F(d(k));其中,F即表示目标寻优函数; 3)若Fc<Fd,则转到步骤4);否则转到步骤5); 4)令a(k+1)=a(k),b(k+1)=d(k),d(k+1)=c(k),Fd=Fc; 令c(k+1)=a(k+1)+(1-T)(b(k+1)-a(k+1)),计算Fc=F(c(k+1)),转到步骤6); 5)令a(k+1)=c(k),c(k+1)=d(k),b(k+1)=b(k),Fc=Fd; 令d(k+1)=b(k+1)-(1-T)(b(k+1)-a(k+1));计算Fd=F(d(k+1)),转到步骤6); 6)令k=k+1;此时若满足b(k)-a(k)>=tol,则返回步骤3)继续迭代计算,直至迭代计算结果收敛;否则停止搜索。 |