FR 3161309 A1 20251017 FR 2403935 A 20240416 La présente invention a pour objet un procédé de contrôle du courant d’un système de batterie (1) comprenant un élément électrochimique de type lithium-ion (10), le procédé comporte les étapes d’estimation d’une concentration de lithium en surface d’une électrode de polarité calculée à partir d’un modèle en éléments finis selon des pas de calcul pour un volume de grain discrétisé de la structure de l’électrode comprenant un volume de surface et au moins un volume interne prédéterminés et le contrôle du courant de l’élément électrochimique (10) comprenant la fourniture d’une première consigne de courant (IDC1) lorsque la concentration de lithium est inférieure à un niveau maximum autorisé et la fourniture d’une deuxième consigne de courant (IDC2) lorsque la concentration de lithium est égale ou supérieure au niveau maximum autorisé. Figure 1. PROCEDE DE CONTROLE D’UN COURANT DE BATTERIE EN FONCTION D’UNE CONCENTRATION DE LITHIUM CALCULEE PAR UN MODELE EN ELEMENTS FINIS Le domaine de l’invention concerne un procédé de contrôle de courant d’un système de batterie, notamment de véhicule électrifié. Les systèmes de batterie de puissance, notamment pour véhicule électrifié, comprennent classiquement un ensemble de cellules électrochimiques montées en série et/ou parallèle dans le système. On désigne par le terme élément de stockage d’énergie une cellule électrochimique ou un assemblage de cellules pouvant être montées en série et/ou en parallèle. Classiquement, en application d’électromobilité, les éléments de stockage d’énergie sont des cellules électrochimiques de type lithium-ion. En fonctionnement, ces cellules sont susceptibles de se dégrader si elles ne sont pas utilisées dans des conditions de fonctionnement nominal, c’est-à-dire dans les limites thermiques, de courants et de tension définies par les constructeurs de cellules. En effet, en cas de température trop basse ou de pilotage de courants trop importants, les électrodes d’une cellule peuvent se saturer d’atomes de lithium en surface, ce qui peut conduire à la formation de dendrites à l’origine de courts-circuits internes. La gestion des courants de charge et décharge des cellules doit être opérées dans des conditions sécuritaires tout en assurant une gestion énergétique optimale et un usage de la batterie le moins contraignant possible, notamment en durée de charge et autonomie de roulage. Le transport de lithium dans une électrode négative en graphite est généralement plus lent que dans les électrodes positives et, par conséquent, est le phénomène limitatif du courant de charge. Ce phénomène s’appelle déposition de lithium ou placage de lithium, désigné aussi par le terme « Lithium Plating » en anglais. Une concentration de lithium trop élevée peut entrainer la formation de dendrites et provoquer une augmentation de la résistance interne, voire un court-circuit. Durant la charge les ions lithium se réduisent en métal lithium sur la surface de la structure en graphite lors de l’échange d’électron. Le taux d’oxydoréduction des ions lithium est proportionnel au courant et la diffusion dans la structure de graphite suit les lois de diffusion. Le taux de diffusion dépend du gradient de concentration en métal lithium et d’un coefficient de diffusion, lequel dépend de l’état de charge, de la température et du vieillissement d’une cellule. La loi de Fick régit la diffusion des ions de lithium dans les différents composants de la batterie, tels que l'électrolyte et les électrodes. Cette loi est bien connue de l’homme du métier dans le domaine des batteries. On peut citer par exemple la référence bibliographique ci-après décrivant des principes d’application de la loi de Fick dans le domaine des cellules lithium-ion : « Wu Musheng, Xu Bo, Ouyang Chuying. Physics of electron and lithium-ion transport in electrode materials for Li-ion batteries. Chinese Physics B, 2016, 25(1): 018206 ». On peut citer en outre le document brevet EP-A1-2266181 décrivant un procédé de charge adaptative se basant sur une technique d’estimation d’une durée de diffusion et d’une concentration de lithium en surface à partir de la durée de diffusion, la capacité de la cellule et le courant de charge mesuré. Ensuite, à partir de la concentration, un courant de charge est calculé. Ce document enseigne de mesurer la durée de diffusion pendant une phase de relaxation à partir de la mesure d’une tension à vide suivant l’application d’un courant continu pendant une durée fixe. Cette technique d’estimation de la concentration de surface en lithium concerne les charges en branchement sur borne, couramment appelées charges « Plug-in ». Cette technique n’est pas donc applicable en freinage récupératif du fait de la dynamique de roulage et de l’impossibilité d’imposer une relaxation des cellules, en particulier dans le cas d’un véhicule à motorisation entièrement électrique. Or, pour améliorer la gestion énergétique, la majorité des véhicules électrifiés est équipée d’un système de freinage récupératif afin d’optimiser la gestion énergétique du véhicule. La gestion du groupe motopropulseur est alors prévue d’une part pour favoriser le freinage électrique récupératif et d’autre part pour piloter des valeurs maximales de courant de charge spécifiques aux situations de freinage tout en respectant les contraintes électrochimiques et sécuritaires liées aux cellules de stockage lithium-ion. On connait de l’état de la technique le document WO-A1-2022237476 qui enseigne une solution de calibration de tables d’impulsions de courant de charge à partir de modèles électrochimiques et thermiques établis par voie expérimentale. Ces modèles électrochimiques comportent également une loi de coefficient de diffusion des ions lithium au niveau des électrodes. Ce document enseigne un principe de commande par impulsion de courant d’une durée donnée, par exemple 10 secondes, dont une condition d’arrêt peut être dépendante de la tension de la cellule par rapport à un seuil de tension haut, d’une valeur de potentiel d’électrode entre l’électrode négative et le séparateur ou d’une température de cellule. Ces tables d’impulsion ont pour objectif de piloter des courants maximums respectant les contraintes électrochimiques des cellules, notamment au regard du dépôt de lithium. On connait en outre le document brevet WO-A1-2022032460 décrivant une solution comportant un procédé de détection de placage de lithium sur une électrode de cellule de batterie comportant la détermination d’une valeur de potentiel en temps réel de l’électrode négative et dans lequel, si ce potentiel devient inférieur ou égal à un seuil de zéro volt, le procédé détecte un placage de lithium. Ce document enseigne l’usage d’un modèle embarqué dans un véhicule pour déterminer ce potentiel en fonction du courant de charge, état de charge et température. Ce document enseigne un principe de commande du courant de charge visant à réduire la valeur de courant en fonction de ce potentiel. L’inconvénient d’un modèle de calcul en temps réel est qu’il est particulièrement énergivore. On connait en outre le document brevet EP-A1-3806270 un dispositif et procédé pour charger une batterie consistant à estimer une concentration en ion Lithium en surface d’une anode à partir d’une technique de réduction de l’ordre du modèle (« Model order reduction» en anglais) afin de réduire la complexité computationnelle des calculs. Dans les documents de l’état de l’art, les tables de contrôle de courant informent de valeurs pouvant être appliquées à partir d’une estimation de l’état de charge instantané et de la température d’une cellule, généralement pour un état relaxé de la cellule. Or, si bien l’estimation de l’état de charge d’une cellule peut être considérée comme l’image du niveau de saturation en lithium du cœur de la structure en graphite de l’électrode, le niveau de saturation de la zone de surface de l’électrode n’est pas toujours identique à celui du cœur de la structure en graphite, notamment en phase de roulage durant laquelle la dynamique des sollicitations électriques en charge et en décharge modifient le niveau de concentration des atomes de lithium en surface. Cela implique que, suite à une charge électrique en freinage récupératif, la surface puisse être à un niveau de saturation déjà critique, supérieur à celui qu’on estime à partir de l’état de charge (SOC pour « State of Charge » en anglais) de la cellule ou, suite à une décharge, puisse être à un niveau inférieur à celui estimé à partir de l’état de charge. Dans le premier cas, il apparait un risque de dépôt de lithium si on applique le courant découlant de la table obtenue par voie expérimentale. Dans le deuxième cas, on limite l’opportunité de récupération d’énergie. Il existe donc un besoin d’optimiser les solutions de contrôle de courant pour les phases de charge et de décharge d’une batterie. Un objectif de l’invention est de pallier les problèmes précités. Un autre objectif de l’invention est de proposer un estimateur en concentration de lithium requérant des ressources de calcul adaptées pour un calculateur embarqué de véhicule automobile. Un autre objectif de l’invention est de mettre en œuvre une solution de contrôle de courant optimisée augmentant les opportunités de commande du courant de charge maximum, en particulier dans les phases de freinage récupératif, tout en limitant le risque de déposition de lithium. Plus précisément, l’invention concerne un procédé de contrôle d’un système de batterie comprenant un élément électrochimique de type lithium-ion, le procédé étant mis en œuvre par une unité de commande dudit système de batterie et comportant la mesure de paramètres électriques instantanés de l’élément électrochimique comprenant la température, le courant et la tension dudit élément électrochimique. Selon l’invention, le procédé comporte en outre les étapes suivantes : - l’estimation d’une concentration de lithium en surface d’une électrode de polarité calculée à partir d’un modèle en éléments finis selon des pas de calcul pour un volume de grain discrétisé de la structure de l’électrode comprenant un volume de surface et au moins un volume interne prédéterminés, l’estimation comprenant à chaque pas de calcul : a) une première sous-étape de calcul d’une quantité réduite de lithium pour le volume de surface sous l’action d’un courant, notée Nredox(t), calculée à partir d’une équation d’oxydoréduction prédéterminée en fonction des paramètres électriques mesurés, b) une deuxième sous-étape de calcul d’une quantité diffusée de lithium depuis le volume de surface vers le volume interne, notée Ndiff(t), calculée à partir d’une équation de diffusion prédéterminée en fonction des paramètres électriques mesurés, c) une troisième sous-étape de calcul de ladite concentration de lithium en surface calculée selon la relation suivante : , CLiZ1(t-1) étant la concentration de surface du pas de calcul précédent et VZ1 étant le volume de surface, - le contrôle du courant de l’élément électrochimique comprenant la fourniture d’une première consigne de courant lorsque la concentration de lithium est inférieure à un niveau maximum autorisé et la fourniture d’une deuxième consigne de courant lorsque la concentration de lithium est égale ou supérieure au niveau maximum autorisé. Le procédé selon l’invention peut comporter les caractéristiques additionnelles suivantes, seules ou en combinaison : - L’équation d’oxydoréduction de lithium se base sur une quantité d’électricité mesurée pendant le pas de calcul en cours et la constante de Faraday selon la relation suivante : , où Nredox (t) est la quantité réduite exprimée en mol pour le pas de calcul en cours, I(t) est le courant instantané mesuré, Pdc est la durée prédéterminée du pas de calcul exprimée en seconde et F la constante de Faraday exprimée en coulombs par mol. - La quantité diffusée depuis le volume de surface vers le volume interne est calculée selon l’équation de diffusion de lithium suivante : , où Ndiff est la quantité diffusée exprimé en mol pendant la durée prédéterminée Pdc du pas de calcul, D est le coefficient de diffusion déterminé à partir d’une table enregistrée en mémoire de l‘unité de commande en fonction des paramètres électriques mesurés et exprimé en µm 2 par seconde, CliZ1(t-1) est la concentration en lithium du volume de surface pour le pas de calcul précédent exprimée en mol/µm 3 , CLiZ2(t-1) est la concentration en lithium du volume interne pour le pas de calcul précédent exprimée en mol/µm 3 , Ldiff est une distance de diffusion prédéterminée exprimée en µm représentative d’une distance entre le volume de surface et le volume interne. - La valeur de la deuxième consigne est calculée à chaque pas de calcul et correspond à une valeur du courant pour laquelle la quantité de lithium diffusée est égale à la quantité réduite de lithium pour le volume de surface lorsque la concentration en lithium du volume de surface atteint le niveau maximum autorisé. - Le procédé comporte en outre, à chaque pas de calcul : - une sous-étape de calcul de valeurs simulées de concentration de lithium à un instant de fin d’une commande simulée d’un profil de courant constant ayant une durée d’application prédéterminée pour, respectivement, la première consigne et une troisième consigne de courant dont la valeur est comprise entre la première consigne et la deuxième consigne, - une sous-étape de sélection de la valeur de courant la plus élevée parmi au moins les première et troisième consignes dont la valeur simulée de concentration est inférieure ou égale au niveau maximum autorisé. - Le contrôle du courant comporte en outre, à chaque pas de calcul : - une sous-étape de calcul, de valeurs simulées de concentration de lithium à un instant de fin d’une rampe de transition graduelle de courant depuis la valeur du courant instantané mesuré vers la deuxième consigne, lesdites valeurs simulées étant calculées en fonction du courant de la rampe de transition, - et une sous-étape d’activation de la rampe de transition lorsqu’une des valeurs simulées atteint le niveau maximum autorisé. - La première consigne est une valeur de courant continu applicable pendant une durée prédéterminée fixe délivrée par une table de courant enregistrée en mémoire de l’unité de commande dépendante des paramètres électriques instantanés. Il est prévu en outre selon l’invention un système de batterie comprenant un élément électrochimique de type lithium-ion et une unité de commande configurée spécifiquement pour la mise en œuvre du procédé de contrôle de courant selon l’un quelconque des modes de réalisation précédents. Il est prévu en outre un véhicule à chaine de traction au moins partiellement électrifiée comportant un tel système de batterie. L’invention prévoit en outre une unité de commande comportant des moyens spécifiquement configurés pour mettre en œuvre le procédé de contrôle de courant selon l’invention, un programme-ordinateur comprenant des instructions qui, lorsque le programme est exécuté par l’unité de commande, conduisent celui-ci à mettre en œuvre l’un quelconque des modes de réalisation du procédé de contrôle de courant, et un support d'enregistrement lisible par ordinateur comprenant des instructions qui, lorsqu'elles sont exécutées par un ordinateur, conduisent celui-ci à mettre en œuvre le procédé de contrôle de courant selon l’invention. L’invention est une solution logicielle mettant en œuvre un estimateur de concentration de lithium en surface d’une structure de matériau actif d’une électrode d’un élément électrochimique de stockage d’énergie. L’estimateur exécute un modèle de calcul en éléments finis où les calculs de la concentration en lithium de la structure de l’électrode sont discrétisés en volumes de concentration homogène durant la durée de chaque pas de calcul. Cet estimateur améliore le contrôle du courant de la batterie et augmente les opportunités de recharge en freinage récupératif tout en protégeant efficacement contre les risques de déposition de lithium. D’autres caractéristiques et avantages de la présente invention apparaitront plus clairement à la lecture de la description détaillée qui suit comprenant des modes de réalisation de l’invention donnés à titre d’exemples nullement limitatifs et illustrés par les dessins annexés, dans lesquels : FIG. 1 représente schématiquement un système de batterie selon l’invention adapté pour la mise en œuvre du modèle en éléments finis de calcul de la concentration en lithium. FIG. 2 représente schématiquement une structure d’électrode illustrant le principe d’oxydoréduction et diffusion des ions lithium à travers le matériau actif de la structure d’électrode. FIG. 3 représente des graphiques de courbes de courant et de valeurs simulées de concentration de lithium pour des profils de courant constants selon l’invention. FIG. 4 représente des graphiques de courbes de courant et de valeurs simulées de concentration de lithium pour une rampe de transition de courant selon l’invention. FIG. 5 représente le procédé de contrôle de courant selon l’invention. FIG. 6 représente un véhicule électrifié prévu pour la mise en œuvre du procédé de contrôle de courant selon l’invention. L’invention s’applique à un système de stockage d’énergie comprenant un élément électrochimique de type lithium-ion et plus précisément à un procédé de contrôle de courant. L’invention trouve une application dans le domaine des véhicules comprenant une chaine de traction au moins partiellement électrifiée, entièrement électrique ou hybrides, de préférence les véhicules automobiles, mais pas seulement tels les aéronefs, tracteur, bicyclettes, navires. La FIG. 1 représente schématiquement un système de stockage d’énergie 1 comprenant un élément de stockage d’énergie 10 comportant au moins une cellule électrochimique 11 conformément à l’invention. L’élément de stockage 10 peut comporter une pluralité de cellules 11 ou groupement de cellules (appelé également cluster de cellules) reliées électriquement en parallèle et/ou en série. Une cellule électrochimique est un accumulateur d’énergie électrique ayant deux bornes et présentant une tension de quelques volts, le plus souvent comprise entre 2,3V et 4,2V, environ. Les cellules sont de type Lithium-ion (un oxyde de Nickel Manganèse Cobalt lithié NMC ou un phosphate lithié de fer LFP peuvent être cités à titre d’exemples de matières actives d’électrode positive). Plus précisément, une cellule Lithium-ion est composée principalement d’une électrode positive poreuse, une électrode négative poreuse, un séparateur et un électrolyte (pouvant être liquide, polymérique ou solide). Le principe de fonctionnement d’une cellule Lithium-ion repose sur l’échange réversible d’ions lithium entre les deux électrodes poreuses. Le système de stockage 1 comporte en outre des moyens de mesure, estimation ou capteur de paramètres électriques P parmi lesquels au moins un capteur de courant 12 apte à mesurer ou estimer une valeur de courant ICR en temps réel traversant l’élément de stockage 10, un capteur de tension et un capteur de température pour des cellules 11 ou chaque cellule 11 de l’élément de stockage d’énergie 10. Le système de stockage 1 comporte en outre une unité de commande 13 (appelée également ECU pour « Electronic Control Unit » ou BMS pour « Battery Management system »). L’unité de commande 13 comporte des informations de calibration de tables de courant continu pour une durée fixe, prenant en entrée des paramètres électriques comprenant l’état de charge et la température d’un élément de stockage d’énergie, d’une cellule ou d’un groupe de cellule. Ces tables de courant continu délivrent des valeurs de courant applicables pendant une durée fixée. Un protocole expérimental réalisé sur des cellules en conditions maitrisées peut être une modélisation numérique à partir des équations de concentration de lithium décrites ci-avant. Un autre protocole expérimental peut être la mesure d’un potentiel de référence de l’électrode négative correspondant à une tension mesurée entre l’électrode négative et une électrode dite de référence, généralement en lithium, connectée à l’électrolyte. Durant une charge, lorsque cette tension atteint 0 volt, cette condition indique que le niveau de saturation est atteint. Une première table délivre des valeurs d’un premier courant de charge continu applicables pendant la première durée au maximum en fonction de paramètres électriques mesurés. Les valeurs de courant sont applicables jusqu’à atteindre un niveau de saturation de dépôt de lithium en surface d’une électrode négative dudit élément de stockage en fin de la première durée. Les paramètres mesurés comprennent l’état de charge et la température d’un élément de stockage, pour une plage de [0%,100%] en points de charge et pour une plage de température [-20°C, 60°C] par exemple. On envisage d’obtenir plusieurs premières tables pour des durées applicables différentes par exemple des durées comprises entre 10 secondes et 180 secondes, par exemple une table pour une durée 10 secondes, une autre table pour une durée 30 secondes et une autre table pour une durée 180 secondes. Ce sont des tables de courant prévues pour une situation dynamique de charge en roulage d’un véhicule, telle le freinage récupératif par exemple. Ainsi, en référence à la FIG. 1 , l’unité de commande 13 comporte, enregistrée en mémoire au moins une première table de courant continu 18 IDC1. L’unité de commande peut comporter au moins deux tables 18 et 19 de courant continu IDC1 et IDC3 respectivement. Par exemple, la table 18 enregistre des valeurs de courants continus associés à une durée d’application de 10 secondes, et la table 19 enregistre des valeurs de courant continu associées à une d’application de 30 secondes. Il peut être prévu une autre table de courant, pour le pilotage d’un courant continu pour une durée d’application supérieure aux tables 18 et 19, c’est-à-dire supérieure à plusieurs minutes, comprise entre plusieurs minutes et plusieurs heures éventuellement. La valeur de courant continu délivrée par la deuxième table correspond à un rythme d’oxydoréduction de lithium proche de celui de la diffusion des atomes dans le cœur de structure d’électrode. Ce deuxième courant peut être piloté constamment, c’est-à-dire jusqu’à l’atteinte de la limite de tension d’une cellule correspondant à la tension dite de « Cut-off voltage » sans que le phénomène de dépôt de lithium apparaisse. Ce sont des valeurs de courant adaptées pour la recharge sur borne par exemple. Cette table délivre des valeurs de courant applicables sans risque de dépôt de lithium en fonction de paramètres d’état de charge et température d’un élément de stockage, pour une plage de [0%,100%] en points de charge et pour une plage de [-20°C, 60°C] par exemple. On notera que la deuxième table peut être remplacée en variante par un module 161 de calcul d’un courant continu IDC2 se basant sur le calcul de la concentration en lithium en temps réel mis en œuvre par l’invention. Ce module de calcul 161 est décrit plus en détail dans la suite de la description. Classiquement, l’unité de commande 13 du système de batterie comporte en outre un estimateur 14 de l’état de vieillissement des cellules fournissant un paramètre SOH représentatif de l’état de vieillissement prévu pour opérer un facteur de correction de courant en fonction du paramètre de vieillissement SOH. Le facteur de correction de courant est prévu pour corriger une valeur du courant IDC1, IDC3 délivrée par les première et deuxième tables 18 et 19 respectivement. Le paramètre d’état de vieillissement est une valeur représentative par exemple d’un nombre de cycles de charge et décharge, d’un nombre d’heures d’usage de la batterie, d’un kilométrage, ou d’un modèle dépendant d’une combinaison des facteurs. L’unité de commande 13 comporte en outre un estimateur 17 de l’état de charge instantané de la batterie 10 en fonction des paramètres P. Par exemple, l’estimateur détermine une valeur d’état de charge en fonction de la tension à vide mesurée, par exemple au démarrage du véhicule. Un objectif de l’invention est de piloter des courants maximums temporairement en phase de charge et pour cela il est nécessaire de connaitre une estimation de la concentration en lithium de la surface de la structure d’une électrode de polarité. L’invention se base sur un principe de modélisation du mécanisme de diffusion de lithium à travers la structure d’une électrode, en particulier la loi de Fick appliquée aux cellules au lithium. A cet effet, l’unité de commande 13 est munie d’un calculateur à circuits intégrés et de mémoires électroniques, le calculateur et les mémoires étant configurés pour exécuter le procédé de contrôle de courant selon l’invention. Mais cela n’est pas obligatoire. En effet, le calculateur pourrait être externe à l’unité de commande de la batterie, tout en étant couplé à cette dernière. L’unité de commande 13, selon l’invention, peut être réalisée sous la forme de modules logiciels (ou informatiques (ou encore « software »)), ou bien de circuits électroniques (ou « hardware »), ou encore d’une combinaison de circuits électroniques et de modules logiciels. L’invention propose un procédé de contrôle de courant se basant sur les propriétés physiques d’oxydoréduction de lithium en surface d’une électrode de polarité et de diffusion de lithium au sein de la structure interne de l’électrode. Plus précisément, le procédé met en œuvre un estimateur 15 d’une concentration de lithium en surface d’une électrode de polarité calculée à partir d’un modèle en éléments finis selon des pas de calcul pour un volume de grain discrétisé de la structure de l’électrode comprenant un volume de surface et au moins un volume interne prédéterminés. Le modèle peut comporter deux volumes internes au plus. Les calculs se réalisent par exemple selon des pas d’une durée d’environ 0,1 seconde, par exemple. Le volume de grain discrétisé correspond à un volume de grain représentatif de la structure du matériau actif. Le volume de grain est discrétisé en un volume cubique comprenant un volume de surface ayant une paroi S de surface d’exposition à l’oxydoréduction du lithium et une longueur d’insertion d’ion lithium L1 et un volume interne ayant une paroi S et une longueur d’insertion d’ion lithium L2. Les volumes sont séparés d’une longueur de diffusion des ions. A titre d’exemple de calibration du modèle, S est une valeur représentative du cumul des surfaces des grains de graphite, proportionnelle à la capacité, égale à 1,6.10 7 cm 2 , L1 est égal à 0,055µm et L2 est égal à 10µm. On considère une longueur de diffusion Ldiff égale à 3,2µm. Le modèle en éléments finis permet de réduire les ressources de calcul nécessaires pour l’estimation afin qu’il puisse être mis en œuvre par un calculateur embarqué d’un système de batterie. Il permet de simplifier le calcul de la concentration en lithium en surface de la structure en particule de graphite. En outre, on prend l’hypothèse que l’estimation du mécanisme de diffusion dans la structure de grain individuelle est représentative du mécanisme de diffusion dans la structure totale du matériau actif. Dans ce mode de réalisation, l’estimateur 15 est configuré pour calculer une quantité réduite de lithium pour le volume de surface sous l’action d’un courant, notée Nredox, calculée à partir d’un module 151 mettant en œuvre une équation d’oxydoréduction prédéterminée en fonction des paramètres électriques mesurés P. Plus précisément, l’équation d’oxydoréduction de lithium se base sur une quantité d’électricité mesurée pendant le pas de calcul en cours et la constante de Faraday selon la relation suivante : , où Nredox (t) est la quantité réduite exprimée en mol pour le pas de calcul en cours, I est le courant instantané mesuré, Pdc est la durée prédéterminée du pas de calcul exprimée en seconde et F la constante de Faraday exprimée en coulombs par mol égale à 96485 A/mol. En outre, l’estimateur 15 est configuré en outre pour calculer une quantité diffusée, notée Ndiff, d’ion lithium depuis le volume de surface prédéterminé vers le volume interne. Plus précisément, la quantité diffusée depuis le volume de surface vers le volume interne est calculée à partir d’un module 152 mettant en œuvre une équation de diffusion de lithium selon la relation suivante : . où Ndiff est la quantité diffusée exprimé en mol pendant la durée prédéterminée Pdc du pas de calcul, D est le coefficient de diffusion déterminé à partir d’une table enregistrée en mémoire de l‘unité de commande en fonction des paramètres électriques mesurés P et exprimé en µm 2 /s, CliZ1(t-1) est la concentration en lithium du volume de surface pour le pas de calcul précédent exprimée en mol/µm 3 , CLiZ2(t-1) est la concentration en lithium du volume interne pour le pas de calcul précédent exprimée en mol/µm 3 , Ldiff est une distance de diffusion prédéterminée exprimée en µm représentative d’une distance entre le volume de surface et le volume interne. Plus précisément, le coefficient de diffusion D est délivré par une table enregistrée en mémoire de l’unité de commande prenant en entrée la valeur du courant mesuré, du paramètre d’état de charge SOC, et de l’état de vieillissement SOH de la batterie. En particulier, l’état de vieillissement peut être utilisé pour appliquer un facteur de correction. On notera que dans le cas où le volume discrétisé comporte deux volumes internes ou plus, la quantité de lithium diffusée s’exprime selon la même formule entre chaque volume interne. L’estimateur 15 est configuré pour calculer en outre la concentration de lithium en surface à partir d’un module 153 selon la relation suivante : , CLiZ1(t-1) étant la concentration de surface du pas de calcul précédent et VZ1 étant le volume de surface avec VZ1=S*L1. On considère qu’au premier pas de calcul t=0s, CliZ1(0) est déterminé en fonction de l’état de charge instantané de la batterie. Ainsi, à l’instant d’initialisation de l’estimateur 15, par exemple au démarrage du véhicule, la valeur de la concentration CliZ1 (0) est délivrée par une table enregistrée en mémoire de l’unité de commande fonction de l’état de charge et de la température des cellules, et éventuellement d’une correction dépendante de l’état de vieillissement de la batterie. Cette table est configurée pour délivrée une valeur de concentration en fonction de valeur de l’état de charge fourni par l’estimateur 17. De même, à l’instant d’initialisation de l’estimateur 15, la valeur de CliZ2(0) est déterminée à partir de la même table enregistrée en mémoire de l’unité de commande. CliZ1(0) =CLiZ2(0). Ensuite, l’équation de diffusion permet d’estimer à chaque pas de calcul la valeur de la concentration de lithium du volume interne à partir d’un module 154 selon la relation suivante : CLiZ2(t-1) étant la concentration de surface du pas de calcul précédent, et VZ2 étant le volume interne avec VZ2=S*L2. Au final, à chaque pas de calcul d’une période prédéterminée Pdc, l’estimateur 15 estime la valeur de concentration en lithium au niveau du volume de surface à partir de la différence entre la quantité de lithium issu du mécanisme l’oxydoréduction et de la quantité diffusée depuis le volume de surface vers le volume interne, et éventuellement entre volumes interne. L’estimateur 15 permet d’estimer l’évolution de l’état de surface d’une structure de matériau actif d’électrode dans des phases de transitions de forte dynamique telle qu’une situation de freinage récupératif durant laquelle le courant de charge est susceptible d’augmenter la vitesse d’insertion d’ion lithium plus rapidement que la vitesse de diffusion au sein de la structure. Le modèle améliore le contrôle courant et la sécurité électrique. En FIG. 2 est illustré schématiquement le mécanisme d’oxydoréduction et de diffusion de lithium dans la structure d’électrode mis en œuvre par le modèle en élément fini selon l’invention. La structure de surface et la structure interne sont modélisées respectivement par le volume de surface ayant un volume VZ1, défini par la surface S et la longueur L1, et une concentration homogène en ion lithium et de plus, par le volume interne ayant un volume VZ2, défini par la surface S et la longueur L2, et une concentration homogène en ion lithium. Le flux d’oxydoréduction Jredox s’opère selon la quantité d’ion lithium réduite en surface et le flux de diffusion Jdiff s’opère entre les deux volumes selon une distance de diffusion Ldiff est entre ces deux volumes. L’état de surface est estimé en calculant la différence entre le flux d’oxydoréduction et le flux de diffusion. Les longueurs des volumes discrétisés L1 et L2, et la longueur de diffusion Ldiff sont des valeurs calibrées et fixes ou peuvent être des valeurs délivrées par des tables enregistrées en mémoire de l’unité de commande prenant en entrée les paramètres d’état de charge SOC et d’état de vieillissement SOH. En particulier, l’état de vieillissement peut être utilisé pour appliquer un facteur de correction. En retour à la FIG. 1 , l’unité de commande 13 comporte un contrôleur de courant 16 comportant une fonction de pilotage délivrant une information CC1 de contrôle de courant dont le pilotage est fonction de l’estimation de la concentration en lithium en surface de la structure d’électrode CLiZ1. L’information CC1 peut être une valeur de consigne de courant, une valeur de courant maximal autorisé à un instant donné ou une valeur de puissance électrique autorisée en charge ou en décharge. Le contrôleur de courant 16 comporte un module 161 délivrant l’information CC1 configuré pour fournir une première consigne de courant lorsque la concentration de lithium est inférieure à un niveau maximum autorisé et la fourniture d’une deuxième consigne de courant lorsque la concentration de lithium est égale ou supérieure au niveau maximum autorisé. Le niveau maximum est une concentration de lithium prédéterminée correspondant à un état de charge de la batterie compris entre 90% et 100% de SOC, par exemple 95% ou 100%. L’information CC1 est délivrée à travers des moyens de communication à d’autres systèmes électriques consommateurs et/ou générateurs d’énergie, par exemple convertisseur de tension, machine motrice électrique, chargeur. La première consigne de courant est la valeur du courant IDC1 délivrée par la table 18. Il peut être prévu que le contrôle de courant 16 soit configuré pour fournir une troisième consigne de courant correspondant à la valeur du courant IDC3 délivré par la table 19 lorsque la concentration de lithium est inférieure au niveau maximum autorisé. La deuxième consigne de courant est délivrée par une table de courant continu enregistrée en mémoire de l’unité de commande délivrant une valeur de courant dépendante de l’état de charge et de la température. En variante, la deuxième consigne de courant est délivrée par le module 162 délivrant une valeur de courant IDC2 calculée à partir de la concentration en surface CLiZ1. Plus précisément, le module 162 est configuré pour calculer la valeur de la deuxième consigne à chaque pas de calcul et correspond à une valeur du courant pour laquelle la quantité de lithium diffusée est égale à la quantité réduite de lithium pour le volume de surface lorsque la concentration en lithium du volume de surface atteint le niveau maximum autorisé. Plus précisément, dans cette situation Nredox(t)=Ndiff(t) pour une concentration de lithium du volume interne CliZ2 égale au niveau de saturation rapportée à 100%, avec : , Ainsi, la valeur de la deuxième consigne de courant IDC2 est obtenue selon la relation suivante : Par ailleurs, selon des variantes de l’invention, l’estimateur 15 peut comporter des modules de calcul de valeurs simulées de concentration de lithium afin de réaliser des projections de concentration pour des situations d’usage de la batterie. Ces projections améliorent le contrôle courant. Une première situation d’usage simulée est l’application de courants continus pendant une durée fixe. Ces projections permettent de vérifier qu’il est possible d’appliquer des valeurs de courant pendant une durée fixe sans atteindre une déposition de lithium en surface d’électrode, notamment pour opérer des freinages récupératifs sans changements de couple de freinage. A cet effet, l’estimateur comporte un module 155 configuré pour mettre en œuvre, à chaque pas de calcul, le calcul de valeurs simulées SCLi1 de concentration de lithium à un instant de fin d’une commande simulée d’un profil de courant constant ayant une durée d’application prédéterminée. Le calcul des valeurs simulées est exécuté pour, respectivement, la première consigne de valeur IDC1 et la troisième consigne de valeur IDC3 dont la valeur est comprise entre la première consigne et la deuxième consigne. En outre, le contrôle de courant 16 est configuré pour sélectionner la valeur de courant la plus élevée IDC1 et IDC3 parmi au moins les première et troisième consignes dont la valeur simulée de concentration est inférieure ou égale au niveau maximum autorisé. Ainsi, à l’instant de déclencher un courant, le contrôle de courant garantit de pouvoir appliquer cette valeur maximale constante pendant la durée associée à cette valeur. Cela évite une baisse indésirable ou inattendue du couple de freinage pendant une phase de freinage récupératif. L’agrément de conduite d’un véhicule est amélioré. En FIG. 3 , un graphique représente les valeurs simulées de concentration en fonction de trois projections dépendantes de trois valeurs de courants continues distinctes. Le graphique supérieur représente une situation de commande de courant où, sur l’axe des ordonnées, les valeurs instantanées du courant Ib sont représentées par un trait plein jusqu’à un instant t1 et des valeurs simulées de différents courants continus IDC1, IDC3 et IDC4 sont représentées en traits en pointillé. Par exemple, IDC1, IDC3 et IDC4 représentent des valeurs de courant applicables pendant une durée fixe de 10s finalisant à l’instant t2. Le courant IDC2 de courant constant est applicable de manière permanente sans atteindre la saturation de la surface et est représenté en double trait. Le graphique inférieur représente les valeurs de concentration en lithium CLi_IDC1, CLi_IDC3 et CLi_IDC4 pour le volume de surface calculées selon l’estimateur de concentration mis en œuvre par l’invention, pour respectivement chaque valeur de courant IDC1, IDC3 et IDC4. A l’instant t1, il peut être observé que parmi les valeurs simulées calculées pour l’instant t2 de finalisation de la durée d’application prédéterminée de 10s, la valeur CLi_IDC1 est supérieure au niveau maximum autorisé CLi_sat. Cette valeur de courant ne peut donc être retenue pour la consigne CC1. La valeur CLi_IDC3 et la valeur CLi_IDC4 sont toutes deux inférieures au niveau maximum autorisé CLi_sat, la fonction de contrôle de courant est configurée pour sélectionner le courant associé IDC3 de plus grande valeur. Une deuxième situation d’usage simulée est l’application d’une rampe de transition graduelle de courant depuis la première consigne vers la deuxième consigne de courant, appelée également rampe de « derating » en anglais. Cette rampe vise à sécuriser la transition de la première consigne vers la deuxième consigne de manière à prévenir un risque de déposition de lithium. A cet effet, l’estimateur comporte un module 156 configuré pour mettre en œuvre, à chaque pas de calcul, de valeurs simulées de concentration de lithium SCLi2 à un instant de fin d’une rampe de transition graduelle de courant depuis la valeur du courant instantané vers la deuxième consigne, lesdites valeurs simulées étant calculées en fonction du courant de la rampe de transition. En outre, le contrôle de courant 16 est configuré pour activer la rampe de transition lorsqu’une des valeurs simulées atteint le niveau maximum autorisé. En FIG. 4 , un graphique représente schématiquement en partie supérieure deux situations de contrôle de courant Iba et Ibb pour lesquelles une rampe de transition depuis le courant instantané Iba ou Ibb, vers la deuxième consigne IDC2. Le courant IDC2 de courant constant est applicable de manière permanente sans atteindre la saturation de la surface et est représenté en double trait. En partie inférieure, les valeurs de concentration en lithium pour le volume de surface calculées selon l’estimateur de concentration mis en œuvre par l’invention, respectivement pour le courant Iba et le courant Ibb, sont représentées par les courbes CLi_Iba et Cli_Ibb. Postérieurement à l’instant t1, les courbes en trait en pointillé représentent la rampe de transition depuis le courant instantané vers la deuxième consigne, la rampe IRPa pour le courant Iba et la rampe IRPb pour le courant Ibb. Les valeurs simulées CLi_RPa et CLi_RPb de la concentration en lithium pour le volume de surface associées respectivement aux rampes IRPa et IRPb sont également représentées en trait en pointillé. Pour la rampe IRPa, les valeurs simulées de concentration en lithium CLI_RPa n’atteignent à aucun instant le niveau maximum autorisé CLi_sat. Le contrôle de courant pour cette situation ne déclenche pas la rampe de transition. Pour la rampe IRPb, les valeurs simulées de concentration en lithium Cli_RPb calculées à l’instant t1 atteignent à l’instant projeté t2 le niveau maximum autorisé CLi_sat. Le contrôle de courant du système de batterie est alors configuré pour activer la rampe de transition à l’instant t1. Cette configuration du contrôle courant permet de prévenir la déposition de lithium pendant l’opération de la rampe de la transition. Il est prévu une variante combinant les valeurs simulées SCLi1 et SCLi2 calculées par les modules 155 et 156 de l’estimateur 15 dans laquelle les valeurs simulées de concentration de la rampe sont calculées à la suite de l’instant de fin de chaque profil de courant constant. Cette variante vise à garantir l’exécution de la transition entre le profil de courant constant et la deuxième consigne IDC2. En FIG. 5 , le procédé de contrôle de courant selon l’invention est représenté. Le procédé de contrôle de courant est mis en œuvre par l’unité de commande du système de stockage d’énergie d’un véhicule électrifié telle que décrite en détail en FIG. 1 . Ce procédé s’applique en phase de roulage d’un véhicule électrifié lorsqu’il est en déplacement à vitesse non nul ou en situation de recharge électrique de la batterie en branchement électrique à une borne de recharge. A une première étape E1, l’unité de commande est démarrée et initialisée. L’unité de commande mesure les paramètres électriques de la batterie, la température, la tension à vide et le courant notamment. Ces mesures permettent de déterminer un état de charge de la batterie. L’unité de commande initialise le modèle en éléments finis à un premier pas de calcul initial, déclenché par exemple par la détection de l’information « KEY-ON » du véhicule. Plus précisément, la valeur de concentration en lithium de chaque volume discrétisé du modèle en éléments finis mis en œuvre par l’estimateur de l’unité de commande est déterminée à partir de la table de valeurs enregistrée en mémoire de l’unité de commande et configurée à cet effet. La table délivre une valeur de concentration en fonction de la température et de l’état de charge pour le volume de surface prédéterminé CLiZ1(0) et pour le volume interne CLIZ2(0). Ensuite, le procédé comporte une deuxième étape E2 de mesure ou estimation des paramètres électriques P de l’élément de stockage de la batterie, les paramètres comportant au moins la température, la tension de chaque élément de stockage de la batterie, le courant échangé I(t), un paramètre d’état de charge, un paramètre d’état de vieillissement, notamment. La mesure est effectuée périodiquement durant le fonctionnement de la batterie, et au moins à chaque pas de calcul. Ensuite, le procédé comporte une troisième étape E3 d’estimation d’une concentration de lithium CLiZ1 en surface d’une électrode de polarité calculée à partir du modèle en éléments finis selon des pas de calcul pour un volume de grain discrétisé de la structure de l’électrode comprenant le volume de surface et au moins le volume interne prédéterminés. Plus précisément, le procédé comporte pour chaque pas de calcul une première sous-étape E31 de calcul d’une quantité réduite de lithium pour le volume de surface sous l’action d’un courant, notée Nredox(t), calculée à partir de l’équation d’oxydoréduction prédéterminée en fonction des paramètres électriques mesurés. L’équation d’oxydoréduction de lithium se base sur une quantité d’électricité mesurée pendant le pas de calcul en cours Pdc et la constante de Faraday selon la relation suivante : . Le procédé comporte en outre pour chaque pas de calcul une deuxième sous-étape de calcul E32 d’une quantité diffusée de lithium depuis le volume de surface vers le volume interne, notée Ndiff(t), calculée à partir de l’équations de diffusion prédéterminée en fonction des paramètres électriques mesurés. L’équation de diffusion de lithium est mise œuvre selon la relation suivante : . Il est rappelé en outre que la concentration du volume interne CLiZ2(t) est calculée selon la relation suivante : Le procédé comporte en outre pour chaque pas de calcul une troisième sous-étape de calcul de ladite concentration de lithium en surface CLiZ1(t) calculée selon la relation suivante : . Ensuite, à chaque pas de calcul, le procédé comporte pour le contrôle de courant du système de batterie, une étape de vérification E4 de la valeur de la concentration de lithium CLiZ1 par rapport au niveau maximum autorisé CLi_sat. Le niveau maximum autorisé peut être une valeur correspondant à un état de charge compris entre 90% et 100% de l’état de charge de la batterie, par exemple 95% ou 100%. Ensuite, le procédé comporte une étape de contrôle du courant E5 de l’élément électrochimique de la batterie comprenant la fourniture de la première consigne de courant lorsque la concentration de lithium est inférieure au niveau maximum autorisé CLi_sat et la fourniture de la deuxième consigne de courant lorsque la concentration de lithium est égale ou supérieure au niveau maximum autorisé CLi_sat. La première consigne de courant est un courant constant IDC1 délivré par une table de courant continu associée à une durée fixe, par exemple le courant IDC1 pour une durée 10 secondes ou le courant continu IDC3 pour une durée fixe 30s. La deuxième consigne de courant est le courant IDC2, soit délivré par une table de courant fonction de l’état de charge et la température mesuré, soit calculé en fonction de la concentration CLiZ1 comme décrit précédemment pour le module 162. De préférence, le procédé prévoit l’application d’une rampe de transition de courant, lorsqu’il est détecté que la concentration Cli_Z1 devient égale ou supérieure au niveau maximum autorisé. Selon des variantes du procédé, il est prévu que le procédé comporte des étapes de calcul des valeurs simulées de concentration de lithium pour les profils de courant constant IDC1, IDC3 et IDC4 décrits précédemment selon l’exemple de la FIG. 3 , et/ou pour la rampe de transition lorsque la limitation de courant est déclenchée. Ces étapes de simulation des concentrations permettent de prévenir une situation de déposition de lithium. Selon ces dernières variantes, le contrôle de courant est en mesure de sélectionner les courants de valeur la plus élevée et qui ont l’avantage d’éviter la déposition de lithium. En outre, ces simulations de concentration permettent d’assurer l’opération de la rampe de transition de courant avant d’atteindre le niveau maximum autorisé de concentration en lithium. En FIG. 6 , un véhicule électrifié 100 comportant un système de batterie 103 selon l’invention est représenté schématiquement. Le véhicule peut être un véhicule hybride ou à motorisation entièrement électrique, il comporte une traction au moins partiellement électrifiée. Le véhicule 100 comporte une machine motrice électrique 101 prévue pour déplacer le véhicule et alimentée par le système de batterie 103. Certaines ou toutes les fonctions du véhicule 100 sont contrôlées par une unité de commande 102. L’unité de commande 102 peut comprendre au moins un processeur qui exécute des instructions stockées dans un support lisible par ordinateur tel qu'une mémoire non volatile. L’unité de commande 102 peut également être de multiples dispositifs informatiques qui contrôlent des composants individuels ou des sous-systèmes du véhicule 100 d'une manière distribuée. Le processeur peut être n'importe quel processeur classique, tel qu'une unité centrale de traitement disponible dans le commerce. En variante, le processeur peut être un dispositif dédié tel qu'un circuit intégré spécifique à une application (ASIC) ou un autre processeur matériel. La machine motrice 101 et le système de batterie 103 peuvent être commandés en fonctionnement de freinage récupératif. Le système de batterie comporte une unité de commande 104 comprenant un calculateur et des mémoires configurés pour la mise en œuvre du procédé de contrôle de courant selon l’invention. En particulier, les mémoires ont un programme enregistré comprenant des instructions qui, lorsque le programme est exécuté par l’unité de commande du système de batterie, conduisent celui-ci à mettre en œuvre l’un quelconque des modes de réalisation du procédé de contrôle selon l’invention. En outre, l’unité de commande 104 est reliée en communication avec des calculateurs du véhicule au moyen d’un bus de communication embarqué 108, par exemple de type CAN (Controler Area Network » en anglais), pour délivrer l’information de contrôle de courant déterminée à partir du procédé selon l’invention. L’unité de commande 104 informe par exemple le groupe motopropulseur d’un courant maximal autorisé déterminé en fonction de la valeur de concentration en lithium calculée à partir du modèle en éléments finis conformément à l’invention. Le courant maximal autorisé est par exemple le courant de charge pilotable en freinage récupératif ou un courant de décharge pour alimenter la machine électrique. La valeur de courant est un courant continu qui est associé à une durée fixe maximal autorisée. Le véhicule comporte en outre un convertisseur de tension 105 de type DC/DC ou AC/DC. Le véhicule comporte en outre un relais ou commutateur électrique 106 prévu pour déconnecter et connecter sélectivement le système de batterie 103 des systèmes électriques du véhicule et, éventuellement d’un boitier de recharge 107. Le boitier de recharge est une interface permettant la connexion d’un câble de recharge à une source d’énergie externe reliée à un réseau d’alimentation électrique étendu fonctionnant en tension alternative. L’invention est décrite dans ce qui précède à titre d’exemple. Il est entendu que la personne de l’art est à même de réaliser différentes variantes de réalisation de l’invention en associant par exemple les différentes caractéristiques ci-dessus prises seules ou en combinaison, sans pour autant sortir du cadre de l’invention. Procédé de contrôle du courant d’un système de batterie (1) comprenant un élément électrochimique de type lithium-ion (10), le procédé étant mis en œuvre par une unité de commande (13) dudit système de batterie (1) et comportant la mesure (E2) de paramètres électriques (P) instantanés de l’élément électrochimique comprenant la température, le courant et la tension dudit élément électrochimique, le procédé étant caractérisé en ce qu’il comporte en outre les étapes suivantes: - l’estimation (E3) d’une concentration de lithium en surface d’une électrode de polarité calculée à partir d’un modèle en éléments finis selon des pas de calcul pour un volume de grain discrétisé de la structure de l’électrode comprenant un volume de surface et au moins un volume interne prédéterminés, l’estimation comprenant à chaque pas de calcul: a) une première sous-étape (E31) de calcul d’une quantité réduite de lithium pour le volume de surface sous l’action d’un courant, notée Nredox(t), calculée à partir d’une équation d’oxydoréduction prédéterminée en fonction des paramètres électriques mesurés, b) une deuxième sous-étape (E32) de calcul d’une quantité diffusée de lithium depuis le volume de surface vers le volume interne, notée Ndiff(t), calculée à partir d’une équations de diffusion prédéterminée en fonction des paramètres électriques mesurés, c) une troisième sous-étape (E33) de calcul de ladite concentration de lithium en surface calculée selon la relation suivante : , CLiZ1(t-1) étant la concentration de surface du pas de calcul précédent et VZ1 étant le volume de surface, - le contrôle (E5) du courant de l’élément électrochimique (10) comprenant la fourniture d’une première consigne de courant (IDC1) lorsque la concentration de lithium est inférieure à un niveau maximum autorisé et la fourniture d’une deuxième consigne de courant (IDC2) lorsque la concentration de lithium est égale ou supérieure au niveau maximum autorisé. Procédé selon la revendication 1 dans lequel l’équation d’oxydoréduction de lithium se base sur une quantité d’électricité mesurée pendant le pas de calcul en cours et la constante de Faraday selon la relation suivante : , où Nredox (t) est la quantité réduite exprimée en mol pour le pas de calcul en cours, I(t) est le courant instantané mesuré, Pdc est la durée prédéterminée du pas de calcul exprimée en seconde et F la constante de Faraday exprimée en coulombs par mol. Procédé selon la revendication 1 ou 2 dans lequel la quantité diffusée depuis le volume de surface vers le volume interne est calculée selon l’équation de diffusion de lithium suivante: , où Ndiff est la quantité diffusée exprimé en mol pendant la durée prédéterminée Pdc du pas de calcul, D est le coefficient de diffusion déterminé à partir d’une table enregistrée en mémoire de l‘unité de commande en fonction des paramètres électriques mesurés et exprimé en µm 2 par seconde, CliZ1(t-1) est la concentration en lithium du volume de surface pour le pas de calcul précédent exprimée en mol/µm 3 , CLiZ2(t-1) est la concentration en lithium du volume interne pour le pas de calcul précédent exprimée en mol/µm 3 , Ldiff est une distance de diffusion prédéterminée exprimée en µm représentative d’une distance entre le volume de surface et le volume interne. Procédé selon l’une quelconque des revendication 1 à 3 dans lequel la valeur de la deuxième consigne est calculée à chaque pas de calcul et correspond à une valeur du courant pour laquelle la quantité de lithium diffusée est égale à la quantité réduite de lithium pour le volume de surface lorsque la concentration en lithium du volume de surface atteint le niveau maximum autorisé. Procédé selon l’une quelconque des revendications 1 à 4 comportant en outre, à chaque pas de calcul : - une sous-étape de calcul de valeurs simulées de concentration de lithium (Cli_IDC1, Cli_IDC3) à un instant de fin d’une commande simulée d’un profil de courant constant ayant une durée d’application prédéterminée pour, respectivement, la première consigne (IDC1) et une troisième consigne de courant (IDC3) dont la valeur est comprise entre la première consigne (IDC1) et la deuxième consigne (IDC2), - une sous-étape de sélection de la valeur de courant la plus élevée parmi au moins les première et troisième consignes (IDC1, IDC3) dont la valeur simulée de concentration est inférieure ou égale au niveau maximum autorisé (CLi_sat). Procédé selon l’une quelconque des revendications 1 à 5 dans lequel le contrôle du courant comporte en outre, à chaque pas de calcul : - une sous-étape de calcul, de valeurs simulées de concentration de lithium (CLi_RPb) à un instant de fin d’une rampe de transition graduelle de courant depuis la valeur du courant instantané mesuré vers la deuxième consigne, lesdites valeurs simulées étant calculées en fonction du courant de la rampe de transition, - et une sous-étape d’activation de la rampe de transition (IRPb) lorsqu’une des valeurs simulées atteint le niveau maximum autorisé (CLi_sat). Procédé selon l’une quelconque des revendications 1 à 6 dans lequel la première consigne est une valeur de courant continu (IDC1) applicable pendant une durée prédéterminée fixe délivrée par une table de courant (18) enregistrée en mémoire de l’unité de commande dépendante des paramètres électriques instantanés (P). Système de batterie (1) comprenant un élément électrochimique de type lithium-ion (10) et une unité de commande (13) configurée spécifiquement pour la mise en œuvre du procédé de contrôle de courant selon l’une quelconque des revendications 1 à 7. Véhicule (100) à chaine de traction au moins partiellement électrifiée comportant un système de batterie (103) selon la revendication 8.