CN110239362B - 一种分布式电驱动车辆多性能优化力矩分配方法 - Google Patents
一种分布式电驱动车辆多性能优化力矩分配方法 Download PDFInfo
- Publication number
- CN110239362B CN110239362B CN201910530350.1A CN201910530350A CN110239362B CN 110239362 B CN110239362 B CN 110239362B CN 201910530350 A CN201910530350 A CN 201910530350A CN 110239362 B CN110239362 B CN 110239362B
- Authority
- CN
- China
- Prior art keywords
- torque
- optimization
- equation
- mpc
- vehicle
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000005457 optimization Methods 0.000 claims abstract description 54
- 238000004364 calculation method Methods 0.000 claims abstract description 22
- 239000013598 vector Substances 0.000 claims description 17
- 239000011159 matrix material Substances 0.000 claims description 14
- 238000005070 sampling Methods 0.000 claims description 8
- 230000004888 barrier function Effects 0.000 claims description 7
- 230000001133 acceleration Effects 0.000 claims description 6
- 238000013461 design Methods 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 5
- 230000000903 blocking effect Effects 0.000 claims description 2
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 claims description 2
- 239000010931 gold Substances 0.000 claims description 2
- 229910052737 gold Inorganic materials 0.000 claims description 2
- 230000002441 reversible effect Effects 0.000 description 7
- 238000011160 research Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 241000764238 Isis Species 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000011423 initialization method Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60L—PROPULSION OF ELECTRICALLY-PROPELLED VEHICLES; SUPPLYING ELECTRIC POWER FOR AUXILIARY EQUIPMENT OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRODYNAMIC BRAKE SYSTEMS FOR VEHICLES IN GENERAL; MAGNETIC SUSPENSION OR LEVITATION FOR VEHICLES; MONITORING OPERATING VARIABLES OF ELECTRICALLY-PROPELLED VEHICLES; ELECTRIC SAFETY DEVICES FOR ELECTRICALLY-PROPELLED VEHICLES
- B60L15/00—Methods, circuits, or devices for controlling the traction-motor speed of electrically-propelled vehicles
- B60L15/20—Methods, circuits, or devices for controlling the traction-motor speed of electrically-propelled vehicles for control of the vehicle or its driving motor to achieve a desired performance, e.g. speed, torque, programmed variation of speed
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/60—Other road transportation technologies with climate change mitigation effect
- Y02T10/72—Electric energy management in electromobility
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Software Systems (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- General Engineering & Computer Science (AREA)
- Operations Research (AREA)
- Mechanical Engineering (AREA)
- Transportation (AREA)
- Power Engineering (AREA)
- Computing Systems (AREA)
- Feedback Control In General (AREA)
- Electric Propulsion And Braking For Vehicles (AREA)
Abstract
本发明公开了一种分布式电驱动车辆多性能优化力矩分配方法,采用基于改进的C/GMRES算法的MPC控制器进行多性能优化力矩分配,包括:MPC控制器以上层控制器给出的纵向总驱动力矩及横摆力矩作为力矩分配控制输入,结合车辆反馈的轮胎相关信息,基于车轮旋转动力学及轮胎滑移率状态更新模型,考虑轮毂电机输出力矩物理特性,以车轮负荷率为稳定性约束前提,以轮胎滑移耗散能最小化为优化目标,建立MPC优化控制问题,并采用改进的C/GMRES算法进行力矩优化分配问题的快速求解,以实现车辆稳定性与经济性的实时优化控制。本发明采用基于改进的C/GMRES算法的MPC控制器进行多性能优化力矩分配优化,在保证车辆稳定性的前提下最小化轮胎滑移耗散能,计算效率高,适于实时应用。
Description
技术领域
本发明涉及电动车辆及车辆动力学控制领域,特别涉及一种分布式电驱动车辆(Distributed Drive Electric Vehicle,DDEV)的多性能优化力矩分配方法。
背景技术
电动汽车(EV)凭借其低排放、快速响应速度及高能量利用率等优点,已成为当前汽车工业的重要发展方向,其中一种特殊的底盘结构便是分布式电驱动车辆。与传统的集中驱动EV不同,DDEV在各个车轮中搭载轮毂电机(In-wheel Motor,IWM)作为主要的动力输出单元,可有效提升车辆稳定性及控制灵活性。然而,由于引入了更多的控制自由度,DDEV实际上是一个过驱动控制***,即控制自由度(力矩分配)大于车辆自由度,这为控制器的设计带来了新的挑战。
广大学者们已对分布式电驱动车辆力矩分配问题开展了大量的研究,其研究问题主要分为两类:安全需求和能耗优化。然而,目前普遍研究的力矩分配目标仅为单一的安全需求或能量优化,极少研究同时考虑***安全及能量进行力矩分配。此外,大部分考虑能量优化的力矩分配方法仅以电机能效为优化目标,而忽略了轮胎滑移耗散能量的影响。
常用的分布式电驱动车辆力矩分配方法有基于规则的分配、菊花链分配、重配置伪逆分配和基于优化的分配。基于规则的分配方法鲁棒性好,易于实车应用,但没有充分利用分布式车辆力矩独立驱动优势,也不能很好地处理执行器及路面条件相关约束;菊花链式分配无法保证车轮充分且均匀地被利用,会产生所谓的“短板效应”,即当一个轮胎已接近约束边界,其余轮胎可能依然处于较低负荷水平;伪逆法不能有选择地优先满足纵向驱动力或者横摆力矩。相比上述方法,基于优化的分配法可以更好地设计优化指标及约束,全面提升力矩分配品质。然而,目前普遍采用的基于优化的方法仅考虑当前反馈的状态量信息,模型动态如驱动器及轮胎动态往往不考虑,这将极大地限制控制性能。模型预测控制(Model Predictive Control,MPC)是一种结合开环优化与闭环反馈的先进算法,可在短期预测时域内以状态更新的形式考虑***未来动态信息,以此进行优化,实现更理想的控制效果。然而,由于需要对未来状态进行优化,MPC控制器优化求解中不容小觑的计算量是一大制约其实际应用的阻碍。
由此可见,上述现有的分布式驱动车辆力矩分配方法,显然仍存在有不便与缺陷,而亟待加以进一步改进。如何设计一种兼顾车辆安全及轮胎滑移耗散能最小化的分布式电驱动车辆力矩分配方法,在保证其控制性能的基础上满足实时计算需求,成为当前业界急需改进的目标。
发明内容
本发明要解决的技术问题是提供一种新形式的基于优化的分布式电驱动车辆力矩分配方法,该方法可同时兼顾车辆稳定性并实现轮胎滑移耗散能最小,在保证其控制性能的基础上满足实时计算需求。
为解决上述技术问题,本发明采用如下技术方案:
一种分布式电驱动车辆多性能优化力矩分配方法,采用基于改进的C/GMRES(延拓/广义最小残差,Continuation/Generalized Minimal Residual)算法的MPC控制器进行多性能优化力矩分配,包括:MPC控制器以上层控制器给出的纵向总驱动力矩及横摆力矩作为力矩分配控制输入,结合车辆反馈的轮胎相关信息,基于车轮旋转动力学及轮胎滑移率状态更新模型,考虑轮毂电机输出力矩物理特性,以车轮负荷率稳定性为约束前提,以轮胎滑移耗散能最小化为优化目标,建立MPC优化控制问题,并采用改进的C/GMRES算法进行力矩优化分配问题的求解。
作为本发明进一步地改进,所述改进的C/GMRES算法中,通过阻碍方程的形式将不等式约束等效转换为目标方程代价,实现不等式约束问题求解。
进一步地,所述改进的C/GMRES算法中,采用基于KKT条件的初始化计算方法求解。
进一步地,所述MPC控制器的控制问题设定过程中将分布式电驱动车辆的前后轴轮毂电机输出力矩约束进行整合。
进一步地,所述轮胎滑移率状态更新模型中的轮胎滑移率更新状态方程表示为:
其中x=[κ1 κ2 κ3 κ4]T为状态量,u=[T1 T2 T3 T4]T为控制量;κi为第i个轮胎的滑移率,Ti为第i个轮毂电机的输出力矩,i=1,2,3,4分别表示前左、前右、后左和后右;
所述MPC控制器的控制问题表示为:
其中,Np表示预测时域长度;为瞬时代价方程;xo为当前时刻***反馈的状态量数值,;ρ为阻碍方程的权重标量;为阻碍方程;x=[κ1 κ2κ3 κ4]T为状态量,u=[T1 T2 T3 T4]T为控制量;κi为第i个轮胎的滑移率,Ti为第i个轮毂电机的输出力矩,i=1,2,3,4分别表示前左、前右、后左和后右;to表示当前采样时刻;为***状态方程。
进一步地,基于所述MPC控制器的控制问题构建哈密尔顿方程,计算预测时域内状态量更新及协状态量;根据庞德里亚金极小值原理最优性条件写出矩阵函数F(U(to),xo(t)),其中xo为当前时刻***反馈的状态量数值,U(to)为当前时刻MPC控制器待优化向量,to表示当前采样时刻;求解计算获得最优解。
进一步地,所述矩阵函数表示为:
其中,表示关于目标方程及***状态方程写出的哈密尔顿方程;Np为预测时域长度;x=[κ1 κ2 κ3 κ4]T为状态量,κi为第i个轮胎的滑移率,Ti为第i个轮毂电机的输出力矩,i=1,2,3,4分别表示前左、前右、后左和后右;λ为协状态向量;xo为当前时刻***反馈的状态量数值;U(to)为当前时刻MPC控制器待优化向量;to表示当前采样时刻;基于xo及U,以有限差分的形式计算即
进一步地,所述改进的C/GMRES求解的初始化计算方法为:
S101、基于车辆反馈的轮胎相关信息,通过获得最优解其中Q=diag(q1 q2 q3 q4)及R=diag(r1 r2r3 r4)均表示对角权重矩阵, Js为轮胎转动惯量,Ttot为车辆总驱动力矩需求,ΔM为车辆外加横摆力矩需求,rw为车轮有效半径,ds为车轮轴间距,为第i个车轮的旋转加速度,i=1,2,3,4分别表示前左、前右、后左和后右;
通过采用上述技术方案,本发明至少具有以下优点:
1、本发明采用基于改进的C/GMRES算法的MPC控制器进行多性能力矩分配优化,在保证控制性能的基础上,该算法计算效率高,适用于实时应用。
2、本发明中将轮胎约束进行合并,将预测时域内任意时刻的优化变量个数减半,且消去了原始问题的等式约束,简化了计算。
3、针对传统C/GMRES算法无法求解不等式约束问题,本发明将不等式约束通过阻碍方程等效转换为优化目标代价,实现不等式约束问题求解。
4、针对C/GMRES算法求解,本发明提出了一种基于KKT条件的初始化计算方法,该方法可快速获得C/GMRES算法初值,可提升C/GMRES的适用性。
附图说明
上述仅是本发明技术方案的概述,为了能够更清楚了解本发明的技术手段,以下结合附图与具体实施方式对本发明作进一步的详细说明。
图1是本发明实施例中分布式电驱动车辆多性能优化力矩分配方法的控制框架图;
图2是本发明实施例中车辆动力学关系图;
图3是本发明实施例中改进的C/GMRES算法快速计算流程图。
具体实施方式
本发明使用改进的C/GMRES算法解决的MPC控制问题,在保证控制性能的基础上,实现快速计算;针对传统C/GMRES算法无法实现不等式约束的问题,采用阻碍方程将不等式约束等效转换为目标方程代价,满足不等式约束边界限制;针对C/GMRES算法初值选取问题,提出了C/GMRES的初值快速计算方法。
下面结合附图和具体实施例对本发明进行详细说明。本实施例以两轴轮毂电机驱动车辆为控制对象,以本发明技术方案为前提进行实施,给出了详细的实施方式,本发明的保护范围包括但不限于两轴车辆。
配合图1所示,本发明的分布式电驱动车辆多性能优化力矩分配方法,采用基于改进的C/GMRES算法的MPC控制器进行多性能优化力矩分配,包括:MPC控制器以上层控制器给出的纵向总驱动力矩及横摆力矩作为力矩分配控制输入,结合车辆反馈的轮胎相关信息,基于车轮旋转动力学及轮胎滑移率状态更新模型,考虑轮毂电机输出力矩物理特性,以车轮负荷率为稳定性约束前提,以轮胎滑移耗散能最小化为优化目标,建立MPC优化控制问题,并采用改进的C/GMRES算法进行力矩优化分配问题的快速求解,以实现车辆稳定性与经济性的实时优化控制。
其中,上层控制器中,车辆的纵横向控制分别由PID(Proportion IntegrationDifferentiation,PID)和线性二次型规划器(Linear Quadratic Regulator,LQR)控制器实现。根据当前车速并给定车辆目标纵向车速,通过PID控制器输出总驱动力矩;根据质心侧偏角、横摆角速度和前轮转角作为LQR控制器的输入信号,输出外加横摆力矩,从而保证车辆横摆稳定性;所述整车控制方法包括但不限于PID控制方法、LQR控制方法、滑模变结构控制方法、模糊控制方法。根据所需的目标纵向总驱动力矩和目标横摆力矩,通过控制电动轮纵向力矩实现所述的整车目标控制力矩的分配。
图2为本发明实施例的车辆动力学关系图。下面建立二自由度车辆模型
系数矩阵为
其中xt=[β γ]T分别表示车辆质心侧偏角和横摆角速度;ut=ΔM为外部横摆力矩;w=δ表示车辆前轮转角;Cf和Cr分别表示轮胎侧偏刚度;la和lb分别表示前后轴到车辆质心的距离;m和vx分别表示车辆质量和纵向车速,Iz表示车辆横摆惯量。
稳态的质心侧偏角βr和横摆角速度γr可由下式表出:
因此可得到偏差的状态更新方程
设立性能指标Jlqr为
其中||·||表示欧几里得二范数,Qlqr和Rlqr分别为对角权重矩阵。通过求解上述方程的黎卡提方程可得出LQR控制器的控制增益向量Klqr,因而有ΔM=Klqre。需要指出,该LQR控制器控制参考选取为βr=0及
其中Kr表示欠转梯度;μ表示路面附着系数;g为重力加速度。
在进行下层MPC控制器设计之前,需要建立***状态更新模型。基于力矩平衡方程且不考虑制动力补偿的,车轮旋转动力学方程可表示为:
其中Js表示轮胎转动惯量;分别表示第i个车轮的角加速度,i=1,2,3,4分别表示左前、右前、左后和右后车轮,下同;Ti表示第i个轮毂电机的输出力矩;rw表示轮胎有效半径;Fxi表示第i个车轮的纵向力。
并定义轮胎纵向滑移率κi为:
其中,wi为第i个车轮的转速,vxi表示第i个车轮质心纵向速度。
其中Cxi为轮胎纵向滑移刚度,由上述方程及车轮旋转动力学方程联立,可通过上一时刻的Fxi和κi逆向估计轮胎纵向滑移刚度:
其中to表示控制器当前采样周期时间节点。上式中的值均是to-1时刻的,其中Ti(to-1)为考虑惯性环节后轮毂电机实际输出力矩,为轮毂电机反馈的旋转加速度,κi(to-1)为根据上式算得的上一时刻滑移率。由于车辆存在多个惯性环节,所以可认为这种简单的估计方法具有合适的精度。
综上,轮胎滑移率更新状态方程可表示为:
其中x=[κ1 κ2 κ3 κ4]T为状态量,u=[T1 T2 T3 T4]T为控制量。
本实施例通过基于改进的C/GMRES算法的MPC控制器实现分布式电动车轮胎力矩分配层轮胎滑移耗散能最小化。在本发明考虑的MPC控制问题中,瞬时代价方程可表示为:
l(x(τ),u(τ))=l1(x(τ),u(τ))+l2(x(τ),u(τ))
其中Ploss=[Fx1κ1vx1 Fx2κ2vx2 Fx3κ3vx3 Fx4κ4vx4]T表示纵向滑移功率的向量;
Q=diag(q1 q2 q3 q4)及R=diag(r1 r2 r3 r4)均表示对角权重矩阵。上式中第一项l1(x(τ),u(τ))的目的是令滑移功率的二次项最小化;l2(x(τ),u(τ))是为了限制轮毂电机输出力矩频繁抖动。
为了保证车辆横摆稳定及纵向驱动控制,控制问题的等式约束定义为
c(x(to),u(τ))=χu(τ)-υ
不等式约束定义为
h(x(to),u(τ))=[(u-umax)T (umin-u)T]T≤0
其中umax=[umax1 umax2 umax3 umax4]T,umaxi表示第i个车轮的最大输出力矩;
umin=[umin1 umin2 umin3 umin4]T,umini表示第i个车轮的最小输出力矩。定义Tmaxi与Tmini分别表示轮毂电机驱动力矩限制,由附着圆理论,轮胎负荷率需满足因而MPC控制器控制量约束定义为
其中η为某裕度系数。在控制中,控制器中使用的模型必定跟实际模型存在一定偏差,此处η取0.85,用以保证实际负荷率不超出1。
传统的C/GMRES算法无法考虑不等式约束,本发明对上述控制问题进行合理简化并有效解决了该不等式约束优化问题。假定等式约束成立,有
通过这一方式,将原方程预测时域内任一τ时刻的优化变量个数从4降为了2,且消去了原始问题的等式约束,简化了计算。
因而,本发明的MPC控制(求解)问题可表示为:
其中,Np表示预测时域长度,xo为当前时刻***反馈的状态量数值,ρ为阻碍方程的权重标量,为阻碍方程,其作用为阻止方程h1和h2到达0,即保证不等式约束成立。特别地,由凸优化理论可知,当ρ→0时,新的优化问题所得到最优解将无限接近于原优化问题的最优解。
根据上述控制问题构建哈密尔顿方程:
其中λ=[λ1 λ2 λ3 λ4]T为协状态向量,上式关于状态量x的偏导可表示为
定义预测时域内优化变量为
U=[T1(to) T2(to) … T1(τ) T2(τ) … T1(to+Np-1) T2(to+Np-1)]T
由于该MPC优化问题为零终端约束,且给定状态初始值xo,可以通过积分的方式数值计算相应的状态量与协状态,
其中
证明:将C/GMRES优化后的U(to)代入若最优解不满足约束,有h1→0,则同理,若h2→0,有此时与PMP最优条件矛盾,且C/GMRES优化的残差亦趋于∞。又由于GMRES算法理论表明:经无穷次迭代后总能收敛于最优解。综上,命题1成立。
由命题1和命题2可知,只要上一时刻控制量向量U(to-Δt)满足不等式约束,阻碍方程可有效抑制不等式约束令U(to)满足不等式约束,从而保证可逆。换句话说,只要U(0)满足不等式约束即可令控制器满足不等式约束与可逆条件。
因而,下面讨论所述改进的C/GMRES的初始化问题并提出一种计算迅速的初始化方法。为了简化初值优化过程,本实施例将算法初始化问题定义为Np=1的MPC控制量优化问题,获得后令其中表示维度为Np′2元素均为1的矩阵,表示克罗内克乘积,表示初始化的最优解。通过这种方法,初始化步骤得到了极大简化且满足命题1和命题2条件。基于Karush-Kuhn-Tucker(KKT)条件进行初值的优化,定义优化方程为
需要指出,为了确保命题1和命题2必定满足,此处对优化约束进行处理,引入某元素很小的列向量其内部元素均取10-10,目的为确保U(0)初始化的解满足不等式约束的。由于上述优化方程不存在等式约束,即不等式与等式约束线性无关,优化问题满足线性无关约束条件(Linearly Independent Constraint Qualification,LICQ),可用KKT条件进行计算。定义拉格朗日方程为
其中ε1为对应下界的拉格朗日乘子向量,ε2为对应上界的拉格朗日乘子向量。
求解步骤如下:
1)首先令ε1和ε2均为零向量,L可表示为,
由于上述优化问题中并不存在化简后同时包含T1和T2的项且是严格凸的,因而经过上述2步计算后的解即为该问题的最优解,也易证此满足KKT最优性条件。以上改进的C/GMRES算法的应用过程可按图3所示流程逐步进行。具体为:判断当前采样时刻是否为改进的C/GMRES算法的第一个采样时刻,若是,则进行“初始化”计算,如不是,则进行“进行期间”的常规计算。
初始化计算包括:
进行期间的常规计算包括:
S201、基于U(t0-Δt),x(t0)以及车辆反馈信息,进行滑移率κi的向前推进计算;
S202、根据κi,采用反向递归的方式计算预测时域内协状态量[λ1λ2λ3λ4]T;
S203、计算最优条件矩阵函数F(U(to),xo(t));
S204、通过有限差分法及GMRES算法,计算最优解U(t0),将U(t0)=U(t0-Δt)+U(t0)Δt设置为力矩分配命令。
综上所述,本发明采用基于改进的C/GMRES算法的MPC控制器进行多性能优化力矩分配,在保证控制性能的基础上,该算法可计算效率高,适用于实时应用;针对传统C/GMRES算法无法求解不等式约束问题,本发明将不等式约束通过阻碍方程等效转换为优化目标代价,实现不等式约束问题求解;针对C/GMRES算法求解,提出了一种基于KKT条件的初始化计算方法,该方法可快速获得C/GMRES算法初值,可提升C/GMRES的适用性。
以上所述,仅是本发明的较佳实施例而已,并非对本发明作任何形式上的限制,本领域技术人员利用上述揭示的技术内容做出些许简单修改、等同变化或修饰,均落在本发明的保护范围内。
Claims (3)
1.一种分布式电驱动车辆多性能优化力矩分配方法,其特征在于,采用基于改进的C/GMRES算法的MPC控制器进行多性能优化力矩分配,包括:
MPC控制器以上层控制器给出的纵向总驱动力矩及横摆力矩作为力矩分配控制输入,结合车辆反馈的轮胎相关信息,基于车轮旋转动力学及轮胎滑移率状态更新模型,考虑轮毂电机输出力矩物理特性,以车轮负荷率为稳定性约束前提,以轮胎滑移耗散能最小化为优化目标,建立MPC优化控制问题,并采用改进的C/GMRES算法进行力矩优化分配问题的求解,其中:
所述改进的C/GMRES算法中,通过阻碍方程的形式将不等式约束等效转换为目标方程代价,实现不等式约束问题求解;所述改进的C/GMRES算法中,采用基于KKT条件的初始化计算方法求解;
所述MPC控制器的控制问题设定过程中将分布式电驱动车辆的前后轴轮毂电机输出力矩约束进行合并;
所述轮胎滑移率状态更新模型中的轮胎滑移率更新状态方程表示为:
其中x=[κ1 κ2 κ3 κ4]T为状态量,u=[T1 T2 T3 T4]T为控制量;κi为第i个轮胎的滑移率,Ti为第i个轮毂电机的输出力矩,i=1,2,3,4分别表示前左、前右、后左和后右;为状态量x关于时间变量t的微分向量;
所述MPC控制器的控制问题表示为:
xo=x(to)
其中,Np表示预测时域长度;为瞬时代价方程;xo为当前时刻***反馈的状态量数值;ρ为阻碍方程的权重标量;为阻碍方程;x=[κ1 κ2 κ3κ4]T为状态量,u=[T1 T2 T3 T4]T为控制量;κi为第i个轮胎的滑移率,Ti为第i个轮毂电机的输出力矩,i=1,2,3,4分别表示前左、前右、后左和后右;to表示当前采样时刻;为***状态方程;Jmpc为MPC控制问题代价,τ为MPC控制问题积分中的微分节点,dτ为MPC控制问题积分中的微分变量,h1和h2分别是由不等式约束方程变换而来的代价值;
基于所述MPC控制器的控制问题构建哈密尔顿方程,计算预测时域内状态量更新及协状态量;根据庞德里亚金极小值原理最优性条件写出矩阵函数F(U(to),xo(t)),其中xo为当前时刻***反馈的状态量数值,U(to)为当前时刻MPC控制器待优化向量,to表示当前采样时刻,t为时间变量;求解计算获得最优解。
2.根据权利要求1所述的分布式电驱动车辆多性能优化力矩分配方法,其特征在于,所述矩阵函数表示为:
其中,表示关于目标方程及***状态方程写出的哈密尔顿方程;Np为预测时域长度;x=[κ1 κ2 κ3 κ4]T为状态量,κi为第i个轮胎的滑移率,Ti为第i个轮毂电机的输出力矩,i=1,2,3,4分别表示前左、前右、后左和后右;λ为协状态向量;xo为当前时刻***反馈的状态量数值;U(to)为当前时刻MPC控制器待优化向量;to表示当前采样时刻;基于xo及U,以有限差分的形式计算即
其中,ψs表示稳态因子,U表示预测时域内优化变量。
3.根据权利要求1所述的分布式电驱动车辆多性能优化力矩分配方法,其特征在于,所述改进的C/GMRES算法求解的初始化计算方法为:
S101、基于车辆反馈的轮胎相关信息,通过获得最优解其中Q=diag(q1 q2 q3 q4)及R=diag(r1r2 r3 r4)均表示对角权重矩阵,Js为轮胎转动惯量,Ttot为车辆总驱动力矩需求,ΔM为车辆外加横摆力矩需求,rw为车轮有效半径,ds为车辆轮间距,为第i个车轮的旋转加速度,i=1,2,3,4分别表示前左、前右、后左和后右;
S102、对于中的元素,若第i个元素大于则令该元素等于若第i个元素小于则令该元素等于其中,ξ为某元素很小的列向量,其内部元素均取10-10;和分别为第i个轮毂电机的最大及最小约束,i=1,2,3,4分别表示前左、前右、后左和后右;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910530350.1A CN110239362B (zh) | 2019-06-19 | 2019-06-19 | 一种分布式电驱动车辆多性能优化力矩分配方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910530350.1A CN110239362B (zh) | 2019-06-19 | 2019-06-19 | 一种分布式电驱动车辆多性能优化力矩分配方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110239362A CN110239362A (zh) | 2019-09-17 |
CN110239362B true CN110239362B (zh) | 2020-10-02 |
Family
ID=67887945
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910530350.1A Active CN110239362B (zh) | 2019-06-19 | 2019-06-19 | 一种分布式电驱动车辆多性能优化力矩分配方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110239362B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111123860B (zh) * | 2019-12-11 | 2020-11-10 | 山东科技大学 | 带有非对称信息的网络***的最优分布式控制方法 |
CN111267834B (zh) * | 2020-02-17 | 2021-03-26 | 北京理工大学 | 一种车辆横摆稳定预测控制方法及*** |
CN111619548B (zh) * | 2020-05-28 | 2021-01-12 | 清华大学 | 车辆驱动防滑控制方法、装置、计算机设备和存储介质 |
CN111703429B (zh) * | 2020-05-29 | 2022-05-10 | 北京理工大学重庆创新中心 | 一种轮毂电机驱动车辆纵向速度估算方法 |
CN112793560B (zh) * | 2020-12-31 | 2022-08-19 | 北京理工大学 | 基于转矩矢量控制的无人驾驶车安全性与操稳性控制方法 |
CN113459829B (zh) * | 2021-07-29 | 2023-02-21 | 南昌智能新能源汽车研究院 | 一种基于路况预测的双电机电动汽车智能能量管理方法 |
CN117125079B (zh) * | 2023-08-28 | 2024-07-05 | 广州汽车集团股份有限公司 | 车辆的控制方法及装置、电子设备、存储介质 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104773169B (zh) * | 2015-04-28 | 2017-05-24 | 吉林大学 | 一种基于轮胎侧偏角的车辆横摆稳定集成控制方法 |
US10209685B2 (en) * | 2015-10-29 | 2019-02-19 | Mitsubishi Electric Research Laboratories, Inc. | Method and apparatus for preconditioned model predictive control |
CN105667343A (zh) * | 2016-03-16 | 2016-06-15 | 吉林大学 | 电动汽车牵引力控制方法 |
CN108556680B (zh) * | 2018-01-17 | 2020-02-14 | 北京理工大学 | 一种用于轮毂电机驱动车辆的稳定性控制方法 |
CN109291932B (zh) * | 2018-10-16 | 2021-03-12 | 浙江东方机电有限公司 | 基于反馈的电动汽车横摆稳定性实时控制装置及方法 |
-
2019
- 2019-06-19 CN CN201910530350.1A patent/CN110239362B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110239362A (zh) | 2019-09-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110239362B (zh) | 一种分布式电驱动车辆多性能优化力矩分配方法 | |
CN109263716B (zh) | 一种四轮毂电机驱动车辆转向的控制方法 | |
CN109849899B (zh) | 一种电动轮汽车电液复合车身稳定控制***及方法 | |
CN111267834B (zh) | 一种车辆横摆稳定预测控制方法及*** | |
Goodarzi et al. | Design of a VDC system for all-wheel independent drive vehicles | |
CN108556680B (zh) | 一种用于轮毂电机驱动车辆的稳定性控制方法 | |
CN112829766B (zh) | 一种基于分布式驱动电动车辆的自适应路径跟踪方法 | |
CN111959486B (zh) | 电机驱动车辆纵横向耦合控制方法、***及存储介质 | |
CN110422052B (zh) | 分布式驱动电动汽车稳定与节能控制*** | |
CN110481343A (zh) | 四轮轮毂电机驱动汽车力矩补偿的组合二阶滑模控制方法 | |
CN112519882A (zh) | 一种车辆参考轨迹跟踪方法及*** | |
CN111158264B (zh) | 面向车载应用的模型预测控制快速求解方法 | |
CN110605973A (zh) | 一种基于分层结构的多轴分布式电驱动车辆操纵稳定性控制方法 | |
CN111142534B (zh) | 一种智能车横纵向综合轨迹跟踪方法及控制*** | |
CN112622875B (zh) | 一种四轮毂电机驱动车辆的下层力矩分配控制方法 | |
CN116872910A (zh) | 无人驾驶车辆轨迹跟踪与横摆稳定性控制方法及*** | |
CN113442906A (zh) | 一种分布式驱动电动汽车横向稳定性分层控制***及方法 | |
CN114148403B (zh) | 一种线控转向***多工况稳定性控制方法 | |
CN109017447B (zh) | 无人驾驶车辆总的驱动力矩输出方法 | |
CN114670828A (zh) | 一种采用复合模型预测控制器的自适应巡航控制方法 | |
Raji et al. | A tricycle model to accurately control an autonomous racecar with locked differential | |
CN115723590A (zh) | 一种轮毂电机驱动汽车的节能转矩矢量控制方法 | |
WO2022106005A1 (en) | Methods for reducing high-speed off-tracking in multi-trailer heavy duty vehicles | |
CN111559361A (zh) | 一种四轮独立驱动车辆转向制动协同控制方法 | |
CN116819972B (zh) | 一种模块化分层式架构的协同控制方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |