METHOD AND DEVICE FOR DETERMINING AN ESTIMATE OF THE TOTAL MASS OF A MOTOR VEHICLE
20190226905 ยท 2019-07-25
Inventors
Cpc classification
International classification
Abstract
A method determines the total mass of an automotive vehicle on the basis of data of a communication network and parameters of the vehicle, in which an estimation of the total laden mass (mv,est) of the vehicle, of the speed of the vehicle (vest) and of the slope of the road (est) is determined at an instant (k) by applying the fundamental equation of dynamics and as a function of the values of the total mass of the vehicle, of the speed of the vehicle and of the slope of the road at a previous instant (k1).
Claims
1.-14. (canceled)
15. A method for determining the total mass of an automotive vehicle on a basis of data of a communication network and of parameters of the vehicle, in which data are retrieved from a database of parameters of the vehicle, the parameters hosted in a computer of the vehicle, comprising the steps of: retrieving setting parameters including execution time of the method for estimating the mass T, covariance matrix R of measurement noise, covariance matrix Q of the modelling indeterminacies with regard to speed and slope, and covariance matrix Q.sub.M of the modelling indeterminacies with regard to the mass; retrieving data originating from the communication network including a speed of the vehicle v, a downshift ratio of a transmission , a level of depression of a brake pedal, a level of use of a retarder, and an estimation of a motor torque provided by the heat engine T.sub.e and of the following mechanical frictions: a motive force F.sub.t(t) generated by the vehicle, determined as a function of the motor torque value T.sub.e, of the downshift ratio of the transmission , of the efficiency coefficient of the transmission .sub.gb, and of a radius of a wheel r.sub.w, an aerodynamic friction force F.sub.a(t), determined as a function of the data retrieved from the database of parameters, a rolling resistance force F.sub.r(t), determined as a function of the mass of the vehicle m.sub.v, of a coefficient of rolling resistance of the tires c.sub.r, and of a slope of a road , and a force induced by gravity F.sub.g(t), determined as a function of the mass of the vehicle m.sub.v, of a gravity constant g and of the slope of the road ; and determining an estimation of a total laden mass m.sub.v,est of the vehicle, of the speed of the vehicle v.sub.est and of the slope of the road .sub.est at an instant k by applying the fundamental equation of dynamics and as a function of the values of the total mass of the vehicle, of the speed of the vehicle and of the slope of the road at a previous instant k1.
16. The method of determination according to claim 15, wherein, in order to determine an estimation of the total laden mass m.sub.v,est of the vehicle, of the speed of the vehicle v.sub.est and of the slope of the road .sub.est, two filters are applied in parallel.
17. The method of determination according to claim 15, further comprising the step of, at the same time as the determining the estimation step, determining an estimation of two covariance matrices P, of dimension 22, expressing a covariance of an error in a speed variable and slope variable and P.sub.M, of dimension 11, expressing a covariance of an error in a mass variable as a function of the covariance matrix Q of the modelling indeterminacies with regard to the speed and the slope and of the covariance matrix Q.sub.M of the modelling indeterminacies with regard to the mass.
18. The method of determination according to claim 17, further comprising the step of, prior to the retrieving the data originating from the communication network step, reinitializing at the initial instant the mass of the vehicle m.sub.v, the speed of the vehicle v, the slope of the road , and the two covariance matrices P, P.sub.M.
19. The method of determination according to claim 15, further comprising the step of, prior to the determining the estimation step, verifying at the instant k that all conditions are fulfilled for undertaking a computation of the estimation of the mass.
20. The method of determination according to claim 19, wherein the conditions are that the speed of the vehicle v is positive and that the position of the brake pedal and the level of use of the retarder are zero.
21. The method of determination according to claim 15, wherein a value of the estimation error (k) is computed by comparing the value of the estimation of the speed v.sub.est determined with the value of the speed v.sub.CAN measured by the communication network.
22. The method of determination according to claim 15, wherein at least one proportionality coefficient L, L.sub.M is computed as a function of the covariance matrices P.sub.est, P.sub.M,est and of the equations of dynamics.
23. The method of determination according to claim 21, wherein the estimations of the speed v.sub.est of the vehicle, of the slope of the road .sub.est and of the mass of the vehicle m.sub.v,est which are determined are corrected as a function of the value of the estimation error (k) and of the proportionality coefficient L.
24. The method of determination according to claim 15, further comprising the step of, after the determining the estimation step, verifying that a journey of the vehicle has terminated, and when the journey of the vehicle has terminated, computing an average value of the mass Avg(m.sub.v) of the vehicle at the end of the journey as a function of the estimation of the mass of the vehicle m.sub.v,est and transmitting the computed average value of the mass Avg(m.sub.v) to a server.
25. A device for determining a total mass of an automotive vehicle on a basis of data of a communication network and of parameters of the vehicle, comprising: a module for processing the data which is able to receive, as input, data of parameters of the vehicle, which parameters are hosted in a computer of the vehicle and comprise an equivalent mass of rolling elements m.sub.r, a coefficient of rolling resistance of the tires c.sub.r, a density of air of the vehicle .sub.a, a frontal surface area of the vehicle A.sub.f, a drag coefficient c.sub.x, and a radius of a wheel r.sub.w of the vehicle, wherein the module for processing the data receives, as input, setting parameters including execution time of the method for estimating the mass T, covariance matrix R of measurement noise, covariance matrix Q of the modelling indeterminacies with regard to speed and slope, and covariance matrix Q.sub.M of the modelling indeterminacies with regard to the mass, and data originating from a communication network including a speed of the vehicle v, a downshift ratio of a transmission , a level of depression of a brake pedal, a level of use of a retarder, and an estimation of a motor torque provided by the heat engine T.sub.e and of the following mechanical frictions: a motive force F.sub.t(t) generated by the vehicle, determined as a function of the motor torque value T.sub.e, of the downshift ratio of the transmission , of the efficiency coefficient of the transmission .sub.gb, and of a radius of a wheel r.sub.w, an aerodynamic friction force F.sub.a(t), determined as a function of the data retrieved from the database of parameters, a rolling resistance force F.sub.r(t), determined as a function of the mass of the vehicle m.sub.v, of a coefficient of rolling resistance of the tires c.sub.r, and of a slope of a road , and a force induced by gravity F.sub.g(t), determined as a function of the mass of the vehicle m.sub.v, of a gravity constant g and of the slope of the road ; and a computation module for computing an estimation of the total laden mass m.sub.v,est of the vehicle, of the speed of the vehicle v.sub.est and of the slope of the road .sub.est at an instant k by applying the fundamental equation of dynamics and as a function of the values of the total mass of the vehicle, of the speed of the vehicle and of the slope of the road at a previous instant k1.
26. The determination device according to claim 25, wherein the computation module comprises two filters disposed in parallel and a module for estimating two covariance matrices P, of dimension 22, expressing a covariance of an error in a speed variable and slope variable and P.sub.M, of dimension 11, expressing a covariance of an error in a mass variable as a function of the covariance matrix Q of the modelling indeterminacies with regard to the speed and the slope and of the covariance matrix Q.sub.M of the modelling indeterminacies with regard to the mass.
27. The determination device according to claim 26, wherein the computation module comprises a correction module for correcting the estimation values determined comprising a module for computing an estimation error (k) comparing the value of the estimation of the speed v.sub.est determined with the value of the speed v.sub.CAN measured by the network and a module for computing a proportionality coefficient L, L.sub.M as a function of the covariance matrices P.sub.est, P.sub.M,est and of the equations of dynamics.
28. The determination device according to claim 25, further comprising an end-of-journey detection module, a module for computing an average value of the mass Avg(m.sub.v) of the vehicle at the end of the journey, and a communication module able to transmit the average value of the mass Avg(m.sub.v) of the vehicle to a server.
Description
[0040] Other aims, characteristics and advantages of the invention will become apparent on reading the following description, given solely by way of nonlimiting example and with reference to the appended drawings in which:
[0041]
[0042]
[0043] A device 10 for determining an estimation of the mass of an automotive vehicle comprises a module for processing the data 20, a computation module 30 for computing the estimation of the mass of the vehicle and of the slope of the road, an end-of-journey detection module 40, a final computation module 50 for finally computing an average value of the mass of the vehicle at the end of the journey, and a communication module 60 for communicating the average value of the mass of the vehicle.
[0044] The module for processing the data 20 receives, as input, data retrieved from a database 21 of parameters of the vehicle, which parameters are hosted in a computer (not represented) of the automotive vehicle. The vehicle parameters transmitted to the processing module comprise: [0045] the equivalent mass of the rolling elements m.sub.r, expressed in kg. This mass is obtained on the basis of the characteristics of the corresponding components. Thus, the equivalent mass of the wheels is obtained on the basis of the characteristics of the wheels and tyres and the equivalent mass of the inertias of the drivetrain (engine, gearbox, transmission shaft . . . ) is obtained on the basis of the characteristics of the heat engine (capacity and number of cylinders). [0046] the coefficient of rolling resistance of the tyres c.sub.r. Since this coefficient is related to the tyres, it evolves over time. In order to obtain the best possible accuracy in the estimation of the mass, provision could be made to estimate this coefficient by a specific estimation method. [0047] the coefficients of the aerodynamic resistance, such as the density of the air of the vehicle .sub.a, the frontal surface area of the vehicle A.sub.f, the drag coefficient c.sub.x. [0048] the radius of the wheels r.sub.w.
[0049] The module for processing the data 20 also receives setting parameters from a settings database 22, such as the execution time of the algorithm for estimating the mass T, the covariance matrix of the measurement noise R, the covariance matrix of the modelling indeterminacies with regard to the speed and the slope and the covariance matrix
.sub.M of the modelling indeterminacies with regard to the mass.
[0050] The data processing module 20 furthermore receives data originating from a communication network 23 well known to the person skilled in the art under the acronym CAN: Controller Area Network. The data transmitted by the CAN network comprise the speed of the vehicle v, the downshift ratio of the gearbox, the level of depression of the brake pedal, the level of use of the retarder and an estimation of the torque provided by the heat engine and of these mechanical frictions.
[0051] The downshift ratio of the gearbox, the motor torque, the mechanical friction torque of the heat engine and the radius of the wheels are used to compute the known so-called control variable u.
[0052] Finally, the module for processing the data 20 can also receive in real time the position of the vehicle in the form of latitude, longitude and altitude as a function of the satellite data received by the global positioning system (not represented).
[0053] The module for processing the data 20 comprises an initialization module for initializing, at the instant t=0, the mass of the vehicle my, the speed of the vehicle v, the slope of the road , as well as two covariance matrices P, of dimension 22, expressing the covariance of the error in the speed variable and slope variable and P.sub.M, of dimension 11, expressing the covariance of the error in the mass variable.
[0054] The module for processing the data 20 comprises a verification module 25 for verifying, at the instant t=k, that all the conditions are fulfilled for undertaking the computation of the estimation of the mass. Stated otherwise, to undertake the computation of the estimation of the mass, the speed of the vehicle v must be strictly positive and the position of the brake pedal and the level of use of the retarder must be zero. As a variant, it would also be possible to envisage conditioning the computation of the estimation of the mass upon the engaging of a ratio of the gearbox.
[0055] If one of these conditions is not satisfied, the module for computing the estimation of the mass will undertake the computation of the mass and of the slope as a function of the data at the instant t=k1.
[0056] The computation module 30 for computing the estimation of the mass of the vehicle and of the slope of the road uses an algorithm with a recursive approach based on a mathematical model of the longitudinal dynamics of the vehicle.
[0057] The equation of dynamics for the vehicle which links the loadings acting on the vehicle to the mass of the vehicle makes it possible to obtain equation Eq. 1:
[0058] With: [0059] m.sub.v, the mass of the vehicle, expressed in kg; [0060] m.sub.r, the equivalent mass of the rolling elements, expressed in kg; [0061] v, the speed of the vehicle, expressed in m.Math.s.sup.2; [0062] F.sub.t(t), the traction force applied to the vehicle, expressed in N; and [0063] F.sub.ext(t), the exterior forces applied to the vehicle, expressed in N.
[0064] The traction force F.sub.t(t) may be written according to the following equation:
F.sub.t(v)=T.sub.e(t).Math.(t).Math..sub.gb.Math.r.sub.w(Eq. 2)
[0065] With: [0066] T.sub.e, the torque of the heat engine, expressed in N.Math.m; [0067] , the downshift coefficient of the transmission; [0068] .sub.gb, the efficiency coefficient of the transmission; and [0069] r.sub.w, the radius of the wheel, expressed in m.
[0070] The set of exterior forces F.sub.ext(t) applied to the vehicle groups together the aerodynamic friction force Fa(t), the rolling resistance force F.sub.r(t) and the force induced by gravity F.sub.g(t).
[0071] The aerodynamic friction force Fa(t) depends on the density of the air of the vehicle .sub.a, on the frontal surface area of the vehicle A.sub.f, on the drag coefficient c.sub.x and on the relative forward speed of the vehicle with respect to the air (vv.sub.wind). The aerodynamic friction force Fa(t) may be written according to the following equation:
F.sub.a(v)=.Math..sub.a.Math.A.sub.f.Math.c.sub.x.Math.(vv.sub.wind).sup.2(Eq. 3)
[0072] The speed relatively of the wind v.sub.wind is not known, in order to compute equation Eq. 3, the assumption is made that this speed is negligible, in so far as it does not have any impact on the accuracy of the estimation of the mass of the vehicle.
[0073] The rolling resistance force F.sub.r(t) depends on the mass of the vehicle m.sub.v, on the coefficient of rolling resistance of the tyres C.sub.r and on the slope of the road , as described in the following equation:
F.sub.r(v)=c.sub.r.Math.m.sub.v.Math.cos()(Eq. 4)
The force induced by gravity F.sub.g(t) depends on the mass of the vehicle m.sub.v, on the gravity constant g and on the slope of the road , as described in the following equation:
F.sub.g()=m.sub.v.Math.g.Math.sin()(Eq. 5)
[0074] The computation module 30 comprises a first module 31 for estimating the speed v.sub.est of the vehicle, the slope of the road .sub.est and the mass of the vehicle m.sub.v,est. The first estimation module 31 comprises a first, so-called Kalman, filter F.sub.1 allowing the filtering of the measurement of the speed and of the estimation of the slope of the road travelled and a second, so-called Kalman, filter F.sub.2 making it possible to estimate the mass of the vehicle. The two filters are executed in parallel so that the information is pooled between the two filters and the variables estimated by the first filter are necessary for the estimation of the second filter, and vice versa.
[0075] The estimation of the speed v.sub.est of the vehicle is computed according to the following equation:
[0076] The estimation of the slope of the road .sub.est is computed according to the following equation:
sin(.sub.est(k))=sin((k1))(Eq. 7)
[0077] The estimation of the mass of the vehicle m.sub.v,est is computed according to the following equation:
m.sub.v,est(k)=m.sub.v(k1)(Eq. 8)
[0078] The first estimation module 31 determines equations 6 to 8 as a function of equations 1 to 5, of the measured variables and of the values of the mass of the vehicle, of the speed of the vehicle and of the slope of the road which were computed at the instant t=k1.
[0079] The computation module 30 also comprises a second estimation module 32 for estimating the two covariance matrices P.sub.est and P.sub.M,is on the basis of the values at the instant t=k1, of equations 1 to 5 and of a so-called Riccati equation whose sought-after variable is the covariance matrix of the estimation errors:
P.sub.est(k)=A(k).Math.P(k1).Math.A(k).sup.T+(Eq. 9)
P.sub.M,est(k)=P.sub.M(k1)+.sub.M(Eq. 10)
[0080] With: [0081] A, the tangent linearized dynamic matrix (linearization of equations Eq. 6 Eq. 7 with respect to the speed and the slope) [0082] , the covariance matrix of the modelling indeterminacies with regard to the speed and the slope and the matrix; and [0083]
, the covariance matrix of the modelling indeterminacies with regard to the mass.
[0084] The computation module 30 furthermore comprises a correction module 33 for correcting the estimation values computed in equations 6 to 8 comprising a module 34 for computing the estimation error (k) comparing the value of the estimation of the speed v.sub.est computed by the first estimation module with the value of the speed v.sub.CAN measured by the CAN network:
(k)=v.sub.est(k)v.sub.CAN(k)(Eq. 11)
[0085] The correction module 33 comprises a module 35 for computing the proportionality coefficients L and L.sub.M, termed gain of the observer, computed on the basis of the covariance matrices P.sub.est and P.sub.M,is and of the equations of dynamics 1 to 5:
L(k)=(L.sub.1(k),L.sub.2(k)).sup.T=P.sub.est(k).Math.C.sup.T(k).Math.(C(k).Math.P.sub.est(k).Math.C.sup.T(k)+R).sup.1(Eq. 12)
L.sub.M(k)=P.sub.M,est(k).Math.C.sub.M.sup.t(k).Math.(C.sub.M(k).Math.P.sub.M,est(k).Math.C.sub.M.sup.T(k)+R).sup.1(Eq. 13)
[0086] With: [0087] L.sub.1, L.sub.2, defined by equation Eq. 12; and [0088] C, C.sub.M, the tangent linearized output matrices (linearization of the output equations of the system with respect to the state.fwdarw.C=(0,1) C.sub.M=Cd.sub.f/d.sub.M, where f=Eq. 6 and Eq. 7)
[0089] With:
P(k)=P.sub.est(k)L(k).Math.C(k).Math.P.sub.est(k)(Eq. 14)
P.sub.M(k)=P.sub.M,est(k)L.sub.M(k).Math.C.sub.M(k).Math.P.sub.M,est(k)(Eq. 15)
[0090] The corrected value of the speed v of the vehicle is computed according to the following equation:
v(k)=v.sub.est(k)+L.sub.1(k).Math.(k)(Eq. 16)
[0091] The corrected value of the slope of the road is computed according to the following equation:
sin((k))=sin(.sub.est(k))+L.sub.2(k).Math.(k)(Eq. 17)
[0092] The corrected value of the mass of the vehicle my is computed according to the following equation:
m.sub.v(k)=m.sub.v,est(k)+L.sub.M(k).Math.(k)(Eq. 18)
[0093] The end-of-journey detection module 40 verifies whether the journey has terminated, for example, with the aid of a means of timeout and of comparing the timeout with a threshold value. If the journey has not terminated, the computation module 30 repeats equations 6 to 18.
[0094] When the journey has terminated, the final computation module 50 computes an average value of the mass Avg(m.sub.v) of the vehicle at the end of the journey. The communication module 60 thereafter transmits the average value of the mass Avg(m.sub.v) of the vehicle to a server (not represented), for example outside the vehicle, for example by telematic means.
[0095] The flowchart represented in
[0096] During a first step 110, the following data are retrieved: the equivalent mass of the rolling elements m.sub.r; the coefficient of rolling resistance of the tyres c.sub.r, the coefficients of the aerodynamic resistance, such as the density of the air of the vehicle .sub.a, of the frontal surface area of the vehicle A.sub.f, of the drag coefficient c.sub.x; the radius of the wheels r.sub.w; and the setting parameters such as the execution time of the algorithm for estimating the mass T, the covariance matrix of the measurement noise R, the covariance matrix of the modelling indeterminacies with regard to the speed and the slope and the covariance matrix
.sub.M of the modelling indeterminacies with regard to the mass.
[0097] By way of nonlimiting example, it is possible to also retrieve, in real time, the position of the vehicle in the form of latitude, longitude and altitude as a function of the satellite data received by the GPS global positioning system (not represented).
[0098] During a second step 120, at the instant t=0, the mass of the vehicle m.sub.v, the speed of the vehicle v, the slope of the road (), as well as two covariance matrices P, of dimension 22, expressing the covariance of the error in the speed variable and slope variable and P.sub.M, of dimension 11, expressing the covariance of the error in the mass variable are initialized.
[0099] During a third step 130, data originating from the CAN communication network are retrieved, such as the speed of the vehicle v, the downshift ratio of the gearbox, the level of depression of the brake pedal, the level of use of the retarder and an estimation of the torque provided by the heat engine and of these mechanical frictions.
[0100] During a fourth step 140, it is verified, at the instant t=k, that all the conditions are fulfilled for undertaking the computation of the estimation of the mass. Stated otherwise, to undertake the computation of the estimation of the mass, the speed of the vehicle v must be strictly positive and the position of the brake pedal and the level of use of the retarder must be zero. As a variant, it would also be possible to envisage conditioning the computation of the estimation of the mass upon the engaging of a ratio of the gearbox.
[0101] If one of these conditions is not satisfied, the module for computing the estimation of the mass will undertake the computation of the mass and of the slope as a function of the data at the instant t=k1. Stated otherwise, the values are pegged at the previous instant k1.
[0102] During a fifth step 151, the estimation of the speed v.sub.est of the vehicle, of the slope of the road .sub.est and of the mass of the vehicle m.sub.v is computed by using an algorithm with a recursive approach based on a mathematical model of the longitudinal dynamics of the vehicle resulting from equations 1 to 5 hereinabove.
[0103] At the same time, during a step 152, the estimation of the two matrices P.sub.est and P.sub.M,is is computed according to equations 9 and 10:
P.sub.est(k)=A(k).Math.P(k1).Math.A(k).sup.T+(Eq. 9)
P.sub.M,est(k)=P.sub.M(k1)+.sub.M(Eq. 10)
[0104] During a sixth step 160, the value of the estimation error (k) is computed by comparing the value of the estimation of the speed v.sub.est computed by the first estimation module with the value of the speed v.sub.CAN measured by the CAN network:
(k)=v.sub.est(k)v.sub.CAN(k)(Eq. 11)
[0105] During a seventh step 170, a proportionality coefficient L, termed gain of the observer, is computed according to the above equations 12 to 15. The proportionality coefficient L computed on the basis of the covariance matrices P.sub.est and P.sub.M,is and of the equations of dynamics 1 to 5 and the value of the vehicle speed v.sub.est computed according to equation 6, the value of the slope of the road .sub.est computed according to equation 7 and the value of the vehicle mass of the vehicle m.sub.v,est computed according to equation 8 are corrected.
[0106] During an eighth step 180, it is verified whether the journey has terminated, for example, with the aid of a means of timeout and of comparing the timeout with a threshold value. If the journey has not terminated, steps 130 to 170 are repeated.
[0107] When the journey has terminated, an average value of the mass Avg(m.sub.v) of the vehicle at the end of the journey is computed, in step 190, and the average value of the mass Avg(m.sub.v) computed is transmitted, in step 200, to a server (not represented), for example outside the vehicle, for example by telematic means.
[0108] Thus, the total mass of a vehicle can be estimated in a reliable and fast manner, doing so even under condition of slope of the road travelled.
[0109] The vehicle can be an automotive vehicle or a train coach or a trolley or a lifting apparatus.
[0110] By mass, by total mass of an automotive vehicle, is meant the mass of the vehicle, including its cargo in the case of a semi-trailer or of articulated vehicles and/or, including its passengers in the case of a cargoless vehicle or of a bus.