CN109733382B - Automobile rollover prevention method based on model predictive control - Google Patents

Automobile rollover prevention method based on model predictive control Download PDF

Info

Publication number
CN109733382B
CN109733382B CN201811557669.5A CN201811557669A CN109733382B CN 109733382 B CN109733382 B CN 109733382B CN 201811557669 A CN201811557669 A CN 201811557669A CN 109733382 B CN109733382 B CN 109733382B
Authority
CN
China
Prior art keywords
ltr
vehicle
automobile
des
front wheel
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
Application number
CN201811557669.5A
Other languages
Chinese (zh)
Other versions
CN109733382A (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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201811557669.5A priority Critical patent/CN109733382B/en
Publication of CN109733382A publication Critical patent/CN109733382A/en
Application granted granted Critical
Publication of CN109733382B publication Critical patent/CN109733382B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Control Of Driving Devices And Active Controlling Of Vehicle (AREA)
  • Regulating Braking Force (AREA)

Abstract

The invention discloses an automobile rollover prevention method based on model predictive control, which comprises the steps of firstly measuring the front wheel corner, the speed, the yaw angular velocity and the transverse position of an automobile; then calculating the transverse load transfer rate of the vehicle; and then comparing the transverse load transfer rate with a preset transverse load transfer rate threshold value, comparing the difference value between the transverse position and a preset reference transverse position with a preset maximum difference threshold value, adjusting the front wheel rotation angle of the vehicle according to the comparison result of the transverse load transfer rate and the preset reference transverse position, and controlling the braking force of the four wheels. The invention considers the different influences of the sprung mass and the unsprung mass on the rollover of the vehicle body, and has more accurate effect on the rollover prevention control effect.

Description

Automobile rollover prevention method based on model predictive control
Technical Field
The invention relates to the technical field of automobile active safety, in particular to an automobile rollover prevention method based on model predictive control.
Background
With the development of economy, the number of automobiles in the world continues to increase, roads become crowded, and the danger of automobile traveling becomes great. Therefore, people pay attention to the comfort and economy of the automobile and look at the safety. Of all automobile safety accidents, there is a class of accidents that have a low incidence, but the probability of fatal injury is quite low. The survey report of the National Highway Traffic Safety Administration (NHTSA) in the united states shows that the number of deaths caused by rollover accidents in 2014 accounts for one third of the number of deaths caused by all traffic accidents in the same year, and the number of deaths caused by rollover accidents is 7659, which is more than that of other vehicle types.
The existing anti-rollover control methods mainly applied to the anti-rollover field include active steering, differential braking, active/semi-active suspension, and a rollover evaluation method based on LTR of a roll angle and a roll angle speed, which is commonly adopted in the aspect of rollover evaluation indexes. The active/semi-active suspension mainly adjusts the damping by changing a hydraulic hole of a suspension damper so as to strengthen the rollover prevention capability of the vehicle.
However, most of the above control methods only focus on the posture of the vehicle body when preventing the vehicle from rolling over, achieve the anti-rolling effect by controlling the vehicle body, and ignore the control effect generated by the vehicle body, and may cause the vehicle to deviate from the original path or to run into another lane to generate secondary collision or other dangers, which violate the intention and the control intention of the driver, so there are places to be improved.
Disclosure of Invention
The invention aims to solve the technical problem of providing an automobile rollover prevention and tracking system based on model predictive control and control thereof aiming at the defects involved in the background technology.
The invention adopts the following technical scheme for solving the technical problems:
a model predictive control-based automobile rollover prevention method comprises the following steps:
step 1), measuring the front wheel corner delta of the vehiclefVehicle speed vrYaw rate βrA lateral position y;
step 2), calculating the lateral load transfer rate LTR of the vehicle:
step 3), the transverse load transfer rate LTR and a preset transverse load transfer rate threshold value LTR are compareddesComparing the horizontal position Y with a preset reference horizontal position YdesComparing the difference value of η with η, wherein η is a preset maximum difference threshold:
when (Y-Y)des) < η and LTR < LTRdesWhen the vehicle body state is not adjusted;
when (Y-Y)des) | is greater than η and LTR < LTRdesWhen the automobile is driven, two front wheels of the automobile are controlled to steer;
when (Y-Y)des) < η and LTR > LTRdesWhen in use, four wheels of the automobile are controlled to brake,
when (Y-Y)des) I > η and LTR > LTRdesMeanwhile, two front wheels of the automobile are controlled to steer, and four wheels of the automobile are controlled to brake.
As a further optimization scheme of the model predictive control-based automobile rollover prevention method, the step 2) comprises the following detailed steps:
step 2.1), according to ay=vrβrCalculating the lateral acceleration a of the vehicley
Step 2.2), force analysis is carried out on the y direction to obtain the roll angle acceleration
Figure BDA0001912401570000021
Figure BDA0001912401570000022
In the formula, m is the total mass of the automobile; m issIs the sprung mass; e is the distance from the sprung mass to the centre of roll;
Figure BDA0001912401570000023
is the roll angular acceleration; k is a radical of1、k2The yaw stiffness of the front wheel and the rear wheel of the vehicle respectively, the distance from the center of mass of the vehicle to the front axle and the rear axle respectively, the yaw angle of the center of mass β, βrThe yaw angular velocity;
step 2.3), carrying out moment balance on the x direction:
Figure BDA0001912401570000024
in the formula IxsThe moment of inertia of the sprung mass of the automobile around the x axis of the vehicle coordinate system; m issIs the sprung mass; g is the acceleration of gravity and the acceleration of gravity,
Figure BDA0001912401570000025
the vehicle body is in a side inclination angle,
Figure BDA0001912401570000026
in order to be the roll rigidity,
Figure BDA0001912401570000027
in order to damp the roll,
Figure BDA0001912401570000028
is the vehicle body roll angle velocity;
step 2.4), carrying out moment balance on the z direction to obtain the yaw angular acceleration
Figure BDA0001912401570000029
Figure BDA00019124015700000210
Wherein, IzThe moment of inertia of the whole vehicle mass around the z axis of a vehicle coordinate system;
Figure BDA00019124015700000211
yaw angular acceleration;
step 2.5), establishing a three-degree-of-freedom vehicle model, and obtaining the roll angle and the roll angle speed of the vehicle in the current state according to the vehicle speed and the rotation angle;
step 2.5), carrying out stress analysis according to the pressure of the tires on the road surface in the rollover process of the automobile:
Figure BDA00019124015700000212
in the formula, FR、FLThe longitudinal counter-force sum of the right wheel and the left wheel, T is the wheel track, mdUnsprung mass, hdThe distance from the unsprung mass center of mass to the ground;
step 2.6) rate of transfer due to transverse load
Figure BDA0001912401570000031
And FR+FLWhen the weight is mg:
Figure BDA0001912401570000032
as a further optimization scheme of the model predictive control-based automobile rollover prevention method, the concrete steps of controlling two front wheels of the automobile to steer in the step 3) are as follows:
step 3.1.1), the electronic control unit determines the reference transverse position Y from YdesCalculating corresponding additional front wheel corner delta by using model prediction control methodf
Step 3.1.2), establishing a nonlinear model of the active front wheel steering controller:
Figure BDA0001912401570000033
in the formula:
Figure BDA0001912401570000034
is the longitudinal acceleration;
Figure BDA0001912401570000035
is the transverse vehicle speed; c. C1,c2Longitudinal stiffness of the front and rear wheels, respectively; s1,s2Respectively the slip rates of the front wheel and the rear wheel;
Figure BDA0001912401570000036
is the lateral acceleration;
Figure BDA0001912401570000037
is the longitudinal vehicle speed;
Figure BDA0001912401570000038
inertial coordinate system lateral velocity;
Figure BDA0001912401570000039
is the longitudinal velocity of the inertial coordinate system; r is a yaw angle;
step 3.1.3), the equation of state is obtained
Figure BDA00019124015700000310
Discretizing the method:
Figure BDA00019124015700000311
wherein
Figure BDA00019124015700000312
u=δfA, B are coefficients;
step 3.1.4), designing an objective function as follows, and calculating an additional front wheel corner:
Figure BDA00019124015700000313
the front wheel steering angle constraint is satisfied: deltaf,min≤δf≤δf,max
Wherein Y (t) is the lateral position actually output by the system at time t; y isc(t) is the lateral position referenced at time t; n is a radical ofpIs the prediction time domain; n is a radical ofcIs the control time domain; u (t) is a control variable output by the active front wheel steering controller at the time t; deltafFor front wheel rotationAngles Q, R are weight matrices of the first term and the second term of the objective function respectively, p is a weight coefficient, and epsilon is a relaxation factor; deltaf,minAt the smallest front wheel angle, deltaf,maxIs the maximum front wheel turning angle;
step 3.1.5), the controller controls two front wheels of the automobile to steer according to the calculated additional front wheel steering angle;
as a further optimization scheme of the model predictive control-based automobile rollover prevention method, the specific steps of controlling four wheels of an automobile to brake in the step 3) are as follows:
step 3.2.1), the electronic control unit transfers the rate LTR according to the LTR and the reference lateral loaddesThe longitudinal force F of the front left tire, the front right tire, the rear left tire and the rear right tire is calculated by adopting an active braking controllerxfl、Fxfr、Fxrl、FxrrThe nonlinear model of the active braking controller is:
Figure BDA0001912401570000041
in the formula, Fyfl、Fyfr、Fyrl、FyrrThe lateral forces of the front left tire, the front right tire, the rear left tire and the rear right tire are respectively;
step 3.2.2) to obtain the equation of state
Figure BDA0001912401570000042
Discretizing the method:
Figure BDA0001912401570000043
wherein
Figure BDA0001912401570000044
U=[Fxfl,Fxfr,Fxrl,Fxrr]TC. D are all coefficients;
step 3.2.3), designing an objective function as follows, and calculating the braking force of four tires:
Figure BDA0001912401570000045
four tire braking force constraints are satisfied: u shapemin≤U≤Umax
In the formula, Y1(t) LTR actually output by the system at time t; y is1c(t) LTR referenced at time tdes;PhIs the prediction time domain; chIs the control time domain; u (t) is a control variable output by the active braking controller at the time t; s, W are weight matrices for the first and second terms of the objective function, respectively, i is a weight coefficient,
Figure BDA0001912401570000046
is a relaxation factor; u shapeminFor minimum tire force, UmaxMaximum tire force;
and 3.2.4), the controller brakes the four wheels according to the calculated braking force of the four tires.
Compared with the prior art, the invention adopting the technical scheme has the following technical effects:
the invention provides an automobile rollover prevention method based on model predictive control, most of the existing rollover prevention systems adopt active steering or differential braking, and in the working process of a controller, other influences brought by the control effect are rarely considered, such as: drive into other lanes, deviating from the driver's intention. The invention can switch modes according to the driving state by adopting active steering, active braking and switching control modes thereof, thereby realizing perfect integration of rollover prevention and tracking paths;
compared with the conventional LTR acquisition method, the rollover evaluation index provided by the invention is more accurate, and has a more accurate effect on the rollover prevention control effect by considering different influences of the sprung mass and the unsprung mass on the rollover of the vehicle body.
Drawings
FIG. 1 is a block diagram of a patented system of the invention;
FIG. 2 is a schematic flow chart of the present invention;
fig. 3 is a flowchart of steering two front wheels.
Detailed Description
The technical scheme of the invention is further explained in detail by combining the attached drawings:
the present invention may be embodied in many different forms and should not be construed as limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art. In the drawings, components are exaggerated for clarity.
The structure of the invention is shown in figure 1, and comprises a sensor, an electronic control unit and an actuating mechanism which are connected with each other; the electronic control unit comprises a state judgment unit, a state identification unit and an actuating mechanism control unit. The state judging unit and the state identifying unit are both connected with a sensor, the state evaluating unit is connected with the state identifying unit through a triggering unit, and the state identifying unit is connected with the actuating mechanism through the actuating mechanism control unit (model predictive control); and the information of the actuating mechanism unit is fed back to the electronic control unit. The executing mechanism is an active front wheel steering mechanism and an active braking mechanism.
As shown in FIG. 2, the invention discloses a model predictive control-based automobile rollover prevention method, which comprises the following steps:
step 1), measuring the front wheel corner delta of the vehiclefVehicle speed vrYaw rate βrA lateral position y;
step 2), calculating the lateral load transfer rate LTR of the vehicle:
step 2.1), according to ay=vrβrCalculating the lateral acceleration a of the vehicley
Step 2.2), force analysis is carried out on the y direction to obtain the roll angle acceleration
Figure BDA0001912401570000051
Figure BDA0001912401570000052
In the formula, m is the total mass of the automobile; ms is the sprung mass; e is the distance from the sprung mass to the centre of roll;
Figure BDA0001912401570000061
is the roll angular acceleration; k is a radical of1、k2The yaw stiffness of the front wheel and the rear wheel of the vehicle respectively, the distance from the center of mass of the vehicle to the front axle and the rear axle respectively, the yaw angle of the center of mass β, βrThe yaw angular velocity;
step 2.3), carrying out moment balance on the x direction:
Figure BDA0001912401570000062
in the formula IxsThe moment of inertia of the sprung mass of the automobile around the x axis of the vehicle coordinate system; m issIs the sprung mass; g is the acceleration of gravity and the acceleration of gravity,
Figure BDA0001912401570000063
the vehicle body is in a side inclination angle,
Figure BDA0001912401570000064
in order to be the roll rigidity,
Figure BDA0001912401570000065
in order to damp the roll,
Figure BDA0001912401570000066
is the vehicle body roll angle velocity;
step 2.4), carrying out moment balance on the z direction to obtain the yaw angular acceleration
Figure BDA0001912401570000067
Figure BDA0001912401570000068
Wherein, IzFor the whole vehicle mass around the vehicle seatThe moment of inertia of the z axis is marked;
Figure BDA0001912401570000069
yaw angular acceleration;
step 2.5), establishing a three-degree-of-freedom vehicle model, and obtaining the roll angle and the roll angle speed of the vehicle in the current state according to the vehicle speed and the rotation angle;
step 2.5), carrying out stress analysis according to the pressure of the tires on the road surface in the rollover process of the automobile:
Figure BDA00019124015700000610
in the formula, FR、FLThe longitudinal counter-force sum of the right wheel and the left wheel, T is the wheel track, mdUnsprung mass, hdThe distance from the unsprung mass center of mass to the ground;
step 2.6) rate of transfer due to transverse load
Figure BDA00019124015700000611
And FR+FLWhen the weight is mg:
Figure BDA00019124015700000612
step 3), the transverse load transfer rate LTR and a preset transverse load transfer rate threshold value LTR are compareddesComparing the horizontal position Y with a preset reference horizontal position YdesComparing the difference value of η with η, wherein η is a preset maximum difference threshold:
when (Y-Y)des) < η and LTR < LTRdesWhen the vehicle body state is not adjusted;
when (Y-Y)des) | is greater than η and LTR < LTRdesWhen the automobile is driven, two front wheels of the automobile are controlled to steer;
when (Y-Y)des) < η and LTR > LTRdesWhen in use, four wheels of the automobile are controlled to brake,
when (Y-Y)des) I > η and LTR > LTRdesWhen the automobile is in use, two front wheels of the automobile are controlled to steer, and four wheels of the automobile are controlled to brake;
as shown in fig. 3, the specific steps of controlling the two front wheels of the automobile to steer are as follows:
step 3.1.1), the electronic control unit determines the reference transverse position Y from YdesCalculating corresponding additional front wheel corner delta by using model prediction control methodf
Step 3.1.2), establishing a nonlinear model of the active front wheel steering controller:
Figure BDA0001912401570000071
in the formula:
Figure BDA0001912401570000072
is the longitudinal acceleration;
Figure BDA0001912401570000073
is the transverse vehicle speed; c. C1,c2Longitudinal stiffness of the front and rear wheels, respectively; s1,s2Respectively the slip rates of the front wheel and the rear wheel;
Figure BDA0001912401570000074
is the lateral acceleration;
Figure BDA0001912401570000075
is the longitudinal vehicle speed;
Figure BDA0001912401570000076
inertial coordinate system lateral velocity;
Figure BDA0001912401570000077
is the longitudinal velocity of the inertial coordinate system; r is a yaw angle;
step 3.1.3), the equation of state is obtained
Figure BDA0001912401570000078
Discretizing the method:
Figure BDA0001912401570000079
wherein
Figure BDA00019124015700000710
u=δfA, B are coefficients;
step 3.1.4), designing an objective function as follows, and calculating an additional front wheel corner:
Figure BDA00019124015700000711
the front wheel steering angle constraint is satisfied: deltaf,min≤δf≤δf,max
Wherein Y (t) is the lateral position actually output by the system at time t; y isc(t) is the lateral position referenced at time t; n is a radical ofpIs the prediction time domain; n is a radical ofcIs the control time domain; u (t) is a control variable output by the active front wheel steering controller at the time t; deltafQ, R are weight matrixes of a first item and a second item of an objective function respectively for front wheel corners, the function is to enable the system to achieve expected effects as soon as possible, p is a weight coefficient, and epsilon is a relaxation factor, so that the problem that no optimal solution exists in the solving process is solved; deltaf,minAt the smallest front wheel angle, deltafmaxIs the maximum front wheel turning angle;
step 3.1.5), the controller controls two front wheels of the automobile to steer according to the calculated additional front wheel steering angle;
the specific steps of controlling four wheels of the automobile to brake are as follows:
step 3.2.1), the electronic control unit transfers the rate LTR according to the LTR and the reference lateral loaddesThe longitudinal force F of the front left tire, the front right tire, the rear left tire and the rear right tire is calculated by adopting an active braking controllerxfl、Fxfr、Fxrl、FxrrThe nonlinear model of the active braking controller is:
Figure BDA0001912401570000081
in the formula, Fyfl、Fyfr、Fyrl、FyrrThe lateral forces of the front left tire, the front right tire, the rear left tire and the rear right tire are respectively;
step 3.2.2) to obtain the equation of state
Figure BDA0001912401570000082
Discretizing the method:
Figure BDA0001912401570000083
wherein
Figure BDA0001912401570000084
U=[Fxfl,Fxfr,Fxrl,Fxrr]TC, D are coefficients;
step 3.2.3), designing an objective function as follows, and calculating the braking force of four tires:
Figure BDA0001912401570000085
four tire braking force constraints are satisfied: u shapemin≤U≤Umax
In the formula, Y1(t) LTR actually output by the system at time t; y is1c(t) LTR referenced at time tdes;PhIs the prediction time domain; chIs the control time domain; u (t) is a control variable output by the active braking controller at the time t; s, W are weight matrices for the first and second terms of the objective function, respectively, the function of which is to bring the system to the desired effect as soon as possible, i being the weight coefficients,
Figure BDA0001912401570000086
the relaxation factor is used for preventing the occurrence of the situation that no optimal solution exists in the solving process; u shapeminFor minimum tire force, UmaxMaximum tire force;
and 3.2.4), the controller brakes the four wheels according to the calculated braking force of the four tires.
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 will be further understood that terms, such as those defined in commonly used dictionaries, should be interpreted as having a meaning that is consistent with their meaning in the context of the prior art and will not be interpreted in an idealized or overly formal sense unless expressly so defined herein.
The above-mentioned embodiments, objects, technical solutions and advantages of the present invention are further described in detail, it should be understood that the above-mentioned embodiments are only illustrative of the present invention and are not intended to limit the present invention, and any modifications, equivalents, improvements and the like made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims (3)

