[go: up one dir, main page]

CN110606079B - Layered control vehicle rollover prevention method and multi-shaft distributed driving vehicle - Google Patents

Layered control vehicle rollover prevention method and multi-shaft distributed driving vehicle Download PDF

Info

Publication number
CN110606079B
CN110606079B CN201910884468.4A CN201910884468A CN110606079B CN 110606079 B CN110606079 B CN 110606079B CN 201910884468 A CN201910884468 A CN 201910884468A CN 110606079 B CN110606079 B CN 110606079B
Authority
CN
China
Prior art keywords
control
active
vehicle
rollover
active suspension
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.)
Expired - Fee Related
Application number
CN201910884468.4A
Other languages
Chinese (zh)
Other versions
CN110606079A (en
Inventor
朱学斌
王成桥
李军求
李志超
万存才
张明星
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Taian Aerospace Special Vehicle Co ltd
Beijing Institute of Technology BIT
Original Assignee
Taian Aerospace Special Vehicle Co ltd
Beijing Institute of Technology BIT
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Taian Aerospace Special Vehicle Co ltd, Beijing Institute of Technology BIT filed Critical Taian Aerospace Special Vehicle Co ltd
Priority to CN201910884468.4A priority Critical patent/CN110606079B/en
Publication of CN110606079A publication Critical patent/CN110606079A/en
Application granted granted Critical
Publication of CN110606079B publication Critical patent/CN110606079B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT 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/00Purposes 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/02Control of vehicle driving stability
    • B60W30/04Control of vehicle driving stability related to roll-over prevention
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT 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
    • B60W40/00Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models
    • B60W40/10Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models related to vehicle motion
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT 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
    • B60W40/00Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models
    • B60W40/10Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models related to vehicle motion
    • B60W40/105Speed
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT 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/00Details 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
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT 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/00Purposes 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/02Control of vehicle driving stability
    • B60W30/04Control of vehicle driving stability related to roll-over prevention
    • B60W2030/043Control of vehicle driving stability related to roll-over prevention about the roll axis
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT 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/00Details 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/0001Details of the control system
    • B60W2050/0019Control system elements or transfer functions
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT 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/00Details 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/0001Details of the control system
    • B60W2050/0019Control system elements or transfer functions
    • B60W2050/0028Mathematical models, e.g. for simulation
    • B60W2050/0031Mathematical model of the vehicle

Landscapes

  • Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Transportation (AREA)
  • Mechanical Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Human Computer Interaction (AREA)
  • Vehicle Body Suspensions (AREA)

Abstract

本发明公开了一种分层控制的车辆防侧翻方法及多轴分布式驱动车辆。采用分层控制结构,上层控制器基于滑模变结构法进行主动侧倾力矩和差动制动的联合控制,选取侧倾角作为防侧翻控制的主要参考指标,以理想侧倾角和实际侧倾角的偏差以及理想横摆角速度和实际横摆角速度的偏差作为输入,输出主动侧倾力矩和差动制动转矩值;中层控制器基于有效集法实现主动侧倾力矩在各个悬架执行机构间的合理分配;最后下层控制器针对因中层控制而可能导致的车轮打滑情况,进行驱动防滑控制,消除垂向力变化对车轮纵向受力的消极影响。本发明采用的分层的控制结构,能够充分发挥两种防侧翻方法的优点,同时结合车辆防滑转控制,进一步提高车辆的行驶稳定性。

Figure 201910884468

The invention discloses a layer-controlled vehicle rollover prevention method and a multi-axle distributed drive vehicle. The layered control structure is adopted, and the upper controller performs joint control of active roll moment and differential braking based on the sliding mode variable structure method. The deviation of inclination angle and the deviation of ideal yaw angular velocity and actual yaw angular velocity are used as input, and output active roll moment and differential braking torque value; the middle-level controller realizes active roll moment in each suspension based on effective set method. Reasonable distribution among the actuators; finally, the lower-level controller performs driving anti-skid control for the wheel slippage that may be caused by the middle-level control, so as to eliminate the negative influence of the vertical force change on the longitudinal force of the wheel. The layered control structure adopted by the present invention can give full play to the advantages of the two anti-rollover methods, and at the same time, combined with the anti-slip control of the vehicle, the driving stability of the vehicle can be further improved.

Figure 201910884468

Description

一种分层控制的车辆防侧翻方法及多轴分布式驱动车辆A layered control vehicle rollover prevention method and multi-axle distributed drive vehicle

技术领域technical field

本发明属于电动汽车主动安全技术领域,具体地涉及一种分层控制的车辆防侧翻方法及多轴分布式驱动车辆。The invention belongs to the technical field of active safety of electric vehicles, and in particular relates to a layer-controlled vehicle rollover prevention method and a multi-axle distributed drive vehicle.

技术背景technical background

随着经济的发展,全世界汽车的数量持续增长,道路变得拥挤,汽车出行的危险变大。据统计,美国2016年全年重型卡车发生侧翻事故共15285起,仅占该车型全年事故总数的3.04%,但侧翻事故中伤亡事故却多达7285件,占总事故量的比例高达47.66%,足以看出重型车辆侧翻伤亡率是相当高的,只要发生侧翻,就会有几乎一半的可能性出现人员伤亡。With the development of the economy, the number of cars around the world continues to grow, the roads become congested, and the danger of car travel increases. According to statistics, there were 15,285 rollover accidents of heavy-duty trucks in the United States in 2016, accounting for only 3.04% of the total number of accidents of this model in the whole year. However, there were as many as 7,285 casualties in rollover accidents, accounting for as high as 47.66%, which is enough to see that the rollover casualty rate of heavy vehicles is quite high. As long as a rollover occurs, there will be almost half of the possibility of casualties.

车辆侧向稳定性是影响车辆行驶安全的重要因素,车辆侧向失稳的情况下易发生侧翻事故,这对于质心高度过高、轴距与质心高度比较小的货运车辆和特种车辆来说显得尤为危险,运载能力强所带来的附加簧上质量意味着车辆操纵稳定性的降低,使得在面临极限工况时易发生侧翻事故。The lateral stability of the vehicle is an important factor affecting the driving safety of the vehicle. In the case of lateral instability of the vehicle, a rollover accident is likely to occur. It is particularly dangerous. The additional sprung mass brought about by the strong carrying capacity means that the handling stability of the vehicle is reduced, making it prone to rollover accidents when faced with extreme working conditions.

专利CN109733382A提出一种基于模型预测控制的汽车防侧翻方法,但该方法只适用于两轴车辆的防侧翻控制;专利CN107499271A提出一种基于电控空气悬架和电控制动系统的客车防侧翻控制系统及方法,该方法通过悬架系统和制动系统协调控制的方式有效地对客车进行侧翻控制,提高客车防侧翻性能,但是对于多轴分布式驱动车辆的主动悬架垂向力的分配以及针对垂向力变化可能导致的车辆滑移情况并未进行稳定性控制。Patent CN109733382A proposes an anti-rollover method based on model predictive control, but this method is only applicable to the anti-rollover control of two-axle vehicles; Patent CN107499271A proposes a passenger car based on electronically controlled air suspension and electronically controlled braking system An anti-rollover control system and method, the method effectively controls the rollover of a passenger car by means of coordinated control of a suspension system and a braking system, and improves the anti-rollover performance of the passenger car. The distribution of the vertical force and the stability control for the possible vehicle slippage caused by the change of the vertical force is not carried out.

发明内容SUMMARY OF THE INVENTION

本发明为解决背景技术中所涉及到的缺陷,提供一种分层控制的车辆防侧翻方法及多轴分布式驱动车辆。本发明采用三层的分层控制方法,针对多轴分布式驱动车辆,通过主动侧倾力矩和差动制动联合控制实现防侧翻,同时实现纵向防滑控制,避免中层主动侧倾力矩分配控制导致的不良影响。In order to solve the defects involved in the background technology, the present invention provides a layered controlled vehicle rollover prevention method and a multi-axle distributed drive vehicle. The present invention adopts a three-layer layered control method, aiming at a multi-axle distributed drive vehicle, realizes anti-rollover through the joint control of active roll moment and differential braking, simultaneously realizes longitudinal anti-skid control, and avoids the active roll moment in the middle layer. Undesirable effects due to distribution control.

本发明的分层控制的车辆防侧翻方法,包括如下步骤:The layered control vehicle rollover prevention method of the present invention comprises the following steps:

测得车辆的方向盘转角、车速,输入给上层防侧翻联合控制系统;Measure the steering wheel angle and speed of the vehicle, and input them to the upper anti-rollover joint control system;

所述上层防侧翻联合控制系统计算主动侧倾力矩和差动制动力矩,并将计算出的主动侧倾力矩作为中层悬架垂向力优化分配系统的输入变量;The upper-layer anti-rollover joint control system calculates the active roll moment and the differential braking moment, and uses the calculated active roll moment as the input variable of the middle-layer suspension vertical force optimization distribution system;

所述中层悬架垂向力优化分配系统实现对于主动悬架执行器的垂向力分配;The middle-level suspension vertical force optimization distribution system realizes the vertical force distribution for the active suspension actuator;

下层驱动与制动防滑控制系统基于差动制动力矩,以理想滑移率为目标,基于滑模变结构法计算驱动转矩。The lower-level driving and braking anti-skid control system is based on the differential braking torque, aiming at the ideal slip rate, and calculating the driving torque based on the sliding mode variable structure method.

优选地,所述上层防侧翻联合控制系统以方向盘转角信号作为输入变量,通过车辆三自由度参考模型计算输出变量理想侧倾角和理想横摆角速度。Preferably, the upper-layer anti-rollover joint control system uses the steering wheel angle signal as an input variable, and calculates the output variables ideal roll angle and ideal yaw rate through a vehicle three-degree-of-freedom reference model.

优选地,基于模糊增益调节的滑模变结构控制算法计算主动侧倾力矩。Preferably, the active roll moment is calculated based on a sliding mode variable structure control algorithm with fuzzy gain adjustment.

优选地,基于滑模变结构控制算法计算差动制动力矩。Preferably, the differential braking torque is calculated based on a sliding mode variable structure control algorithm.

优选地,主动悬架执行器的垂向力分配具体为对主动防侧倾力矩在前桥主动悬架执行器与中后桥主动悬架执行器之间分配比例,以及左侧主动悬架执行器与右侧主动悬架执行器之间分配比例的控制。Preferably, the vertical force distribution of the active suspension actuators is specifically the distribution ratio of the active anti-roll moment between the front axle active suspension actuators and the middle and rear axle active suspension actuators, and the left active suspension Control of proportional distribution between the actuator and the right active suspension actuator.

优选地,基于有效集法实现对于主动悬架执行器的垂向力分配。Preferably, the vertical force distribution to the active suspension actuators is achieved based on the active set method.

优选地,所述基于有效集法实现对于主动悬架执行器的垂向力分配,具体包括如下步骤:Preferably, the realization of the vertical force distribution for the active suspension actuator based on the active set method specifically includes the following steps:

设置初始时刻总控制转矩在前桥与中后桥间、左侧轮与右侧轮之间均为各轮处执行器平均分配;Set the total control torque between the front axle and the middle and rear axles, and between the left and right wheels at the initial moment, all of which are equally distributed by the actuators at each wheel;

确定有效集法的设计变量X的约束条件,所述约束条件考虑车轮在良好路面的附着条件、横向载荷转移率要求、垂向动态力对车身的拉压要求、主动侧倾力矩的输出要求、以及主动悬架执行器的输出能力;Determine the constraints of the design variable X of the active set method, which consider the adhesion conditions of the wheels on the good road surface, the lateral load transfer rate requirements, the tension and compression requirements of the vertical dynamic force on the body, and the output requirements of the active roll moment. , and the output capability of the active suspension actuator;

在迭代可行点周围确定所有约束中的有效约束集合;Determine a valid set of constraints among all constraints around an iterative feasible point;

经过多次迭代得到在上述约束下的目标函数的最优解,所述目标函数为各轮胎平均相对动载,得到所述主动悬架执行器的最优垂向力分配比例,计算左前桥、右前桥、左中后桥和右中后桥的主动悬架执行器所需的垂向控制力。After several iterations, the optimal solution of the objective function under the above constraints is obtained, the objective function is the average relative dynamic load of each tire, the optimal vertical force distribution ratio of the active suspension actuator is obtained, and the left front axle, The vertical control force required by the active suspension actuators of the right front axle, left middle rear axle and right middle rear axle.

优选地,述计算主动侧倾力矩具体包括如下步骤:Preferably, the calculation of the active roll moment specifically includes the following steps:

计算侧倾角偏差e及其偏差率

Figure BDA0002206265440000021
其计算过程如下:Calculate the roll angle deviation e and its deviation rate
Figure BDA0002206265440000021
Its calculation process is as follows:

e=φ0e=φ 0

Figure BDA0002206265440000022
Figure BDA0002206265440000022

选用指数趋近律,设计滑模变结构控制器的线性切换函数s;The exponential reaching law is used to design the linear switching function s of the sliding mode variable structure controller;

Figure BDA0002206265440000023
Figure BDA0002206265440000023

式中,c为常数,ε为切换系数,k表示滑模切换增益,且c>0,ε>0,k>0;In the formula, c is a constant, ε is the switching coefficient, k is the sliding mode switching gain, and c>0, ε>0, k>0;

sgn表示符号函数,其表达式为:sgn represents a symbolic function whose expression is:

Figure BDA0002206265440000024
Figure BDA0002206265440000024

最终的滑模控制器设计为;The final sliding mode controller is designed as;

Figure BDA0002206265440000025
Figure BDA0002206265440000025

式中,

Figure BDA0002206265440000026
由基于39DOF车辆模型的上层防侧翻控制系统的控制方程得到:In the formula,
Figure BDA0002206265440000026
The control equation of the upper-layer anti-rollover control system based on the 39DOF vehicle model is obtained:

Figure BDA0002206265440000027
Figure BDA0002206265440000027

其中,确定切换系数ε的方法为:模糊控制器中输入变量侧倾角偏差e的论域为[0,15],输出变量切换增益ε的论域为[0,1.5];侧倾角偏差e、切换增益ε的语言值集合都为{ZO,PM,PB},表示自然语言中的{零,正中,正大};所述输入变量和输出变量隶属度函数选用三角-S型混合隶属度函数,解模糊方法选用面积重心法。优选地,所述下层驱动与制动防滑控制系统以20%的滑移率对应的旋转角速度为参考目标。Among them, the method of determining the switching coefficient ε is as follows: the domain of discourse of the input variable roll angle deviation e in the fuzzy controller is [0, 15], and the domain of discourse of the output variable switching gain ε is [0, 1.5]; The language value sets of the switching gain ε are all {ZO, PM, PB}, representing {zero, positive, positive} in natural language; the membership function of the input variable and output variable is a triangular-S type mixed membership function, The defuzzification method adopts the area centroid method. Preferably, the lower-level driving and braking anti-skid control system takes the rotation angular velocity corresponding to the slip ratio of 20% as the reference target.

本发明还涉及一种多轴分布式驱动车辆,其特征在于采用如前述分层控制的车辆防侧翻方法。The present invention also relates to a multi-axle distributed drive vehicle, which is characterized by adopting the above-mentioned layered control vehicle rollover prevention method.

本发明具有如下有益的效果:The present invention has the following beneficial effects:

1.采用分层稳定性控制框架,中层控制系统可实现上层控制系统计算出的主动侧倾力矩在多轴分布式驱动车辆的主动悬架间的优化分配,下层控制系统针对中层垂向力变化现象与中层垂向力变化可能导致的车辆滑移情况,进行纵向防滑控制以避免中层分配控制的不良影响,实现车辆在行驶过程中综合性能要求。1. Using a layered stability control framework, the middle-layer control system can realize the optimal distribution of the active roll moment calculated by the upper-layer control system among the active suspensions of the multi-axle distributed drive vehicle, and the lower-layer control system can target the middle-layer vertical force For the vehicle slippage that may be caused by the change phenomenon and the change of the vertical force of the middle layer, the longitudinal anti-skid control is carried out to avoid the adverse effect of the distribution control of the middle layer, and the comprehensive performance requirements of the vehicle during the driving process are realized.

2.采用主动侧倾力矩和差动制动的联合控制,可以充分发挥两种控制方法的优势,并解决单一方法使用所带来的不足之处。2. The joint control of active roll moment and differential braking can give full play to the advantages of the two control methods and solve the shortcomings brought by the use of a single method.

具体来讲,当车辆侧倾姿态偏移较小时,侧倾控制不足以对横摆参数造成影响,这时首先采用主动侧倾力矩系统进行控制,充分发挥悬架控制直接有效的优势,而在车身侧倾较为严重,或者主动侧倾力矩控制明显改变了车辆横摆角速度变化规律时,差动制动系统参与工况,一方面通过制动减速进一步提高车辆的抗侧倾能力,另一方面通过差动附加横摆力矩的作用,补偿侧倾力矩控制在横摆方向带来的消极影响。Specifically, when the vehicle roll attitude deviation is small, the roll control is not enough to affect the yaw parameters. At this time, the active roll moment system is firstly used for control to give full play to the direct and effective advantages of suspension control, while the When the body roll is more serious, or the active roll moment control obviously changes the vehicle yaw rate change law, the differential braking system participates in the working condition. On the one hand, the negative influence of the control of the rolling moment in the yaw direction is compensated by the action of the differential additional yaw moment.

3.采用基于模糊控制的方法确定主动侧倾力矩滑模变结构控制中的增益系数,相较与传统滑模控制系统中的常数切换增益,能随着控制效果的改变而不断变化,能够一定程度上的改善控制过程中的抖动现象。3. The method based on fuzzy control is used to determine the gain coefficient in the active roll moment sliding mode variable structure control. Compared with the constant switching gain in the traditional sliding mode control system, it can change continuously with the change of the control effect. To a certain extent, the jitter phenomenon in the control process is improved.

4.取侧倾角和横摆角速度作为防侧翻控制参考指标,与选取横向载荷转移率LTR作为指标相比,能够根据控制参考指标的偏差,来直接实现对于车辆的侧倾稳定性控制。4. Taking the roll angle and yaw rate as the reference index of anti-rollover control, compared with selecting the lateral load transfer rate LTR as the index, the roll stability control of the vehicle can be directly realized according to the deviation of the control reference index.

附图说明Description of drawings

图1为本发明分层防侧翻协调控制系统原理示意图1 is a schematic diagram of the principle of the layered anti-rollover coordination control system of the present invention

图2为本发明基于模糊-滑模控制的车辆防侧翻控制原理框图Fig. 2 is a block diagram of the vehicle rollover prevention control principle based on fuzzy-sliding mode control of the present invention

图3为本发明有效集方法进行主动侧倾力矩优化分配流程示意图FIG. 3 is a schematic diagram of the optimal allocation process of active roll moment in the active set method of the present invention.

具体实施方式Detailed ways

下面结合附图对本发明的技术方案做进一步的详细说明:Below in conjunction with accompanying drawing, the technical scheme of the present invention is described in further detail:

如图1所示,本发明所述分层控制结构,左上角虚线框内为上层防侧翻联合控制系统,采用双闭环的控制结构,计算出主动侧倾力矩和差动制动力矩用于车辆防侧翻控制,然后将计算出的主动侧倾力矩作为中层悬架垂向力优化分配系统的输入变量,将主动侧倾力矩分配给各个主动悬架系统,最后下层驱动与制动防滑控制系统针对中层悬架垂向力分配优化控制可能导致的车轮打滑现象,以滑移率为控制目标进行驱动防滑控制。本发明的控制过程具体包含以下步骤:As shown in Figure 1, the layered control structure of the present invention, the upper left corner of the dashed box is the upper-layer anti-rollover joint control system, using a double closed-loop control structure to calculate the active roll moment and the differential braking moment. For the vehicle anti-rollover control, the calculated active roll moment is used as the input variable of the vertical force distribution system of the middle suspension, and the active roll moment is distributed to each active suspension system. The dynamic anti-skid control system optimizes and controls the wheel slip phenomenon that may be caused by the vertical force distribution of the mid-level suspension, and uses the slip rate to control the driving anti-skid control. The control process of the present invention specifically comprises the following steps:

步骤1),测得车辆的方向盘转角δs、车速v、横摆角速度r、侧倾角φ,然后输入给上层防侧翻联合控制系统;Step 1), measure the steering wheel angle δ s , vehicle speed v, yaw rate r, roll angle φ of the vehicle, and then input them to the upper-layer anti-rollover joint control system;

