Brevet FR2803664A1 — publié le 20010713

La présente invention est relative aux batteries d'accumulateurs électrochimiques au lithium. Plus particulièrement, l'invention concerne un procédé pour déterminer l'état de charge d'une batterie d'accumulateurs au lithium pendant son utilisation notamment dans un véhicule automobile de type hybride. Dans de telles applications en particulier, il est important de connaître à tout moment 'état de charge de la batterie, afin de pouvoir le maintenir, par une régulation appropriée, à une valeur prédéterminée qui est fixée en général à la moitié de la charge totale de la batterie. De la sorte, pendant que le véhicule automobile circule, il est possible de maîtriser aussi bien les phases motrices pendant lesquelles de l'énergie est puisée dans la batterie que les phases de freinage pendant lesquelles de l'énergie peut être réinjectée dans la batterie. Pour assurer cette maîtrise, il est souhaitable de pouvoir déterminer l'état de charge de la batterie avec une précision inférieure à 1%. On connaît essentiellement deux méthodes de détermination de l'état de charge d'une batterie d'accumulateurs. La première dont on trouve une description notamment dans US 5 352 968 est une méthode estimative qui consiste à mesurer l'état de charge à partir plusieurs mesures de tension. Ces mesures de tension sont corrigées afin de tenir compte de la résistance interne de la batterie. L'état de charge de la batterie peut ensuite être déterminé par interpolation à l'aide de cartographies de la tension à vide en fonction de l'état de charge de la température. Cependant, cette méthode statique permet pas une évaluation précise de l'état de charge au cours d'un cycle dynamique de la batterie. De plus, la tension à vide est une grandeur qui évolue souvent trop lentement en fonction de l'etat de charge pour permettre une mesure satisfaisante. US 5 444 378 décrit une autre méthode dite coulométrique de détermination de l'état de charge d'une batterie. Dans ce cas, on pose: dans laquelle EC est l'état de charge, Q la charge fournie par la batterie et Q",aX la charge maximale que batterie peut fournir. La charge Q est alors déterminée mesurant le courant I en l'intégrant par rapport au temps selon la formule. k étant un coefficient de rendement de la batterie. Cette méthode serait satisfaisante si on pouvait se contenter de déterminer l'état de charge sur une période de mesure brève, ce qui n'est pas le cas dans l'application qui vient d'être mentionnée plus haut. Si on applique cette méthode sur une période longue, les erreurs mesure du courant sont accumulées. Par conséquent, pour disposer dans ce cas d'un résultat satisfaisant, il faudrait périodiquement recaler la mesure en partant d'un état de charge connu, par exemple la pleine charge. C'est donc une méthode difficile à mettre en oeuvre dans le cadre d'une régulation qui doit être utilisée dans une voiture automobile hybride. En outre, le coefficient de rendement doit être déterminé expérimentalement et il faut le corriger dans le temps. FR 2 740 877 propose de déterminer l'état de charge d'une batterie en combinant les deux méthodes connues qui vient d'être brièvement décrites ci-dessus. Cependant, cette méthode mixte est difficile à mettre en oeuvre, nécessite également l'établissement expérimental de coefficients de rendement et la correction de ceux-ci dans le temps, mais surtout elle ne permet pas de tenir compte d'un phénomène propre aux batteries au lithium qui est celui de la relaxation. La présente invention a ainsi pour but de fournir un procédé de détermination de l'état de charge d'une batterie appliqué spécifiquement aux batteries d'accumulateurs au lithium qui permette, en partant de la technique décrite dans FR 2 740 877, d'améliorer la précision de la détermination de l'état de charge, en particulier par la prise en compte du phénomène de relaxation. L'invention a donc pour objet un procédé de détermination de l'état de charge une batterie d'accumulateurs au lithium consistant: a) à déterminer l'état de charge estimatif de la batterie par une méthode coulométrique basée sur l'intégration du courant extérieur échangé la batterie avec l'extérieur, et b) à recaler par une méthode estimative basée sur la loi d'Ohm, l'état de charge estimatif pour obtenir un état charge recalé, procédé étant caractérisé en ce qu il consiste en outre: c) à modéliser le phénomène de relaxation de la batterie d'une part en répartissant sa matière active en des particules microscopiques présentant une zone dite "de charge réservoir" contenant du lithium disponible instantanément pour le fonctionnement de la batterie, et une zone dite "de charge effective" contenant du lithium immédiatement disponible pour ce fonctionnement, et d'autre part en considérant que le flux d'échange entre les deux zones est lié à l'écart entre les états de charge des deux zones, et d) à corriger l'opération de recalage définie sous b) en tenant compte du flux échange ainsi modélisé pour faire intervenir le phénomene de relaxation dans la détermination de l'état de charge. Ainsi, par une correction continue du recalage de la valeur de l'état de charge estimatif, on arrive à déterminer de façon dynamique l'état de charge de la batterie avec une très bonne précision, le résultat de la mesure pouvant ainsi être utilisé dans un processus de régulation au cours duquel l'état de charge peut être ajusté constamment à une valeur de consigne. Le procédé selon l'invention est ainsi particulièrement approprié pour être applique à la gestion temps réel d'une batterie d'accumulateurs alimentant un véhicule automobile hybride. L'invention peut également présenter les caractéristiques complémentaires suivantes: - l'opération de correction de recalage est effectuée en calculant un terme correctif représentant l'écart entre un état de charge réservoir et un état de charge effectif respectivement desdites zones réservoir et de charge effective; - ledit terme correctif est choisi proportionnel audit écart; - ledit état de charge réservoir et ledit état de charge effectif sont tous deux recalés par application de la loi d'Ohm en fonction de l'état de charge effectif; - pour l'application de la loi d'Ohm en vue du recalage desdits états de charge, le résultat de l'intégration dudit courant extérieur, est filtré de manière à tenir compte du comportement dynamique de la batterie; - le filtrage est effectué en appliquant audit courant un coefficient de filtrage déterminé expérimentalement et un filtrage du premier ordre; - pour le recalage desdits états de charge, l'application de la loi d'Ohm consiste à faire la somme d'un terme représentant la tension à vide en fonction de l'état de charge effectif et d'un terme représentant le produit dudit courant extérieur par la résistance interne en fonction de l'état de charge effectif; - pour le recalage desdits états de charge, l'application de la loi d'Ohm consiste à ajouter à ladite somme un terme supplémentaire représentant le produit dudit courant filtré par une résistance interne complémentaire en fonction de l'état de charge effectif; - ledit terme supplémentaire est déduit dudit terme représentant le produit dudit courant extérieur par la résistance externe en fonction de l'état de charge effectif, multipliée par un coefficient de proportionnalité déterminé expérimentalement; - les termes impliqués dans l'application de la d'Ohm en fonction de l'état de charge effectif sont également établis en fonction de la température de la batterie; - l'intégration du courant est effectuée introduisant deux gains différents à savoir un gain pour le courant de charge et un gain pour le courant de décharge; - l'intégration du courant extérieur, est effectuée en prenant en compte un état de saturation à la charge et d'un état de saturation à la décharge, l'intégration de courant étant désactivée lorsque l'état de charge de la batterie s'approche respectivement desdits états de saturation; - la température de la batterie est calculée par l'intermédiaire d'une équation de chaleur faisant intervenir les échanges thermiques avec l'extérieur, les pertes Joules dues au courant extérieur et les pertes Joule dues au courant filtré; - le filtrage du courant extérieur est fait à partir de l'expression: dans laquelle i est le courant filtré, I le courant échangé par la batterie avec l'extérieur a est un coefficient de filtrage; - la température de ladite batterie est calculée à partir l'expression: dans laquelle T est la température de la batterie, text est la temperature ambiante, m est la masse de la batterie, c est la chaleur massique de la batterie, k est un coefficient d'échange thermique avec l'extérieur, R est la résistance interne, r est la résistance interne complémentaire, i est le courant filtré et I est le courant échangé par la batterie avec l'extérieur; - les états de charge réservoir et effectif sont calculés en résolvant les équations différentielles suivantes: et dans lesquelles - Sr est l'état de charge réservoir, - X, est un coefficient de proportionnalité, - Se est l'état de charge effectif, - kr et ke sont des constantes rapportant les grandeurs macroscopiques définissant la batterie et son fonctionnement à l'échelle microscopique desdites particules, - Umes est la tension mesurée aux bornes de la batterie, - Uo est la tension à vide de la batterie, - T est la température de la batterie, - R est la résistance interne, - Imes est le courant échangé par la batterie avec 'extérieur, - y est un facteur de proportionnalité, - i est le courant filtré, - g- et g+ sont des gains scalaires, - h désigne une fonction de saturation, - v est la valeur minimale de l'état de charge effectif, et - p est la valeur maximale de 'état de charge effectif. D'autres caractéristiques et avantages de la présente invention apparaîtront au cours de la description qui va suivre, donnée uniquement à titre d'exemple et faite en se référant aux dessins annexés sur lesquels: - la figure 1 est une vue schématique de la structure d'une batterie d'accumulateurs au lithium à laquelle le procédé selon l'invention peut être appliqué; - la figure 2 montre schematiquement la modélisation d'une particule active d' batterie au lithium; - la figure 3 est un schéma équivalent électrique une batterie au lithium; - la figure 4 est un graphe représentant les états de saturation d'une batterie au lithium; et - la figure 5 est un schéma illustrant les calculs effectués pour mettre en #uvre le procédé selon l'invention. figure 1 montre le schéma de la structure d' batterie au lithium. Celle-ci comprend une anode 1 à masse active 2, un séparateur 3 par lequel transitent des ions lithium + et une cathode 4 à masse active 5. Sur cette figure, on suppose que la batterie est en décharge alimente une charge 6 (un moteur de véhicule automobile hybride par exemple) par un circuit 7 dans lequel les électrons circulent selon les flèches f1 et le courant extérieur I dans le sens de la flèche f2. La batterie est enfermee dans une cuve 8 remplie d'électrolyte. La batterie représentée est du type à insertion dans laquelle l'anode et la cathode ont une structure poreuse enrichie de la matière active 2 et 5 respectivement. Selon l'invention, le phénomène de relaxation de la batterie est modélisé en considérant que sa matière active 2, 5 formée de particules sphériques 9 permettant de réaliser des réactions de stockage et de déstockage d'ions lithium+ Ces particules ont un rayon de l'ordre de la dizaine de microns. La figure 2 représente schématiquement l'une de ces particules sphériques 9. La modélisation selon l'invention consiste à décrire chaque particule comme étant formée d'une zone sphérique interne Zl et d'une zone en calotte sphérique Z2 dans lesquelles se produisent des phénomènes de diffusion radiale. La calotte extérieure matérialisant la zone Z2 contient le lithium immédiatement disponible pour le fonctionnement de la batterie, tandis que la sphère Z1 interne contient du lithium en tant que charge réservoir qui ne participe pas directement au fonctionnement de la batterie. On comprend que les réactions de stockage et de déstockage à partir de la zone Z2 s'accompagnent selon le cas de la libération ou de la capture d'électrons qui sont à l'origine du courant de charge ou de décharge de la batterie. Afin de prendre en compte la diffusion du lithium à l'interieur de chaque particule active 9, on definit les états de charge suivants: l'état de charge effectif Se (t) qui représente le taux d'utilisation de la zone extérieure Z2 de la particule 9. Il s'agit ici du lithium contribuant à l'échange avec l'exterieur de la particule; l'état de charge réservoir Sr(t) qui représente le taux d'utilisation de la zone interne Z1 de la particule 9. suppose également qu'à tout instant, 'ensemble des particules d'une électrode (anode et cathode stockent une quantité identique de lithium, répartie d' manière identique à l'intérieur de ces particules. Les flux d'échange de lithium entre les zones Z1 et Z2 et avec l'extérieur sont définis de la façon suivante: - j, (t) est le flux de lithium échangé la zone Z2 avec l'électrolyte, et jr(t) est le flux de lithium échangé entre les deux zones, étant entendu que le signe de ces grandeurs représente le sens flux considéré. On suppose que le flux de lithium s'écoulant entre les deux zones de chaque particule 9 est proportionnel à l'écart entre les états de charge dans les deux zones selon l'expression: j, (t) C', (t) Se (t)) ( dans laquelle J, est un facteur de proportionnalité. La valeur du courant filtré i de la batterie est obtenue en mesurant le courant externe I qui est filtré à l'aide d'un filtre du premier ordre mettant en oeuvre l'équation suivante: dans laquelle î est la valeur estimée (dans suite de la description, une notation semblable sera utilisée également pour toute autre valeur estimée), a étant coefficient de filtrage du premier ordre. La tension U de la batterie est déterminée à partir de la tension à vide Uo, des courants filtré i et externe I et de la résistance interne complémentaire r et interne R selon l'expression suivante: U(t) Uo(Se,T) R(se,T).I(t) r(.se,T).l(t) (5) dans laquelle T est la température de la batterie, et Un, R et r sont définis en fonction de l'état de charge effectif Se et de la température T, à partir de cartographies. La température T de la batterie peut être calculée à partir d'un bilan de puissance que l'on peut décrire par l'expression suivante: dans laquelle m est la masse de la batterie en g, c est la chaleur massique en J/g/ C, k est un coefficient d'échange thermique avec l'extérieur en W/ C, et T,Xt est la température ambiante en C. Dans cette expression, la tension à vide et les résistances interne et interne complémentaire sont déterminées comme dans l'expression précédente. A l'aide des différents éléments que l'on vient d'énoncer, on peut mettre en ceuvre les opérations de calcul nécessaires pour effectuer le recalage de la détermination estimative coulométrique de l'état de charge et la correction de ce recalage, en résolvant le système d'équations différentielles suivantes: et dans laquelle - Y est un facteur proportionnalité entre la résistance interne R la résistance interne complémentaire r (voir le schéma équivalent de la batterie figurant à la figure 3) selon l'expression r(Se,T) _ rR (S?, T) ; - kr et ke sont des constantes attachées respectivement aux zones Z1 et Z2 de chaque particule 9 de masse active, et permettant de rapporter les grandeurs macroscopiques qui sont le courant I, la tension U et les résistances R et r à 'échelle microscopique des particules. Ces deux équations différentielles mettent donc en oeuvre l'opération estimative coulométrique alliée à un recalage par application de loi d'Ohm (deuxième terme de l'expression 7 et quatrième terme de l'expression 8), ainsi qu'une correction de recalage par le premier terme des deux expressions afin de permettre la prise en compte de la relaxation. En outre, l'expression (8) contient deux termes qui rendent compte de l'état de saturation de la batterie. Lorsqu'on applique un courant de charge à une batterie pleine, le flux de lithium vers l'électrolyte j, doit être nul pour des raisons physiques évidentes. De même, lorsqu'on applique un courant de décharge à une batterie vide, le flux j, doit egalement être nul. Par conséquent pour tenir compte de effets de saturation, il vient: jn(t) lu - h(p Se())- h(t'I- +,u; -'rf,),(t)-V). h(t@+ dans laquelle g- et sont des gains scalaires strictement positifs, h(x) est une fonction de saturation représentée sur la figure 4 et v p sont respectivement des valeurs minimale et maximale l'état de charge effectif. IxI- et IxI+ désignent respectivement les parties négatives et positives de x: IxI_= min (x, 0) IxI =max (x, 0) Le graphique de la figure 4 représente une fonction de saturation h(x) adéquate qui est tracée en fonction de l'expression suivante: dans laquelle F,,, 1. La figure 5 représente un schéma fonctionnel d'un dispositif de calcul de l'état de charge d'une batterie au lithium mettant en oeuvre le procédé suivant l'invention et conduisant à résoudre équations différentielles (7) et (8). Ce dispositif peut être réalisé sous forme d'un circuit matériel ou être mis en oeuvre par un calculateur tournant sur un logiciel correspondant écrit pour remplir les fonctions décrites. tel logiciel pourra être mis au point par les spécialistes ayant pris connaissance de la présente description. Le dispositif de calcul comprend plusieurs blocs distincts indiqués en traits mixtes. Un premier 10 de blocs est destiné à calculer la valeur du courant i conformément à l'équation (4). Un second bloc 11 du dispositif de calcul est destiné à calculer la tension U conformément à l'expression Un troisième bloc 12 permet de calculer la température T de la batterie conformément à l'expression (6) . Un quatrième bloc 13 est conçu pour calculer les termes de saturation à la charge et à la décharge la batterie conformément aux deuxième et troisième termes de l'expression (8). Un cinquième bloc 14 est destiné, en utilisant les résultats calculés par les blocs 10 à 13, à déterminer les valeurs Se et Sr en résolvant les équations différentielles (7) et (8). Le bloc 10 comprend un sommateur 10a auquel est appliquée la valeur mesurée Imes du courant échangé par la batterie avec l'extérieur. La sortie de l'additionneur est appliquée à un amplificateur lob de gain a qui fournit sa valeur de sortie à un intégrateur 10c délivrant la valeur i. Celle-ci est réinjectée dans le sommateur 10a en étant retranchée de la valeur bloc 11 est destiné à calculer la valeur de la tension U. I1 comporte une première table cartographique lla dans laquelle sont mémorisées des valeurs de la résistance R en fonction de la valeur d'état de charge Se et la température T de la batterie, fournies respectivement par les blocs 14 et 12. De cette table est extraite la valeur R(Se, T). Celle-ci est appliquée à un multiplicateur 11b dans laquelle elle est multipliée la constante y. La valeur de sortie r (Se, T) ce multiplicateur représente donc la résistance interne complementaire r en fonction de l'état de charge Se et de la température T. valeur de sortie R (SE, T) de la table cartographique 11a est appliquée à un multiplicateur llc qui effectue la multiplication de cette valeur avec valeur In,,,. La valeur de sortie r (S,, T) du multiplicateur llb appliquée à un multiplicateur lld qui la multiplie par la valeur i calculée dans le bloc 10. Une table cartographique lle contient des valeurs Uo en fonction de l'état de charge Se et de la température T. I1 en est extrait la valeur Uo(SE, T). Celle-ci est appliquée, de concert avec les valeurs de sortie des multiplicateurs 11c et 11d, à un sommateur llf qui délivre ainsi la valeur de tension U sur sa sortie. Celle-ci est connectée à un autre sommateur llg qui retranche la valeur U de la valeur Unes, c'est à dire celle de la tension les bornes de la batterie. valeur calculée par le sommateur 11g est envoyée au bloc 14 dans lequel elle est d'abord appliquée à deux multiplicateurs 14a et 14b pour être multipliée par constantes kr et k.. Les sorties des multiplicateurs 14a et 14b delivrent donc respectivement les valeurs du deuxième terme de l'expression (7) et du quatrième terme de l'expression (8). Le bloc 12 comprend deux élévateurs au carré 12a et 12b qui calculent respectivement les valeurs I,52 et i2. sortie de l'élévateur au carré 12a est connectée à multiplicateur 12c dans laquelle la valeur I;nes2 est multipliée par la valeur R(SE, T) issue de la table 11a. La sortie de l'élévateur au carré 12b est connectée à un multiplicateur 12d dans laquelle la valeur i2 est multipliée par la valeur r(Se,T) issue du multiplicateur 11b. La valeur Te,t représentant la température ambiante est appliquée à un sommateur 12e dans laquelle la valeur de sortie T du bloc 12 est retranchée de la valeur Te,,t. La valeur calculée dans le sommateur 12e est appliquée à un multiplicateur 12f dans laquelle elle est multipliée par la constante k d'échange thermique de la batterie avec l'extérieur. Les valeurs issues des multiplicateurs 12c, 12d et 12f sont additionnées dans un sommateur 12g, puis la valeur somme ainsi calculée est appliquée à un multiplicateur dans laquelle elle est divisée par le produit des valeurs c qui représentent respectivement la masse de la batterie et la chaleur massique de cette dernière. La valeur calculée dans le multiplicateur 12h est intégrée par rapport au temps dans un intégrateur 12i qui delivre la valeur de température T. Le bloc 13 calcule, selon le signe du courant i (charge ou décharge de la batterie), les valeurs des deuxième et troisième termes de l'expression (8). A cet effet, il comprend deux sélecteurs 13a et 13b qui sélectionnent une voie de charge ou une voie de décharge selon que la polarité du courant i est négative ou positive. Les sélecteurs 13a et 13b transmettent, la valeur du courant i respectivement à des multiplicateurs 13c et 13d ou elle est multipliée selon le cas par la valeur h(1 Se) ou h(Se). Les résultats de ces calculs sont appliqués à deux autres multiplicateurs 13e et 13f qui les multiplient par les gains scalaires respectifs g - et g+, puis nouveau résultat est transmis par l'intermédiaire sommateur 13g au bloc 14. Outre les multiplicateurs 14a et 14b déjà décrits dernier comprend un sommateur 14c qui additionne la valeur issue du multiplicateur 14a à une valeur multiplicative engendrée dans un multiplicateur 14d. Ce dernier est relié par son entrée à un sommateur 14e dont il multiplie valeur de sortie par la constante X, (voir expression (3) . Le sommateur 14e est chargé de retrancher la valeur Sr calculée de la valeur Se également calculee dans le bloc 14. La sortie du sommateur 14c est connectée à un intégrateur 14f qui effectue une intégration par rapport au temps de la valeur calculée par le sommateur 14c pour délivrer la valeur Sn. Celle-ci est également appliquée à sommateur 14g dans laquelle est retranchée la valeur Le résultat calculé dans ce sommateur 14g est appliqué à un multiplicateur 14h dans lequel il est multiplié par le facteur @. La sortie du multiplicateur 14h connectée à un sommateur 14i qui additionne la valeur multipliée par le multiplicateur 14h à la valeur de sortie du bloc 13 et à la valeur de sortie du multiplicateur 14b. Le résultat calculé dans le sommateur 14i est appliqué à un intégrateur 14j qui l'intègre par rapport au temps pour délivrer la valeur Se. On comprend désormais que le procédé selon l'invention permet de réaliser une mesure en temps réel de l'état de charge d'une batterie au lithium en tenant compte du phenoméne de relaxation propre à ces batteries et ce avec précision satisfaisante, suffisante pour permettre l'utilisation d'une telle batterie dans un véhicule hybride Par une modélisation du comportement statique et dynamique de la batterie, le procédé selon l'invention permet de prendre en compte les phénomènes de pertes de capacité, notamment pour des courants de décharge élevés, et la saturation de l'état de charge, à batterie pleine ou vide. Le procédé permet en outre de résoudre facilement la limitation du courant, de la tension et de la puissance délivrés par la batterie. REVENDICATIONS 1. Procédé de détermination de l'état de charge d'une batterie d'accumulateurs au lithium consistant: a) à déterminer l'état de charge estimatif de la batterie par une méthode coulométrique basée sur l'intégration du courant extérieur (I) échangé la batterie avec l'extérieur, et b) à recaler par une méthode estimative basee sur la loi Ohm, l'état de charge estimatif (Se) pour obtenir un état charge recalé, ce procédé étant caractérisé en ce qu'il consiste en outre: c) à modéliser le phénomène de relaxation la batterie une part en répartissant sa matière active en des particules microscopiques (9) présentant une zone (Z2) dite "de charge réservoir" contenant du lithium non disponible instantanément pour le fonctionnement la batterie, une zone (Z1) dite "de charge effective" contenant lithium immédiatement disponible pour ce fonctionnement, et d'autre part en considérant que flux d'échange entre les deux zones (Z1 et Z2) est lié à 'ecart entre les états de charge (Se et Sr) des deux zones, d) a corriger l'opération de recalage définie sous b) en tenant compte du flux d'échange ainsi modélise pour faire intervenir le phénomène de relaxation dans la détermination de l'état de charge (Se). 2. Procédé suivant la revendication 1, caractérisé en ce que l'opération de correction de recalage est effectuée en calculant un terme correctif (X, (Sr-Se) ) représentant l'écart entre un état de charge réservoir (Sr) et un état de charge effectif (Se) respectivement desdites zones réservoir et de charge effective (Z1 et Z2). Procédé suivant la revendication 2, caractérisé en que ledit terme correctif (k. (Sr-Se) est choisi proportionnel audit écart. Procédé suivant l'une quelconque des revendications 1 à 3, caractérisé en ce que ledit état de charge réservoir (Sr) et ledit état de charge effectif (Se) sont tous deux recalés par application de la loi d'Ohm en fonction de l'état de charge effectif (Se). Procédé suivant l'une quelconque des revendications 1 à 4, caractérisé en que pour l'application de la loi d'Ohm en vue du recalage desdits états de charge (Se et Sr). le résultat de 'intégration dudit courant extérieur (I), est filtré de manière à tenir compte du comportement dynamique de la batterie; 6. Procédé suivant la revendication 5, caractérisé en ce que le filtrage est effectué en appliquant audit courant (I) un coefficient de filtrage (a) déterminé expérimentalement et un filtrage du premier ordre. 7. Procédé suivant la revendication 6, caractérisé en ce que, pour le recalage desdits états de charge (Se, Sr), l'application de la loi d'Ohm consiste à faire la somme d'un terme représentant la tension à vide (Uc) en fonction de l'état de charge effectif (Se) et d'un terme représentant le produit dudit courant extérieur (I) par la résistance interne (R) en fonction de l'état de charge effectif (Se) . 8. Procédé suivant la revendication 7, caractérisé en ce que, pour le recalage desdits états de charge (Se et Sr), l'application de la loi d'Ohm consiste à ajouter à ladite somme un terme supplémentaire représentant le produit dudit courant filtré (i) par une résistance interne complémentaire (r) en fonction de l'état de charge effectif (Se) Procédé suivant la revendication 8, caractérisé en ce que ledit terme supplémentaire est déduit dudit terme représentant le produit dudit courant extérieur (I) par la résistance interne (R) en fonction de l'état de charge effectif (Se) multiplié par un coefficient de proportionnalité (y) déterminé expérimentalement. 0. Procédé suivant l'une quelconque des revendications 7 à 9, caractérisé en ce que les termes impliqués dans l'application de la loi d'Ohm en fonction de l'état de charge effectif (Se) sont également établis en fonction de la température (T ) de la batterie. 1. Procédé suivant l'une quelconque des revendications précédentes, caractérisé en ce que l'intégration du courant filtré (i) est effectuee en introduisant deux gains différents à savoir un gain (h(1- Se)) pour le courant de charge et un gain h(Se)) pour le courant de décharge. 12. Procédé suivant la revendication 11, caractérisé en ce que l'intégration dudit courant filtré (' est effectuée en prenant en compte un état de saturation à la charge et d'un état de saturation à la décharge, l'intégration de courant étant désactivée lorsque l'état de charge (Se) de la batterie s'approche respectivement desdits états de saturation. 13. Procédé suivant l'une quelconque des revendications précédentes, caractérisé en ce la température (T) de la batterie est calculée par l'intermédiaire d'une équation de chaleur faisant intervenir les échanges thermiques avec l'extérieur, les pertes Joules dues au courant extérieur (I) et les pertes Joule dues au courant filtré (i). 14. Procédé suivant l'une quelconque des revendications précédentes, caractérisé en ce que le filtrage du courant extérieur (I) est fait ' partir de l'expression: dans laquelle i est le courant filtré, I le courant échangé par la batterie avec l'extérieur a est un coefficient de filtrage. 15. Procédé suivant la revendication 14, caractérisé en ce que la température (T) de ladite batterie est calculée à partir de l'expression: dans laquelle T est la température de la batterie, text est la température ambiante, m est la masse de la batterie, c est la chaleur massique de la batterie, k est un coefficient d'échange thermique avec l'extérieur, R est la résistance interne, r est la résistance interne complémentaire, i est le courant filtré et I est le courant extérieur. 16. Procédé suivant la revendication 15, caractérisé en ce que les états de charge réservoir et effectif (Sr et Se) sont calculés en résolvant les équations différentielles suivantes: et dans lesquelles - est l'état de charge réservoir, - est un coefficient de proportionnalité, - est l'état de charge effectif, k, et ke sont des constantes rapportant grandeurs macroscopiques définissant la batterie et fonctionnement à l'échelle microscopique desdites particules, Umes est la tension mesurée aux bornes de la batterie, - Uo est la tension à vide de la batterie, - T est la température de la batterie, - R est la résistance interne, - In,es est le courant échangé par la batterie avec l'extérieur, - y est un facteur de proportionnalité, - i est le courant filtré, - g- et g+ sont des gains scalaires, - h désigne une fonction de saturation, - v est la valeur minimale de l'état de charge effectif, et - p est la valeur maximale de l'état de charge effectif.