CN111695197A - Highly-reliable dynamic estimation method for rollover threshold value of tank car - Google Patents

Highly-reliable dynamic estimation method for rollover threshold value of tank car Download PDF

Info

Publication number
CN111695197A
CN111695197A CN202010444574.3A CN202010444574A CN111695197A CN 111695197 A CN111695197 A CN 111695197A CN 202010444574 A CN202010444574 A CN 202010444574A CN 111695197 A CN111695197 A CN 111695197A
Authority
CN
China
Prior art keywords
rollover
threshold value
test
tank car
air bag
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.)
Granted
Application number
CN202010444574.3A
Other languages
Chinese (zh)
Other versions
CN111695197B (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.)
Southeast University
Original Assignee
Southeast University
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 Southeast University filed Critical Southeast University
Priority to CN202010444574.3A priority Critical patent/CN111695197B/en
Publication of CN111695197A publication Critical patent/CN111695197A/en
Application granted granted Critical
Publication of CN111695197B publication Critical patent/CN111695197B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M17/00Testing of vehicles
    • G01M17/007Wheeled or endless-tracked vehicles
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Abstract

The invention provides a highly reliable dynamic estimation method for a tank car rollover threshold value. The method provides a rollover characterization parameter of the lateral transfer rate of the air bag pressure, and meets the requirement of rollover risk redundancy judgment; the tank car rollover threshold value calibration device is built by utilizing the high-precision inertia measurement unit, the pressure sensor, the wheel force sensor and the rollover prevention frame, and can be applied to a real-vehicle test; and fitting a functional relation between the rollover threshold and the influence factors thereof through real vehicle test data and the SVR to realize dynamic estimation of the rollover threshold.

Description