所述上层防侧翻联合控制系统的控制过程包括以下步骤:The control process of the upper-layer anti-rollover joint control system includes the following steps:

步骤2),以方向盘转角信号δs作为输入变量,通过车辆三自由度参考模型计算输出变量理想侧倾角φ0和理想横摆角速度r0Step 2), taking the steering wheel angle signal δs as the input variable, and calculating the output variables ideal roll angle φ 0 and ideal yaw rate r 0 through the vehicle three-degree-of-freedom reference model;

以八轴分布式驱动车辆为例,步骤2)的具体计算过程如下:Taking an eight-axis distributed drive vehicle as an example, the specific calculation process of step 2) is as follows:

以侧倾角φ、侧向速度vy、横摆角速度r及侧倾角速度

Figure BDA0002206265440000041
作为状态变量,即
Figure BDA0002206265440000042
Figure BDA0002206265440000043
方向盘转角信号δs作为输入变量u,以理想侧倾角φ0和理想横摆角速度r0为输出变量y;With roll angle φ, lateral velocity v y , yaw angular velocity r and roll angular velocity
Figure BDA0002206265440000041
as a state variable, i.e.
Figure BDA0002206265440000042
Figure BDA0002206265440000043
The steering wheel angle signal δ s is used as the input variable u, and the ideal roll angle φ 0 and the ideal yaw rate r 0 are used as the output variable y;

