FR 3010797 A1 20150320 FR 1358978 A 20130918 DOMAINE TECHNIQUE AUQUEL SE RAPPORTE L'INVENTION La présente invention concerne de manière générale la gestion des cycles de charge et décharge d'une batterie d'accumulateurs. Elle concerne plus particulièrement un procédé d'estimation du vieillissement d'au moins une cellule de batterie d'accumulateurs, comportant des étapes : a) d'acquisition d'une tension aux bornes de ladite cellule et d'une intensité du courant traversant ladite cellule, et b) de calcul d'une capacité maximale et d'un niveau de charge de ladite cellule en fonction de la tension et de l'intensité du courant acquises à l'étape a). L'invention s'applique particulièrement avantageusement dans les véhicules automobiles équipés d'un moteur électrique alimenté par une batterie d'accumulateurs dite de traction. ARRIERE-PLAN TECHNOLOGIQUE De manière bien connue, la puissance électrique que peut fournir une batterie d'accumulateurs diminue au cours d'un cycle de décharge. La quantité maximale d'énergie que peut stocker la batterie d'accumulateurs diminue quant à elle progressivement au cours de la durée de vie de cette batterie. Afin de tirer au mieux partie de la puissance électrique disponible dans la batterie pour alimenter un moteur électrique de véhicule automobile, afin de prédire l'autonomie résiduelle du véhicule avant une nouvelle charge de batterie, et afin de prévoir à quel moment il sera nécessaire de changer toute ou partie de la batterie, on cherche à déterminer deux paramètres particuliers de la batterie, que sont son niveau de charge noté SOC (acronyme anglais de "State Of Charge") et sa capacité maximale notée Q. Le niveau de charge de la batterie, qui s'exprime généralement en pourcentage, traduit le niveau de chargement de la batterie entre un niveau de charge minimal où la batterie n'est plus utilisable, et un niveau de charge maximal. La capacité maximale de la batterie, qui s'exprime généralement en ampère-heure, permet de connaître la durée pendant laquelle la batterie peut fournir un courant électrique d'un ampérage donné. Cette capacité se dégrade au cours du temps, en fonction notamment de l'historique de température de la batterie et de son historique de cycles de charge et décharge. On connaît alors du document FR2971854 un procédé d'estimation de ce niveau de charge et de cette capacité maximale, utilisant un filtre de Kalman étendu dont la seule relation de base est constituée par l'équation suivante : dSOC I = avec I le courant traversant la batterie d'accumulateurs. dt Q' Les inconvénients majeurs de ce procédé d'estimation sont, d'une première part, son manque de précision et de fiabilité, d'une seconde part, sa non-détection d'une cellule défectueuse au sein de la batterie d'accumulateurs, et, d'une troisième part, son erreur d'interprétation concernant des variations réversibles de capacité considérées à tort comme étant dues au vieillissement. OBJET DE L'INVENTION Afin de remédier aux inconvénients précités de l'état de la technique, la présente invention propose un procédé d'estimation tel que défini dans l'introduction, dans lequel, à l'étape b), le calcul est réalisé par résolution : - d'un système d'équations discrétisées correspondant à une modélisation de ladite cellule sous la forme d'un circuit électrique comportant, en série, une source de tension parfaite, une résistance, un premier couple de résistance et de capacité branchées en parallèle, et au moins un second couple de résistance et de capacité branchées en parallèle, et - d'une équation discrétisée d'estimation de la variation de la capacité maximale de ladite cellule. Selon l'invention, la cellule est donc modélisée par un circuit électrique. La source de tension parfaite modélise alors le potentiel électrochimique de la cellule. La différence de potentiels aux bornes de cette source de tension parfaite dépend alors du niveau de charge de la cellule. La résistance modélise la chute de tension induite par la connectique de la cellule et la résistance interne de la cellule. Dans le premier couple de résistance et de capacité branchées en parallèle, la résistance modélise le phénomène de transfert de charge et la capacité représente le phénomène de double couche aux interfaces électrode/électrolyte. On appelle « phénomène de transfert de charge » le courant résultant des réactions d'oxydation et de réduction aux électrodes au sein de la cellule. On appelle « phénomène de double couche » le phénomène de séparation des charges ayant lieux à l'interface entre l'électrode et l'électrolyte. Les autres couples de résistance et de capacité branchées en parallèle modélisent le phénomène de diffusion des ions lithium dans la cellule, c'est-à-dire la polarisation due au déplacement des ions dans l'électrolyte et à l'intérieur des électrodes. Ainsi, les états de charge et capacités maximales sont calculées de manière fiable, au moyen d'un algorithme représentatif de la réalité électrochimique d'une cellule de batterie, qui converge rapidement vers des valeurs proches de la réalité. Par ailleurs, grâce à l'invention, il est possible d'estimer la capacité maximale de chaque cellule de la batterie d'accumulateurs. Dès lors, une comparaison de ces différentes capacités maximales permettra de détecter une cellule défectueuse au sein de la batterie d'accumulateurs. D'autres caractéristiques avantageuses et non limitatives du procédé d'estimation conforme à l'invention sont les suivantes : - ledit système d'équations discrétisées s'exprime sous la forme suivante : SOC (k)= SOC (k - 1) + (1+ 3,0).1(k -1).T O 1(k -1).T, Vdr(k)- (1 & 1)+ TU 1(k -1).T Tyk)= (1 1) + , R avec = 1 C (k) = E,q(SOC (k)) + RQ.1(k)+17,(k)+ ,=1 (k) avec Te une période d'échantillonnage considérée, k un indice incrémental, SOC le niveau de charge de ladite cellule, I l'intensité du courant traversant ladite cellule, 4Q étant un paramètre relatif à une capacité nominale en début de vie et à la capacité maximale instantanée de ladite cellule, Voeu I la tension aux bornes du circuit électrique, Vdi la tension aux bornes du premier couple de résistance et capacité, dont les valeurs sont respectivement notées Rct et Cdi, et V1 la tension aux bornes du second couple de résistance et capacité, dont les valeurs sont respectivement notées R1 et C1 ; - ladite équation discrétisée d'estimation de la variation de la capacité maximale de ladite cellule s'exprime en fonction d'un paramètre relatif à une capacité nominale en début de vie et à la capacité maximale instantanée de ladite cellule ; - ladite équation discrétisée d'estimation de la variation de la capacité maximale de ladite cellule s'exprime sous la forme AO(k) = AQ(k -1), k étant un indice incrémental ; - ladite équation discrétisée d'estimation de la variation de la capacité maximale de ladite cellule s'exprime sous la forme AQ(k) = g (SOC (k AQ(k - 0, 1(k - 0, T(k -1)) , avec k un indice incrémental, g une fonction représentative de la dynamique de vieillissement de la cellule, SOC le niveau de charge de ladite cellule, I l'intensité du courant traversant ladite cellule, T une température mesurée sur ladite cellule ; - à l'étape b), le calcul est réalisé en résolvant lesdites équations à l'aide d'un observateur d'état dont le vecteur d'état comporte le niveau de charge de ladite cellule, une variable représentative de la perte de capacité de ladite cellule, la tension aux bornes du premier couple de résistance et capacité, et la tension aux bornes de chaque second couple de résistance et capacité ; - à l'étape b), le calcul est réalisé en résolvant lesdites équations à l'aide d'un observateur d'état dont le vecteur d'entrée comprend la tension aux bornes du circuit électrique ; - il est prévu une étape supplémentaire c) de calcul d'un indicateur d'une perte d'autonomie de ladite cellule, en fonction de la capacité maximale de ladite cellule calculée à l'étape b) ; - à l'étape c), la valeur dudit indicateur est modifiée seulement si l'estimation de la capacité maximale de ladite cellule est considérée comme 25 convergée ; - à l'étape c), la valeur dudit indicateur est modifiée seulement si la température mesurée de ladite batterie d'accumulateurs est comprise dans la plage de température de caractérisation en début de vie de ladite batterie d'accumulateurs ; 30 - à l'étape c), la valeur dudit indicateur d'état de santé est modifiée seulement si le niveau de charge de ladite batterie d'accumulateurs est compris dans une plage de caractérisation en début de vie de ladite batterie d'accumulateurs ; - à l'étape c), la valeur dudit indicateur d'état de santé est modifiée seulement si la tension acquise à l'étape a) est comprise dans une plage de tension de caractérisation en début de vie de ladite batterie d'accumulateurs ; et - à l'étape c), la valeur dudit indicateur d'état de santé est modifiée seulement si l'intensité acquise à l'étape a) est comprise dans une plage d'intensité de caractérisation en début de vie de ladite batterie d'accumulateurs. DESCRIPTION DETAILLEE D'UN EXEMPLE DE REALISATION La description qui va suivre en regard des dessins annexés, donnés à titre d'exemples non limitatifs, fera bien comprendre en quoi consiste l'invention et comment elle peut être réalisée. Sur les dessins annexés : - la figure 1 est une vue schématique d'un véhicule automobile équipé d'une batterie d'accumulateurs comprenant une pluralité de cellules, et d'une unité de pilotage adaptée à mettre en oeuvre un procédé d'estimation du vieillissement de ces cellules ; - la figure 2 est un schéma électrique modélisant l'une des cellules de la batterie d'accumulateurs de la figure 1 ; et - la figure 3 est un graphique illustrant, en pointillés, les variations de l'estimation de la capacité maximale de la cellule de la figure 2, en trait interrompu, les variations d'un indicateur intermédiaire d'état de santé, et, en trait continu, les variations d'un indicateur d'état de santé filtré. En préliminaire, dans la suite de cet exposé, on entendra par « batterie d'accumulateurs » un élément apte à stocker de l'énergie électrique lorsqu'il est alimenté en courant par un réseau électrique extérieur, puis à restituer ultérieurement cette énergie électrique. Une telle batterie d'accumulateurs peut être par exemple de type électrochimique (par exemple Lithium-ion) ou d'un type différent (par exemple condensateur). Sur la figure 1, on a représenté très schématiquement un véhicule automobile 1. Ce véhicule automobile 1 est ici un véhicule électrique. Il comporte alors un moteur électrique 10 prévu pour entraîner en rotation ses roues motrices 20. En variante, il pourrait s'agir d'un véhicule hybride, comportant un moteur à combustion interne et un moteur électrique pour l'entraînement de ses roues motrices. Comme le montre la figure 1, le véhicule automobile 1 comporte une batterie d'accumulateurs, appelée batterie de traction 40. Cette batterie de traction 40 est ici exclusivement prévue pour alimenter le moteur électrique 10 en courant. En variante, cette batterie de traction pourrait être également prévue pour alimenter en courant différents appareils électriques consommateurs de courant, tels que le système de direction assistée, le système de climatisation... Cette batterie de traction 40 comporte un boîtier extérieur 43 duquel émergent deux bornes, une positive 44 et l'autre négative 45 connectées au moteur électrique 10 via un boîtier électronique de puissance (non représenté). La batterie de traction 40 comporte également une pluralité de cellules 41 qui sont logées dans le boîtier extérieur 43 et qui sont ici branchées en série entre les deux bornes positive 44 et négative 45. En variante, on pourrait prévoir que les cellules soient connectées par couples en parallèles, et que ces couples de cellules soient branchés en série entre les deux bornes positive et négative. Le nombre de cellules 41 utilisées est déterminé de telle sorte que le moteur électrique 10 puisse développer un couple et une puissance suffisants pour propulser le véhicule automobile 1 pendant une durée prédéterminée. Une cellule de batterie de traction délivrant habituellement une tension de l'ordre de 3 à 5 V, le nombre de cellules 41 est alors calculé de telle manière que la tension aux bornes de la batterie de traction 40 puisse atteindre 400 V. Chaque cellule 41 de la batterie de traction 40 sera alors ici surveillée indépendamment des autres cellules 41. On caractérise à cet effet ici chaque cellule 41 par trois paramètres particuliers appelés niveau de charge SOC, capacité maximale Q- et indicateur d'état de santé SOH_E. Le niveau de charge SOC s'exprime en pourcentage. Il traduit l'état de chargement de la cellule 41 considérée, entre un niveau de charge minimal où la batterie n'est plus utilisable (0 %), et un niveau de charge maximal (100 %). La capacité maximale Q- s'exprime en ampère-heure. Elle indique la durée pendant laquelle la cellule 41 peut fournir un courant électrique d'un ampérage donné. Au moment de la fabrication de la cellule 41, cette capacité maximale Q, est généralement égale ou légèrement inférieure à la capacité nominale Qn pour laquelle la cellule a été conçue. Puis, elle se dégrade au cours 3010 797 7 du temps, en fonction notamment de l'historique de température de la cellule 41 et de son historique de cycles de charge et décharge. L'indicateur d'état de santé SOH_E s'exprime en pourcentage. Il fournit une information quant au niveau de vieillissement de la cellule 41 considérée. En 5 règle générale, au moment de la fabrication de la cellule 41, l'indicateur d'état de santé SOH_E est égal ou légèrement inférieur à 100%, puis il diminue avec le temps, selon l'utilisation qui est faite de la cellule 41. Il s'avère alors essentiel, et c'est d'ailleurs l'objet de la présente invention, de déterminer précisément ces trois paramètres afin de surveiller le 10 vieillissement de la batterie de traction 40 en vue de prévenir toute panne du véhicule automobile 1. Pour cela, le véhicule automobile 1 comporte un calculateur 30, qui est ici représenté comme étant indépendant de la batterie de traction 40. En variante, ce calculateur pourrait être intégré à la batterie. Selon une autre variante, le 15 calculateur pourrait faire partie intégrante de l'unité globale de pilotage du moteur électrique. Ici, comme le montre la figure 1, le calculateur 30 comprend un processeur (CPU), une mémoire vive (RAM), une mémoire morte (ROM), des convertisseurs analogiques-numériques (A/D), et différentes interfaces d'entrée et 20 de sortie. Grâce à ses interfaces d'entrée, le calculateur 30 est adapté à recevoir de différents capteurs des signaux d'entrée relatifs à la batterie de traction 40. Dans sa mémoire vive, le calculateur 30 mémorise ainsi en continu : - la température T de la batterie de traction 40, mesurée ici à l'aide d'un 25 capteur de température situé dans le boîtier extérieur 43 de la batterie de traction 40, - l'intensité I du courant débité par la batterie de traction 40, à l'aide d'un capteur de courant connecté électriquement entre les cellules 41 et la borne négative 45 de la batterie de traction 40, et 30 - la tension Voeu aux bornes de chaque cellule 41, à l'aide de voltmètres connectés électriquement aux bornes des cellules 41 de la batterie de traction 40. Grâce à un logiciel mémorisé dans sa mémoire morte, le calculateur 30 est adapté à déterminer, en fonction des valeurs mesurées, le niveau de charge SOC, la capacité maximale Q,, et l'indicateur d'état de santé SOH_E de chaque cellule 41. Enfin, grâce à ses interfaces de sortie, le calculateur 30 est adapté à transmettre ces données à l'unité globale de pilotage du moteur électrique 10. Ici, le fonctionnement du logiciel sera exposé en relation avec une seule cellule 41 de la batterie de traction 40. En pratique, il sera mis en oeuvre de la même manière sur chacune des autres cellules de la batterie de traction 40. Au démarrage, le calculateur 30 met en oeuvre une opération d'initialisation au cours de laquelle il affecte des « valeurs initiales estimées » à différents paramètres, notamment au niveau de charge SOC, à la capacité maximale Q- et à l'indicateur d'état de santé SOH_E. Ces valeurs initiales peuvent par exemple être choisies égales aux valeurs calculées au précédent cycle de fonction du moteur électrique 10. Au premier démarrage du moteur, elles peuvent par exemple être choisies ainsi : - SOC = 100 %, - = Qr,, et - SOH_E = 100%. Le calculateur 30 met ensuite en oeuvre un algorithme en quatre étapes, répétées de manière récurrente à chaque pas de temps (le pas de temps considéré étant ici noté k). La première étape est une étape d'acquisition des paramètres de la cellule 41 considérée de la batterie de traction 40. La seconde étape est une étape de calcul de la capacité maximale Qr et du niveau de charge SOC de la cellule 41. La troisième étape est une étape de validation des données calculées. La quatrième étape est une étape de calcul de l'indicateur d'état de santé SOH_E de la cellule 41 et d'un indicateur d'état de santé SOH_E de la batterie de traction 40. Au cours de la première étape, le calculateur 30 acquiert les valeurs de la tension Voeu I aux bornes de la cellule 41 considérée, de l'intensité I du courant traversant cette cellule 41, et de la température T de la batterie de traction 40. La seconde étape, qui consiste à calculer la capacité maximale Q- et le niveau de charge SOC de la cellule 41, est ici opérée au moyen d'un observateur d'état. Un tel observateur d'état est utilisé pour reconstruire, à partir des mesures, les variables internes d'un système dynamique. L'observateur, à partir d'un modèle embarqué de la cellule et de la valeur du courant d'entrée prédit la tension de cette cellule. Il compare ensuite cette prédiction à la mesure de tension de la cellule. La différence entre tension prédite et tension mesurée est utilisée pour adapter les états internes du modèle et les faire converger de manière à annuler l'écart entre tension prédite et tension mesurée. Ici, selon une caractéristique particulièrement avantageuse de l'invention, le calcul de la capacité maximale C), et du niveau de charge SOC de la cellule 41 est réalisé par résolution : - d'un système d'équations discrétisées correspondant à une modélisation de la cellule 41 sous la forme d'un circuit électrique 42 comportant, en série, une source de tension parfaite Eeq, une résistance RQ, un premier couple de résistance Re et de capacité Cdi branchées en parallèle, et au moins un second couple de résistance R1 et de capacité C; branchées en parallèle, et - d'une équation discrétisée d'estimation de la variation de la capacité maximale Q- de la cellule 41. Théoriquement, ces équations ont été élaborées de la manière suivante. Comme le montre la figure 2, la cellule 41 est modélisée sous la forme d'un circuit électrique 42 comportant des composants idéaux. Dans ce modèle, la source de tension parfaite Eeq modélise le potentiel électrochimique de la cellule. La différence de potentiels aux bornes de cette source de tension parfaite Eeq dépend alors directement du niveau de charge SOC de la cellule 41. La résistance RQ modélise la chute de tension induite par la connectique de la cellule 41 et la résistance interne de la cellule 41. Dans le premier couple de résistance Rd et de capacité Cd branchées en parallèle, la résistance Re modélise le phénomène de transfert de charge et la capacité Cdi représente le phénomène de double couche aux interfaces électrode/électrolyte. On note Vdi la tension aux bornes de ce premier couple. Il est ici prévu N couples de résistance R1 et de capacité C1 branchées en parallèle, avec i allant de 1 à N. Ces N couples modélisent le phénomène de diffusion des ions lithium dans la cellule 41 et dans les électrodes de la cellule 41. On note V1 la tension aux bornes du ième couple. N peut par exemple être choisi égal à 2. La mise en équation de ce circuit électrique 42 équivalent conduit au système d'équations suivant : dSOC I dt I7 di V 'di dt R,C Cd, dI' 17 I ' + , pour i = 1... iNT L dt R,C C VCell = Eeq (SOC )+ Rç2.1 +V cldl+2_, r--1V L'équation d'estimation de la variation de la capacité maximale Q- de la cellule 41 est quant à elle obtenue sur la base du constat suivant : tandis que le niveau de charge SOC évolue rapidement au cours d'un cycle de charge ou de décharge, que la tension Vdi peut varier de OV à quelques dixièmes de volts en quelques secondes, et que la polarisation de diffusion y N V peut passer de OV à quelques dixièmes de volts en quelques minutes, la capacité maximale Q- de la cellule 41 varie très lentement (en moyenne, de l'ordre de 20 % en 5 ans). De ce fait, l'hypothèse de considérer une variation de capacité maximale Q- nulle entre deux pas de temps est raisonnable. On écrit alors : 614Q = 0 , avec AC) =- 1 . dt Qr On notera ici que l'utilisation de ce paramètre 4Q plutôt que la capacité maximale Q- améliorera le fonctionnement de l'algorithme et évitera des problèmes numériques pouvant entraîner sa divergence. La discrétisation de ces équations permet d'obtenir le système suivant : SOC (k) = SOC (k - 1) + (1+ 40).1(k -1).T, Q' 40(k) = .40(k - 1) T 1(k -1).T, (k ) - (1 (k 1) + RaC C dl V(k)=(1- e ),I'(k 1) +1(k -1).T, avec = 1... N R C V,,11(k) = E ,' (SOC (k)) + RÇ.2.1(k ) + (k)+ (k) Dans ce système d'équations, Te est la période d'échantillonnage. L'observateur d'état ici utilisé pour obtenir à chaque pas de temps k une bonne évaluation de la capacité maximale Q- et du niveau de charge SOC de la cellule 41 est un filtre de Kalman étendu. En pratique, le calculateur 30 met en oeuvre une pluralité de calculs basés sur ce filtre de Kalman étendu afin d'obtenir des estimations de la capacité maximale Q, et du niveau de charge SOC de la cellule 41. Ces calculs sont connus de l'homme du métier et ne seront donc pas décrits ici en détail. On pourra, pour plus de détail, se référer par exemple à l'ouvrage « Optimal State Estimation » de Dan Simon, publié aux éditions Wiley. En résumé, cet observateur d'état utilise un vecteur d'état xk, un vecteur d'entrée uk, un vecteur de sortie yk, un bruit d'état wk, et un bruit de mesure vk. Il considère alors le système stochastique suivant : fk -I Xk -1 - I (I) k - I ) Yk = h, (x,'v,,) k(o,k) V ,(0,R,) Dans ce système, Rk est la matrice de covariance de bruit de mesure et Qk la matrice de covariance de bruit d'état. En pratique, le vecteur d'état xk est défini ainsi : SOC (k)- AQ(k) Vcii(k) V1(k) VN(k) - Le vecteur de sortie yk est défini ainsi : Yk = Vcell(k) Le vecteur d'entrée uk est défini ainsi : uk =1(k) Les fonctions fk-1 et hk sont obtenues à l'aide du système d'équations discrétisées mentionné supra. Le calculateur 30 linéarise alors le système d'équations en calculant les matrices jacobiennes suivantes : Xk = - afk ax F i54-1 Lk-1 x-k_, Il estime alors la matrice de covariance d'erreur d'estimation de la variable d'état Pk et la variable d'état xk de la manière suivante : P, = Fk IPkIFkTI+ LkIQkILI = f, 1(x k_i,u k_1,0) Le calculateur 30 linéarise alors ces équations en calculant les matrices jacobiennes suivantes : H h_ ah ax ah av - - Enfin, le calculateur 30 met à jour le gain de Kalman Kk, estime la variable d'état xk et la matrice de covariance d'erreur d'estimation de la variable d'état Pk de la manière suivante : K k PT (H ,P, I I ,T T + AI kR,AI , )1 = , H, xk = +K JyA - h(x,'u A,0)] = (1 - pkII k )1'; Le calculateur 30 obtient ainsi les valeurs de la capacité maximale Q- et du niveau de charge SOC de la cellule 41. La précision de la modélisation permet le choix d'un paramétrage 15 dynamique de l'observateur conduisant à une convergence rapide. La troisième étape est une étape de validation des données calculées. Cette étape consiste à vérifier, d'une part, si la valeur de la capacité maximale Q- calculée a suffisamment convergé pour être utilisable, et, d'autre part, si les conditions d'utilisation de la cellule 41 sont proches des conditions de 20 sa caractérisation en début de vie (lorsque la capacité maximale Q, de la cellule était considérée égale à la capacité nominale Qn). C'est-à-dire que nous disposons effectivement de données de comparaison permettant de conclure sur la perte irréversible de capacité. Pour vérifier que l'algorithme a suffisamment convergé, le calculateur 30 détermine la différence entre la dernière capacité ao maximale Q,-(k) calculée, et la capacité maximale Q,-(k-k') calculée quelques pas de temps auparavant, puis il estime que la capacité maximale Q- a suffisamment convergé si cette différence est inférieure à un seuil SO prédéterminé. Ici, le calculateur 30 considère que la capacité maximale Q- a 5 suffisamment convergé si : Q,(k)-Q,(k-10) I Sur la figure 3, on a représenté les convergences de l'estimation de la capacité maximale Qr, de l'indicateur intermédiaire de vieillissement SOH_Eint, et de l'indicateur d'état de santé SOH_E. On observe sur cette figure que dans un premier temps (partie A), la capacité maximale Q, n'a pas encore convergé, si bien que l'indicateur de validité est maintenu égal à 0. De ce fait, les indicateurs de vieillissement SOH_Ebr et SOH_E ne varient pas. Dans un second temps (partie B), le calculateur 30 estime que la capacité maximale Qr a convergé et affecte donc la valeur 1 à l'indicateur de validité Ô. L'indicateur intermédiaire de vieillissement SOH_E1 varie alors brusquement pour converger vers la valeur du rapport Q,-/Qn. Grâce au filtre, l'indicateur d'état de santé SOH_E varie quant à lui moins brusquement, pour converger progressivement vers la valeur de ce rapport (partie C). Enfin, le calculateur 30 détermine l'indicateur d'état de santé SOH_E de la batterie de traction 40. Pour cela, le calculateur 30 affecte à l'indicateur d'état de santé SOH_E de la batterie de traction 40 la valeur de l'indicateur d'état de santé SOH E de la cellule 41 qui a le plus vieilli, selon la formule suivante : SOH_Ebatterie = mirl(SOH_Ecellules). A ce stade, le calculateur 30 peut aussi déterminer de manière fiable si l'une des cellules 41 est défectueuse. Pour cela, il compare les valeurs des indicateurs de vieillissement SOH E de chaque cellule avec une valeur attendue. Cette valeur attendue peut par exemple être constituée par la moyenne des indicateurs de vieillissement SOH E des cellules 41. Le calculateur 30 en déduit un dysfonctionnement si cet écart dépasse une valeur prédéterminée, par exemple égale à 30%. La présente invention n'est nullement limitée au mode de réalisation décrit et représenté, mais l'homme du métier saura y apporter toute variante 5 conforme à son esprit. Ainsi, on pourrait par exemple prévoir que la température soit mesurée à l'aide d'un capteur placé contre la face externe du boîtier extérieur de la batterie de traction. On pourrait aussi prévoir d'utiliser plusieurs capteurs de température afin d'obtenir une meilleure estimation de la température au niveau de la cellule 10 considérée. Selon une autre variante, on aurait pu prévoir de considérer que l'indicateur d'état de santé de la cellule soit constitué par sa capacité maximale. Encore en variante, on aurait pu utiliser un paramètre différent de 4Q pour calculer le niveau de charge SOC et la capacité maximale Q-. Ainsi, on aurait 15 par exemple pu considérer un paramètre 0 égal au rapport de la capacité nominale Q, sur la capacité maximale Q- de la cellule. Alors, l'équation discrétisée d'estimation de la variation de la capacité maximale Q- de la cellule 41 se serait exprimée sous la forme suivante : 13(k) = 13(k-1), avec p(k)=Q,/Q,(k). 20 Selon une autre variante, on aurait pu utiliser le paramètre 4Q, en estimant toutefois que ce dernier varie légèrement entre deux pas de temps, en fonction par exemple de l'un ou l'autre des paramètres que sont le niveau de charge SOC, l'intensité I et la température T de la cellule 41 considérée. Dans cette variante, l'équation discrétisée d'estimation de la variation de 25 la capacité maximale de la cellule se serait alors exprimée sous la forme suivante : AO(k) = g (SOC (k - 0, AQ(k - 1), I(k -1).T (k -1)), avec g une fonction représentative de la dynamique de variation de la capacité de la cellule. REVENDICATIONS1. Procédé d'estimation du vieillissement d'au moins une cellule (41) de batterie d'accumulateurs (40), comportant des étapes : a) d'acquisition d'une tension (Voeu') aux bornes de ladite cellule (41) et d'une intensité (I) du courant traversant ladite cellule (41), et b) de calcul d'une capacité maximale (C),) et d'un niveau de charge (SOC) de ladite cellule (41) en fonction de la tension (Voeu) et de l'intensité (I) acquises à l'étape a), caractérisé en ce que, à l'étape b), le calcul est réalisé par résolution : - d'un système d'équations discrétisées correspondant à une modélisation de ladite cellule (41) sous la forme d'un circuit électrique (42) comportant, en série, une source de tension parfaite une résistance (R1-2), un premier couple de résistance (R0) et de capacité (Cd1) branchées en parallèle, et au moins un second couple de résistance (R) et de capacité (C1) branchées en parallèle, et - d'une équation discrétisée d'estimation de la variation de la capacité maximale (Q,-) de ladite cellule (41). 2. Procédé d'estimation selon la revendication précédente, dans lequel ledit système d'équations discrétisées s'exprime sous la forme suivante : SOC (k)= SOC (k - 1) + (1+ 3,Q).1(k O 1(k -1).TeT 17d1(k) - (1 ).Vai(k 1) + 1?,,C C V,(k)= (1 ).V (k 1) + 1(k -1).T , avec I = 1... N C avec : - Te une période d'échantillonnage considérée, - k un indice incrémental, - SOC le niveau de charge (SOC) de ladite cellule (41), - I l'intensité (I) du courant traversant ladite cellule (41), - 4Q une fonction d'une capacité nominale (Qn) en début de vie de lacellule (41) et la capacité maximale (0,-) instantanée de la cellule (41), - Voeu l la tension aux bornes du circuit électrique (42), - Vdi la tension aux bornes du premier couple de résistance (Rd) et capacité (Cd1), dont les valeurs sont respectivement notées Re et Cdi, et - V, la tension aux bornes du second couple de résistance (R) et capacité (C1), dont les valeurs sont respectivement notées R1 et C. 3. Procédé d'estimation selon l'une des revendications précédentes, dans lequel ladite équation discrétisée d'estimation de la variation de la capacité maximale (0,-) de ladite cellule (41) s'exprime en fonction d'un paramètre (4Q) relatif à une capacité nominale (Q,-,) en début de vie de ladite cellule (41) et à la capacité maximale (Q,-) instantanée de ladite cellule (41). 4. Procédé d'estimation selon la revendication 3, dans lequel ladite équation discrétisée d'estimation de la variation de la capacité maximale (Q,) de ladite cellule (41) s'exprime sous la forme suivante : AO(k) = AO(k -1), k étant un indice incrémentai. 5. Procédé d'estimation selon la revendication 3, dans lequel ladite équation discrétisée d'estimation de la variation de la capacité maximale (Q,) de ladite cellule (41) s'exprime sous la forme suivante : AQ (k) = g (SOC (k - 1). AQ (k - 1), 1(k -1),T (k - 1)) , avec : - k un indice incrémental, - g une fonction représentative de la dynamique de variation de la capacité de la cellule (41), - SOC le niveau de charge (SOC) de ladite cellule (41), - I l'intensité (I) du courant traversant ladite cellule (41), - T une température mesurée sur ladite cellule (41). 6. Procédé d'estimation selon l'une des revendications précédentes, dans lequel, à l'étape b), le calcul est réalisé en résolvant lesdites équations à l'aide d'un observateur d'état dont le vecteur d'état (xk) comporte : - le niveau de charge (SOC) de ladite cellule (41), - la capacité maximale (Q) de ladite cellule (41), - la tension (Vd1) aux bornes du premier couple de résistance (Re) et capacité (Cd1), et - la tension (V1) aux bornes de chaque couple i de résistance (R) etcapacité (C1). 7. Procédé d'estimation selon l'une des revendications précédentes, dans lequel, à l'étape b), le calcul est réalisé en résolvant lesdites équations à l'aide d'un observateur d'état dont le vecteur de sortie (yk) comprend la tension (Vali) aux bornes du circuit électrique (42). 8. Procédé d'estimation selon l'une des revendications précédentes, dans lequel il est prévu une étape supplémentaire c) de calcul d'un indicateur d'état de santé (SOH_E) de ladite cellule (41), en fonction de la capacité maximale (Q,-) de ladite cellule (41) calculée à l'étape b). 9. Procédé d'estimation selon la revendication précédente, dans lequel, à l'étape c), la valeur dudit indicateur d'état de santé (SOH_E) est modifiée seulement si la variation de la capacité maximale (Qr) de ladite cellule (41) est inférieure à un seuil prédéterminé. 10. Procédé d'estimation selon l'une des deux revendications précédentes, dans lequel, à l'étape c), la valeur dudit indicateur d'état de santé (SOH_E) est modifiée seulement si la température mesurée (T) de ladite batterie d'accumulateurs (40) est comprise dans une plage de température de caractérisation en début de vie (Tmin,carac 11. Procédé d'estimation selon l'une des trois revendications précédentes, dans lequel, à l'étape c), la valeur dudit indicateur d'état de santé (SOH_E) est modifiée seulement si le niveau de charge (SOC) de ladite batterie d'accumulateurs (40) est compris dans une plage de caractérisation en début de vie (S0Cmin,carac 12. Procédé d'estimation selon l'une des quatre revendications précédentes, dans lequel, à l'étape c), la valeur dudit indicateur d'état de santé (SOH_E) est modifiée seulement si la tension (Voeu') acquise à l'étape a) est comprise dans une plage de tension de caractérisation en début de vie (Vmin,carac 13. Procédé d'estimation selon l'une des cinq revendications précédentes, dans lequel, à l'étape c), la valeur dudit indicateur d'état de santé (SOH_E) est modifiée seulement si l'intensité (I) acquise à l'étape a) est comprise dans une plage d'intensité de caractérisation en début de vie (Imin,carac