Brevet FR3031814A1 — publié le 20160722

PROCEDE D'IDENTIFICATION DE LA COURBE DE TENSION A VIDE D'UNE CELLULE ELECTRIQUE EN VIEILLISSEMENT Le domaine de l'invention concerne un procédé 5 d'identification de la courbe d'une tension à vide d'une cellule d'accumulation d'énergie électrique en vieillissement. Les véhicules hybrides et électriques puisent leur énergie de traction dans des stockeurs électriques embarqués. Une gestion efficiente de l'énergie électrique est primordiale 10 pour réduire le coût du véhicule. Cela permet de dimensionner le stockeur pour correspondre au plus proche des besoins électriques du véhicule et de configurer le dispositif de pilotage du stockeur électrique (communément appelé « Battery Management System » en anglais) avec des lois de pilotage 15 efficientes et préservant les propriétés électriques des stockeurs durant l'utilisation du véhicule. Plus précisément, le stockeur électrique d'un véhicule automobile est généralement composé de plusieurs cellules d'accumulation d'énergie électrique. On définit dans la suite 20 de la description le stockeur électrique d'un véhicule comme l'ensemble des cellules reliées entre elles pour obtenir la tension et la capacité désirée pour la fourniture d'énergie du véhicule. Une cellule d'accumulation d'énergie est caractérisée 25 par des indicateurs de fonctionnement, notamment une résistance interne, une tension à ses bornes, des limites de tension de fonctionnement, une capacité, une tension à vide (communément appelée « Open Circuit Voltage » en anglais), un état de charge (communément appelé « State of Charge » en 30 anglais), et la température de fonctionnement. Par extension, le stockeur électrique est caractérisé par ces mêmes indicateurs. Le dispositif de pilotage du stockeur électrique a pour tâche principale de surveiller les indicateurs de fonctionnement électrique et thermique du stockeur électrique pour assurer un usage sécurisé et un fonctionnement optimal du véhicule. Parmi ces indicateurs de fonctionnement, la courbe de tension à vide est d'une grande importance car elle est indispensable pour déterminer l'état de charge du stockeur électrique ainsi que d'autres indicateurs du dispositif de pilotage. La tension à vide d'un accumulateur électrique correspond à la tension existante aux bornes de l'accumulateur électrique en absence de tout courant électrique à ses bornes et au bout d'un temps suffisamment long de sorte que la tension n'évolue plus. Ainsi, on définit généralement la courbe de tension à vide d'un accumulateur comme la représentation de la tension à vide pour l'ensemble des états de charge de l'accumulateur entre un état de charge minimum et un état de charge maximum. Par exemple, l'énergie récupérable par le stockeur électrique dans les limites de tension de fonctionnement maximale et minimale est déterminée à partir de l'indicateur de fonctionnement de la courbe de tension à vide. On comprend donc l'importance de cet indicateur pour obtenir une bonne efficience du pilotage du stockeur électrique. Or, l'indicateur de la courbe de tension à vide, ainsi que d'autres indicateurs du stockeur électrique, évolue au cours du temps sous l'effet du vieillissement. Il s'ensuit une dégradation de l'aptitude du stockeur à accumuler et délivrer de l'énergie électrique. Ceci provoque une baisse de performance du véhicule liée à une perte de son rendement. Dans les pires cas, une mauvaise détection du vieillissement, par conséquent de l'état de charge ou de la puissance disponible, peut conduire à une immobilisation du véhicule, voire des risques d'incendie ou d'explosion. Parallèlement à ces contraintes techniques, une mauvaise gestion de ces indicateurs conduit à des surcoûts à l'usage tant pour l'utilisateur, en cas de changement précoce du stockeur électrique, que pour le constructeur qui voit le coût se répercuter sur la prestation de la garantie proposée. Par ailleurs, les paramètres de l'indicateur de la courbe de tension à vide du stockeur électrique d'un véhicule sont propres à chaque cellule du stockeur. Ceci provient des dispersions de fabrication des cellules électrochimiques mais aussi du vieillissement différencié des cellules. Celles-ci ont en effet une histoire propre et donc un vieillissement individuel. Il est par conséquent nécessaire de mettre à jour l'indicateur de la courbe de tension à vide de chaque cellule pour proposer un dispositif de pilotage du stockeur performant. On connaît de l'état de la technique les documents 20 nommés ci-après qui décrivent des propositions techniques pour identifier l'indicateur de la courbe de tension à vide d'un stockeur électrique. Le brevet français FR2964464B1 propose une méthode d'identification de la courbe de tension à vide reposant sur 25 la mesure de la capacité lors du vieillissement du stockeur. Cette mesure peut par exemple être réalisée en garage ou encore lors de la recharge du véhicule par l'utilisateur. Il est également supposé que la courbe de tension à vide ainsi définie est invariante avec le vieillissement à condition de 30 réactualiser la capacité. L'inconvénient majeur de la méthode réside dans la mesure de la capacité. Chaque cellule ayant sa propre courbe de tension à vide il faut mesurer la capacité pour chacune d'elles ce qui nécessite, pour toutes les cellules d'un stockeur électrique, de charger et décharger complètement toutes les cellules. Le temps de mise en oeuvre de la mesure de la capacité peut donc être long et induire des dépenses énergétiques liées à la décharge complète des cellules. Même si une charge et décharge complète est possible dans le cas d'un véhicule électrique ou hybride disposant d'une interface de rechargement du stockeur (véhicule dit « plug-in » en anglais pour interface de rechargement) la méthode reste néanmoins difficile car il faut pouvoir équilibrer en tension toutes les cellules du stockeur à la tension minimale et la tension maximale. Par ailleurs, dans le cas d'un véhicule hybride classique il n'y a pas d'interface de rechargement, ce qui complique la mesure de la capacité. Enfin, l'hypothèse faite sur la non dépendance au vieillissement de la courbe de tension à vide par rapport l'état de charge n'est pas justifiée pour tous les types d'électrochimie. Le brevet français FR2966250B1 a pour objectif de palier les défauts majeurs du brevet décrit dans le paragraphe précédent. La méthode ne cherche plus à mesurer directement la capacité lors du vieillissement du stockeur. La méthode proposée consiste à accéder à la capacité via une cartographie propre à chaque cellule en reliant différentes valeurs de capacité à différentes valeurs de résistance interne. Le défaut majeur de cette méthode est qu'elle suppose la connaissance au préalable du lien à tout niveau de vieillissement entre la valeur de la résistance interne d'une cellule et la valeur de sa capacité réelle. Or, l'existence d'un lien entre la résistance d'une cellule et sa capacité n'est pas justifiée. Il n'a pas été prouvé qu'il y ait une corrélation entre ces deux valeurs. Par ailleurs, même si un lien existe, cette cartographie introduit une imprécision d'estimation. De plus, l'hypothèse faite sur la non dépendance au vieillissement de la courbe de tension à vide par rapport l'état de charge, n'est toujours pas valable pour tous les types d'électrochimie. L'article intitulé « Battery open-circuit voltage estimation by a method of statistical analysis », publié en 2006 par Snihir et accessible dans la publication « Journal of Power Sources, 159(2), pp.1484-1487 », divulgue une estimation de la courbe de tension à vide par approche statistique. La méthode propose d'obtenir la courbe de tension à vide comme la limite, lorsque le courant tend vers zéro, d'un modèle donnant la tension de la cellule en réponse à une sollicitation de type courant constant. Il y est considéré le cas d'une cellule subissant plusieurs cycles de charge et décharge complets pour différentes valeurs de courant décroissantes. Cette méthode présente l'inconvénient d'une procédure de calcul très lourde à implémenter en temps réel. Si la procédure de caractérisation, qui consiste à faire subir au stockeur plusieurs cycles de charge et décharge sous des courants différents, n'est réalisée qu'initialement, on ne voit alors pas apparaître la prise en compte du vieillissement. De plus, si la procédure de caractérisation est réalisée à chaque identification, la méthode devient alors très longue. On ne voit alors plus l'avantage par rapport à une méthode ne nécessitant qu'un seul cycle de charge ou décharge complet. De plus, la méthode serait inapplicable aux véhicules hybrides, quasi dépourvus de cycles de charge/décharge complets. L'article intitulé « Analysis and prediction of the open circuit potential of lithium-ion cells », publié en 2013 par Schmidt et accessible dans la publication « Journal of Power Sources, 239, 696-704 », divulgue une première méthode de paramétrage de l'évolution de la tension à vide en fonction du vieillissement et une deuxième méthode d'identification de la courbe de tension à vide reposant sur la connaissance du potentiel d'équilibre de chaque électrode. Ces méthodes présentent l'inconvénient de la nécessité de connaitre le potentiel d'équilibre de chaque électrode, cette courbe n'étant en général pas fournie par le constructeur du stockeur et encore moins pour chaque cellule. En outre, la mise en oeuvre en temps réel est très difficile. Au regard des documents décrits dans l'état de la technique et des inconvénients de chaque méthode, il existe donc un besoin de proposer une méthode d'identification de la courbe de tension à vide du stockeur électrique pouvant être mise en oeuvre en temps réel et de façon autonome par le calculateur de supervision d'un véhicule. En outre, Il est recherché une méthode d'identification de la courbe de tension à vide pouvant être opérée en cours de vie du véhicule par un diagnostic autonome afin d'adapter les fonctions de commandes du dispositif de pilotage du stockeur électrique. Plus précisément, l'invention concerne un procédé d'identification du vieillissement d'une courbe de tension à vide d'une cellule d'accumulation d'énergie électrique. Selon l'invention, le procédé comprend les étapes successives suivantes : - la détermination d'un polynôme représentant une courbe de tension à vide de référence de la cellule, - le pilotage de la cellule pour maintenir au moins des premier et deuxième paliers de tension à vide à au moins des premier et deuxième états de charge respectivement, - la mesure des première et deuxième tensions à vide pour les premier et deuxième paliers respectivement et d'une 30 première variation de charge de la cellule entre les premier et deuxième paliers, - la détermination d'un coefficient de perte a de la capacité de la cellule à partir d'au moins un premier facteur de perte al calculé selon une fonction de vieillissement du polynôme spécifiquement appliquée avec les première et deuxième tensions à vide et la première variation de charge, - le calcul de la courbe de tension à vide selon la fonction de vieillissement paramétrée avec le coefficient de perte a et appliquée à l'ensemble du polynôme. Avantageusement, la fonction de vieillissement 10 correspond à la formule suivante : SOC - S°Cloou. avec : - U représentant la courbe de tension à vide, - k, un indice de coefficient constituant le polynôme 15 parmi N coefficients du polynôme, - Dkt un coefficient du polynôme et correspondant un indice k, - SOC un état de charge de la cellule, - SOC100%, l'état de charge maximum de la cellule, 20 -a, le coefficient de perte de capacité de la cellule. De préférence, pour le procédé : - le pilotage comprend en outre le maintien à une pluralité de troisièmes paliers de tension à vide, - la mesure comprend en outre la mesure d'une pluralité 25 de troisièmes tensions à vide pour les troisièmes paliers et d'une pluralité de deuxièmes variations de charge de la cellule entre chacun des troisièmes paliers et entre le deuxième palier et au moins un des troisièmes paliers, - la détermination du coefficient de perte de la capacité de la cellule comprend une outre une sélection parmi le premier facteur de perte al et une pluralité de deuxièmes facteurs de pertes at, chaque deuxième facteur de perte at étant calculé identiquement au premier facteur de perte OEl en fonction des troisièmes tensions à vide et des deuxièmes variations de charge. Selon une première variante du procédé, pour la détermination du coefficient de perte a, le calcul du premier facteur de perte al comprend : - le calcul d'une première variation d'état de charge de la cellule entre les première et deuxième tensions à vide 15 mesurées, la première tension à vide mesurée étant la tension à vide pour l'état de charge maximum de la cellule, - le calcul de l'état de charge de la cellule pour la deuxième tension à vide en soustrayant la première variation d'état de charge à l'état de charge maximum de la cellule, 20 - le calcul du premier facteur de perte al selon la tension à vide mesurée et la fonction de vieillissement appliquée pour l'état de charge et pour la deuxième tension à vide de la cellule. Selon une deuxième variante du procédé, pour la 25 détermination du coefficient de perte a, le calcul du premier facteur de perte al comprend: - le calcul d'une deuxième variation d'état de charge de la cellule entre les première et deuxième tensions à vide mesurée, et d'une troisième variation d'état de charge de la 30 cellule entre les tensions à vide correspondante du polynôme, - le calcul du premier facteur de perte al en fonction du rapport entre la deuxième variation et la troisième variation. Selon une variante de déroulement du procédé, le 5 pilotage de la cellule pour maintenir le premier, deuxième et les troisièmes paliers de tension à vide est opéré dans un même cycle de recharge de la cellule. De préférence, le polynôme correspond à la courbe de tension à vide pour un état neuf de la cellule. 10 Avantageusement, au moins les étapes du procédé successives à l'étape de détermination du polynôme sont exécutées lorsque la cellule est installée dans un équipement hôte et lorsque celui-ci fonctionne au moyen de l'énergie électrique fournie par la dite cellule. 15 L'invention prévoit également, un dispositif de pilotage d'un stockeur électrique, le dit stockeur électrique comprenant une pluralité de cellules électriques. Le dispositif comprend : un moyen de calcul opérant le procédé 20 d'identification du vieillissement de la courbe de tension à vide conformément à l'une quelconque des variantes décrites précédemment pour chacune des dites cellules du stockeur, un moyen de détermination des indicateurs de fonctionnement du stockeur électrique, par exemple une courbe 25 de tension à vide du stockeur et un état de charge, à partir des courbes de tension à vide de chacune des dites cellules calculées par le dit procédé d'identification. L'invention prévoit un véhicule automobile comprenant un stockeur électrique, une machine de traction électrique et 30 un dispositif de pilotage du stockeur électrique. Avantageusement, le dispositif de pilotage est conforme au dispositif de pilotage décrit précédemment. Grâce à l'invention, le procédé d'identification est réalisable sur tout type de cellule électrique de technologie 5 lithium-ion sans le besoin de recharger ou décharger complétement la cellule. La connaissance de la courbe de tension initiale, par exemple obtenue par les données de constructeur permet au moyen de la fonction de vieillissement de calculer une courbe de tension à vide d'une cellule en 10 vieillissement par la mesure de deux paliers de tension à vide. Plusieurs paliers permettront d'identifier un vieillissement plus précis. D'autres caractéristiques et avantages de la présente invention apparaîtront plus clairement à la lecture de la 15 description détaillée qui suit de modes de réalisation de l'invention donnés à titre d'exemples nullement limitatifs et illustrés par les dessins annexés, dans lesquels : - la figure 1 représente un schéma d'un équipement électronique comprenant un stockeur électrique composé de 20 cellules d'accumulation d'énergie électrique ; - la figure 2 est un diagramme représentant le procédé de d'identification du vieillissement d'une cellule électrique selon l'invention ; - la figure 3A est un graphique représentant les 25 données mesurées de la cellule électrique lors de l'exécution d'un premier mode de réalisation du procédé d'identification selon l'invention ; - la figure 3B est un graphique représentant les données mesurées de la cellule électrique lors de l'exécution 30 d'un deuxième mode de réalisation du procédé d'identification selon l'invention ; - la figure 4 est un graphique représentant un polynôme de la courbe de tension à vide à neuf de la cellule et les mesures de tension à vide de la cellule lors de l'exécution du deuxième mode de réalisation du procédé, le graphique ayant en abscisse l'état de charge de la cellule et en ordonnée la tension à vide de la cellule ; - la figure 5 représente une formule de la fonction de vieillissement de la courbe de tension à vide à neuf de la cellule et un mode opératoire pour la sélection du coefficient 10 de perte a de capacité de la cellule ; - la figure 6 est un graphique représentant un polynôme de la courbe de tension à vide à neuf de la cellule et la courbe de tension à vide de la cellule en vieillissement. L'invention s'applique à tout équipement électronique, 15 comprenant un stockeur électrique formé d'une ou plusieurs cellules d'accumulation d'énergie électrique. L'invention s'applique de préférence à une cellule électrique de type lithium-ion qui est une technologie requérant un pilotage sécurisé pour éviter les effets de surchauffe ou d'explosion. 20 L'invention propose un procédé d'identification du vieillissement de la courbe de tension à vide d'une cellule électrique dans le but d'améliorer la gestion de l'énergie et proposer un pilotage sécurisé s'adaptant au vieillissement de la cellule électrique. 25 Le procédé selon l'invention peut être exécuté par tout type d'équipement électronique disposant d'un moyen de pilotage du stockeur électrique. Il s'applique donc par exemple à un équipement électronique portable de type téléphone mobile ou ordinateur, à un véhicule automobile 30 électrique ou hybride électrique et ainsi qu'à une installation électrique de grande dimension comprenant un stockeur électrique. Dans la suite de la description, l'invention sera décrite pour un véhicule automobile électrique (ou hybride électrique). On notera que l'une quelconque des variantes du procédé décrit par la suite peut s'appliquer à tout dispositif électronique comprenant un stockeur électrique. L'application à un véhicule automobile n'est pas limitative de la portée de l'invention. La figure 1 est un schéma simplifié représentant un véhicule automobile 11. Le véhicule automobile comporte un stockeur électrique 12 constitué de plusieurs cellules 121 d'accumulation d'énergie électrique de type lithium-ion. Le nombre de cellules 121 est adapté selon les besoins électriques du véhicule 11 et n'est pas un facteur limitatif de l'invention. Le véhicule automobile comprend également un dispositif de pilotage 13 du stockeur électrique 12. Par simplification de représentation, le dispositif de pilotage 13 est le calculateur de supervision du véhicule qui commande également des équipements électroniques 14 et 15. Ces équipements sont par exemple une machine électrique de traction 14 associée à son onduleur 15 et un convertisseur de tension pour un réseau électrique du véhicule. Ce réseau alimente par exemple le réseau de bord. On notera que ce mode de réalisation n'est pas limitatif et d'autres équipements électriques du véhicule non représentés sur la figure 1 sont pilotés par le calculateur de supervision et alimentés par le stockeur électrique 12, ou alimentent le stockeur électrique 12, par exemple un chargeur embarqué qui à partir du réseau électrique 220V recharge le stockeur électrique. Le dispositif de pilotage 13 du stockeur 12 peut également être un calculateur décentralisé au calculateur de supervision et dédié aux fonctions de gestion de l'alimentation électrique et aux fonctions de pilotage du stockeur 12 et de chaque cellule électrique 121. Pour la mise en oeuvre du procédé d'identification de vieillissement de la courbe de tension à vide d'une cellule, 5 le dispositif de pilotage comprend des moyens de mesures de divers indicateurs de fonctionnement de chacune des cellules. Le dispositif de pilotage 13 comprend notamment des capteurs de tension et de courant de chacune des cellules électriques 121. Le dispositif de pilotage 13 comprend également un ou 10 plusieurs calculateurs aptes à calculer des indicateurs de fonctionnement en fonction des données mesurées. Par exemple, partir des données de la variation de la charge d'une cellule, les calculateurs déterminent des variations d'état de charge des cellules électriques et du stockeur dans son 15 intégralité en fonction d'un état de charge de référence. En outre, le dispositif de pilotage 13 comprend des moyens de pilotage des cellules 121 pour les commander en tension et en courant. Ainsi chaque cellule peut être commandée par des fonctions de gestion d'énergie du véhicule 20 pour délivrer de l'énergie électrique à un équipement électronique du véhicule ou pour être rechargée. Les fonctions de gestion sont les fonctions de stratégies énergétiques commandant les différents modes de roulage d'un véhicule hybride. Il s'agit des commandes mettant en fonctionnement la 25 machine électrique de traction pour récupérer de l'énergie cinétique ou pour fournir un couple aux roues. Il peut s'agir également de commandes pour une machine électrique d'aide au démarrage d'un moteur thermique du véhicule. On définit dans la description l'état de charge d'une 30 cellule 121 comme le niveau de charge en énergie électrique de la cellule 121, et l'état de charge du stockeur 12 comme le niveau de charge de la totalité des cellules 121. La variation de charge d'une cellule est obtenue par le calcul de l'intégration entre deux instants donnés de la variation du courant mesuré d'une cellule. Plus précisément, le stockeur électrique 12 comprend 5 une pluralité de cellules électriques 121 agencées pour fournir une capacité d'énergie totale du stockeur 12. Une cellule comprend deux électrodes, une cathode et une anode, appelées bornes dans la suite de la description. La tension à vide d'une cellule est la différence de tension stabilisée 10 entre ses bornes lorsqu'aucun courant ne circule. On dit que la cellule est en état d'équilibre. La littérature, notamment la publication de Karthikeyan et al. Divulgué en 2008 et intitulé "Thermodynamic model development for lithium intercalation electrodes », numéro 15 185, page 1398-1407, a démontré que le vieillissement d'une cellule est dépendant du taux d'insertion en lithium des bornes d'une cellule. On notera que le taux d'insertion en lithium des bornes d'une cellule est également dépendant de l'état de charge 20 d'une cellule électrique 121. On peut donc affirmer que le vieillissement peut être modélisé à partir d'une modélisation d'une courbe de tension à vide de référence, correspondant à la courbe de tension à vide à l'état neuf, c'est à dire en sortie du constructeur, et d'un paramétrage de coefficient de 25 perte de capacité de la cellule et une définition spécifique de l'état de charge. On prend également l'hypothèse que les conditions de température durant l'exécution du procédé d'identification du vieillissement de la courbe de tension à vide n'impactent pas 30 le procédé. Dans la suite de la description, on prend l'hypothèse que la courbe de tension à vide de référence est modélisée par une interpolation polynomiale CO. A partir du polynôme CO, d'une fonction de vieillissement 50 paramétrée par un coefficient de perte de capacité de la cellule a, il est possible d'identifier la courbe de tension à vide de la cellule après le vieillissement suivant l'instant de référence. La courbe de tension à vide identifiée par le procédé est un polynôme Cl. La figure 6 représente une modélisation Cl de la courbe de tension à vide d'une cellule après vieillissement et une modélisation CO du polynôme de la courbe de tension à vide de référence. On appelle a le coefficient de perte de la capacité de la cellule suivant le vieillissement de la cellule après un instant de référence, correspondant ici à l'instant pour lequel la courbe de tension à vide de référence a été déterminée. On prend également l'hypothèse que la variation de la tension à vide de l'anode de la cellule 121 est invariante au vieillissement sur les plages de valeur de fonctionnement d'une cellule électrique. Cette hypothèse est prouvée dans l'article publié par Schmidt et mentionné dans l'état de la technique décrivant la variation des tensions à vide en fonction des taux d'insertion en lithium de chaque borne. La fonction de vieillissement 50 de la cellule 121 exploitée dans le cadre du procédé de l'invention par le dispositif de pilotage du stockeur 12 est fonction du polynôme CO de la courbe de tension à vide de référence et est alors décrite par la formule suivante, représentée en figure 5: aCSO C - SOCLOC.:!4, avec : - U représentant la courbe de tension à vide pour les différents états de charge de la cellule, - k, un indice de coefficient constituant le polynôme CO parmi N coefficients du polynôme CO, - Dkt un coefficient du polynôme CO et correspondant à un indice k compris entre 0 et N, - SOC un état de charge de la cellule 121, - SOCl00%, l'état de charge maximum de la cellule 121, - a, le coefficient de perte de capacité de la cellule. La figure 2 représente le procédé d'identification de la courbe de tension à vide de la cellule 121. Il comprend les étapes successives ci-après pouvant être réalisées selon deux modes de calcul distincts du coefficient de perte de capacité a de la cellule. Tout d'abord, le procédé comprend une première étape de détermination 21 du polynôme CO représentant une courbe de tension à vide de référence de la cellule 121. Le polynôme CO peut être obtenu selon deux méthodes. Une première méthode est la mesure à l'instant de référence, c'est à dire pour un état de sortie d'usine du constructeur du stockeur 12, de la tension à vide pour chacune des cellules 121 du stockeur 12 pour l'ensemble des états de charge de la cellule 121. Une deuxième méthode est l'utilisation des mesures de référence fournies par le constructeur du stockeur 12. A la suite des deux méthodes d'obtention du polynôme CO modélisant la courbe de tension à vide d'une cellule, les facteurs composants le polynôme CO sont mémorisés dans une mémoire accessible par le dispositif de pilotage 13 du stockeur. Le polynôme CO de chacune des cellules est ainsi mémorisé et exploitable par le dispositif de pilotage. Lors du fonctionnement du véhicule, le procédé comprend ensuite une deuxième étape de pilotage 22 de la cellule pour maintenir au moins des premier et deuxième paliers de tension à vide Pl, P2 à au moins des premier et deuxième états de charge SOC11, SOC21 respectivement. Les figures 3A et 3B représentent le pilotage de la cellule lors de la mise en oeuvre de la deuxième étape de pilotage 22. Plus précisément, les figures 3A et 3B représentent le pilotage de quatre paliers de tension Pl, P2, P3, P4 pour quatre états de charge SOC11, SOC21, SOC31 et SOC41 distincts de la cellule 121 lors d'une phase de recharge de la cellule 121. On notera qu'un palier de tension est défini par une phase d'équilibre de la cellule lorsque celle-ci est relaxée. La cellule est en équilibre lorsque le courant qui la traverse est très faible, sensiblement nul, et lorsque la tension ne varie pas. Dans le cadre du véhicule automobile, le dispositif de pilotage 13 du stockeur électrique 12 ne sollicite ni en charge ni en décharge le stockeur 12 et donc les cellules 121. Une sollicitation est opérée lorsqu'un chargeur au réseau électrique est en fonctionnement (s'il s'agit d'un véhicule disposant d'une prise de branchement) ou lorsqu'un équipement 14 ou 15 (machine électrique ou réseau de bord) consomme de l'énergie électrique ou délivre de l'énergie électrique. Ainsi, la deuxième étape 22 du procédé peut être commandée dans une première variante lors d'une phase de recharge du stockeur 13 sur le réseau électrique (pour le cas d'un véhicule disposant d'une prise de branchement), ou dans une deuxième variante lors des phases de vie entre deux arrêts de la machine électrique 14 (pour le cas d'un véhicule hybride électrique). Une troisième variante commande durant le roulage du véhicule les paliers de tension Pl, P2, P3, P4 à des instants lorsqu'aucune sollicitation n'est exercée sur le stockeur 12. Par ailleurs, le pilotage d'une séquence de paliers est commandé lorsque le dispositif de pilotage 13 reçoit ou émet une demande d'identification de la courbe de tension à vide d'une cellule, par exemple sur requête extérieure au véhicule ou sur demande d'un programme de diagnostic du stockeur 12 interne au véhicule. D'autres séquences de pilotages comportant seulement deux paliers sont possibles, préférentiellement à l'état de charge minimum de la cellule et à l'état de charge maximum de la cellule. Cependant, l'identification de la courbe de tension à vide de la cellule en vieillissement sera préférée avec une séquence d'une pluralité de paliers successifs entre l'état de charge minimum et l'état de charge maximum de la cellule 121. Ainsi, selon une variante du procédé comportant de multiples paliers de tension à vide, la deuxième étape de pilotage de la cellule comporte en outre le maintien à une pluralité de troisièmes paliers de tension à vide intermédiaires entre le palier correspondant à un état de charge minimum P2 et le palier correspondant à un état de charge maximum Pl. Ensuite, le procédé d'identification comprend une troisième étape de mesure 23 des première et deuxième tensions vide Ull, U21 pour les premier et deuxième paliers respectivement et d'une première variation de charge de la cellule entre les premier et deuxième paliers Ull et U21. Selon la variante du procédé comportant de multiples 5 paliers Pl, P2, P3 et P4, la troisième étape de mesure 23 comporte la mesure d'une pluralité de troisièmes tensions à vide U31 et U41 pour les troisièmes paliers P3, P4 respectivement et une pluralité de deuxièmes variations de charge de la cellule 121 entre chacun des troisièmes paliers 10 P3 et P4 et entre le deuxième palier et au moins un des troisièmes paliers. On notera que les variations de charge de la cellule 121 entre les paliers Pl, P2, P3 et P4 sont obtenues par la mesure de la variation de l'intensité entre les paliers de 15 tension et par le calcul de l'intégrale de la variation du courant. En outre, le procédé d'identification comporte une quatrième étape de détermination 24 d'un coefficient de perte a de la capacité de la cellule 121 à partir d'au moins un 20 premier facteur de perte al calculé selon la fonction de vieillissement 50 du polynôme CO spécifiquement appliquée avec les première et deuxième tensions à vide Ull, U21 et la première variation de charge. Cette méthode de calcul est justifiée par les hypothèses suivantes : 25 - le polynôme CO est strictement croissant, l'application de la fonction de vieillissement 50 pour une tension à vide mesurée permet d'obtenir la relation suivante entre les états de charge correspondant SOCtl= CXt*SOCto+SOCrrta, (1- Ot) r avec SOCtl correspondant à l'état de 30 charge pour une mesure de tension à vide mesurée Utl, SOCto l'état de charge pour le polynôme CO pour une tension à vide Uto égale à Utl, SOCmax l'état de charge de la cellule maximum et 0t le coefficient de perte pour la mesure de tension à vide mesurée, - la tension à vide est identique quel que soit le 5 vieillissement de la cellule pour l'état de charge maximum de la cellule dans la plage de fonctionnement. En outre, le procédé propose deux méthodes de calcul des facteurs de perte al, out pouvant être réalisées par le dispositif de pilotage 13. 10 Une première méthode, représentée sur la figure 2 et la figure 3A, pour le couple de palier Pl, P2 comprend : - le calcul 241 d'une première variation 31 d'état de charge de la cellule 121 entre les première et deuxième tensions à vide mesurées Ull, U21, la première tension à vide 15 mesurée Ull étant la tension à vide pour l'état de charge maximum SOC11 de la cellule, - le calcul 242 de l'état de charge SOC21 de la cellule 121 pour la deuxième tension à vide U21 en soustrayant la première variation d'état de charge 31 à l'état de charge 20 maximum SOC11 de la cellule, - le calcul 243 du premier facteur de perte al selon la tension à vide mesurée U21 et la fonction de vieillissement 50 appliquée pour l'état de charge SOC21 et pour la deuxième tension à vide U21 de la cellule. 25 Grâce à un couple de paliers de tension à vide (P1, P2 par exemple), on déduit la tension à vide mesurée U21 et son état de charge SOC21, la tension à vide pour le polynôme CO de la courbe de tension à vide de référence et l'état de charge correspondant SOC20. De ces valeurs on déduit un premier 30 facteur de perte de la capacité de la cellule. Une deuxième méthode de calcul d'un facteur de perte pour le couple de palier Pl, P2 est représenté par la figure 2, la figure 3B et la figure 4. La figure 3B représente les étapes du procédé d'identification pour la mesure des tensions 5 à vide Ull, U21, U31 et U41 des paliers Pl, P2, P3 et P4. Les états de charge SOC11, SOC21, SOC31 et SOC41 de chaque palier Pl, P2, P3 et P4 sont calculées en fonction de la variation de courant lors du pilotage des paliers. Cette méthode de calcul exploite le fait que la variation de l'état de charge est 10 différente à chaque état de vieillissement d'une cellule. Ainsi, pour la détermination 24, le calcul d'un facteur de perte al (pour le couple de palier Pl, P2) comprend: - le calcul 244 de la variation 34 d'état de charge de la cellule entre les première et deuxième tensions à vide 15 mesurée Ull, U21, et de la variation 37 d'état de charge de la cellule entre les tensions à vide correspondante U10, U20 du polynôme CO, (avec U20=U21 et U10=U11), - le calcul 245 du facteur de perte al en fonction du rapport entre la variation 34 et la variation 37. 20 Du fait de l'application de la fonction de vieillissement 50 et en considérant les états de charge d'une cellule en vieillissement et les états de charge du polynôme CO de la courbe de tension à vide de référence, deux états de charge de même tension à vide sont liés par la relation 25 suivante : SOC100."K SOCc°, étant un état de charge de la cellule pour le polynôme CO, SOC, étant un état de charge de la cellule en vieillissement calculé lors de l'étape 24, SOCl00%, étant l'état de charge maximum de la cellule. Ainsi, on peut affirmer pour le couple de palier de tension à vide P1 et P2 que la variation d'état de charge 34 est égale à la variation d'état de charge 37 divisée par le facteur de perte al. Et, connaissant la variation des états de charge du polynôme CO et les variations d'états de charges mesurés, on calcule le coefficient de perte al. De préférence, le procédé d'identification comporte une pluralité de paliers de tension Pl, P2, P3, P4, et la quatrième étape de détermination 24 du coefficient de perte a comporte en outre une sélection parmi le premier facteur de perte al et parmi une pluralité de deuxièmes facteurs de perte at, chaque deuxième facteur de perte at étant calculé identiquement au premier facteur de perte al et en fonction des troisièmes tensions à vide U31, U41 et des deuxièmes variations de charge mesurées durant la troisième étape de mesure 23. La sélection peut être réalisée selon plusieurs variantes de méthode sélective. Selon une première méthode, représentée en figure 5, la sélection est une fonction de calcul de la moyenne 51. Une moyenne de l'ensemble des facteurs de perte obtenus 01_ r out est calculée selon la formule suivante : Avec a le coefficient de perte, M le nombre de paliers de tension piloté et out le facteur de perte calculé à partir d'un couple de mesures de palier, t étant l'indice du facteur de perte parmi les M-1 facteurs de perte calculés. Dans une variante privilégiée, la méthode sélective pour déterminer le coefficient de perte consiste à utiliser la 5 méthode des moindres carrées et à utiliser un critère J adapté à la méthode de calcul des facteurs de perte ai, Selon la première méthode de calcul représentée par la figure 3A, le critère à minimiser est la différence entre la tension à vide mesurée et la tension à vide calculée avec la 10 pluralité de facteurs de perte at et l'application de la fonction de vieillissement 50 pour un état de charge donné est le suivant: S C, Avec, Ut la mesure de la tension à vide mesurée pour un 15 état de charge de la cellule étudié, Dkneuf le coefficient du polynôme correspondant à l'état de charge de la cellule étudié, SOCt l'état de charge de la cellule étudié correspondant à la mesure de tension à vide Ut et a le facteur de perte propre à la courbe de tension à vide Cl. 20 Selon la deuxième méthode de calcul, représentée par la figure 3B et la figure 4, à chaque variation d'état de charge entre le palier de tension à vide minimum et un palier de tension à vide, le critère à minimiser est la différence entre les variations d'état de charge mesurées ASOC(Ui;Ut) et les 25 variations d'état de charge calculées avec plusieurs valeurs d'un facteur de perte ASOC(Ui;Ut;a) . Selon cette méthode de calcul, le critère à minimiser peut être représenté par la formule suivante : ASOC(U Ut 3SOC U1; - t-2 On notera que la méthode des moindres carrés est une méthode de sélection parmi d'autres et ne limite en aucun cas la portée de l'invention. Par ailleurs, une fois que le coefficient de perte de la capacité de la cellule est déterminé on peut calculer la courbe de tension à vide de la cellule en vieillissement. Le procédé comporte une cinquième étape de calcul 25 de la courbe de tension à vide Cl selon la fonction de vieillissement 50 paramétrée avec le coefficient de perte a déterminée lors de l'étape 24 et appliquée à l'ensemble des coefficients du polynôme CO de la courbe de tension à vide de référence. La figure 6 est un graphique ayant en abscisse les états de charge de la cellule et en ordonnée les tensions à vide. Le graphique représente le polynôme CO de la courbe de tension à vide de référence et la courbe de tension à vide Cl obtenu par le calcul pour chaque coefficient du polynôme CO en appliquant la fonction de vieillissement 50 paramétrée avec le coefficient de perte de capacité de la cellule. La justesse du procédé a été validée par des mesures sur une cellule d'accumulation d'énergie électrique après deux ans de vieillissement. Le procédé présente de nombreux avantages cités ci-après. Il suffit d'identifier une courbe de tension à vide de 25 référence avant l'utilisation de la cellule électrique (obtenue par les données du constructeur ou une procédure de mesure initiale). Le procédé est opérable dans tout type d'équipement disposant d'un dispositif de pilotage de la cellule. En effet, les étapes 22, 23, 24, 25 du procédé successives à l'étape de détermination 21 du polynôme CO sont exécutées lorsque la cellule 121 est installée dans un équipement hôte 11 et lorsque celui-ci fonctionne au moyen de l'énergie électrique fournie par la dite cellule 121. Dans le cas d'un stockeur, la courbe de tension à vide de chaque cellule peut être identifiée. En outre, l'invention prévoit un dispositif de pilotage d'un stockeur électrique comprenant un moyen de calcul opérant le procédé d'identification du vieillissement de la courbe de tension à vide pour chacune des dites cellules du stockeur, et un moyen de détermination des indicateurs de fonctionnement du stockeur électrique, par exemple une courbe de tension à vide du stockeur et un état de charge, a partir des courbes de tension à vide de chacune des dites cellules calculées par le dit procédé d'identification. Il n'est pas nécessaire de réaliser un processus de charge et décharge complet de la cellule car les mesures peuvent être réalisées à partir de deux paliers de tension à vide seulement. Cependant, la méthode privilégie le pilotage de plusieurs paliers afin d'obtenir une mesure plus précise. La capacité réelle estimée ne nécessite pas de charge ni décharge complète. Grâce à l'identification du vieillissement de la courbe de tension à vide, notamment dans un véhicule électrique ou hybride, l'invention permet une amélioration de la précision et de la fiabilité de l'indicateur de niveau d'énergie dans le stockeur et la précision de pilotage par le dispositif de pilotage et une amélioration de la durée de vie du stockeur. II permet aussi de réduire le dimensionnement du stockeur, et ainsi le coût, et d'améliorer le diagnostic du stockeur. REVENDICATIONS1. Procédé d'identification du vieillissement d'une courbe de tension à vide (C1) d'une cellule d'accumulation 5 d'énergie électrique (121), caractérisé en ce qu'il comprend les étapes successives suivantes : - la détermination (21) d'un polynôme (CO) représentant une courbe de tension à vide de référence de la cellule (121), le pilotage (22) de la cellule pour maintenir au 10 moins des premier et deuxième paliers de tension à vide (P1, P2) à au moins des premier et deuxième états de charge (SOC11, SOC21) respectivement, - la mesure (23) des première et deuxième tensions à vide (U11, U21) pour les premier et deuxième paliers 15 respectivement et d'une première variation de charge de la cellule entre les premier et deuxième paliers (U11, U21), - la détermination (24) d'un coefficient de perte a de la capacité de la cellule à partir d'au moins un premier facteur de perte al calculé selon une fonction de 20 vieillissement (50) du polynôme (CO) spécifiquement appliquée avec les première et deuxième tensions à vide (U11, U21) et la première variation de charge, - le calcul (25) de la courbe de tension à vide (Cl) selon la fonction de vieillissement (50) paramétrée avec le 25 coefficient de perte a et appliquée à l'ensemble du polynôme (CO). 2. Procédé selon la revendication 1, caractérisé en ce que la fonction de vieillissement (50) correspond à la formule suivante : SO C10013/, 5° ri o D,3{: avec : - U représentant la courbe de tension à vide (Cl), - k, un indice de coefficient constituant le polynôme (CO) parmi N coefficients du polynôme (CO), Dkt un coefficient du polynôme (CO) et correspondant à un indice k, - SOC un état de charge de la cellule (121), - SOC100%, l'état de charge maximum de la cellule (121), -a, le coefficient de perte de capacité de la cellule. 3. Procédé selon la revendication 1 ou 2, caractérisé en ce que : - le pilotage (22) comprend en outre le maintien à une pluralité de troisièmes paliers de tension à vide (P3, P4), - la mesure (23) comprend en outre la mesure d'une pluralité de troisièmes tensions à vide (U31, U41) pour les troisièmes paliers (P3, P4) et d'une pluralité de deuxièmes variations de charge de la cellule (121) entre chacun des troisièmes paliers (P3, P4) et entre le deuxième palier et au moins un des troisièmes paliers, - la détermination (24) du coefficient de perte a de la capacité de la cellule comprend une outre une sélection (51) parmi le premier facteur de perte al et une pluralité de deuxièmes facteurs de pertes at, chaque deuxième facteur de perte at étant calculé identiquement au premier facteur de perte al en fonction des troisièmes tensions à vide (U31, U41) et des deuxièmes variations de charge. 4. Procédé selon l'une quelconque des revendications 1 3, caractérisé en ce que, pour la détermination (24) ducoefficient de perte cx, le calcul du premier facteur de perte al comprend : - le calcul (241) d'une première variation (31) d'état de charge de la cellule entre les première et deuxième tensions à vide mesurées (U11, U21), la première tension à vide mesurée (U11) étant la tension à vide pour l'état de charge maximum (SOC11) de la cellule, - le calcul (242) de l'état de charge (SOC21) de la cellule pour la deuxième tension à vide (U21) en soustrayant 10 la première variation d'état de charge (31) à l'état de charge maximum (SOC11) de la cellule, - le calcul (243) du premier facteur de perte al selon la tension à vide mesurée (U21) et la fonction de vieillissement (50) appliquée pour l'état de charge (SOC21) et 15 pour la deuxième tension à vide (U21) de la cellule. 5. Procédé selon l'une quelconque des revendications 1 3, caractérisé en ce que, pour la détermination (24) du coefficient de perte cx, le calcul du premier facteur de perte al comprend: 20 - le calcul (244) d'une deuxième variation (34) d'état de charge de la cellule entre les première et deuxième tensions à vide mesurée (U11, U21), et d'une troisième variation (37) d'état de charge de la cellule entre les tensions à vide correspondante (U10, U20) du polynôme (CO), 25 - le calcul (245) du premier facteur de perte al en fonction du rapport entre la deuxième variation (34) et la troisième variation (37). 6. Procédé selon l'une quelconque des revendications 3 5, caractérisé en ce que le pilotage de la cellule pour 30 maintenir le premier (P1), deuxième (P2) et les troisièmes (P3, P4) paliers de tension à vide est opéré dans un même cycle de recharge de la cellule (121). 7. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que le polynôme (CO) correspond à la courbe de tension à vide pour un état neuf de la cellule (121). 8. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'au moins les étapes (22, 23, 24, 25) du procédé successives à l'étape de détermination (21) du polynôme (CO) sont exécutées lorsque la cellule (121) est installée dans un équipement hôte (11) et lorsque celui-ci fonctionne au moyen de l'énergie électrique fournie par la dite cellule (121). 9. Dispositif de pilotage (13) d'un stockeur électrique (12), le dit stockeur électrique comprenant une pluralité de cellules électriques (121), caractérisé en ce qu'il comprend : - un moyen de calcul opérant le procédé d'identification du vieillissement de la courbe de tension à vide conformément à l'une quelconque des revendications 1 à 8 pour chacune des dites cellules (121) du stockeur (12), - un moyen de détermination des indicateurs de fonctionnement du stockeur électrique (12), par exemple une courbe de tension à vide du stockeur et un état de charge, partir des courbes de tension à vide de chacune des dites cellules (121) calculées par le dit procédé d'identification. 10. Véhicule automobile comprenant un stockeur électrique, une machine de traction électrique et un dispositif de pilotage du stockeur électrique, caractérisé en ce que le dispositif de pilotage est conforme à la revendication 9.30