建立三自由度整车状态方程如下:The three-degree-of-freedom vehicle state equation is established as follows:

Figure BDA0002206265440000044
Figure BDA0002206265440000044

其中,A=-F-1G,B=F-1HC=[1 0 1 0]Where, A=-F -1 G, B=F -1 HC=[1 0 1 0]

Figure BDA0002206265440000045
Figure BDA0002206265440000045

Figure BDA0002206265440000046
Figure BDA0002206265440000046

式中,m、ms分别表示整车质量和悬挂质量;where m and m s represent the vehicle mass and suspension mass, respectively;

vx、vy分别表示车辆的纵向和侧向车速;v x and v y represent the longitudinal and lateral speeds of the vehicle, respectively;

φ、

Figure BDA0002206265440000047
分别表示侧倾角、侧倾角速度、侧倾角加速度;φ,
Figure BDA0002206265440000047
respectively represent the roll angle, roll angular velocity, and roll angular acceleration;

r、

Figure BDA0002206265440000048
表示车辆绕z轴的横摆角速度、横摆角加速度;r.
Figure BDA0002206265440000048
Indicates the yaw angular velocity and yaw angular acceleration of the vehicle around the z-axis;

Ixz表示车辆对ox与oz轴的惯性积;I xz represents the inertia product of the vehicle on the ox and oz axes;