1. A model predictive control-based automobile rollover prevention method is characterized by comprising the following steps:
step 1.1), measuring the front wheel corner delta of the vehiclefVehicle speed vrYaw rate βrA lateral position y;
step 1.2), calculating the lateral load transfer rate LTR of the vehicle:
step 1.3), the transverse load transfer rate LTR and a preset transverse load transfer rate threshold value LTR are compareddesComparing the horizontal position Y with a preset reference horizontal position YdesComparing the difference value of η with η, wherein η is a preset maximum difference threshold:
when (Y-Y)des) < η and LTR < LTRdesWhen the vehicle body state is not adjusted;
when (Y-Y)des) | is greater than η and LTR < LTRdesWhen the automobile is driven, two front wheels of the automobile are controlled to steer;
when (Y-Y)des) < η and LTR > LTRdesWhen the four wheels of the automobile are controlled to carry outThe brake is carried out by the brake device,
when (Y-Y)des) I > η and LTR > LTRdesWhen the automobile is in use, two front wheels of the automobile are controlled to steer, and four wheels of the automobile are controlled to brake;
step 1.4), the step 1.3) LTR is obtained by the following formula, and the detailed steps are as follows:
step 1.4.1), according to ay=vrβrCalculating the lateral acceleration a of the vehicley
Step 1.4.2), force analysis is carried out on the y direction to obtain the roll angle acceleration
Figure FDA0002406875070000011
Figure FDA0002406875070000012
In the formula, m is the total mass of the automobile; m issIs the sprung mass; e is the distance from the sprung mass to the centre of roll;
Figure FDA0002406875070000013
is the roll angular acceleration; k is a radical of1、k2The yaw stiffness of the front wheel and the rear wheel of the vehicle respectively, the distance from the center of mass of the vehicle to the front axle and the rear axle respectively, the yaw angle of the center of mass β, βrThe yaw angular velocity;
step 1.4.3), performing moment balance on the x direction:
Figure FDA0002406875070000014
in the formula IxsThe moment of inertia of the sprung mass of the automobile around the x axis of the vehicle coordinate system; m issIs the sprung mass; g is the acceleration of gravity and the acceleration of gravity,
Figure FDA0002406875070000015
the vehicle body is in a side inclination angle,
Figure FDA0002406875070000016
in order to be the roll rigidity,
Figure FDA0002406875070000017
in order to damp the roll,
Figure FDA0002406875070000018
is the vehicle body roll angle velocity;
step 1.4.4), carrying out moment balance on the z direction to obtain the yaw angular acceleration
Figure FDA0002406875070000019
Figure FDA00024068750700000110
Wherein, IzThe moment of inertia of the whole vehicle mass around the z axis of a vehicle coordinate system;
Figure FDA00024068750700000111
yaw angular acceleration;
step 1.4.5), establishing a three-degree-of-freedom vehicle model, and obtaining the roll angle and the roll angle speed of the vehicle in the current state according to the vehicle speed and the rotation angle;
step 1.4.6), carrying out stress analysis according to the pressure of tires on the road surface in the rollover process of the automobile:
Figure FDA0002406875070000021
in the formula, FR、FLThe longitudinal counter-force sum of the right wheel and the left wheel, T is the wheel track, mdUnsprung mass, hdThe distance from the center of mass of the unsprung mass to the ground, and h is the height from the center of mass to the ground;
step 1.4.7) due to lateral load transfer rate
Figure FDA0002406875070000022
And FR+FLWhen the weight is mg:
Figure FDA0002406875070000023
2. the model predictive control-based automobile rollover prevention method according to claim 1, wherein the specific steps of controlling two front wheels of the automobile to steer in the step 1.3) are as follows:
step 2.1), the electronic control unit determines the reference transverse position Y from YdesCalculating corresponding additional front wheel corner delta by using model prediction control methodf
Step 2.2), establishing a nonlinear model of the active front wheel steering controller:
Figure FDA0002406875070000024
in the formula:
Figure FDA0002406875070000025
is the longitudinal acceleration;
Figure FDA0002406875070000026
is the transverse vehicle speed; c. C1,c2Longitudinal stiffness of the front and rear wheels, respectively; s1,s2Respectively the slip rates of the front wheel and the rear wheel;
Figure FDA0002406875070000027
is the lateral acceleration;
Figure FDA0002406875070000028
is the longitudinal vehicle speed;
Figure FDA0002406875070000029
inertial coordinate system lateral velocity;
Figure FDA00024068750700000210
is the longitudinal velocity of the inertial coordinate system; r is a yaw angle;
step 2.3), obtaining an equation of state
Figure FDA00024068750700000211
Discretizing the method:
Figure FDA00024068750700000212
wherein
Figure FDA00024068750700000213
u=δfA, B are coefficients, X is the longitudinal position of the inertial frame, and Y is the transverse position of the inertial frame;
step 2.4), designing an objective function as follows, and calculating an additional front wheel corner:
Figure FDA0002406875070000031
the front wheel steering angle constraint is satisfied: deltaf,min≤δf≤δf,max
Wherein Y (t) is the lateral position actually output by the system at time t; y isc(t) is the lateral position referenced at time t; n is a radical ofpIs the prediction time domain; n is a radical ofcIs the control time domain; u (t) is a control variable output by the active front wheel steering controller at the time t; deltafQ, R are weight matrixes of a first term and a second term of an objective function respectively for front wheel turning angles, p is a weight coefficient, and epsilon is a relaxation factor; deltaf,minAt the smallest front wheel angle, deltaf,maxIs the maximum front wheel turning angle;
and 2.5), controlling two front wheels of the automobile to steer by the controller according to the calculated additional front wheel steering angle.
3. The model predictive control-based vehicle rollover prevention method according to claim 1, wherein the specific steps of controlling four wheels of the vehicle to brake in the step 1.3) are as follows:
step 3.1), the electronic control unit transfers the rate LTR according to the LTR and the reference transverse loaddesThe longitudinal force F of the front left tire, the front right tire, the rear left tire and the rear right tire is calculated by adopting an active braking controllerxfl、Fxfr、Fxrl、FxrrThe nonlinear model of the active braking controller is:
Figure FDA0002406875070000032
in the formula, Fyfl、Fyfr、Fyrl、FyrrThe lateral forces of the front left tire, the front right tire, the rear left tire and the rear right tire are respectively;
step 3.2), obtaining an equation of state
Figure FDA0002406875070000033
Discretizing the method:
Figure FDA0002406875070000034
wherein
Figure FDA0002406875070000035
U=[Fxfl,Fxfr,Fxrl,Fxrr]TC, D are coefficients;
step 3.3), designing an objective function as follows, and calculating the braking force of four tires:
Figure FDA0002406875070000036
four tire braking force constraints are satisfied: u shapemin≤U≤Umax
In the formula, Y1(t) LTR actually output by the system at time t; y is1c(t) LTR referenced at time tdes;PhIs the prediction time domain; chIs the control time domain; u (t) isthe control variable output by the active brake controller at the moment t; s, W are weight matrices for the first and second terms of the objective function, respectively, i is a weight coefficient,
Figure FDA0002406875070000041
is a relaxation factor; u shapeminFor minimum tire force, UmaxMaximum tire force;
and 3.4), the controller brakes the four wheels according to the calculated braking forces of the four tires.
CN201811557669.5A 2018-12-19 2018-12-19 Automobile rollover prevention method based on model predictive control Active CN109733382B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811557669.5A CN109733382B (en) 2018-12-19 2018-12-19 Automobile rollover prevention method based on model predictive control

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811557669.5A CN109733382B (en) 2018-12-19 2018-12-19 Automobile rollover prevention method based on model predictive control

