CN112622875B - 一种四轮毂电机驱动车辆的下层力矩分配控制方法 - Google Patents
一种四轮毂电机驱动车辆的下层力矩分配控制方法 Download PDFInfo
- Publication number
- CN112622875B CN112622875B CN202011218555.5A CN202011218555A CN112622875B CN 112622875 B CN112622875 B CN 112622875B CN 202011218555 A CN202011218555 A CN 202011218555A CN 112622875 B CN112622875 B CN 112622875B
- Authority
- CN
- China
- Prior art keywords
- control
- vehicle
- wheel
- torque
- tire
- 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 63
- 238000012546 transfer Methods 0.000 claims abstract description 10
- 238000011217 control strategy Methods 0.000 claims description 23
- 230000008569 process Effects 0.000 claims description 11
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000005096 rolling process Methods 0.000 claims description 9
- 238000012937 correction Methods 0.000 claims description 6
- 230000001133 acceleration Effects 0.000 claims description 5
- 238000005457 optimization Methods 0.000 claims description 5
- 238000013459 approach Methods 0.000 claims description 2
- 230000009977 dual effect Effects 0.000 claims description 2
- 239000000126 substance Substances 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 claims 1
- 150000001875 compounds Chemical class 0.000 claims 1
- 238000004364 calculation method Methods 0.000 abstract description 9
- 230000000694 effects Effects 0.000 description 11
- 238000010586 diagram Methods 0.000 description 8
- 230000008859 change Effects 0.000 description 7
- 238000004088 simulation Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 230000002411 adverse Effects 0.000 description 2
- 210000001503 joint Anatomy 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 230000001687 destabilization Effects 0.000 description 1
- 230000000368 destabilizing effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000036039 immunity Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000013598 vector Substances 0.000 description 1
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
- B60W30/00—Purposes of road vehicle drive control systems not related to the control of a particular sub-unit, e.g. of systems using conjoint control of vehicle sub-units
- B60W30/02—Control of vehicle driving stability
-
- 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/32—Control or regulation of multiple-unit electrically-propelled vehicles
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
- B60W50/00—Details of control systems for road vehicle drive control not related to the control of a particular sub-unit, e.g. process diagnostic or vehicle driver interfaces
-
- 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
- B60L2240/00—Control parameters of input or output; Target parameters
- B60L2240/40—Drive Train control parameters
- B60L2240/42—Drive Train control parameters related to electric machines
- B60L2240/423—Torque
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
- B60W50/00—Details of control systems for road vehicle drive control not related to the control of a particular sub-unit, e.g. process diagnostic or vehicle driver interfaces
- B60W2050/0001—Details of the control system
- B60W2050/0019—Control system elements or transfer functions
- B60W2050/0028—Mathematical models, e.g. for simulation
- B60W2050/0031—Mathematical model of the vehicle
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
- B60W2720/00—Output or target parameters relating to overall vehicle dynamics
- B60W2720/40—Torque distribution
Landscapes
- Engineering & Computer Science (AREA)
- Transportation (AREA)
- Mechanical Engineering (AREA)
- Automation & Control Theory (AREA)
- Human Computer Interaction (AREA)
- Power Engineering (AREA)
- Control Of Driving Devices And Active Controlling Of Vehicle (AREA)
- Electric Propulsion And Braking For Vehicles (AREA)
Abstract
一种四轮毂电机驱动车辆(4MIDEV)的下层力矩分配控制方法,相对于现有技术,提出了一种下层双级最优转矩分配控制方法。上层控制层获得了目标控制力及力矩后,通过下层双级最优转矩分配控制器将其以转矩指令形式分配到各驱动轮。第一级分配控制计算成本低、实时性较好,可在良好路面附着条件下实现车辆最大的稳定性裕度;第二级分配控制引入权重系数考虑各轮不同的附着情况,以保证车辆在附着条件较差尤其是道路附着系数发生改变时的操纵稳定性。最后,基于轮胎力最优解,采用经典滑模控制方法根据各轮滑转/移率调节电动轮输出转矩,完成最终的转矩分配,实现对车辆的操纵稳定性控制,有效改善了4MIDEV面对不同道路时的操纵稳定性与实时性。
Description
技术领域
本发明涉及分布式驱动车辆控制领域,具体涉及一种四轮毂电机驱动车辆(4MIDEV) 的下层力矩分配控制方法,用于实现操纵稳定性控制。
背景技术
四轮毂电机驱动车辆四个驱动轮的驱动/制动转矩可以独立控制,其驱动冗余相对于现有技术更有利于实现转向过程中的操纵稳定性控制。车辆的操纵稳定性对于行驶机动性及安全性有重要意义,是车辆动力学控制中重点研究内容之一。4MIDEV的操纵性稳定性控制策略通常被设计成分层控制结构,上层运动控制器根据驾驶员模型指令产生目标控制力及力矩,下层力矩分配控制器根据目标控制力及力矩产生各驱动轮转矩指令。然而,现有的4MIDEV操纵稳定性下层力矩分配控制还主要存在以下几个问题:
1、在目前研究中,转矩分配控制往往忽略了由地面附着不同带来的各轮附着条件不一致的问题,在转矩分配过程中,每个车轮的轮胎附着利用率通常被认为是同等重要,即它们的权重系数是相等的,这就有可能导致某个车轮超出附着极限,而其他车轮仍保持在较高附着裕度的情况。2、最优分配的约束,除来自上级控制层的广义力(矩)的约束外,还有摩擦椭圆约束、电机物理约束、滑移率约束等。其中,摩擦椭圆约束与轮胎纵向力和侧向力以及车轮处路面附着系数相关,是典型的非线性不等式约束,会加大最优问题求解的难度和时间,对控制的实时性造成不利影响。
发明内容
为克服上述现有技术中所存在的技术问题,本发明提供了一种四轮毂电机驱动车辆的下层力矩分配控制方法,具体包括以下步骤:
首先,构建双级控制结构,首先根据上层运动控制器输出的目标控制力及力矩,采用第一级分配控制求得轮胎纵向力最优解判断是否满足地面附着约束;若满足,则输出若不满足,则采用第二级分配控制,引入权重系数考虑各轮处不同的附着情况,求得轮胎纵向力最优解并输出;
其次,构建经典滑模控制,根据各轮滑转/移率,以及上一步输出的轮胎纵向力最优解,求得电动轮输出转矩,完成转矩的分配。
本发明还提供了一种四轮毂电机驱动车辆的综合控制方法,具体包括以下步骤:
步骤1、针对四轮毂电机驱动车辆,建立七自由度车辆动力学模型;
步骤2、建立单点预瞄驾驶员模型;
步骤3、建立具有多层级结构的操纵稳定性双级最优转矩分配控制策略;
其中,步骤3具体包括:
步骤3.1、建立运动参考层,获得期望质心侧偏角与期望横摆角速度;
步骤3.2、设计上层运动控制器,根据横摆角速度期望值与纵向速度期望值获得期望纵向力与期望横摆力矩;
步骤3.3、设计下层双级力矩分配控制器,求解轮胎纵向力的最优解,基于滑模控制方法设计滑转/移率控制器,根据各轮滑转/移率与轮胎纵向力最优解调节电动轮输出转矩,完成执行器对轮胎力的实现,从而达成最终的转矩分配。
进一步地,综合控制方法中所述步骤1中所建立的七自由度动力学模型,采用以下形式:
其中,矩阵Bx、By分别为:
式中,Vx、Vy和γ分别表示车辆的纵向速度、侧向速度及横摆角速度,和则分别是它们的导数。Fxij与Fyij代表各轮胎所受的纵向力和侧向力,其中i∈{f,r}表示前后,j∈{l,r}表示左右。δf为前轮转角。m为车辆质量,a与b分别表示车辆的前轴距和后轴距,d为轮距的一半,Iz表示横摆转动惯量。
车轮旋转运动作为电动轮特性分析中重要部分,需要对其进行准确分析。忽略轮胎侧向力的影响,对单个车轮进行受力分析,其旋转运动方程可以描述为:
其中,Iwij表示车轮旋转惯量,Rij为车轮滚动半径,为车轮旋转角加速度,Tij为轮毂电机驱动/制动转矩,Tbij为机械制动转矩。Tfij为滚动阻力矩,与垂向力Fzij和轮胎拖矩Δij有关,Tfij=FzijΔij=FzijfRij,其中f为路面摩擦系数。
车轮的所受垂向力则可以描述为:
在车辆实际行驶过程中,车轮并不常处于纯滚动状态,其轮心速度和滚动速度不相等时,就将发生滑转/滑移。滑转/移率λij用以描述车轮滑转/滑移程度,定义为:
其中,vwij代表各轮纵向轮心速度,Rij为各轮有效滚动半径,ωij为车轮转动角速度。
进一步地,所述步骤2所述驾驶员模型,采用如下形式:
xOy与XOY分别是车辆坐标系和大地坐标系,车辆任意时刻位置在上述坐标系中的坐标(x,y)与(X,Y)存在以下关系:
在实际行车过程中,驾驶员通常选择前方道路一点(或一段)作为参照。假定驾驶员预瞄时间为T,则t时刻期望轨迹上的参照点为A(x(t+T),f(x(t+T))),f(x)为关于x 的函数,表示期望轨迹参照点的纵坐标。而经过T时间后,车辆将运动至 B(x(t+T),y(t+T))。考虑到预瞄时间T一般较小,在预瞄过程中可以将车速vx简化为恒定值,于是有:
进一步地,所述步骤3中建立的具有多层级结构的操纵稳定性双级最优转矩分配控制策略具有如下形式:
基于DYC方法,发明了一种具有多层级结构的操纵稳定性双级最优转矩分配控制策略,对下层转矩分配控制器设置双级分配策略。其主要包括运动参考层、上层运动控制器及下层转矩分配控制器;
上层运动控制器由速度跟随控制和横摆力矩控制组成,根据运动参考层生成的期望值,生成纵向牵引力和横摆力矩虚拟控制指令,可根据路面附着自适应改变控制方法,分别实现对车辆速度的控制以及横摆角速度和质心侧偏角的综合控制。
下层转矩分配控制器基于最优控制方法,设计为双级分配控制结构,通过分配控制级的自适应切换及对车轮滑转/移率的调节,实现了对各驱动轮转矩的双级最优转矩分配控制,在满足车辆期望纵向力及横摆力矩的同时,保证了车辆在不同附着条件下最大的稳定性裕度并减小了求解计算成本。
进一步地,所述步骤3.1中建立的运动参考层,具有以下形式:
为更加直接和清楚地反映驾驶员的转向意图,运动参考层将上述车辆动力学模型简化为二自由度作为参考模型。通过对二自由度模型的分析,可以获得质心侧偏角和横摆角速度的期望值:
考虑到所得的横摆角速度与质心侧偏角期望值,出现在轮胎的线性工作区,如果路面附着系数较低,轮胎将进入非线性工作区,此时必须对横摆角速度和质心侧偏角的期望值上限进行限制,即
其中,μ为附着系数。
综上,通过运动参考层可以获得横摆角速度和质心侧偏角的期望值:
进一步地,所述步骤3.2中建立的上层运动控制器,具有以下形式:
控制目标纵向力可以实现车辆纵向速度对期望纵向车速的跟随控制。以车辆纵向速度为控制变量,根据实际纵向车速与人为设定的期望车速之间的偏差,利用PID控制方法控制车辆期望纵向驱动力Fx_des,以满足车辆对期望车速跟随。实际纵向车速可以通过驾驶员对加速踏板和制动踏板的输入信号估计得出。
横摆力矩控制选取质心侧偏角β和横摆角速度γ作为控制变量,减小其实际值与运动参考层输出期望值的偏差,考虑到质心侧偏角与横摆角速度相互相耦合,并且质心侧偏角不便直接测量且估计精度有限,本文通过对质心侧偏角的检测利用模糊控制方法设计修正系数K,以修正名义横摆角速度,即γ0=Kγ。又考虑到质心侧偏角变化率的影响,通过β-Method理论和相平面法,设计综合控制量作为模糊控制器的输入。其中,B1是与附着系数有关的系数。此外,由于路面附着系数不同时,质心侧偏角对稳定性的影响不同,将附着系数μ也引入修正系数的设计中,将其作为模糊控制器的另一个输入。表1示出了修正系数K的模糊规则,其中负大(NB),负中(NM),负小(NS),零(Z),正小(PS),正中(PM),正大(PB)为定义的模糊集。
表1修正系数K模糊规则
为适应系统的非线性并尽可能减少计算成本,本文采用模糊PID控制方法对修正名义横摆角速度γ0进行控制。模糊PID控制器的输入分别是修正名义横摆角速度γ0和横摆角速度期望值γ_des的偏差,以及将偏差进行微分处理得到的偏差率,输出是横摆力矩期望值Mz_des。表2示出了模糊PID控制器模糊规则,其中Δkp、Δki、Δkd分别是比例、积分、微分环节系数修正量。
表2模糊PID控制器模糊规则
进一步地,所述步骤3.3中建立的下层双级力矩分配控制器,具有以下形式:
上层与下层制层之间的关系可以描述为:
v=Bu
其中,v表示所需的目标控制力及力矩Fd,v=[Fx_des Mz_des]T;u代表各轮胎纵向力,u=[Fxfl Fxfr Fxrl Fxrr]T;B为系数矩阵,可描述为:
由于采样周期一般较短,可以忽略转向角δf在此期间的变化,系数矩阵B可视为常数矩阵。
考虑到来自地面附着的限制,轮胎纵向力还需要满足摩擦圆约束:
式中μij及Fzij分别表示各轮胎处路面附着系数及所受垂向载荷。
v和u分别是一个2维和4维的向量,且u中各项相互独立。则对于同一目标v,往往存在多组满足约束条件的可行解u,这也反映了目标车辆是一个过驱动系统。不同的u 代表不同可行的轮胎纵向力分配方式,由此影响的系统稳定性裕度、轮胎磨损程度、整车能量利用率等也会不一样。针对上述情况,为获取最佳的轮胎纵向力分配u*,可构造相关目标函数J(u),求解摩擦圆约束下该目标J(u)的最优解。于是,轮胎力对目标控制力及力矩的实现,可归纳为一个带有约束的优化问题。
最优轮胎纵向力u*在执行器上的实现可通过车轮转旋转运动方程分析,若忽略滚动阻力的影响,执行器与轮胎力纵向力之间的关系可以描述为:
其中Twij为电动轮的输出转矩,Twij∈[-Tijmax-Tbijmax,Tijmax]。在这主要考虑轮毂电机驱动转矩Tij的影响,机械制动转矩Tbij主要在制动工况时对其起补充作用,不展开讨论。可见,当车轮滑转/移率较小时,可以忽略轮速在采样时间内的变化,电动轮输出转矩与轮胎力存在清晰关系:
Twij=FxijRij
而当滑转/移率较大时,为保证车辆操纵稳定性则需要根据滑转/移率对电动轮输出转矩进行调节。
因此,下层转矩分配的两个子问题,可以归结为一个带有约束条件的优化问题及一个滑转/移率控制问题。
考虑到路面附着的影响及算法实现的实时性,轮胎力对Fx_des的实现基于最优分配控制算法,设计为双级控制结构,其中第一级分配控制计算成本低、实时性较好,可在良好路面附着条件下实现车辆最大的稳定性裕度;第二级分配控制引入权重系数考虑各轮处不同的附着情况,在保证车辆较大稳定性裕度的同时确保低附着处的车轮具有较大的附着裕度,以改善车辆在附着条件较差尤其是道路附着系数发生改变时的操纵稳定性。下面分别对其层级切换准则及各级控制层进行阐述。
轮胎纵向力需要满足摩擦圆的约束。摩擦圆约束与轮胎纵向力和侧向力相关,是典型的非线性不等式约束,并且由魔术公式计算轮胎侧向力,过程复杂、计算成本高,若将摩擦圆约束直接考虑到最优问题求解当中,会大大加大求解的难度和时间,对控制的实时性造成不利影响。此外,车辆失稳工况往往出现在轮胎侧向力和摩擦圆的非线性区域。
于是,本文将摩擦圆约束作为层级切换准则,不使其直接参与最优问题的求解,而是用以判断实时性较好的第一级分配控制得到的轮胎力最优解是否满足地面附着约束:若满足,则控制电动轮输出转矩实现完成转矩分配;若不满足,则继续第二级分配控制,求得轮胎力最优解再进行转矩控制,以在满足上层目标控制力及力矩的同时,保证车辆在不同附着条件下较大的稳定性裕度。
考虑到车辆操纵稳定性控制是由当前时刻的车辆状态,施加对下一时刻的预控制,也即当前求得的轮胎力分配方式是下一个时刻所需实现的。为区别表示,设当前时刻为t,则下一时刻需实现的轮胎纵向力为Fxij(t+1),摩擦圆约束可改写为:
简化轮胎模型可在保证一定精度的同时,简化对轮胎侧向力的估算过程,便于检验分配算法中轮胎力是否超出附着极限。其中,估算的侧向力可描述为:
式中,αij是各个轮胎的侧偏角,Kαij为各轮侧偏刚度,其中i∈{f,r},j∈{l,r},分别代表前后、左右。
由上式可以看出,ρij越大,轮胎纵向力越接近地面附着极限,一方面降低了轮胎纵向力的输出储备,影响车辆动力性;另一方面由于轮胎侧向力在摩擦圆中所占比例相应减小,易导致车辆侧滑,降低了车辆的稳定性裕度。尤其当ρij=1时,轮胎侧向力为0,已经无法保证车辆转向时的操纵稳定性。相反地,减小ρij可以获得更大的车辆稳定性裕度及轮胎纵向力储备。因此,选取各轮轮胎附着利用率平方和构建优化目标函数如下:
此外,轮胎纵向力Fxij需满足七自由度车辆动力学模型,以实现目标控制力及力矩Fd。若作小角度假设,可简化为:
为减小计算成本、提高求解效率,结合轮胎附着利用率消去任意两个变量,可得:
综上,第一级最优分配控制方程为:
第二级分配控制中同样采用轮胎附着利用率ρij描述车辆的稳定性裕度,但考虑到车辆在实际行驶过程中,各车轮处的附着条件不总是一样的,如果将轮胎附着利用率ρij简单相加作为性能目标函数,可能出现各轮附着裕度差别较大的情况。此时,若遇见道路附着系数发生改变等特殊工况,则可能出现处于低附着条件的车轮失稳,而处于高附着条件的车轮仍有较大稳定裕度的情况。
于是,第二级分配控制在第一级分配控制的性能目标函数J1的基础上,引入权重系数ψij考虑各轮附着情况,构建新的目标函数:
J′2=||Γ1u||2
定义ψij为:
通过引入权重系数ψij可以使得在求解车辆最大稳定性裕度的同时,低附着条件下的车轮具有较大的附着裕度,以避免道路附着系数改变时某个车轮超出附着极限。
此外,考虑到上下级动力学关系作为等式约束时,其约束强度过强,将Bu=v转换为||Wv(Bu-v)||2作为逼近误差惩罚项,则有:
J2=||Γ1u||2+ξ||Wv(Bu-v)||2
式中,ξ为惩罚项的权重系数,其值越大,误差逼近项的控制优先级就越高。 Wv=diag(WFxd,WMzd),表示权重矩阵,其与性能目标函数中牵引力和横摆力矩的控制优先级有关。对于低附着路面,可以增大横摆力矩的控制权重WMzd,减小牵引力的控制权重WFxd,以保证车辆良好操纵稳定性。
综上,第二级最优分配控制方程为:
将上式进一步整理,有
但当车轮滑转/移率较大时,为保证车辆操纵稳定性,需要根据滑转/移率对电动轮输出转矩进行调节,以实现期望的轮胎纵向力。
在不同附着系数的道路上,车轮最佳滑转/移率是不一样的。道路附着系数越小,最佳滑转/移率越小。若通过附着系数辨识,由“滑转/滑移率-附着系数”曲线确定最佳滑转/ 移率,再进行转矩控制,其过程复杂、计算成本过高,不利于控制的实时性。考虑到当轮胎附着利用率较小时,轮胎纵向力较小,一定程度上避免了车轮打滑。在此,针对各路面统一取车轮最佳滑转/移率λ*为0.2,以便于后续控制。
vwij代表各轮纵向轮心速度,Rij为各轮有效滚动半径。
采用经典滑模控制对滑移率控制器进行设计,定义滑模面如下:
取带饱和函数的等速趋近率,有:
可得控制律Twij为:
综合滑转/移率较小的情况,最终可得电动轮输出转矩Twij为:
其中Twij∈[-Tijmax-Tbijmax,Tijmax]。
上述本发明所提供的方法,相对于现有技术至少具有以下有益效果:
(1).本发明设计了下层控制器轮胎力双级分配控制结构,以获取轮胎力的最优解。其中第一级分配控制计算成本低、实时性较好,可在良好路面附着条件下实现车辆最大的稳定性裕度。
(2).第二级分配控制引入权重系数考虑各轮不同的附着情况,以保证车辆在附着条件较差尤其是道路附着系数发生改变时的操纵稳定性。
(3)本发明采用滑模控制方法根据各轮滑移/转率调节电动轮输出转矩,以完成最终的转矩分配,考虑到了不同附着系数对滑转/移率的影响,在满足车辆目标纵向力及附加横摆力矩的同时,保证了车辆在不同附着条件下较大的稳定性裕度。
附图说明
图1为本发明所提供方法的执行过程框图;
图2为四轮毂电机驱动车辆七自由度动力学模型示意图;
图3为车轮旋转运动模型;
图4为驾驶员预瞄模型示意图;
图5为四轮毂电机驱动车辆二自由度动力学模型示意图;
图6为上层运动模糊PID控制器示意图;
图7为下层转矩分配控制器的控制逻辑框图;
图8为轮胎力分配的双级控制结构示意图;
图9为基于Simulink-Carsim仿真实例的控制系统模型;
图10为Carsim对接路面设置图;
图11为对接路面下双移线工况仿真结果;
图12为Carsim对开路面设置图;
图13为对开路面下双移线工况仿真结果;
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明所提供的四轮毂电机驱动车辆下层力矩分配控制方法,在图1中示出了其如何用于实现此种车辆的操纵稳定性力矩分配控制全过程。该方法具体包括以下步骤:
步骤1、针对四轮毂电机驱动车辆,建立七自由度车辆动力学模型;
步骤2、建立单点预瞄驾驶员模型;
步骤3、建立具有多层级结构的操纵稳定性双级最优转矩分配控制策略;
其中,步骤三具体包括:
步骤3.1、建立运动参考层,获得期望质心侧偏角与期望横摆力矩;
步骤3.2、设计上层运动控制器,根据横摆角速度期望值与纵向速度期望值获得期望纵向力与期望横摆力矩;
步骤3.3、设计下层双级力矩分配控制器,求解轮胎力的最优解,基于滑模控制方法设计滑转/移率控制器,根据各轮滑转/移率与轮胎力最优解调节电动轮输出转矩,完成执行器对轮胎力的实现,从而达成最终的转矩分配。
图9-图13示出了基于本发明所提供方法的一个具体实例,采用Matlab/Simulink仿真实验,表1示出了所设置车辆及电机参数。
为整体评估操纵稳定性双级最优转矩分配控制策略的有效性,本文在不同车速在对接道路、对开双移线道路上进行驾驶员闭环工况仿真,分析提出的双级最优转矩分配控制策略(简称“双级控制”)相较于其他稳定性控制策略的控制效果及可靠性。对比策略包括:1)简称为“速度控制”策略,其上级运动控制层仅控制纵向车速而不对侧向运动进行控制,下级转矩控制层采用转矩平均分配方法且不对滑移率进行控制;2)简记为“一般控制”策略,其上级运动控制层采用经典滑模控制方法实现对横摆角速度、质心侧偏角及车速的控制,下层转矩分配则仅采用文中第一级分配控制结构。
图10示出了Carsim中对接路面设置。
图11(a)描述了在不同操纵稳定性控制策略下对接路面上车辆的行驶轨迹。可以看出,各控制策略在附着系数改变前均可以保证车辆对期望轨迹良好的跟踪效果,而当附着系数改变后(距起始位置105m处,约8s时),车辆行驶轨迹相比于期望轨迹出现了不同程度的偏移:一般控制策略下车辆轨迹侧向偏移量最大,但其可在一定时间内逐渐抑制由附着条件变化带来的扰动,保证车辆的侧向稳定;速度控制下的车辆开始出现逐渐发散的左右摆动行驶,趋于失稳状态;相比之下,双级控制可针对道路附着条件自适应调节,能在较短时间内抑制扰动,车辆的侧向偏移量最小,呈现最好的轨迹保持控制效果。
图11(b)-(d)分别展示了各控制策略下对接路面上车辆的速度、横摆角速度及质心侧偏角的响应。可以看到,当路面附着系数发生改变后,各运动状态量均出现了程度不一的扰动。其中,双级控制采用改进滑模控制方法较快地抑制了扰动,并在参考值动态调节的基础上实现了各量对其期望值良好的跟踪控制;而一般控制下各量与期望值的偏差均较大,抑制扰动时间略长;速度控制下各量出现发散的情况,运动控制效果最差。由图11(e)所示的车辆质心侧偏角及其变化率的相平面可见,双级控制与速度控制下的各相平面曲线分别离原点最近与最远,这同样表明双级控制策略的操纵稳定控制效果最好,而速度控制策略控制效果最差,一般控制策略介于二者之间。
图11(f)-(g)分别表示各控制策略下对接路面上电动轮的输出转矩。速度控制输出转矩整体较小,但因无法形成附加横摆力矩而过度依赖驾驶员转向操作,在工况后期出现逐渐发散的波动。双级控制相较于一般控制策略可自适应调节运动控制参考值及转矩分配权重,在附着系数变化时侧重于对质心侧偏角的抑制,并降低了前轴的转矩输出,其输出转矩的最大值及平均值均较小。
图12示出了Carsim中对开路面设置。
图13(a)及图13(b)分别描述了在不同操纵稳定性控制策略下对开路面上车辆的行驶轨迹及速度响应。由于一侧道路附着条件良好,驾驶员仍可以通过修正方向盘转角实现对车辆行驶方向的控制,此时三种控制策略均能较好地保证车辆对期望轨迹与速度的跟随,其中双级控制下车辆侧向偏移量及速度波动最小,如Carsim实际仿真图13(a) 及图13(b)所示,呈现了最好的控制效果。
图13(c)与13(d)分别展示了对开路面上车辆的横摆角速度及质心侧偏角的响应,图13(e)为车辆质心侧偏角及其变化率的相平面图。可以看到,当车辆进入对开路面后,各控制策略下横摆角速度及质心侧偏角均出现了程度不一的扰动,而双级控制下二者受扰动最小且最快得到了抑制,表现出较强的抗扰性。另一方面,双级控制下的相平面曲线离原点最近,同样表明了其操纵稳定控制效果最好;相比之下,一般控制策略的控制效果次之,速度控制的效果最次。
从图13(f)-(g)可以看出,在附着系数变化后,一般控制下输出转矩总保持在电机的峰值转矩,轮胎纵向力输出储备小,外侧车轮滑转率过大;双级控制可自适应调节转矩分配控制层级,利用转矩分配权重降低了外侧车轮(图中1和3号车轮)的转矩输出,输出转矩峰值及平均值均较小。
应理解,本发明实施例中各步骤的序号的大小并不意味着执行顺序的先后,各过程的执行顺序应以其功能和内在逻辑确定,而不应对本发明实施例的实施过程构成任何限定。
尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由所附权利要求及其等同物限定。
表3车辆及电机参数
Claims (8)
2.根据权利要求1所述的方法,其特征在于,所述的双级控制结构的构建包括以下步骤:
所述第一级分配控制构建优化目标函数J1如下:
式中,Fx_des和Mz_des分别为上层运动控制器输出的目标控制力及目标控制力矩,Fxij代表各轮胎所受的纵向力,其中i∈{f,r}表示前后,j∈{l,r}表示左右,μij及Fzij分别表示各轮胎处路面附着系数及所受垂向载荷,d为轮距的一半;
所述地面附着约束具体为摩擦圆约束:
式中,Fyij代表各轮胎所受的侧向力;
所述第二级分配控制构建优化目标函数J2如下:
J2=||Γ1u||2+ξ||Wv(Bu-v)||2
式中,u代表各轮胎纵向力,u=[Fxfl Fxfr FxrlFxrr]T,ξ为惩罚项的权重系数,Wv=diag(WFd,WMzd),表示权重矩阵,WMzd为横摆力矩的控制权重,WFd为牵引力的控制权重;B为系数矩阵,a表示车辆的前轴距,δf为前轮转角,v表示上层运动控制器输出的目标控制力及力矩Fd,v=[Fx_des Mz_des]T;
将上式进一步整理,有
4.一种四轮毂电机驱动车辆的综合控制方法,具体包括以下步骤:
步骤1、针对四轮毂电机驱动车辆,建立七自由度车辆动力学模型;
步骤2、建立单点预瞄驾驶员模型;
步骤3、建立具有多层级结构的操纵稳定性双级最优转矩分配控制策略;
其特征在于,步骤3具体为:
步骤3.1、建立运动参考层,获得期望质心侧偏角与期望横摆角速度;
步骤3.2、设计上层运动控制器,根据期望横摆角速度、期望纵向速度获得期望纵向力与期望横摆力矩作为目标控制力及力矩;
步骤3.3、设计下层双级转矩分配控制器求解轮胎纵向力的最优解并调节电动轮输出转矩,所述下层双级转矩分配控制器执行如权利要求1-3之一所述的下层力矩分配控制方法。
6.根据权利要求5所述的方法,其特征在于,所述步骤2中所述单点预瞄驾驶员模型,指在实际行车过程中,驾驶员选择前方道路一点或一段作为参照,当驾驶员预瞄时间为T,则t时刻期望轨迹f上的参照点为A(x(t+T),f(x(t+T))),而经过T时间后,车辆将运动至B(x(t+T),y(t+T)),预瞄时间T较小因此预瞄过程中车速vx简化为恒定值,于是有:
其中车辆坐标系xOy与大地坐标系XOY的转换方式为,车辆任意时刻位置在两坐标系中的坐标(x,y)与(X,Y)存在以下关系:
8.根据权利要求7所述的方法,其特征在于,所述上层运动控制器具有以下形式:
根据实际纵向车速与人为设定的期望车速之间的偏差,利用PID控制方法控制车辆期望纵向力Fx_des,以满足车辆对期望车速跟随;
横摆力矩控制选取质心侧偏角β和横摆角速度γ作为控制变量,其中利用模糊控制方法设计修正系数K,用对质心侧偏角的检测得到系数K以修正名义横摆角速度,即γ0=Kγ;采用模糊PID控制方法对修正名义横摆角速度进行控制,模糊PID控制器的输入分别是修正名义横摆角速度γ0和横摆角速度期望值γ_des的偏差,以及将偏差进行微分处理得到的偏差率,输出是横摆力矩期望值Mz_des。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011218555.5A CN112622875B (zh) | 2020-11-04 | 2020-11-04 | 一种四轮毂电机驱动车辆的下层力矩分配控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011218555.5A CN112622875B (zh) | 2020-11-04 | 2020-11-04 | 一种四轮毂电机驱动车辆的下层力矩分配控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112622875A CN112622875A (zh) | 2021-04-09 |
CN112622875B true CN112622875B (zh) | 2022-07-19 |
Family
ID=75302951
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011218555.5A Active CN112622875B (zh) | 2020-11-04 | 2020-11-04 | 一种四轮毂电机驱动车辆的下层力矩分配控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112622875B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113246749B (zh) * | 2021-07-14 | 2021-11-05 | 北京理工大学 | 一种多轮/多轴独立驱动电动汽车多目标动力学控制方法 |
CN113479184B (zh) * | 2021-07-30 | 2023-03-17 | 南京航空航天大学 | 一种机械弹性电动轮车辆纵横向协同控制方法 |
CN113997799A (zh) * | 2021-10-26 | 2022-02-01 | 武汉理工大学 | 轮毂电机车辆驱动方法、系统、装置及存储介质 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103303157B (zh) * | 2013-06-19 | 2015-04-08 | 电子科技大学 | 一种四轮驱动电动汽车转矩分配方法 |
CN105109477B (zh) * | 2015-09-09 | 2017-06-06 | 北京理工大学 | 一种轮毂电机驱动车辆的转矩分配方法 |
GB2566500B (en) * | 2017-09-15 | 2020-04-22 | Detroit Electric Ev Tech Zhejiang Limited | Spin-control system for electrical vehicles |
CN108189705B (zh) * | 2017-12-11 | 2021-01-15 | 江苏大学 | 一种兼顾节能和稳定的分布式驱动电动汽车控制方法 |
CN109263716B (zh) * | 2018-10-10 | 2020-06-19 | 北京理工大学 | 一种四轮毂电机驱动车辆转向的控制方法 |
CN109927706A (zh) * | 2019-03-04 | 2019-06-25 | 北京精密机电控制设备研究所 | 一种轮毂电机驱动的四轮驱动电动汽车控制方法 |
CN110228462B (zh) * | 2019-05-17 | 2021-03-26 | 吉林大学 | 四轮轮毂电机驱动电动汽车横摆稳定性控制方法 |
CN111806427B (zh) * | 2020-06-08 | 2021-10-08 | 北京理工大学 | 一种四轮毂电机驱动车辆的综合控制方法 |
-
2020
- 2020-11-04 CN CN202011218555.5A patent/CN112622875B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN112622875A (zh) | 2021-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112622875B (zh) | 一种四轮毂电机驱动车辆的下层力矩分配控制方法 | |
CN111806427B (zh) | 一种四轮毂电机驱动车辆的综合控制方法 | |
CN111890951B (zh) | 智能电动汽车轨迹跟踪与运动控制方法 | |
CN109017778B (zh) | 四轮独立驱动车辆的期望路径主动转向控制方法 | |
CN109849899B (zh) | 一种电动轮汽车电液复合车身稳定控制系统及方法 | |
CN109263716B (zh) | 一种四轮毂电机驱动车辆转向的控制方法 | |
CN111002840B (zh) | 一种分布式驱动电动汽车的容错控制方法 | |
CN109017760B (zh) | 车辆期望轨迹跟踪方法、装置及滚动时域优化算法 | |
CN110605973B (zh) | 一种基于分层结构的多轴分布式电驱动车辆操纵稳定性控制方法 | |
WO2007074716A1 (ja) | 車両の制御装置 | |
CN112572411B (zh) | 一种考虑轮胎侧偏特性的车辆底盘协调控制方法及系统 | |
WO2007074717A1 (ja) | 車両の制御装置 | |
CN110239362B (zh) | 一种分布式电驱动车辆多性能优化力矩分配方法 | |
WO2007074714A1 (ja) | 車両の制御装置 | |
CN112373459B (zh) | 一种四轮毂电机驱动车辆上层运动状态控制方法 | |
CN109017759B (zh) | 期望路径车辆横摆控制方法 | |
CN112051851B (zh) | 极限工况下电动四驱车辆的自主漂移控制方法及系统 | |
CN110605975B (zh) | 一种多轴分布式电驱动车辆转矩分配集成控制器和控制方法 | |
CN113002528A (zh) | 四轮轮毂电机驱动电动汽车稳定性协调控制方法和系统 | |
CN113221257B (zh) | 考虑控制区域的极限工况下车辆横纵向稳定控制方法 | |
CN109017804B (zh) | 力矩分配控制器为车辆各轮毂电机分配驱动力矩的方法 | |
CN109017446B (zh) | 期待路径车辆纵向速度追踪控制方法及装置 | |
CN113978263A (zh) | 一种驱动轮防滑与转矩优化融合的电动汽车稳定性控制方法 | |
CN113002324A (zh) | 一种四轮独立驱动和独立转向电动汽车电子差速系统 | |
CN109017447B (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 |