Ix表示车辆绕x轴的转动惯量;I x represents the moment of inertia of the vehicle around the x-axis;

Iz表示车辆绕z轴的转动惯量;I z represents the moment of inertia of the vehicle around the z-axis;

Li表示第i轴到质心的距离;Li represents the distance from the i -th axis to the centroid;

Ci表示第i轴悬架侧倾角刚度;C i represents the roll angle stiffness of the i-th axle suspension;

h表示悬挂质量重心到侧倾轴线的距离;h represents the distance from the center of gravity of the suspension to the roll axis;

Ksi、Csi分别表示第i个悬架的刚度和阻尼。K si and C si represent the stiffness and damping of the ith suspension, respectively.

上层防侧翻控制系统的控制方程如下所示建立:The control equation of the upper anti-rollover control system is established as follows:

以xc=[r ωx]T为目标状态变量,其中ωx、r分别为上层防侧翻联合控制系统中主动侧倾力矩控制和差动制动控制所对应的状态变量。以uc=[FX MX]T作为目标控制输入变量,其中FX表示差动制动力矩,MX表示主动侧倾力矩;Take x c =[r ω x ] T as the target state variable, where ω x and r are the state variables corresponding to active roll moment control and differential braking control in the upper-layer anti-rollover joint control system, respectively. Take uc = [F X M X ] T as the target control input variable, where F X represents the differential braking torque, and M X represents the active roll moment;

基于待分析车辆模型得到非线性状态方程,得到控制方程;The nonlinear state equation is obtained based on the vehicle model to be analyzed, and the control equation is obtained;

优选地基于39DOF车辆模型得到非线性状态方程,得到的控制方程如下:The nonlinear state equation is preferably obtained based on the 39DOF vehicle model, and the obtained control equation is as follows:

Figure BDA0002206265440000051
Figure BDA0002206265440000051

式中,u、v、vz分别表示车辆的纵向、侧向和垂向车速;In the formula, u, v , vz represent the longitudinal, lateral and vertical speeds of the vehicle, respectively;

ωx、ωy、r分别表示车身绕车辆坐标系x轴、y轴和z轴的旋转角速度;ω x , ω y , and r represent the rotational angular velocity of the vehicle body around the x-axis, y-axis and z-axis of the vehicle coordinate system, respectively;

Fx、Fy、Fz分别表示车辆在X、Y、Z坐标的纵向、侧向和垂向受力;F x , F y , and F z represent the longitudinal, lateral and vertical forces of the vehicle at the X, Y, and Z coordinates, respectively;

Mz、Mx、My分别代表车辆在z、x、y轴方向所受的横摆力矩;M z , M x , and M y represent the yaw moment of the vehicle in the directions of the z, x, and y axes, respectively;

M、mb分别表示整车和车身质量;M and m b represent the vehicle and body mass, respectively;

Ix、Iy、Iz分别代表车辆绕x轴、y轴和z轴的转动惯量;I x , I y , and I z represent the moment of inertia of the vehicle around the x-axis, the y-axis and the z-axis, respectively;

αs和f分别表示车道坡度和滚动阻力系数; αs and f are the lane gradient and rolling resistance coefficient, respectively;

CD和A分别表示空气阻力系数和迎风面积;C D and A represent the air resistance coefficient and the windward area, respectively;

φ和ψ分别表示车身的侧倾和俯仰角;φ and ψ represent the roll and pitch angles of the body, respectively;

he表示簧上质量质心至侧倾中心距离。h e is the distance from the sprung mass center to the roll center.

步骤3),基于模糊增益调节的滑模变结构控制算法计算主动侧倾力矩;Step 3), calculating active roll moment based on the sliding mode variable structure control algorithm of fuzzy gain adjustment;

如图2所示,基于模糊增益调节的滑模变结构控制算法计算主动侧倾力矩的原理框图.As shown in Figure 2, the principle block diagram of calculating the active roll moment based on the sliding mode variable structure control algorithm based on fuzzy gain adjustment.

步骤3)具体步骤如下:Step 3) The specific steps are as follows:

步骤3.1.1),计算侧倾角偏差e及其偏差率

Figure BDA0002206265440000052
其计算过程如下:Step 3.1.1), calculate the roll angle deviation e and its deviation rate
Figure BDA0002206265440000052
Its calculation process is as follows:

e=φ0e=φ 0

Figure BDA0002206265440000053
Figure BDA0002206265440000053

步骤3.1.2),选用指数趋近律,设计滑模变结构控制器的线性切换函数s;Step 3.1.2), select the exponential reaching law to design the linear switching function s of the sliding mode variable structure controller;

Figure BDA0002206265440000054
Figure BDA0002206265440000054

式中,c为常数,ε为切换系数,k表示滑模切换增益,且c>0,ε>0,k>0;In the formula, c is a constant, ε is the switching coefficient, k is the sliding mode switching gain, and c>0, ε>0, k>0;

sign表示符号函数,其表达式为:sign represents a symbolic function whose expression is:

Figure BDA0002206265440000055
Figure BDA0002206265440000055

步骤3.1.3),最终的滑模控制器设计为;Step 3.1.3), the final sliding mode controller is designed as;

Figure BDA0002206265440000061
Figure BDA0002206265440000061

式中,

Figure BDA0002206265440000062
由上述基于39DOF车辆模型的上层防侧翻控制系统的控制方程得到:In the formula,
Figure BDA0002206265440000062
It is obtained from the above control equation of the upper-layer anti-rollover control system based on the 39DOF vehicle model:

Figure BDA0002206265440000063
Figure BDA0002206265440000063

其中,采用模糊控制法确定步骤3.1.2)的切换系数ε,具体步骤如下:Among them, the fuzzy control method is used to determine the switching coefficient ε in step 3.1.2), and the specific steps are as follows:

步骤3.2),模糊控制器中输入变量侧倾角偏差e的论域为[0,15],输出变量切换增益ε的论域为[0,1.5];侧倾角偏差e、切换增益ε的语言值集合都为{ZO,PM,PB},表示自然语言中的{零,正中,正大};所述输入变量和输出变量隶属度函数选用三角-S型混合隶属度函数,解模糊方法选用面积重心法。Step 3.2), the domain of universe of the input variable roll angle deviation e in the fuzzy controller is [0, 15], and the domain of discourse of the output variable switching gain ε is [0, 1.5]; the language values of the roll angle deviation e and switching gain ε The sets are all {ZO, PM, PB}, which represent {zero, positive, positive} in natural language; the membership functions of the input variables and output variables use the triangular-S type mixed membership function, and the defuzzification method uses the area center of gravity Law.

步骤4),基于滑模变结构控制算法计算差动制动力矩FxStep 4), calculate the differential braking torque F x based on the sliding mode variable structure control algorithm;

步骤4)具体步骤如下:Step 4) The specific steps are as follows:

步骤4.1),计算横摆角速度偏差er及其偏差率

Figure BDA0002206265440000064
其计算过程如下:Step 4.1), calculate the yaw rate deviation er and its deviation rate
Figure BDA0002206265440000064
Its calculation process is as follows:

