FR 3074916 A1 20190614 FR 1762093 A 20171213 Arrière-plan de l'invention L'invention concerne le domaine général de l'aéronautique. Elle se rapporte plus particulièrement à la surveillance d'un réseau de retour de courant utilisé par les équipements électriques 5 embarqués dans un aéronef tel que par exemple dans un avion. L'alimentation électrique des équipements électriques embarqués à bord des aéronefs (aussi désignés par « charges électriques ») est traditionnellement assurée par une source électrique, reliée aux appareils via un câble d'alimentation, et par un réseau de retour 10 de courant assurant une liaison basse impédance entre la source et les équipements électriques de l'aéronef. Une telle liaison basse impédance assure une chute de tension faible entre la source électrique et les diverses charges électriques embarquées dans l'aéronef, et par conséquent moins de pertes. 15 Traditionnellement, la fonction de retour de courant (ainsi qu'éventuellement d'autres fonctions complémentaires comme par exemple la métallisation, la protection au choc de foudre, la protection contre les champs rayonnés de forte intensité, la référence de potentiel, etc.) était assurée intrinsèquement par le fuselage et la structure 20 métalliques des aéronefs. Toutefois, depuis quelques années, l'aéronautique voit l'essor des matériaux composites qui sont de plus en plus utilisés dans les structures primaires et secondaires des aéronefs. Ces matériaux composites ont de très bonnes propriétés mécaniques, mais ne supportent pas la circulation de courants ni les échauffements que ceux-ci 25 peuvent engendrer. Il a donc été nécessaire de prévoir pour les aéronefs en matériaux composites, un réseau électrique complémentaire permettant d'assurer pour les équipements électriques embarqués à bord de ces aéronefs la fonction de retour de courant. Des exemples de tels réseaux 30 sont par exemple le réseau ESN (pour Electrical Structural Network) utilisé par Airbus et le réseau CRN (pour Current Return Network) utilisé par Boeing. Ces réseaux sont des réseaux basse impédance connectés en parallèle de la peau et de la structure composite de l'avion et qui permettent d'assurer les fonctions décrites précédemment. 35 Le réseau CRN est un réseau dédié aux fonctions électriques évoquées ci-dessus formé de câbles et de barres profilées interconnectés, 2 formant un réseau maillé. Le réseau ESN quant à lui est un réseau maillé et redondé utilisant des éléments combinant à la fois des fonctionnalités mécaniques et électriques et faisant partie de la structure primaire et secondaire de l'aéronef, des protections des câbles et des harnais, des tresses, des câbles, etc. Les technologies d'interconnexion utilisées actuellement pour connecter les éléments du réseau entre eux ne sont pas adaptées à la réalisation de contacts électriques fiables et viables. La vérification de chaque contact du point de vue électrique est en outre impossible en pratique. Les résistances de contact au niveau des interconnexions du réseau sont par conséquent assez mal maitrisées, dépendant de nombreux paramètres et sont amenées à évoluer dans le temps. Ainsi, dans l'exemple d'une jonction entre deux éléments réalisés au moyen de boulons, la résistance de contact dépend notamment de la force de serrage des boulons, du type de boulons, de l'état des surfaces, etc. Le réseau ESN se présente donc comme un réseau maillé de grande taille dans lequel les circulations des courants sont pilotées par les résistances de contact entre les différents éléments structurels du réseau. On comprend bien dès lors qu'il est nécessaire de surveiller plus particulièrement l'état des jonctions du réseau ESN réalisant l'interconnexion entre les différents éléments structurels sur lesquels s'appuie le réseau. Le document EP 2 653 859 décrit une solution permettant de surveiller une jonction flexible ou électromécanique entre deux éléments structurels d'un réseau ESN. Cette solution s'appuie sur l'utilisation d'un transformateur de courant monté en abaisseur de tension pour injecter, le temps du contrôle, un courant alternatif de forte intensité dans un conducteur électrique mis en parallèle de la jonction à contrôler. Cette solution repose sur l'hypothèse que le courant injecté passe en grande partie par la jonction que l'on cherche à surveiller. Or la plupart des résistances de contact au niveau des jonctions ne peuvent pas être mesurées directement, car elles sont connectées en parallèle au reste de la structure de l'aéronef. Une mesure de résistance donne en effet la valeur de la résistance de contact souhaitée en parallèle avec la résistance du reste du réseau, sans possibilité de discriminer entre les deux. 3 Il existe donc un besoin d'un procédé de surveillance facile à mettre en oeuvre et permettant de surveiller efficacement les éléments d'un réseau de retour de courant dans un aéronef, et plus particulièrement les jonctions utilisées par ce réseau pour l'interconnexion 5 des éléments structurels de l'aéronef sur lesquels il s'appuie. Objet et résumé de l'invention La présente invention répond notamment à ce besoin en proposant un procédé de surveillance d'un réseau électrique maillé de 10 retour de courant d'un aéronef, ledit aéronef comprenant au moins une structure en matériau composite, ledit réseau électrique comprenant une pluralité d'éléments conducteurs, ledit procédé de surveillance comprenant : - une étape d'identification d'au moins un élément critique du réseau, 15 pour au moins un couple source-charge comprenant une source et une charge de l'aéronef connectées au réseau ; - une étape d'estimation paramétrique d'une valeur de la caractéristique électrique de chaque élément critique, cette étape utilisant un modèle numérique modélisant un comportement électrique du réseau en 20 fonction de ladite caractéristique électrique et des mesures électriques réalisées au moyen de capteurs positionnés dans le réseau et permettant d'évaluer ledit comportement électrique, un capteur étant positionné dans le réseau pour chaque élément critique sur au moins un élément du réseau ou sur au moins un chemin du réseau formé 25 d'une pluralité d'éléments du réseau influencé par une perturbation de la caractéristique électrique de cet élément critique ; - une étape d'évaluation, pour ledit au moins un couple source-charge, d'un critère fonctionnel du réseau prédéfini, cette étape d'évaluation utilisant la valeur de la caractéristique électrique estimée pour chaque 30 élément critique ; et - une étape de détermination d'un état du réseau électrique en tenant compte du critère fonctionnel évalué pour ledit au moins un couple source-charge. Corrélativement, l'invention concerne aussi un système de 35 surveillance d'un réseau électrique maillé de retour de courant d'un aéronef, ledit aéronef comprenant au moins une structure en matériau 4 composite, ledit réseau électrique comprenant une pluralité d'éléments conducteurs, ce système de surveillance comprenant : - un dispositif d'identification, apte à identifier pour au moins un couple source-charge comprenant une source et une charge de l'aéronef 5 connectées au réseau, au moins un élément critique au regard d'un état de fonctionnement du réseau (normal ou dégradé, ou d'un état de connectivité) parmi les éléments du réseau ; - une pluralité de capteurs positionnés dans le réseau et aptes à réaliser des mesures électriques ; 10 - un module d'estimation paramétrique, configuré pour estimer une valeur de la caractéristique électrique de chaque élément critique, ce module utilisant un modèle numérique modélisant un comportement électrique du réseau en fonction de ladite caractéristique électrique et des mesures électriques réalisées au moyen desdits capteurs et 15 permettant d'évaluer le comportement électrique, un capteur étant positionné dans le réseau pour chaque élément critique sur au moins un élément du réseau ou sur au moins un chemin du réseau formé d'une pluralité d'éléments du réseau influencé par une perturbation de la caractéristique électrique de cet élément critique ; 20 - un module d'évaluation, configuré pour évaluer pour ledit au moins un couple source-charge, un critère fonctionnel du réseau prédéfini, ce module utilisant la valeur de la caractéristique électrique estimée pour chaque élément critique ; et - un module de détermination, configuré pour déterminer un état de 25 fonctionnement du réseau électrique en tenant compte du critère fonctionnel évalué pour ledit au moins un couple source-charge considéré. L'invention vise aussi un procédé d'identification d'éléments critiques d'un réseau électrique maillé de retour de courant d'un aéronef, 30 cet aéronef comprenant au moins une structure en matériau composite, le réseau électrique comprenant une pluralité d'éléments conducteurs. Ce procédé d'identification comprend : - une étape de modélisation du réseau électrique au moyen d'un graphe comprenant une pluralité de noeuds et une pluralité de branches reliant 35 ces noeuds, chaque noeud du graphe correspondant à un point de connexion entre deux éléments du réseau, chaque branche entre deux 5 noeuds représentant un élément du réseau, ledit graphe étant représenté par une matrice reflétant la topologie du réseau et comprenant des valeurs d'une caractéristique électrique des éléments du réseau représentés par les branches du graphe ; 5 - une étape de détermination des valeurs propres et des vecteurs propres de la matrice représentant le graphe ; - une étape d'estimation, à partir des vecteurs propres et des valeurs propres de la matrice, et pour chaque élément d'une pluralité d'éléments du réseau représenté par une branche du graphe, d'une 10 sensibilité d'une propriété électrique déterminée entre une source et une charge connectées au réseau électrique à une perturbation de la valeur de la caractéristique électrique de cet élément ; et - une étape de détermination, parmi les éléments du réseau, d'éléments critiques au regard d'un état de fonctionnement du réseau, un élément 15 du réseau étant déterminé comme critique lorsque la sensibilité de la propriété électrique estimée pour cet élément dépasse un seuil de criticité prédéterminé. L'invention vise également un dispositif d'identification de points de connexion critiques dans un réseau électrique maillé de retour de 20 courant d'un aéronef, ledit aéronef comprenant au moins une structure en matériau composite, ledit réseau électrique comprenant une pluralité d'éléments conducteurs, ledit dispositif d'identification comprenant : - un module de modélisation, configuré pour modéliser le réseau électrique au moyen d'un graphe comprenant une pluralité de noeuds 25 et une pluralité de branches reliant ces noeuds, chaque noeud du graphe correspondant à un point de connexion entre deux éléments du réseau, chaque branche entre deux noeuds représentant un élément du réseau, ledit graphe étant représenté par une matrice reflétant la topologie du réseau et comprenant des valeurs d'une caractéristique 30 électrique des éléments du réseau représentés par les branches du graphe ; - un module de détermination, configuré pour déterminer des valeurs propres et des vecteurs propres de la matrice représentant le graphe ; - un module d'estimation, configuré pour estimer à partir des vecteurs 35 propres et des valeurs propres de la matrice, et pour chaque élément du réseau représenté par une branche du graphe, d'une sensibilité 6 d'une propriété électrique déterminée entre une source et une charge connectées au réseau électrique à une perturbation de la valeur de la caractéristique électrique de cet élément ; et - un module de détermination, configuré pour déterminer parmi les 5 éléments du réseau, des éléments critiques au regard d'un état de fonctionnement du réseau, un élément étant déterminé comme critique lorsque la sensibilité de la propriété électrique estimée pour cet élément dépasse un seuil de criticité prédéterminé. Le procédé (respectivement le dispositif) d'identification selon 10 l'invention est utilisé préférentiellement par le procédé (respectivement par le système) de surveillance pour identifier les éléments critiques du réseau électrique de retour de courant. L'invention propose une méthode permettant de surveiller efficacement et simplement un réseau de retour de courant embarqué 15 dans un aéronef. Cette méthode repose sur une identification préalable des éléments critiques du réseau qui sont susceptibles notamment d'affecter la connectivité du réseau et sa capacité à remplir les fonctions pour lesquelles il a été dimensionné. De cette sorte, il est possible de simplifier 20 la surveillance du réseau en se limitant à la surveillance de ces éléments critiques. L'identification des éléments critiques proposée par l'invention a une complexité maîtrisée : elle s'appuie sur une représentation du réseau de retour de courant par un graphe, particulièrement bien adaptée pour 25 un tel réseau, et nécessite un unique calcul numérique consistant à déterminer le spectre de la matrice représentant le graphe (i.e. ses valeurs propres et ses vecteurs propres). A partir des valeurs propres et des vecteurs propres de la matrice représentant le graphe, il est alors possible, via une simple étude perturbative de ces valeurs propres et de ces 30 vecteurs propres, d'en déduire quels sont les éléments critiques du réseau à surveiller plus particulièrement pour s'assurer du fonctionnement normal de celui-ci. Ces éléments critiques sont identifiés en observant la sensibilité d'une propriété électrique du réseau (par exemple une chute de tension ou d'intensité dans le réseau) à des variations (perturbations) des 35 caractéristiques électriques des éléments du réseau. En particulier, si cette sensibilité dépasse un seuil prédéterminé lors d'une perturbation appliquée 7 à un élément donné, cet élément est considéré comme critique, car les variations de ses caractéristiques électriques sont susceptibles d'entrainer une défaillance du réseau. Dans un mode privilégié de réalisation : 5 - la caractéristique électrique considérée est une admittance de l'élément du réseau ; - la matrice du graphe est une matrice Laplacienne ; et - la propriété électrique est une tension électrique entre la source et la charge considérée (ou une chute de tension entre ces deux éléments). 10 Il convient de noter que la modélisation sous forme de graphe, ainsi que la décomposition spectrale sur laquelle s'appuie l'invention et l'étude perturbative du spectre de la matrice du graphe appliquée par l'invention permettent avantageusement d'obtenir des relations analytiques directes entre les perturbations appliquées aux 15 caractéristiques électriques des éléments du réseau et des composants du spectre du graphe, pour en déterminer ensuite la sensibilité des dites propriétés électriques du réseau face aux perturbations. Il en résulte un gain de temps substantiel par rapport aux méthodes classiques utilisant des calculs circuit et une complexité réduite de l'invention. 20 L'identification des éléments critiques du réseau permet de simplifier la surveillance du réseau, en se focalisant sur les éléments critiques. La surveillance proposée par l'invention permet de vérifier la connectivité des charges électriques embarquées dans l'aéronef et le respect de critères fonctionnels imposés au réseau. 25 Ce critère fonctionnel est par exemple une valeur maximale d'une chute de tension constatée entre chaque source et chaque charge, l'état du réseau électrique étant déterminé comme satisfaisant (ou normal) si cette chute de tension est, pour chaque couple source-charge, inférieure à ladite valeur maximale. 30 Un autre critère fonctionnel peut être une variation constatée entre la caractéristique électrique estimée pour un élément critique et une valeur nominale de cette caractéristique électrique, l'état du réseau électrique étant déterminé comme défaillant si pour au moins un élément critique du réseau pour au moins un couple source-charge, cette variation 35 est supérieure à un seuil d'alerte prédéterminé. 8 Bien entendu d'autres critères fonctionnels peuvent être considérés. L'invention permet ainsi de détecter facilement les dégradations d'éléments du réseau de retour de courant. Elle permet en outre de 5 localiser facilement les éléments dégradés et de faciliter les opérations de maintenance sur le réseau. Dans un mode particulier de réalisation, le procédé de surveillance comprend en outre une étape préliminaire de sélection, parmi une pluralité de charges de l'aéronef connectées au réseau électrique, 10 d'au moins une charge considérée comme essentielle au regard du fonctionnement de l'aéronef, les étapes d'identification, d'estimation paramétrique, et d'évaluation étant appliquées uniquement aux couples source-charge comprenant ladite au moins une charge essentielle sélectionnée. 15 Cette sélection préliminaire permet de réduire la complexité de la surveillance mise en oeuvre. Elle a en effet un impact sur le nombre d'éléments critiques considérés lors de la surveillance qui dépend des charges électriques prises en compte. Dans un mode particulier de réalisation, lors de l'étape 20 d'estimation paramétrique, on utilise autant de capteurs que d'éléments critiques identifiés, chaque capteur étant associé à un élément critique et placé sur l'élément du réseau ou sur le chemin du réseau le plus influencé par une perturbation de la valeur de caractéristique électrique de cet élément critique. 25 Ceci permet de limiter les mesures nécessaires pour estimer les caractéristiques techniques des éléments critiques lors de la surveillance. On s'assure dans ce mode de réalisation de disposer de mesures pertinentes pour évaluer de façon fiable ces caractéristiques techniques tout en limitant la complexité liée à la mise en oeuvre de ces mesures. Les 30 inventeurs ont constaté par le biais de tests que ce mode de réalisation conduit à un bon compromis entre complexité et fiabilité de la surveillance opérée, et permet d'optimiser le placement des capteurs dans le réseau. Bien entendu, il est possible d'utiliser davantage de capteurs et de réaliser davantage de mesures pour améliorer la précision de la 35 surveillance. 9 Dans un mode particulier de réalisation, l'étape d'estimation paramétrique met en oeuvre une technique d'optimisation dite des points intérieurs sur une pluralité d'itérations. Une telle technique est particulièrement efficace lorsque de 5 nombreux paramètres doivent être optimisés. Toutefois d'autres techniques d'optimisation peuvent être utilisées en variante. Dans un autre mode de réalisation, le procédé de surveillance comprend une étape préalable de caractérisation initiale des caractéristiques électriques à partir du graphe modélisant le réseau, de la 10 matrice représentant le graphe et de mesures électriques réalisées dans le réseau. Ainsi, la modélisation proposée par l'invention permet également la caractérisation initiale du réseau électrique de retour de courant et en particulier des caractéristiques électriques de ses éléments critiques 15 (résistances de contact). Brève description des dessins D'autres caractéristiques et avantages de la présente invention ressortiront de la description faite ci-dessous, en référence aux dessins 20 annexés qui en illustrent un exemple de réalisation dépourvu de tout caractère limitatif. Sur les figures : - la figure 1 représente, dans son environnement, un système de surveillance d'un réseau électrique de retour de courant conforme à l'invention ; 25 - la figure 2 représente, sous forme d'ordinogramme, les principales étapes d'un procédé d'identification d'éléments critiques du réseau de retour de courant, conforme à l'invention, tel qu'il est mis en oeuvre par le système de surveillance de la figure 1; - les figures 3A et 36 illustrent un premier exemple de réseau et sa 30 représentation par un graphe pouvant être utilisée lors du procédé d'identification représenté sur la figure 2; - les figures 4A et 4B illustrent un second un exemple de réseau et sa représentation par un graphe pouvant être utilisée lors du procédé d'identification représenté sur la figure 2 ; 10 - la figure 5 représente les principales étapes d'un procédé de surveillance conforme à l'invention mis en oeuvre par le système de surveillance de la figure 1; - la figure 6 représente un exemple de technique d'estimation 5 paramétrique mise en oeuvre lors du procédé de surveillance ; et - les figures 7 à 9 illustrent un exemple d'application de l'invention. Description détaillée de l'invention La figure 1 représente, dans son environnement, un système 10 de surveillance 1 d'un réseau électrique maillé 2 d'un aéronef, conforme à l'invention dans un mode particulier de réalisation. Dans l'exemple envisagé ici, l'aéronef est un avion comprenant des structures internes ou externes en matériaux composites. A titre illustratif, on envisage ici que le fuselage de l'avion est en matériau 15 composite, comprenant éventuellement une couche métallique (treillis métallique en surface des matériaux composites). Le fuselage est connecté aux autres éléments conducteurs du réseau par des boulons ou des rivets (qui sont eux aussi des éléments conducteurs). Le réseau électrique 2 considéré est un réseau de retour de 20 courant utilisé à bord de l'avion. Un tel réseau électrique assure, de façon connue en soi, le retour des courants fonctionnels des différentes charges électriques vers les différentes sources électriques de l'avion, mais peut également remplir une ou plusieurs des fonctions suivantes : - métallisation et protection contre les défauts électriques ; 25 - protection contre les effets indirects de la foudre ; - protection contre les champs rayonnés de forte intensité : - référence de potentiel pour les équipements et les électroniques embarqués à bord de l'avion ; et - protection des personnes à bord. 30 Il convient de noter qu'un réseau d'alimentation électrique d'un avion comporte plusieurs sources électriques (ex. sources de tension) : au moins un générateur électrique par moteur, un générateur électrique auxiliaire (APU - Auxiliary Power Unit), des batteries, un générateur associé à une turbine déployée en cas d'urgence (RAT - Ram Air Turbine). 35 A terme, d'autres sources d'énergie électriques pourraient être envisagées (piles à combustibles, générateurs en bout d'aile, etc.). 11 Par souci de simplification dans la suite de la description on considère une unique source électrique et plusieurs charges dans l'avion. Dans l'exemple envisagé ici, le réseau électrique 2 est un réseau ESN (Electrical Structural Network) tel que décrit par exemple dans le 5 document EP 2 653 859 précédemment cité. Il s'agit d'un réseau basse impédance, connecté en parallèle du fuselage en matériau composite de l'avion, et reliant les différentes charges électriques de l'avion aux différentes sources électriques. Une liaison basse impédance assure avantageusement une chute (i.e. différence) de tension faible entre 10 chaque charge et chaque source électrique de l'avion, et de ce fait moins de pertes. Le réseau électrique 2 est maillé et redondé, chaque élément du réseau électrique étant relié à au moins deux autres éléments du réseau afin de respecter notamment les normes de sécurité en vigueur. Il 15 s'appuie sur les composants métalliques de l'avion et comprend une pluralité d'éléments conducteurs de différentes natures, à savoir des éléments structurels constitués notamment de cadres métalliques, de traverses métalliques, de rails de sièges métalliques, du sol de cabine, etc., ainsi que de tout autre élément conducteur connecté en parallèle, 20 tels que des chemins de câbles métalliques (raceways en anglais), des protections électromagnétiques, des tresses, des câbles ou des barres, ou tout élément de jonction réalisant l'interconnexion ou la connexion entre ces différents éléments. Les jonctions entre les différents éléments conducteurs du réseau électrique 2 sont réalisées par des éléments de 25 fixation tels que des rivets, des boulons ou des vis. Les contacts entre les éléments du réseau électrique 2 peuvent être en conséquence des contacts complexes, réalisés à l'aide de plusieurs rivets, boulons ou vis, ou des contacts simples, réalisés à l'aide d'un seul rivet, boulon ou vis. Les contacts complexes concernent généralement les connexions 30 électromécaniques entre des éléments structurels du réseau électrique 2, tandis que les contacts simples sont généralement réalisés entre les autres éléments du réseau ou entre les autres éléments du réseau et des éléments structurels. On note que les segments d'éléments conducteurs du réseau 35 définis de part et d'autre le cas échéant des différentes jonctions 12 constituent eux-mêmes des éléments conducteurs du réseau au sens de l'invention. Quelle que soit la jonction considérée entre deux éléments du réseau (structurels ou pas), il apparaît, au niveau du contact une 5 résistance de contact. Les résistances induites par ces contacts électriques ne sont pas négligeables et impactent les circulations de courant dans le réseau électrique 2 et son état de connectivité. Dans l'exemple envisagé ici, pour faciliter la compréhension de l'invention, on s'attache uniquement au retour des courants continus 10 assuré par le réseau électrique 2. Le réseau électrique 2 est supposé ici être un réseau électrique résistif, chaque élément du réseau pouvant être assimilé à une admittance (définie comme l'inverse de l'impédance). Comme nous avons fait l'hypothèse d'être en courant continu, les parties inductives et capacitives de l'impédance sont nulles, donc l'impédance est 15 égale à la résistance et par conséquent, l'admittance est l'inverse de la résistance. Toutefois, ces hypothèses ne sont pas limitatives en soi et l'invention s'applique bien entendu à d'autres types de réseaux, et à d'autres caractéristiques électriques permettant de caractériser les 20 éléments du réseau. Ainsi, la prise en compte des courants alternatifs ou impulsionnels peut être envisagée. Cela implique un calcul fréquentiel, donc la prise en compte des inductances et des capacités des éléments du réseau, ainsi que des couplages inductifs et des capacitifs entre les différents éléments du réseau. La prise en compte des aspects 25 fréquentiels complexifie les modèles, mais ne change pas la démarche présenté dans cette étude. Le système de surveillance 1 est avantageusement apte à surveiller et à déterminer l'état du réseau électrique 2 lorsque celui-ci est en fonctionnement ou excité par des courants judicieusement choisis. Par 30 état du réseau, on entend ici son état de fonctionnement (par exemple, état normal/satisfaisant ou état défaillant) et sa capacité à assurer les fonctions primaires pour lesquelles il est prévu et dimensionné, et plus particulièrement sa fonction de retour des courants continus entre les charges et les sources électriques de l'aéronef. Cet état est évalué par le 35 système de surveillance 1 au regard d'un critère fonctionnel prédéfini représentatif du fonctionnement du réseau (normal ou dégradé et/ou de 13 son état de connectivité). Dans l'exemple envisagé ici, le critère fonctionnel considéré est une chute de tension maximale entre les sources et les charges de l'aéronef : autrement dit, toutes les chutes de tension constatées entre chaque source et chaque charge de l'aéronef doivent 5 être inférieures à une valeur maximale prédéterminée. Le critère fonctionnel est considéré ici comme respecté (et l'état du réseau est considéré comme satisfaisant) si cette chute de tension maximale est inférieure à la chute de tension maximale imposée par le système d'ordre supérieur. Par système d'ordre supérieur, on entend tout système 10 électrique utilisant tout ou une partie du réseau électrique 2 pour assurer le bon fonctionnement d'un équipement électrique en régime normal de fonctionnement ou en régime dégradé. Dans la suite de la description, à titre illustratif, on considère une chute de tension maximale de 0,5 V. Bien entendu, d'autres critères fonctionnels peuvent être 15 envisagés en variante par le système de surveillance 1 pour juger de l'état du réseau électrique 2, comme par exemple un critère fonctionnel portant sur une intensité de courant, sur une variation des résistances de contact au sein du réseau, etc. Pour surveiller l'état du réseau électrique 2, le système de 20 surveillance 1 s'appuie sur diverses entités. Plus particulièrement, il comprend : - un dispositif d'identification 3 conforme à l'invention, apte à identifier parmi les différents éléments du réseau 2, un ou plusieurs éléments dits critiques. Par élément critique, on entend ici un élément 25 susceptible de dégrader fortement la connectivité des charges de l'aéronef en cas de dysfonctionnement ou de dégradation. Dans l'exemple envisagé ici, on considère que les éléments critiques du réseau électrique 2 se situent au niveau des contacts assurés par les jonctions du réseau électrique 2, et non au niveau des éléments 30 structurels du réseau à proprement parler qui sont considérés comme fiables. L'identification de ces éléments critiques permet de simplifier la surveillance du réseau 2, qui est alors limitée à la surveillance de ces éléments critiques ; - une pluralité de capteurs 4 positionnés dans le réseau électrique 2 35 permettant de réaliser des mesures électriques. De tels capteurs sont par exemple des capteurs de champ (électrique ou magnétique), de 14 courant, de tension, etc. classiquement utilisés pour réaliser des mesures électriques ; et une pluralité de modules comprenant notamment : o un module d'estimation paramétrique 5: ce module est configuré 5 pour estimer les valeurs des caractéristiques électriques des éléments critiques. On s'intéresse ici aux valeurs des admittances des éléments critiques. Le module d'estimation paramétrique 5 utilise à cet effet, comme détaillé davantage ultérieurement, un modèle numérique modélisant le comportement électrique du 10 réseau 2 en fonction des admittances de ses éléments et, des mesures électriques réalisées au moyen des capteurs 4 dûment positionnés dans le réseau pour permettre l'évaluation du comportement électrique de celui-ci lorsqu'il est en fonctionnement normal, en régime dégradé ou alimenté sur un schéma précis ; 15 o un module d'évaluation 6: ce module est configuré pour évaluer pour différents couples source-charge, le critère fonctionnel auquel on s'attache pour déterminer l'état du réseau (chute de tension maximale inférieure à 0.5V ici). Il utilise à cet effet les valeurs des caractéristiques électriques des éléments critiques estimés par le 20 module d'estimation 5 ; et o un module de détermination 7: ce module est configuré pour déterminer un état du réseau électrique 2 en tenant compte du critère fonctionnel évalué par le module d'évaluation 6 sur l'ensemble des couples source-charge considérés. 25 Les fonctions des modules 5-7 du système de surveillance sont détaillées davantage ultérieurement. Dans le mode de réalisation décrit ici, les modules 5-7 sont des modules logiciels définis par des instructions d'un ou de plusieurs programmes d'ordinateurs stockés dans la mémoire morte d'un ordinateur 30 8, et lisibles par un processeur de cet ordinateur. Cet ordinateur peut être embarqué par exemple à bord de l'aéronef. Dans le mode de réalisation décrit ici, le dispositif d'identification 3 selon l'invention est également intégré dans l'ordinateur 8 et s'appuie sur les moyens matériels de cet ordinateur. Conformément à l'invention, le 35 dispositif d'identification 3 comprend divers modules logiciels ici, définis par des instructions de programmes d'ordinateurs, lisibles par le 15 processeur de l'ordinateur 8 et stockées dans sa mémoire morte. Ces instructions, lorsqu'elles sont lues par le processeur de l'ordinateur 8, conduisent à l'exécution d'un procédé d'identification selon l'invention. Les différents modules du dispositif d'identification 3 5 comprennent notamment : - un module de modélisation 3A: ce module est configuré pour modéliser le réseau électrique 2 au moyen d'un graphe G comprenant une pluralité de noeuds et une pluralité de branches reliant ces noeuds. Comme détaillé davantage ultérieurement, chaque noeud du graphe 10 correspond à un point de connexion entre deux éléments du réseau : il peut s'agir d'un point de connexion entre deux éléments structurels du réseau ou entre un élément structurel du réseau et une jonction utilisée pour relier cet élément structure à un autre élément structurel du réseau, etc. Chaque branche entre deux noeuds représente donc un 15 élément du réseau, qu'il s'agisse d'un élément structurel (ou d'un segment d'un élément structurel du réseau), d'une jonction assurant une connexion entre deux éléments structurels ou de tout autre élément conducteur du réseau. Le graphe est représenté par une matrice reflétant la topologie du réseau et comprenant des valeurs 20 d'une caractéristique électrique des éléments du réseau représentés par les branches du graphe. Le graphe G est représenté, au niveau du module de modélisation 3A, par une matrice notée dans la suite L reflétant la topologie du réseau 2 (c'est-à-dire ses éléments conducteurs et les interconnexions entre ses différents éléments 25 conducteurs), et comprend les valeurs de caractéristiques électriques des éléments du réseau, à savoir ici les valeurs de leurs admittances ; - un module de détermination 3B: ce module est configuré pour déterminer les valeurs propres et les vecteurs propres de la matrice L; - un module d'estimation 3C : ce module est configuré pour estimer, à 30 partir des vecteurs propres et des valeurs propres de la matrice L, et pour chaque élément d'une pluralité d'éléments du réseau représentés par une branche du graphe, la sensibilité d'une propriété électrique du réseau déterminée entre une source et une charge de l'aéronef à une perturbation de la valeur de la caractéristique électrique comprise dans 35 la matrice L de cet élément. Par sensibilité d'une propriété électrique, on entend ici une variation de cette propriété électrique en réponse 16 aux perturbations appliquées aux caractéristiques électriques des éléments du réseau. Dans l'exemple envisagé ici, la caractéristique électrique perturbée est l'admittance de chaque élément et la propriété électrique en question est une tension entre la source et la charge 5 (désignée dans la suite de la description compte tenu du contexte par « chute » de tension entre la source et la charge). Elle est choisie en liaison avec le critère fonctionnel auquel le système de surveillance s'intéresse pour déterminer l'état du réseau 2 (à savoir ici la chute de tension maximale entre toutes les sources et les charges de l'aéronef). 10 Toutefois, en variante d'autres propriétés peuvent être estimées (par exemple une intensité, une résistance entre la source et la charge, etc.) ; et - un module de détermination 3D: ce module est configuré pour déterminer parmi les éléments du réseau 2, des éléments critiques au 15 regard de la propriété électrique considérée (tension ici) pour la source et la charge auxquelles on s'intéresse. Le module de détermination 3D identifie un élément comme critique lorsque la sensibilité à la perturbation appliquée à la valeur de l'admittance de cet élément dépasse un seuil de criticité prédéterminé noté THR. Un tel élément 20 critique est en effet susceptible de dégrader fortement la connectivité des charges de l'aéronef en cas de dysfonctionnement ou de dégradation. Dans l'exemple envisagé ici, THR = 2% de la chute de tension maximale autorisée (0.5V ici) . Les fonctions des différents modules du dispositif d'identification 25 3 sont décrites plus en détail maintenant en référence aux étapes du procédé d'identification selon l'invention. Comme mentionné précédemment, l'invention offre une solution efficace permettant la surveillance de l'état d'un réseau électrique 2. Cette surveillance vise à vérifier le respect par le réseau électrique 2 d'un critère 30 fonctionnel prédéfini du réseau. Comme indiqué précédemment, dans l'exemple envisagé ici, le critère fonctionnel est considéré comme respecté (et l'état du réseau est jugé normal ou satisfaisant par le système de surveillance 1) si la chute de tension maximale aux bornes de tous les couples source-charge de l'aéronef est inférieure à une chute de tension 35 maximale prise égale à 0.5V lors du fonctionnement normal du réseau 2. 17 La surveillance mise en oeuvre par l'invention se concentre avantageusement sur un nombre limité d'éléments du réseau 2, à savoir sur les éléments identifiés comme critiques par le dispositif d'identification 3 pour chaque couple source-charge de l'aéronef considéré lors de la 5 surveillance. En effet, comme indiqué précédemment, ces éléments critiques sont les éléments du réseau 2 qui sont susceptibles de dégrader le plus fortement la connectivité des charges de l'aéronef en cas de dysfonctionnement ou de dégradation, et donc qui risquent de rendre le réseau électrique 2 de retour de courant non fonctionnel. L'identification 10 de ces éléments critiques permet de simplifier la surveillance du réseau 2. Nous allons décrire, en référence à la figure 2, les principales étapes mises en oeuvre par le dispositif d'identification 3 pour identifier ces éléments critiques du réseau 2. Ces éléments critiques peuvent varier d'une charge à l'autre. Ils sont donc identifiés pour chaque charge de 15 l'aéronef considérée lors de la surveillance du réseau électrique de retour de courant. Pour mieux comprendre le principe de l'invention, il est parfois détaillé dans la suite les hypothèses et les développements mathématiques sur lesquels se sont appuyés les inventeurs pour 20 configurer le dispositif d'identification 3 en vue d'identifier les éléments critiques du réseau 2. Conformément à l'invention, pour identifier les éléments critiques du réseau 2, le dispositif d'identification 3 s'appuie sur une modélisation numérique du réseau électrique. Plus précisément, le 25 dispositif d'identification 3, via son module de modélisation 3A, modélise le réseau électrique 2 au moyen d'un graphe complexe G (étape E10). De façon connue, un graphe est constitué d'un ensemble de noeuds reliés par des branches (ou liens). L'utilisation d'un graphe pour modéliser un réseau maillé tel que le réseau électrique 2 est donc 30 particulièrement bien adaptée. Chaque noeud du graphe complexe G correspond à un point de connexion entre deux éléments conducteurs du réseau électrique 2. Chaque branche du graphe représente un élément conducteur du réseau électrique 2 (élément structurel, segment d'élément structurel ou tout 35 autre élément conducteur faisant partie du réseau électrique 2) . Chaque connexion (jonction) nécessitant au moins une vis est représentée par 18 deux noeuds dans le graphe, ayant les mêmes coordonnées, et une branche de longueur 0 et de poids égal à la valeur de la résistance de contact équivalente. A titre illustratif, pour mieux comprendre comment un tel 5 graphe G est construit, les figures 3A et 3B représentent respectivement une portion de réseau électrique EX1 et le graphe complexe G(EX1) modélisant cette portion de réseau électrique. Plus spécifiquement, la figure 3A représente différents éléments conducteurs du réseau EX1 comprenant notamment deux 10 éléments structurels E1 et E2 et une tresse T. Divers points de connexion existent entre ces éléments conducteurs et sont caractérisés par une jonction J(E1,T) entre l'élément structurel E1 et la tresse T, une jonction 3(E2,T) entre l'élément structurel E2 et T, et une jonction J(E1,E2) entre les deux éléments structurels 31 et J2. On suppose que ces jonctions sont 15 réalisées par exemple au moyen d'une vis. On désigne par segments, les portions des éléments structurels situés de part et d'autre de ces jonctions. Ainsi, chacun des éléments structurels El et E2 est composé respectivement de plusieurs segments, à savoir des segments Bi, B2 et B3 pour l'élément El, et des segments B4, B5 et B6 pour l'élément E2. On 20 note que chaque segment constitue en soi un élément conducteur du réseau au sens de l'invention. La figure 3B représente le graphe G(EX1) modélisant la portion de réseau électrique EX1 illustrée à la figure 3A. Comme mentionné précédemment, chaque connexion nécessitant au moins une vis est 25 représentée par 2 noeuds dans le graphe, ayant les mêmes coordonnées, et une branche de longueur 0 et de poids égal à la valeur de la résistance de contact équivalente. Ainsi, qu'il s'agisse d'une connexion entre un élément de structure et une tresse (comme les jonctions J(E1,T) et J(E2,T) ou tout autre élément conducteur du réseau électrique 2, ou bien 30 une connexion entre deux éléments de structure (comme la jonction J(E1,E2)), les « points » de connexion sont représentés de la même manière. Les branches du graphe correspondent soit à des éléments du réseau, soit à des segments du réseau (comme dans l'exemple illustré sur la figure 3B). Les poids des branches du réseau (B1 à B7 à titre 35 d'exemple) qui ne sont pas des résistances de contact sont égaux aux valeurs des résistances des segments des éléments respectifs. 19 De cette sorte, le graphe complexe G modélise non seulement l'interconnexion des éléments structurels et des autres éléments conducteurs connectés en parallèle sur lesquels s'appuie le réseau mais également les contacts établis au niveau de chaque point de connexion 5 entre ces différents éléments. Un autre exemple est illustré aux figures 4A et 4B. La figure 4B représente un graphe complexe G(EX2) modélisant le réseau électrique résistif EX2 représenté à la figure 4A. Le réseau électrique EX2 comprend 35 résistances identiques R qui sont 10 modélisées par 35 branches B sur le graphe G(EX2) reliant entre eux 25 noeuds. Les 25 noeuds N du graphe G(EX2) représentent les interconnexions entre les résistances R du réseau électrique 2ex. Chaque branche du graphe G(EX2) est caractérisée par une admittance. Le graphe G peut être représenté par une matrice L reflétant la 15 topologie du réseau électrique 2 ainsi que les caractéristiques électriques de ses éléments. Dans l'exemple envisagé ici, comme mentionné précédemment, on s'intéresse comme caractéristiques électriques aux admittances des éléments du réseau 2. Chaque branche du graphe G est donc associée à une admittance qu'il s'agisse d'une admittance associée à 20 un élément structurel du réseau (égale dans ce cas à l'inverse de son impédance, c'est-à-dire dans l'exemple résistif envisagé ici à l'inverse de sa résistance), ou d'une admittance correspondant à une résistance de contact au niveau d'une jonction. La matrice L utilisée par le module de modélisation 3A pour 25 représenter le graphe G est ici la matrice Laplacienne du graphe G modélisant le réseau électrique 2. Les composantes L(u,v) de cette matrice de dimensions nxn où n désigne le nombre de noeuds du graphe G sont définies par, pour u et v désignant des noeuds du graphe : L(u,y = - y', si u et y sont adjacents, 1 yrie Si u = V, tEnel(u) 0 sinon. 30 où yuv désigne l'admittance associée à la branche reliant les noeuds u et y, et nei(u) désigne l'ensemble des noeuds voisins du graphe G directement connectés au noeud u. 20 En variante, une autre matrice représentant le graphe G peut être envisagée, comme par exemple une matrice d'adjacence. Suite à cette modélisation, le module de détermination 3B du dispositif d'identification 3 réalise une analyse spectrale du graphe G et 5 notamment de sa matrice Laplacienne L. Cette analyse spectrale consiste à extraire les valeurs propres et les vecteurs propres de la matrice Laplacienne L (étape E20). Les valeurs propres et les vecteurs propres de la matrice L sont aussi connus sous le nom de spectre de la matrice Laplacienne L. 10 Pour déterminer le spectre de la matrice Laplacienne L, le module de détermination 3B utilise des techniques de diagonalisation de la matrice L connues en soi et non décrites en détail ici. De façon connue, les n valeurs propres 8=(01,021...,en) de la matrice Laplacienne L de dimensions nxn sont les n racines du polynôme 15 caractéristique : det(L -.0 * J) = où 3 désigne la matrice identité, et det le déterminant. La matrice Laplacienne étant symétrique et semi-définie positive, toutes ses valeurs propres sont réelles et positives. On note de plus que la matrice Laplacienne d'un graphe a au moins une valeur propre égale à 0, la 20 multiplicité de cette valeur propre nulle étant égale au nombre de noeuds formant des clusters connectés dans le graphe. Dans la suite de la description, on suppose que les valeurs propres de la matrice L vérifient : O 0 s 02 s ...s O' 25 Les vecteurs propres droits te, et gauches çb, associés à la ième valeur propre ei sont définis par : Lyi = Os!, et OiL =0 ,çb, Seuls les vecteurs propres droits sont utilisés dans la suite de la description. On utilise ici la notation suivante, pour désigner la matrice regroupant les vecteurs propres droits de la matrice L: 30 --[tgi 1/12 - - 1/1' On suppose en outre que les vecteurs propres droits sont orthogonaux entre eux soit : .v =o si ie j, i=1...n, j=1,...,n 35 et qu'ils sont normalisés soit : 21 =1, ce qui implique que IIPIP = . Ces propriétés sont aussi valables pour les vecteurs propres gauches. 5 Le spectre de la matrice Laplacienne L est ensuite utilisé par le module d'estimation 3C du dispositif d'identification 3 pour estimer pour chaque élément d'une pluralité d'éléments du réseau représenté par une branche du graphe G, la sensibilité d'une propriété électrique déterminée entre une source et une charge connectées au réseau électrique, à une 10 perturbation de la valeur de la caractéristique électrique de cet élément du réseau (étape E30). Dans le mode de réalisation décrit ici, le module d'estimation 3C considère plus particulièrement comme propriété électrique du réseau, la tension (ou chute de tension) Vsc entre une source S et une charge C d'un 15 couple source-charge de l'aéronef connecté au réseau électrique 2, lorsque celui-ci est fonctionnel. Pour déterminer les équations appliquées par le module d'estimation 3C pour estimer la sensibilité de la tension Vc, les inventeurs ont conduit le raisonnement suivant. 20 La matrice Laplacienne L représentant le graphe du réseau électrique 2 peut s'écrire sous la forme : L =IP Air où à désigne la matrice diagonale contenant les valeurs propres de L. Par ailleurs, par définition, la matrice Laplacienne L relie les potentiels électriques notés V=(Vi,...,Vn) et les courants nodaux notés 25 I=(Ii,...,In) du graphe G (i.e. considérés en chaque noeud du graphe) via la relation : LV=I Il en découle que : v où A. désigne la matrice pseudo-inverse de Moore-Penrose de A. La matrice A n'étant pas inversible, on considère en outre ici 30 une relation mathématique de régularisation entre les potentiels nodaux \ft, à savoir : Eivi =0 Le potentiel électrique V; d'un noeud i peut donc s'écrire : 22 tek (i)vik (i))1 Eq.1 ok ' En outre, si l'on considère une source de courant électrique S et une charge C de l'aéronef connecté au réseau électrique 2, la résistance entre la source S et la charge C lorsque l'on fait transiter entre eux un 5 courant IO est : -Vc Rsc - OÙ Vs désigne le potentiel de la source S et Vc le potentiel de la charge C. En remplaçant Vs et Vc par la décomposition indiquée à l'équation (Eq. 1) fournie ci-dessus on obtient : I, I=- ( k (C)W k i) j) 0 k k 10 En faisant l'hypothèse que tous les courants nodaux étant nuls exceptés ceux des noeuds S et C et que : /, = -/, =10 on obtient pour la résistance Rsc In ( (s)2 - 2ÇV (S)1//k (c)+vk (C) 2 e ) Ok (tek (s) - tek (c))2 0, 15 Il apparaît donc au regard des équations (Eq. 1) (en remplaçant dans l'équation le noeud i par la source S et par la charge C respectivement) et (Eq. 2) que la tension Vsc = Vs-Vc entre la source S et la charge C et la résistance électrique Rsc entre ces deux noeuds peuvent être déterminées à partir du spectre de la matrice L. Les inventeurs ont 20 donc eu l'idée judicieuse de réaliser une étude perturbative du spectre de la matrice L en vue de prévoir les variations des propriétés électriques du réseau électrique 2 (à savoir notamment ici de la tension Vsc), lorsque des perturbations sont appliquées aux composantes de la matrice Laplacienne L du graphe G modélisant le réseau électrique 2. Ces perturbations sont 25 liées à la modification de la structure du réseau 2 (par exemple en raison d'une dégradation d'un contact entre deux éléments structurels du réseau). On suppose ici que les perturbations de la matrice Laplacienne sont assez faibles afin de rester dans le domaine linéaire. Eq.2 23 Dans la suite de la description, on note LO la matrice Laplacienne du graphe G avant application des perturbations (LO correspond à la matrice L considérée à l'étape E20 décrite précédemment). Toutes les notations relatives à cette matrice introduites 5 précédemment (par exemple pour désigner les valeurs propres et les vecteurs propres de la matrice) sont conservées et indexées par 0. L'équation reliant les ième valeur propre et vecteur propre de la matrice Laplacienne LO est : L0 le = 00, voi Lorsque le réseau 2 est perturbé (par exemple si une impédance 10 de branche est augmentée à cause de la détérioration de la connexion entre les deux noeuds reliés par la branche), la matrice Laplacienne représentant le réseau 2 ainsi que son spectre sont modifiés. On note dans la suite L la matrice Laplacienne caractérisant le réseau électrique 2 après perturbation. 15 La relation entre les ième valeur propre et vecteur propre de la matrice Laplacienne perturbée L est : vf, Oiw, Si la perturbation est faible, il est possible de calculer les valeurs propres et les vecteurs propres de la matrice Laplacienne L du réseau perturbé à partir du spectre de la matrice Laplacienne initiale LO, en 20 linéarisant la relation indiquée ci-dessus. On note AL, Atpi et Mi les variations respectives appliquées à la matrice LO, et aux vecteurs propres droits tpi et aux valeurs propres ei lorsque le réseau 2 est perturbé. Plus particulièrement, il peut être démontré que les valeurs propres Di et les vecteurs propres tpi de la matrice Laplacienne L perturbée 25 peuvent s'écrire sous la forme suivante : 0, = 0,0 + ii0it à"' i,ii, j_ n Oki AL oi e oz Ok - eok Ces équations ont été utilisées judicieusement par les inventeurs pour déterminer les relations appliquées par le module d'estimation 3C du dispositif d'identification 3 pour estimer la sensibilité des potentiels nodaux Vi à une modification de l'admittance d'une branche du graphe. Il s'agit, 30 en calculant cette sensibilité entre la source S et la charge C considérées pour différentes branches du graphe G (ou de manière équivalente, pour 24 différents éléments du réseau 2), de déterminer quelles sont les branches qui, lorsque leur admittance est perturbée, conduisent à une sensibilité accrue (i.e. supérieure à un seuil prédéterminé). Ces branches désignent les éléments critiques du réseau 2 auxquels il convient de porter plus 5 particulièrement attention lors de la détermination de l'état du réseau 2. Plus précisément, lorsque une admittance non nulle Ayuv est ajoutée à l'admittance d'une branche quelconque reliant des noeuds u et v, il en découle une modification AL de la matrice Laplacienne vérifiant : 0 - - - 0 0 - 3Y'i, - AY. AL = 0 0 0 - - AY UV 4145, 0 (les composantes non nulles +/-Ayuy correspondent aux composantes se 10 trouvant à l'intersection de la u-ième ligne et de la v-ième colonne). Les valeurs propres de la matrice Laplacienne LO sont alors modifiées de la façon suivante : 0, = Il résulte de cette relation que la sensibilité (i.e. variation en réponse à la perturbation Ayuv appliquée à la branche reliant les noeuds u 15 et y) des valeurs propres de la matrice Laplacienne à la modification de l'admittance de la branche entre les noeuds u et v est : a ei , De manière similaire, les vecteurs propres du réseau perturbé vérifient la relation suivante : te; teoi .4',, 51, n [( li/O- (U)(gOk (11) - Y/0k (y))+ (Y)(Vok (y) t v0 (u))) I 001- Ook Ooi - Ok La sensibilité des vecteurs propres de la matrice Laplacienne à 20 la modification de l'admittance de la branche entre les noeuds u et v est donc donnée par : a Yfi woi (11)(Wok lgok (4+ Woi(V)(Wok (y) - leok (11))),,, Ok Eq.4 - Ook Ooi - Ook (42 - l'0(u)- w 0, Eq.3 Dans le mode de réalisation décrit ici, le dispositif d'identification 3 utilise les sensibilités des valeurs propres et des vecteurs propres de la 25 matrice Laplacienne du graphe G données respectivement par les équations (Eq.3) et (Eq.4) pour déterminer la sensibilité de la tension entre la source S et la charge C considérées à la variation de l'admittance de chaque branche du graphe. 5 Plus précisément, le potentiel électrique Vo, du noeud i avant perturbations s'écrit : '(i)te k(j) 01, 0 )Ji j-1 -- °Ok où Ij désigne l'intensité nodale au noeud j. Après une petite variation de la matrice Laplacienne, ce potentiel électrique devient : (i)Wk (i) Ok 10 La modification du potentiel électrique résultant d'une petite perturbation appliquée sur une admittance est donc : Il [I1 (;;; t;\ ; ;); ;\ A ;;;; ; ;;; ;\;;; ;\ (Y" Ok VI 1-1Wk le/XY Okl J 1) + kli I Y' Okl' Okl J J)1 AVi = - Onk + A Ok Onk 4- A Ok gOk I- soit, en négligeant les termes du second ordre : A vi [ ( Igk (i)Vok (i) + A j)Wok (0) A00,0* (i)Vok (J) )1 (00/, + Aok ) Onk(Ook +A9k) La sensibilité du potentiel électrique V; du noeud i à la variation d'admittance de la branche du réseau 2 située entre les noeuds u et v 15 vérifie donc l'équation (Eq.5) suivante : If avi _--- ()Yin, j".1 , atek(i),' r 00kvk Ok \JI ' Ok('/ , Okl- YOk \Jok ây', 2 (Ou Ay ', 031 ae A partir de cette équation, la sensibilité de la tension entre la source S et la charge C (autrement dit la sensibilité de la chute de tension mesurée entre ces deux noeuds) à la variation d'admittance d'une branche du réseau peut être aisément calculée par le module d'estimation 3C à 20 partir de la relation : ay, ay',, ay ',, aym, 26 dans laquelle les sensibilités des potentiels électriques Vs et Vc sont évaluées par le second module d'estimation 3C à partir de l'équation (Eq. 5) indiquée ci-dessus appliquée au noeud source S et au noeud charge C, et les sensibilités des vecteurs propres et des valeurs propres de la 5 matrice L sont calculées à partir des équations (Eq.3) et (Eq.4). Le module d'estimation 3C peut être configuré pour calculer la sensibilité de la tension à une variation d'admittance de chaque branche du graphe G. En variante, on peut faire l'hypothèse que seules les jonctions électromécaniques sont susceptibles d'affecter fortement la 10 connectivité du réseau (en raison par exemple d'un boulon desserré, d'une dégradation du contact assurée par la jonction) afin de limiter les calculs : le module d'estimation 3C peut dans ce cas se limiter aux variations des admittances (i.e. des résistances de contact) des branches représentant ces jonctions électromécaniques, les autres admittances étant gardées 15 fixes égales à leurs valeurs initiales. Il convient de noter que la dégradation des autres éléments du réseau correspond à une dégradation structurelle qui peut être détectée par d'autres moyens. A partir des sensibilités estimées par le module d'estimation 3C, le module de détermination 3D du dispositif d'identification 3 est en 20 mesure de déterminer (i.e. d'identifier) les éléments critiques du réseau 2 pour au moins un couple source-charge connecté au réseau (étape E40). Plus spécifiquement, un élément du réseau 2 est identifié comme critique par le module de détermination 3D lorsque la sensibilité de la tension Vsc entre la source S et la charge considérées à la perturbation Ayuv appliquée 25 à la valeur de l'admittance de cet élément critique (c'est-à-dire ici à la branche du graphe représentant cet élément) dépasse un seuil de criticité THR prédéterminé. En relation avec le critère fonctionnel considéré par le système de surveillance, le seuil de criticité THR est pris égal ici à 2% de la chute de tension maximale admissible entre la source électrique et la 30 charge. Ainsi, tous les éléments du réseau pour lesquels la sensibilité de la tension Vsc dépasse le seuil de criticité THR sont considérés comme critiques par le module de détermination 3D. Dans le mode de réalisation décrit ici, les autres éléments du réseau modélisés par une branche sont considérés comme fixes et leurs admittances sont prises égales à leur 35 valeur initiale (valeur considérée dans la matrice Laplacienne LO avant perturbation). 27 La surveillance de l'état du réseau électrique 2 peut donc se concentrer sur les éléments du réseau identifiés comme critiques à l'étape E40, pour chaque source électrique et chaque charge de l'avion connectées au réseau électrique 2. On note que le procédé d'identification 5 fonctionne quel que soit le seuil de criticité THR choisi. Le choix de ce seuil de criticité impacte le compromis entre la précision des résultats obtenue et la rapidité d'exécution du procédé de surveillance. Ainsi, un seuil de criticité pris égal à 0 impose la prise en compte lors de la surveillance de tous les éléments du réseau électrique 2. 10 La figure 5 représente, dans un mode particulier de réalisation, les principales étapes mises en oeuvre lors de cette surveillance par le système de surveillance 1. Ces étapes correspondent aux étapes du procédé de surveillance selon l'invention. La surveillance proposée par l'invention comprend, outre 15 l'identification des éléments critiques du réseau à surveiller, plusieurs étapes qui visent à vérifier la connectivité des charges de l'avion et le respect du critère fonctionnel selon lequel les chutes de tension constatées entre les différentes charges et les différentes sources de l'avion sont inférieures à une chute de tension maximale admissible fixée 20 ici à 0,5 V. Dans le mode de réalisation décrit ici, pour simplifier la surveillance du réseau électrique 2, on considère une étape préliminaire de classement des différentes charges électriques de l'avion connectées au réseau électrique 2 en fonction de leur importance au regard du 25 fonctionnement de l'avion. Cette étape est destinée à déterminer quelles sont les charges connectées au réseau électrique 2 qui sont critiques, c'est-à-dire déterminantes ou encore essentielles pour le bon fonctionnement de l'avion. Ce classement peut être établi à partir d'informations fournies par le systémier ou l'intégrateur aéronautique. Il 30 permet au système de surveillance 1 de sélectionner, parmi la pluralité de charges de l'avion connectées au réseau électrique, un nombre réduit de charges à prendre en compte lors de la surveillance du réseau électrique (étape F10). Plus précisément, seules les charges les plus déterminantes désignées dans la suite par « essentielles » au regard du bon 35 fonctionnement de l'avion sont sélectionnées par le système de surveillance 1. Ceci permet de limiter la complexité de mise en oeuvre de 28 la surveillance, celle-ci n'étant conduite qu'au regard des charges sélectionnées. A partir des charges essentielles sélectionnées, le système de surveillance 1 identifie, par l'intermédiaire du dispositif d'identification 3, 5 pour chaque couple source S - charge C comprenant une charge essentielle, les éléments critiques du réseau électrique 2 à considérer (étape F20). Cette identification est réalisée en appliquant à la source S et à chaque charge C le procédé d'identification selon l'invention, précédemment décrit en référence à la figure 2. Il convient de noter que 10 des éléments du réseau peuvent être critiques pour des charges essentielles distinctes. Le système de surveillance 1, par l'intermédiaire de son module d'estimation paramétrique 6, estime ensuite la valeur de l'admittance de chaque élément critique identifié par le dispositif d'identification 3 pour les 15 différents couples source-charge considérés (étape F30). A cet effet, le module d'estimation paramétrique 6 utilise un modèle numérique MOD modélisant le comportement électrique du réseau électrique 2 en fonction de l'admittance de ses éléments, et des mesures électriques réalisées au moyen des capteurs 4 permettant d'évaluer le 20 comportement du réseau électrique 2 et de le comparer au comportement modélisé par le modèle numérique. Dans le mode de réalisation décrit ici, le modèle numérique MOD considéré par le module d'estimation paramétrique 5 pour représenter le comportement du réseau électrique en fonction des 25 admittances de ses éléments est le même modèle que celui utilisé par le dispositif d'identification, à savoir que dans le réseau 2, les potentiels électriques nodaux V sont reliés aux courants nodaux par l'intermédiaire des admittances des branches du graphe, selon la relation : Y'/=J (MOD) où Y désigne la matrice des admittances des branches du graphe 30 (autrement dit des éléments du réseau), cette matrice étant identique à la matrice Laplacienne. Ce modèle définit les équations circuits du réseau électrique 2. L'estimation paramétrique vise à déterminer les composantes de la matrice Y via une technique itérative. Cette technique est illustrée à la 35 figure 6. Elle consiste de façon générale à exciter le réseau électrique 2 29 au moyen de courants nodaux continus (dans une ou plusieurs configurations), et à mesurer les chutes de tension résultant de l'excitation du réseau électrique au moyen de capteurs dûment positionnés dans le réseau électrique 2. Les chutes de tension correspondant aux courants 5 d'excitation considérés sont par ailleurs estimées au moyen du modèle numérique MOD, et comparer aux chutes de tension mesurées. Les admittances sont estimées en minimisant la différence constatée entre les chutes de tension estimées et mesurées au moyen d'un algorithme d'optimisation. Cette optimisation peut être éventuellement répétée au 10 cours d'une pluralité d'itérations pour permettre l'estimation de toutes les admittances des éléments critiques, les admittances déterminées au cours d'une itération étant considérées comme connues à l'itération suivante. La technique d'estimation paramétrique mise en oeuvre par le module d'estimation 5 s'appuie donc sur des mesures électriques (ex. 15 mesures de potentiels, de chutes de tension et/ou de courants comme détaillé ultérieurement) réalisées au moyen des capteurs 4 au sein du réseau électrique 2. Le placement des capteurs 4 pour réaliser ces mesures est un paramètre important à prendre en compte qui impacte à la fois la complexité du procédé de surveillance et sa précision. 20 Pour déterminer (et optimiser) les emplacements des capteurs 4, les inventeurs ont considéré que plus la chute de tension mesurée est influencée par un élément critique du réseau (autrement dit par une branche du graphe), plus l'information obtenue grâce à la mesure de cette chute de tension est utile pour estimer l'admittance de cet élément 25 critique. Ils ont donc décidé de l'emplacement des capteurs 4 en fonction des tensions (ou des courants) influencées par un ou plusieurs éléments critiques identifiés à l'étape F20. Il convient de noter que le nombre de capteurs 4 à placer dans le réseau électrique 2 résulte également d'un compromis entre complexité 30 et précision de la surveillance mise en oeuvre : en effet, un faible nombre de capteurs fournit un faible nombre de mesures, ce qui se traduit par un gain de temps en maintenance du réseau électrique 2 ; à l'inverse un grand nombre de capteurs fournissant un grand nombre de mesures permet d'améliorer la précision de l'estimation des admittances des 35 éléments critiques. 30 Dans le mode de réalisation décrit ici, il est choisi d'utiliser autant de capteurs que d'éléments critiques identifiés dans le réseau électrique 2 à l'étape F30. Pour chaque admittance à estimer, le chemin le plus influencé entre la source S et la charge C par les perturbations 5 appliquées à cette admittance est identifié dans le réseau électrique 2. Ce chemin peut comprendre un ou plusieurs éléments du réseau (ex. éléments structurels, segments d'éléments structurels et autres éléments conducteurs de connexion), et est donc représenté dans le graphe par une ou plusieurs branches. Il peut être déterminé au moyen d'un procédé 10 similaire à celui utilisé pour identifier les éléments critiques du réseau électrique 2, c'est-à-dire en tenant compte des sensibilités des tensions aux perturbations appliquées aux admittances des éléments du réseau entre les deux noeuds, ces sensibilités étant évaluées à partir de l'équation (Eq. 5) et des équations reflétant la sensibilité des valeurs propres et des 15 vecteurs propres de la matrice Laplacienne L du graphe. Dans ce mode de réalisation, la sensibilité est déterminée sur tous les chemins source S - charges C possibles, mais il est également possible, en considérant que certains éléments du réseau électrique 2 peuvent être inaccessibles ou difficilement accessibles, de déterminer les 20 sensibilités des mesures accessibles et de déterminer les meilleures mesures répondant au besoin. On note que dans l'exemple mentionné ici, on s'attache aux sensibilités des tensions aux extrémités des chemins et des mesures de tensions sont réalisées au moyen des capteurs 4. En variante, il est 25 possible de considérer des mesures de courant de branches ou des mesures de chutes de tension et des mesures de courant de branches. Les capteurs 4 sont placés aux extrémités de chaque chemin identifié pour chaque élément critique du réseau 2. On note en outre que pour obtenir des mesures pertinentes 30 pour l'estimation des admittances des éléments critiques, il est préférable de s'assurer que lors de l'excitation du réseau électrique pour réaliser les mesures par les capteurs 4, les éléments critiques du réseau soient traversés par un courant. On choisira donc des points d'injection des courants d'excitation en conséquence. Plusieurs configurations (par 35 exemple un nombre entier K supérieur à 1) peuvent être envisagées à cet 31 effet. Ce choix peut être réalisé dès la conception du réseau électrique 2 et reste figé. Une fois les capteurs 4 dûment placés dans le réseau électrique, l'estimation paramétrique des admittances des éléments critiques réalisée 5 par le module d'estimation paramétrique 5 est mise en oeuvre de la façon suivante, en s'appuyant sur la représentation par le graphe complexe G du réseau électrique 2. Plus spécifiquement, conformément au modèle numérique MOD utilisé par le module d'estimation paramétrique 5, les potentiels nodaux et 10 les courants nodaux sont liés par le système d'équations (Eq.6) suivant : *v dV (Eq.6) où - M Inc désigne la matrice d'incidence du graphe complexe G. Cette matrice reflète les interconnexions entre les éléments du réseau 15 électrique 2, autrement dit entre les branches du graphe complexe G représentant ces éléments ; - W est la matrice des poids du graphe G; il s'agit d'une matrice diagonale dans laquelle chaque élément représente l'admittance d'une branche du réseau. On note que selon les éléments critiques identifiés 20 lors de l'étape F30, tout ou partie des admittances consignées sont à estimer, les autres étant supposées connues ; - V est un vecteur représentant les potentiels nodaux ; - I est un vecteur représentant les courants nodaux injectés pour la mesure ; 25 - dV désigne un vecteur contenant les chutes de tensions mesurées par les capteurs 4; et - Cmes désigne la matrice dite de combinaison, représentant la matrice d'incidence des mesures ; le nombre de colonnes de cette matrice est égale au nombre de mesures effectuées, tandis que le nombre de 30 lignes est égale à la dimension n du graphe G. La matrice Cmes est construite de manière analogue à la matrice d'incidence sauf qu'elle ne décrit plus une branche entre deux noeuds voisins, mais un chemin (qui est la chute de tension mesurée) entre 2 noeuds non voisins. 32 A titre illustratif, pour mieux comprendre la construction de cette matrice, on considère un exemple basé sur le réseau EX2 illustré à la figure 4A et modélisé par le graphe G(EX2) à la figure 4B (pour un point d'injection du courant d'excitation situé en noeud 11 et une sortie 5 du courant au noeud 15). On suppose qu'on réalise trois mesures de chutes de tension entre les noeuds 6 et 9, entre les noeuds 15 et 24 et entre les noeuds 3 et 9 respectivement. La matrice des combinaisons en résultant est donnée par. 0 0 0 \ 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 -1 0 -1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -1 0 0 0 0 10 Le système d'équations (Eq.6) peut s'écrire de façon équivalente sous la forme : 33 f 01,11 ,L) dV= g(V ,W , I) où f désigne une fonction non linéaire et g désigne une fonction linéaire. On désigne par K le nombre de jeux de mesures expérimentales réalisés sur le réseau électrique 2 au moyen des capteurs 4 au cours d'une même itération : chaque jeu de mesures correspond par exemple à un 5 point d'injection du courant dans le réseau différent, et comprend les mesures réalisées par les différents capteurs 4 associés à chaque élément critique du réseau. L'estimation paramétrique réalisée par le module d'estimation 5 consiste alors à estimer les admittances des éléments critiques 10 correspondant le mieux aux K jeux de mesures effectuées. On désigne dans la suite par West l'estimation de la matrice des poids des éléments du réseau 2 : cette matrice comprend les admittances estimées des éléments critiques et les admittances connues des autres éléments du réseau 2. Le module d'estimation paramétrique 5 estime à 15 partir de l'estimation West des admittances des éléments du réseau 2 et en utilisant le système d'équations (Eq.6), le vecteur dVest des chutes de tension constatées dans le réseau électrique 2 pour chaque jeu de mesures. Ce vecteur vérifie les équations suivantes : dV = g(Ves'W I) Puis il calcule, pour chaque jeu de mesures indexé par k, 20 k=1,...,K, l'erreur d'estimation notée Ek entre les chutes de tensions réelles notées clVr:_eelk mesurées sur le réseau et les chutes de tensions estimées dVestk, soit : ek c1V,'IA -dV''*(V'') La valeur optimale de West est alors obtenue par le module d'estimation paramétrique 5 en minimisant l'erreur quadratique entre les 25 chutes de tensions réelles et estimées sur l'ensemble des K jeux de mesures réalisés, soit : min y ek 2 = min y (d17,,,eik dV''k(V,'))2 0 - f Ksi 'West >1) { 34 Le module d'estimation paramétrique 5 considère par ailleurs, pour la minimisation une contrainte physique imposant que les admittances soient positives. L'estimation paramétrique mise en oeuvre par le module 5 d'estimation paramétrique 5 revient à résoudre un problème d'optimisation dont la formulation est explicitée ci-dessous : min (dV1, - dlles,,(11))2 avec M *Wes, *M' *ve.' -1 wi Pour réaliser cette optimisation, le module d'estimation paramétrique 5 peut avoir recours à différentes techniques d'optimisation. Dans le mode de réalisation décrit ici, compte-tenu du nombre de 10 variables à estimer (chutes de tension et admittances des éléments critiques, les autres étant supposées connues), le module d'estimation paramétrique 5 utilise un algorithme des points intérieurs, connu de l'homme du métier et non détaillé ici. Cet algorithme est décrit par exemple dans l'article de Wâchter et Biegler intitulé «On the 15 implementation of an interior-point ?liter line-search algorithm for large- scale nonlinear programming », 2006. Le module d'estimation paramétrique 5 peut par ailleurs utiliser l'algorithme des points intérieurs en combinaison avec une stratégie permettant d'éviter une convergence de l'algorithme vers des minima 20 locaux. Une telle stratégie peut consister par exemple à appliquer plusieurs fois l'algorithme en modifiant son point de départ (i.e. les valeurs initiale des admittances des branches des contacts critiques). D'autres algorithmes ayant les mêmes propriétés peuvent être utilisés en variante par le module d'estimation paramétrique 5, comme par 25 exemple l'algorithme de Levenberg-Marquardt. Dans un autre mode de réalisation de l'invention, le module d'estimation paramétrique 5 estime les admittances des éléments critiques à partir de mesures de chutes de tensions dVreel et à partir de mesures de courant .reel- Dans ce cas, le problème d'optimisation à résoudre par le 30 module d'estimation paramétrique 5 s'écrit : 35 .1 min (dVre', - dV,',(1/'())2 + y(reeç_ i,', (v,' , y_))2 avec M *W *M *Vest = I IIIC evt 111C - 0 Pour réaliser cette estimation paramétrique, le module d'estimation paramétrique 5 doit non seulement d'estimer les chutes de tension dVest et les poids des branches W mais également les courants nodaux 'est. 5 On note que le module d'estimation paramétrique 5 peut mettre en oeuvre une technique d'estimation paramétrique itérative (en réitérant les étapes décrites ci-dessus), notamment pour garantir que tous les éléments critiques dont l'admittance doit être estimée sont bien parcourus par un courant lors de la réalisation des mesures, et les admittances 10 estimées sont fiables. Chaque itération peut correspondre à un point d'injection distinct du courant sur le réseau 2 pour réaliser les mesures de sorte à permettre l'estimation de toutes les admittances. Les valeurs d'admittance estimées lors d'une itération sont préférentiellement figées et considérées comme connues à l'itération suivante. 15 Par ailleurs, il convient de noter que l'estimation paramétrique qui vient d'être décrite peut également être utilisée pour déterminer les admittances des éléments non critiques du réseau 2 lors d'une étape préliminaire. Ces admittances requièrent une unique détermination puisqu'elles sont ensuite considérées comme invariantes au cours du 20 procédé de surveillance. Pour réaliser cette détermination, il suffit de considérer que toutes les admittances de la matrice L.est doivent être déterminées (il peut dès lors être nécessaire de considérer un plus grand nombre de capteurs pour réaliser des mesures pertinentes pour l'ensemble des éléments du réseau). Des tests conduits par les inventeurs 25 sur des réseaux simples en utilisant cette méthodologie ont montré que des valeurs très proches des valeurs exactes des admittances peuvent être obtenues aisément. L'estimation des admittances des éléments critiques du réseau 2 réalisée par le module d'estimation 5 permet alors au module d'évaluation 30 6 du système de surveillance 1 d'évaluer sur l'ensemble des couples source - charge C considérés, un critère fonctionnel du réseau prédéfini représentatif du fonctionnement du réseau (étape F40). Dans le mode de réalisation décrit ici, le critère fonctionnel évalué est la chute de tension 36 maximale entre chaque charge C et la source S, afin de déterminer si cette chute de tension maximale est conforme à la chute de tension maximale autorisée par la règlementation de 0.5V. Pour évaluer ce critère fonctionnel, le module d'évaluation 6 5 définit un modèle estimé du réseau électrique à partir d'une matrice Laplacienne estimée Lest. Il convient de noter que ce modèle fournit une estimation des paramètres influents du réseau (i.e. des admittances des éléments critiques du réseau) sur les chutes de tensions entre chaque source et 10 chaque charge essentielle mais ne garantit pas les valeurs des admittances des autres éléments du réseau. Un modèle dégradé de la totalité du réseau peut être réalisé en estimant l'ensemble des admittances des éléments du réseau ; une telle estimation requiert toutefois d'installer un grand nombre de capteurs 4 dans le réseau 2 pour 15 disposer des mesures nécessaires à cette estimation. Cette démarche peut être utilisée également pour l'estimation initiale des admittances des éléments du réseau 2 de retour de courant (et plus particulièrement des admittances ou résistances de contact correspondant aux jonctions électromécaniques qui créent des résistances de contact) dès lors que les 20 valeurs de ces admittances ne sont pas connues (les valeurs des admittances des éléments structurels peuvent être obtenues du systémier ou de l'intégrateur aéronautique). La caractérisation initiale des résistances de contact (ou plus généralement des admittances du réseau 2) peut être ainsi mise en oeuvre en reprenant les étapes F10, F20, F30 25 (appliquées à plusieurs reprises) et une étape de validation des admittances estimées en comparant les chutes de tension source-charge estimées via le modèle dérivé des admittances estimées aux chutes de tension source-charge mesurées. Puis, à partir du courant injecté dans le réseau 2 par chaque 30 charge essentielle C considérée lors du fonctionnement normal de l'avion (ce courant est connu des spécifications de l'intégrateur ou peut être aisément mesuré), le module d'évaluation 6 évalue les chutes de tensions entre la source S et chaque charge C essentielle grâce à l'équation fonctionnelle suivante intégrant les équations circuits : 35 L*V=Ifif 37 où V désigne dans cette équation fonctionnelle le vecteur des potentiels nodaux en fonctionnement, et Ifct le vecteur des courants nodaux injectés par les charges essentielles. Puisque la matrice Laplacienne Lest d'un graphe n'est pas 5 inversible, les potentiels nodaux sont calculés en multipliant le vecteur des courants nodaux par la pseudo-inverse de Moore-Penrose de la matrice Laplacienne L. La chute de tensions entre chaque source et chaque charge essentielle peut alors être obtenue à partir des potentiels nodaux estimés (différence des potentiels entre la source S et chaque charge 10 essentielle C) et est comparée par le module d'évaluation 6 au seuil 0.5V. Le système de surveillance 1. est alors en mesure, à partir des critères fonctionnels évalués pour chaque couple source charge considéré, de déterminer l'état de fonctionnement du réseau électrique 2 (étape F50). 15 Si toutes les chutes de tensions constatées entre tous les couples source-charge considérés (c'est-à-dire comprenant les charges « essentielles » au bon fonctionnement de l'avion déterminées à l'étape F10) sont inférieures à 0.5V, alors l'état (de fonctionnement) du réseau électrique 2 de retour de courant est considéré par le module de 20 détermination 7 du système de surveillance 1 comme satisfaisant ou normal. Sinon, le module de détermination 7 détermine que le réseau électrique 2 est dans un état dégradé ou défaillant et nécessite une maintenance. Les éléments dégradés du réseau peuvent être détectés et 25 localisés aisément au vu des différentes quantités évaluées au cours des étapes F20 à F40 et des éléments critiques identifiés, ce qui permet d'envisager une maintenance, une réparation ou un remplacement de ces éléments et d'anticiper un dysfonctionnement plus important du réseau de retour de courant. 30 Dans le mode de réalisation décrit ici, le critère fonctionnel évalué pour déterminer l'état du réseau de retour de courant est la chute de tension maximale entre les sources et les charges de l'avion. Toutefois, d'autres critères fonctionnels peuvent être envisagés, en remplacement ou en plus de celui-ci. Par exemple, on peut envisager également comme 35 critère fonctionnel une valeur maximale de variation constatée entre l'admittance estimée pour un élément critique et une valeur nominale de 38 cette admittance (ou de toute autre caractéristique électrique de l'élément), l'état du réseau électrique étant déterminé comme défaillant si pour au moins un élément critique du réseau pour au moins un couple source-charge considéré, cette variation maximale est supérieure à un 5 seuil d'alerte prédéterminé (par exemple 50%). Les valeurs nominales (ou initiales) des admittances peuvent être connues, ou être déterminées via une estimation paramétrique comme décrit à l'étape F40. L'invention permet ainsi de surveiller aisément l'état de fonctionnement d'un réseau électrique de retour de courant. 10 Nous allons maintenant décrire brièvement à titre illustratif un exemple d'application de l'invention. On considère ici le réseau de retour de courant modélisé à la figure 7. Sur cette figure, le noeud source est représenté par un cercle et les noeuds charges (au nombre de 7) par des carrés. On note que la plupart des noeuds représentés sur la figure 6 qui 15 modélisent un contact entre deux éléments structurels du réseau « cachent » en fait deux noeuds du graphe en raison du contact établi entre les éléments structurels au moyen de jonctions (mais qui n'apparaissent pas en perspective dans le mode de représentation adopté sur cette figure). 20 Dans cet exemple, on fait l'hypothèse que seules les admittances de contact peuvent être dégradées, les admittances des autres éléments du réseau étant considérées comme fixes et connues. De plus les valeurs initiales des admittances de contact sont supposées connues. Il est par ailleurs fait l'hypothèse que les admittances de contact 25 peuvent diminuer de maximum 70% de leur valeur initiale (autrement dit, les résistances de contacts peuvent varier entre leur valeur initiale et 333 % de cette valeur). Dans l'exemple envisagé ici, on considère que le graphe modélisant le réseau de retour de courant illustré à la figure 7 comprend 30 244 noeuds et 284 admittances. Par souci de simplification les noeuds du graphe sont numérotés dans la suite de la description. Comme indiqué précédemment, on considère 7 charges distinctes situées respectivement aux noeuds 108, 109, 112, 114, 119, 127 et 136, et qui injectent un courant identique de 100 A dans le réseau de 35 retour de courant. La source (de puissance est), quant-à-elle, située au noeud 212. Les chutes de tensions entre la source et les charges en 39 fonctionnement normal (i.e. quand le réseau est dans un état satisfaisant) sont données ci-dessous dans la Table 1). On remarque que le réseau de retour de courant est bien dimensionné car le critère fonctionnel est respecté pour toutes les charges (chutes de tensions entre la source et les 5 charges inférieures à 0,5 V). Charge Charge Charge Charge Charge Charge Charge 1 2 3 4 5 6 7 (noeud (noeud (noeud (noeud (noeud (noeud (noeud 108) 109) 112) 114) 119) 127) 136) Chute de tension 0,4014 0,3709 0,3657 0,2768 0,2092 0,2415 0,3955 (V) Table 1: Chutes de tensions entre la source et les charges en fonctionnement normal 10 Dans l'exemple envisagé ici, on suppose qu'un classement préliminaire des charges de sorte à identifier les charges essentielles au bon fonctionnement de l'avion est mis en oeuvre (cf. étape F10 précédemment décrite), et que seules les charges situées aux noeuds 112, 15 127 et 136 sont considérées comme essentielles. Autrement dit, seules les chutes de tensions relatives à ces charges sont surveillées. En application du procédé de surveillance selon l'invention, les éléments critiques du réseau sont identifiés entre chaque couple source-charge considéré. A cet effet, comme détaillé à l'étape F20, les matrices 20 de sensibilité de chaque chute de tension entre la source et une charge aux variations des admittances sont évaluées grâce à l'étude perturbative du spectre du graphe. Les sensibilités sont calculées uniquement pour les admittances des contacts puisque les admittances des autres branches du réseau sont considérées fixées. Dans cet exemple, le seuil de criticité est 25 fixé à -2 % de modification des chutes de tensions. Les sensibilités des chutes de tensions sont calculées et données à la figure 8, pour les différentes charges considérées. On peut remarquer sur cette figure que pour la chute de tension entre la source et la charge 3, il y a 9 branches critiques dont la sensibilité 30 est inférieure à -2 %. Cela signifie que si l'une des admittances de contact correspondant à ces 9 branches critiques diminue de 100 % alors la chute de tension entre la source et la charge augmenterait de plus de 2 °/0. Pour la chute de tension entre la source et la charge 6 on trouve 12 branches critiques et pour celle entre la source et la charge 7 il y en a 8. Cependant, il convient de noter que certaines branches critiques 5 sont communes aux chutes de tensions constatées pour les trois charges. On peut par conséquent limiter à la surveillance à seulement 16 branches critiques pour monitorer les chutes de tensions aux bornes des 3 charges essentielles. Les 16 branches critiques (i.e. éléments critiques du réseau) à surveiller sont les contacts correspondant aux noeuds 11, 12, 13, 26, 28, 10 41, 43, 56, 58, 71, 73, 77, 82, 83, 100 et 101. Dans la suite, les autres admittances de contact sont considérées comme fixes et égales à leur valeur initiale (même si dans le réseau réel, ces admittances peuvent être amenées à varier). Les branches correspondant aux éléments critiques à surveiller 15 (désignées dans la suite par « branches critiques ») sont désormais connues et il convient maintenant de déterminer les mesures électriques qui doivent être réalisées dans le réseau pour fournir les bonnes informations à l'estimation paramétrique. Pour cela les matrices de sensibilités issues de l'étude 20 perturbative du spectre du graphe (étape F20 et E40) sont de nouveau utilisées. On recherche grâce à ces matrices les chemins du réseau les plus influencés par les branches critiques. On suppose ici qu'une mesure de chute de tension est réalisée pour chaque branche critique. La figure 9 donne une représentation graphique de la matrice de sensibilité des 25 chemins du réseau à la variation de l'admittance de la lère branche critique correspondant au noeud 11. On remarque que le chemin entre les noeuds 11 et 148 est le chemin le plus influencé par le contact. On peut noter que ce chemin correspond au chemin entre les bornes du contact 11 car la modification de l'admittance du contact modifie en premier temps la 30 chute de tension de celui-ci. Pour toutes les autres branches critiques les chemins le plus influencés sont aussi les chemins entre les extrémités des branches correspondant aux éléments critiques. Les capteurs sont donc placés aux bornes des contacts critiques afin de réaliser 16 mesures de chutes de tensions. De plus afin de garantir l'excitation des branches 35 critiques, des injections de courant sont réalisées entre leurs bornes. 41 On note qu'il n'est pas toujours possible d'installer des capteurs aux bornes des contacts puisque la structure de l'avion ne permet pas l'accès à tous les éléments du réseau lors des maintenances. Dans ce cas des contraintes techniques peuvent être rajoutées aux matrices de 5 sensibilité afin de ne prendre en compte que les noeuds accessibles (ceci peut influencer la sélection du chemin le plus influencé par chaque branche critique). Une fois les capteurs placés, les 16 mesures de tensions sont effectuées. Dans l'exemple envisagé ici, ces mesures sont réalisées à 10 partir d'une simulation numérique dans laquelle les admittances des contacts sont dégradées, leurs admittances ont été choisies aléatoirement entre la moitié de leur valeur initiale et leur valeur initiale. L'estimation paramétrique décrite à l'étape F30 est ensuite réalisée. Il s'agit de minimiser l'erreur quadratique entre les mesures sur 15 le réseau et les chutes de tensions obtenues par estimation des admittances des branches critiques. Du fait de l'identification des éléments critiques réalisée, les autres admittances de contact sont considérées comme égales à leur valeur initiale dans l'optimisation. Le problème d'optimisation est résolu 15 fois avec différents points initiaux afin d'éviter 20 la convergence vers un minimum local. L'erreur quadratique obtenue entre les 16 chutes de tensions mesurées sur le réseau réel et celles estimées est de 3,08*10-7. Cette solution est en moyenne obtenue après 95 itérations en 1 min et 15 secondes. Les estimations des admittances des branches critiques ainsi 25 que les erreurs avec les admittances réelles sont données dans la Table 2 reproduite ci-dessous : Erreur 42 ontact 11 Contact 12 Contact 13 Contact 26 Contact 28 Contact 41 Contact 43 Contact 56 Contact 58 Contact 71 Contact 73 Contact 77 Contact 82 Contact 83 Contact 100 Contact 101 b4S,2 640,9 0,66 922,6 983,7 -6,21 317,2 315,7 0,46 453,4 463,3 -2,15 614,9 623,3 -1,34 614,2 611,7 0,41 526,2 522,9 0,64 509,0 508,2 0,16 640,3 639,5 0,13 628,1 625,2 0,47 642,9 640,6 0,36 370,2 370,2 0,00 835,1 823,3 1,44 974,3 990,0 -1,58 52 1,6 622,1 -0,09 890,8 890,2 0,07 Table 2 Comparaison entre les admittances estimées réelles des 16 branche.s critiques. Il apparait dans cette table que l'erreur quadratique est très faible. Les admittances estimées sont pertinentes pour reconstituer les chutes de tensions mesurées sur le système dégradé réel. On remarque en outre que les emplacements des capteurs sont bien définis puisque la minimisation de l'erreur quadratique permet d'obtenir une bonne estimation des admittances du réseau. On note cependant que l'estimation 10 du noeud 12 est légèrement moins bonne avec une erreur de -6,21 %. Ceci peut s'expliquer par l'incertitude existant sur les admittances des éléments (et plus particulièrement des contacts) non critiques du réseau. Pour améliorer cette estimation, un autre seuil de critieité peut être considéré afin de prendre en compte un nombre d'éléments critiques plus 15 importants. L'estimation paramétrique a donc permis d'estimer correctement les admittances des éléments critiques du réseau. Afin de vérifier le critère fonctionnel de chute de tensions entre la source et les charges essentielles, le modèle estimé du réseau dégradé 20 est construit à partir des acim ttances de contact critiques estimées lors de l'estimation paramétriques. Les admittances des contacts non critiques et des autres éléments sont égales à leur valeur initiale. 43 Les potentiels de l'ensemble des noeuds du réseau de retour de courant lors de son fonctionnement normal sont estimés au moyen du modèle fonctionnel : La Table 3 reproduite ci-dessous reporte les chutes de tensions estimées entre la source et les 3 charges essentielles Table 3 : Comparaison entre chutes de tensions estimées et réelles du réseau dégradé 10 D'après les estimations des admittances des branches critiques, le réseau de retour de courant est fonctionnel pour les 3 charges essentiel;es. Cependant les chutes de tensions des charges 3 et 7 sont proches de la limite de 0,5 V. Les chutes de tensions réelles permettent la 15 même conclusion. Pour les charges essentielles, l'erreur d'estimation maximale enregistrée est de 1,1%, ce qui demeure extrêmement précis. Cet exemple montre bien que la surveillance proposée par l'invention permet d'obtenir une bonne. approximation de l'évolution des chutes de tensions entre la source et les charges essentielles lors de la 20 dégradation du réseau de retour de courant. La précision de l'estimation peut être améliorée en diminuant le seuil de criticité et donc en estimant davantage de branches critiques. htite de tension ée ( 0,3957 0,4001 -1,1 0,2684 0,2704 -0,7 0,5931 0,5969 -0,6 REVENDICATIONS1. Procédé d'identification d'éléments critiques d'un réseau électrique (2) maillé de retour de courant d'un aéronef, ledit aéronef comprenant au moins une structure en matériau composite, ledit réseau électrique comprenant une pluralité d'éléments conducteurs, ledit procédé d'identification comprenant : - une étape de modélisation (E10) du réseau électrique au moyen d'un graphe comprenant une pluralité de noeuds et une pluralité de branches reliant ces noeuds, chaque noeud du graphe correspondant à un point de connexion entre deux éléments du réseau, chaque branche entre deux noeuds représentant un élément du réseau, ledit graphe étant représenté par une matrice reflétant la topologie du réseau et comprenant des valeurs d'une caractéristique électrique des éléments du réseau représentés par les branches du graphe ; - une étape de détermination (E20) des valeurs propres et des vecteurs propres de la matrice représentant le graphe ; - une étape d'estimation (E30), à partir des vecteurs propres et des valeurs propres de la matrice, et pour chaque élément d'une pluralité d'éléments du réseau représenté par une branche du graphe, d'une sensibilité d'une propriété électrique déterminée entre une source et une charge connectées au réseau électrique à une perturbation de la valeur de la caractéristique électrique de cet élément ; et - une étape de détermination (E40), parmi les éléments du réseau, d'éléments critiques au regard d'un état de fonctionnement du réseau, un élément du réseau étant déterminé comme critique lorsque la sensibilité de la propriété électrique estimée pour cet élément dépasse un seuil de criticité prédéterminé. 2. Procédé d'identification selon la revendication 1 dans lequel : - la caractéristique électrique est une admittance de l'élément du réseau ; - la matrice du graphe est une matrice Laplacienne ; et - la propriété électrique est une tension électrique entre la source et la charge. 45 3. Procédé de surveillance d'un réseau électrique maillé de retour de courant d'un aéronef, ledit aéronef comprenant au moins une structure en matériau composite, ledit réseau électrique comprenant une pluralité d'éléments conducteurs, ledit procédé de surveillance 5 comprenant : - une étape d'identification (F20) d'au moins un élément critique du réseau, pour au moins un couple source-charge comprenant une source et une charge de l'aéronef connectées au réseau, en utilisant un procédé d'identification selon la revendication 1 ou 2; 10 - une étape d'estimation (F30) paramétrique d'une valeur de la caractéristique électrique de chaque élément critique, cette étape utilisant un modèle numérique modélisant un comportement électrique du réseau en fonction de ladite caractéristique électrique et des mesures électriques réalisées au moyen de capteurs positionnés dans 15 le réseau et permettant d'évaluer ledit comportement électrique, un capteur étant positionné dans le réseau pour chaque élément critique sur au moins un élément du réseau ou sur au moins un chemin du réseau formé d'une pluralité d'éléments du réseau influencé par une perturbation de la caractéristique électrique de cet élément critique ; 20 - une étape d'évaluation (F40), pour ledit au moins un couple source-charge, d'un critère fonctionnel du réseau prédéfini, cette étape d'évaluation utilisant la valeur de la caractéristique électrique estimée pour chaque élément critique ; et - une étape de détermination (F50) d'un état du réseau électrique en 25 tenant compte du critère fonctionnel évalué pour ledit au moins un couple source-charge. 4. Procédé de surveillance selon la revendication 3 dans lequel le critère fonctionnel évalué est une valeur maximale d'une chute 30 de tension constatée entre chaque source et chaque charge, l'état du réseau électrique étant déterminé comme satisfaisant si cette chute de tension est, pour chaque couple source-charge, inférieure à ladite valeur maximale. 35 5. Procédé de surveillance selon la revendication 3 dans lequel le critère fonctionnel évalué est une variation constatée entre la 46 caractéristique électrique estimée pour un élément critique et une valeur nominale de cette caractéristique électrique, l'état du réseau électrique étant déterminé comme défaillant si pour au moins un élément critique du réseau pour au moins un couple source-charge, cette variation est 5 supérieure à un seuil d'alerte prédéterminé. 6. Procédé de surveillance selon l'une quelconque des revendications 3 à 5 comprenant en outre une étape préliminaire (F10) de sélection, parmi une pluralité de charges de l'aéronef connectées au 10 réseau électrique, d'au moins une charge considérée comme essentielle au regard du fonctionnement de l'aéronef, les étapes d'identification, d'estimation paramétrique, et d'évaluation étant appliquées uniquement aux couples source-charge comprenant ladite au moins une charge essentielle sélectionnée. 15 7. Procédé de surveillance selon l'une quelconque des revendications 3 à 6 dans lequel lors de l'étape d'estimation paramétrique, on utilise autant de capteurs que d'éléments critiques identifiés, chaque capteur étant associé à un élément critique et placé sur l'élément ou sur le 20 chemin du réseau le plus influencé par une perturbation de la valeur de caractéristique électrique de cet élément critique. 8. Procédé de surveillance selon l'une quelconque des revendications 3 à 7 dans lequel l'étape d'estimation paramétrique met en 25 oeuvre une technique d'optimisation dite des points intérieurs sur une pluralité d'itérations. 9. Procédé de surveillance selon l'une quelconque des revendications 3 à 8 comprenant une étape préalable de caractérisation 30 initiale des caractéristiques électriques à partir du graphe modélisant le réseau, de la matrice représentant le graphe et de mesures électriques réalisées dans le réseau. 10. Dispositif (3) d'identification de points de connexion 35 critiques dans un réseau électrique maillé de retour de courant d'un aéronef, ledit aéronef comprenant au moins une structure en matériau 47 composite, ledit réseau électrique comprenant une pluralité d'éléments conducteurs, ledit dispositif d'identification comprenant : un module de modélisation (3A), configuré pour modéliser le réseau électrique au moyen d'un graphe comprenant une pluralité de noeuds 5 et une pluralité de branches reliant ces noeuds, chaque noeud du graphe correspondant à un point de connexion entre deux éléments du réseau, chaque branche entre deux noeuds représentant un élément du réseau, ledit graphe étant représenté par une matrice reflétant la topologie du réseau et comprenant des valeurs d'une caractéristique 10 électrique des éléments du réseau représentés par les branches du graphe ; un module de détermination (3B), configuré pour déterminer des valeurs propres et des vecteurs propres de la matrice représentant le graphe ; 15 un module d'estimation (3C), configuré pour estimer à partir des vecteurs propres et des valeurs propres de la matrice, et pour chaque élément du réseau représenté par une branche du graphe, d'une sensibilité d'une propriété électrique déterminée entre une source et une charge connectées au réseau électrique à une perturbation de la 20 valeur de la caractéristique électrique de cet élément ; et un module de détermination (3D), configuré pour déterminer parmi les éléments du réseau, des éléments critiques au regard d'un état de fonctionnement du réseau, un élément étant déterminé comme critique lorsque la sensibilité de la propriété électrique estimée pour cet 25 élément dépasse un seuil de criticité prédéterminé. 11. Système de surveillance (1) d'un réseau électrique maillé de retour de courant d'un aéronef, ledit aéronef comprenant au moins une structure en matériau composite, ledit réseau électrique comprenant une 30 pluralité d'éléments conducteurs, ledit système de surveillance comprenant : un dispositif d'identification (3) selon la revendication 10, apte à identifier pour au moins un couple source-charge comprenant une source et une charge de l'aéronef connectées au réseau, au moins un 35 élément critique au regard du fonctionnement du réseau parmi les éléments du réseau ; 48 - une pluralité de capteurs (4) positionnés dans le réseau et aptes à réaliser des mesures électriques ; - un module d'estimation paramétrique (5), configuré pour estimer une valeur de la caractéristique électrique de chaque élément critique, ce 5 module utilisant un modèle numérique modélisant un comportement électrique du réseau en fonction de ladite caractéristique électrique et des mesures électriques réalisées au moyen desdits capteurs et permettant d'évaluer le comportement électrique, un capteur étant positionné dans le réseau pour chaque élément critique sur au moins 10 un élément du réseau ou sur au moins un chemin du réseau formé d'une pluralité d'éléments du réseau influencé par une perturbation de la caractéristique électrique de cet élément critique ; - un module d'évaluation (6), configuré pour évaluer pour ledit au moins un couple source-charge, un critère fonctionnel du réseau prédéfini, ce 15 module utilisant la valeur de la caractéristique électrique estimée pour chaque élément critique ; et - un module de détermination (7), configuré pour déterminer un état de fonctionnement du réseau électrique en tenant compte du critère fonctionnel évalué pour ledit au moins un couple source-charge 20 considéré.