Highly-reliable dynamic estimation method for rollover threshold value of tank car
Technical Field
The invention relates to a rollover threshold value estimation method, in particular to a highly reliable dynamic estimation method for a rollover threshold value of a tank car, and belongs to the technical field of vehicle safety.
Background
The dangerous goods transportation safety situation of our country is severe, and the rollover is one of the main accident patterns. Due to the large loading capacity and high transportation efficiency, the tank truck is a main carrier for road transportation of dangerous goods. However, the whole tank truck has large mass and high mass center and is disturbed by liquid, so that the tank truck is easy to turn on one side when turning and changing lanes, dangerous goods are leaked, burnt and exploded, and property loss, environmental pollution, ecological damage and casualties are caused. Therefore, the research on the tank car rollover prevention and control method has important significance on road safety.
At present, a rollover prevention and control method is based on a fixed rollover threshold value, such as lateral acceleration of 0.4g, and the like, however, the rollover threshold value of a tank truck is dynamically changed under different driving behaviors, different liquid filling ratios and different road slopes, so that the single fixed rollover threshold value is difficult to adapt to the actual rollover prevention and control requirements of the tank truck. According to data research, the acquisition way of the tank car rollover threshold value is mainly based on a multidimensional simulation test, and the reasons are as follows: 1. the method is characterized by lacking of a highly reliable calibration device for the tank car rollover threshold value, 2 lacking of a typical rollover scene library suitable for the tank car, and 3 lacking of a method capable of dynamically estimating the rollover threshold value according to the tank car running condition.
Disclosure of Invention
The invention provides a highly reliable dynamic estimation method for a tank car rollover threshold value, aiming at the problem that a single fixed rollover threshold value is difficult to meet the actual rollover prevention and control requirements of a tank car. The method can dynamically estimate the rollover threshold value according to the current driving state of the tank car, and is beneficial to improving the accuracy of rollover early warning.
In order to achieve the above purpose, the invention provides the following technical scheme:
the method comprises the following steps: side turning characteristic parameters and factors influencing a side turning threshold are determined and a side turning threshold calibration device is built
The rollover characterizing parameters are selected as a side inclination angle tau, a lateral acceleration theta and an air bag pressure transverse transfer rate eta, factors influencing the rollover threshold value of the tank car are selected as a vehicle speed v, a whole vehicle mass m and a liquid filling ratio lambda, and the eta calculation formula is as follows:
Figure BDA0002505277430000021
in the formula (1), FlsIs the pressure of the left air bag of the s-th axle of the tank car, FrsThe pressure of an air bag at the right side of the s th axle of the tank car, s is the number of the axles, and s is 1,2, …, e and e are the total number of the axles of the tank car;
the rollover threshold value calibration device comprises a high-precision inertia measurement unit, a plurality of pressure sensors, two wheel force sensors and two rollover prevention frames, wherein the high-precision inertia measurement unit is arranged close to the center of mass of the tank car; the high-precision inertial measurement unit acquires a roll angle tau and a lateral acceleration theta, and the wheel force sensor acquires a vertical force G of a left wheel of a final shaftlAnd right wheel vertical force GrAcquiring the average value of information of two wheel speed sensors on a last shaft through a vehicle body CAN bus as a vehicle speed v, wherein the data output frequencies of the sensors are the same, and the vehicle mass m and the liquid filling ratio lambda are obtained by static measurement in advance;
step two: tank car rollover threshold calibration test under typical rollover scene carried out in closed test field
The division of each scene element of the tank car rollover scene library is shown as the following table:
Figure BDA0002505277430000022
after the scene elements are arranged and combined, 30 rollover scenes exist; the tank truck rollover threshold calibration test under a typical rollover scene specifically comprises the following steps:
substep 1: the calibration test is carried out on a dry and solid road surface, and the peak value adhesion coefficient of the road surface is not less than 0.9; checking the reliability of each sensor in the calibration device, the safety of the anti-rollover frame and whether potential safety hazards exist in a checking test field;
substep 2: sequentially setting driving behaviors, liquid filling ratios and longitudinal gradients in a closed test field based on the determined rollover scene, drawing test tracks of J steering and double shifting lines on the test field, and statically measuring the mass m and the liquid filling ratio lambda of the whole vehicle before the test is started;
and substep 3, keeping constant vehicle speed to drive according to a track, gradually increasing for 2km/h by taking 32km/h as an initial speed until the vertical force of the tire on the side of the final shaft is 0 or the side of the rollover prevention support lands in the continuous 5-time test process to finish the rollover threshold calibration test under the current scene, and storing the 5-time test roll angle tau, the lateral acceleration theta, the transverse air bag pressure transfer rate η, the vehicle speed v, the vehicle mass m, the liquid filling ratio lambda and the vertical force G of the left wheel of the final shaftlLast right wheel vertical force GrThe data of (a);
substep 4: repeating the substep 2 and the substep 3, completing threshold value calibration tests under 30 rollover scenes and storing data;
step three: processing data and dynamically calibrating rollover thresholds in different typical rollover scenarios
After 30 rollover scene calibration tests are completed, 5 times of test data are stored in each scene, and 150 times of test data are totally stored, and the roll angle tau, the lateral acceleration theta, the air bag pressure transverse transfer rate η, the vehicle speed v and the vertical force G of the left wheel of the final axle, which are acquired in each test, are filtered by adopting a weighted average valuelLast right wheel vertical force GrProcessing the data; take lateral acceleration as an example, θtFor a set of data of lateral acceleration in a certain test, t is 1,2, …, c, c is the data quantity of the lateral acceleration in the test, and the weighted mean filtering processing method is as follows:
when t is 1 and t is c,
Figure BDA0002505277430000031
when t is 2 and t is c-1,
Figure BDA0002505277430000032
when t is more than or equal to 3 and less than or equal to c-2,
Figure BDA0002505277430000033
the roll angle tau, the air bag pressure transverse transfer rate η, the vehicle speed v and the vertical force G of the left wheel of the final axle of the testlLast right wheel vertical force GrFiltering the data by the same method;
after the data of each item of 150 tests are processed by the method, the roll state H is calculated, and the formula is as follows:
Figure BDA0002505277430000034
in the formula (2), the reaction mixture is,
Figure BDA0002505277430000035
the vertical force of the left wheel of the final axle after the weighted average filtering processing,
Figure BDA0002505277430000036
the method comprises the following steps of calibrating a rollover threshold value for the vertical force of the right wheel of the last axle after weighted mean filtering treatment by using the side-tipping state H of the tank car, and specifically comprises the following steps:
substep 1: uniformly displaying the data of the roll angle, the lateral acceleration, the transverse transfer rate of the air bag pressure, the vehicle speed and the roll state of a certain test on a time axis, namely the curves of all the data have the same starting time and ending time; the whole vehicle mass and the liquid filling ratio in the test are constant values;
substep 2: searching a point at which the value H of the side-rolling state reaches 0.9 for the first time from the initial moment, recording the values of the vehicle speed, the whole vehicle mass, the liquid filling ratio, the side-rolling angle, the lateral acceleration and the air bag pressure transverse transfer rate at the moment, and recording the values as
Figure BDA0002505277430000041
Wherein
Figure BDA0002505277430000042
The rollover threshold for this trial;
substep 3: repeating the above steps, recording the first arrival of the roll state H value in 150 testsThe values of the vehicle speed, the vehicle mass, the liquid filling ratio, the roll angle, the lateral acceleration and the lateral air bag pressure transfer rate at 0.9 are recorded
Figure BDA0002505277430000043
Step four: fitting a functional relation between the rollover threshold and the influence factors thereof by using the SVR
Binding sample data
Figure BDA0002505277430000044
Respectively fitting functions f of the roll angle threshold value, the vehicle speed, the whole vehicle mass and the liquid filling ratio by using SVRτ(x) Function f of lateral acceleration threshold, vehicle speed, vehicle mass and charge ratioθ(x) Function f of lateral transfer rate threshold of air bag pressure and vehicle speed, vehicle mass and liquid filling ratioη(x);
Step five: tank car rollover threshold dynamic estimation based on fitted functional relation
When the tank car runs, the speed of the car is read through the CAN bus of the car body
Figure BDA0002505277430000045
Static measurement of vehicle mass in advance
Figure BDA0002505277430000046
And liquid filling ratio
Figure BDA0002505277430000047
To obtain
Figure BDA0002505277430000048
Respectively calculate
Figure BDA0002505277430000049
And
Figure BDA00025052774300000410
Figure BDA00025052774300000411
is a roll angle control threshold value that is,
Figure BDA00025052774300000412
is a lateral acceleration control threshold value and is,
Figure BDA00025052774300000413
controlling a threshold value for the lateral transfer rate of the air bag pressure;
Figure BDA00025052774300000414
is a roll angle pre-warning threshold value,
Figure BDA00025052774300000415
a threshold value is pre-warned for the lateral acceleration,
Figure BDA00025052774300000416
and (4) early warning threshold value for the transverse transfer rate of the air bag pressure.
Compared with the prior art, the invention has the following advantages and beneficial effects:
1. the invention provides a rollover characterization parameter of the lateral transfer rate of the air bag pressure, meets the requirement of redundant judgment of rollover danger, and improves the reliability of rollover prevention and control;
2. the rollover threshold value calibration and calibration device can accurately acquire data of various characterization parameters of the tank car when rollover occurs in real time, and can be applied to real-vehicle tests;
3. the rollover scene library summarized by the invention covers typical scenes of rollover when the tank truck actually runs;
4. according to the invention, the functional relation between the rollover threshold and the influence factors thereof is fitted through real vehicle test data and the SVR, so that the dynamic estimation of the rollover threshold is realized.
Drawings
FIG. 1 is a general design scheme diagram of tank truck rollover threshold value calibration method and device
FIG. 2 is a structural diagram of a tank car rollover threshold value calibration device
FIG. 3 is a left-hand steering test trace diagram of a tank truck J with a radius of 45.7m
FIG. 4 is a track diagram of a double-line-shifting test of a tank car
Detailed Description
The technical solutions provided by the present invention will be described in detail below with reference to specific examples, and it should be understood that the following detailed description is only illustrative and not intended to limit the scope of the present invention.
The invention provides a highly reliable dynamic estimation method for a tank car rollover threshold value, which comprises the steps of firstly determining rollover characterization parameters and factors influencing the rollover threshold value, building a rollover threshold value calibration device, then carrying out tank car rollover threshold value calibration tests under typical rollover scenes in a closed test field, then processing data and dynamically calibrating rollover threshold values under different typical rollover scenes, further fitting a functional relation between the rollover threshold value and the influencing factors by utilizing SVR (singular value representation) and finally realizing dynamic estimation of the tank car rollover threshold value based on the fitted functional relation. The invention provides a rollover characterization parameter of the lateral transfer rate of the air bag pressure for a tank car using an air bag separated air suspension, and meets the requirement of redundant judgment of rollover dangers; the tank car rollover threshold value calibration device is built by utilizing the high-precision inertia measurement unit, the pressure sensor, the wheel force sensor and the rollover prevention frame, and can be applied to a real-vehicle test; summarizing a typical rollover scene library suitable for the tank car by referring to a plurality of vehicle stability test standards in combination with an actual rollover scene of the tank car; and fitting a functional relation between the rollover threshold and the influence factors thereof through real vehicle test data and the SVR to realize dynamic estimation of the rollover threshold. The general design scheme of the invention is shown in figure 1, and the specific steps comprise:
the method comprises the following steps: side turning characteristic parameters and factors influencing a side turning threshold are determined and a side turning threshold calibration device is built
The roll angle, the lateral acceleration and the transverse load transfer rate are common rollover characterizing parameters and can visually reflect the roll stability of the vehicle. However, the load measurement of the wheels is difficult in a high-speed operation state, and for a tank car using an air bag separated air suspension, because the corresponding relation exists between the air bag pressure at the wheels and the load of the wheels and the air bag distribution is bilaterally symmetrical, a rollover characteristic parameter of the air bag pressure transverse transfer rate is provided, so that the rollover characteristic parameters are selected as a roll angle tau, a lateral acceleration theta and an air bag pressure transverse transfer rate eta. The η calculation formula is:
Figure BDA0002505277430000061
in the formula (1), FlsIs the pressure of the left air bag of the s-th axle of the tank car, FrsThe air bag pressure on the right side of the second axle of the tank truck, s is the axle number, and s is 1,2, …, e and e are the total axle number of the tank truck.
When the curvature radius of the road is fixed, the larger the vehicle speed and the whole vehicle mass are, the larger the centrifugal force of the tank car is, and when the moment of the self weight of the vehicle on the wheels is not enough to overcome the centrifugal force, the vehicle can turn over. Under the non-full load state of the tank car, because of the continuous change of the motion state of the tank car and the liquidity of liquid, the liquid in the tank body is easy to shake, and additional force and moment can be generated on the side wall of the tank body, so that the side-tipping stability of the tank car is reduced, and the side tipping of the vehicle is induced. Therefore, the factors influencing the rollover threshold value of the tank car are selected as the vehicle speed v, the mass m of the whole vehicle and the liquid filling ratio lambda.
The side-turning threshold calibration device comprises a high-precision inertia measurement unit, a plurality of pressure sensors and two wheel force sensors (the introduction and the function of the wheel force sensors are shown in a reference document- -Yanhuawen, a Bluetooth-based wheel force data transmission system design [ D ]]Zhenjiang, university of Jiangsu science and technology, 2013) and two anti-rollover frames. The high-precision inertia measurement unit is installed at a position close to the center of mass of the tank car, the pressure sensors are installed at the vent valves of all air bags of the air suspension of the tank car, the wheel force sensors are installed on wheels on two sides of the last shaft of the tail of the tank car, and the anti-rollover frame is installed on two sides of the tank car as shown in figure 2 (the two-shaft tank car is shown in the figure, and the number of the pressure sensors is increased for the three-shaft tank. The high-precision inertial measurement unit acquires a roll angle tau and a lateral acceleration theta, and the wheel force sensor acquires a vertical force G of a left wheel of a final shaftlAnd the vertical force G of the right wheel of the final axler. As more and more tank cars are equipped with electronic systems such as an anti-lock braking system (ABS), etc., wheel speed sensors have been installed in vehicles, and the average of information of wheel speed sensors of two wheels (i.e., non-steered wheels) on the last axle is collected as a vehicle speed v through a vehicle body CAN bus. The mass m and the liquid filling ratio lambda of the whole vehicle are obtained by static measurement in advance,the data output frequency of each sensor is the same, and the collected data are ensured to be in one-to-one correspondence.
Step two: tank car rollover threshold calibration test under typical rollover scene carried out in closed test field
After determining the rollover characterization parameters, the factors influencing the rollover threshold value and the calibration device, designing a tank car typical rollover scene library, and developing a rollover threshold value calibration test in a closed test field.
The tank truck rollover scene library needs to consider three scene elements of driving behavior, liquid filling ratio and longitudinal gradient (downhill), and each scene element is divided as shown in the following table. Selecting one quantitative parameter from the driving behavior element, the liquid filling ratio element and the longitudinal gradient element, and carrying out permutation and combination to obtain different rollover scenes. To distinguish from sideslip, the low adhesion coefficient is not considered for the moment.
Figure BDA0002505277430000071
The J-steering test trajectory setting refers to steering test regulations in GB/T6323-2014 automobile steering stability test method, as shown in FIG. 3. The double-lane-changing test trajectory is set according to the requirements specified in ISO 3888-2, part 2, namely obstacle avoidance, of a passenger car, a test lane for lane abrupt change operation, and is shown in FIG. 4. When the scene elements are combined in an arrangement mode, in consideration of test safety, the J-turn test is not supported by the longitudinal steep slope of < 7%, so that 30 rollover scenes are available, namely 4 (J-turn) × 5 (liquid filling ratio) × 1 (longitudinal gradient < 3%) +1 (double shift line) × 5 (liquid filling ratio) × 2 (longitudinal gradient).
Through the real vehicle test, the numerical changes of rollover representation parameters and factors influencing the rollover threshold value of the tank car in the process from safe driving to rollover occurrence under different rollover scenes are recorded. The method comprises the following specific steps:
substep 1: the calibration test is carried out on a dry and solid road surface, and the peak value adhesion coefficient of the road surface is not less than 0.9; the reliability of each sensor in the calibration device, the safety of the anti-rollover frame and the potential safety hazard of a test site need to be checked;
substep 2: sequentially setting driving behaviors, liquid filling ratios and longitudinal gradients in a closed test field based on a determined rollover scene, drawing test tracks of J steering and double shifting lines on the test field by using striking colors, and statically measuring the mass m and the liquid filling ratio lambda of the whole vehicle before the test starts;
and substep 3, referring to the speed setting in a J steering test of GB/T38185 plus 2019 electronic stability control system performance requirements and test method for commercial vehicles, keeping a constant vehicle speed as much as possible, driving the vehicle according to a track, gradually increasing by 2km/h by taking 32km/h as an initial speed until the vertical force of the tire on the last shaft side is 0 in the continuous 5-time test process or the side of the rollover prevention support lands to finish the rollover threshold calibration test under the current scene, and storing the 5-time test roll angle tau, the lateral acceleration theta, the lateral transfer rate of the air bag pressure η, the vehicle speed v, the vehicle mass m, the liquid filling ratio lambda and the vertical force G of the left wheel on the last shaftlLast right wheel vertical force GrThe data of (a);
substep 4: repeating the substep 2 and the substep 3, completing threshold value calibration tests under 30 rollover scenes and storing data;
step three: processing data and dynamically calibrating rollover thresholds in different typical rollover scenarios
In order to further improve the accuracy and reliability of the data, the roll inclination angle tau, the lateral acceleration theta, the air bag pressure transverse transfer rate η, the vehicle speed v and the vertical force G of the left wheel of the final axle collected by each test by adopting the weighted average filtering practical in engineeringlLast right wheel vertical force GrThe data of (A) is processed (the mass m and the liquid filling ratio lambda of the whole vehicle in a single test are constant values, and the data do not need to be processed). Take lateral acceleration as an example, θtFor a set of data of lateral acceleration in a certain test, t is 1,2, …, c, c is the data quantity of the lateral acceleration in the test, and the weighted mean filtering processing method is as follows:
when t is 1 and t is c,
Figure BDA0002505277430000081
when t is 2 and t is c-1,
Figure BDA0002505277430000082
when t is more than or equal to 3 and less than or equal to c-2,
Figure BDA0002505277430000083
the roll angle tau, the air bag pressure transverse transfer rate η, the vehicle speed v and the vertical force G of the left wheel of the final axle of the testlLast right wheel vertical force GrThe data is filtered in the same way.
After the data of each item of 150 tests are processed by the method, the roll state H is calculated, and the formula is as follows:
Figure BDA0002505277430000084
in the formula (2), the reaction mixture is,
Figure BDA0002505277430000085
the vertical force of the left wheel of the final axle after the weighted average filtering processing,
Figure BDA0002505277430000086
the vertical force of the right wheel of the final axle after the weighted mean filtering processing. The method comprises the following steps of calibrating a rollover threshold value by using a tank car side-tipping state H:
substep 1: uniformly displaying the data of the roll angle, the lateral acceleration, the transverse transfer rate of the air bag pressure, the vehicle speed and the roll state of a certain test on a time axis, namely the curves of all the data have the same starting time and ending time; the whole vehicle mass and the liquid filling ratio in the test are constant values;
substep 2: searching a point at which the value H of the roll state reaches 0.9 (the value can be changed to be 0.9 according to different requirements of the roll control degree) for the first time from the initial moment, recording the values of the vehicle speed, the whole vehicle mass, the liquid filling ratio, the roll angle, the lateral acceleration and the air bag pressure transverse transfer rate at the moment, and recording the values as
Figure BDA0002505277430000091
Wherein
Figure BDA0002505277430000092
The rollover threshold for this trial;
substep 3: repeating the steps, recording the values of the vehicle speed, the whole vehicle mass, the liquid filling ratio, the roll angle, the lateral acceleration and the air bag pressure transverse transfer rate when the value H of the roll state in 150 tests reaches 0.9 for the first time, and recording the values as the values
Figure BDA0002505277430000093
Step four: fitting a functional relation between the rollover threshold and the influence factors thereof by using the SVR
The support Vector machine is developed for general estimation and prediction problems, and can be used for performing function regression fitting, namely SVR (support Vector regression), besides using the support Vector machine classification. SVR has the advantages of low computational complexity and good generalization performance, by minimizing the prediction error, finding a function that can well approximate the training examples, and maximizing the flatness of the function when the error is minimized reduces the risk of fitting.
Respectively fitting functions f of the roll angle threshold value, the vehicle speed, the whole vehicle mass and the liquid filling ratio by using SVRτ(x) Function f of lateral acceleration threshold, vehicle speed, vehicle mass and charge ratioθ(x) Function f of lateral transfer rate threshold of air bag pressure and vehicle speed, vehicle mass and liquid filling ratioη(x) With fτ(x) An example is the SVR fitting function method.
For the training set (x)i,yi),
Figure BDA0002505277430000094
The expression pattern of SVR is:
Figure BDA0002505277430000095
in equation (3), w is an adjustable weight vector;
Figure BDA0002505277430000096
representing a non-linear mapping, mapping the input quantities to a higher dimensional feature space; b represents a bias; the superscript T denotes transposing the matrix.
Introducing relaxation variables under the condition of controlling fitting accuracyiAnd
Figure BDA0002505277430000097
fτ(x) The objective function solved by the expression is:
Figure BDA0002505277430000101
in the formula (4), C is a regularization constant. The Lagrange multiplier mu is introduced to be more than or equal to 0,
Figure BDA0002505277430000102
α≥0,
Figure BDA0002505277430000103
obtained by the lagrange multiplier method:
Figure BDA0002505277430000104
order to
Figure BDA0002505277430000105
To w, b,iAnd
Figure BDA0002505277430000106
has a partial derivative of zero:
Figure BDA0002505277430000107
bringing formula (6) into formula (5) yields the dual problem of SVR:
Figure BDA0002505277430000108
obtained by the formula (7)
Figure BDA0002505277430000109
According to
Figure BDA00025052774300001010
To obtain w*Finally, find fτ(x) Comprises the following steps:
Figure BDA00025052774300001011
the specific use of SVR and the calculation of b can be referenced (ZhouZhihua. machine learning [ M)]Beijing, Qinghua university Press 2016: 133-137). Finding f by the same methodθ(x) And fη(x)。
Step five: tank car rollover threshold dynamic estimation based on fitted functional relation
When the tank car runs, the speed of the car is read through the CAN bus of the car body
Figure BDA0002505277430000111
Static measurement of vehicle mass in advance
Figure BDA0002505277430000112
And liquid filling ratio
Figure BDA0002505277430000113
To obtain
Figure BDA0002505277430000114
Respectively calculate
Figure BDA0002505277430000115
And
Figure BDA0002505277430000116
the rollover prevention and control part comprises an early warning part and a control part, the rollover early warning part is used for reminding a driver of safe driving when the vehicle has small rollover danger, and the rollover control part is used for controlling the vehicle to have large rollover dangerAnd executing the operation to prevent the rollover. Because the rollover early warning and rollover control are applied to rollover scenes with different risk degrees, the rollover threshold value is divided into an early warning threshold value and a control threshold value. Thus, it is possible to provide
Figure BDA0002505277430000117
Is a roll angle control threshold value that is,
Figure BDA0002505277430000118
is a lateral acceleration control threshold value and is,
Figure BDA0002505277430000119
controlling a threshold value for the lateral transfer rate of the air bag pressure;
Figure BDA00025052774300001110
is a roll angle pre-warning threshold value,
Figure BDA00025052774300001111
a threshold value is pre-warned for the lateral acceleration,
Figure BDA00025052774300001112
and (4) early warning threshold value for the transverse transfer rate of the air bag pressure.