er=r0-rer = r 0 -r

Figure BDA0002206265440000065
Figure BDA0002206265440000065

步骤4.2),选用指数趋近律,设计滑模变结构控制器的线性切换函数s;Step 4.2), select the exponential reaching law to design the linear switching function s of the sliding mode variable structure controller;

Figure BDA0002206265440000066
Figure BDA0002206265440000066

式中,c为常数,ε为切换系数,k表示滑模切换增益,且c>0,ε>0,k>0;In the formula, c is a constant, ε is the switching coefficient, k is the sliding mode switching gain, and c>0, ε>0, k>0;

sgn表示符号函数,其表达式为:sgn represents a symbolic function whose expression is:

Figure BDA0002206265440000067
Figure BDA0002206265440000067

步骤4.3),最终的滑模控制器设计为;Step 4.3), the final sliding mode controller is designed as;

Figure BDA0002206265440000068
Figure BDA0002206265440000068

式中,g′、f′由上述基于39DOF车辆模型的上层防侧翻控制系统的控制方程得到:In the formula, g′ and f′ are obtained from the above control equations of the upper-layer anti-rollover control system based on the 39DOF vehicle model:

Figure BDA0002206265440000069
Figure BDA0002206265440000069

Figure BDA00022062654400000610
Figure BDA00022062654400000610

在前述上层防侧翻联合控制系统作用之下,防止车身侧翻所需的主动侧倾力矩经由模糊控制-滑模变结构控制计算得到,该主动侧倾矩能够使得车辆侧倾角能按照参考轨迹的变化规律实现安全的侧向运动。Under the action of the above-mentioned upper-layer anti-rollover joint control system, the active roll moment required to prevent the body from rolling over is calculated through fuzzy control-sliding mode variable structure control. The active roll moment can make the vehicle roll angle according to the reference The changing law of the trajectory realizes safe lateral movement.

总主动侧倾力矩对车辆系统上的控制最终是通过主动悬架执行机构来实现的,故而采用中层悬架垂向力分配优化控制系统实现对于主动悬架执行器的垂向力分配控制。The control of the total active roll moment on the vehicle system is finally realized by the active suspension actuator, so the vertical force distribution control of the active suspension actuator is realized by the vertical force distribution optimization control system of the mid-level suspension.

本发明将整车主动悬架执行机构简化到各个车轮处并依照车轮位置进行分区,主要分为支撑驾驶员舱的左前桥执行区和右前桥执行区,以及承载货物的左中后桥执行区和右中后桥执行区。前桥执行区包括前述左前桥执行区和右前桥执行区,中后桥执行区包括前述左中后桥执行区和右中后桥执行区,左侧执行区包括前述左前桥执行区和左中后桥执行区,右侧执行区包括前述前右桥执行区和右中后桥执行区。The present invention simplifies the vehicle active suspension actuator to each wheel and divides it according to the position of the wheel. It is mainly divided into the left front axle implementation area and the right front axle implementation area that support the driver's cabin, and the left middle rear axle implementation area that carries the cargo. and the execution area of the right middle rear axle. The front axle execution area includes the foregoing left front axle execution area and the right front axle execution area, the middle and rear axle execution area includes the foregoing left middle rear axle execution area and the right middle rear axle execution area, and the left execution area includes the foregoing left front axle execution area and left middle execution area. Rear axle execution area, right execution area includes the aforementioned front right axle execution area and right middle rear axle execution area.

故而对于主动悬架执行器的垂向力分配可以转化为对主动防侧倾力矩在前桥执行区与中后桥执行区之间分配比例,以及左侧执行区与右侧执行区之间分配比例的控制。即于主动悬架执行器的分配控制可以转化为对主动防侧倾力矩在前桥主动悬架执行器与中后桥主动悬架执行器之间分配比例,以及左侧主动悬架执行器与右侧主动悬架执行器之间分配比例的控制。Therefore, the vertical force distribution for the active suspension actuator can be transformed into the distribution ratio of the active anti-roll moment between the front axle execution area and the middle and rear axle execution areas, as well as between the left and right execution areas. Control of distribution ratio. That is, the distribution control of the active suspension actuator can be transformed into the distribution ratio of the active anti-roll moment between the active suspension actuator of the front axle and the active suspension actuator of the middle and rear axle, and the active suspension actuator of the left side. Controls proportional distribution with the right active suspension actuator.

在主动侧倾力矩控制的作用下,一方面车身的侧向姿态可以得到有效的调节,另一方面车轮处的动态垂向力会发生改变,具体来讲,悬架受压侧压力增大,受拉侧拉力增大。故而从提高行车安全性的角度考虑,以减小车轮垂向相对动态载荷为目标进行优化分配控制。Under the action of active roll moment control, on the one hand, the lateral posture of the body can be effectively adjusted, on the other hand, the dynamic vertical force at the wheels will change. Specifically, the pressure on the compressed side of the suspension increases. , the tension on the tension side increases. Therefore, from the perspective of improving the driving safety, the optimal distribution control is carried out with the goal of reducing the vertical relative dynamic load of the wheel.

所述中层悬架垂向力分配优化系统的控制过程包括以下步骤:The control process of the vertical force distribution optimization system for the mid-level suspension includes the following steps:

步骤5),以步骤3)计算出的主动侧倾力矩作为中层控制器的输入,基于有效集法对各个主动悬架执行器进行主动防侧倾力矩优化分配;Step 5), using the active roll moment calculated in step 3) as the input of the middle-level controller, and based on the active set method, perform an optimal distribution of the active anti-roll moment for each active suspension actuator;

基于上述的主动悬架执行机构分区,设计变量X为主动侧倾力矩优化分配比例:Based on the above active suspension actuator partition, the design variable X is the optimal distribution ratio of active roll moment:

X=[kx1,kx2]T X=[k x1 ,k x2 ] T

其中,kx1为前桥主动悬架执行器、与中后桥主动悬架执行器之间的力矩分配比例,kx2为左侧主动悬架执行器与右侧主动悬架执行器之间的力矩分配比例。Among them, k x1 is the torque distribution ratio between the front axle active suspension actuator and the middle rear axle active suspension actuator, and k x2 is the ratio between the left active suspension actuator and the right active suspension actuator. Torque distribution ratio.

如图3所示的中层悬架垂向力分配优化系统的有效集法流程原理图,其具体包括以下步骤:Figure 3 shows the schematic diagram of the effective set method flow chart of the vertical force distribution optimization system of the middle suspension, which specifically includes the following steps:

步骤5.1),设置仿真初始点为X=[0.25,0.5],表示初始时刻总控制转矩在前桥与中后桥间、左侧轮与右侧轮之间均为各轮处执行器平均分配。Step 5.1), set the initial point of the simulation as X=[0.25, 0.5], indicating that the total control torque at the initial moment is the average of the actuators at each wheel between the front axle and the middle and rear axle, and between the left and right wheels. distribute.

步骤5.2),确定设计变量X的约束条件。所述约束条件考虑车轮在良好路面的附着条件、横向载荷转移率要求、垂向动态力对车身的拉压要求、主动侧倾力矩的输出要求、主动悬架执行器的输出能力等,其数学表述简写形式为

Figure BDA0002206265440000071
Figure BDA0002206265440000072
具体表达如下:Step 5.2), determine the constraints of the design variable X. The constraint conditions consider the adhesion conditions of the wheels on a good road surface, the lateral load transfer rate requirements, the tension and compression requirements of the vertical dynamic force on the body, the output requirements of the active roll moment, the output capability of the active suspension actuator, etc. The short form of mathematical expression is
Figure BDA0002206265440000071
Figure BDA0002206265440000072
The specific expression is as follows:

