Method and device for analyzing the energy expenditure distribution of a motor vehicle

Abstract

A method for analyzing the distribution of energy expenditures of a motor vehicle from data from a communications network and from parameters of the vehicle includes steps in which the energy expenditures of the vehicle over a journey are calculated, the said energy expenditures are analyzed by comparing them with at least one model of the vehicle simulating the same journey, an energy balance report is formulated on the basis of the analysis of the energy expenditures and of the fuel consumption and the said energy balance report is communicated to an external server.

Claims

1. A device-implemented method for analyzing the distribution of energy expenditures of a motor vehicle from data from a communications network and from parameters of the vehicle, the method comprising: calculating energy expenditures of the vehicle over a journey; analyzing the energy expenditures by comparing them with a predetermined model of the vehicle simulating the same journey; formulating an energy balance report on the basis of the analysis of the energy expenditures and of the fuel consumption; telematically communicating the energy balance report to an external server, wherein the external server receives a plurality of the energy balance reports, corresponding to a respective plurality of journeys of the vehicle, from which a check on performance of the vehicle is made; gathering by a data processing module, prior to the calculating energy expenditures step, data from the communications network by the data processing module receiving a data transmission via the communications network, the data comprising data selected from the group consisting of the volume of fuel injected q.sub.fuel, the indicated combustion engine torque T.sub.ind, the engine speed Ω.sub.eng, the combustion engine mechanical friction torque T.sub.fric, the vehicle speed ν, the gearbox gear ratio engaged γ.sub.gear, the gradient of the roadway α, the longitudinal acceleration of the vehicle {dot over (ν)}, the position of the brake pedal P.sub.brake and the torque of the retarder T.sub.ret; gathering by the data processing module, prior to the calculating energy expenditures step, data transmitted from a database of vehicle parameters housed in a computer of the vehicle or on a remote server, the parameters being selected from the group consisting of the radius of the wheels r.sub.w, the mass of the vehicle m, the equivalent mass of the running elements m.sub.r, the rolling resistance coefficient for the tires c.sub.r, the vehicle air density ρ.sub.a, the frontal area of the vehicle A.sub.f, the coefficient of drag c.sub.x, the lower calorific value PCI, the density of the ρ.sub.fuel, the axle ratio γ.sub.axle, the gearbox efficiency η.sub.GB, the transmission efficiency μ.sub.TRANS, the inertia of the primary shaft J.sub.prim and the inertia of the running elements of the engine J.sub.eng; and gathering by the data processing module, prior to the calculating energy expenditures step, the engine speed at idle Ω.sub.idle, the mean power approximated by the auxiliary elements P.sub.aux, a map H.sup.Q(x, y) of fuel consumption, and a map H.sup.T(x) of the physical limits of the engine, wherein, during the step of calculating the energy expenditures: the energy E.sub.fuel produced by burning the fuel is calculated on the basis of the lower calorific value PCI, of the volume of fuel injected q.sub.fuel, and of the density of the fuel ρ.sub.fuel, the indicated mechanical energy E.sub.ind is calculated on the basis of the indicated combustion engine torque T.sub.ind and of the engine speed Ω.sub.eng, the energy E.sub.fric associated with internal friction within the engine is calculated on the basis of the combustion engine mechanical friction torque T.sub.fric and of the engine speed Ω.sub.eng, the energy E.sub.idle expended during the engine idling phases is calculated on the basis of the combustion engine mechanical friction torque T.sub.fric, of the engine speed Ω.sub.eng, and of the engine idling speed Ω.sub.idle, the energy E.sub.aux expended by the auxiliary elements is calculated on the basis of the mean power approximated by the auxiliary elements P.sub.aux, the mechanical energy E.sub.pwt.sup.+,wo/losses on the output side of the engine shaft is calculated on the basis of the indicated combustion engine torque T.sub.ind, of the vehicle speed ν, of the combustion engine mechanical friction torque T.sub.fric, and of the gearbox gear ratio engaged γ.sub.gear, the mechanical energies E.sub.pwt.sup.+,w/losses and E.sub.pwt.sup.−,w/losses with respect to the wheel are calculated on the basis of the mechanical energy E.sub.pwt.sup.+,wo/losses on the output side of the engine shaft, of the radius of the wheels r.sub.w, of the axle ratio γ.sub.axle, of the gearbox efficiency η.sub.GB, and of the transmission efficiency η.sub.TRANS, the potential energy E.sub.g is calculated on the basis of the mass of the vehicle m, of the vehicle speed ν, and of the gradient of the roadway α, the kinetic energy E.sub.inertia is calculated on the basis of the vehicle speed ν, of the longitudinal acceleration of the vehicle {dot over (ν)}, of the mass of the vehicle m, of the equivalent mass of the running elements m.sub.r, of the axle ratio γ.sub.axle, of the gearbox efficiency η.sub.GB, of the transmission efficiency η.sub.TRANS, of the inertia of the primary shaft J.sub.prim, of the inertia of the running elements of the engine J.sub.eng, of the gearbox gear ratio engaged and of the radius of the wheels, the energy E.sub.aero.sup.wo/wind associated with the aerodynamics of the vehicle is calculated on the basis of the vehicle air density ρ.sub.a, of the frontal area of the vehicle A.sub.f, of the coefficient of drag c.sub.x, and of the vehicle speed ν, and the energy E.sub.r associated with the rolling resistance is calculated on the basis of the mass of the vehicle m, of the rolling resistance coefficient for the tires c.sub.r and of the gradient of the roadway α.