Publications (2)

Publication Number Publication Date
CN109733382A CN109733382A (en) 2019-05-10
CN109733382B true CN109733382B (en) 2020-07-07

Family

ID=66360759

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811557669.5A Active CN109733382B (en) 2018-12-19 2018-12-19 Automobile rollover prevention method based on model predictive control

Country Status (1)

Country Link
CN (1) CN109733382B (en)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110126915B (en) * 2019-05-17 2020-11-17 中国农业大学 Active steering control-based active rollover prevention control method and system for wheeled tractor
CN111027132B (en) * 2019-08-08 2022-08-23 中国第一汽车股份有限公司 Vehicle roll control method, device, equipment and storage medium
CN110626353B (en) * 2019-09-09 2020-09-08 武汉理工大学 Vehicle dangerous state early warning method based on roll risk index
CN110562243B (en) * 2019-09-11 2021-07-27 武汉理工大学 Automobile rollover prevention method and device and rollover prevention automobile
CN110471428B (en) * 2019-09-18 2021-05-07 吉林大学 Path tracking method based on variable pre-aiming distance and speed constraint of model
CN111674387B (en) * 2020-05-27 2022-08-05 南京航空航天大学 Method for generating novel rollover early warning index based on derivative iterative prediction
CN111873991B (en) * 2020-07-22 2022-04-08 中国第一汽车股份有限公司 Vehicle steering control method, device, terminal and storage medium
CN112406854B (en) * 2020-11-25 2022-05-10 东风越野车有限公司 Method for controlling side-tipping stability of wheel hub motor-driven off-road vehicle
CN113370798B (en) * 2021-05-27 2022-09-13 武汉理工大学 Differential torsion rollover prevention control method and device for hub motor driven vehicle
CN113806958A (en) * 2021-09-26 2021-12-17 上汽通用五菱汽车股份有限公司 Anti-roll control method, device and storage medium based on MPC algorithm
CN113879282A (en) * 2021-11-09 2022-01-04 杭州云栖智能汽车创新中心 Automatic-driving vehicle rollover prevention control method
CN114890345B (en) * 2022-04-26 2024-05-28 安徽合力股份有限公司 Forklift stability control method and control system

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3736999B2 (en) * 1999-09-06 2006-01-18 本田技研工業株式会社 Vehicle rollover judgment method
EP1873020A1 (en) * 2006-06-29 2008-01-02 Delphi Technologies, Inc. Integrated vehicle crash sensing system and method
CN106080553A (en) * 2016-07-13 2016-11-09 南京航空航天大学 A kind of four-wheel steering automobile anti-rollover control system merging speed change and method
CN106427957A (en) * 2015-08-11 2017-02-22 比亚迪股份有限公司 Stabilization control system and stabilization control method based on four-wheel drive for electric vehicle, as well as electric vehicle
CN107176216A (en) * 2017-06-15 2017-09-19 石家庄铁道大学 Heavy-duty car anti-rollover system
CN107571706A (en) * 2017-09-22 2018-01-12 合肥工业大学 A kind of fork truck anti-rollover control method
CN107745709A (en) * 2017-09-26 2018-03-02 湖北文理学院 Preventing vehicle rollover pre-warning and control method, system and hardware-in-loop simulation method
CN108680364A (en) * 2018-03-29 2018-10-19 南京航空航天大学 A kind of vehicle side turning evaluation index and evaluation method
CN108773376A (en) * 2018-05-07 2018-11-09 南京航空航天大学 A kind of the automobile objective layered Collaborative Control and optimization method of fusion driving intention

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3736999B2 (en) * 1999-09-06 2006-01-18 本田技研工業株式会社 Vehicle rollover judgment method
EP1873020A1 (en) * 2006-06-29 2008-01-02 Delphi Technologies, Inc. Integrated vehicle crash sensing system and method
CN106427957A (en) * 2015-08-11 2017-02-22 比亚迪股份有限公司 Stabilization control system and stabilization control method based on four-wheel drive for electric vehicle, as well as electric vehicle
CN106080553A (en) * 2016-07-13 2016-11-09 南京航空航天大学 A kind of four-wheel steering automobile anti-rollover control system merging speed change and method
CN107176216A (en) * 2017-06-15 2017-09-19 石家庄铁道大学 Heavy-duty car anti-rollover system
CN107571706A (en) * 2017-09-22 2018-01-12 合肥工业大学 A kind of fork truck anti-rollover control method
CN107745709A (en) * 2017-09-26 2018-03-02 湖北文理学院 Preventing vehicle rollover pre-warning and control method, system and hardware-in-loop simulation method
CN108680364A (en) * 2018-03-29 2018-10-19 南京航空航天大学 A kind of vehicle side turning evaluation index and evaluation method
CN108773376A (en) * 2018-05-07 2018-11-09 南京航空航天大学 A kind of the automobile objective layered Collaborative Control and optimization method of fusion driving intention

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
融合助力转向功能的新型主动转向***LQG控制策略;赵万忠等;《中国机械工程》;20140215;第25卷(第3期);417-421页 *