Figure BDA0002206265440000073
Figure BDA0002206265440000073

其中,Fxi表示每个车轮所受纵向力,Fzsdi表示每个车轮处的垂向总载荷,μ表示路面附着系数,Fzli表示左侧执行区第i轴车轮处所受垂向载荷,Fzri表示右侧执行区第i轴车轮处所受垂向载荷,ΔFzi表示每个车轮处的垂向力,Fmax表示单个车轮垂向力极限值,Kx1、Kx1则分别表示分配系数的边界限制,受到主动悬架执行器输出能力的影响,n表示总轴数,LTRmax为最大横向载荷转移率;Among them, F xi represents the longitudinal force on each wheel, F zsdi represents the vertical total load at each wheel, μ represents the road adhesion coefficient, F zli represents the vertical load on the wheel of the i-th axis in the left execution area, F zri represents the vertical load on the wheel of the i-th axis in the right execution area, ΔF zi represents the vertical force at each wheel, F max represents the limit value of the vertical force of a single wheel, K x1 , K x1 respectively represent the distribution The boundary limit of the coefficient is affected by the output capacity of the active suspension actuator, n represents the total number of axles, and LTR max is the maximum lateral load transfer rate;

步骤5.3),在迭代可行点周围确定所有约束中的有效约束集合;Step 5.3), determine the effective set of constraints in all constraints around the iterative feasible point;

步骤5.4),经过多次迭代得到在上述约束下的目标函数的最优解X,目标函数为各轮胎平均相对动载,其表达式如下:Step 5.4), after several iterations, the optimal solution X of the objective function under the above constraints is obtained, and the objective function is the average relative dynamic load of each tire, and its expression is as follows:

Figure BDA0002206265440000081
Figure BDA0002206265440000081

其中,Fzdi表示车辆每个车轮处的垂向动态载荷,Fzsdi表示每个车轮处的垂向总载荷,m表示轮胎总个数。Among them, F zdi represents the vertical dynamic load at each wheel of the vehicle, F zsdi represents the vertical total load at each wheel, and m represents the total number of tires.

步骤5.5),由步骤5.4)计算得到所述主动悬架执行器的最优垂向力分配比例,计算左前桥、右前桥、左中后桥和右中后桥的主动悬架执行器所需的垂向控制力ΔFzfl、ΔFzfr、ΔFzrl、ΔFzrr,并进行主动悬架垂向力控制,各主动悬架垂向力的具体计算过程如下:Step 5.5), the optimal vertical force distribution ratio of the active suspension actuator is calculated from step 5.4), and the required amount of active suspension actuators for the left front axle, right front axle, left middle rear axle and right middle rear axle is calculated. The vertical control forces ΔF zfl , ΔF zfr , ΔF zrl , and ΔF zrr are used to control the vertical force of the active suspension. The specific calculation process of the vertical force of each active suspension is as follows:

Figure BDA0002206265440000082
Figure BDA0002206265440000082

式中,B为轮距;In the formula, B is the wheelbase;

通过对主动侧倾力矩的中层优化分配,得到了整车主动悬架的最优力矩分配比例,由于主动悬架的控制对于车轮垂向力的影响,车轮相对动载减小,垂向安全性得到提高。在垂向力因控制分配而发生变化的同时,考虑到轮胎垂向力对纵向力取值的限制作用,车轮的纵向输出能力会不可避免地发生改变。针对因中层垂向力控制而可能导致的车轮打滑情况,下层驱动与制动防滑控制系统以理想滑移率为目标,基于滑模变结构法进行驱动防滑控制。Through the middle-level optimal distribution of the active roll moment, the optimal moment distribution ratio of the vehicle's active suspension is obtained. Due to the influence of the active suspension control on the vertical force of the wheel, the relative dynamic load of the wheel is reduced, and the vertical safety is Sex is improved. When the vertical force changes due to the control distribution, the longitudinal output capacity of the wheel will inevitably change considering the limiting effect of the vertical force of the tire on the value of the longitudinal force. Aiming at the possible wheel slip caused by the vertical force control of the middle layer, the lower driving and braking anti-skid control system takes the ideal slip rate as the target, and performs the driving anti-skid control based on the sliding mode variable structure method.

所述下层驱动与制动防滑控制系统的控制过程包括以下步骤:The control process of the lower driving and braking anti-skid control system includes the following steps:

步骤6),以理想滑移率为目标,基于滑模变结构法计算驱动转矩TwStep 6), with the ideal slip rate as the target, the driving torque Tw is calculated based on the sliding mode variable structure method.

步骤6.1),以20%的滑移率对应的旋转角速度w0为参考目标,计算旋转角速度偏差ew及其偏差率

Figure BDA0002206265440000083
其计算过程如下:Step 6.1), take the rotation angular velocity w 0 corresponding to 20% slip rate as the reference target, calculate the rotation angular velocity deviation e w and its deviation rate
Figure BDA0002206265440000083
Its calculation process is as follows:

ew=w0-we w =w 0 -w

Figure BDA0002206265440000084
Figure BDA0002206265440000084

步骤6.2),选用指数趋近律,设计滑模变结构控制器的线性切换函数s;Step 6.2), select the exponential reaching law to design the linear switching function s of the sliding mode variable structure controller;

Figure BDA0002206265440000085
Figure BDA0002206265440000085

式中,c为常数,ε为切换系数,k表示滑模切换增益,且c>0,ε>0,k>0;In the formula, c is a constant, ε is the switching coefficient, k is the sliding mode switching gain, and c>0, ε>0, k>0;

sgn表示符号函数,其表达式为:sgn represents a symbolic function whose expression is:

Figure BDA0002206265440000086
Figure BDA0002206265440000086

步骤6.3),最终的滑模控制器设计为;Step 6.3), the final sliding mode controller is designed as;

Figure BDA0002206265440000091
Figure BDA0002206265440000091

式中,

Figure BDA0002206265440000092
由单个轮胎控制所需状态方程得到:In the formula,
Figure BDA0002206265440000092
The equation of state required for the control of a single tire is obtained:

Figure BDA0002206265440000093
Figure BDA0002206265440000093

Figure BDA0002206265440000094
Figure BDA0002206265440000094

其中,控制输出u′=Tw,状态变量x=ω,

Figure BDA0002206265440000095
Among them, the control output u '=Tw, the state variable x=ω,
Figure BDA0002206265440000095

本技术领域技术人员可以理解的是,除非另外定义,这里使用的所有术语(包括技术术语和科学术语)具有与本发明所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是,诸如通用字典中定义的那些术语应该被理解为具有与现有技术的上下文中的意义一致的意义,并且除非像这里一样定义,不会用理想化或过于正式的含义来解释。It will be understood by those skilled in the art that, unless otherwise defined, all terms (including technical and scientific terms) used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. It should also be understood that terms such as those defined in general dictionaries should be understood to have meanings consistent with their meanings in the context of the prior art and, unless defined as herein, are not to be taken in an idealized or overly formal sense. explain.

以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。The specific embodiments described above further describe the purpose, technical solutions and beneficial effects of the present invention in further detail. It should be understood that the above descriptions are only specific embodiments of the present invention, and are not intended to limit the present invention. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included within the protection scope of the present invention.

Claims (9)

