9 89420 1 PROCEDE DE CALCUL AUTOMATIQUE D'UNE MASSE DE SUIES PRESENTES DANS UN FILTRE À PARTICULES [0001] L'invention concerne un procédé de calcul automatique d'une masse de suies présentes dans un filtre à particules. L'invention concerne également un procédé de régénération du filtre à particules. L'invention concerne enfin un calculateur électronique programmé pour mettre en oeuvre le procédé de calcul. [0002] On connaît un filtre à particules pour épurer des gaz d'échappement émis par un moteur à combustion interne d'un véhicule automobile. Un tel filtre à particules est configuré pour piéger des particules contenues dans ces gaz d'échappement, avant que ces gaz ne soient rejetés dans l'environnement. Ces particules comprennent par exemple des suies issues d'une combustion de carburant dans le moteur. Lors du fonctionnement du véhicule, ces suies piégées s'accumulent dans le filtre à particules, et doivent être éliminées lors d'opérations de régénération du filtre à particules. Pour éviter d'endommager le filtre à particules, la masse de suies doit être maintenue en deçà d'une certaine valeur. [0003] Typiquement, cette masse de suies est estimée à l'aide d'une cartographie, en fonction de paramètres physiques mesurés, tels que des débits et des pressions des gaz d'échappement traversant le filtre à particules. Il est nécessaire d'intégrer numériquement plusieurs paramètres avant de les appliquer à l'entrée de la cartographie. [0004] Cette méthode présente de nombreux inconvénients. Notamment, l'utilisation de cette cartographie nécessite une étape de calibration préalable, longue et compliquée à mettre en place. De plus, le temps de calcul de l'intégration numérique préalable à l'utilisation de la cartographie dépend des valeurs des paramètres physiques mesurés. En particulier, pour une cartographie donnée, des valeurs de paramètres mesurés se trouvant en dehors d'une certaine gamme de valeurs peuvent conduire à des temps de calcul élevés et le modèle peut s'avérer sensible à une non linéarité entre les paramètres mesurés. [0005] II existe donc un besoin pour un procédé de calcul d'une masse de suies présentes dans un filtre à particules, qui présente une exécution rapide, une précision accrue et une calibration simplifiée. [0006] L'invention vise à résoudre un ou plusieurs de ces inconvénients. L'invention porte ainsi sur un procédé de calcul automatique d'une masse de suies présentes dans un filtre à particules d'une ligne d'échappement du moteur à combustion interne d'un véhicule automobile, ce procédé comprenant une succession d'itérations de calcul espacées dans le temps, chaque itération permettant le calcul de la masse mt de suies instantanée en fonction de la valeur de la masse mt_At de suies calculée lors d'une itération précédente, dans lequel chaque itération comprend : -la mesure d'une différence de pression instantanée de gaz d'échappement APmes entre une entrée et une sortie du filtre à particules ; -le calcul, par un calculateur électronique, de la différence de pression An -,calc, la différence de pression An étant calculée par une fonction polynomiale -,calc d'un débit volumique Q, dont des coefficients sont des fonctions rationnelles, c'est-à-dire de rapports de fonctions polynômes, de la masse de suies, en fonction du débit volumique Q instantané des gaz d'échappement et de la masse mt_At de suies calculée lors de l'itération précédente ; -le calcul de la masse mt de suies instantanée par intégration numérique de la fonction F(AP mes - AP calc(nit-At)) , où F est une fonction telle que les dérivées dF/dm et dF/dQ sont nulles. [0007] Selon une variante, l'intégration comprend le calcul de la quantité mt -int_3, = F (A P ni' - A P.i, (mt_axAt où At est un intervalle de temps entre deux itérations consécutives. [0008] Selon une variante, l'intervalle de temps At séparant deux itérations (21, 22) consécutives est compris entre 10ms et 500ms. [0009] Selon encore une variante, la fonction F est une fonction mathématique qui, à toute valeur x réelle, associe la valeur K*x, où K est une constante. [0010] Selon une autre variante, la valeur du coefficient K est comprise entre 0,01 et 100. [0011] Selon encore une autre variante, chaque itération comprend la mémorisation de la valeur de la masse mt de suies instantanée calculée lors de ladite itération. [0012] Selon une variante, le procédé comprend, lors d'une première exécution du procédé, une étape d'initialisation de la valeur de la masse de suies à une valeur prédéfinie. [0013] Selon une autre variante, la différence de pression est calculée par la fonction AP./, (in) = aAi (m) x Q2 + Pm (m) x Q+ Yci(m)xQ + aA2 x Q2 dans laquelle Q est le débit volumique des gaz d'échappement, aA2 est un coefficient réel et les coefficients aA1(m) , 13B1(m) , yC1(m) sont des fonctions rationnelles monotones de la masse de suie m. [0014] L'invention porte par ailleurs sur un procédé de régénération d'un filtre à particules, comprenant : 2 9 89420 3 -le calcul d'une masse de suies présentes dans le filtre à particules ; -la régénération du filtre à particules, lorsque la masse de suies calculée est supérieure ou égale à un seuil ; -la masse de suies est calculée selon un procédé de calcul automatique tel 5 que défini ci-dessus. [0015] L'invention porte en outre sur un calculateur électronique pour la mise en oeuvre d'un procédé de calcul automatique tel que défini ci-dessus, le calculateur électronique est programmé pour exécuter une succession d'itérations de calcul espacées dans le temps, chaque itération permettant le 10 calcul de la masse mt de suies instantanée en fonction de la valeur de la masse mt-At de suies calculée lors d'une itération précédente, chaque itération comprenant : -la mesure d'une différence de pression instantanée de gaz d'échappement Apmes entre une entrée et une sortie du filtre à particules ; 15 -le calcul de la différence de pression Apcalc, la différence de pression Apcalc étant calculée par une fonction polynomiale d'un débit volumique Q, dont des coefficients sont des fonctions rationnelles de la masse de suies, en fonction du débit volumique Q instantané des gaz d'échappement et de la masse mt-At de suies calculée lors de l'itération précédente ; 20 -le calcul de la masse mt de suies instantanée par intégration numérique de la fonction F(AD mes - A P calc(int-At)) , où F est une fonction telle que les dérivées dF/dm et dF/dQ sont nulles. [0016] D'autres caractéristiques et avantages de l'invention ressortiront clairement de la description qui en est faite ci-après, à titre indicatif et 25 nullement limitatif, en référence aux dessins annexés, dans lesquels : - la figure 1 illustre un véhicule automobile comprenant un filtre à particules ; - la figure 2 est un organigramme d'un procédé de calcul d'une masse de suies présentes dans le filtre à particules de la figure 1 ; - la figure 3 est un organigramme d'une itération du procédé de la figure 2 ; 30 - la figure 4 est un organigramme d'un procédé de régénération du filtre à particules de la figure 1. [0017] Dans la suite de cette description, les fonctions et les caractéristiques bien connues de l'homme de métier ne sont pas décrites en détail. [0018] L'invention propose un procédé de calcul d'une masse de suies 35 piégée par un filtre à particules. [0019] La figure 1 représente un véhicule automobile 2, comprenant : -un moteur 4 à combustion interne, configuré pour mouvoir le véhicule 2 ; -une ligne d'échappement 6, incluant un filtre à particules 8 ; -un calculateur électronique 10. [0020] La ligne d'échappement 6 est reliée fluidiquement à un orifice 7 d'évacuation de gaz d'échappement du moteur 4. Le filtre à particules 8 est configuré pour épurer des gaz d'échappement circulant dans la ligne 6. [0021] Le filtre 8 comprend : -une entrée et une sortie, adaptées pour permettre l'entrée et la sortie de gaz d'échappement issus de l'orifice 7 ; -une pluralité de canaux, s'étendant entre l'entrée et la sortie. [0022] Ces canaux sont configurés pour permettre l'écoulement des gaz d'échappement entre l'entrée et la sortie. Chaque canal comprend une extrémité d'entrée et une extrémité de sortie, respectivement, reliée fluidiquement à l'entrée et la sortie du filtre 8. Ces canaux sont séparés les uns des autres par des parois poreuses. Ces parois sont configurées pour filtrer les gaz d'échappement et, notamment, pour piéger des particules de suies contenues dans ces gaz. Ici, les canaux sont identiques entre eux. [0023] Le calculateur 10 est programmé pour exécuter un procédé de calcul de la masse de suies contenues dans le filtre 8. A cet effet, le calculateur 10 est configuré pour recevoir des signaux représentatifs de pressions des gaz d'échappement en entrée et en sortie du filtre 8. Ces pressions sont ici mesurées, respectivement, par des capteurs de pression 12 et 14. [0024] Le fonctionnement du calculateur 10 va maintenant être décrit, en référence au procédé de la figure 2. [0025] Ce procédé comprend une succession d'itérations de calcul 21, 22 espacées dans le temps. Chaque itération 21, 22 permet de calculer une valeur instantanée d'une masse de suies contenues dans le filtre 8, en fonction d'une valeur de la masse de suies calculée lors d'une itération précédente. Pour simplifier la figure 2, seules deux itérations 21, 22 sont illustrées. Ces itérations 21, 22 sont identiques ; aussi, dans la suite de cette description, seule l'itération 22 sera décrite en détail. [0026] Avantageusement, le procédé comprend une étape d'initialisation 20 de la masse de suies à une valeur prédéfinie mo, représentative d'une masse de suies connue dans le filtre 8. La masse de suies est ici initialisée, préalablement à la première utilisation du filtre 8, à une valeur mo par exemple inférieure à cinq grammes et, de préférence, égale à zéro. [0027] La figure 3 représente plus en détail l'itération de calcul 22. Cette itération 22 comprend : -une étape 30 d'acquisition d'une valeur de la masse mt_At de suies calculée lors d'une itération immédiatement précédente ; -une étape 32 de mesure d'une différence de pression instantanée de gaz d'échappement Apmes entre l'entrée et la sortie du filtre 8 ; 2 9 89420 5 -une étape 34 de calcul de la différence de pression Apcalc, calculée en fonction de la valeur de la masse mt_At de suies calculée lors de l'itération précédente ; -une étape 36 de calcul d'une valeur de masse mt de suies instantanée, en 5 fonction des différences de pression mesurée et calculée. [0028] La valeur initiale de la masse acquise lors de l'étape 30 est, par exemple, la valeur de la masse mt_At de suies calculée lors de l'itération précédente. [0029] Lors de l'étape 32, la mesure de la différence de pression instantanée 10 APmes des gaz d'échappement est ici réalisée en mesurant les pressions de gaz d'échappement en entrée et en sortie du filtre 8, respectivement, à l'aide des capteurs de pression 12 et 14. [0030] Au cours de l'étape 34, la différence de pression An -r-calc des gaz entre cette entrée et cette sortie est calculée. Ce calcul est réalisé grâce à une 15 modélisation du filtre 8 en fonction de paramètres physiques des gaz d'échappement. La différence de pression An est modélisée par une fonction polynomiale du débit volumique des d'échappement. Les coefficients de cette fonction polynomiale sont des fonctions rationnelles de la masse de suies. Cette différence de pression An -r-calc est calculée en fonction du 20 débit volumique Q instantané des gaz d'échappement et de la masse mt_At de suies calculée lors de l'itération précédente. Ce débit volumique Q est par exemple mesuré ou calculé par un estimateur. [0031] Ici, cette différence de pression An est exprimée en fonction d'une masse m de suies, par la formule suivante : 25 AP'/, (m) = aAi (m) Q2+ Am (m) Q+2/ci(m)xQ+ aA2 x Q2 dans laquelle Q est le débit volumique instantané des gaz d'échappement, aA2 est un coefficient réel et les coefficients aAl(m), PB1(m), Yci(m) sont des fonctions rationnelles monotones dépendantes d'une masse de suie m. [0032] Cette formulation particulière de Apcalc(m) découle d'un modèle 30 physique de la perte de charge des gaz d'échappement lors de leur circulation entre l'entrée et la sortie du filtre 8. Ce modèle physique permet notamment de calculer An (m) en fonction d'une masse m de suies contenues dans le filtre 8. [0033] Dans ce modèle, la masse m de suies contenues dans le filtre 8 est 35 considérée comme formant une couche d'épaisseur E(m), déposée sur des parois des canaux du filtre 8. L'épaisseur E(m) dépend de la masse m, par le truchement de la fonction : E(m)= m x 2 1x ^ 100- %Résidu Pd orNbCANAUX 4x(LFAp - LpLuG )X 100 X LciN 2 9 89420 6 dans laquelle m est exprimée en kilogrammes, NbcANAUX est le nombre de canaux du filtre 8 ; (LFAP - LpLuG) est une longueur du filtre 8 sur laquelle les gaz sont filtrés ; LciN est la largeur d'un canal, mesurée sur son extrémité d'entrée ; 0 .-depot est la masse volumique de la couche de suies et %Résidu est 5 un taux d'encrassement du filtre 8. Ce taux d'encrassement %Résidu est, par exemple, défini comme étant le rapport entre la surface des parois recouverte par des cendres (ce qui reste des suies après une régénération) et la surface totale des parois. Ce taux %Résidu est compris entre 0% et 100%. [0034] Ici, le nombre de canaux NbcANAUX est inférieur ou égal à 5000. La 10 longueur (LFAP - LPLUG) est comprise entre 0,1m et 1m et, de préférence, entre 0,2m et 0,5m. La largeur LciN est inférieure ou égale à lcm, 0,5cm ou 0,2cm. La masse volumique n 1- depot de la couche de suies est inférieure ou égale à 500 kg/m3 ou 200 kg/m3. [0035] Les coefficients am (m) et aA2 modélisent des pertes de charges des 15 gaz d'échappement, respectivement, lors de leur entrée et de leur sortie du filtre 8, du fait de changements de section. Ces coefficients sont donnés, respectivement, par les formules suivantes : «Ai (m) = Alx P gaz 20 P gaz A2 x ceA2(m) = 4 LCOUT où Al et A2 sont des constantes sans unités, pgaz est la masse volumique des gaz d'échappement et LcouT est la largeur d'un canal, mesurée sur l'extrémité 25 de sortie de ce canal. [0036] Ici, les constantes Al et A2 sont choisies pour assurer l'homogénéité des unités entre des grandeurs physiques intervenant dans le calcul de la différence de pressions. Al et A2 sont compris entre 10-14 et 10-13. La largeur l_cOUT est inférieure ou égale à 1 cm, 0,5cm ou 0,2cm. 30 [0037] Le coefficient 13131(m) modélise des pertes de charge des gaz d'échappement lors de leur écoulement le long des canaux du filtre 8, selon un modèle basé sur l'équation de Bernoulli. Ce coefficient 13,31(m) est ici donné par la formule suivante : r 35 /3131(m) = (LC/N 2XE(M))4 B1 x,te )x LFAP L pLUG X / 2 ( + LCOUT 4 COUT LCOUT + (I, cm - 2x E kin ff )) 2 T 100 - %Résidu , LFAP X LPLUG I \\ 2 100 (LaN - 2 E (m ll ,i x (LOIN 2 xE(m)/4 LCOUT 2 + (Law 2 xE(m))2 2 LCOUT où B1 est une constante sans unité. Ici, B1 est choisie pour assurer l'homogénéité des unités entre des grandeurs physiques intervenant dans le calcul de la différence de pressions. B1 est, par exemple, compris entre 10-6 et 105 [0038] Le coefficient yci(m) modélise des pertes de charge des gaz d'échappement lors de leur filtrage par les parois poreuses du filtre 8, selon un modèle basé sur la formule de Darcy. Ce coefficient yci(m) est ici donné par la formule suivante : 1 du (T )x (ePAROI + E(m))Yci(m) = x Cl Kp (E (m))x (LFAp - LpL')x (L cm, - 2 x E(m)) où ePAROI est l'épaisseur des parois poreuses ; Kp(E(m)) est un coefficient de perte de charge dépendant de la perméabilité des suies, calculé en fonction de l'épaisseur E(m) ; C1 est une constante sans unité et p(T) est la viscosité dynamique instantanée des gaz d'échappement, calculée en fonction d'une valeur instantanée de la température T des gaz d'échappement dans le filtre 8. Cette température T est par exemple mesurée dans la ligne d'échappement par un capteur de température, non représenté ou calculée à partir d'un estimateur. [0039] Ici, l'épaisseur ePAROI est comprise entre 10-4m et 10-3m. La constante C1 est choisie pour assurer l'homogénéité des unités entre des grandeurs physiques intervenant dans le calcul de la différence de pressions. C1 est, par exemple, compris entre 10-7 et 10-6. [0040] Le coefficient Kp(E(m)) est, par exemple, donné par la formule suivante : Kp (E (m)) = (e PAROI + E(M))X Ksic X K sooT X C pe ea (E(M)) K SOOT X C pe ea (E (m))x e pARoi + K , x E (m) où Ksour, Ksic sont les perméabilités intrinsèques des matériaux constituant, respectivement, la suie et la paroi d'un canal ; Cpermea(E(m)) est un coefficient de perméabilité de la couche de suies. [0041] Ici, la perméabilité Ksour est comprise entre 10-15m2 et 10-13m2 ; la perméabilité Ksic est comprise entre 10-13m2 et 10-11m2. [0042] Le coefficient Cpermea(E(m)) est ici calculé en fonction de l'épaisseur E(m) de suies, par la formule suivante : E(m) eoX +c0 Dpart Cpermea(E(m)) = doxE(m)+1 si E(m) 0.2 e Dp,t Dpart si E(m) < 0.2 Dpart a0 x E(m) + bo Dpart OÙ Dpart est une constante et ao, bo, co, do et eo sont des constantes sans unité. [0043] Ici, ces constantes sont choisies empiriquement pour définir le coefficient de perméabilité. Dpart est compris entre 10-6m et 10-4m. ao est compris entre 1 et 10. bo est compris entr0e 0,05 et 0,3. co est compris entre -5 et -0,5. do est compris entre 1 et 10. eo est compris entre -10-1 et -10-3. [0044] Avantageusement, l'étape 34 comprend la détermination de la masse volumique instantanée pgaz des gaz d'échappement. Par exemple, cette masse volumique pgaz est donnée par la formule pgaz = Dm/Q, où Dm est le débit massique des gaz d'échappement. Ce débit massique Dm est, par exemple, inférieur à 1000kg/h. [0045] L'étape 34 comprend également la détermination de la viscosité dynamique instantanée p(T) des gaz d'échappement, en fonction d'une valeur mesurée de la température T. Par exemple, cette viscosité p(T), exprimée en Pa.s, est donnée par la formule : ,u(T)= -1,07e-6 x T2 +5,009e-3 x T + 0,44 où la température T est exprimée en Kelvin. Cette température T est, par exemple, inférieure à 2500 K. [0046] Ainsi, avec ce choix de modèle physique, le filtre 8, permet de modéliser avec précision la circulation des gaz, quelle que soit la structure du filtre 8. [0047] Lors de l'étape 36, la masse mt de suies instantanée est calculée, à partir des valeurs des différences de pressions Apcaic et Apmes. Cette masse mt est calculée par intégration numérique de l'équation suivante : dm dt = F(AP.' - AP'ic(mt_A M où F est une fonction telle que les dérivées dF/dm et dF/d0 sont nulles. [0048] Cette intégration comprend ici le calcul de la quantité suivante : nit - m,-A, = F (A P., - APc.,(mt_At)))< Ar où At est un intervalle de temps entre deux itérations consécutives. L'intervalle At est, par exemple, au moins égal à 10ms, ou 25ms. At est, par exemple, au plus égal à 500ms ou 400ms. De préférence, At est égal à 100ms. [0049] La masse de suies est calculée à partir d'un modèle physique du filtre 35 8, qui relie analytiquement une masse m de suies à une différence de pression Apcab(m) de gaz d'échappement entre une entrée et une sortie du filtre 8. L'utilisation de ce modèle permet de simplifier le calcul de la masse de suies par rapport à la lecture d'une cartographie nécessitant plusieurs intégrations en entrée.30 [0050] L'intégration d'une quantité dépendante de la différence entre An et Apmes(m) permet d'assurer un calcul rapide de la masse de suies, ce que ne permettrait pas une inversion littérale du modèle physique An (m) En effet, les coefficients am(m), PB1(m) et Nice) sont des fonctions rationnelles en m de degré élevé. Leur inversion en temps réel, pour le calcul littéral, à chaque instant, de la masse m en fonction de Q et de Apmes, nécessiterait un temps de calcul élevé et demanderait des ressources informatiques importantes. [0051] Dans cet exemple, la fonction F est une fonction mathématique qui, à toute valeur x réelle, associe la valeur K*x, où K est une constante ne dépendant pas du débit volumique Q ni de la masse mt de suies. Ce coefficient K est, par exemple, supérieur ou égal à 0,01m.s ou 100m.s. [0052] Ainsi, la constante K joue le rôle d'un gain dans un système de régulation à boucle fermée, ce qui permet d'améliorer la vitesse de convergence du calcul de la masse mt et de réduire la sensibilité de ce calcul à des conditions initiales. [0053] Ici, la masse calculée lors de l'étape 36 est mémorisée par le calculateur 10. Ainsi, lorsque le procédé est redémarré à la suite d'une interruption, par exemple après un arrêt de fonctionnement du moteur, le procédé peut, lors de l'étape 30 suivant ce redémarrage, acquérir la dernière valeur calculée de la masse de suies. [0054] Avantageusement, lors d'une étape 40, illustrée à la figure 4, le filtre 8 est régénéré si la masse de suies calculée est supérieure ou égale à un seuil. [0055] De nombreux autres modes de réalisation sont possibles. [0056] La détermination de la viscosité dynamique instantanée p(T) peut être omise. La détermination de la masse volumique instantanée pgaz des gaz d'échappement peut également être omise. [0057] En variante, la mesure des pressions de gaz d'échappement en entrée et en sortie du filtre 8 est remplacée par l'acquisition de valeurs de pressions fournies par un estimateur. [0058] La masse de suies peut également être initialisée à une valeur m0 postérieurement à une régénération du filtre 8. FIELD OF THE INVENTION The invention relates to a method for automatically calculating a mass of soot present in a particulate filter. The invention also relates to a method of regeneration of the particulate filter. The invention finally relates to an electronic calculator programmed to implement the calculation method. [0002] A particulate filter is known for purifying exhaust gases emitted by an internal combustion engine of a motor vehicle. Such a particulate filter is configured to trap particles contained in these exhaust gases before these gases are released into the environment. These particles include, for example, soot from fuel combustion in the engine. During the operation of the vehicle, these trapped soot accumulate in the particulate filter, and must be removed during regeneration operations of the particulate filter. To prevent damage to the particulate filter, the soot mass must be kept below a certain value. Typically, this mass of soot is estimated using a map, according to measured physical parameters, such as flow rates and pressures of the exhaust gas passing through the particle filter. It is necessary to numerically integrate several parameters before applying them to the input of the cartography. This method has many disadvantages. In particular, the use of this mapping requires a preliminary calibration step, long and complicated to set up. In addition, the calculation time of the digital integration prior to the use of the map depends on the values of the physical parameters measured. In particular, for a given map, measured parameter values lying outside a certain range of values can lead to high computation times and the model may be sensitive to non-linearity between the measured parameters. There is therefore a need for a method of calculating a mass of soot present in a particle filter, which has a fast execution, increased accuracy and simplified calibration. The invention aims to solve one or more of these disadvantages. The invention thus relates to a method for automatically calculating a mass of soot present in a particulate filter of an exhaust line of the internal combustion engine of a motor vehicle, this method comprising a succession of iterations of computation spaced in time, each iteration allowing the calculation of the mass and instantaneous soot as a function of the value of the mass mt_At of soot computed during a previous iteration, in which each iteration comprises: the measurement of a difference instantaneous pressure of exhaust gas APmes between an inlet and an outlet of the particulate filter; the calculation, by an electronic calculator, of the pressure difference An -, calc, the pressure difference An being calculated by a polynomial function -, calc of a volume flow rate Q, whose coefficients are rational functions, and that is, polynomial function ratios, soot mass, as a function of the instantaneous flow rate Q of the exhaust gas and the mass mt_At of soot computed during the previous iteration; the calculation of the instantaneous mass of soot by numerical integration of the function F (AP mes - AP calc (nit-At)), where F is a function such that the derivatives dF / dm and dF / dQ are zero. According to a variant, the integration comprises the calculation of the quantity mt -int_3, = F (AP n '- A Pi, (mt_axAt where At is a time interval between two consecutive iterations. [0008] According to a variant the interval of time At between two consecutive iterations (21, 22) is between 10ms and 500ms. According to another variant, the function F is a mathematical function which, at any real value x, associates the value K * x, where K is a constant According to another variant, the value of the coefficient K is between 0.01 and 100. According to yet another variant, each iteration comprises the storage of the value of the instantaneous mass of soot calculated during said iteration According to a variant, the method comprises, during a first execution of the method, a step of initializing the value of the mass of soot to a predefined value. According to another variant, the pressure difference is driven by AP./, (in) = aAi (m) x Q2 + Pm (m) x Q + Yci (m) xQ + aA2 x Q2 where Q is the volume flow of the exhaust gas, aA2 is a real coefficient and coefficients aA1 (m), 13B1 (m), yC1 (m) are monotonic rational functions of soot mass m. The invention furthermore relates to a method for regenerating a particle filter, comprising: the calculation of a mass of soot present in the particulate filter; the regeneration of the particulate filter, when the calculated mass of soot is greater than or equal to a threshold; the mass of soot is calculated according to an automatic calculation method as defined above. The invention also relates to an electronic computer for implementing an automatic calculation method as defined above, the electronic computer is programmed to execute a succession of calculation iterations spaced in time each iteration allowing the calculation of the instantaneous mass of soot as a function of the value of the soot mass mt-At calculated during a previous iteration, each iteration comprising: the measurement of an instantaneous pressure difference of exhaust gas Apmes between an inlet and an outlet of the particulate filter; The calculation of the pressure difference Apcalc, the pressure difference Apcalc being calculated by a polynomial function of a volume flow rate Q, whose coefficients are rational functions of the soot mass, as a function of the instantaneous volume flow rate Q of the exhaust gas and mt-At soot mass calculated during the previous iteration; The calculation of the mass of instantaneous soot by numerical integration of the function F (AD mes -A P calc (int-At)), where F is a function such that the derivatives dF / dm and dF / dQ are zero. Other features and advantages of the invention will emerge clearly from the description which is given below, by way of indication and in no way limiting, with reference to the accompanying drawings, in which: - Figure 1 illustrates a vehicle automobile comprising a particulate filter; FIG. 2 is a flowchart of a method for calculating a soot mass present in the particle filter of FIG. 1; FIG. 3 is a flow diagram of one iteration of the method of FIG. 2; FIG. 4 is a flowchart of a regeneration method of the particle filter of FIG. 1. In the remainder of this description, the functions and characteristics well known to those skilled in the art are not described. in detail. The invention proposes a method for calculating a soot mass 35 trapped by a particulate filter. Figure 1 shows a motor vehicle 2, comprising: an internal combustion engine 4, configured to move the vehicle 2; an exhaust line 6, including a particulate filter 8; an electronic computer 10. The exhaust line 6 is fluidly connected to an exhaust gas discharge port 7 of the engine 4. The particulate filter 8 is configured to purify exhaust gases flowing in line 6. [0021] The filter 8 comprises: an inlet and an outlet adapted to allow the entry and the exit of exhaust gas issuing from the orifice 7; a plurality of channels, extending between the input and the output. These channels are configured to allow the flow of exhaust gas between the inlet and the outlet. Each channel comprises an inlet end and an outlet end, respectively, fluidly connected to the inlet and the outlet of the filter 8. These channels are separated from each other by porous walls. These walls are configured to filter the exhaust gas and, in particular, to trap soot particles contained in these gases. Here, the channels are identical to each other. The computer 10 is programmed to execute a method of calculating the mass of soot contained in the filter 8. For this purpose, the computer 10 is configured to receive signals representative of the exhaust gas pressures at the input and in the These pressures are here measured, respectively, by pressure sensors 12 and 14. The operation of the computer 10 will now be described, with reference to the method of FIG. 2. [0025] This method comprises a succession of calculation iterations 21, 22 spaced apart in time. Each iteration 21, 22 makes it possible to calculate an instantaneous value of a mass of soot contained in the filter 8, as a function of a value of the mass of soot calculated during a previous iteration. To simplify FIG. 2, only two iterations 21, 22 are illustrated. These iterations 21, 22 are identical; also, in the remainder of this description, only the iteration 22 will be described in detail. Advantageously, the method comprises a step 20 of initializing the soot mass to a predefined value mo, representative of a soot mass known in the filter 8. The mass of soot is here initialized, prior to the first use of the filter 8, at a value mo for example less than five grams and preferably equal to zero. FIG. 3 represents in greater detail the calculation iteration 22. This iteration 22 comprises: a step 30 of acquiring a value of the mass mt_At of soot computed during an immediately preceding iteration; a step 32 of measuring an instantaneous pressure difference of the exhaust gas Apmes between the inlet and the outlet of the filter 8; A step 34 of calculating the pressure difference Apcalc, calculated as a function of the value of the mass mt_At of soot computed during the previous iteration; a step 36 of calculating a mass value of instantaneous soot, as a function of the measured and calculated pressure differences. The initial value of the mass acquired during step 30 is, for example, the value of the mass mt_At soot computed during the previous iteration. In step 32, the measurement of the instantaneous pressure difference 10 AP of the exhaust gas is here carried out by measuring the exhaust gas pressures at the inlet and at the outlet of the filter 8, respectively, at Using the pressure sensors 12 and 14. In step 34, the pressure difference Δn -r-calc of the gases between this input and this output is calculated. This calculation is carried out by modeling the filter 8 as a function of physical parameters of the exhaust gases. The pressure difference An is modeled by a polynomial function of the volume flow rate of the exhaust. The coefficients of this polynomial function are rational functions of the soot mass. This difference in pressure An -r-calc is calculated as a function of the instantaneous flow rate Q of the exhaust gas and of the mass mt_At of soot computed during the previous iteration. This volume flow Q is for example measured or calculated by an estimator. Here, this difference in pressure An is expressed as a function of a mass m of soot, by the following formula: AP '/, (m) = aAi (m) Q2 + Am (m) Q + 2 / ci (m) xQ + aA2 x Q2 where Q is the instantaneous flow rate of the exhaust gases, aA2 is a real coefficient and the coefficients aAl (m), PB1 (m), Yci (m) are monotonic rational functions dependent on a mass of soot m. This particular formulation of Apcalc (m) derives from a physical model 30 of the pressure drop of the exhaust gases during their circulation between the inlet and the outlet of the filter 8. This physical model makes it possible in particular to calculate An (m) as a function of a mass m of soot contained in the filter 8. In this model, the mass m of soot contained in the filter 8 is considered to form a layer of thickness E (m) , deposited on the walls of the channels of the filter 8. The thickness E (m) depends on the mass m, by means of the function: E (m) = m × 2 × 100-% Pd or NbCANAUX residual 4x (LFAp - LpLuG) X 100 X LciN 2 9 89420 6 where m is expressed in kilograms, NbcANAUX is the number of channels of the filter 8; (LFAP - LpLuG) is a length of the filter 8 on which the gases are filtered; LciN is the width of a channel, measured on its input end; The deposit is the density of the soot layer and the residue is a fouling rate of the filter 8. This fouling rate% Residue is, for example, defined as being the ratio between the covered wall surface area. by ashes (what remains soot after regeneration) and the total surface of the walls. This% Residual rate is between 0% and 100%. Here, the number of NbcANAUX channels is less than or equal to 5000. The length (LFAP - LPLUG) is between 0.1m and 1m and, preferably, between 0.2m and 0.5m. The LciN width is less than or equal to 1cm, 0.5cm or 0.2cm. The density n-1 of the soot layer is less than or equal to 500 kg / m3 or 200 kg / m3. The coefficients am (m) and aA2 model the pressure losses of the exhaust gases, respectively, at their entry and at their exit from the filter 8, due to section changes. These coefficients are given, respectively, by the following formulas: "Ai (m) = Alx P gas 20 P gas A2 x ceA2 (m) = 4 LCOUT where Al and A2 are constants without units, pgaz is the density of the gases Exhaust and LcouT is the width of a channel, measured on the output end of this channel. Here, the constants Al and A2 are chosen to ensure the homogeneity of the units between physical quantities involved in the calculation of the pressure difference. Al and A2 are from 10-14 to 10-13. The width l_cOUT is less than or equal to 1 cm, 0.5 cm or 0.2 cm. The coefficient 13131 (m) models pressure losses of the exhaust gases as they flow along the channels of the filter 8, according to a model based on the Bernoulli equation. This coefficient 13.31 (m) is given here by the following formula: ## EQU1 ## (LC / N 2XE (M)) 4 B1 x, te) x LFAP L pLUG X / 2 (+ LCOUT 4 COST LCOUT + (I, cm - 2x E kin ff)) 2 T 100 -% Residue, LFAP X LPLUG I \\ 2 100 (LaN - 2E (m ll, ix (LOIN 2 xE (m) / 4 LCOUT 2 + (Law 2 xE (m)) 2 2 LCOUT where B1 is a constant without unit Here, B1 is chosen to ensure the homogeneity of units between physical quantities involved in the calculation of the pressure difference. for example, between 10-6 and 105 The coefficient yci (m) models exhaust pressure losses of the exhaust gases when they are filtered by the porous walls of the filter 8, according to a model based on the Darcy formula. This coefficient yci (m) is here given by the following formula: 1 of (T) x (ePAROI + E (m)) Yci (m) = x Cl Kp (E (m)) x (LFAp - LpL ') x (L cm, - 2 x E (m)) where ePAROI is the thickness of the porous walls, Kp (E (m)) is a coefficient of pressure loss depending on the permeability of the substrates. ies, calculated according to the thickness E (m); C1 is a unitless constant and p (T) is the instantaneous dynamic viscosity of the exhaust gas, calculated as a function of an instantaneous value of the temperature T of the exhaust gas in the filter 8. This temperature T is for example measured in the exhaust line by a temperature sensor, not shown or calculated from an estimator. Here, the thickness ePAROI is between 10-4m and 10-3m. The constant C1 is chosen to ensure the homogeneity of the units between physical quantities involved in the calculation of the pressure difference. C1 is, for example, between 10-7 and 10-6. The coefficient Kp (E (m)) is, for example, given by the following formula: ## EQU1 ## M)) K SOOT XC pe ea (E (m)) xe pARoi + K, x E (m) where Ksour, Ksic are the intrinsic permeabilities of the materials constituting, respectively, the soot and the wall of a channel; Cpermea (E (m)) is a coefficient of permeability of the soot layer. Here, the permeability Ksour is between 10-15m2 and 10-13m2; Ksic permeability is between 10-13m2 and 10-11m2. The coefficient Cpermea (E (m)) is here calculated as a function of the thickness E (m) of soot, by the following formula: E (m) eoX + c0 Dpart Cpermea (E (m)) = doxE (m) +1 if E (m) 0.2 e Dp, t Start if E (m) <0.2 Start a0 x E (m) + bo Start where Ù Start is a constant and ao, bo, co, do and eo are constant without unity. Here, these constants are chosen empirically to define the coefficient of permeability. Departure is between 10-6m and 10-4m. ao is from 1 to 10, and is from 0.05 to 0.3. co is between -5 and -0.5. do is between 1 and 10. eo is between -10-1 and -10-3. Advantageously, step 34 comprises the determination of the instantaneous density pgaz of the exhaust gas. For example, this density pgaz is given by the formula pgaz = Dm / Q, where Dm is the mass flow rate of the exhaust gas. This mass flow Dm is, for example, less than 1000 kg / h. Step 34 also comprises determining the instantaneous dynamic viscosity p (T) of the exhaust gas as a function of a measured value of the temperature T. For example, this viscosity p (T), expressed as Pa.s, is given by the formula:, u (T) = -1,07e-6 x T2 + 5,009e-3 x T + 0,44 where the temperature T is expressed in Kelvin. This temperature T is, for example, less than 2500 K. Thus, with this choice of physical model, the filter 8 makes it possible to accurately model the flow of gases, whatever the structure of the filter 8. [0047] In step 36, the instantaneous mass of soot is calculated from the values of the pressure differences Apcaic and Apmes. This mass mt is calculated by numerical integration of the following equation: ## EQU1 ## where F is a function such that the derivatives dF / dm and dF / d0 are zero. This integration here includes the calculation of the following quantity: nit-m, -A, = F (AP, - APc., (Mt_At))) <Ar where At is a time interval between two consecutive iterations. The interval At is, for example, at least equal to 10ms, or 25ms, At is, for example, at most equal to 500ms or 400ms, Preferably At is equal to 100ms, the soot mass is calculated to from a physical model of the filter 8, which analytically links a mass m of soot to a pressure difference Apcab (m) of exhaust gas between an inlet and an outlet of the filter 8. The use of this model allows to simplify the calculation of the mass of soot compared to the reading of a cartography requiring several input integrations. [0050] The integration of a quantity endante of the difference between An and Apmes (m) makes it possible to ensure a quick calculation of the mass of soot, which would not allow a literal inversion of the physical model An (m) Indeed, the coefficients am (m), PB1 (m) and Nice) are rational functions in m of high degree. Their inversion in real time, for the literal calculation, at any moment, of the mass m as a function of Q and Apmes, would require a high computation time and would require important computer resources. In this example, the function F is a mathematical function which, at any real value x, associates the value K * x, where K is a constant that does not depend on the volume flow rate Q nor the mass mt of soot. This coefficient K is, for example, greater than or equal to 0.01m.s or 100m.s. Thus, the constant K plays the role of a gain in a closed-loop control system, which makes it possible to improve the convergence speed of the calculation of the mass mt and to reduce the sensitivity of this calculation to initial conditions. Here, the mass calculated during step 36 is stored by the computer 10. Thus, when the process is restarted following an interruption, for example after an engine running stoppage, the method can, in step 30 following this restart, acquire the last calculated value of the mass of soot. Advantageously, during a step 40, illustrated in Figure 4, the filter 8 is regenerated if the calculated soot mass is greater than or equal to a threshold. Many other embodiments are possible. The determination of the instantaneous dynamic viscosity p (T) can be omitted. The determination of the instantaneous density of the exhaust gases can also be omitted. In a variant, the measurement of the exhaust gas pressures at the inlet and the outlet of the filter 8 is replaced by the acquisition of pressure values provided by an estimator. The mass of soot can also be initialized to a value m0 subsequent to a regeneration of the filter 8.