FR 2980052 A1 20130315 FR 1158205 A 20110914 STABILISATION D'UN RESEAU ELECTRIQUE DC L'invention se rapporte à un procédé de stabilisation d'un réseau ou sous-réseau électrique DC à bus de tension continue, et par exemple d'un réseau ou sous-réseau HVDC (initiales de High Voltage Direct Current - courant continu haute tension). Ce réseau peut être embarqué dans un moyen de transport tel qu'un aéronef ou une automobile ou alimenter un bâtiment tel qu'une usine. Dans la présente demande, on entend par sous-réseau, un réseau local appartenant à un réseau électrique plus important. Un réseau ou sous-réseau DC comprend typiquement une source de tension continue alimentant plusieurs charges électriques, telles que des convertisseurs (du type par exemple onduleur, actionneur, source de stockage (super-capacité), machine synchrone à aimants permanents, etc.), qui sont branchées en parallèle aux bornes de la source de tension et qui sont chacune destinées à recevoir une consigne de courant ou de puissance. Chaque charge peut donc correspondre à un équipement différent et n'ayant pas forcément de lien de communication avec les autres charges, et peut avoir son contrôle actif de charge indépendant. Un réseau ou sous-réseau DC se différencie d'un réseau AC distribué par l'importance accrue de la contrainte liée à la stabilité du réseau à courant continu et de ses charges. L'ensemble des charges et le réseau lui-même doivent être normalement dimensionnés pour que leurs impédances respectives permettent d'avoir un réseau stable. Un réseau DC est instable à un point de fonctionnement si le signal de tension ou de courant aux bornes de sa charge oscille et a une amplitude qui augmente au cours du temps. Il est connu de stabiliser un réseau ou une charge de manière active ou passive. La stabilisation passive d'une charge d'un réseau consiste à dimensionner son impédance de manière à satisfaire les critères de 30 stabilité du réseau. Ceci se traduit par des surdimensionnements importants des éléments passifs (résistances et capacités) du réseau pour permettre la stabilisation du réseau suivant les différentes configurations. Cette solution n'est pas toujours satisfaisante car elle peut entraîner une augmentation de la taille et de la masse du réseau, ce qui présente un problème important lorsque ce réseau doit être embarqué, en particulier dans un aéronef. Par ailleurs, il est possible que ce surdimensionnement ne respecte pas certaines contraintes prévues initialement comme une fréquence de coupure d'un filtre ou autre. Par ailleurs, l'addition d'une charge à un tel réseau stabilisé de manière passive, peut entraîner une instabilité du réseau, même si la charge ajoutée au réseau est elle-même stable. Ceci montre que la stabilité des charges seules n'assure pas la stabilité du réseau comportant ces charges, l'interconnexion des charges sur un même bus pouvant conduire à l'instabilité de celui-ci. Une autre manière passive d'assurer la stabilité d'un réseau est de s'assurer, via l'intégration de saturations sur les consignes de référence des charges, que celles-ci ne consomment pas plus d'une certaine puissance sur le réseau. La stabilisation active d'une charge sur un réseau dédié consiste à modifier sa consigne de commande afin d'augmenter la stabilité du réseau. Cependant, cette solution est limitée car elle s'applique à un réseau simple 20 comportant une seule charge. Il est également connu de stabiliser un réseau à plusieurs charges de façon centralisée. Le réseau comprend alors un organe de stabilisation qui va acquérir des informations de l'ensemble du réseau et générer un vecteur contenant les signaux stabilisants de chaque charge et ainsi 25 stabiliser l'ensemble du réseau. Cependant, toutes les solutions de stabilisation active d'une charge ou d'un réseau électrique de la technique actuelle ne sont pas entièrement satisfaisantes car on a constaté qu'en cas de défaut d'une des charges du réseau ou en cas de reconfiguration du réseau, la stabilisation du réseau 30 n'est plus assurée. L'invention a notamment pour but d'éviter ces inconvénients de façon simple, efficace et économique. Elle a pour objet un procédé de stabilisation active et décentralisée d'un réseau ou sous réseau DC. Elle propose à cet effet un procédé de stabilisation d'un réseau électrique DC, ce réseau comportant une source de tension continue alimentant des charges électriques qui sont branchées en parallèle aux bornes de la source de tension et qui sont chacune destinées à recevoir une consigne de courant ou de puissance, caractérisé en ce que le réseau est stabilisé par régulation des consignes appliquées aux charges au moyen d'une impédance virtuelle de stabilisation générée aux bornes de chaque charge, ces impédances virtuelles étant dimensionnées pour stabiliser le réseau dans différents points de fonctionnement désirés et dans différentes configurations données du réseau comprenant l'état où au moins une charge est inactive ou en panne et l'état où la stabilisation d'une charge est inactive, chaque impédance virtuelle étant générée au moyen d'une boucle de régulation non linéaire agissant sur la consigne de la charge correspondante. Selon une autre caractéristique de l'invention, chaque boucle de régulation utilise une loi non linéaire de régulation générant un signal de sortie du type K.f(x), x étant une variable de commande de la charge et K étant un coefficient de correction propre à cette charge, les coefficients de correction (K1, K2, ..., Kn) étant déterminés en effectuant les étapes consistant à : a) définir un modèle mathématique du réseau électrique, ce modèle comportant des équations comprenant les coefficients de correction (K1, K2, ..., Kn) précités pour assurer la stabilité des charges et du réseau ; b) définir des contraintes à respecter pour que le réseau reste stable, notamment lors de reconfiguration du réseau ou de défaut d'au moins une charge ; et c) déterminer un vecteur optimal de coefficients de correction au moyen d'un algorithme d'optimisation destiné à vérifier le respect des contraintes précitées pour un vecteur donné et à calculer une fonction objectif de stabilisation avec ce vecteur, puis à répéter ces étapes jusqu'à ce qu'un vecteur satisfasse aux contraintes à respecter. Le procédé selon l'invention permet d'assurer la stabilité d'un réseau ou sous-réseau DC par la stabilisation des charges de ce réseau ou sous-réseau, et ce pour des points de fonctionnement définis quelle que soit la configuration du réseau. La stabilisation de chaque charge prend en compte l'impact de chaque charge sur la stabilisation du réseau global. Le procédé selon l'invention permet de moins ou de ne pas surdimensionner les éléments passifs du réseau en créant une impédance (par exemple une capacité ou une résistance) virtuelle en amont de chaque charge permettant la stabilisation des charges et de l'ensemble du réseau. Cette impédance virtuelle a pour but d'amortir les perturbations éventuellement générées sur le réseau ou aux bornes de la charge et est « dimensionnée » de manière à stabiliser le réseau dans les différents points de fonctionnement désirés, dans les différentes configurations réseau données, et dans les cas où au moins une charge est en défaut (c'est-à-dire qu'elle est désactivée ou en panne). L'impédance virtuelle générée aux bornes de chaque charge est obtenue par la régulation de la consigne de commande de la charge au moyen d'une loi non-linéaire et d'un coefficient de correction prédéterminé. Un vecteur comportant les coefficients de correction optimaux de toutes les charges est déterminé pour que les contraintes précitées soient respectées et qu'une fonction objectif soit réalisable. L'étape a) peut par exemple consister à définir un modèle non linéaire du réseau et à le linéariser autour d'un point de fonctionnement ou à exprimer son jacobien. Les contraintes à respecter à l'étape b) peuvent comprendre des plages de valeurs tolérées pour les coefficients de correction et le maintien de la stabilisation en cas de reconfiguration du réseau, de défaut de stabilisation d'au moins une charge, et de variations ou imprécisions dans les valeurs paramétriques propres au réseau. La prise en compte des variations ou imprécisions dans les valeurs paramétriques propres au réseau permet de définir une marge de robustesse et de sécurité sur la stabilité. Avantageusement, le modèle mathématique comprend des équations caractérisant le réseau en l'absence d'une ou plusieurs charges et/ou en l'absence de la stabilisation d'une ou plusieurs charges. La stabilisation est ainsi robuste aux pannes et aux reconfigurations du réseau sans identification de la configuration en temps réel et sans modification des paramètres de régulation /stabilisation utilisés. L'algorithme d'optimisation utilisé à l'étape c) est par exemple destiné à comparer le calcul de la fonction objectif avec un vecteur au calcul de cette fonction avec un autre vecteur dont les coefficients de correction sont les plus faibles possibles, puis à répéter ces étapes jusqu'à ce que le vecteur considéré minimise la fonction objectif définie. Celle-ci est choisie de telle sorte à orienter les efforts de stabilisation sur le réseau selon la volonté d'un concepteur (par exemple, minime sur les actionneurs, maximale sur les sources de stockage). N'importe quel algorithme d'optimisation sous contraintes peut être utilisé (génétique, méthode du gradient, etc.). La présente invention concerne également un réseau électrique du type DC, comportant une source de tension alimentant des charges électriques qui sont branchées en parallèle aux bornes de la source de tension et qui sont chacune destinées à recevoir une consigne de courant ou de puissance, caractérisé en ce que chaque charge est stabilisée par la régulation de sa consigne au moyen d'un bloc de stabilisation ayant une entrée reliée à des moyens de mesure de la tension aux bornes de la charge et une sortie reliée à une borne d'application de la consigne de la charge par un sommateur relié également à des moyens d'émission d'une consigne. Selon l'invention, un bloc de stabilisation local est implanté dans la boucle de régulation de chaque charge. Le fait que le réseau comprenne plusieurs blocs de stabilisation locaux entraîne une diminution des câbles nécessaires aux transmissions d'informations (mesures, référence du signal stabilisant) par rapport à l'intégration d'un bloc de stabilisation unique dans le réseau (pour une stabilisation centralisée) qui doit être relié à un grand nombre de capteurs pour ramener toutes les mesures à ce bloc. L'intégration de blocs de stabilisation locaux au réseau permet également de diminuer le nombre de capteurs mis en jeu car seules les tensions aux bornes des charges sont nécessaires. Le bloc de stabilisation extrait la composante haute fréquence de la tension, c'est à dire qu'il doit filtrer sa composante continue. A titre d'exemple, ce filtrage peut être réalisé à l'aide d'un filtre passe-bas et d'un soustracteur, l'entrée du filtre étant reliée aux moyens de mesure précités et sa sortie étant reliée au soustracteur, qui est destiné à soustraire le signal de sortie du filtre au signal provenant des moyens de mesure. En variante, on peut utiliser des filtres passe-haut ou passe-bande. La présente invention concerne encore un moyen de transport, tel qu'un aéronef, caractérisé en ce qu'il comprend un réseau électrique stabilisé tel que décrit ci-dessus. L'invention sera mieux comprise et d'autres détails, avantages et caractéristiques de l'invention apparaîtront à la lecture de la description suivante faite à titre d'exemple non limitatif, en référence aux dessins annexés dans lesquels : - la figure 1 représente une architecture d'un réseau électrique DC comportant une source de tension continue et trois charges ; - la figure 2 représente un réseau électrique DC simplifié comprenant une source de tension et une seule charge qui est stabilisée activement par la création d'une impédance virtuelle Z à ses bornes ; - la figure 3 représente un réseau électrique DC similaire à celui de la figure 1 et comprenant des blocs de stabilisation locaux destinés à créer virtuellement des impédances pour stabiliser les charges ; - les figures 4a et 4b représentent chacune un bloc de stabilisation d'une charge d'un réseau DC ; - les figures 5 et 6 sont des schémas blocs représentant des étapes du procédé selon l'invention ; - la figure 7 est un graphe illustrant l'évolution de la tension de bus DC sur les trois charges du réseau de la figure 1 en réponse à un échelon de puissance, lorsque ces trois charges ne sont pas stabilisées ; et - les figures 8 et 9 sont des graphes illustrant l'évolution de la tension de bus DC sur les trois charges du réseau de la figure 3 en réponse à un échelon de puissance, lorsque ces trois charges sont stabilisées dans le cas de la figure 8, et lorsque seulement deux charges sont stabilisées, la fonction stabilisante d'une charge étant en défaut dans le cas de la figure 9. On se réfère d'abord à la figure 1 qui représente un réseau électrique DC 10 qui peut être embarqué dans un moyen de transport tel qu'un aéronef, ce réseau comprenant une source 11 de tension continue dont les bornes de sortie sont reliées en parallèle aux bornes d'entrée de trois charges différentes. La première charge 12 comprend une machine synchrone à aimants 25 permanents 14 (MSAP ou PMSM, acronyme de Permanent Magnet Synchronous Machine) alimentée par un onduleur triphasé 16 régulé en couple via une commande vectorielle d-q classique. La deuxième charge 18 comprend un convertisseur DC/DC 20 débitant sur une simple résistance 22, l'ensemble étant régulé en 30 puissance de sorte que l'on dispose ainsi d'une charge consommant une puissance parfaitement maîtrisée sur le réseau. La troisième charge 24 est une super-capacité (SC) qui est connectée au réseau via un autre convertisseur DC/DC 26 bidirectionnel en courant. L'ensemble est également régulé en puissance. On peut ainsi contrôler les échanges de puissance entre la SC et le réseau. La tension du bus DC est vn et les tensions aux bornes des trois charges sont respectivement v1, v2 et vn. Le réseau ou sous-réseau DC selon l'invention est par exemple un boitier du type EBAC, un sous réseau ELS (Landing Gear system) ou EFCS (Flight Control System), un boitier de distribution primaire type PRIMES, etc. Les charges du réseau 10 sont suffisamment bien asservies pour que l'on puisse les considérer comme fonctionnant à puissance constante. La tension d'entrée Ve est aussi supposée parfaitement constante. La stabilité d'une charge peut être étudiée individuellement avant de traiter le cas où les trois charges sont simultanément connectées sur le réseau. La stabilité de chaque charge peut être étudiée par l'étude des tracés fréquentiels de Nyquist de la FTBO des systèmes pour les configurations en cascade (1 charge). La stabilité du réseau global (3 charges) est appréciée par l'étude des zéros du dénominateur d'une fonction de transfert, qui est pour chaque charge la fonction de transfert du rapport de la tension vs aux bornes de la charge (vi, v2 ou vn) sur la tension d'entrée Ve (vs/ Ve). Les résultats obtenus peuvent être confirmés par simulation, comme on le verra dans ce qui suit. Comme décrit précédemment, les configurations électriques des réseaux peuvent être instables si la consigne de commande (par exemple de puissance) des charges dépasse un certain seuil. Cette limite de puissance est définie par la structure du réseau ainsi que son contenu, c'est à dire par la valeur des paramètres qui le composent. Il est connu que la valeur des éléments passifs comme les capacités, les résistances ou les inductances varient en fonction de plusieurs paramètres qui ne sont pas toujours maîtrisables et/ou constants. On peut prendre en exemple les variations de température ou le vieillissement qui agissent sur les propriétés électriques des systèmes. Tout ceci va donc entraîner une modification des propriétés de stabilité du système pouvant l'amener dans un état d'instabilité. Afin de contrer ce phénomène non désiré, l'invention propose de définir des commandes adaptées qui généreront des signaux dits stabilisants, qui viendront se superposer aux références des charges et permettront d'assurer totalement ou d'augmenter la stabilité du système. Ces signaux sont avantageusement nuls en régime permanent dans le but de ne pas perturber ou modifier les consignes souhaitées par un utilisateur ou un point de fonctionnement original. La modification des paramètres de consigne des charges entraîne une modification du point de fonctionnement du réseau. Autrement dit, pour chaque point de fonctionnement du réseau sont définis des paramètres de consigne des différentes charges. Par exemple, on définit un point de fonctionnement correspondant à iq = 7A, P2 = 500W, Psi = OW. Le procédé selon l'invention permet de moins ou de ne pas surdimensionner les éléments passifs du réseau en créant une impédance Z (capacité, résistance ou autre) virtuelle en amont de chaque charge permettant la stabilisation des charges et de l'ensemble du réseau (figure 2). Cette impédance virtuelle a pour but d'amortir les perturbations éventuellement générées sur le réseau ou aux bornes de la charge et est « dimensionnée » de manière à stabiliser le réseau dans les différents points de fonctionnement désirés, dans les différentes configurations réseau données, et dans les cas où au moins une charge est en défaut (c'est-à-dire qu'elle est désactivée ou est en panne). Selon l'invention, l'impédance de stabilisation de chaque charge est générée à partir d'un bloc de stabilisation 28 (figure 3) qui agit directement sur la boucle de consigne de la charge, pour réguler cette consigne et assurer la stabilité notamment en : - ajustant les variations de la consigne entre les différents points de fonctionnement demandés ; - assurant la stabilité du réseau lorsqu'au moins l'un des blocs de stabilisation n'est plus actif ; et - assurant la stabilité du réseau lorsqu'au moins l'une des charges est déconnectée (reconfiguration du réseau). Le bloc de stabilisation 28 de la figure 4a permet de générer une résistance virtuelle aux bornes d'une charge. Le bloc 28 comprend un moyen 30 de multiplication de la tension vs mesurée aux bornes de la charge par un gain ou coefficient de correction K (équivalent à 1/R) directement proportionnel à l'amplitude du signal stabilisant, pour produire un signal U. L'entrée d'un filtre passe bas 32 de fréquence de coupure wc est reliée à la sortie du moyen 30 et reçoit donc le signal U. Le signal de sortie X3 du filtre 32 est soustrait au signal U au moyen d'un soustracteur 34. Le signal de sortie P, du soustracteur 34 est additionné au signal P' de la consigne de la charge par un sommateur 36 dont la sortie est reliée à la charge. La variante de réalisation du bloc de stabilisation 28' de la figure 4b permet de générer une capacité virtuelle C, aux bornes d'une charge. Ce bloc 28' comprend également un filtre passe bas 32 de fréquence de coupure wc qui reçoit en entrée la tension vs du bus DC et dont la sortie est notée vsif. Un soustracteur 34 permet de retrancher le signal filtré vsif au signal vs. Le signal de sortie du soustracteur 34 correspond à une oscillation haute fréquence qui est multipliée par la fréquence de pulsation du filtre passe bas par un moyen 35 de multiplication de gain wc. Le signal de sortie du moyen 35 et le signal vsif sont transmis à un moyen 37 supplémentaire de multiplication pour produire un signal régulé de puissance pv, en fonction de la capacité virtuelle C, à générer aux bornes de la charge. Le modèle mathématique pour étudier la stabilité de la tension du bus DC d'une charge telle que représentée en figure 2 (pour lequel Ps = P so est la puissance d'entrée ou de consigne de la charge pour un point de fonctionnement veo, v ; - P so = Vso-ie0 étant supposée être constante), peut être exprimé de la façon suivante : L f - -ccliti e = -rf-ie-vs+Ve d P (1) C--VS - - v = i - , dt 8 evs dans lequel Lf représente une inductance, ie représente le courant délivré par la source de tension Ve, rf représente une résistance, vs représente la tension aux bornes de la charge, Ve représente la tension, C ou Cf représente une capacité qui forme avec l'inductance Lf un filtre LC. Pour simplifier ce modèle, il est considéré que la commande de la charge est efficace et donc que les variations de la tension du bus DC ont peu ou pas d'influence sur les courants commandés de la charge. L'analyse de la stabilité du signal de la charge représentée en figure 2 peut être effectuée en étudiant les pôles du modèle mathématique localement linéarisé autour d'un point de fonctionnement. Ceci amène aux critères suivants qui donnent la puissance de charge maximale admissible pour un point de fonctionnement stable : 517 (2) dans lequel vs, est la tension aux bornes de la charge pour le point de fonctionnement considéré. Dans le cas où la consigne de puissance (Pso) de la charge 30 atteindrait Pso max, cette charge ne serait plus stable. Pso max représente donc un seuil qui ne doit pas être atteint par la consigne Pso. Comme il ressort des équations ci-dessus, Pso max est directement proportionnel à la capacité C. L'augmentation de cette capacité, par l'addition d'une capacité virtuelle supplémentaire aux bornes de la charge, permet d'augmenter le 35 seuil Pso max- En pratique, cette capacité virtuelle est générée en injectant une puissance additionnelle P, à la consigne de puissance Ps, cette puissance additionnelle répondant à la formule suivante et représentant le signal de stabilisation de la charge délivré par le bloc de stabilisation (lorsque P, = 0, la charge n'est pas stabilisée) : (1«.. = v I y 8 dans lequel C, est la capacité virtuelle générée en amont de la charge par le bloc de stabilisation, et i, est le courant traversant cette capacité. Pour connaître l'influence du bloc de stabilisation sur l'impédance, il 15 est nécessaire d'intégrer le courant supplémentaire i' aux bornes de la charge, dans le modèle précité. Pour cela, le modèle intègre une dérivée filtrée de la mesure de la tension du bus DC (vs), ce modèle étant du type : - r Xi = - L. f 1.'80 ,} - , = + X2 + - 30 Ces équations représentent le modèle mathématique du réseau de la figure 2, stabilisé par un bloc de stabilisation, dans lequel les variables d'état sont : 35 (5) 40 avec les paramètres suivants du point de fonctionnement (e o, vso) (4) 25 1 j-S = -1-}.A - "..v....es.. = - (6) , . f et dans lequel, pour un filtre de premier ordre, vsif est une variable d'état additionnelle. Pour stabiliser une charge, on peut aussi proposer la loi de stabilisation suivante, dans laquelle K est le coefficient de correction à définir pour stabiliser la charge : U = K.vs2. La loi de stabilisation intègre un amortissement du signal et prend en compte le fait que les oscillations non amorties sont destinées à disparaître. Du fait que cet amortissement est seulement nécessaire autour d'un point de fonctionnement, il doit être suivi d'un filtre approprié, par exemple un filtre passe bas, comme représenté aux figures 4a et 4b. Le réseau de la figure 3 est modélisé et linéarisé à un point de fonctionnement donné. Ce réseau correspond au réseau de la figure 1 dans lequel des blocs de stabilisation selon l'invention ont été ajoutés, chaque bloc étant destiné à réguler la consigne de commande de la charge à laquelle il est associé de façon à la stabiliser au sein du réseau. Le modèle mathématique linéaire du réseau DC peut être défini par la matrice suivante (Msain), dans laquelle les paramètres K1, K2 et Ksc représentent respectivement les coefficients de correction de la première, de la deuxième et de la troisième charge. 0 o 0 t. 0 0 0 2K2 Une fonction objectif est déterminée pour stabiliser le réseau en fonction de contraintes qui seront décrites plus en détail dans ce qui suit. Cette fonction objectif peut être exprimée par l'équation suivante : 10 Les contraintes vont garantir la stabilité du système. Elles seront imposées aux parties réelles des valeurs propres X des matrices d'état du système dans ses différentes configurations. Les contraintes étant valables localement autour d'un point de fonctionnement, on étudiera la stabilité de tous les points de fonctionnement correspondant à un domaine D donné 15 par deux intervalles de puissance, tel qu'il est défini par la relation suivante : 49 25 et ceci pour les cas suivants : (a) au moins deux des stabilisateurs sont actifs sur les trois (possibilité de panne sur l'un des trois) ; (b) déconnexion de l'une des branches du système (déconnexion 30 automatique par sécurité). Dans les scénarios de défaut, on va donc prendre en compte l'évènement et ce qu'il implique sur le système par l'algorithme chargé de l'optimisation qui, s'il parvient à converger, donnera un dimensionnement des blocs de stabilisation qui assure la stabilité de l'ensemble des cas 35 considérés. On peut rappeler qu'assurer la stabilité locale pour tous les points de fonctionnement d'un domaine, ne correspond pas à assurer la stabilité large signal dans ce domaine. D'autres contraintes à prendre en compte concernent les cas de reconfiguration du réseau, par exemple en cas de défaut sur une branche du réseau. La matrice Mdéfautt suivante est une modélisation mathématique du réseau dans le cas où une branche de la première charge est ouverte, par exemple au point 01 en figure 3. 0 o 0 0 u E2 2 0 o 0 3," 0 0 0 La matrice Mdéf't2 suivante est une modélisation mathématique du 25 réseau dans le cas où une branche de la deuxième charge est ouverte, par exemple au point 02 en figure 3. ---- I 1 -1 L Ii I L1 0 î - .. ,,i. si) . 0 (- 1 ;)'-- :J f '1 (.-,,, 0 2A,.. ,.A1.,-.,,,, 0 o 0 0 2Ki1'103.,,i o!^. T. 0 0 0 0 o 0 o f z . La matrice Mdéfautsc suivante est une modélisation mathématique du réseau dans le cas où une branche de la troisième charge est ouverte, par exemple au point 03 en figure 3. (ï) t- 0 L o P-.., ,.-- o 0 ) (1) ,-,-- 0 u L. i\ .- O () ( ) ) 0 L2 0 () ''' ' 5 (t) 1.6 0 7-7- 0 0 ï ' O 0 `--- fï 0 0 1 -7--- 0 0 1.. 0 7,7 1 - (..j 0 0 r 0 0 --...5.'i Mi Viow; O 0 , C .2 V 7..t.F (2 O f) -tu, 2 K -, i''' -,T, '1.., -, Afin d'intégrer un critère de robustesse, on va imposer aux parties réelles des valeurs propres des matrices M -sain, Mdefautl, Mdefautsc et Mdefaut2 d'être inférieures à une valeur strictement négative notée MR. Ceci assurera une marge de sécurité pour prendre en compte les éventuelles variations et imprécisions sur les valeurs paramétriques. Finalement, les contraintes à respecter se traduisent par la relation 15 25 Toute autre contrainte exprimée sous forme d'équations peut être prise en compte dans l'algorithme. La linéarisation du système permet de définir les valeurs propres du système, en jouant sur les différentes valeurs de K, les parties réelles de chaque pôle sont rendues négatives (ce qui garantit la stabilité dans le domaine proche des points de linéarisation), et sont avantageusement rendues inférieures à -5 pour augmenter les marges de robustesse du système. Ensuite, en fonction par exemple des différents points de fonctionnement désirés, des différentes configurations possibles, et de critères de minimisation de l'influence des performances dynamiques sur les charges, les différentes valeurs de K répondant à l'ensemble de ces critères sont déterminées pour assurer la stabilité du réseau. La figure 5 représente très schématiquement les étapes du procédé selon l'invention de détermination des coefficients de correction qui sont appliqués aux consignes des charges pour les stabiliser au sein du réseau DC. Cette figure illustre les contraintes qui peuvent être prises en compte pour la définition des coefficients de correction. Dans l'exemple représenté, les contraintes prises en compte sont : a) des intervalles tolérés pour les coefficients K des blocs de stabilisation (bloc 60) ; b) une marge de robustesse sur la stabilité de chaque charge (bloc 62); c) la gestion des reconfigurations du réseau (bloc 64); et d) la gestion de la perte d'un des blocs de stabilisation (bloc 66). Le modèle du réseau établi ci-dessus permet de définir la fonction objectif à minimiser et d'apprécier l'impact de chaque charge sur la stabilisation (bloc 68). La modélisation du réseau, de ses points de fonctionnements, de ses modes dégradés tolérés souhaités permettent de définir un ensemble d'équations permettant de définir un unique jeu de valeurs de coefficients (K1 à Kn). L'intérêt de cette unique solution permet de faire en sorte que les valeurs définies sont robustes et uniques pour tous les cas souhaités et pris en compte dans la modélisation et le système d'équations. Un algorithme d'optimisation est ensuite utilisé pour itérer et trouver les meilleures valeurs de K répondant aux critères d'optimisation (blocs 70 et 72). L'algorithme de calcul et d'optimisation n'est pas spécifique, et peut par exemple être la fonction fmincon de matlab qui résout et optimise l'ensemble des équations. Ceci permet de définir les valeurs des coefficients de correction K précités. La figure 6 illustre schématiquement les étapes du procédé d'optimisation selon l'invention. Un utilisateur détermine tout d'abord un domaine D dans lequel les coefficients de correction K seront déterminés. L'algorithme d'optimisation va définir dans ce domaine un vecteur X, comportant l'ensemble des coefficients de correction des charges (bloc 80). Dans le cas présent d'un réseau DC à trois charges, X, = [K1, K2 et Ksc]. L'algorithme d'optimisation va dans un premier temps vérifier que toutes les contraintes précitées sont respectées pour le vecteur X, (bloc 82). Si les contraintes ne sont pas respectées, l'algorithme va définir un nouveau vecteur (84, 80). Si, au contraire, toutes les contraintes sont respectées. L'algorithme va calculer la fonction objectif avec le vecteur X, (bloc 86) puis va comparer le résultat obtenu à un résultat de cette même fonction lorsque le vecteur utilisé est le plus petit possible, c'est-à-dire que les coefficients de correction qui le composent sont les plus petits possibles (bloc 88). Si la fonction objectif fobj (X0) est supérieure ou égale à fobj (branches 90, 84), l'algorithme va calculer un nouveau vecteur et réitérer les étapes des blocs 80, 82 et 86. Lorsque la fonction objectif fobj (X0) est inférieure à fobj (Xmin) et que X0 est égal à Xmin, c'est-à-dire que le vecteur est le plus petit possible (bloc 92), l'optimisation est terminée (branche 94) et le vecteur identifié peut être considéré comme le vecteur d'optimisation (X0 = X,,,, = X0pt - bloc 96), ce vecteur d'optimisation comportant les coefficients de correction K1, K2 et Ksc optimaux pour la stabilisation de l'ensemble des charges et du réseau DC. L'efficacité de la stabilisation du réseau DC complet est vérifiée par simulation et expérimentation. Pour cela, un modèle Simulink du réseau peut être développé. Les courbes des figures 7 à 9 ont été obtenues par expérimentation. La figure 7 est un graphe illustrant l'évolution de la tension de bus DC (vi, v2 et vn) sur les trois charges du réseau de la figure 1 en réponse à un échelon de puissance, lorsque les trois charges ne sont pas stabilisées. On constate que, suite à l'échelon de puissance, les oscillations des signaux ne se stabilisent pas. Au contraire, leurs amplitudes augmentent. Lorsque les valeurs propres de la fonction de transfert vn/Ve de chaque charge sont étudiées, on constate qu'au moins une partie des zéros du dénominateur de cette fonction de transfert ont une partie réelle positive, ce qui prouve l'instabilité du système. La figure 8 est un graphe illustrant l'évolution de la tension de bus DC (vi, v2 et vn) sur les trois charges du réseau de la figure 1 en réponse à un échelon de puissance, lorsque ces trois charges sont stabilisées par le procédé selon l'invention. Au contraire de la figure précédente, les oscillations des signaux se stabilisent et leurs amplitudes diminuent au cours du temps, en réponse à un échelon de puissance. Le réseau DC est donc stabilisé. Lorsque les valeurs propres de la fonction de transfert vn/Ve de chaque charge sont étudiées, on constate que les zéros du dénominateur de cette fonction de transfert sont tous à partie réelle négative, ce qui prouve la stabilité du réseau quand les trois branches fonctionnent simultanément. La figure 9 est un graphe correspondant à la figure 8 et pour lequel le bloc de stabilisation de la charge SC est mis en défaut et n'assure donc plus sa stabilisation. On constate toutefois que les oscillations des signaux se stabilisent et leurs amplitudes diminuent au cours du temps, en réponse à un échelon de puissance. Le réseau DC reste donc stabilisé malgré la désactivation d'un bloc de stabilisation. REVENDICATIONS1. Procédé de stabilisation d'un réseau électrique (10) DC, ce réseau comportant une source de tension continue (11) alimentant des charges électriques (12, 18, 24) qui sont branchées en parallèle aux bornes de la source de tension et qui sont chacune destinées à recevoir une consigne de courant ou de puissance, caractérisé en ce que le réseau est stabilisé par régulation des consignes appliquées aux charges au moyen d'une impédance virtuelle de stabilisation générée aux bornes de chaque charge, ces impédances virtuelles étant dimensionnées pour stabiliser le réseau dans différents points de fonctionnement désirés et dans différentes configurations données du réseau comprenant l'état où au moins une charge est inactive ou en panne et l'état où la stabilisation d'une charge est inactive, chaque impédance virtuelle étant générée au moyen d'un bloc de stabilisation d'une boucle de régulation non linéaire qui agit sur la consigne de la charge correspondante et qui est apte à générer un signal de stabilisation pv(t) destiné à être superposé à la consigne de la charge pour augmenter le seuil de cette consigne au-delà duquel le réseau est instable, le signal de stabilisation étant du type pv(t) = vs(t).Cv.dvs/dt ou pv(t) = K.vs2- X3 dans lesquelles vs représente la tension aux bornes de la charge, C, et K sont des paramètres définissant l'impédance virtuelle de stabilisation et X3 est le signal de sortie d'un filtre de fréquence de coupure wc qui reçoit en entrée le signal K.vs2, plus le paramètre C, ou K est important et plus le signal de stabilisation est important, ce qui augmente le seuil de la consigne de la charge et garantit la stabilité du réseau. 2. Procédé selon la revendication 1, caractérisé en ce que chaque boucle de régulation utilise une loi non linéaire de régulation générant un signal de sortie du type K.f(x), x étant une variable de commande de lacharge et K étant un coefficient de correction propre à cette charge, les coefficients de correction (Ki, K2, ..., Kn) étant déterminés en effectuant les étapes consistant à : a) définir un modèle mathématique du réseau électrique, ce modèle comportant des équations comprenant les coefficients de correction (Ki, K2, ..., Kn) précités pour assurer la stabilité des charges et du réseau ; b) définir des contraintes à respecter pour que le réseau reste stable, notamment lors de reconfiguration du réseau ou de défaut d'au moins une charge ; et c) déterminer un vecteur (X0) optimal de coefficients de correction au moyen d'un algorithme d'optimisation destiné à vérifier le respect des contraintes précitées pour un vecteur donné et à calculer une fonction objectif de stabilisation avec ce vecteur, puis à répéter ces étapes jusqu'à ce qu'un vecteur dont les coefficients de correction sont les plus faibles possibles satisfasse aux critères. 3. Procédé selon la revendication 2, caractérisé en ce que les contraintes à respecter comprennent des plages de valeurs tolérées pour les coefficients de correction et le maintien de la stabilisation en cas de reconfiguration du réseau, de défaut de stabilisation d'au moins une charge, et de variations ou imprécisions dans les valeurs paramétriques propres au réseau. 4. Procédé selon la revendication 2 ou 3, caractérisé en ce que le 25 modèle mathématique comprend des équations caractérisant le réseau en l'absence d'une ou plusieurs charges et/ou en l'absence de la stabilisation d'une ou plusieurs charges. 5. Procédé selon l'une des revendications 2 à 4, caractérisé en ce que l'algorithme d'optimisation est destiné à comparer le calcul de la fonction 30 objectif avec un vecteur au calcul de cette fonction avec un autre vecteur dont les coefficients de correction sont les plus faibles possibles, puis àrépéter ces étapes jusqu'à ce que le vecteur considéré soit identique au vecteur dont les coefficients sont les plus faibles possibles. 6. Réseau électrique (10) du type DC, comportant une source de tension continue (11) alimentant des charges électriques (12, 18, 24) qui sont branchées en parallèle aux bornes de la source de tension et qui sont chacune destinées à recevoir une consigne de courant ou de puissance, caractérisé en ce que le réseau est stabilisé par régulation des consignes appliquées aux charges au moyen de blocs de stabilisation (28) implantés dans des boucles de régulation non linéaire, chaque bloc de stabilisation étant destiné à générer une impédance virtuelle de stabilisation à une charge et ayant une entrée reliée à des moyens de mesure de la tension aux bornes de la charge correspondante et une sortie reliée à une borne d'application de la consigne de la charge par un sommateur (36) relié également à des moyens d'émission d'une consigne. 7. Réseau selon la revendication 6 caractérisé en ce que le bloc de stabilisation (28) comprend un filtre passe bas (32) et un soustracteur (34), l'entrée du filtre étant reliée aux moyens de mesure précités et sa sortie étant reliée au soustracteur qui est destiné à soustraire le signal de sortie du filtre du signal provenant des moyens de mesure. 8. Moyen de transport, tel qu'un aéronef, caractérisé en ce qu'il comprend un réseau électrique (10) selon la revendication 6 ou 7.