1.一种分层控制的车辆防侧翻方法,其特征在于包含如下步骤:1. a vehicle rollover prevention method of layered control is characterized in that comprising the steps: 测得车辆的方向盘转角、车速,输入给上层防侧翻联合控制系统;Measure the steering wheel angle and speed of the vehicle, and input them to the upper anti-rollover joint control system; 所述上层防侧翻联合控制系统计算主动侧倾力矩和差动制动力矩,并将计算出的主动侧倾力矩作为中层悬架垂向力优化分配系统的输入变量;The upper-layer anti-rollover joint control system calculates the active roll moment and the differential braking moment, and uses the calculated active roll moment as the input variable of the middle-layer suspension vertical force optimization distribution system; 所述中层悬架垂向力优化分配系统实现对于主动悬架执行器的垂向力分配,具体包括如下步骤:The vertical force distribution system for the middle-level suspension realizes the vertical force distribution for the active suspension actuator, which specifically includes the following steps: 设置初始时刻总控制转矩在前桥与中后桥间、左侧轮与右侧轮之间均为各轮处执行器平均分配;Set the total control torque between the front axle and the middle and rear axles, and between the left and right wheels at the initial moment, all of which are equally distributed by the actuators at each wheel; 确定有效集法的设计变量X的约束条件,所述约束条件考虑车轮在良好路面的附着条件、横向载荷转移率要求、垂向动态力对车身的拉压要求、主动侧倾力矩的输出要求、以及主动悬架执行器的输出能力;Determine the constraints of the design variable X of the active set method, which consider the adhesion conditions of the wheels on the good road surface, the lateral load transfer rate requirements, the tension and compression requirements of the vertical dynamic force on the body, and the output requirements of the active roll moment. , and the output capability of the active suspension actuator; 在迭代可行点周围确定所有约束中的有效约束集合;Determine a valid set of constraints among all constraints around an iterative feasible point; 经过多次迭代得到在上述约束下的目标函数的最优解,所述目标函数为各轮胎平均相对动载,得到所述主动悬架执行器的最优垂向力分配比例,计算左前桥、右前桥、左中后桥和右中后桥的主动悬架执行器所需的垂向控制力;After several iterations, the optimal solution of the objective function under the above constraints is obtained, the objective function is the average relative dynamic load of each tire, the optimal vertical force distribution ratio of the active suspension actuator is obtained, and the left front axle, The vertical control force required by the active suspension actuators of the right front axle, left middle rear axle and right middle rear axle; 下层驱动与制动防滑控制系统基于差动制动力矩,以理想滑移率为目标,基于滑模变结构法计算驱动转矩。The lower-level driving and braking anti-skid control system is based on the differential braking torque, aiming at the ideal slip rate, and calculating the driving torque based on the sliding mode variable structure method. 2.如权利要求1所述的方法,其特征在于:所述上层防侧翻联合控制系统以方向盘转角信号作为输入变量,通过车辆三自由度参考模型计算输出变量理想侧倾角和理想横摆角速度。2. The method according to claim 1, characterized in that: the upper-layer anti-rollover joint control system takes the steering wheel angle signal as an input variable, and calculates the output variables ideal roll angle and ideal yaw angular velocity through a vehicle three-degree-of-freedom reference model . 3.如权利要求1所述的方法,其特征在于:基于模糊增益调节的滑模变结构控制算法计算主动侧倾力矩。3 . The method of claim 1 , wherein the active roll moment is calculated based on a sliding mode variable structure control algorithm with fuzzy gain adjustment. 4 . 4.如权利要求1所述的方法,其特征在于:基于滑模变结构控制算法计算差动制动力矩。4. The method of claim 1, wherein the differential braking torque is calculated based on a sliding mode variable structure control algorithm. 5.如权利要求1至4任一项所述的方法,其特征在于:主动悬架执行器的垂向力分配具体为对主动防侧倾力矩在前桥主动悬架执行器与中后桥主动悬架执行器之间分配比例,以及左侧主动悬架执行器与右侧主动悬架执行器之间分配比例的控制。5. The method according to any one of claims 1 to 4, characterized in that: the vertical force distribution of the active suspension actuator is specifically for the active anti-roll moment between the front axle active suspension actuator and the middle and rear. The control of the distribution ratio between the axle active suspension actuators and the distribution ratio between the left active suspension actuator and the right active suspension actuator. 6.如权利要求1至4任一项所述的方法,其特征在于:基于有效集法实现对于主动悬架执行器的垂向力分配。6. The method according to any one of claims 1 to 4, wherein the vertical force distribution for the active suspension actuator is realized based on an active set method. 7.如权利要求1至4任一项所述的方法,其特征在于:所述计算主动侧倾力矩具体包括如下步骤:7. The method according to any one of claims 1 to 4, wherein the calculating the active roll moment specifically comprises the following steps: 计算侧倾角偏差e及其偏差率
Figure FDA0002571026120000021
其计算过程如下:
Calculate the roll angle deviation e and its deviation rate
Figure FDA0002571026120000021
Its calculation process is as follows:
e=φ0e=φ 0
Figure FDA0002571026120000022
Figure FDA0002571026120000022
其中,φ为侧倾角,φ0为理想侧倾角;Among them, φ is the roll angle, and φ 0 is the ideal roll angle; 选用指数趋近律,设计滑模变结构控制器的线性切换函数s;The exponential reaching law is used to design the linear switching function s of the sliding mode variable structure controller;
Figure DEST_PATH_IMAGE001
Figure DEST_PATH_IMAGE001
式中,c为常数,ε为切换系数,k表示滑模切换增益,且c>0,ε>0,k>0;In the formula, c is a constant, ε is the switching coefficient, k is the sliding mode switching gain, and c>0, ε>0, k>0; sgn表示符号函数,其表达式为:sgn represents a symbolic function whose expression is:
Figure FDA0002571026120000024
Figure FDA0002571026120000024
最终的滑模控制器设计为;The final sliding mode controller is designed as;
Figure DEST_PATH_IMAGE002
Figure DEST_PATH_IMAGE002
式中,
Figure FDA0002571026120000026
由基于39DOF车辆模型的上层防侧翻控制系统的控制方程得到:
In the formula,
Figure FDA0002571026120000026
The control equation of the upper-layer anti-rollover control system based on the 39DOF vehicle model is obtained:
Figure FDA0002571026120000027
Figure FDA0002571026120000027
其中,
Figure FDA0002571026120000028
为侧倾角加速度,确定上述切换系数ε采用模糊控制法,具体为模糊控制器中输入变量侧倾角偏差e的论域为[0,15],输出变量切换系数ε的论域为[0,1.5];侧倾角偏差e、切换系数ε的语言值集合都为{ZO,PM,PB},表示自然语言中的{零,正中,正大};所述输入变量和输出变量隶属度函数选用三角-S型混合隶属度函数,解模糊方法选用面积重心法。
in,
Figure FDA0002571026120000028
For the roll angle acceleration, the fuzzy control method is used to determine the above switching coefficient ε. Specifically, the domain of the input variable roll angle deviation e in the fuzzy controller is [0, 15], and the domain of the output variable switching coefficient ε is [0, 1.5. ]; the language value sets of the roll angle deviation e and the switching coefficient ε are all {ZO, PM, PB}, representing {zero, middle, positive} in natural language; the membership functions of the input variables and output variables are selected from triangular- S-shaped mixed membership function, the defuzzification method adopts the area centroid method.
8.如权利要求1至4任一项所述的方法,其特征在于:所述下层驱动与制动防滑控制系统以20%的滑移率对应的旋转角速度为参考目标。8. The method according to any one of claims 1 to 4, wherein the lower-layer driving and braking anti-skid control system takes a rotation angular velocity corresponding to a slip ratio of 20% as a reference target. 9.一种多轴分布式驱动车辆,其特征在于采用如权利要求1-8任一项所述的分层控制的车辆防侧翻方法。9. A multi-axle distributed drive vehicle, characterized in that the layered control vehicle rollover prevention method according to any one of claims 1-8 is adopted.
CN201910884468.4A 2019-09-18 2019-09-18 Layered control vehicle rollover prevention method and multi-shaft distributed driving vehicle Expired - Fee Related CN110606079B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910884468.4A CN110606079B (en) 2019-09-18 2019-09-18 Layered control vehicle rollover prevention method and multi-shaft distributed driving vehicle

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910884468.4A CN110606079B (en) 2019-09-18 2019-09-18 Layered control vehicle rollover prevention method and multi-shaft distributed driving vehicle

