FR 2931553 A1 20091127 FR 0853350 A 20080522 La présente invention se rapporte à un procédé et un système d'estimation d'un critère de charge subie par un composant structurel d'un aéronef, en particulier un avion de transport civil, à un procédé et un système correspondant d'aide à la détection d'un atterrissage dur dudit aéronef, ainsi qu'à un aéronef comprenant de tels systèmes. Par charge, on entend, dans le cadre de la présente invention, tout effort que supporte la structure de l'aéronef, et en particulier les composants analysés. Les "charges", dans leur sens générique, regroupent toute force, moment, facteur de charge (ou accélération), contraintes structurelles. Dans un but évident de sécurité, une inspection de la structure de l'aéronef, ainsi que la réparation d'éventuelles détériorations, doivent être effectuées dès lors que l'aéronef est soumis à des charges importantes, ce qui est notamment le cas lors d'un atterrissage dur (ou "hard landing" selon la terminologie anglo-saxonne largement utilisée dans le domaine). Un atterrissage dur se caractérise par un facteur de charge vertical (ou accélération) et une vitesse verticale de l'aéronef importants lors de l'atterrissage. L'inspection évoquée ci-dessus est généralement déclenchée par le pilote, en fonction de ce qu'il a ressenti au moment de l'atterrissage, lorsqu'il juge que cet atterrissage a pu provoquer des dommages réels de la structure de l'aéronef. Cependant, les réactions perçues au niveau du poste de pilotage traduisent parfois mal le niveau de charge réel appliqué sur l'ensemble de l'aéronef, et en particulier sur chaque composant que l'on souhaite surveiller. Des systèmes automatisant la mesure de charges de dimensionnement d'un aéronef ont alors été proposés. Ainsi, la solution du document EP-A1-0 657 723 décrit un système comparant une charge mesurée à un emplacement, à une valeur de charge admissible déterminée à partir d'une loi d'évolution prédéterminée et d'un ensemble de paramètres mesurés. Le résultat de la comparaison permet alors d'indiquer si une inspection doit être effectuée, et éventuellement les parties de l'aéronef devant être inspectées. Néanmoins, cette solution souffre de nombreux inconvénients. D'une part, elle requiert l'incorporation, dans l'aéronef et en sus des capteurs classiquement utilisés, de capteurs spécifiques pour mesurer les charges subies aux divers emplacements d'intérêt. De tels capteurs peuvent être nombreux selon le nombre d'éléments à surveiller. Ils constituent, en outre, un coût financier conséquent pour équiper toute une flotte d'avions, un coût de maintenance des capteurs (pour remédier à d'éventuelles dérives), ainsi qu'un coût en poids supplémentaire, pour chacun des avions, proportionnel au nombre d'emplacements à analyser. D'autre part, la réalisation parfois millimétrique des aéronefs rend l'intégration de tels capteurs très complexe sans requérir des études onéreuses et approfondies, voire impossible dans les régions inaccessibles. Il en résulte que l'inspection physique de ces parties non équipées par un opérateur doit toujours être maintenue après chaque atterrissage. Cette situation est alors préjudiciable en termes de coûts et de temps de maintenance de l'aéronef. La présente invention entend remédier aux inconvénients de l'art antérieur en proposant un procédé et un système associé permettant d'estimer plus simplement les charges subies en de nombreux emplacements de la structure de l'aéronef. A cet effet, l'invention vise notamment un procédé d'estimation d'un critère de charge subie par un composant structurel d'un aéronef, comprenant une étape de mesures de paramètres dudit aéronef, caractérisé en ce qu'il comprend une étape de calcul d'au moins un critère de charge dudit composant structurel à l'aide d'au moins un réseau de neurones recevant en entrée des valeurs déterminées en fonction desdits paramètres mesurés. De façon avantageuse, les réseaux de neurones permettent de représenter simplement des modèles très complexes. Ces modèles sont ainsi créés et utilisés à moindres coûts, notamment de traitements, ce qui permet de réaliser le procédé selon l'invention à bord de l'aéronef. On utilise ainsi ces réseaux pour représenter les modèles, généralement non connus, liant les valeurs mesurées et différents emplacements structurels de l'aéronef. On comprend ici que le réseau de neurones correspond à toute représentation ou structure mathématique pour l'exploitation de poids et seuils qui auront été déterminés au préalable, notamment par apprentissage du réseau à partir d'au moins un jeu d'apprentissage associant des valeurs des paramètres mesurés à des valeurs de critère de charge. La puissance de modélisation des réseaux de neurones permet également de se satisfaire des capteurs classiques équipant les avions, sans en installer de nouveaux, afin de déterminer les efforts en tout point de la structure de l'aéronef. De tels capteurs sont, par exemple, ceux prévus pour l'enregistrement des mesures d'un enregistreur numérique de données de vol ou DFDR ("Digital Flight Data Recorder" selon la terminologie anglo-saxonne). Classiquement, il s'agit de façon non exhaustive des accélérations et vitesses du corps rigide de l'avion, de différents angles de représentation de l'avion dans l'espace et des vitesses associées à ces angles, de mesures dites statiques telles que la position du centre de gravité de l'avion, le carburant à bord et le poids brut de l'avion. Ainsi, l'on n'a pas besoin de prévoir un capteur de mesure de charge pour chaque composant à surveiller. En outre, sans l'ajout de nouveaux capteurs, il est aisé d'intégrer la présente invention aux avions déjà en service. Dans un mode de réalisation, on prévoit une étape préalable de formation dudit réseau de neurones, cette étape préalable comprenant une première étape d'exécution d'une pluralité de simulations en faisant varier au moins un, généralement un ensemble, desdits paramètres mesurés, la simulation déterminant une valeur simulée d'une charge subie par ledit composant, et une deuxième étape d'apprentissage dudit réseau de neurones à l'aide desdites simulations, ledit réseau de neurones recevant en entrées ledit au moins un paramètre. Grâce à la simulation, généralement assistée par ordinateur, il est possible d'obtenir des mesures cohérentes pour des composants inaccessibles dans un avion et notamment en tout point de l'avion. Corrélativement, l'invention concerne un système d'estimation d'un critère de charge subie par un composant structurel d'un aéronef, comprenant des moyens de mesures de paramètres dudit aéronef, caractérisé en ce qu'il comprend au moins un réseau de neurones agencé pour recevoir en entrée des valeurs déterminées en fonction desdits paramètres mesurés et pour fournir en sortie au moins un critère de charge dudit composant structurel. Un tel système selon l'invention peut comprendre des moyens se rapportant aux caractéristiques du procédé présenté ci-dessus. L'invention concerne également un procédé d'aide à la détection d'un 10 atterrissage dur d'un aéronef, caractérisé en ce qu'il comprend: - la détermination d'un instant d'impact dudit aéronef sur une piste d'atterrissage à partir desdits paramètres mesurés; - l'estimation d'une pluralité desdits paramètres audit instant d'impact déterminé; 15 - l'estimation d'au moins un critère de charge d'un composant structurel dudit aéronef selon le procédé vu ci-dessus, dans lequel le réseau de neurones reçoit en entrée lesdits paramètres estimés. Ainsi, lors de l'atterrissage, des critères de charge en divers emplacements de la structure d'aéronef sont calculés au moment de l'impact de 20 cet aéronef avec la piste, c'est-à-dire lorsque les contraintes d'atterrissage sur la structure sont les plus importantes. En pratique, on opère de la sorte pour plusieurs composants structurels, c'est-à-dire en utilisant plusieurs réseaux de neurones configurés en conséquences, recevant chacun tout ou partie des paramètres estimés. Ces 25 composants peuvent par exemple être la voilure, le train d'atterrissage, le fuselage, les mâts de réacteur. Grâce à la présente invention, on dispose d'informations sur les charges subies par un ensemble de composants structurels de l'aéronef. Il devient alors aisé, par exemple par simple comparaison, d'évaluer si ces 30 composants nécessitent une inspection. Ces critères de charge calculés permettent également au pilote de conforter ou infirmer son analyse par ressenti, pour fournir un rapport de vol. De façon optionnelle, le procédé peut comprendre des étapes se rapportant aux caractéristiques des procédé et système d'estimation d'un critère de charge tels que présentés ci-dessus. Dans un mode de réalisation de l'invention, ledit au moins un critère de charge calculé est comparé à au moins une première valeur seuil de sorte à déterminer une opération d'inspection d'un premier type à opérer sur ledit composant structurel, suite audit atterrissage. Cette première valeur seuil peut représenter une première capacité structurale du composant, par exemple une limite entre l'utilisation du composant dans le cadre d'une déformation élastique et l'utilisation du composant dans le cadre d'une déformation plastique. Il est ainsi possible de prévoir plusieurs valeurs seuil délimitant autant de types d'opérations d'inspection à prévoir. En particulier, si ledit au moins un critère de charge calculé est supérieure à ladite première valeur seuil, on compare ledit critère de charge calculé à une deuxième valeur seuil de sorte à déterminer une opération d'inspection d'un deuxième type à opérer sur ledit composant structurel. Cette deuxième valeur seuil peut représenter une autre capacité structurale du composant à partir de laquelle des dommages affectant la sécurité de l'aéronef peuvent être créés. Des mesures de sécurité telles que le blocage de l'aéronef au sol jusqu'à analyse de celui-ci par l'avionneur peuvent être alors prévues. Ainsi, grâce à ces dispositions, on peut obtenir rapidement et automatiquement une indication des composants structurels à inspecter selon le degré d'inspection à effectuer. La maintenance et l'entretien de l'avion s'en trouvent facilités. On peut également prévoir qu'une pluralité de critères de charge d'une pluralité respective de composants structurels est calculée et que le procédé comprend un étape d'indication audit pilote, par exemple par l'affichage d'une information selon laquelle l'atterrissage a été dur, lorsque la comparaison d'au moins un critère de charge est telle qu'une opération d'inspection doit être effectuée. Cela permet au pilote de confirmer ou non son analyse par ressenti. Dans un mode de réalisation de l'invention, ladite détermination d'un instant d'impact comprend le calcul de l'instant correspondant au premier impact de l'aéronef sur la piste, ledit calcul de l'instant étant réalisé à l'aide de paramètres mesurés représentatifs de la vitesse d'au moins une roue dudit aéronef. Au contraire des solutions connues basant la détermination d'impact à l'aide de mesures d'enfoncement d'amortisseurs des trains d'atterrissage, ce mode de réalisation permet d'obtenir une détermination plus précise de l'instant d'impact compatible avec le calcul d'un critère de charge. Notamment, ce calcul peut comprendre la comparaison de la vitesse de l'au moins une roue dudit aéronef avec une valeur seuil. Par exemple, elle comprend une régression linéaire à partir de deux mesures de vitesse de ladite au moins une roue. En particulier, ladite détermination comprend, en outre, le calcul d'au moins un autre instant d'impact, postérieur audit instant de premier impact, ledit calcul de l'autre instant étant réalisé à partir d'un paramètre mesuré d'accélération verticale dudit aéronef. Alors que la détermination de l'impact à l'aide de la vitesse des roues n'est plus possible dès le premier impact réalisé, cette caractéristique permet d'assurer par un autre biais une détermination de l'instant d'impact également plus précise que dans les solutions connues et compatible avec la précision requise pour le calcul d'un critère de charge. Selon une caractéristique particulière, on calcule un incrément d'accélération verticale aux voisinages d'au moins deux instants représentatifs d'au moins deux impacts et on sélectionne l'instant présentant le plus grand incrément. On simplifie ainsi les calculs ultérieurs en se concentrant sur l'impact ayant le plus grand mouvement vertical de l'aéronef, c'est-à-dire celui susceptible d'engendrer le plus de dommages. De façon classique, les capteurs utilisés dans les avions effectuent des mesures discrètes à des vitesses d'échantillonnage définies. Ainsi, on prévoit que ladite pluralité de paramètres est estimée par interpolation linéaire audit instant d'impact desdits paramètres mesurés. On accroît ainsi la précision du calcul final du critère de charge. Selon une caractéristique de l'invention, on envisage que lesdits paramètres sont mesurés et enregistrés à compter de l'instant où ledit aéronef passe sous une altitude prédéfinie par rapport à ladite piste d'atterrissage lors d'une phase d'approche. En combinaison, on prévoit de cesser ces mesures et enregistrements dès l'instant où la vitesse de l'aéronef, roulant sur ladite piste, passe sous une vitesse sol seuil. En variante, les mesures et enregistrements cessent après une durée prédéterminée, par exemple 30 secondes. On s'assure ainsi de disposer des mesures pour aussi bien le premier impact de l'aéronef sur la piste que pour les éventuels rebonds ultérieurs de l'aéronef. Corrélativement, l'invention a également trait à un système d'aide à la détection d'un atterrissage dur d'un aéronef, caractérisé en ce qu'il 10 comprend: - un moyen de détermination d'un instant d'impact dudit aéronef sur une piste d'atterrissage à partir desdits paramètres mesurés; - un moyen d'estimation d'une pluralité desdits paramètres audit instant d'impact déterminé; 15 - un système d'estimation d'au moins un critère de charge d'un composant structurel dudit aéronef selon le système ci-dessus, dans lequel le réseau de neurones reçoit en entrée lesdits paramètres estimés. De façon optionnelle, le système peut comprendre des moyens se rapportant aux caractéristiques du procédé d'aide à la détection d'un 20 atterrissage dur ci-dessus. L'invention vise également un aéronef comprenant un système d'estimation d'un critère de charge ou d'aide à la détection d'un atterrissage dur tel que présenté ci-dessus. De façon optionnelle, l'aéronef peut comprendre des moyens se 25 rapportant aux caractéristiques de systèmes présentées ci-dessus. D'autres particularités et avantages de l'invention apparaîtront encore dans la description ci-après, illustrée par les dessins ci-joints, dans lesquels : - la figure 1 représente une vue générale d'un tracé d'aéronef lors de 30 son atterrissage sur piste; - la figure 2 représente schématiquement un exemple d'architecture de système selon la présente invention, embarqué dans l'aéronef de la figure 1; - la figure 3 illustre la structure générale d'un réseau de neurones; - la figure 4 représente, sous forme d'un logigramme, des étapes d'un exemple de traitement selon l'invention; - la figure 5 représente, sous forme d'un logigramme, des étapes du traitement des paramètres mesurés lors du processus de la figure 4; - la figure 6 illustre un exemple d'organisation de trains d'atterrissage illustrant une étape du processus de la figure 5; - la figure 7 illustre le calcul d'un instant de premier toucher da l'avion de la figure 1 sur la piste; - la figure 8 illustre le calcul d'un instant de toucher ultérieur de l'avion sur la piste correspondant à un rebond de celui-ci; - la figure 9 illustre l'étape de calcul d'un critère de charge mis en oeuvre lors du processus de la figure 4; et - la figure 10 illustre l'étape de calcul d'un critère de charge mis en oeuvre lors du processus de la figure 4. Sur la figure 1, on a représenté un avion 10 lors de son atterrissage sur une piste 12 d'aéroport. Sur cet exemple, la trajectoire de l'avion 10 montre qu'il subit un rebond et présente donc deux impacts 14 et 16 avec le sol de la piste avant de s'immobiliser en fin de piste. L'avion 10 est en communication avec une station 18, ici prévue à l'aéroport, pour échanger des données. L'avion 10 embarque à son bord un système 20 pour la mise en oeuvre de l'invention, un exemple de ce système étant plus détaillé sur la figure 2. La finalité de ce système 20 est de fournir un rapport renseignant de la sévérité de l'atterrissage, éventuellement détaillé en précisant différents composants structurels potentiellement altérés par l'atterrissage. En référence à la figure 2, ce système 20 pour déterminer la sévérité de l'atterrissage comprend: - une pluralité de capteurs 22; répartis à divers emplacements de l'avion 10; - un système informatique 24 comprenant un unité de traitement 26, une mémoire 28 mémorisant un fichier 30 de paramètres mesurés dans l'avion 10 et plusieurs définitions de réseaux de neurones 32, ainsi qu'une interface 33 avec des moyens d'entrées et de sorties, type écran et clavier, prévus dans le cockpit de l'avion 10; et - un moyen de communication 34, type radio ou système ACARS ("Aircraft Communications Addressing and Reporting System" selon la terminologie anglo-saxonne), permettant de communiquer des données avec la station 18. Les capteurs 22; sont des capteurs classiquement utilisés dans les avions pour diverses mesures de paramètres dits "paramètres avion". Ceux-ci correspondent aux paramètres enregistrés de manière standard pour l'enregistreur numérique de données de vol DFDR équipant les avions. Ici, on prévoit la mesure de paramètres d'accélération et de vitesse du corps rigide de l'avion (par exemple la vitesse au sol Vx(t)), la mesure de différents angles permettant la représentation de l'avion dans l'espace (par exemple l'angle de roulis cp(t)), la mesure des vitesses associées à ces mêmes angles, la mesure de paramètres considérés comme "statiques" au regard de l'atterrissage tels que la position du centre de gravité de l'avion 10, le carburant à bord et le poids brut de l'avion. D'autres paramètres sont classiquement mesurés et leurs capteurs associés peuvent également être prévus sans affecter l'invention, par exemple la position des volets, de la gouverne, ... On utilise des capteurs présentant une fréquence d'échantillonnage définie (At étant la période d'échantillonnage) pour mesurer les paramètres variant dynamiquement. Les paramètres relativement stable tout au long de l'atterrissage ne sont mesurés qu'une seule fois lors du début d'enregistrement comme expliqué plus tard; c'est notamment le cas de l'estimation du carburant à bord, du poids de l'avion ainsi que de celui du centre de gravité. Le système informatique 24 appartient à un système global de surveillance des conditions de l'avion, autrement appelé ACMS ("Aircraft Condition Monitoring System" selon la terminologie anglo-saxonne). Le fichier 30 est un fichier numérique stockant les paramètres mesurés par les capteurs 22; au fur et à mesure de leurs acquisitions, ces paramètres étant éventuellement horodatés pour une exploitation ultérieure plus aisée. Ce fichier 30 peut, en fin d'atterrissage, être transmis, via les moyens radio 34, à la station au sol 18 pour effectuer une analyse similaire ou plus poussée que celle réalisée à bord de l'avion et décrite ci-après dans un exemple. Chaque définition de réseau de neurones 32 correspond à une modélisation mathématique pour le calcul d'un charge subie par un composant structurel de l'aéronef, type une voilure 100, une portion de fuselage 102 ou un mât moteur 104. On dispose d'autant de définitions 32 qu'il y a de composants structurels identifiés à surveiller lors de l'atterrissage. Ces définitions 32 comprennent des coefficients (poids et seuils) permettant de créer, de manière logicielle dans l'unité de traitement 26, des réseaux de neurones, ici de type réseau de neurones linéaire entièrement connecté à plusieurs couches intermédiaires, ici deux, et à rétroaction ("twohidden layer linear feed-forward fully-connected neural network"). Les réseaux de neurones correspondent à des structures mathématiques mises en oeuvre généralement par une programmation logicielle, utilisant les poids et seuils déterminés lors d'un apprentissage préalable de la structure. Une illustration partielle d'un tel réseau 36 est proposée en figure 3, disposant d'une couche d'entrée 38 recevant des entrées x;, des deux couches intermédiaires 40 définies par les poids a;,; et bi,k et par les seuils ao,i et bo,k, d'une couche de sortie 42 définie par ses coefficients ck et fournissant une valeur de sortie s (les valeurs seuils sont ici unitaires pour chacun des neurones des couches d'entrées 38 et de sortie 42). Un tel réseau 36 est muni d'une fonction d'activation f(x), par exemple 1/(1+x) à chaque neurone. Les couches intermédiaires 40 sont notamment pourvues chacune de 20 à 30 neurones, par exemple 25 neurones. Comme illustré par la figure 4, le présent exemple comprend généralement une première phase d'acquisition des paramètres de l'avion 50 en relevant notamment les mesures des capteurs 22;, une deuxième phase de traitement de ces paramètres acquis 52, puis le calcul 54 d'un critère de charge pour un composant structurel (100, 102, 104) de l'avion 10, notamment à l'aide du réseau de neurones 36 associé à ce composant structurel (tandis que d'autres paramètres peuvent être traités directement comme on le verra par la suite), et enfin la détermination 56 de la sévérité de l'atterrissage au moins au regard du critère de charge calculé. En sortie de traitement, on obtient une indication 58, par exemple unique, qui est affichée dans le cockpit de l'avion 10 pour aider le pilote à confirmer son ressenti de l'atterrissage aux fins d'établir un rapport sur la sévérité de l'atterrissage. Cette indication 58 peut également prendre la forme d'un rapport informatique automatiquement transmis par exemple aux équipes de maintenance. En variante, cette indication 58 peut préciser la situation de sévérité au regard de chacun des composants analysés, comme évoquée par la suite. Plus en détail, la première phase 50 d'acquisition comprend l'enregistrement des mesures effectuées dans le fichier 30. Cette acquisition est déclenchée lors que l'avion 10 est en approche et passe sous un plafond prédéterminé 60, comme illustré par la figure 1, ici l'altitude de 100 pieds (soit 300 m environ). Les acquisitions et enregistrements des mesures sont poursuivis jusqu'à ce que l'avion 10 ait définitivement atterri. Pour ce faire, on cesse ces opérations après une période de 30 secondes jugée suffisante. En variante, on peut prévoir de cesser ces opérations dès lors que la vitesse au sol de l'avion 10 passe sous une valeur prédéterminée 62, ici 20 noeuds (environ 37 km.h-l), permettant d'avoir entre 30 et 60 secondes d'enregistrement de paramètres pour le traitement ultérieur 19. Pour la suite de la description, on utilise l'unité de "noeud" largement diffusée dans le domaine aéronautique. L'homme du métier saura aisément convertir ces mesures en km.h-' par l'équivalent 1 noeud = 1,85 km.h-l. Le traitement 52 des données ainsi acquises est déclenché, quant à lui, lorsque la vitesse au sol de l'avion 10 passe également sous une valeur prédéterminée, par exemple la même valeur que celle mettant fin aux acquisitions. En variante, on peut retenir une valeur plus faible permettant de commencer le traitement ultérieurement, par exemple à l'arrêt de l'avion au parking. Le traitement 52 vise à obtenir les paramètres pertinents pour l'évaluation des critères de charge des composants surveillés au moment de l'impact le plus fort. Comme illustré par la figure 5, ce traitement 52 comprend une première étape 70 de prétraitement au cours de laquelle on vérifie que tous les paramètres mesurés attendus pour la suite de l'analyse sont présents. Puis on vérifie que la fréquence d'échantillonnage correspond à celle attendue. Enfin, on s'assure que les valeurs acquises des paramètres satisfont des plages de valeurs prévues pour ces paramètres. A défaut un message d'erreur peut être remonté au système et au pilote. A l'étape 72, on détermine l'instant to de premier impact 14 de l'avion 10 sur la piste 12. En effet, plusieurs impacts successifs peuvent avoir lieu, par exemple lorsque l'avion 10 rebondit sur la piste 12. D'une façon générale, on utilise la vitesse des roues des trains d'atterrissage pour déterminer cet instant. On illustre cette étape à l'aide de la figure 6 représentant les trains d'atterrissage d'un avion A380. 22 roues sont réparties sur cinq trains. Seules les roues R1 à R16 de la figure sont munies de tachymètres pour mesurer leurs vitesses respectives. On récupère donc les vitesses linéaires cl) des roues R1 à R16 enregistrées dans le fichier 30. On détermine l'instant t; où les vitesses mesurées dépassent une valeur seuil A1, ici 15 noeuds. En notant At la période d'échantillonnage des mesures sur les roues, ici 1/16 s, on calcule pour chacune des roues R1 à R16, l'instant t tel que ((t;) Al, ce qui représente sensiblement l'instant du toucher de la roue i sur la piste 12. A partir de cette valeur t;, on calcule trois valeurs tom, toi2, toi3 par interpolation, telles qu'illustrées sur la figure 7. Puis on choisit comme instant toi de toucher de la roue i sur la piste, la valeur maximale obtenue: to; = max { tom, toi2, toi3}. Le toucher to; le plus précoce pour l'ensemble des roues correspond à l'instant to de premier impact de l'avion 10 sur la piste 12: to = min { toi}iE[1;16]. On a ainsi interpolé la valeur to de l'instant d'impact à partir des mesures enregistrées. A l'étape 74, on calcule les paramètres correspondant au moment de l'impact 14 déterminé. Ici, on s'intéresse aux paramètres dynamiques qui nécessitent une acquisition dans le temps, au contraire de paramètres peu variables tels que le poids de l'avion, la carburant à bord ou le centre de gravité acquis une fois pour toute au début de l'enregistrement 60. On interpole les paramètres dynamiques à l'instant to à partir des mesures enregistrées, ici par interpolation linéaire entre les deux échantillons mémorisés au voisinage de to. On obtient ainsi des accélérations et vitesses interpolées du corps rigide de l'avion (par exemple les accélérations latérales et verticales à to: nyo et nzo), des angles interpolés de représentation dans l'espace de l'avion, des vitesses angulaires interpolées correspondantes, tous à l'instant to du toucher.. On détermine ensuite la plus grande accélération verticale nzomax subie sensiblement au moment et postérieurement de l'impact, comme suit: nzomax = max { nz(t) } lorsque te [to; to + 1 seconde]. On déduit alors l'incrément de facteur de charge verticale subie: Anzo = (nzomax - nzo) De façon similaire, on détermine l'incrément de facteur de charge latérale subie: Anyo = (I ny(t) - nyo I) lorsque te [to; to + 1 seconde]. L'ensemble des paramètres correspondant à l'instant to d'impact peut être mémorisé sous forme d'un vecteur d'une dimension. Un vecteur K1 peut être appliqué pour multiplier ce vecteur de paramètres afin d'appliquer, le cas échéant, un coefficient correcteur. C'est notamment le cas lorsque les essais pour élaborer les réseaux de neurones 36 sont réalisés dans des conditions légèrement différentes que lors de l'exploitation de l'avion, par exemple avec des fréquences d'échantillonnage différentes. Par défaut, ce vecteur est composé uniquement de 1. On calcule également lors de cette étape 74, un coefficient de dérotation ed-derot correspondant à la valeur maximale d'abaissement du nez de l'avion 10 dans la deuxième phase de l'atterrissage (6d(t) étant négatif lors de l'abaissement du nez), c'est-à-dire après le toucher initiale 14 à to. Ce calcul est effectué à l'aide de la vitesse angulaire de tangage 6d(t) comme suit: on détermine t2 tel que 6d(t2) On détermine si une telle fenêtre existe dans laquelle: max(nz(t)) û min(nz(t)) > A3, par exemple A3 = 0,7g; auquel cas un rebond a été identifié. On tient compte du rebond uniquement si cet incrément max(nz(t)) û min(nz(t)) > Anzo. En effet, ici on cherche à déterminer l'instant où l'accélération verticale subie par l'avion 10 est la plus grande. On évacue alors les rebonds présentant un incrément inférieur à celui du premier impact to sur la piste 12. Si plusieurs rebonds ont été identifiés, on ne conserve que celui présentant un incrément d'accélération verticale Anz le plus important. Dans le cas où un rebond 16 est identifié et conservé (sortie "oui" de l'étape 76 sur la figure 5), on exécute l'étape 78 consistant à déterminer les valeurs des paramètres à l'instant de ce rebond. A cette étape, on calcule l'accélération verticale minimale nz1min obtenue dans cette fenêtre temporelle identifiée correspond à l'instant t1 d'impact de l'avion 10 sur la piste 12 lors du rebond 16 identifié. On calcule également l'accélération verticale maximale nzimax dans cette fenêtre (l'instant associé est noté timax) de sorte à définir l'incrément à cet impact de rebond: Anz1 = (nzimax - nz1min)• A l'étape 78 toujours, on calcule de façon similaire à l'étape 74 les paramètres à l'instant t1 de l'impact de rebond. On calcule également l'incrément de charge latérale Any1 = (1 ny(t) ù ny1 1) lorsque te [ti; t1 + 1 seconde], où 1 x l est la valeur absolue de x. Un facteur K2 similaire à K1 est éventuellement utilisé pour appliquer des coefficients correcteurs à ces paramètres calculés. On poursuit à l'étape 80. Si aucun rebond n'est identifié ou conservé (sortie "non" de l'étape 76), on passe directement à l'étape 80. Pour la suite de la description, on considère qu'un rebond 16 a été identifié et conservé. Les calculs ultérieurs sont alors basés en partie sur les valeurs indexées par "1", c'est-à-dire correspondantes à l'impact t1 de rebond 16. On comprend aisément qu'en l'absence de rebond pertinent, ces valeurs indexées par "1" sont à remplacer par celles correspondantes au premier impact 14 de l'avion sur la piste, c'est-à-dire indexées par "0". L'étape 80 a pour objectif de distinguer les calculs à effectuer en fonction de premiers critères représentatifs de la sévérité de l'impact de l'avion sur la piste. Deux critères rentrent en ligne de compte: d'une part le poids brut Pavion de l'avion, et l'incrément d'accélération verticale Anz1 subie par l'avion. On qualifie l'atterrissage comme faiblement sévère lorsque: - Anz1 A5, par exemple A6 = 0,5g. Dans ce cas, on considère que seuls les trains d'atterrissage peuvent avoir été endommagés. Le traitement est alors plus simple que pour un atterrissage jugé plus sévère comme décrit plus loin. Ainsi, ici, on détermine à l'aide d'une table prédéfinie et selon les différents anglés interpolés à l'instant t1, les trains d'atterrissage à inspecter suite à l'atterrissage. Dans les cas d'atterrissage faiblement sévère non couverts par les cellules de cette table, aucune inspection de composants de l'avion n'est à prévoir; l'atterrissage n'est pas qualifié de "dur". Cette information déterminée peut être remontée au pilote. Le traitement se termine alors. On s'intéresse maintenant plus précisément aux cas où l'on considère l'atterrissage comme potentiellement sévère, c'est-à-dire lorsque: Anz1 > A4 et Pavion A6 et Pavion > A5. On poursuit alors à l'étape 82 par le calcul d'un coefficient d'adhérence Sr utilisé par la suite comme entrée de certains réseaux de neurones dans l'évaluation des critères de charge. Ce coefficient Sr caractérise une contribution des efforts dus aux roues. Cette étape 82 comprend tout d'abord la détermination du temps tsu de mise en rotation des roues, calculé par exemple comme la valeur minimale du temps de mise en rotation de chaque roue: tsu = min {tsu_i} i E [1; 16]. Chaque temps tsu i de mise en rotation d'une roue Ri est déterminé comme le temps entre toi précédemment calculé lors de l'étape 72 ci-dessus et un instant tsi interpolé au voisinage où la vitesse angulaire coi dépasse une valeur seuil. L'étape 82 comprend ensuite la détermination du coefficient d'adhérence ou de frottement Sr. A cette fin, on utilise un premier réseau de neurones qui prend, en entrées, certains des paramètres avion précédemment introduits ainsi que les valeurs calculées Anzo et tsu. La configuration du réseau de neurones est généralement effectuée, de façon préalable à son intégration dans l'avion, par l'estimation des poids et seuils adaptés à définir un modèle mathématique liant les paramètres ci-dessus à la valeur Sr souhaitée en sortie. Ce processus d'élaboration comprend par exemple un ensemble de simulations, par exemple plusieurs milliers de simulations d'atterrissage, dans lesquelles ont fait varier les paramètres d'entrées ci-dessus. Plusieurs essais d'apprentissage peuvent être réalisés en jouant sur le nombre de neurones par couche et/ou le nombre de couches. On détermine ensuite, avec les jeux de valeurs pour valider le modèle, un écart type représentatif de l'erreur engendrée par le modèle. On conserve alors la configuration de réseau de neurones présentant l'écart type minimal. On peut donc obtenir des architectures de réseaux de neurones différentes selon les modèles mathématiques et les simulations réalisées. En appliquant les valeurs d'entrée au réseau de neurones (aux positions d'entrées respectives prévues lors de l'apprentissage) prévu pour le calcul du coefficient d'adhérence, on obtient une valeur Sr que l'on compare à deux valeurs seuils: 0,8 et A8=0,95. Si 0,8 Sr A8, alors Sr = 0,8; permettant de limiter la valeur de l'adhérence à celles classiquement utilisées Si Sr Ainsi, à l'instant t1, on détermine la vitesse minimale des roues que l'on compare à une valeur seuil A9, par exemple 40 noeuds, à partir de laquelle on estime que l'adhérence n'engendre plus d'effort: Si min {c (ti)} ; E [1; 16] > A9, alors Sr = 0 car les roues sont déjà en rotation Sinon Sr est inchangé. Cette limitation par rapport à A9 est optionnelle et peut être omise pour diminuer les traitements à effectuer. A la fin de l'étape 82, on dispose de l'ensemble des paramètres calculés ou interpolés (tous issus des paramètres avion enregistrés) pour procéder au calcul à proprement parler de critères de charge sur différents composants structurels de l'avion 10. On rappelle ici que l'on est dans le cas où un atterrissage potentiellement sévère a été détecté. On passe alors à l'étape 54. En référence à la figure 9, cette étape 54 comprend trois ensembles de calcul de critères de charge ou de sévérité de l'impact, calculs qui peuvent être réalisés successivement ou en parallèle: charge sur le train d'atterrissage du nez (90), charge sur les autres trains d'atterrissage à l'aide de l'accélération latérale (92), et charges sur tout type de composant structurel à surveiller, incluant les trains d'atterrissage si souhaité (94). Lors de l'étape 90, on calcule à l'aide d'un réseau de neurones un premier paramètre Cmod, caractérisant la vitesse de l'avion à l'impact considéré. Ce réseau de neurones reçoit en entrée les valeurs calculées ou interpolées précédemment. Puis, on procède aux calculs suivants de sorte à générer un critère de charge du train avant en fonction du poids de l'avion Pavion, de Cmod, de l'angle de dérotation Od derot. En variante, cette deuxième sous-étape peut être intégrée au réseau de neurones utilisé pour le calcul de Cmod sans difficulté particulière. On obtient ainsi un premier critère de charge du train d'atterrissage de nez à l'aide d'un réseau de neurones. Lors de l'étape 92, on détermine, de façon similaire à l'étape 80, à l'aide d'une table, si une intervention sur un ou plusieurs trains d'atterrissage est nécessaire. L'information est remontée directement au pilote qui la consigne dans son rapport. Lors de l'étape 94, on détermine le critère de charge pour chacun des composants structurels à surveiller. N'importe quel composant structurel à l'intérieur de l'avion peut être concerné, un réseau de neurones correspondant étant alors prévu comme indiqué précédemment. Ainsi, on injecte en entrée xi du réseau de neurones correspondant à un composant à surveiller, un ensemble de paramètres calculés précédemment. Ici, on n'utilise, de façon avantageuse, que des paramètres "avion" issus des capteurs classiques et estimées au moment du toucher concerné par les traitements décrits ci-dessus. Ce réseau ayant été configuré comme expliqué ci-dessus, il fournit une valeur Ci de critère de charge subie par ledit composant structurel. A titre d'exemple, on liste ci-après quelques exemples de composants dont il est utile de déterminer les critères de charge subie: - charge verticale des trains d'atterrissage des ailes ou du fuselage; - charge de bielle de recul des trains d'atterrissage du fuselage; - moment de torsion des ailes; - moment de flexion du fuselage; - moment de torsion du fuselage; et - charge verticale du fuselage. Les critères de charge ainsi obtenus en sortie des différents réseaux de neurones ne sont pas corrélés entre eux, au contraire des paramètres d'entrées qui correspondent généralement à l'instant de l'impact considéré. Une fois ces critères de charge obtenus, on passe à l'étape 96 déterminant la sévérité de l'atterrissage sur chacun des composants surveillés. Comme illustré par la figure 10, lors de cette étape, on compare le critère de charge C; obtenu à deux valeurs seuils prédéfinies Sil et Si2 (spécifiques au composant considéré) comme caractéristiques des capacités structurales du composant. Cette comparaison peut notamment être mise en oeuvre à l'aide d'un ratio entre le critère calculé et par exemple la première valeur seuil Sil . Si le critère de charge obtenu est inférieur aux deux valeurs seuils, aucun message n'est émis concernant le composant associé. Si le critère de charge obtenu est supérieur au premier des deux seuils, un message d'un premier type, ici qu'une inspection du composant doit être effectuée, est émis. Cela correspond à un composant ayant potentiellement subi une altération détectable. Dans ce cas, on vérifie également si le critère de charge est supérieur au deuxième seuil. Dans l'affirmative, le composant doit être considéré comme potentiellement sérieusement endommagé. Un message indiquant qu'une inspection et des vérifications détaillées doivent être effectuées est émis. Généralement, une réparation comme indiquée dans les manuels de réparation structurelle prévus est effectuée. Il peut également être indiqué de contacter le constructeur de l'avion. Les traitements précédents à base de réseaux de neurones ne nécessitent pas de puissances de calcul considérables. Ainsi, les messages peuvent être émis alors même que l'avion arrive à son emplacement de parking. Ainsi, le pilote reçoit l'ensemble des messages dont il tient compte dans son rapport de vol. Ces messages lui permettent notamment de confirmer ou infirmer son ressenti de l'atterrissage pour qualifier ou non celui-ci de "dur". Ces messages peuvent être également sauvegardés ou envoyés, vis les moyens radio 34, à la station 18. Cette dernière dispose ainsi de façon détaillée de l'ensemble des composants à inspecter. Ainsi, au contraire des nombreuses solutions de l'art antérieur où tous les composants devaient être inspectés en cas d'atterrissage dur, ici l'inspection et la maintenance sont guidées et largement réduites. La présente invention peut ainsi s'appliquer en vue de la détermination d'un indicateur de sévérité d'atterrissage aussi bien pour les atterrissages classiques, que ceux avec rebond(s) ou avec dérapage du à une accélération latérale. Les exemples qui précèdent ne sont que des modes de réalisation de l'invention qui ne s'y limite pas. REVENDICATIONS1. Procédé d'estimation d'un critère de charge (C;) subie par un composant structurel (100, 102, 104) d'un aéronef (10), comprenant une étape de mesures (50) de paramètres dudit aéronef, caractérisé en ce qu'il comprend une étape de calcul (54) d'au moins un critère de charge (C;) dudit composant structurel à l'aide d'au moins un réseau de neurones (36) recevant en entrée (x;) des valeurs déterminées en fonction desdits paramètres mesurés. 2. Procédé d'aide à la détection d'un atterrissage dur d'un aéronef (10), caractérisé en ce qu'il comprend: - la détermination (72, 76) d'un instant d'impact (14, 16) dudit aéronef sur une piste (12) d'atterrissage à partir desdits paramètres mesurés; - l'estimation (74, 78) d'une pluralité desdits paramètres audit instant d'impact déterminé (to, ti); - l'estimation (54) d'au moins un critère de charge d'un composant structurel dudit aéronef selon le procédé de la revendication 1, dans lequel le réseau de neurones (36) reçoit en entrée lesdits paramètres estimés. 3. Procédé selon la revendication précédente, dans lequel ledit au moins un critère de charge (C;) calculé est comparé (96) à au moins une première valeur seuil (Sil, Si2) de sorte à déterminer une opération d'inspection d'un premier type à opérer sur ledit composant structurel, suite audit atterrissage. 4. Procédé selon la revendication 3, dans lequel une pluralité de critères de charge d'une pluralité respective de composants structurels est calculée et le procédé comprend un étape d'indication (58) au pilote de l'aéronef lorsque la comparaison (96) d'au moins un critère de charge est telle qu'une opération d'inspection doit être effectuée. 5. Procédé selon l'une des revendications 2 à 4, dans lequel ladite détermination d'un instant d'impact comprend le calcul (72) de l'instant (to) correspondant au premier impact de l'aéronef sur la piste à l'aide de paramètres mesurés représentatifs de la vitesse (ci)) d'au moins une roue (R;) dudit aéronef,et le calcul d'au moins un autre instant d'impact (t1), postérieur audit instant de premier impact à partir d'un paramètre mesuré d'accélération verticale (nz) dudit aéronef (10). 6. Procédé selon la revendication précédente, dans lequel on calcule un incrément d'accélération verticale (Anzo, Onze) aux voisinages d'au moins deux instants (to, ti) représentatifs d'au moins deux impacts (14, 16) et on sélectionne l'instant présentant le plus grand incrément. 7. Procédé selon l'une des revendications 2 à 6, dans lequel lesdits paramètres sont mesurés et enregistrés à compter de l'instant où ledit aéronef passe sous une altitude prédéfinie par rapport à ladite piste d'atterrissage lors d'une phase d'approche. 8. Système d'estimation d'un critère de charge (C;) subie par un composant structurel d'un aéronef, comprenant des moyens de mesures (22;) de paramètres dudit aéronef, caractérisé en ce qu'il comprend au moins un réseau de neurones (36) agencé pour recevoir en entrée des valeurs déterminées en fonction desdits paramètres mesurés et pour fournir en sortie au moins un critère de charge dudit composant structurel. 9. Système (20) d'aide à la détection d'un atterrissage dur d'un aéronef, caractérisé en ce qu'il comprend: - un moyen de détermination d'un instant d'impact dudit aéronef sur une piste d'atterrissage à partir desdits paramètres mesurés; - un moyen d'estimation d'une pluralité desdits paramètres audit instant d'impact déterminé; - un système d'estimation d'au moins un critère de charge d'un composant structurel dudit aéronef selon le système de la revendication 8, dans lequel le réseau de neurones reçoit en entrée lesdits paramètres estimés. 10. Aéronef (10) comprenant un système (20) d'estimation d'un critère de charge ou d'aide à la détection d'un atterrissage dur selon l'une des revendications 8 et 9.