Claims (1)

1. A highly reliable dynamic estimation method for a tank car rollover threshold value is characterized by comprising the following specific steps:
the method comprises the following steps: side turning characteristic parameters and factors influencing a side turning threshold are determined and a side turning threshold calibration device is built
The rollover characterizing parameters are selected as a side inclination angle tau, a lateral acceleration theta and an air bag pressure transverse transfer rate eta, factors influencing the rollover threshold value of the tank car are selected as a vehicle speed v, a whole vehicle mass m and a liquid filling ratio lambda, and the eta calculation formula is as follows:
Figure FDA0002505277420000011
in the formula (1), FlsIs the pressure of the left air bag of the s-th axle of the tank car, FrsIs the air bag pressure on the right side of the s axle of the tank carForce, s is the axle number, s is 1,2, …, e, e is the total axle number of the tank car;
the rollover threshold value calibration device comprises a high-precision inertia measurement unit, a plurality of pressure sensors, two wheel force sensors and two rollover prevention frames, wherein the high-precision inertia measurement unit is arranged close to the center of mass of the tank car; the high-precision inertial measurement unit acquires a roll angle tau and a lateral acceleration theta, and the wheel force sensor acquires a vertical force G of a left wheel of a final shaftlAnd right wheel vertical force GrAcquiring the average value of information of two wheel speed sensors on a last shaft through a vehicle body CAN bus as a vehicle speed v, wherein the data output frequencies of the sensors are the same, and the vehicle mass m and the liquid filling ratio lambda are obtained by static measurement in advance;
step two: tank car rollover threshold calibration test under typical rollover scene carried out in closed test field
The division of each scene element of the tank car rollover scene library is shown as the following table:
Figure FDA0002505277420000012
after the scene elements are arranged and combined, 30 rollover scenes exist; the tank truck rollover threshold calibration test under a typical rollover scene specifically comprises the following steps:
substep 1: the calibration test is carried out on a dry and solid road surface, and the peak value adhesion coefficient of the road surface is not less than 0.9; checking the reliability of each sensor in the calibration device, the safety of the anti-rollover frame and whether potential safety hazards exist in a checking test field;
substep 2: sequentially setting driving behaviors, liquid filling ratios and longitudinal gradients in a closed test field based on the determined rollover scene, drawing test tracks of J steering and double shifting lines on the test field, and statically measuring the mass m and the liquid filling ratio lambda of the whole vehicle before the test is started;
substep 3: keeping constant speed and running according to the track, and gradually increasing by 2km/h with 32km/h as initial speed until continuingIn the process of 5 times of continuous tests, the vertical force of the tire on the side of the last shaft is 0 or the side of the rollover prevention support lands to finish the rollover threshold value calibration test under the current scene, and the roll angle tau, the lateral acceleration theta, the lateral transfer rate of the air bag pressure η, the vehicle speed v, the whole vehicle mass m, the liquid filling ratio lambda and the vertical force G of the left wheel of the last shaft in the 5 times of tests are storedlLast right wheel vertical force GrThe data of (a);
substep 4: repeating the substep 2 and the substep 3, completing threshold value calibration tests under 30 rollover scenes and storing data;
step three: processing data and dynamically calibrating rollover thresholds in different typical rollover scenarios
After 30 rollover scene calibration tests are completed, 5 times of test data are stored in each scene, and 150 times of test data are totally stored, and the roll angle tau, the lateral acceleration theta, the air bag pressure transverse transfer rate η, the vehicle speed v and the vertical force G of the left wheel of the final axle, which are acquired in each test, are filtered by adopting a weighted average valuelLast right wheel vertical force GrProcessing the data; take lateral acceleration as an example, θtFor a set of data of lateral acceleration in a certain test, t is 1,2, …, c, c is the data quantity of the lateral acceleration in the test, and the weighted mean filtering processing method is as follows:
when t is 1 and t is c,
Figure FDA0002505277420000021
when t is 2 and t is c-1,
Figure FDA0002505277420000022
when t is more than or equal to 3 and less than or equal to c-2,
Figure FDA0002505277420000023
the roll angle tau, the air bag pressure transverse transfer rate η, the vehicle speed v and the vertical force G of the left wheel of the final axle of the testlLast right wheel vertical force GrFiltering the data by the same method;
after the data of each item of 150 tests are processed by the method, the roll state H is calculated, and the formula is as follows:
Figure FDA0002505277420000031
in the formula (2), the reaction mixture is,
Figure FDA0002505277420000032
the vertical force of the left wheel of the final axle after the weighted average filtering processing,
Figure FDA0002505277420000033
the method comprises the following steps of calibrating a rollover threshold value for the vertical force of the right wheel of the last axle after weighted mean filtering treatment by using the side-tipping state H of the tank car, and specifically comprises the following steps:
substep 1: uniformly displaying the data of the roll angle, the lateral acceleration, the transverse transfer rate of the air bag pressure, the vehicle speed and the roll state of a certain test on a time axis, namely the curves of all the data have the same starting time and ending time; the whole vehicle mass and the liquid filling ratio in the test are constant values;
substep 2: searching a point at which the value H of the side-rolling state reaches 0.9 for the first time from the initial moment, recording the values of the vehicle speed, the whole vehicle mass, the liquid filling ratio, the side-rolling angle, the lateral acceleration and the air bag pressure transverse transfer rate at the moment, and recording the values as
Figure FDA0002505277420000034
Wherein
Figure FDA0002505277420000035
The rollover threshold for this trial;
substep 3: repeating the steps, recording the values of the vehicle speed, the whole vehicle mass, the liquid filling ratio, the roll angle, the lateral acceleration and the air bag pressure transverse transfer rate when the value H of the roll state in 150 tests reaches 0.9 for the first time, and recording the values as the values
Figure FDA0002505277420000036
Step four: fitting a functional relation between the rollover threshold and the influence factors thereof by using the SVR
Binding sample data
Figure FDA0002505277420000037
Respectively fitting functions f of the roll angle threshold value, the vehicle speed, the whole vehicle mass and the liquid filling ratio by using SVRτ(x) Function f of lateral acceleration threshold, vehicle speed, vehicle mass and charge ratioθ(x) Function f of lateral transfer rate threshold of air bag pressure and vehicle speed, vehicle mass and liquid filling ratioη(x);
Step five: tank car rollover threshold dynamic estimation based on fitted functional relation
When the tank car runs, the speed of the car is read through the CAN bus of the car body
Figure FDA0002505277420000038
Static measurement of vehicle mass in advance
Figure FDA00025052774200000315
And liquid filling ratio
Figure FDA00025052774200000310
To obtain
Figure FDA00025052774200000311
Respectively calculate
Figure FDA00025052774200000312
And
Figure FDA00025052774200000313
Figure FDA00025052774200000314
is a roll angle control threshold value that is,
Figure FDA0002505277420000041
is a lateral acceleration control threshold value and is,
Figure FDA0002505277420000042
controlling a threshold value for the lateral transfer rate of the air bag pressure;
Figure FDA0002505277420000043
is a roll angle pre-warning threshold value,
Figure FDA0002505277420000044
a threshold value is pre-warned for the lateral acceleration,
Figure FDA0002505277420000045
and (4) early warning threshold value for the transverse transfer rate of the air bag pressure.
CN202010444574.3A 2020-05-23 2020-05-23 Highly-reliable dynamic estimation method for rollover threshold value of tank car Active CN111695197B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010444574.3A CN111695197B (en) 2020-05-23 2020-05-23 Highly-reliable dynamic estimation method for rollover threshold value of tank car

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010444574.3A CN111695197B (en) 2020-05-23 2020-05-23 Highly-reliable dynamic estimation method for rollover threshold value of tank car