Publications (2)

Publication Number Publication Date
CN110606079A CN110606079A (en) 2019-12-24
CN110606079B true CN110606079B (en) 2020-11-17

Family

ID=68892229

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910884468.4A Expired - Fee Related CN110606079B (en) 2019-09-18 2019-09-18 Layered control vehicle rollover prevention method and multi-shaft distributed driving vehicle

Country Status (1)

Country Link
CN (1) CN110606079B (en)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111645696B (en) * 2020-04-28 2021-09-10 武汉理工大学 Method for identifying complex off-road working condition of distributed driving off-road vehicle
CN111731268B (en) * 2020-06-01 2022-09-30 南京航空航天大学 Vehicle rollover prevention control method considering different driver characteristics
CN112172789B (en) * 2020-09-25 2021-07-02 河南坐骑科技有限公司 Anti-rollover control method for active-side-tipping vehicle
CN112172540B (en) * 2020-09-27 2021-12-31 上海工程技术大学 Rollover prevention comprehensive control method for distributed drive electric automobile
CN113147735B (en) * 2021-04-15 2022-05-06 浙江理工大学 Differential braking/driving coordination rollover prevention control system and control method thereof
CN113072012B (en) * 2021-04-25 2022-06-07 合肥工业大学 An anti-rollover control method for a counterbalanced forklift
CN113370798B (en) * 2021-05-27 2022-09-13 武汉理工大学 A differential torsional rollover prevention control method and device for in-wheel motor-driven vehicles
CN114604234B (en) * 2022-04-12 2024-12-27 中国第一汽车股份有限公司 Vehicle rollover prevention control method, system, computer readable storage medium and vehicle
CN114802310B (en) * 2022-06-08 2025-02-18 中国第一汽车股份有限公司 Vehicle trajectory tracking control method, device, electronic device and storage medium
CN115972840B (en) * 2022-09-30 2025-06-24 合肥工业大学 Hydraulic interconnection suspension system and control method thereof
CN116902108B (en) * 2023-09-14 2023-12-05 北京理工大学 High-speed steering active roll control method and system applied to wheel-legged vehicles
CN117485325B (en) * 2024-01-02 2024-03-19 中国重汽集团济南动力有限公司 A multi-axis distributed electric drive vehicle steering control method and vehicle

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8712639B2 (en) * 2006-08-30 2014-04-29 Ford Global Technologies Integrated control system for stability control of yaw, roll and lateral motion of a driving vehicle using an integrated sensing system to determine longitudinal velocity
DE102008014104A1 (en) * 2008-03-13 2008-10-23 Daimler Ag Controlling handling of automobile, e.g. to anticipate adverse weather conditions, adjusts suspension components as function of road friction
US20130060423A1 (en) * 2010-05-14 2013-03-07 Lord Corporation Land vehicles and systems with controllable suspension systems
FR3003834B1 (en) * 2013-03-26 2015-04-24 Univ Blaise Pascal Clermont Ii METHOD OF REDUCING THE RISK OF TILTING OF A MOTOR VEHICLE EQUIPPED WITH A PILOTABLE SUSPENSION SYSTEM
US10543849B2 (en) * 2016-03-31 2020-01-28 Airgo Ip, Llc Vehicle operational diagnostics and condition response system
US10384672B1 (en) * 2016-05-11 2019-08-20 Apple Inc. Vehicle stability control system
CN106394560B (en) * 2016-09-09 2018-08-28 山东理工大学 A kind of wheel hub driving electric automobile chassis coordinated control system
US11110913B2 (en) * 2016-11-18 2021-09-07 Polaris Industries Inc. Vehicle having adjustable suspension
WO2018173303A1 (en) * 2017-03-24 2018-09-27 株式会社ショーワ Control device and suspension device
CN107512262B (en) * 2017-08-14 2018-12-21 吉林大学 A kind of vehicle stability control system tire force distribution method when the driving limited space for execution
CN109649376B (en) * 2017-10-11 2021-01-26 郑州宇通客车股份有限公司 Rollover prevention active intervention control method and system
CN108146431B (en) * 2017-12-22 2020-08-07 天津职业技术师范大学 Automobile rollover stability integrated control method under limit braking condition
CN108891221A (en) * 2018-07-24 2018-11-27 山东大学 A kind of active suspension system and its working method based on mode energy distribution method
CN109606352B (en) * 2018-11-22 2020-06-26 江苏大学 Vehicle path tracking and stability coordination control method
CN109895579A (en) * 2019-03-04 2019-06-18 岭南师范学院 A kind of braking coupling Active suspension anti-side turns over integrated control device and its control method

Also Published As

Publication number Publication date
CN110606079A (en) 2019-12-24

Similar Documents

Publication Publication Date Title
CN110606079B (en) Layered control vehicle rollover prevention method and multi-shaft distributed driving vehicle
CN112644457B (en) A distributed driving vehicle steering stability control system and its control method
CN109733382B (en) Automobile rollover prevention method based on model predictive control
CN113002324B (en) Electronic differential system of four-wheel independent driving and independent steering electric automobile
WO2020187259A1 (en) Safety monitoring method and system for autonomous vehicle, and motion control system
CN110605973B (en) A multi-axis distributed electric vehicle handling stability control method based on layered structure
CN109606133A (en) Torque vector control method for distributed drive electric vehicle based on double-layer control
CN109552312A (en) Intact stability model predictive control method
CN112224036B (en) Four-wheel driving torque distribution method and system for distributed driving electric vehicle
CN111391595A (en) Model Predictive Control Method for Vehicle Rollover Prevention Active Tilting
CN110901630A (en) Method and system for controlling lateral stability of semi-trailer train
CN114771503A (en) Automobile transverse stability control method based on nonsingular terminal sliding mode control
CN114801781A (en) Fuzzy-sliding mode composite control system of four-wheel drive AFS/DYC integrated control system
Liang et al. Integration of active tilting control and full-wheel steering control system on vehicle lateral performance
CN119017884B (en) Vehicle body posture estimation and active control method thereof
CN113449378A (en) Transverse load transfer rate calculation method based on vehicle suspension deformation
CN115431790A (en) Cooperative control method of AFS and DYC for 8-wheel distributed electric drive vehicle
CN115230679A (en) Modeling method of vehicle steering control based on the proportionality of front and rear wheel angle state feedback input
CN114148319B (en) A Rollover Early Warning and Control Method Based on Quadratic Predictive Zero Moment Point
CN114560011B (en) System and method for improving dead time safety of wheel steering clamp of heavy multi-axle vehicle
CN116767201A (en) Commercial vehicle lane keeping control method considering roll stability
CN118928364A (en) Active rollover prevention control method for heavy commercial vehicles
CN117141461A (en) A collaborative control method of active front wheel steering and direct yaw moment based on vehicle road conditions
CN108839656B (en) Method for Determining Driving Moment of Multi-axis Distributed Drive Articulated Bus
CN114194035B (en) Balanced arm type six-wheel independent drive intelligent chassis torque distribution method

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20201117

CF01 Termination of patent right due to non-payment of annual fee