Also Published As

Publication number Publication date
CN109733382A (en) 2019-05-10

Similar Documents

Publication Publication Date Title
CN109733382B (en) Automobile rollover prevention method based on model predictive control
CN106585625B (en) Four-wheel steering vehicle rollover prevention system and control method thereof
CN106945670B (en) Anti-rollover system for automobiles and control method based on driver&#39;s input prediction
Odhams et al. Active steering of a tractor–semi-trailer
US8718872B2 (en) Vehicle attitude controller
CN101288081B (en) Vehicular yaw stabilizing system
US9636965B2 (en) Suspension system
US6684140B2 (en) System for sensing vehicle global and relative attitudes using suspension height sensors
US7647148B2 (en) Roll stability control system for an automotive vehicle using coordinated control of anti-roll bar and brakes
US7063334B2 (en) Vehicle stability system using active tilting mechanism
US8855856B2 (en) Vehicle roll control method using controllable friction force of MR dampers
CN110606079B (en) Layered control vehicle rollover prevention method and multi-shaft distributed driving vehicle
JP4882848B2 (en) Integrated vehicle behavior controller
GB2414815A (en) Determining lateral velocity and yaw rate of a vehicle
US10059344B2 (en) Vehicle behavior control apparatus
GB2425189A (en) Method and system of controlling a vehicle system
CN101184641A (en) Vehicle damping force control apparatus
CN107662468B (en) Vehicle roll motion safety H for active suspension2/H∞Controller design method
CN109591537A (en) A kind of automotive semi-active suspension control system and method
EP1391330B1 (en) An anti-roll or anti-yaw suspension device for vehicles
Sathishkumar et al. Trajectory control for tire burst vehicle using the standalone and roll interconnected active suspensions with safety-comfort control strategy
US6298294B1 (en) Device for controlling suspension shock absorbers of vehicles based upon phantom substitute therefor
US7085639B2 (en) System and method for characterizing the road bank for vehicle roll stability control
Liang et al. Integration of active tilting control and full-wheel steering control system on vehicle lateral performance
JP4515240B2 (en) Control device for suspension system for vehicle

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