Publications (2)

Publication Number Publication Date
CN111695197A true CN111695197A (en) 2020-09-22
CN111695197B CN111695197B (en) 2023-04-11

Family

ID=72477444

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010444574.3A Active CN111695197B (en) 2020-05-23 2020-05-23 Highly-reliable dynamic estimation method for rollover threshold value of tank car

Country Status (1)

Country Link
CN (1) CN111695197B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112389450A (en) * 2020-11-20 2021-02-23 交通运输部公路科学研究所 Dynamic estimation method for overall vehicle mass of tank car based on plate spring deformation

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103213582A (en) * 2013-04-18 2013-07-24 上海理工大学 Anti-rollover warning control method based on vehicle roll angle estimation
CN108099919A (en) * 2017-11-09 2018-06-01 珠海格力电器股份有限公司 Preventing vehicle rollover method for early warning, device, storage medium and vehicle
CN109368076A (en) * 2018-08-31 2019-02-22 淮阴工学院 A kind of flow container vehicle tank anti-rollover control system and control method
CN109515368A (en) * 2018-10-29 2019-03-26 淮阴工学院 A kind of tank truck side overturning crash cushioning device and its control method
CN109878509A (en) * 2019-03-11 2019-06-14 东南大学 Monoblock type tank car Multi-source Information Fusion rollover method for early warning based on fuzzy logic

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103213582A (en) * 2013-04-18 2013-07-24 上海理工大学 Anti-rollover warning control method based on vehicle roll angle estimation
CN108099919A (en) * 2017-11-09 2018-06-01 珠海格力电器股份有限公司 Preventing vehicle rollover method for early warning, device, storage medium and vehicle
CN109368076A (en) * 2018-08-31 2019-02-22 淮阴工学院 A kind of flow container vehicle tank anti-rollover control system and control method
CN109515368A (en) * 2018-10-29 2019-03-26 淮阴工学院 A kind of tank truck side overturning crash cushioning device and its control method
CN109878509A (en) * 2019-03-11 2019-06-14 东南大学 Monoblock type tank car Multi-source Information Fusion rollover method for early warning based on fuzzy logic

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112389450A (en) * 2020-11-20 2021-02-23 交通运输部公路科学研究所 Dynamic estimation method for overall vehicle mass of tank car based on plate spring deformation

