FR 2949565 A1 20110304 FR 0904159 A 20090902 Domaine de l'invention La présente invention concerne une méthode pour estimer les caractéristiques d'un système électrochimique, de type batterie, qui ne sont pas directement mesurables. La méthode permet de gérer une batterie électrochimique, notamment pendant son fonctionnement dans un véhicule hybride ou électrique. La méthode selon l'invention permet de simuler le comportement thermique et/ou électrique d'une batterie. La méthode peut également être utile hors ligne, notamment pour dimensionner une batterie et optimiser les stratégies de gestion énergétique et thermique en fonction de l'application visée. La batterie électrochimique est l'un des composants les plus critiques d'un véhicule hybride ou électrique. Le dimensionnement de cet organe au sein de l'architecture véhicule repose sur une simulation de la batterie qui vise à représenter précisément son comportement électrique et thermique à différents niveaux de sollicitation dynamique. Le fonctionnement correct du véhicule s'appuie sur un système intelligent de gestion de la batterie (BMS) qui s'occupe de faire fonctionner la batterie au meilleur compromis entre les différents niveaux de sollicitation dynamique. La connaissance précise et fiable de l'état de charge (SoC), de l'état de santé (SoH) et de l'état thermique (T) est indispensable pour le BMS. L'état de charge d'une batterie (SoC) est sa capacité disponible (exprimée comme pourcentage de sa capacité nominale). La connaissance du SoC permet d'estimer combien de temps la batterie pourra continuer à fournir de l'énergie à un courant donné ou jusqu'à quand elle pourra en absorber. Cette information conditionne le fonctionnement du véhicule entier et notamment la gestion de l'énergie entre ses composants. Pendant la vie d'une batterie, ses performances tendent à se dégrader graduellement à cause des variations physiques et chimiques qui ont lieu avec l'usage jusqu'à rendre la batterie inutilisable. L'état de santé (SoH), qui est la capacité disponible après une recharge (exprimée en Ah), est donc une mesure du point qu'a été atteint dans le cycle de vie de la batterie. L'état thermique (T) d'une batterie conditionne ses performances, car les réactions chimiques et phénomènes de transport mis en jeu dans les systèmes électrochimiques sont thermiquement activés. L'état thermique initial est lié à la température extérieure au véhicule qui peut être amené à fonctionner sur une large plage de températures, typiquement entre -40°C et +40°C. L'état thermique en service évolue en fonction de la sollicitation en charge et décharge de la batterie, de son design et de son environnement. L'estimation plus précise et fiable du SoC, du SoH et de l'état thermique T a comme conséquence plusieurs avantages. Cette estimation permet d'éviter que le superviseur du véhicule ne soit amené à se comporter d'une façon trop prudente dans l'utilisation du potentiel d'énergie de la batterie ou inversement. Elle permet aussi d'éviter des surdimensionnements de sûreté des batteries, donc de sauvegarder du poids à bord et, en conséquence, du carburant consommé; elle permet également de réduire le coût total du véhicule. Un estimateur correct constitue donc une garantie pour une exploitation efficace et sûre de la capacité de la batterie dans toute la plage de fonctionnement du véhicule. Objet de l'invention La présente invention concerne une méthode de simulation du comportement et d'estimation de l'état interne d'un système électrochimique rechargeable (de type batterie), à partir d'un modèle mathématique à paramètres concentrés (modèle zéro-dimensionnel) de la batterie qui peut de manière avantageuse être exécuté de façon synchrone au fonctionnement de la batterie elle-même (temps réel). En particulier, pour reconstruire l'état interne de la batterie à partir des mesures facilement obtenues par des moyens classiques, la méthode permettra notamment d'estimer l'état de charge (SoC), l'état de santé (SoH) et l'état thermique (T) d'une batterie électrochimique, qui sont les caractéristiques internes les plus intéressantes pour des applications concernant les véhicules hybrides et électriques, mais qui ne sont pas des valeurs directement mesurables. La méthode selon l'invention peut être utilisée pour: - le contrôle embarqué et la gestion d'énergie d'un véhicule hybride ou électrique, - un simulateur du comportement électrique et thermique d'une batterie, par exemple dans le cadre du dimensionnement d'une batterie pour une application donnée, d'une calibration, d'une optimisation, ou d'une validation des stratégies de gestion énergétique et thermique de la batterie, - un chargeur/déchargeur de batterie. L'invention concerne également un système intelligent de gestion énergétique et thermique de la batterie pendant son fonctionnement. Art antérieur La méthode d'estimation du SoC désignée par les termes "coulombcounting" ou "book-keeping" est connue dans le domaine technique, mais elle conduit à des erreurs d'estimation en négligeant des phénomènes, notamment 15 comme l'autodécharge. La mesure de la tension à vide comme indicateur du SoC est également connue. L'usage d'autres indicateurs, par exemple l'estimation d'une résistance interne, (documents US6191590B1 et EP1835297A1) est une méthode également connue. 20 Ces deux dernières méthodes sont caractérisées par le fait que le SoC est préalablement associé à une ou plusieurs quantités mesurables ou facilement estimables, à travers des cartographies statiques ou des dépendances fonctionnelles analytiques. Pourtant, ces dépendances sont en réalité beaucoup plus complexes que ce qui est pris en compte dans le BMS, ce qui 25 amène souvent à des erreurs dans l'estimation du SoC. Une méthode potentiellement plus prometteuse est basée sur la mesure par spectroscopie d'impédance (SIE) d'une grandeur physique paramétrée par le SoC. Par exemple, la demande de brevet US-2007/0090843 propose de déterminer par SIE la fréquence f� associée à la transition 30 capacitive/inductive. Une corrélation entre la fréquence f� et le SoC est présentée pour une batterie au plomb, ainsi que pour les batteries Ni-Cd et Ni- MH. Une approche similaire à la précédente est basée sur la modélisation des spectres SIE par circuits électriques équivalents, dont les composants sont paramétrés par l'état de charge (SoC), comme décrit dans le brevet US6778913B2 de la société Cadex Electronics qui permettent le développement d'un testeur de batteries automobiles Spectro CA-12 basé sur la spectroscopie d'impédance électrochimique à multifréquences pour le couple acide-plomb. Les spectres SIE sont approchés par des circuits électriques équivalents, et l'évolution des composants est paramétrée par l'état de charge (SoC). De même, le brevet US6037777 de K. S. Champlin détermine l'état de charge et autres propriétés de batteries par la mesure des parties réelles et imaginaires de l'impédance ou admittance complexes pour des batteries au plomb ou d'autres systèmes. Une approche alternative se base sur des modèles mathématiques de batterie, afin d'utiliser des techniques d'estimation connues dans d'autres domaines. D'une manière générale, il existe deux grandes catégories de modèle du comportement électrique et thermique des batteries permettant de simuler les systèmes électrochimiques : d'une part les modèles électrochimiques, inspirés des travaux de Newman, reposant sur la connaissance des réactions chimiques et phénomènes physico-chimiques qui ont lieu à l'échelle microscopique dans la cellule de batterie ; d'autre part les modèles à circuit électrique équivalent, constitués d'éléments électriques tels que des résistances, capacités, inductances, agencés en série et/ou parallèle, visant à représenter au mieux la dynamique du comportement d'une batterie. La modélisation du comportement de batterie par analogie électrique est la plus courante, car il s'agit de modèles intuitifs, à paramètres concentrés (i.e. dépendants seulement du temps), donc peu gourmands en temps de calcul. Les éléments électriques, le plus souvent identifiés par une mesure physique, sont paramétrés par l'état de charge SoC, l'état de santé SoH, la température T, et la valeur du courant (données sous forme de cartographies). La demande de brevet US2007/0035307A1 décrit notamment une méthode pour estimer les variables d'état et les paramètres d'une batterie à partir des données de service (tension U, courant I, T), utilisant un modèle mathématique de batteries. Le modèle mathématique est caractérisé par le fait qu'il comprend une pluralité de sous-modèles mathématiques, et permet une réponse plus rapide. Les sous-modèles sont des modèles de type circuit électrique équivalent, dit modèles RC, associés à des gammes de fréquences restreintes. L'utilisation de modèles RC est décrite également dans le brevet EP880710 (Philips), la description des phénomènes électrochimiques et physiques aux électrodes et dans l'électrolyte servant de support au développement du modèle RC, la température de la batterie étant simulée par le modèle, afin de gagner en précision, par rapport à une mesure externe. Dans les modèles du type RC, l'état de charge (SoC) est toujours introduit uniquement pour paramétrer d'autres variables. Il n'est aucunement fait mention du SoC lui-même comme variable électrochimique. La précision de la prédiction des modèles par analogie électrique dépend des cartographies établies expérimentalement, avec un nombre d'essais limités. Aussi, ces modèles rencontrent souvent des difficultés à rendre compte précisément de la dynamique de batterie en condition d'usage nominal et extrême, sous courant fort ou nul, alors que les phénomènes de relaxation dans la batterie s'échelonnent durant plusieurs heures avant le retour à l'équilibre thermodynamique. Dans les modèles électrochimiques, les équations de la cinétique des principales réactions électrochimiques complétées par le bilan des masses et des charges à l'échelle d'un élément constituent un système d'équations algébro-différentielles, dont la résolution donne en chaque instant la concentration des espèces actives. L'état de charge (SoC) est défini à partir des concentrations des espèces actives aux électrodes. Une autre méthode d'estimation du SoC connue dans la littérature ([Gu, White, etc.]) repose sur la description mathématique des réactions d'un système électrochimique. Le SoC est calculé à partir de variables d'état du système. Cette description s'appuie sur des bilans de matière, charge, énergie, ..., ainsi que sur des corrélations semi empiriques. Dans la demande de brevet déposée le 28 Mars 2008 sous le numéro FR-08/01.709, est proposée une "METHODE POUR ESTIMER LES CARACTERISTIQUES NON MESURABLES D'UN SYSTEME ELECTROCHIMIQUE" qui repose sur la description mathématique des réactions d'un système électrochimique permettant d'obtenir un modèle de référence, puis d'en déduire un modèle réduit du comportement de la batterie par dérivation du modèle de référence. Mais les phénomènes de redistribution des charges et relaxation dans la batterie, au delà de quelques dizaines de secondes, sont toujours négligés. De manière surprenante, il a été trouvé qu'une méthode comprenant une description mathématique des réactions d'un système électrochimique directement sous forme d'un modèle OD (c'est-à-dire reposant sur une approximation �zéro-dimensionnelle�, détaillée ci-après), permet de rendre compte du comportement dynamique d'une batterie sous courant fort et sous courant nul, avec une précision accrue, de façon synchrone au fonctionnement de la batterie, et permet également d'estimer plus finement les caractéristiques internes non mesurables d'un système électrochimique. Description de l'invention Ainsi, l'invention concerne une méthode d'estimation de l'état interne d'un système électrochimique rechargeable dans lequel : o on dispose d'au moins un signal d'entrée d'au moins un paramètre représentatif d'une quantité physique du système électrochimique rechargeable, o on établit un modèle électrochimique à paramètres concentrés (OD ou � zéro-dimensionnel �) du système comprenant au moins : ^ une représentation mathématique de la cinétique des réactions électrochimiques ayant lieu aux interfaces entre chacune des électrodes et l'électrolyte, ladite représentation prenant en compte les concentrations d'interface, • une représentation mathématique de l'accumulation spatiale de capacité de double couche des charges à chaque électrode, ^ une représentation mathématique de la redistribution des charges à chacune des électrodes, ^ une représentation mathématique de la diffusion des charges ioniques de l'électrolyte au travers des différents compartiments poreux, ^ un bilan de matière dans toutes les phases du système, ^ un bilan énergétique du système pour calculer la température de la batterie, les concentrations des espèces et les autres variables étant homogènes au sein de chacune des régions du système électrochimique, o on génère au moins un signal de sortie dérivé du calcul. On peut distinguer la concentration interfaciale et la concentration moyenne C=f(t) pour chaque espèce active au sein de chaque région du système électrochimique. On peut décrire le potentiel d'équilibre thermodynamique de chaque électrode par une relation mathématique thermodynamique, par exemple les équations de Nernst, Margules, Van Laar, Redlich-Kister, ou analytique, par exemple: polynomiale, exponentielle, ou autre. Le signal d'entrée d'au moins un paramètre représentatif d'une quantité physique du système électrochimique rechargeable peut être l'intensité I, et/ou le potentiel V, et/ou la température T, et/ou la puissance électrique demandée. On peut recueillir comme signal de sortie, au moins l'un des signaux suivants: le potentiel, l'état de charge, l'état de santé, la température du 20 système électrochimique. On peut recueillir comme signal de sortie: la tension aux bornes du système électrochimique et la température du système électrochimique. On peut recueillir comme signal de sortie: l'état de charge, l'état de santé, la température du système électrochimique. 25 L'invention concerne également un système de gestion intelligent d'une batterie électrochimique comprenant: o un moyen d'entrée connecté à un moyen de mesure sur la batterie pour recevoir une valeur d'entrée d'au moins un paramètre représentatif d'une quantité physique de la batterie, 30 o un moyen de traitement pour générer au moins un signal de sortie d'au moins une caractéristique calculée par la méthode selon l'invention décrite ci-dessus, o un moyen d'information/commande pour présenter des informations sur la quantité physique de la batterie et/ou commander la charge/décharge et/ou le refroidissement de la batterie en réponse au signal de sortie du moyen de traitement et/ou de comparaison. Système de gestion selon l'invention dans lequel le moyen de traitement peut comprendre un filtre récursif. Le système de gestion peut être utilisé pour le contrôle et la gestion d'énergie d'un véhicule hybride ou électrique, ou d'une installation de fourniture d'énergie photovoltaïque ou éolien. Le système de gestion peut être utilisé pour le contrôle et la gestion d'un chargeur/déchargeur d'un véhicule hybride ou électrique, ou d'une installation de fourniture d'énergie photovoltaïque ou éolien. Également, l'invention concerne un simulateur de l'état thermique et/ou électrique d'une batterie comprenant : o un moyen d'entrée pour recevoir une valeur d'entrée d'au moins un paramètre représentatif d'une quantité physique d'une batterie, o un moyen de traitement pour générer au moins une caractéristique de sortie calculée par la méthode décrite ci-dessus. On peut utiliser la méthode selon l'invention pour un simulateur de 20 spectroscopie d'impédance. On peut utiliser ladite méthode pour simuler le comportement thermique et électrique d'une batterie de traction de véhicule hybride ou électrique, ou d'une installation de fourniture d'énergie photovoltaïque ou éolien. On peut utiliser ladite méthode pour le dimensionnement d'une batterie. Description détaillée de l'invention Le modèle électrochimique OD utilisé dans la méthode selon l'invention, dit modèle à paramètres concentrés, repose sur l'hypothèse que les 30 concentrations des espèces et les autres variables sont homogènes dans chacune des régions du système électrochimique correspondant typiquement 25 aux électrodes, au séparateur, et au compartiment destiné à recueillir les espèces gazeuses. C'est l'approximation homogène zéro-dimensionnelle (OD). Dans l'approximation homogène zéro-dimensionnelle, les variables dépendent donc uniquement du temps. Par ailleurs une approche pseudo-1D est utilisée à l'intérieur des grains de matière active constitutifs des matériaux d'électrode pour prendre en compte les aspects de diffusion en phase solide. Le modèle électrochimique OD de la méthode selon l'invention peut calculer les variations dans le temps de toutes les variables électrochimiques internes de 10 la batterie, et en particulier de l'état de charge. Comme une des entrées du modèle est le courant aux bornes de la batterie, les cas simulés dépendent du choix de cette dernière variable. Les grandeurs utilisables comme signal d'entrée du modèle sont, dans le cas 15 d'une batterie électrochimique : l'intensité I, la température T, le potentiel V, ou la puissance électrique demandée au système de stockage. De manière avantageuse, on distingue la concentration interfaciale et la concentration moyenne C=f(t) pour chaque espèce active dans chaque région du système électrochimique. 20 De manière avantageuse, on décrit le potentiel d'équilibre thermodynamique de chaque électrode par une relation mathématique thermodynamique (Nernst, Margules, Van Laar, Redlich-Kister) ou analytique (polynomiale, exponentielle...) On peut recueillir comme signal de sortie, le potentiel, et/ou l'état de 25 charge, et/ou l'état de santé, et/ou la température du système électrochimique. Avantageusement pour une application de la méthode à un simulateur de batterie, on recueille comme signaux de sortie: la tension aux bornes du système électrochimique et la température du système électrochimique. 30 Avantageusement pour une application de la méthode à un estimateur de l'état d'une batterie, on recueille comme signaux de sortie: l'état de charge, l'état de santé et la température du système électrochimique. L'invention concerne également un système de gestion intelligent d'une batterie électrochimique (notamment dénommé: "BMS" (Battery Management System)) comprenant: o un moyen d'entrée connecté à un moyen de mesure sur la batterie pour recevoir une valeur d'entrée d'au moins un paramètre représentatif d'une quantité physique de la batterie; o un moyen de traitement pour générer au moins un signal de sortie d'au moins une caractéristique calculée par la méthode utilisant le modèle électrochimique OD selon l'invention; o un moyen d'information/commande pour présenter des informations sur la quantité physique de la batterie et/ou commander la charge/décharge et/ou le refroidissement de la batterie en réponse au signal de sortie du moyen de traitement et/ou de comparaison. Le moyen de traitement peut comprendre un filtre récursif (par exemple de type Kalman). Ledit système de gestion est utilisable pour le contrôle embarqué et la gestion d'énergie d'un véhicule hybride. L'invention comprenant ledit système de gestion concerne également un chargeur/déchargeur de batterie. L'invention concerne en outre un simulateur de batterie comprenant : o un moyen d'entrée pour recevoir une valeur d'entrée d'au moins un paramètre représentatif d'une quantité physique d'une batterie; o un moyen de traitement pour générer au moins une caractéristique de sortie calculée par la méthode selon l'invention. Le simulateur de batterie permet notamment de simuler le comportement thermique et électrique de la batterie. L'invention concerne également un simulateur de spectroscopie d'impédance utilisant la méthode selon l'invention. La méthode selon l'invention permet la mise en oeuvre d'un procédé de 30 dimensionnement et/ou de conception d'une batterie. L'invention concerne également un simulateur du système véhicule hybride ou électrique comportant une batterie de traction, utilisant la méthode selon l'invention d'estimation des caractéristiques internes de la batterie. Description des figures: Les figures 1 à 11 illustrent un mode de réalisation de l'invention à titre non limitatif. La figure 1 illustre une représentation schématique du modèle selon la méthode, où les références ont la signification suivante : EBV: Équations de Butler-Volmer (équation 2) qui décrivent la cinétique des réactions électrochimiques; (1): Iapp est le courant appliqué à la cellule; (2): JZ est la densité de courant pour la réaction chimique z; (3): Jz+k est la densité de courant pour la réaction chimique z+k; (4): Jz+k+1 est la densité de courant pour la réaction chimique z+k+1; (5): Jn est la densité de courant pour la réaction chimique n; BCh: bilan de charge -représente les équations du bilan de charge électrique; (6): (bpos est le surpotentiel à l'interface du matériau d'électrode positive; (7): tneg est le surpotentiel à l'interface du matériau d'électrode négative; BMa: bilan de matière -représente les équations du bilan de matière et de conservation de la masse; (8): [X]°pos est la concentration de surface de l'espèce X à l'électrode positive; (9): [X]°neg est la concentration de surface de l'espèce X à l'électrode négative; (10): [X] z,pos,max est la concentration maximale de l'espèce X dans l'électrode positive; (11): [X]z,neg,max, est la concentration maximale de l'espèce X dans l'électrode négative; (12): V; est le potentiel aux bornes de la cellule électrochimique; (13): SOC;est l'état de charge de la batterie (State of Charge); Ech: état de charge -erprésente l'ensemble des équations donnant l'état de charge du système; Esanté: état de santé -représente l'ensemble des équations donnant l'état de santé du système; BEn: bilan d'énergie -représente l'ensemble des équations donnant l'évolution de la température en fonctions des flux thermique Phi générés lors du fonctionnement de la batterie; (14): T est la température du système; (15): SOH est l'état de santé de la batterie (State of Health). Le courant aux bornes de la cellule est considéré comme une entrée du modèle, tandis que la tension est une de ses sorties. Les signaux d'entrée, courant et température, sont représentatifs de quantités physiques mesurées sur la batterie. Des moyens de traitement basés sur les équations de Butler Volmer, le bilan des charges, le bilan de matière et le bilan énergétique calculent l'état de la batterie sur la base des signaux d'entrée et génèrent des signaux de sortie dérivés du calcul, comme le potentiel, l'état de charge, l'état de santé et la température. La figure 2 est une représentation schématique d'une cellule de batterie Ni-MH, où MH-el désigne l'électrode négative poreuse à base d'hydrure métallique, Ni-el l'électrode positive poreuse à base de Nickel, ReG le compartiment réserve de gaz, Sep le séparateur isolant électriquement les deux électrodes, Col les deux collecteurs de courant, et x la direction prévalente. Pour garantir la conduction ionique entre les deux électrodes lorsque il y a un flux de courant, les électrodes et le séparateur sont imprégnés d'une solution alcaline concentrée. Le gaz (oxygène ou Hydrogène) qui peut se libérer pendant la charge de la batterie se concentre dans un espace commun ReG au-dessus des cellules. La figure 3 est une représentation schématique d'une cellule de batterie Li-ion, où Neg désigne l'électrode négative poreuse à base de composés carbonés, LiMO2 l'électrode positive poreuse à base d'oxydes métalliques, Sep le séparateur isolant électriquement les deux électrodes, Col les collecteurs de courant, et x la direction prévalente. Pour garantir la conduction ionique entre les deux électrodes lorsque il y a un flux de courant, les électrodes et le séparateur sont imprégnés d'un électrolyte organique concentré en sels de lithium. La figure 4 représente un schéma du filtre de Kalman qui est appliqué à une cellule électrochimique selon la méthode de l'invention, avec X: état interne calculé par l'estimateur, U: entrée , Y: sortie, F: variation de l'état interne selon le modèle. La figure 5 représente un schéma de fonctionnement de l'algorithme d'estimation du SoC, avec Sph: système physique, M: modèle, FNL: filtre non linéaire, Est: estimateur, U: entrées mesurées, Y: sorties mesurées, Ye: sorties calculées par le modèle, Xe: état interne calculé par l'estimateur, F: variation de l'état interne selon le modèle, L: gain en sortie du filtre non linéaire. Les figures 6a,b,c,d montrent un exemple de prédiction en tension du modèle réduit selon l'invention à partir d'une batterie Ni-MH 6,5 Ah fabriquée par la société Panasonic, sollicitée à température ambiante suivant un régime de courant dynamique correspondant à un profil routier. Les résultats simulés sont comparés aux mesures expérimentales et aux simulations suivant l'art antérieur. La figure 6a donne l'évolution de la vitesse (en ordonnée) en fonction du temps T (en abscisse) sur un profil de roulage donné. La figure 6b donne la variation du courant batterie en fonction du profil de roulage. La figure 6c montre la comparaison entre les tensions simulées à partir d'un modèle V1 de l'art antérieur et du modèle V2 selon l'invention, et les mesures expérimentales. La figure 6d est un agrandissement de la figure 6c, entre l'intervalle de temps 850 et 1250. Les figures 7a,b,c,d montrent un exemple de prédiction en tension du modèle réduit selon l'invention d'une batterie Li-ion 2,3Ah de Al23s, sollicitée à différentes températures (0°, 10°, 20°, 30° respectivement pour les figures 7a, 7b, 7c et 7d) suivant un régime de courant dynamique selon un profil HPPC. Les résultats simulés par le modèle V2 selon l'invention, sont comparés aux mesures expérimentales. Les figures 8a,b,c,d montrent un exemple de prédiction en température du modèle réduit selon l'invention d'une batterie Ni-MH 6,5Ah de Panasonic, sollicitée à température ambiante suivant un régime de courant dynamique selon un profil HPPC. Les résultats simulés par V2 sont comparés aux mesures expérimentales sur les figures 8a, 8b et 8c. La figure 8a montre l'évolution de la tension de la batterie en fonction du temps entre 0% et 110% de SoC. La figure 8b est un agrandissement d'une phase de la figure 8a, sur une phase de pics de puissance en décharge/charge à 50% de SoC suivie d'une relaxation (à courant nul). Les pics de charge et décharge correspondent à des puissances d'environ +5,5 kW et -7,5 kW. La figure 8c montre l'évolution de la température de la batterie en fonction du temps. la figure 8d donne l'évolution de l'état de charge sur le profil HPPC HC. Les figures 9a,b,c montrent un exemple de prédiction en température du modèle V2 selon l'invention d'une batterie Li-ion 2,3Ah de Al23s, sollicitée à différentes températures suivant un régime de courant dynamique selon un profil HPPC. Les résultats simulés par V2 sont comparés aux mesures expérimentales. Les figures 9a, 9b et 9c correspondent respectivement à des températures de -20°, 0°, et 30°C. La figure 10 est un schéma fonctionnel d'un simulateur de véhicule hybride utilisant la méthode d'estimation des caractéristiques internes selon l'invention, dans lequel DP représente le profil de roulage, CS la stratégie de contrôle, V le véhicule, GB la boite de vitesse, M la motorisation, B la batterie, CU l'unité de contrôle, G le générateur, IC le système d'injection, FC la consommation de carburant, et D l'affichage. La figure 11 montre un exemple d'une courbe de spectroscopie d'impédance électrochimique simulée à partir de la méthode selon l'invention, représentant la partie imaginaire de l'impédance Imag(Z) (en ordonnée) en fonction de la partie réelle de l'impédance Real(Z) (en abscisse). Modèle électrochimique OD: Comme décrit précédemment, le modèle électrochimique OD, dit à paramètres concentrés, repose sur l'hypothèse que les concentrations des espèces et les autres variables sont homogènes dans chacune des régions du système électrochimique (par exemple de la cellule de batterie) correspondant typiquement aux électrodes, au séparateur, et au compartiment destiné à recueillir les espèces gazeuses. C'est ce qui est dénommé : l'approximation homogène zéro-dimensionnelle (OD). Les équations permettant la mise en oeuvre du modèle zéro-dimensionnel utilisé dans la méthode selon l'invention sont explicitées ci-après. Des réactions électrochimiques ont lieu aux interfaces entre les électrodes et l'électrolyte. L'électrode positive est le siège de réactions électrochimiques de réduction des espèces oxydantes, pendant la décharge, alors que l'électrode négative est le siège de réactions d'oxydation des espèces réductrices. L'équation généralisée des réactions électrochimiques entre des espèces réactives R et produites P est présentée dans la relation (1) : M' -4 y VrR -1- ER (1) z 1 4---car r' '1 La cinétique des réactions électrochimiques peut être décrite par les 20 équations de Butler-Volmer, dont la forme générale pour la réaction générique "z" est : Jins,z(t,T) = Jz,o n [P.]V 1=1 t ll[P f .%=1 N 0 11[Ri] exp(a0X,,nKî~r(t))ù 'N1 ~io exp(ùare,/ ,znK~lz(t)) exp [Ri lref i=1 Eaz 1 1 R T T 0 Y (2) où z = Y'z,electrode ù eq ,z ù Y'layer,z (3) dans laquelle JIns,Z est la densité de courant de transfert de charge (insertion), Jz,O est la densité de courant d'échange électronique, [X]° est la concentration de l'espèce X produite (P) ou consommée (R) à la surface d'une particule de matériau d'électrode, 1b(t) est la surtension d'électrode définie comme la différence de potentiel électrique à l'interface immédiate électrode/électrolyte entre (Dz,electrode le potentiel à l'électrode considérée, Ueq,z le potentiel d'équilibre thermodynamique et ~Z un terme de surtension supplémentaire due à la présence possible d'une couche d'accumulation des charges à la surface des grains de matière active. K est une fonction de la température T et TO est la température de référence (298K) (K = F/RT), et Ea,z est l'énergie d'activation de la réaction électrochimique considérée. La forme du potentiel d'équilibre thermodynamique Ueq,z peut prendre la forme d'expressions mathématiques complexes, décrites par les lois de Nernst, de Margules, de Van Laar ou bien de Redlich-Kister. Ces expressions sont paramétrées en fonction du taux d'avancement de la réaction électrochimique au sein de l'électrode et peuvent aussi prendre la forme de fonctions mathématiques polynomiales, exponentielles, hyperboliques. Les phénomènes d'hystérésis aux électrodes, le cas échéant, sont pris en compte par l'équation (4) : dU eq,z dt (Ueqz>°D ù Ueq,z )k(l app )1'app app k (I )= kc , I aPP > 0 kd ,Iap, 0 Ueq ,z,d ,Iapp 15 dans laquelle lllelectrode est une grandeur fonctionnelle caractéristique du matériau d'électrode et proportionnelle à une longueur de diffusion, Ds,, est le coefficient de diffusion de l'espèce X dans le matériau d'électrode. Bilan des charges 20 - Représentation mathématique de l'accumulation de double couche électrique des charges à chacune des électrodes dans le modèle électrochimique OD: En prenant en compte l'accumulation des charges à chaque électrode, la méthode selon l'invention étend le domaine de validité du modèle 25 électrochimique OD aux cas non-stationnaires, où des effets capacitifs aux interfaces solide-liquide sont importants comme par exemple "l'effet double couche". L'équation (7) retranscrit la prise en compte de la capacité de double couche dans le bilan des charges : d oz,electrode = rr l dt gelectrode \electrode Iapp ù Jz J OÙ O,electrode est une grandeur fonctionnelle caractéristique de l'électrode considérée où se déroule la réaction électrochimique z faisant intervenir l'espèce active X. Iapp est le courant appliqué au système de stockage et �electrode une grandeur fonctionnelle caractéristique de l'électrode considérée. JZ est la densité de courant totale traversant l'électrode considérée et est la somme du courant d'insertion Jlns,z et de courant de réactions secondaires Jsecondaire,z ou parasites Jpara,z de type processus de vieillissement. - Représentation mathématique de la redistribution des charges à chacune des électrodes dans le modèle électrochimique OD: Le phénomène de redistribution des charges dans la profondeur de 15 chaque électrode est pris en compte dans le modèle électrochimique OD général à l'aide des équations (8)-(10): (7) d0~ z, electrode dt gelectrode (Iapp Jz ) ù helectrode VRp,electrode (8) gelectrode \I app ' Jz / gelectrode \Iapp ù electrode `jz J d V Rp ,electrode _ d0{D z, electrode dt dt kelectrode VRp,electrode dans lesquelles gelectrode, est une fonction mathématique linéaire faisant intervenir le courant appliqué Iapp, ainsi que la densité de courant de transfert 20 de charge et où VRp,electrode est le sur potentiel de redistribution des charges dans la profondeur de l'électrode considérée. helectrode, gelectrode, Eelectrode, et kelectrode sont des grandeurs fonctionnelles caractéristiques du matériau d'électrode considérée. (9) (10) - Représentation mathématique de la diffusion ionique dans la porosité à chacun des compartiments du système dans le modèle électrochimique OD: Le modèle électrochimique OD générique retranscrit le phénomène de diffusion des charges de l'électrolyte dans la porosité des différents compartiments du système électrochimique en intégrant au bilan électrique global des termes de surtensions de concentrations qui prennent la forme de l'équation différentielle (11) : V i E N, d rici = /l1(T)(Iapp ù 0 (T )ici ) dt dans laquelle Ài et coi sont des fonctions mathématiques exponentielles caractéristiques du système de stockage et reflétant les aspects diffusifs à courts et longs termes. Ces fonctions dépendent de la température pour rendre compte du caractère classique thermiquement activé du processus de diffusion des espèces en présence. Généralement, deux termes de surtensions suffisent à rendre compte des phénomènes diffusifs. - Vieillissement des systèmes de stockage d'énergie: Le vieillissement des systèmes de stockage d'énergie se caractérise par une perte de capacité (perte d'énergie) accompagnée d'une augmentation de la résistance interne et des surtensions (perte de puissance), souvent causées par une ou plusieurs réactions parasites de vieillissement. La cinétique des réactions de vieillissement peut être décrite par une expression mathématique analogue à l'équation (2), donnée dans l'équation (12) : Tpaa,~(t,T) = o 'N~ PA exp (aox.znKr1z(t))- 1N [R;1 exp (ù ares ~nK,,,(t)) l [' l [Ri J, j=1 dans laquelle Jpara,z est la densité de courant de transfert de charge parasite, Jpara,z,o est la densité de courant d'échange électronique de la réaction parasite z. E"'( 1 1 \ exp - - R ,T Tä 1 (12) La perte de capacité se traduit par une diminution de la concentration maximale des espèces stockées dans le matériau d'électrode selon l'équation (13) : x electrode para ,z (13) dans laquelle Xelectrode est une grandeur fonctionnelle caractéristique du matériau d'électrode La perte de puissance se traduit par une augmentation de la résistance interne du système provoquée notamment par la croissance de couches résistives en surfaces des particules de matières actives selon l'équation (14) : d dt R layer � X layer 1 para , z (14) dans laquelle )(layer est une grandeur fonctionnelle caractéristique de la composition chimique des couches parasites formées en surface des électrodes par réactivité avec l'électrolyte. Le signe +/- dépend du mécanisme physico- chimique de développement de la couche de surface. - Bilan thermique: La température de la cellule peut être calculée comme sortie du bilan d'énergie. D'une part le flux de chaleur interne %en généré par l'activité de la 20 cellule électrochimique est donné par : (t) I (t) dUeq ~ . (t)) A(z~ ù v(t)l(t) (15) où le terme (Ueq,ref,z-V) peut être associé aux pertes irréversibles pour chaque réaction électrochimique z, tandis que le terme de génération réversible T dUeq,ref,z/dT est directement lié aux variations d'entropie dues aux v eq,r ef, réactions électrochimiques. D'autre part le flux transféré vers l'ambiant à température Ta, (Pt, est donné par la loi de Fourier : où h est un coefficient de transfert thermique associée aux phénomènes de convection et radiation, et ACeä est la surface de la cellule. Le flux thermique net à travers la batterie, cp, peut être facilement calculé comme différence entre les flux internes et externes, i.e. (p =Yen - (ptra• La quantité de chaleur stockée dans la batterie, obtenue par intégration du flux de chaleur dans le temps, permet de calculer la température de la batterie suivant les relations : (Ptra(t) = hAcell(T(t) ù Ta), (16) dT( n'lcellCp Ott) = lPgen (t) - 1Ptra(t) (17) où Cp est la capacité thermique spécifique de la cellule et MCeä sa masse. - Bilan électrique global: Le modèle électrochimique OD générique est complété par un bilan global 15 du potentiel électrique sur la cellule : N V(t,T) = electrode ,pos (t, T) ù electrode ,neg (t,T) + E 1l ci (t, T) + %S2 (t,T ) (18) i=1 dans laquelle V(t,T) est la tension aux bornes de la cellule, roc; sont des termes de surtension de concentration et tin est une surtension ohmique faisant intervenir la résistance interne du système, résultant des conductivités des phases solide et liquide. La résistance ohmique est une grandeur 20 fonctionnelle caractéristique du système qui prend la forme d'une fonction mathématique exponentielle. En résumé, le modèle électrochimique OD de la méthode selon l'invention illustré sur la figure 1 comprend les équations (1 à 18), soit un total a minima 25 en fonction des réactions électrochimiques prises en considération de 16 équations, pour les 16 variables J&, ..., Jn, il, ..., 11näLXl°pos [Xl°neg 'pos, CDneg,Vrpos, Vrneg, OIpos, Atneg,ni,nn V, T. Les autres grandeurs apparaissant dans les équations qui constituent la méthode sont traitées comme des paramètres à calibrer. - Définition de l'état de charge selon l'invention L'état de charge de la cellule dans la méthode selon l'invention, q(t), est donné par la concentration d'une des espèces réactives X selon la relation (19) [X Imax (19) dans laquelle y et b sont des grandeurs fonctionnelles caractéristiques des matériaux d'électrodes. Ce calcul se distingue nettement du calcul connu dans l'art antérieur dit 15 "coulomb counting", qui donne dq(t) I(t) dt Qmax La relation entre Xmax et Qmax est donnée par Q max = KF [X ]max (21) dans laquelle F est la constante de Faraday, K est une grandeur fonctionnelle caractéristique de la géométrie de l'électrode limitante. 20 L'estimation de q s'appuie donc sur l'estimation de X, alors que cette variable n'est pas directement mesurable à partir d'une batterie, en particulier à bord du véhicule. 10 SOC (t) = y + 8 CX l ` (20) Exemples d'applications aux technologies Ni-MH et Li-ion ) Cas d'une batterie Ni-MH Dans le cas d'une batterie Ni-MH, représentée schématiquement sur la figure 2, les espèces actives sont l'oxyhydroxyde de Nickel NiOOH, l'hydroxyde de nickel Ni(OH)2, l'hydrure métallique MH, l'oxygène 02 partiellement dissous 5 dans l'électrolyte en équilibre avec la phase gazeuse. Les réactions électrochimiques à l'électrode positive sont, pendant la décharge: NiOOH+H2O+e- -- Ni(OH)2+OH- (22) 0.502+H2O+e- -p 2OH- tandis qu'à l'électrode négative: (23) MH+OH- -+ H2O+e-+M 2OH- -- 0.502+H2O+2e- En appliquant l'équation (2) aux quatre réactions (22)-(23), avec 10 z=1,...4, on obtient: J1(t) = 1 (co(t)) Ce e0.5•K r t J ric)'_ Cn,ref ~'r.rrf C Cn,max Cn (t) 0.5~K rai (t) }: eX Cn,max ù.Cnn,ref / l ((t) l J3(t) J3,o P Ce e0.5•K 113(t) Cm,re f Ce,re f dans lesquelles en est la concentration des protons dans l'électrode positive (hydroxyde de nickel), ce est la concentration de l'électrolyte i.e. des 1/2 e1.5•K•r/a(t) Ca(t) eù0.5.K•rtz(t) CQTefJ 23 _e-0.5•K•113(4l . Cxp CTTe 1 2 , Cùo.5•K•Ca(t) (24) (25) (26) (27) ions OH co est la concentration de l'oxygène dans l'électrode négative et la même variable barrée est la concentration interfaciale d'oxygène, en équilibre avec la phase gazeuse, cm est la concentration de l'hydrogène dans l'électrode négative (matériau métallique), tandis que les indices "ref" et "max" se réfèrent aux valeurs de référence et maximales, respectivement ; enfin p représente l'ordre de réaction. Les surpotentiels riz sont définis comme ~11(t) = Opos(t)-Ueq,ref,1(t) 12(t) = (Dpos(t)-Ueq,ref,2(t) 13(t) _ (Dneg(t)-Ueq,ref,3(t) (28) 14(t) = (Dneg(t)-Ueq,ref,4(t) où (Dpos et (reneg sont les différences de potentiel solide-liquide à l'électrode positive et négative, respectivement. Des expressions pour Ueq,ref,l, Ueq,ref,2, Ueq,ref,3 et Ueq,ref,4 en fonction de la température sont connues. Par exemple, on connaît Ueq,re 1= K(T) (26= 1)+U+(T-Tp T (28a) où 0 est le SoC, kl une constante et U1 est caractérisé par sa dérivée par rapport à la température T. Pareillement, on connaît les expressions paramétrées : L eq,re dU~ dT (28b) pour z = 2,...,4. 2) Cas d'une batterie Li-ion Dans le cas d'une batterie Li-ion, les espèces actives sont des oxydes de métaux pour l'électrode positive et des Composés Carbonés, des Métaux, ou des oxydes de Métaux pour l'électrode négative. Une représentation schématique d'une cellule Li-ion est donnée sur la Fig. 3. Les réactions électrochimiques à l'électrode positive sont, pendant la charge, Lil_xMO2 + xe- + xLi+ ----> LiMO2 (29) tandis qu'à l'électrode négative, en prenant l'exemple d'un composé carboné, Li yC6 --> yLi + + 6C + ye (30) En appliquant l'équation (2) aux deux réactions (29)-(30), avec z=1,...2, 10 on obtient: F _ Cpos xexaoxpoS 77pos \ mspos=i0,in.~pos C posrefC C poSrnax pos x C X eX -eq'osL Ip°s ex L a,pos 1 1 (31) Cposmax--C po.~ref Ce,ref RT R \T T C C " " ( E /1 negnax eg x e x ex) edne eg ex, a,neg (32) ùnegnaxCneg'ef Cref RT /) R \ T 1 J) dans lesquelles cpos est la concentration des ions Li+ dans l'électrode positive (Oxyde métallique), ce est la concentration de l'électrolyte i.e. des ions Li+, et la même variable barrée est la concentration interfaciale en porteurs de charge, Cneg est la concentration des ions Li+ dans l'électrode 15 négative (matériau d'insertion généralement à base de graphite), tandis que les indices "ref" et "max" se réfèrent aux valeurs de référence et maximales, respectivement. Les surpotentiels ri, sont définis comme 77 pos = T pos ù v eq, pos ù RSLI, pos X j pos 1/neg = rneg ù v eq,neg ù RSLI,neg X J neg (33)10 Des expressions pour Ueq,pos, et Ueq,neg en fonction de la température sont connues: Ueq pos =2.8-1.33x+0.40exp(-254x)-0.14tanh(36.5(x-0.033))+(T-To) dUP°S (33a) dT Ueq,neg =0.12+2.5exp(-70x)-exp(50(x-1.04))+(T-To)dUneg (33b) dT Le tableau suivant récapitule les valeurs numériques nécessaires à la 5 calibration du modèle réduit suivant l'invention : - pour un élément prismatique de batterie Ni-MH de capacité 6,5Ah, commercialisé par Panasonic ; - pour un élément cylindrique de batterie Li-ion de chimie LFP/graphite de capacité 2,3Ah, commercialisé par Al23system. T=298K Ni-MH Li-ion Polarité et matériau Electrode (+) Electrode (-) Electrode (+) Electrode (-) d'électrode NiOOH MH LiFePO4 Carbone J01,ins 2.1 0 0.28 0 302,secondaire 8.0e-08 0 0 0 J03,ins 0 5.4 0 13.75 J04,secondaire 0 1.5e-011 0 0 kc 0.0035 0 0 0 kd 0.0025 0 0 0 I3electrode 1.0516 0.4834 1.03e+003 5.65 Welectrode 4.75e-010 0 3.73e-014 6.85e-012 D, le-014 le-015 2.92e-020 2.50e-013 A,electrode 0.0222 0.2059 0.74 0.028 Eelectrode 0.0903 -0.3098 0.0033 -0.39 helectrode, 1.5428e-004 0 0.0023 0.0023 eelectrode, 0.0020 -0.0638 0.0024 -0.0112 Selectrode, 0.0903 - 0.3098 307.09 -2.5620 kelectrod 1.7094e-004 0 0.0018 0.0021 À1 1.51e-005 1.44e-011 c,~1 176.06 1.81e+009 Rint 9.0e-04 9.0e-03 y 100 102.11 Ô -100 -126.58 Acell 0.0083 0.0063 Mcell 0.166 0.07 Cp 2000 1750 Présentation du filtre récursif La méthode utilise avantageusement un filtre récursif pour estimer l'état du système dynamique à partir des mesures disponibles, dont un schéma est proposé sur la figure 4. Des caractéristiques notoires de ce problème d'estimation sont le fait que les mesures sont affectées par du bruit, et le fait que le système modélisé selon la méthode est fortement non-linéaire. Un filtre récursif utilisé de préférence dans la méthode sera le filtre de Kalman étendu connu par l'Homme de l'art. Selon le modèle de la méthode, le vecteur d'état de la cellule électrochimique de batterie (figure 4) s'écrit x = {cn, cm, po, Atpos, A neg, T}, dont la première composante est liée a l'état de charge à estimer par l'équation (19). Les mesures disponibles sont la tension aux bornes de la cellule et la température de la batterie, qui représentent la sortie y du modèle, et le courant Iapp aux bornes qui représente l'entrée u du modèle. Selon la méthode connue du filtre récursif, les équations du modèle sont réorganisées en: x f (x, u y = h(x, u)I, Algorithme d'estimation du SoC (34) La méthode prévoit alors une étape (M dans la figure 5) où le modèle fournit le vecteur des variations f (F dans la figure 5) et la sortie calculée y (Ye dans la figure 5) selon l'équation (34). Ensuite, ces deux variables sont manipulées par une deuxième étape (Est dans la figure 5) qui reconstruit l'état Xe à partir de F, Ye, et de la mesure Y. L'algorithme d'estimation utilise ainsi la sortie d'une troisième étape (FNL dans la figure 5) qui fournit la variable L en fonction de l'état reconstruit, des caractéristiques du système électrochimique (selon le modèle de la méthode) et des caractéristiques du bruit qui affecte les mesures. L'étape FNL peut être réalisée avec une méthode connue par l'Homme de l'art, par exemple le filtre de Kalman étendu. Le modèle réduit représente directement l'état de charge comme une variable d'état du modèle. Au contraire, les méthodes connues utilisent des modèles dits "à circuit électrique équivalent", où l'état de charge n'est pas une variable dynamique du modèle, mais une variable exogène, en fonction de laquelle d'autres variables dynamiques ou statiques sont paramétrées. On comprend qu'un système de gestion d'une batterie électrochimique (BMS) utilisant la méthode suivant l'invention pour estimer les caractéristiques non mesurables de la batterie à partir des mesures couramment disponibles fournira en sortie une information plus fiable et plus précise. Estimation de l'état de santé selon l'invention Le modèle réduit selon l'invention est basé sur des paramètres physiques du système, et non sur des paramètres globaux équivalents comme les modèles RC connus dans l'art antérieur. Cette propriété facilite l'estimation du vieillissement et donc de l'état de santé de la batterie. En fait, les méthodes utilisées pour l'estimation de l'état de charge, qui se basent sur des observateurs d'état, peuvent être étendues pour comprendre une adaptation lente des paramètres du modèle. Cette extension est connue dans l'art antérieur pour plusieurs applications différentes. Dans le cas de l'invention, les mêmes signaux qui circulent dans la Figure 4 peuvent être utilisés aussi pour cette extension adaptative. Les variations estimées des paramètres du modèle réduit serviront à détecter d'éventuelles variations macroscopiques dans le comportement de la batterie, et donc des altérations de ses performances, ce qui est communément entendu comme "vieillissement". La reconnaissance et la quantification du vieillissement, à travers les variations relatives des paramètres, donne l'estimation recherchée de l'état de santé. Simulateur du comportement électrique et thermique d'une batterie Le modèle électrochimique OD de la méthode selon l'invention peut calculer les variations dans le temps de toutes les variables électrochimiques internes de la batterie, et en particulier de l'état de charge. Comme l'entrée du modèle est le courant aux bornes de la batterie, les cas simulés dépendent du choix de cette dernière variable. Par exemple, on peut représenter une charge ou une décharge contrôlée à courant constant, ou à courant variable selon un profil fixé, ou encore à courant variable en fonction de la tension. Ce dernier cas est représentatif des conditions de sollicitation de la batterie dans un véhicule, où le courant imposé à la batterie dépend de la tension, selon les caractéristiques des composants électriques associés (électronique de puissance, moteur(s) électrique(s), etc.). Des résultats typiques de prédiction du comportement électrique obtenus par un simulateur de batterie utilisant les modèles suivant l'invention sont présentés sur la figure 6 pour le cas de la batterie Ni-MH de l'exemple 1 et sur la figure 7 pour le cas de la batterie Li-ion de l'exemple 2. Dans les deux cas, la confrontation des résultats du modèle OD de la méthode selon l'invention aux résultats expérimentaux souligne la précision du rendu du comportement dynamique obtenu. Sur figure 6, la confrontation aux résultats donnés par un modèle réduit de l'art antérieur souligne les bénéfices de l'utilisation du modèle OD suivant l'invention. La présence du bilan d'énergie dans le modèle OD de la méthode selon l'invention permet de simuler l'évolution thermique du système, couplé avec l'évolution de l'état électrique donné par l'équation (18). Des résultats typiques de prédiction du comportement thermique d'une batterie à l'aide d'un simulateur utilisant les modèles suivant l'invention sont présentés sur la figure 8 pour le cas de la batterie Ni-MH de l'exemple 1 et sur la figure 9 pour le cas de la batterie Li-ion de l'exemple 2. Par conséquence, la méthode selon l'invention peut ainsi servir pour le dimensionnement de la batterie, la définition, la calibration et la validation des stratégies de management électrique et thermique, enfin l'optimisation des systèmes de gestion thermique, qui doivent nécessairement équiper la batterie elle-même. En fait, les flux thermiques générés et la température de la batterie sont des variables d'entrée pour ces systèmes, qui ont pour but le réglage de ces flux et de cette température autour des valeurs admissibles. La représentation des transitoires thermiques permet ainsi de synthétiser et valider les stratégies de contrôle et optimisation associés aux systèmes de gestion thermiques. Ces stratégies peuvent ainsi profiter de la présence d'un modèle réduit pendant leur utilisation en ligne, pour avoir des estimations de certaines variables qui ne sont pas mesurables (températures dans de points spécifiques, flux thermiques, etc.), ou qui sont mesurables, mais avec des temps de réponse des capteurs associés trop lents. Autres utilisations des modèles OD de la méthode Simulateur du système véhicule Le modèle OD selon l'invention est aussi utile comme aide au dimensionnement des chaînes de traction pour véhicules hybrides. Un exemple de simulateur du système véhicule hybride intégrant un modèle de batterie est donné figure 10. Typiquement, ces applications ont besoin de modèles de comportement d'une batterie à paramètres concentrés capables de simuler le comportement dynamique d'une batterie de traction plus efficacement et fidèlement que des modèles de type circuit électrique équivalent. Simulateur de spectroscopie d'impédance Le modèle OD de la méthode selon l'invention permet également de reproduire les essais expérimentaux de spectroscopie d'impédance, afin de prévoir les relations entre ces mesures et l'état de charge interne de la batterie. L'équation du potentiel est alors modifiée pour tenir compte des effets inductifs dus aux connections entre les cellules et avec les terminaux. V(t) _ A~pp$(t) ù A eg(t)+ Ri, I�t) + Ldtt) (35) Cette potentialité de la méthode selon l'invention est illustrée sur la 30 figure 11. Procédé de dimensionnement pour la production d'une batterie Tout procédé de production d'une batterie qui s'appuiera sur un simulateur du comportement électrique et thermique d'une batterie tirera avantageusement parti du modèle OD de la méthode selon l'invention, pour son temps de calcul minimisé, sa fiabilité et sa précision sur la prédiction des caractéristiques d'une batterie. Ainsi, on peut mettre en oeuvre un procédé de fabrication d'une batterie en dimensionnant la batterie avec la méthode selon l'invention. REVENDICATIONS1) Méthode d'estimation de l'état interne d'un système électrochimique rechargeable dans lequel : o on dispose d'au moins un signal d'entrée d'au moins un paramètre représentatif d'une quantité physique du système électrochimique rechargeable, o on établit un modèle électrochimique à paramètres concentrés (OD ou � zéro-dimensionnel �) du système comprenant au moins : ^ une représentation mathématique de la cinétique des réactions électrochimiques ayant lieu aux interfaces entre chacune des électrodes et l'électrolyte, ladite représentation prenant en compte les concentrations d'interface, une représentation mathématique de l'accumulation spatiale de capacité de double couche des charges à chaque électrode, une représentation mathématique de la redistribution des charges à chacune des électrodes, • une représentation mathématique de la_diffusion des charges ioniques de l'électrolyte au travers des différents compartiments poreux, ^ un bilan de matière dans toutes les phases du système, un bilan énergétique du système pour calculer la température de la batterie , les concentrations des espèces et les autres variables étant homogènes au sein de chacune des régions du système électrochimique , o on génère au moins un signal de sortie dérivé du calcul. 2) Méthode selon la revendication 1 dans laquelle on distingue la concentration interfaciale et la concentration moyenne C=f(t) pour chaque espèce active au sein de chaque région du système électrochimique. 30 3) Méthode selon l'une des revendications précédentes dans laquelle on décrit le potentiel d'équilibre thermodynamique de chaque électrode par une25relation mathématique thermodynamique (Nernst, Margules, Van Laar, Redlich-Kister) ou analytique (par exemple: polynomiale, exponentielle). 4) Méthode selon l'une des revendications précédentes dans laquelle le signal d'entrée d'au moins un paramètre représentatif d'une quantité physique du système électrochimique rechargeable est l'intensité I, et/ou le potentiel V, et/ou la température T, et/ou la puissance électrique demandée. 5) Méthode selon l'une des revendications précédentes dans laquelle on recueille comme signal de sortie le potentiel, et/ou l'état de charge, et/ou l'état de santé, et/ou la température du système électrochimique. 6) Méthode selon l'une des revendications précédentes dans laquelle on recueille comme signal de sortie la tension aux bornes du système électrochimique et la température du système électrochimique. 7) Méthode selon l'une des revendications 1 à 5 dans laquelle on recueille comme signal de sortie l'état de charge, l'état de santé et la température du système électrochimique. 8) Système de gestion intelligent d'une batterie électrochimique comprenant oun moyen d'entrée connecté à un moyen de mesure sur la batterie pour recevoir une valeur d'entrée d'au moins un paramètre représentatif d'une 25 quantité physique de la batterie; oun moyen de traitement pour générer au moins un signal de sortie d'au moins une caractéristique calculée par la méthode selon l'une des revendications 1 à 7. oun moyen d'information/commande pour présenter des informations sur 30 la quantité physique de la batterie et/ou commander la charge/décharge et/ou le refroidissement de la batterie en réponse au signal de sortie du moyen de traitement et/ou de comparaison.20 9) Système de gestion selon la revendication 8 dans lequel le moyen de traitement comprend un filtre récursif. 10) Utilisation du système de gestion selon l'une des revendications 8 ou 9 pour le contrôle et la gestion d'énergie d'un système électrochimique rechargeable d'un véhicule électrique ou hybride, ou d'une installation de fourniture d'énergie photovoltaïque ou éolien. 11) Utilisation du système de gestion selon l'une des revendications 8 ou 9 pour le contrôle et la gestion d'un chargeur/déchargeur. 12) Simulateur de l'état thermique et/ou électrique d'une batterie comprenant : o un moyen d'entrée pour recevoir une valeur d'entrée d'au moins un paramètre représentatif d'une quantité physique d'une batterie; o un moyen de traitement pour générer au moins une caractéristique de sortie calculée par la méthode selon l'une des revendications 1 à 7. 13) Utilisation de la méthode selon l'une des revendications 1 à 7, pour un simulateur de spectroscopie d'impédance. 14) Utilisation de la méthode selon l'une des revendications 1 à 7 pour simuler le comportement thermique et électrique d'une batterie de traction de véhicule hybride ou électrique, ou d'une installation de fourniture d'énergie photovoltaïque ou éolien. 15) Utilisation de la méthode selon l'une des revendications 1 à 7 pour le dimensionnement d'une batterie.30