2. The method according to claim 1, wherein, during the analyzing energy expenditures step, resistive forces are analyzed by comparing the calculated energy expenditures with a model of the vehicle, supplied by a database, simulating the same journey under standard conditions in which the gradient is zero, and the fuel consumption is analyzed by comparing a measured fuel consumption with an estimate of the fuel consumption, derived from the model of the vehicle.

3. The method according to claim 2, wherein, during the analysis of the resistive forces, the measured energy E.sub.meas.sup.res associated with the resistive forces is calculated on the basis of the mechanical energy E.sub.pwt.sup.+,wo/losses on the output side of the engine shaft, of the kinetic energy E.sub.inertia, and of the potential energy E.sub.g, the modelled energy E.sub.model.sup.res associated with the resistive forces is calculated on the basis of the energy E.sub.aero.sup.wo/wind associated with the aerodynamics and of the energy E.sub.r associated with the rolling resistance, and the energies E.sub.meas.sup.res, E.sub.model.sup.res associated with resistive forces are compared.

4. The method according to claim 2, wherein, during the analysis of the fuel consumption, a first quantity Q.sub.eng of fuel injected is calculated on the basis of the engine speed Q.sub.eng and of a map H.sup.Q(x, y) of fuel consumption.

5. The method according to claim 4, wherein, during the analysis of the fuel consumption, a second quantity Q.sub.f of fuel injected is calculated on the basis of the speed profile and of the gradient profile and of an engine map H.sup.T(x) that connects engine speed with effective torque.

6. The method according to claim 5, wherein, during the analysis of the fuel consumption, a third quantity Q.sub.f,sd of fuel injected is calculated on the basis of the speed profile and of the gradient profile under standard conditions in which the gradient is zero, and a fourth quantity Q.sub.f,ν,sd of fuel injected is calculated on the basis of the speed profile and of the gradient profile under standard conditions in which the gradient is zero and on a model of the vehicle simulating the same journey under standard conditions in which the gradient is zero and the third and fourth quantities of fuel are compared with the first and second quantities of fuel.

7. The method according to claim 1, wherein the energy expenditures are calculated only if the position of the brake pedal P.sub.brake and the torque of the retarder T.sub.ret are zero and no gearshift is in progress.

Description

BRIEF DESCRIPTION OF THE FIGURES

(1) Other objectives, features 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:

(2) FIG. 1 schematically depicts the device for analysing the distribution of the energy expenditures of a vehicle according to the invention; and

(3) FIG. 2 illustrates the steps of the method for analysing the distribution of the energy expenditures of a vehicle according to the invention.

DETAILED DESCRIPTION

(4) A device 10 for analysing the distribution of the energy expenditures of a motor vehicle comprises a data processing module 20, a module 30 for formulating an energy balance report on the motor vehicle, a module 40 for detecting the end of a journey and a module 50 for communicating the energy balance report to a server (not depicted), for example external to the vehicle, for example using telematic means.

(5) The data processing module 20 receives as input data coming from a communications network 21 better known to those skilled in the art by the abbreviation “CAN: Controller Area Network”. The data transmitted by the CAN network comprise the volume of fuel injected q.sub.fuel, the indicated combustion engine torque T.sub.ind, the engine speed Ω.sub.eng, the combustion engine mechanical friction torque T.sub.fric, the vehicle speed ν, the gearbox gear ratio engaged γ.sub.gear, the local gradient of the roadway α, the longitudinal acceleration of the vehicle {dot over (ν)}, the position of the brake pedal P.sub.brake and the torque of the retarder T.sub.ret.