Also Published As

Publication number Publication date
CN111695197B (en) 2023-04-11

Similar Documents

Publication Publication Date Title
CN103407451B (en) A kind of road longitudinal and additional forces method of estimation
CN102112354B (en) Road surface friction coefficient estimating device and road surface friction coefficient estimating method
CN106740873B (en) Rollover early warning system and early warning method thereof
CN100562735C (en) Quantitative analysis method about the steering behaviour of steering stability of vehicles/tires
US20090177346A1 (en) Dynamic estimation of vehicle inertial parameters and tire forces from tire sensors
US6904351B1 (en) Operating a vehicle control system
CN108819950B (en) Vehicle speed estimation method and system of vehicle stability control system
CN111645698B (en) Self-adaptive estimation method for rollover threshold value of heavy-duty vehicle
CN104354700A (en) Vehicle parameter on-line estimation method based on unscented Kalman filtering
CN102165300A (en) Method and device for determining a center of gravity of a motor vehicle
CN105774458A (en) Method For Controlling Suspension System
CN1882464B (en) Vehicle motion-realizing method and device
CN111695197B (en) Highly-reliable dynamic estimation method for rollover threshold value of tank car
WO2022134929A1 (en) Method and apparatus for determining mass of vehicle, and device and medium
Tandy et al. The true definition and measurement of oversteer and understeer
CN105579824B (en) Method for simulating turning driving
CN111645670B (en) Heavy-duty vehicle roll state detection method based on support vector machine
CN111695196B (en) Tank car side-tipping state prejudging method based on vehicle body kinematics information
Ghosh et al. Sideslip angle estimation of a Formula SAE racing vehicle
CN115406669A (en) Multi-axle special vehicle rollover index optimization method
Zheng et al. A Heavy Tractor Semi-Trailer Stability Control Strategy Based on Electronic Pneumatic Braking System HIL Test
US20210221432A1 (en) Systems and methods for real-time monitoring of vehicle inertia parameter values using lateral dynamics
Ehlers et al. State and parameter estimation in a semitrailer for different loading conditions only based on trailer signals
JPH11115720A (en) Estimation device of road surface friction coefficient
Chakravarthy Development of a steer axle tire blowout model for tractor semitrailers in TruckSim

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