Procédé et dispositif d'estimation de la masse d'un véhicule automobile La présente invention concerne un procédé et un dispositif d'estimation de la masse totale d'un véhicule automobile. La connaissance de la masse totale d'un véhicule automobile est nécessaire au bon fonctionnement de nombreux dispositifs embarqués dans le véhicule, notamment les dispositifs de gestion de freinage ou de gestion de boîte de vitesse automatique.The present invention relates to a method and a device for estimating the total mass of a motor vehicle. Knowledge of the total mass of a motor vehicle is necessary for the proper functioning of many devices embedded in the vehicle, including brake management devices or automatic transmission management.
Par ailleurs, pour les véhicules utilitaires offrant une capacité de chargement importante, une estimation de la masse du véhicule permet à l'utilisateur de garantir le respect du poids total en charge autorisé. Il est donc souhaitable d'obtenir rapidement une estimation fiable de la masse du véhicule, même lorsque le véhicule est engagé sur une pente. Aujourd'hui l'estimation de la masse d'un véhicule automobile peut se faire de différentes façons, par exemple par mesure directe de la masse, par estimation à l'aide de capteurs d'assiettes et à partir de la raideur des amortissements. Toutefois, ces méthodes sont chères et ne permettent pas d'obtenir une estimation de la masse précise dû à la dispersion de fabrication et au vieillissement des raideurs des suspensions. Il existe également des procédés d'estimation de la masse d'un véhicule automobile par application du principe fondamental de la dynamique à un seul instant. On peut se référer à cet égard au document US 6, 249, 735 qui décrit un procédé d'estimation de l'état de chargement d'un véhicule en fonction du couple du moteur mesuré et de l'accélération du véhicule calculée par approximation discrète de la dérivée de la vitesse et par filtrage. Un tel procédé ne permet pas d'obtenir une estimation de la masse suffisamment fiable et robuste. On peut également se référer au document FR 2 857 090 qui décrit un procédé d'estimation de la masse reposant sur un algorithme de moindres carrés récursifs en utilisant l'accélération du véhicule. Toutefois, un tel procédé d'estimation ne permet pas d'estimer rapidement la masse du véhicule automobile. De plus, les procédés d'estimation de la masse d'un véhicule automobile par application du principe fondamental de la dynamique à un seul instant sont particulièrement imprécis et lourds à mettre en oeuvre car ils nécessitent une estimation de nombreux paramètres du véhicule et se basent sur des hypothèses qui ne sont pas toujours vérifiées, telles qu'une vitesse du vent nulle, une déclivité nulle, etc.Moreover, for commercial vehicles offering a large load capacity, an estimate of the weight of the vehicle allows the user to guarantee compliance with the total permissible laden weight. It is therefore desirable to quickly obtain a reliable estimate of the mass of the vehicle, even when the vehicle is engaged on a slope. Today the estimation of the mass of a motor vehicle can be done in different ways, for example by direct measurement of the mass, by estimation with the help of plate sensors and from the stiffness of depreciation. However, these methods are expensive and do not allow to obtain an estimate of the precise mass due to the manufacturing dispersion and the aging of the stiffness of the suspensions. There are also methods for estimating the mass of a motor vehicle by applying the fundamental principle of momentum dynamics. Reference may be made in this regard to document US Pat. No. 6,249,735 which describes a method for estimating the loading state of a vehicle as a function of the torque of the engine measured and the acceleration of the vehicle calculated by discrete approximation. of the derivative of the speed and by filtering. Such a method does not make it possible to obtain a sufficiently reliable and robust estimate of the mass. Reference can also be made to document FR 2 857 090 which describes a mass estimation method based on a recursive least squares algorithm using vehicle acceleration. However, such an estimation method does not make it possible to rapidly estimate the mass of the motor vehicle. In addition, the methods for estimating the mass of a motor vehicle by applying the fundamental principle of the dynamics at a single instant are particularly imprecise and cumbersome to implement because they require an estimation of many parameters of the vehicle and are based on on assumptions that are not always verified, such as zero wind speed, zero gradient, etc.
Par ailleurs, dans certains véhicules actuels, l'utilisateur doit appuyer sur un bouton de commande lorsque le véhicule est chargé. Le bouton de commande permet de configurer le véhicule dans un mode « chargé » ou « non chargé » pour activer un programme correspondant de pilotage du calculateur de la transmission automatique. Toutefois, l'implantation d'un tel bouton de commande est relativement coûteuse. La présente invention a donc pour objet de remédier à ces inconvénients. Le but de l'invention est donc de fournir un procédé et un dispositif permettant d'obtenir rapidement une information sur la masse totale du véhicule automobile, sans nécessiter un appui sur un bouton de commande par l'utilisateur du véhicule. L'invention a pour objet le pilotage du calculateur d'un dispositif embarqué, tel que par exemple une transmission automatique, un dispositif d'estimation de la pression des pneus, etc..., en fonction de la masse du véhicule automobile. L'invention concerne tout type d'engin muni de roues. Dans un mode de réalisation, l'invention concerne un procédé d'estimation de la masse totale d'un véhicule automobile dans lequel on calcule la différence entre les accélérations sans freinage d'une roue du véhicule à partir de l'équation fondamentale de la dynamique appliquée à deux instants temporels et l'on estime la masse du véhicule par un algorithme de moindres carrés récursifs à facteur d'oubli à partir de cette différence d'accélération.Moreover, in some current vehicles, the user must press a command button when the vehicle is loaded. The control button is used to configure the vehicle in a "loaded" or "unloaded" mode to activate a corresponding program for controlling the automatic transmission computer. However, the implementation of such a button is relatively expensive. The present invention therefore aims to overcome these disadvantages. The object of the invention is therefore to provide a method and a device for quickly obtaining information on the total mass of the motor vehicle without requiring a button on the control button by the user of the vehicle. The object of the invention is to control the computer of an on-board device, such as, for example, an automatic transmission, a device for estimating tire pressure, etc., depending on the mass of the motor vehicle. The invention relates to any type of vehicle equipped with wheels. In one embodiment, the invention relates to a method for estimating the total mass of a motor vehicle in which the difference between the braking accelerations of a vehicle wheel is calculated from the fundamental equation of the vehicle. dynamic applied at two time instants and the mass of the vehicle is estimated by a recursive least squares algorithm with forgetting factor from this acceleration difference.
Ainsi, on obtient une estimation fiable et rapide de la masse du véhicule par l'utilisation d'un algorithme simple d'utilisation et ne nécessitant pas un stockage important de données grâce au facteur d'oubli.Thus, a reliable and rapid estimation of the mass of the vehicle is obtained by using an easy-to-use algorithm and not requiring a large data storage thanks to the forgetting factor.
Avantageusement, on récupère les données comprenant l'inertie de la roue, l'accélération de la roue à deux instants, la résultante du couple moteur ramenée à la roue, le rayon de la roue et l'accélération longitudinale du véhicule à deux instants. Dans un mode de réalisation, on filtre la valeur de la masse estimée par l'algorithme de moindres carrés récursifs à facteur d'oubli. La résultante du couple moteur ramenée à la roue peut être déterminée par une cartographie. Le facteur d'oubli est, par exemple, compris entre 0,99 et 0,9959. Selon un deuxième aspect, l'invention concerne un dispositif d'estimation de la masse totale d'un véhicule automobile, comportant une unité de commande électronique destinée à calculer la masse totale du véhicule automobile et à piloter un dispositif embarqué dans le véhicule. Le dispositif d'estimation comprend un module de calcul de la différence entre les accélérations sans freinage d'une roue du véhicule à partir de l'équation fondamentale de la dynamique appliquée à deux instants temporels et un module d'estimation de la masse du véhicule par un algorithme de moindres carrés récursifs à facteur d'oubli à partir de cette différence d'accélération. Avantageusement, le dispositif d'estimation comprend un module de récupération des données comprenant l'inertie de la roue, l'accélération de la roue à deux instants, la résultante du couple moteur ramenée à la roue, le rayon de la roue et l'accélération longitudinale du véhicule à deux instants. Le dispositif d'estimation peut comprendre un filtre passe bas destiné à filtrer la valeur de la masse totale estimée par le module d'estimation.Advantageously, the data are collected including the inertia of the wheel, the acceleration of the wheel at two moments, the resultant of the engine torque brought back to the wheel, the radius of the wheel and the longitudinal acceleration of the vehicle at two times. In one embodiment, the value of the estimated mass is filtered by the forgetting factor recursive least squares algorithm. The resultant motor torque brought back to the wheel can be determined by mapping. The forgetting factor is, for example, between 0.99 and 0.9959. According to a second aspect, the invention relates to a device for estimating the total mass of a motor vehicle, comprising an electronic control unit for calculating the total mass of the motor vehicle and driving a device embedded in the vehicle. The estimating device comprises a module for calculating the difference between the accelerations without braking of a wheel of the vehicle from the fundamental equation of the dynamics applied at two time instants and a module for estimating the mass of the vehicle. by a recursive least squares algorithm with forgetting factor from this acceleration difference. Advantageously, the estimating device comprises a data recovery module comprising the inertia of the wheel, the acceleration of the wheel at two moments, the resultant of the motor torque brought back to the wheel, the radius of the wheel and the longitudinal acceleration of the vehicle at two moments. The estimation device may comprise a low-pass filter for filtering the value of the total mass estimated by the estimation module.
D'autres buts, caractéristiques et avantages de l'invention apparaîtront à la lecture de la description suivante, donnée uniquement à titre d'exemple non limitatif, et faite en référence aux dessins annexés, sur lesquels : - la figure 1 représente schématiquement le dispositif d'estimation de la masse totale d'un véhicule selon l'invention; et - la figure 2 illustre les étapes du procédé d'estimation de la masse totale d'un véhicule selon l'invention. Le dispositif d'estimation de la masse totale M d'un véhicule automobile (non représenté), selon la figure 1, comporte une unité de commande électronique 1 destinée à piloter un calculateur 2, par exemple de la transmission automatique (non représentée), de l'estimation de la pression des pneumatiques ou de tout autre dispositif du véhicule automobile nécessitant la connaissance de la masse totale du véhicule. L'unité de commande 1 comprend un module de traitement des données 3, un module de calcul 4 de la différence entre les accélérations d'une roue à deux instants temporels et un module d'estimation de la masse 5 par un algorithme des moindres carrés récursifs à facteur d'oubli. Le module de traitement des données 3 reçoit en entrée des données comprenant : Jr, l'inertie de la roue, exprimée en kg.m2 ; 6.,.(t), l'accélération de la roue à deux instants ti et t2, exprimée en rad. s_2 ; Cr(t), la résultante du couple moteur ramenée à la roue, exprimé en N.m; R , le rayon de la roue, exprimé en m ; et a, l'accélération de la roue à deux instants ti et t2, exprimée en m.s-2. L'équation du bilan des couples au niveau d'une des roues non motrices déterminée par l'application du principe fondamental de la dynamique permet d' obtenir l' équation Eq. 1 correspondant l'application du principe fondamental de la dynamique en rotation d'une roue avant équivalente lors d'une phase d'accélération sans freinage : JYOY (t) = Cr(t) - Fx(t).R - Cf(t) (Eq. 1) Avec : Jr, l'inertie de la roue, exprimée en kg.m2 ; 6.,.(t), l'accélération de la roue correspondant à la dérivée de la vitesse de rotation de la roue, exprimée en rad.s-2 ; Cr(t), la résultante du couple moteur ramenée à la roue, exprimé en N.m ; Fx(t), la force de la route sur la roue, exprimée en N ; R , le rayon de la roue, exprimé en m ; et Cf(t), le couple lié au frottement, exprimé en N.m. L'équation du bilan des forces au niveau de la roue avant équivalente lors d'une phase d'accélération sans freinage est déterminée par le principe fondamental de la dynamique en translation, selon l'équation suivante : M .a(t) = Fx(t)- Fext(t) (Eq. 2) Avec : M, la masse du véhicule, exprimé en kg ; a, l'accélération longitudinale du véhicule, exprimée en m.s-2 ; et Fext(t), l'effort extérieur appliqué sur le véhicule, exprimé en N En remplaçant, dans l'équation 2, l'expression de la force de la route sur la roue Fx(t) calculée à l'équation 1, on obtient, à deux instants t1 et t2, les équations 3 et 4 suivantes: Jrd (t2) = Cr(t2) - (M .a(t2)- Fext(t2)).R - Cf(t2) (Eq. 3) Jré.,. (h) = Cr(ti) - (M Fext(ti)).R - Cf(ti) (Eq. 4) Le module de traitement des données 3 calcule les équations 3 et 4 ci-dessus par application des principes fondamentaux de la dynamique à deux instants ti et t2.Other objects, 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: FIG. 1 schematically represents the device estimating the total mass of a vehicle according to the invention; and FIG. 2 illustrates the steps of the method for estimating the total mass of a vehicle according to the invention. The device for estimating the total mass M of a motor vehicle (not shown), according to FIG. 1, comprises an electronic control unit 1 intended to control a computer 2, for example of the automatic transmission (not represented), the estimation of the pressure of the tires or any other device of the motor vehicle requiring knowledge of the total mass of the vehicle. The control unit 1 comprises a data processing module 3, a calculation module 4 of the difference between the accelerations of a wheel with two time instants and a mass estimation module 5 by a least squares algorithm. recursive forgetting factors. The data processing module 3 receives input data including: Jr, the inertia of the wheel, expressed in kg.m2; 6.,. (T), the acceleration of the two-moment wheel ti and t2, expressed in rad. s_2; Cr (t), the resultant of the motor torque brought back to the wheel, expressed in Nm; R, the radius of the wheel, expressed in m; and a, the acceleration of the two-moment wheel t1 and t2, expressed in m.s-2. The equation of the torque balance at one of the non - driving wheels determined by the application of the fundamental principle of dynamics allows to obtain the equation Eq. 1 corresponding to the application of the fundamental principle of the dynamics in rotation of an equivalent front wheel during an acceleration phase without braking: JYOY (t) = Cr (t) - Fx (t) .R - Cf (t) ) (Eq. 1) With: Jr, the inertia of the wheel, expressed in kg.m2; 6.,. (T), the acceleration of the wheel corresponding to the derivative of the speed of rotation of the wheel, expressed in rad.s-2; Cr (t), the resultant of the motor torque brought back to the wheel, expressed in Nm; Fx (t), the force of the road on the wheel, expressed in N; R, the radius of the wheel, expressed in m; and Cf (t), the frictional torque, expressed in N.m. The equilibrium equilibrium equation for the equivalent front wheel during an acceleration phase without braking is determined by the fundamental principle of the dynamics in translation, according to the following equation: M .a (t) = Fx (t) - Fext (t) (Eq. 2) With: M, the mass of the vehicle, expressed in kg; a, the longitudinal acceleration of the vehicle, expressed in m.s-2; and Fext (t), the external force applied to the vehicle, expressed in N By replacing, in equation 2, the expression of the force of the road on the wheel Fx (t) calculated in equation 1, at two instants t1 and t2, we get the following equations 3 and 4: Jrd (t2) = Cr (t2) - (M .a (t2) - Fext (t2)) R - Cf (t2) (Eq. 3) Jré.,. (h) = Cr (ti) - (M Fext (ti)) R - Cf (ti) (Eq.4) The data processing module 3 calculates equations 3 and 4 above by applying the fundamental principles of the two-moment dynamic ti and t2.
Le module de calcul 4 calcule la différence entre les équations 3 et 4. En prenant pour hypothèse que les deux instants t2 et ti sont suffisamment proches pour que les efforts extérieurs Fext(t), Cf(t) appliqués sur le véhicule ne varient pas significativement, c'est-à-dire que Fe),(t2)-F,,,(0-, 0 et Cf(t2)-Cf(t1),--- 0, on obtient l'équation 5: JrAé, = ACr - RM.Aa (Eq. 5) Où: Ad., - 6., (t2) -6.,(t1) ACr -Cr(t2)-Cr(t1), provient, par exemple, d'une cartographie Aa= a(t2) - a(t1) On peut réécrire l'équation 5 sous la forme : (OC, - JrA0) (Eq. 6) AaM = (1) Y En appliquant, au module d'estimation 5, un algorithme basé sur la méthode des moindres carrés récursifs à facteur d'oubli k, on recalcule une nouvelle estimation des paramètres à chaque instant t sur la base de l'estimation précédente. On pose : R(k) = f(k) = Exk-to (t) T (t) (Eq.7) t=1 k (Eq.8) Xk-t 0(0 .Y(t) t=1 Avec : (1)T la transposée de cl), À, un facteur de pondération qui prend toujours une valeur positive inférieure à 1, par exemple compris entre 0,99 et 0,9959. À, est également appelé « facteur d'oubli » et peut être issu de cartographie en fonction du nombre de points acquis, qui peut être mise à jour en cas d'arrêt et/ou ouverture de porte (descente ou montée de passager). En remplaçant ces valeurs dans l'équation 6, on en déduit l'équation suivante : Û(k)- R(k)-1 f (k) (Eq. 9) Avec : , la masse estimée par la méthode des moindres carrés Les variables R(k) et f(k) peuvent être mis sous forme récursive sans aucune approximation selon les équations suivantes : R(k) = X/r?(k -1)+0(k).0T(k) (Eq.10) f(k) = X(k) f (k -1) +(1) (k).y(k) (Eq.11) En remplacement les équations 10 et 11 dans l'équation 9, on en déduit : M (k) = R(k)-1k(k). f (k -1) + ) + (I) (k) . y (k)] = R(k)-1k(k).R(k - 1) .û (k -1) + (I) (k) . y (k) M(k) = R(k)-11.[R(k) - (k).OT (k)lû (k -1) + (k) .y (k).1 SoitÛ(k) = û (k -1) + R(k)-1 (k)i.y(k) - T (k) (k -1) .1 (Eq.12) Appliquer la méthode des moindres carrés récursifs à facteur d'oubli se traduit par calculer les trois équations suivantes à chaque nouvel échantillon k : Û(k)=Û(k -1)+ L(k)i.y(k)- 0(k) û (k - 1)] (Eq.13) P (k -1)0(k) L(k) - (Eq. 14) (X +0 (k)P (k -1» (k)) [1- L(k» (k)]P (k -1) P (k) (Eq.15) X Avec P (k) = R(k)-1 En effet, la solution consistant à calculer la masse M par la méthode des moindres carrés récursifs classique a l'avantage d'être optimale et rapide. Toutefois, tous les échantillons k ont la même importance. En d'autres termes, si après une heure de roulage, la masse du véhicule varie, par exemple lorsqu'une personne descend du véhicule, il faudra attendre au moins une heure de roulage pour que les nouveaux points compensent les anciens. L'attribution d'un poids aux échantillons k par la méthode connue dans la littérature comme des moindres carrés récursifs à facteur d'oubli permet d'attribuer un poids à un échantillon d'autant plus important que l'échantillon est proche dans le passé. Le module d'estimation 5 transmet ensuite la valeur de la masse estimée à un filtre passe bas 6 avant d'être transmise au calculateur 2.The calculation module 4 calculates the difference between the equations 3 and 4. Assuming that the two instants t2 and ti are sufficiently close so that the external forces Fext (t), Cf (t) applied to the vehicle do not vary. significantly, that is to say that Fe), (t2) -F ,,, (0-, 0 and Cf (t2) -Cf (t1), --- 0, we obtain equation 5: JrAé , = ACr - RM.Aa (Eq.5) Where: Ad., - 6., (t2) -6., (T1) ACr -Cr (t2) -Cr (t1), comes, for example, from a map Aa = a (t2) - a (t1) We can rewrite equation 5 in the form: (OC, - JrA0) (Eq.6) AaM = (1) Y By applying, to the estimation module 5 , an algorithm based on the recursive least squares method with forgetting factor k, we recalculate a new estimate of the parameters at each instant t on the basis of the previous estimate: R (k) = f (k) = Exk-to (t) T (t) (Eq.7) t = 1k (Eq.8) Xk-t 0 (0.Y (t) t = 1 With: (1) T the transpose of cl) , At, a weighting factor that always takes a p-value ositive less than 1, for example between 0.99 and 0.9959. To, is also called "forgetting factor" and can be derived from mapping based on the number of points acquired, which can be updated in the event of stopping and / or opening of the door (descent or passenger climb). By replacing these values in equation 6, we deduce the following equation: Û (k) - R (k) -1 f (k) (Eq.9) With: the mass estimated by the least squares method The variables R (k) and f (k) can be put into recursive form without any approximation according to the following equations: R (k) = X / r? (K -1) + 0 (k) .0T (k) ( Eq.10) f (k) = X (k) f (k -1) + (1) (k) .y (k) (Eq.11) Replacing equations 10 and 11 in equation 9, deduce from this: M (k) = R (k) -1k (k). f (k -1) +) + (I) (k). y (k)] = R (k) -1k (k) .R (k - 1) .t (k -1) + (I) (k). y (k) M (k) = R (k) -11. [R (k) - (k) .OT (k) û (k -1) + (k) .y (k) .1 Let ((k ) = û (k -1) + R (k) -1 (k) iy (k) - T (k) (k -1) .1 (Eq.12) Apply the least-squares recursive least squares method oblivion is calculated by calculating the following three equations for each new sample k: Û (k) = Û (k -1) + L (k) y (k) - 0 (k) û (k - 1)] (Eq. 13) P (k -1) 0 (k) L (k) - (Eq 14) (X + O (k) P (k -1) (k)) [1- L (k "(k)] P (k -1) P (k) (Eq.15) X With P (k) = R (k) -1 Indeed, the solution consisting of calculating the mass M by the recursive least squares method advantage of being optimal and fast, however, all the samples k have the same importance, ie if after one hour of driving, the mass of the vehicle varies, for example when a person gets off the vehicle, it will be necessary to wait at least one hour for the new points to compensate for the old ones The allocation of a weight to the samples k by the method known in the literature as recursive least squares with forgetting factor makes it possible to assign a weight to a sample all the more important since the sample is close in the past. The estimation module 5 then transmits the value of the estimated mass to a low pass filter 6 before being transmitted to the computer 2.
L'organigramme représenté sur la figure 2 illustre le procédé d'estimation de la masse totale d'un véhicule automobile. Lors d'une première étape 10, on récupère les données comprenant l'inertie de la roue Jr, l'accélération de la roue à deux instants 6.,.(ti),6,.(t2), la résultante du couple moteur ramenée à la roue Cr(ti) Cr(t2), le rayon de la roue R et l'accélération longitudinale du véhicule à deux instant a(ti), a(t2). Lors d'une deuxième étape 20, on calcule les équations 1 et 2 par application du principe fondamental de la dynamique en rotation et en translation à deux instants temporels ti, t2 et on calcule l'accélération sans freinage d'une roue du véhicule 6.,.(0,6,.(t2) respectivement à chacun des deux instants ti, t2, selon les équations 3 et 4. Lors d'une troisième étape 30, on calcule l'équation 5 correspondant à la différence Aé.,. entre les accélérations sans freinage 6,.(0,6,.(t2)calculées précédemment, et en prenant pour hypothèse que les instants t2 et t1 sont suffisamment proches pour que les efforts extérieurs Fext(t), Cf(t) appliqués sur le véhicule ne varient pas significativement.The flowchart shown in FIG. 2 illustrates the method for estimating the total mass of a motor vehicle. In a first step 10, the data comprising the inertia of the wheel Jr, the acceleration of the two-moment wheel 6, (ti), 6, (t2), the resultant of the driving torque are recovered. brought back to the wheel Cr (ti) Cr (t2), the radius of the wheel R and the longitudinal acceleration of the vehicle at two times a (ti), a (t2). In a second step 20, equations 1 and 2 are calculated by applying the fundamental principle of the dynamics in rotation and in translation at two time instants t 1, t 2 and the acceleration without braking of a wheel of the vehicle 6 is calculated. . (0.6,. (T2) respectively at each of the two instants ti, t2, according to equations 3 and 4. In a third step 30, the equation 5 corresponding to the difference Ae, is calculated. between the accelerations without braking 6, (0.6,. (t2) calculated previously, and assuming that the instants t2 and t1 are sufficiently close so that the external forces Fext (t), Cf (t) applied on the vehicle do not vary significantly.
On calcule ensuite à l'étape 40, l'estime de la masse M du véhicule par un algorithme de moindres carrés récursif à facteur d'oubli à partir de cette différence d'accélération Enfin, à l'étape 50, on filtre la valeur de la masse M estimée par un filtre passe bas 6. Ainsi, la masse totale du véhicule peut être estimée de manière fiable et rapide.Then, in step 40, the estimate of the mass M of the vehicle is calculated by a recursive least-squares algorithm with forgetting factor from this difference in acceleration Finally, in step 50, the value of mass M estimated by a low pass filter 6. Thus, the total mass of the vehicle can be estimated reliably and quickly.