(6) The data processing module 20 also receives vehicle parameters derived from a database 22 of vehicle parameters housed in a computer (not depicted) of the motor vehicle or on a remote server. The vehicle parameters transmitted to the processing module 20 comprise: the radius of the wheels r.sub.w, the mass of the vehicle m, the equivalent mass of the running gear m.sub.r. This mass is obtained from the characteristics of the corresponding components. Thus, the equivalent mass of the wheels is obtained from characteristics of the wheels and tyres, and the equivalent mass of the inertias of the drivetrain (engine, gearbox, transmission shaft . . . ) is obtained from the characteristics of the combustion engine (cylinder capacity and number of cylinders. the rolling resistance coefficient for the tyres c.sub.r. Because this coefficient is connected to the tyres, it evolves over time. In order to obtain the best possible accuracy on the mass estimate, provision could be made for this coefficient to be estimated using a special purpose estimation method. the coefficients of the aerodynamic resistance, such as the vehicle air density ρ.sub.a, the frontal area of the vehicle A.sub.f, the coefficient of drag c.sub.x. the lower calorific value PCI, the density of the fuel ρ.sub.fuel, the axle ratio γ.sub.axle, the gearbox efficiency η.sub.GB, the transmission efficiency η.sub.TRANS, the constant due to gravity g, the inertia of the primary shaft J.sub.prim and the inertia of the running elements of the engine J.sub.eng.

(7) The data processing module 20 also receives three calibration constants: the engine speed at idle Ω.sub.idle, a constant regarding the level of uncertainty around the idling speed δΩ, and a constant representing the mean power approximated by the auxiliary elements in a heavy duty vehicle P.sub.aux. These constants are, for example, calibrated according to the knowledge of a person skilled in the art. Thus, the constant P.sub.aux is fixed between 2 kW and 6 kW, the constant Ω.sub.idle is fixed at around 600 rpm, and the constant δΩ is fixed at around 200 rpm. The data processing module 20 also receives a map H.sup.Q(x, y) of fuel consumption, and a map H.sup.T(x) of the engine linking engine speed to effective torque.

(8) Finally, by way of nonlimiting example, the data processing module 20 may also receive, in real-time, the position of the vehicle in the form of a latitude, a longitude and an altitude on the basis of satellite data received by the global positioning system (not depicted) so as to determine the gradient of the roadway α.

(9) The module 30 for formulating an energy balance report for the vehicle comprises a module 32 for calculating the energies at various stages along the drivetrain, in order therefrom to deduce a distribution of the energy expenditures of the vehicle.

(10) The calculation module 32 calculates the energy E.sub.fuel produced by burning the fuel, the indicated mechanical energy E.sub.ind, the energy E.sub.fric associated with the internal friction of the engine, the energy E.sub.idle expended during engine idling phases, the energy E.sub.aux expended by the auxiliary elements, the mechanical energy E.sub.pwt.sup.+,wo/losses on the output side of the engine shaft, the mechanical energies E.sub.pwt.sup.+,w/losses and E.sub.pwt.sup.−,w/losses with respect to the wheel, the potential energy E.sub.g, the kinetic energy E.sub.inertia, the energy E.sub.aero.sup.wo/wind associated with the aerodynamics of the vehicle and the energy E.sub.r associated with the rolling resistance. The exponent “wo/wind” is an abbreviation for the expression “without wind” which means to say that the equations are implemented without the presence of a wind sensor.

(11) The energy E.sub.fuel produced by burning the fuel can be written according to the following equation:

(12) $\begin{array}{cc}{E}_{\mathrm{fuel}}=\mathrm{PCI}.Math.{\rho }_{\mathrm{fuel}}{\int }_0^{\mathrm{Tend}}{q}_{\mathrm{fuel}}\left(t\right).Math.\mathrm{dt}& \left(\mathrm{Eq}.1\right)\end{array}$

(13) Where:

(14) PCI is the lower calorific value or “lower heating value” corresponding to the amount of heat released by the complete combustion of a unit of fuel, expressed in MJ/kg;

(15) ρ.sub.fuel is the density of the fuel, expressed in kg/L;

(16) q.sub.fuel is the volume of fuel injected, expressed in L/s; and

(17) Tend is the end of the calculation of the distribution of energy expenditures.

(18) The indicated mechanical energy E.sub.ind can be written according to the following equation:

(19) $\begin{array}{cc}{E}_{\mathrm{ind}}={\int }_0^{\mathrm{Tend}}{T}_{\mathrm{ind}}\left(t\right).Math.{\Omega }_{\mathrm{eng}}\left(t\right).Math.\mathrm{dt}& \left(\mathrm{Eq}.2\right)\end{array}$

(20) Where:

(21) T.sub.ind is the indicated combustion engine torque, expressed in N.Math.m; and

(22) Ω.sub.eng is the engine speed, expressed in rad.Math.s.sup.−1.

(23) The energy E.sub.fric associated with the internal friction of the engine can be written according to the following equation:

(24) $\begin{array}{cc}{E}_{\mathrm{fric}}={\int }_0^{\mathrm{Tend}}{T}_{\mathrm{fric}}\left(t\right).Math.{\Omega }_{\mathrm{eng}}\left(t\right).Math.\mathrm{dt}& \left(\mathrm{Eq}.3\right)\end{array}$

(25) Where:

(26) T.sub.fric is the combustion engine mechanical friction torque, expressed in N.Math.m.

(27) In order not to account twice for the energy associated with the internal friction of the engine, it is desirable to make a distinction between the energy E.sub.frix.sup.w/idle associated with the internal friction of the engine during engine idling phases, and the energy E.sub.frix.sup.wo/idle associated with the internal friction of the engine during engine idling phases outside of the engine idling phases:

(28) $\begin{array}{cc}{E}_{\mathrm{frix}}^{w/\mathrm{idle}}={\int }_0^{\mathrm{Tend}}{T}_{\mathrm{fric}}\left(t\right).Math.{\Omega }_{\mathrm{eng}}\left(t\right).Math.\mathrm{dt}.Math.s.t..Math.{\Omega }_{\mathrm{eng}}-{\Omega }_{\mathrm{idle}}.Math.<\mathrm{\delta \Omega }& \left(\mathrm{Eq}.4\right)\end{array}$

(29) When |Ω.sub.eng−Ω.sub.idle|<δΩ

(30) $\begin{array}{cc}{E}_{\mathrm{frix}}^{\mathrm{wo}/\mathrm{idle}}={\int }_0^{\mathrm{Tend}}{T}_{\mathrm{fric}}\left(t\right).Math.{\Omega }_{\mathrm{eng}}\left(t\right).Math.\mathrm{dt}& \left(\mathrm{Eq}.5\right)\end{array}$

(31) When |Ω.sub.eng−Ω.sub.idle|≥δΩ

(32) Where:

(33) Ω.sub.idle is the engine idling speed expressed in rad.Math.s.sup.−1; and

(34) δΩ is a constant relating to the level of uncertainty around the idling speed, expressed in rad.Math.s.sup.−1.

(35) The energy E.sub.idle expended during the engine idling phases can be written according to the equation Eq. 4:

(36) $E_{\mathrm{idle}}=E_{\mathrm{frix}}^{w/\mathrm{idle}}={\int }_0^{\mathrm{Tend}}{T}_{\mathrm{fric}}\left(t\right).Math.{\Omega }_{\mathrm{eng}}\left(t\right).Math.\mathrm{dt}$

(37) The energy E.sub.aux expended by the auxiliary elements can be written according to the following equation:

(38) $\begin{array}{cc}{E}_{\mathrm{aux}}=P_{\mathrm{aux}}{\int }_0^{\mathrm{Tend}}t.Math.\mathrm{dt}& \left(\mathrm{Eq}.6\right)\end{array}$

(39) Where:

(40) P.sub.aux is a constant representing the mean power approximated by the auxiliary elements in a heavy duty vehicle, given by a database, expressed in kW.

(41) The mean power P.sub.aux can be broken down into the auxiliary power needed for vehicle traction, notably needed by the air compressor, the power steering, the alternator, into auxiliary power associated with comfort, notably the air conditioning, the heating etc., and into the auxiliary power specific to the power take-off systems, notably when the vehicle is stationary.

(42) The mechanical energy at the wheel E.sub.pwt.sup.+,wo/losses corresponds to the energy of the powerplant, and can be written according to the following equation:

(43) $\begin{array}{cc}{E}_{\mathrm{pwt}}^{+,\mathrm{wo}/\mathrm{losses}}=\frac{{\gamma }_{\mathrm{axle}}}{r_w}{\int }_0^{\mathrm{Tend}}\left(T_{\mathrm{ind}}\left(t\right)-T_{\mathrm{fric}}\left(t\right)\right).Math.v\left(t\right).{\gamma }_{\mathrm{gear}}\left(t\right).Math.\mathrm{dt}& \left(\mathrm{Eq}.7\right)\end{array}$

(44) When T.sub.ind(t)>0

(45) Where:

(46) γ.sub.axle is the axle ratio, which is dimensionless;

(47) r.sub.w is the radius of the wheel, expressed in m;

(48) ν is the speed of the vehicle, expressed in m.Math.s.sup.−1;

(49) γ.sub.gear is the gearbox gear ratio engaged, which is dimensionless.

(50) The mechanical energies E.sub.pwt.sup.+,w/losses and E.sub.pwt.sup.−,w/losses with respect to the wheel, can be written according to the following equations:

(51) $\begin{array}{cc}{E}_{\mathrm{pwt}}^{+,w/\mathrm{losses}}={\eta }_{\mathrm{GB}}.Math.{\eta }_{\mathrm{Trans}}.Math.E_{\mathrm{pwt}}^{+,\mathrm{wo}/\mathrm{losses}}& \left(\mathrm{Eq}.8\right)\\ E_{\mathrm{pwt}}^{-,w/\mathrm{losses}}=\frac{1}{{\eta }_{\mathrm{GB}}{\eta }_{\mathrm{Trans}}}.Math.E_{\mathrm{pwt}}^{+,\mathrm{wo}/\mathrm{losses}}& \left(\mathrm{Eq}.9\right)\end{array}$

(52) Where:

(53) η.sub.GB is the gearbox efficiency, which is dimensionless; and

(54) η.sub.TRANS is the transmission efficiency, which is dimensionless.

(55) The potential energy E.sub.g can be written according to the following equation:

(56) 0$\begin{array}{cc}{E}_g=m.g.{\int }_0^{\mathrm{Tend}}v\left(t\right).Math.\sin \left(\alpha \left(t\right)\right).Math.\mathrm{dt}& \left(\mathrm{Eq}.10\right)\end{array}$

(57) Where:

(58) m is the mass of the vehicle, expressed in kg;

(59) g is the constant due to gravity, expressed in N/kg; and

(60) α is the gradient of the roadway, expressed in rad.

(61) It may also be beneficial to determine the positive potential energy stored up during the journey:

(62) $\begin{array}{cc}{E}_g^{+}={\int }_0^{\mathrm{Tend}}v\left(t\right).Math.F_g\left(t\right).Math.\mathrm{dt}& \left(\mathrm{Eq}.11\right)\end{array}$

(63) When sin(α(t))>0

(64) The force induced by gravity F.sub.g(t) is determined in the remainder of the description using equation Eq.20.

(65) The kinetic energy E.sub.inertia can be written according to the following equation:

(66) $\begin{array}{cc}{E}_{\mathrm{intertia}}=E_{\mathrm{intertia}}^v+E_{\mathrm{inertia}}^{\mathrm{eng}}& \left(\mathrm{Eq}.12\right)\\ \mathrm{Where}\text{:}& \\ E_{\mathrm{inertia}}^v=\left(m+m_r\right).Math.{\int }_0^{\mathrm{Tend}}\stackrel{.}{v}\left(t\right).Math.v\left(t\right).Math.\mathrm{dt}\mathrm{when}v\left(t\right)>v_0& \left(\mathrm{Eq}.13\right)\\ E_{\mathrm{inertia}}^{\mathrm{eng}}=\frac{{\eta }_{\mathrm{eng}}.Math.{\eta }_{\mathrm{Trans}}.Math.\left(J_{\mathrm{prim}}+J_{\mathrm{eng}}\right).Math.{\gamma }_{\mathrm{axle}}^2.Math.}{r_w^2}{\int }_0^{\mathrm{Tend}}\stackrel{.}{v}\left(t\right).Math.v\left(t\right).Math.{\gamma }_{\mathrm{gear}}^2\mathrm{dt}\mathrm{when}v\left(t\right)>v_0& \left(\mathrm{Eq}.14\right)\end{array}$

(67) Where:

(68) m.sub.r is the equivalent mass of the running elements, expressed in kg;

(69) {dot over (ν)} is the longitudinal acceleration of the vehicle, expressed in m.Math.s.sup.−2;

(70) J.sub.prim is the inertia of the primary shaft, expressed in kg.Math.m.sup.2;

(71) J.sub.eng is the inertia of the running elements of the engine, expressed in kg.Math.m.sup.2; and

(72) ν.sub.0 is a threshold constant so as not to consider very low speeds, ν.sub.0 is, for example, equal to 1 km/h.

(73) The energy E.sub.aero.sup.w/wind associated with the aerodynamics of the vehicle can be written according to the following equation:

(74) $\begin{array}{cc}{E}_{\mathrm{aero}}^{w/\mathrm{wind}}=\frac{1}{2}.Math.{\rho }_a.Math.A_f.Math.c_x.Math.{\int }_0^{\mathrm{Tend}}{\left(v\left(t\right)-v_{\mathrm{wind}}\left(t\right)\right)}^2.Math.v\left(t\right).Math.\mathrm{dt}& \left(\mathrm{Eq}.15\right)\end{array}$

(75) Where:

(76) ρ.sub.a is the density of the air, expressed in kg.Math.m.sup.−3;

(77) A.sub.f is the frontal area of the vehicle, expressed in m.sup.2;

(78) c.sub.x is the coefficient of drag, which is dimensionless; and

(79) ν.sub.wind is the wind speed, expressed in m.Math.s−1.

(80) The wind speed is information that is not always available at vehicle level. In one exemplary embodiment, this wind speed is obtained by external means not detailed here and which do not form part of the present invention.

(81) However, the wind speed ν.sub.wind can be neglected, notably in cases in which the information is not available, and in that case, equation Eq. 15 can be written:

(82) $E_{\mathrm{aero}}^{w/\mathrm{wind}}=\frac{1}{2}.Math.{\rho }_a.Math.A_f.Math.c_x.Math.{\int }_0^{\mathrm{Tend}}{v}^3\left(t\right).Math.\mathrm{dt}$

(83) The energy E.sub.r associated with the rolling resistance can be written according to the following equation:

(84) $\begin{array}{cc}{E}_r=c_r.Math.m.Math.{\int }_0^{\mathrm{Tend}}v\left(t\right).Math.\cos \left(\alpha \left(t\right)\right).Math.\mathrm{dt}& \left(\mathrm{Eq}.16\right)\end{array}$

(85) Where:

(86) c.sub.r is the rolling resistance coefficient for the tyres, which is dimensionless.

(87) The module 30 for formulating an energy balance report on the vehicle further comprises a module 34 for analysing the energy expenditures of the vehicle as determined by the energy calculation module 32.

(88) The module 34 for analysing the energy expenditures comprises a module 36 for analysing the resistive forces that compares the energy flows calculated in Equations 1 to 16 with a physical model of the vehicle, provided by a database 35, simulating the same journey. Reference may be made to Equations 25 to 32 below for the dynamic models of vehicle fuel consumption.

(89) The vehicle dynamics equation which connects the forces acting on the vehicle with the mass of the vehicle yields Equation Eq. 17:
(m.sub.ν+m.sub.r).Math.{dot over (ν)}(t)=F.sub.t(t)−F.sub.ext(t)−F.sub.dist(t)  (Eq. 17)

(90) Where:

(91) m.sub.v is the mass of the vehicle, expressed in kg;

(92) m.sub.r is the equivalent mass of the running elements, expressed in kg;

(93) {dot over (ν)} is the longitudinal acceleration of the vehicle, expressed in m.Math.s.sup.−2;

(94) F.sub.t(t) is the traction force applied to the vehicle, expressed in N;

(95) F.sub.ext(t) is the external forces applied to the vehicle, expressed in N; and

(96) F.sub.dist(t) is the force associated with the modelling errors.

(97) The set of external forces F.sub.ext(t) applied to the vehicle combines the aerodynamic friction force F.sub.a(t), the rolling resistance force F.sub.r(t) and the force induced by gravity F.sub.g(t).

(98) The aerodynamic friction force F.sub.a(t) is dependent on the vehicle air density ρ.sub.a, on the frontal area of the vehicle A.sub.f, on the coefficient of drag c.sub.x and on the relative speed at which the vehicle makes progress with respect to the air (ν−ν.sub.wind) The aerodynamic friction force F.sub.a(t) can be written according to the following equation:
F.sub.a(ν)=½.Math.ρ.sub.a.Math.A.sub.f.Math.c.sub.x.Math.(ν−ν.sub.wind).sup.2  (Eq. 18)

(99) The relative wind speed ν.sub.wind is unknown and so to calculate Equation Eq.3, it is assumed that this speed is negligible, insofar as it has no impact on the accuracy of the energy expenditures analysis.

(100) The rolling resistance force F.sub.r(t) is dependent on the mass of the vehicle m.sub.ν, on the rolling resistance coefficient for the tyres c.sub.r and on the gradient of the roadway α as described in the following equation:
F.sub.r(ν)=c.sub.r.Math.m.sub.ν.Math.cos(α)  (Eq. 19)

(101) The force induced by gravity F.sub.g(t) is dependent on the mass of the vehicle on the gravity constant g and on the gradient of the roadway α as described in the following equation:
F.sub.g(α)=m.sub.ν.Math.g.Math.sin(α)  (Eq. 20)

(102) From an energy standpoint, Equation Eq. 17 can be written:
E.sub.pwt−E.sub.inertia−E.sub.g=E.sub.r+E.sub.aero.sup.w/wind  (Eq.21)

(103) Because the braking torque and the torque of the retarder are not available, the energies are calculated only when the engine torque is positive, namely under the following condition:
T.sup.+(t)=[ν(t)>ν0;T.sub.ret(t)>1,P.sub.brake(t)=0;γ.sub.gear(t)>0.1;T.sub.ret(t)=0]

(104) Where:

(105) P.sub.brake is the position of the brake pedal, which is dimensionless; and

(106) T.sub.ret is the torque of the retarder, expressed as a %.

(107) Thus, this is used to deduce the following equation:

(108) $\begin{array}{cc}{E}_{\mathrm{pwt}}^{{+}^{T^{+}}}-E_{\mathrm{inertia}}^{T^{+}}-E_g^{T^{+}}=E_r^{T^{+}}+E_{\mathrm{aero}}^{w\text{/}{\mathrm{wind}}^{T^{+}}}& \left(\mathrm{Eq}.22\right)\end{array}$

(109) The measured energy E.sub.meas.sup.res associated with resistive forces can be written according to the following equation:

(110) $\begin{array}{cc}{E}_{\mathrm{meas}}^{\mathrm{res}}=E_{\mathrm{pwt}}^{{+}^{T^{+}}}-E_{\mathrm{inertia}}^{T^{+}}-E_g^{T^{+}}=E_r^{T^{+}}& \left(\mathrm{Eq}.23\right)\end{array}$

(111) The modelled energy E.sub.model.sup.res associated with resistive forces can be written according to the following equation:

(112) $\begin{array}{cc}{E}_{\mathrm{model}}^{\mathrm{res}}=E_r^{T^{+}}+E_{\mathrm{aero}}^{w\text{/}{\mathrm{wind}}^{T^{+}}}& \left(\mathrm{Eq}.24\right)\end{array}$

(113) The module 36 for analysing the resistive forces compares the measured energy with that from the model: E.sub.meas.sup.res−E.sub.model.sup.res.

(114) If the two energies are equal: E.sub.meas.sup.res−E.sub.model.sup.res, that means that the model is correctly reproducing the measurements. This also indicates that the influence of the wind ν.sub.wind is small.

(115) If E.sub.meas.sup.res>E.sub.model.sup.res, then the resistive forces in actual fact have been more significant than the model had calculated. That may be explained by the fact that one of the constants c.sub.x or c.sub.r is higher than the value indicated in the model or that the wind on the journey has been of appreciable strength without being compensated for over the entirety of the journey.

(116) If E.sub.meas.sup.res<E.sub.model.sup.res, then the resistive forces in actual fact have been less significant than the model had calculated. That may be explained by the fact that one of the constants c.sub.x or c.sub.r is lower than the value indicated in the model.

(117) There are various sources of error that may be present when calculating the energy breakdown, notably delays on measurements coming from various sensors or inherent errors associated with the modelling of the vehicle. It is thus preferable to analyse the distribution of energy expenditures over a plurality of journeys.

(118) It may also be emphasized that once the energy balance report has been created and communicated to an external server, it is possible to analyse not the absolute value of the difference between the modelled energy and the measured energy, but how this difference evolves from one journey to another, the analysis being performed using the one same model over the plurality of journeys.

(119) In addition to the analysis of the resistive forces that enable vehicle parameter validation, an analysis of the fuel consumption is also made.

(120) The module 34 for analysing the energy expenditures for this purpose comprises a module 38 for analysing the fuel consumption by comparing the measured consumption with consumptions taken from various models derived from databases simulating the same journey under standard conditions. Reference may be made to Equations 25 to 32 below for the dynamic models of vehicle fuel consumption.

(121) Thus, it is possible to normalize the journey and the vehicle so as to be able subsequently to compare the performance of several vehicles and rank them according to their fuel consumption against a standard usage.

(122) The module 38 for analysing the fuel consumption calculates the quantity .sub.eng of fuel injected from the engine speed and from a consumption map (x,y):

(123) $\begin{array}{cc}{Q}_{\mathrm{eng}}=\frac{1}{\rho }{\int }_0^{\mathrm{Tend}}{\Omega }_{\mathrm{end}}\left(t\right).Math.H^Q({\Omega }_{\mathrm{eng}}\left(t\right),\left(T_{\mathrm{ind}}\left(t\right)-T_{\mathrm{fric}}\left(t\right)\right).Math.\mathrm{dt}\mathrm{When}{\Omega }_{\mathrm{eng}}\left(t\right)>200\frac{\pi }{180}& \left(\mathrm{Eq}.25\right)\end{array}$

(124) Where:

(125) (x,y) is a map of the engine linking engine speed and effective torque with a quantity of fuel injected.

(126) The module 38 for analysing the fuel consumption also calculates the quantity .sub.f of fuel injected from the speed profile and from the gradient profile:

(127) 0$\begin{array}{cc}{Q}_f=\frac{1}{\rho }{\int }_0^{\mathrm{Tend}}{\Omega }_{\mathrm{eng}}\left(t\right).Math.H^Q\left({\Omega }_{\mathrm{eng}}\left(t\right),h\left(t\right)\right).Math.\mathrm{dt}& \left(\mathrm{Eq}.26\right)\\ \mathrm{When}{\Omega }_{\mathrm{eng}}\left(t\right)>200\frac{\pi }{180}& \\ \mathrm{Where}\text{:}& \\ h\left(t\right)={\mathrm{SAT}}_{\underset{\_}{\mathrm{Te}}}^{\stackrel{\_}{\mathrm{Te}}}.Math.\left(\frac{r_w}{{\gamma }_{\mathrm{gear}}{\lambda }_{\mathrm{axle}}}\left(F_{\mathrm{inertia}}+F_r+F_{\mathrm{aero}}+F_g\right)\right)& \left(\mathrm{Eq}.27\right)\\ \stackrel{\_}{\mathrm{Te}}=\max \left(H^T.Math.\left({\Omega }_{\mathrm{eng}}\left(t\right)\right)\right)& \left(\mathrm{Eq}.28\right)\\ \underset{\_}{\mathrm{Te}}=\min \left(H^T.Math.\left({\Omega }_{\mathrm{eng}}\left(t\right)\right)\right)& \left(\mathrm{Eq}.29\right)\end{array}$

(128) H.sup.T(x) is a map of the engine linking engine speed with effective torque.

(129) The module for analysing the fuel consumption then calculates the quantity .sub.f,sd of fuel injected from the speed profile and from the gradient profile under standard conditions, namely conditions in which α=0:

(130) $\begin{array}{cc}{Q}_{f,\mathrm{sd}}=\frac{1}{\rho }{\int }_0^{\mathrm{Tend}}{\Omega }_{\mathrm{eng}}\left(t\right).Math.H^Q\left({\Omega }_{\mathrm{eng}}\left(t\right),h_{\mathrm{sd}}\left(t\right)\right).Math.\mathrm{dt}& \left(\mathrm{Eq}.30\right)\\ \mathrm{When}{\Omega }_{\mathrm{eng}}\left(t\right)>200\frac{\pi }{180}& \\ \mathrm{Where}\text{:}& \\ h_{\mathrm{sd}}\left(t\right)={\mathrm{Sat}}_{\underset{\_}{\mathrm{Te}}}^{\stackrel{\_}{\mathrm{Te}}}.Math.\left(\frac{r_w}{{\gamma }_{\mathrm{gear}}{\lambda }_{\mathrm{axle}}}\left(F_{\mathrm{inertia}}+F_r+F_{\mathrm{aero}}\right)\right)& \left(\mathrm{Eq}.31\right)\end{array}$

(131) The module for analysing the fuel consumption also calculates the quantity .sub.f,ν,sd of fuel injected from the speed profile and from the gradient profile under standard conditions, and on a normalized vehicle, namely conditions in which α=0:

(132) $\begin{array}{cc}{Q}_{f,v,\mathrm{sd}}=\frac{1}{\rho }{\int }_0^{\mathrm{Tend}}{\Omega }_{\mathrm{eng}}\left(t\right).Math.H^Q\left({\Omega }_{\mathrm{eng}}\left(t\right),h_v{,}_{\mathrm{sd}}\left(t\right)\right).Math.\mathrm{dt}& \left(\mathrm{Eq}.32\right)\\ \mathrm{When}{\Omega }_{\mathrm{eng}}\left(t\right)>200\frac{\pi }{180}& \\ \mathrm{Where}\text{:}& \\ h_{v,\mathrm{sd}}\left(t\right)={\mathrm{Sat}}_{\underset{\_}{\mathrm{Te}}}^{\stackrel{\_}{\mathrm{Te}}}.Math.\left(\frac{r_w}{{\gamma }_{\mathrm{gear}}{\lambda }_{\mathrm{axle}}}\left(\left(F_{\mathrm{inertia}}+F_r\right).Math.\frac{m_0}{m}+F_{\mathrm{aero}}\right)\right)& \left(\mathrm{Eq}.33\right)\end{array}$

(133) m is the mass of the vehicle; and

(134) m.sub.0 is the mass of the normalized vehicle. What is considered here is a vehicle of constant mass, such as an 18 T vehicle, so as to compare, over a similar journey, two vehicles of different masses with respect to a similar reference.

(135) The fuel consumptions calculated in Equations Eq. 25 to Eq. 32 are then compared with each other in order to determine whether the consumption of the vehicle is normal.

(136) The gearshift occurrence needs to be zero during all the calculation operations.

(137) Once the energy balance report has been formulated, and the end of the journey has been detected, it is transmitted to an external server by telematic means.

(138) It is thus possible to group, over several journeys, and for the one same vehicle, the energy balance report formulated for each journey and to make a reliable check on the performance of the vehicle.

(139) It is also possible to group the energy balance reports over several journeys, and for a plurality of motor vehicles, in order to detect, with in the one same fleet of vehicles, those which have a high fuel consumption, and those journeys which require a high fuel consumption, for example because of changes in level, or because of a strong wind presence.

(140) The flow diagram depicted in FIG. 2 illustrates the method 60 for analysing the energy distribution of a motor vehicle.

(141) During a first step 61, data is collected, the data comprising: the volume of fuel injected q.sub.fuel, the indicated combustion engine torque T.sub.ind, the engine speed Ω.sub.eng, the combustion engine mechanical friction torque T.sub.ind, the vehicle speed ν, the gearbox gear ratio engaged γ.sub.gear, the gradient of the roadway α, the longitudinal acceleration of the vehicle {dot over (ν)}, the position of the brake pedal P.sub.brake and the torque of the retarder T.sub.ret, the radius of the wheels r.sub.w, the mass of the vehicle m, the equivalent mass of the running elements m.sub.r, the rolling resistance coefficient for the tyres c.sub.r, the aerodynamic resistance coefficients, such as the vehicle air density ρ.sub.a, the frontal area of the vehicle A.sub.f, the coefficient of drag c.sub.x, the lower calorific value PCI, the density of the fuel ρ.sub.fuel, the axle ratio γ.sub.axle, the gearbox efficiency η.sub.GB, the transmission efficiency η.sub.TRANS, the inertia of the primary shaft J.sub.prim and the inertia of the running elements of the engine J.sub.eng.

(142) Three calibration constants are also gathered: the engine speed at idle Ω.sub.idle, a constant regarding the level of uncertainty around the idling speed δΩ, and a constant representing the mean power approximated by the auxiliary elements in a heavy-duty vehicle P.sub.aux, as well as a map H.sup.Q (x,y) of fuel consumption and a map H.sup.T(x) of the engine linking engine speed with effective torque.

(143) By way of nonlimiting example, it is also possible, in real-time, to collect the position of the vehicle in the form of a latitude, a longitude and an altitude on the basis of satellite data received by the global positioning system GPS (not depicted) so as to determine the gradient of the roadway α.

(144) During a second step 62, the following energies are calculated: the energy E.sub.fuel produced by burning the fuel, the indicated mechanical energy E.sub.ind, the energy E.sub.fric associated with the internal friction of the engine, the energy E.sub.idle expended during engine idling phases; the energy E.sub.aux expended by the auxiliary elements, the mechanical energy E.sub.pwt.sup.+,wo/losses on the output side of the engine shaft, the mechanical energies E.sub.pwt.sup.+,w/losses and E.sub.pwt.sup.−,w/losses with respect to the wheel, the potential energy E.sub.g, the kinetic energy E.sub.inertia, the energy E.sub.aero.sup.wo/wind associated with the aerodynamics of the vehicle and the energy E.sub.r associated with the rolling resistance.

(145) During a third step 63, the energy expenditures of the vehicle, as determined during the second step 62, are analysed by analysing, in step 64, the resistive forces by comparing the energy flows calculated in Equations 1 to 16 with a physical model of the vehicle, and by analysing, in step 66, the fuel consumption by comparing the measured fuel consumption with the fuel consumptions taken from models 67 derived from databases simulating the same journey under standard conditions, namely without gradient. Thus, it is possible to normalize the journey and the vehicle so as to be able subsequently to compare the performance of several vehicles and rank them according to their fuel consumption against a standard usage.

(146) During a fourth step 68, an energy balance report is formulated for the vehicle, indicating the fuel consumption, and making it possible to detect a potential anomaly between the measurements and the model (the presence of wind, a significant difference in the coefficient of drag . . . )

(147) By virtue of the invention, it is possible to normalize the journey and the vehicle so as to be able subsequently to compare the performance of several vehicles and rank them according to their fuel consumption against a standard usage. It is also possible to ensure the validity of the parameters of the vehicle model and to rank all of the vehicles of the one same fleet according to their fuel consumption.