FR 2615641 A1 19881125 FR 8707048 A 19870520 La présente invention concerne un procédé d'élaboration d'un modèle statistique pour prédéterminer la charge de travail d'un pilote d'aéronef, notamment d'un avion de type commercial, au moyen des paramètres objectifs de vol. L'invention vise aussi, & titre de produit industriel, le modèle statistique élaboré par le procédé considéré. L'invention vise également diverses applications de ce modèle, notamment & la certification de la composition de l'équipage d'aéronefs, au développement de nouveaux aéronefs et & la formation, au perfectionnement, au contr8le des performances et au suivi du degré d'adaptation à un nouvel aéronef de pilotes d'aéronefs. L'invention concerne encore un dispositif pour la mise en oeuvre du procédé susvisé. La charge de travail imposée aux pilotes des avions commerciaux est une source de préoccupation majeure pour les constructeurs d'avions et les compagnies aériennes pour au moins deux raisons: il est généralement reconnu que des niveaux élevés de charge de travail entraînent une augmentation de la probabilité d'erreur du pilote et une diminution de sa capacité à faire- face aux situations d'urgence. D'autre part, la possibilité de démontrer que les charges de travail d'un équipage ne dépassent pas un niveau acceptable est une partie essentielle du processus pour obtenir une certification de type par les autorités compétentes. L'accord n'est pas unanime en ce qui concerne la définition du terme "charge de travail". Même si la plupart des chercheurs font une différence entre la charge de travail "subjective" et la charge de travail "effective", les écrits traitant de ce sujet contiennent encore de nombreuses références vagues & la "charge de travail" pour lesquelles il n'est pas parfaitement clair de déterminer s'il s'agit de la charge mentale subjective ou de la charge de travail effectif. De multiples recherches sur la charge de travail se sont concentrées sur les pilotes et l'environnement relativement complexe qui existe dans un poste de pilotage et au cours de vols dans des zones de trafic aérien de densité variable élevée à faible. Cependant, on a fourni peu de données sur des pilotes commerciaux en vol réel. L'étude approfondie de l'Armée de l'Air américaine dans l'Air Force Aerospace Medical Research Laboratory (AFAMRL) s'est limitée aux avions et systèmes militaires (Shingledecker 1983). Lorsque des pilotes civils ont été utilisés, les expériences ont eu lieu sur simulateurs ou installations de laboratoire analogues. La recherche sur la charge de travail s'est intéressée particulièrement à deux principaux domaines. Le premier implique l'élaboration d'une théorie de charge de travail ou d'un modèle qui caractériserait avec précision le fonctionnement du système humain de traitement de l'information et décrirait la nature des capacités de ce système qui limite les performances (Shingledecker, 1983). Aucun modèle réalisé n'a encore fait l'unanimité. Certains chercheurs mettent l'accent sur l'environnement extérieur comme source essentielle de charge de travail. D'autres donnent comme cause déterminante majeure de la charge de travail l'aptitude du pilote et tentent de fixer la réserve de capacité comme une mesure de charge de travail. D'autres encore se fondent sur la nature subjective de la charge de travail, suivant la théorie selon laquelle les pilotes doivent "se sentir" stressés pour l'ftre. Enfin, un groupe de chercheurs a envisagé le concept de charge de travail comme une construction de l'esprit déterminée par des facteurs multiples et n'ayant aucune existence propre. Dans cette optique, les impératifs de tâche, les opérations mentales et les impressions subjectives se combinent toutes avec une multitude de variables d'état comme la formation, le niveau de fatigue, la motivation et la condition physique pour produire dans l'individu une capacité réelle à réaliser 261564-1 des performances (O'Donnel and Enggemeier, 1985). Le second domaine majeur dans lequel s'est effectuée la recherche concerne la mesure de la charge de travail, de quelque manière qu'elle soit définie, et la mise au point de techniques de mesures de comportement physiologiques et subjectives. Un grand nombre d'études a été consacré & la mise au point de trois types de mesures de charge de travail. Cependant aucune proposition significative n'a été faite en ce qui concerne une mesure de la charge de travail des pilotes de transport aérien commercial au cours de vols "typiques". Les techniques connues d'étude du comportement pour mesurer la charge de travail consistent le plus souvent à donner au pilote une tache secondaire à exécuter et à mesurer le niveau de rendement pour cette tache si le scénario de vol se modifie. On suppose que le niveau de rendement en ce qui concerne la tache secondaire se détériorera en fonction de la charge de travail impliquée par la tâche principale de vol. Baser les estimations de charge de travail sur le comportement m me des pilotes présente plusieurs avantages théoriques. D'abord ces mesures sont parfaitement objectives et il est donc très peu probable qu'elles soient influencées par des facteurs extérieurs comme la motivation ou la fausse interprétation du pilote. Ensuite, il serait rationnel de les utiliser puisqu'elles sont issues de performances humaines. En dépit de leurs possibilités, il a été constaté, dans le cadre de l-a présente invention, que les techniques d'étude du comportement entraînent plusieurs problèmes majeurs, en particulier si l'on envisage de les utiliser dans des conditions réelles de vol. D'abord, l'introduction d'une tache secondaire nécessite souvent une instrumentation supplémentaire. Ensuite, elle peut gner l'exécution de la tâche principale du pilote et diminuer ainsi la sécurité. En fait, une mesure valable de la charge de travail par l'étude du comportement doit empiéter sur la tâche principale. Enfin, un certain degré d'entraînement à la tâche d'étude du comportement doit exister pour que cette étude soit utilisée d'une manière fiable. A cause de ces restrictions, la recherche sur des mesures de charge de travail par l'étude du comportement a été limitée essentiellement au laboratoire ou aux simulateurs. Différentes mesures physiologiques ont été également proposées et essayées comme mesures de charge de travail. Elles vont de la mesure de la variabilité du rythme cardiaque (arythmie) jusqu'à des essais neurophysiologiques extrêmement complexes comme ceux étudiés par 1'AFAMRL (O'Donnell, 1983). Ces mesures physiologiques requièrent toutes des équipements spéciaux qui ne pourraient pas, d'une façon pratique, &tre montés dans un poste de pilotage et utilisés pendant un vol réel. On a également examiné le rapport entre le rythme cardiaque et la charge de travail. et des mesures de variabilité de rythme cardiaque ont fait l-objet d'un travail considérable (Mulder and Mulder- Hajonides Van Der Meulen. 1973). On s'accorde généralement à reconnaitre que la variabilité du rythme cardiaque, correctement mesurée, a de grandes possibilités, en tant qu'indice d'une certaine composante de la charge de travail d'un pilote. La présente invention a précisément mis à profit le fait que le rythme cardiaque pouvait se mesurer en, vol avec facilité et discrétion à l'aide d'un matériel portatif relativement bon marché. On a également étudié l'application des mesures subjectives de charge de travail des pilotes dans le cadre à la fois de la recherche et de la certification. Cooper et Harper (1969), à titre d'exemple, ont ainsi mis au point une échelle utilisable pour évaluer les qualités de vol d'un avion. Il est relativement facile d'appliquer des mesures subjectives sans pratiquement aucune incidence sur les opérations de vol. Ces mesures ont donc été utilisées en particulier au cours de processus de certification d'avions pour démontrer que la cotation des pilotes sur la charge de travail imposée par un nouvel avion reste dans des limites acceptables. Cependant m me si les mesures subjectives sont faciles à appliquer. à enregistrer et à analyser, elles représentent encore une charge supplémentaire dans le processus de certification et sont sujettes à des déformations possibles. Si les valeurs subjectives données par les pilotes au cours des essais de certification pouvaient Etre prévues d'une façon précise en utilisant un modèle approprié, le besoin de recueil de cotation serait supprimé. L'une des idées à la base de la présente invention est la suivante: comme les mesures subjectives réelles intègrent la plupart des facteurs déterminants d'environnement, psychologiques et physiologiques de la charge de travail. tout modèle d'évaluation subjective de charge de travail devrait bénéficier de l'introduction des paramètres physiologiques, de comportement et d'environnement. Un but de la présente invention est de proposer un procédé d'élaboration d'un modèle que l'on puisse mettre en oeuvre pratiquement tant en vol réel que Simulé pour déterminer de manière aussi précise que possible la charge de travail subjective des pilotes d'aéronef, à partir d'un nombre réduit de données représentatives d'un ou plusieurs vols réels d'aéronefs (notamment de l'avion AIRBUS A 310) et mesurables pendant ces vols sans incidence néfaste sur la charge de travail du pilote et sur les conditions de sécurité de vol. Un autre but encore de l'invention est de faire en sorte que le modèle précité puisse Etre utilisé pour simplifier les procédures de certification d'aéronef et en tant qu'aide au développement de nouveaux aéronefs, ou bien encore puisse tre applicable à la sélection, la formation ou le perfectionnement des pilotes d'aéronef. L'invention vise également à fournir un dispositif permettant une saisie fiable et commode des diverses données servant à l'élaboration et à l'utilisation du modèle en question. La présente invention vise ainsi un procédé 6 156 4 1 d'élaboration d'un modèle statistique pour prédéterminer la charge de travail d'un pilote d'aéronef, notamment d'un avion de type commercial, en fonction de paramètres - objectifs de vol, dans lequel au cours d'au moins un vol, réel ou simulé, de l'aéronef, on demande au pilote de l'aéronef, dans des phases et conditions du vol et dans des conditions de pilotage déterminées constituant les paramètres de vol, de donner une cotation de sa charge de travail au moyen d'une échelle préétablie, on consigne cette cotation en mme temps que les paramètres de vol correspondants et on mesure simultanément la valeur de données significatives desdits paramètres de vol, puis on établit finalement une corrélation statistique entre les valeurs cotées de la charge de travail d'une part et les données précitées et les paramètres de vol d'autre part. En établissant une corrélation statistique entre les valeurs cotées de charge de travail et des données indicatives de la phase de vol et des conditions de vol et de pilotage, la présente invention permet d'élaborer un modèle capable de fournir une estimation objective de la charge de travail d'un pilote d'aéronef à partir de données objectives de vols dont la mesure en vol n'influe pas sur les conditions de pilotage et sur la sécurité du vol. Un modèle pouvant Etre mis en oeuvre tant au cours de vols réguliers qu'au cours de vols simulés, sans exigences déraisonnables d'instrumentation et sans gqner le pilotage de l'aéronef et recouvrant une large gamme de scénarios de vols normaux ou de pannes, peut en particulier aider à simplifier le long et coûteux processus de vols d'essais spéciaux visant à démontrer qu'un nouvel aéronef peut Etre piloté avec un niveau de charge de travail acceptablement bas. De mme, un modèle répondant aux exigences scientifiques de validité et de fiabilité et pratique & utiliser présentera un intér2t considérable s'il est appliqué au processus de développement d'un nouvel aéronef car il permetLra de déterminer artificiellement les conséquences des choix conceptuels sur la charge de travail des pilotes. Suivant une caractéristique avantageuse de l'invention, le procédé précité est caractérisé en ce qu'on relie la charge de travail prédéterminée aux paramètres de vol par une relation du type: Yi= AI Xil + A2 Xi2 + Ai... + j Xij +... + An Xin ± B1 + B2 +... + Bk +... + Bm + C dans laquelle - Yi est une variable représentative de la charge de travail prédéterminée; lés Xij sont des variables représentatives de données significatives d'un premier ensemble de paramètres de vol, - les constantes Aji Bk et C sont des coefficients multiplicatifs ou additifs caractéristiques du modèle précité, la valeur des coefficients Bk étant indicative d'un second ensemble de paramètres de vol, et en ce qu'on calcule les coefficients Ai, Bk et C en appliquant sur la base de la relation précitée, les méthodes statistiques dites de régression linéaire multiple et d'analyse de la variance à partir de valeurs cotées par le pilote de sa charge de travail et des données précitées et des paramètres de vol mesurés ou relevés au cours des vols d'essai. La validité de la relation statistique précitée entre les cotations de charge de travail faites par le pilote et les paramètres de vol objectifs a été vérifiée par la demanderesse. Celle-ci a comparé les valeurs prédéterminées de charge de travail calculées à partir de la formule précitée en utilisant des données recueillies lors de l'ensemble des vols d'essai réels regroupant des conditions de v- et de pilotage variées. Les valeurs calculées de charge de travail ont été comparées aux valeurs cotées par le pilote lors des vols précités et une excellente correspondance a pu tre établie entre les valeurs cotées et les valeurs calculées. 2 6 1 5 6 4 1 Suivant l'invention également, le modèle obtenu selon le procédé précité est caractérisé par la relation suivante: Yi = AI Xil + A2 Xi2 +... + Aj Xij + + An Xin + B1 + B2 +.. + Bk +... Bm + C dans laquelle: - Yi est une variable représentative de la charge de travail prédéterminée, - les Xij sont des variables représentatives de données significatives d'un premier ensemble de paramètres de vol, - les constantes Ai, Bk et C sont des coefficients multiplicatifs ou additifs caractéristiques du modèle, la valeur des coefficients Bk étant indicative d'un deuxième ensemble de paramètres de vol. De préférence, la charge de travail prédéterminée PRi s'exprime suivant une échelle de cotation comportant six indices de 2 à 7! en fonction de la variable Yi par la relation: PRi = 3 x 10Oî 1, les variables Xij sont indicatives de données de vol relatives à l'aéronef et au rythme et à la variabilité cardiaques des pilotes tandis que les paramètres dont dépend la valeur des coefficients Bk sont les conditions de vol et de pilotage de l'aéronef ainsi que la phase de vol à l'instant correspondant à la relevée des données précitées. En élaborant un modèle mettant en oeuvre des données et/ou paramètres de vol aisément accessibles ou mesurables pendant des opérations de vol normales ou anormales, on obtient un outil de travail pratique à utiliser. Au moyen d'un équipement de calcul approprié, le modèle suivant l'invention pourra ainsi tre utilisé pour 261i641 fournir des cotations objectives de charge de travail instantanées en vol ou sur simulateur. Suivant l'invention également, le dispositif pour la mise en oeuvre du procédé précité est caractérisé en ce- qu'il comprend des moyens pour consigner les estimations de charge de travail des membres d'équipage, des moyens pour interroger sélectivement les membres d'équipage quant & leur charge de travail et consigner simultanément la phase de vol et la condition de vol correspondantes, l'ensemble des moyens précites étant relié à des moyens de traitement et d'enregistrement des données fournies par ledit ensemble de moyens en synchronisme avec une base de temps de référence et des moyens de consignation et/ou d'enregistrement des paramètres de vol en synchronisme avec ladite base de temps de référence. De préférence, le dispositif précité comprend en outre des moyens de mesure et d'enregistrement de l'électrocardiogramme des membres d'équipage en synchronisme avec la base de temps précitée. Un tel dispositif est ainsi particulièrement bien adapté & la mise en oeuvre du procédé précité. Selon d'autres aspects, l'invention vise une application du modèle précité au développement d'un nouvel aéronef dans laquelle on utilise la courbe Yi = f(Xij) pour simuler l'influence d'une nouvelle caractéristique de l'aéronef sur la charge de travail de ses pilotes, une application de ce modèle à la formation, au perfectionnement, au contr8le des performances et au suivi du degré d'adaptation à un nouvel aéronef d'un pilote d'aéronef, dans laquelle on évalue les performances dudit pilote en calculant sa charge de travail prédéterminée dans une situation de vol déterminée, et son application à la certification de la composition de l'équipage d'un nouvel aéronef, caractérisée en ce qu'au cours d'au moins un vol d'essai du nouvel aéronef, on mesure dans différentes conditions et phases de vol représentatives des règles de certification, les paramètres de vol dont sont indicatives 2 6 1 5 6 4 1 o les variables Xiji, on calcule les valeurs correspondantes de la variable Yi à partir desquelles on déduit les valeurs correspondantes de charge de travail relativement à une échelle déterminée, ces valeurs de charge de travail étant ensuite comparées à celles d'un aéronef de référence précédemment certifié. D'autres caractéristiques et avantages de l'invention ressortiront encore de la description qui va suivre. 1i0 Aux dessins annexés donnés à titre d'exemples non limitatifs: - la figure 1 est un organigramme illustrant schématiquement le procédé conforme à l'invention, - la figure 2 est un schéma synoptique du dispositif pour la mise en oeuvre du procédé selon l'invention, - les figures 3A et 3B sont des vues de la face avant du boîtier "observateur" utilisé pour la cotation de la charge de travail, - la figure 4 est une vue de la face avant du boîtier "membre d'équipage" utilisé pour la cotation de la charge de travail, - la figure 5 est une vue schématique partielle de la cabine de pilotage montrant l'implantation du boîtier de la figure 5, - la figure 6 est une courbe illustrant l'enregistrement de l'électrocardiogramme des membres d'équipage, - la figure 6A est une courbe illustrant un exemple de régression sur 64 battements, - les figures 7 et 8 sont des histogrammes illustrant la variation de la fréquence cardiaque des membres d'équipage selon les phases de vol, - la figure 9 est un organigramme illustrant de manière plus détaillée le procédé selon l'invention, - la figure 10 est un organigramme illustrant la validation du modèle suivant l'invention selon le concept 261'5641 des "moitiés partagées". et - les figures 11 à 15 sont des courbes comparatives illustrant la correspondance entre les valeurs de charge de travail cotées subjectivement et celles calculées au moyen du modèle selon l'invention. Les tableaux I à VI insérés dans le mémoire descriptif doivent Etre considérés comme faisant partie intégrante de cette description ainsi que les légendes apparaissant sur les organigrammes des figures 1, 9 et 10. Dans la présente description, donnée à titre non limitatif, on va décrire plus particulièrement un procédé d'élaboration d'un modèle pour prédéterminer la charge de travail d'un pilote d'avion de transport commercial tel que l'AIRBUS A 510. Dans l'exposé détaillé de l'invention qui va tre maintenant présenté en référence aux dessins et aux tableaux, on décrira d'abord la nature et le mode de saisie des données utilisées pour l'élaboration du modèle puis les mesures ou variables dérivées de ces données. On procédera ainsi pour chaque ensemble de données concerné. La procédure d'obtention d'une certification de type pour un avion du genre précité implique l'exécution d'un grand nombre de vols d'essai destines en particulier à déterminer si la charge de travail mentale subjective imposée à ses pilotes (pilote et copilote) par le pilotage de ce nouvel avion est acceptable. Durant ces vols! on demande à l'équipage de coter sa charge de travail à des instants prédéterminés, au moyen d'une échelle préétablie. Ces vols d'essai étant longs et coûteux, la demanderesse a décidé de développer un modèle mathématique capable de prédéterminer la cotation des pilotes de manière fiable et objective. L'étape préliminaire du procédé d'élaboration d'un modèledu type précité consiste à déterminer s'il existe une corrélation statistique entre les cotations faites par les pilotes durant le vol et des performances de l'avion qui soient directement mesurables. Au cours d'une étude préliminaire réalisée par la demanderesse, l'existence d'une telle corrélation a été établie. L'étape suivante consiste à collecter suffisamment de données et à calculer le modèle en appliquant des techniques statistiques appropriées. Dans l'exemple, cette seconde étape a été mise en oeuvre en utilisant les données collectées lors de la campagne de certification de l'avion AIRBUS A 310-200. Ces données constituent en effet une excellente base pour l'élaboration du modèle. Ainsi, les 60 vols de la campagne de certification A 310 procurent une abondance de données parmi lesquelles des cotations subjectives de charge de travail. des mesures de performances de l'avion (ci-après dénommées données avion) et les enregistrements de la fréquence cardiaque (électrocardiogramme) des 14 pilotes avant participé à cette campagne. Ainsi, comme on le décrira plus loin, des variables indicatives de données avion et physiologiques et pouvant tre reliées à la charge de travail, pourront être calculées. En outre, les scénarios de vol utilisés pendant la campagne de lrA 310 couvrent une gamme de conditions de vol normales, anormales et d'urgence relativement large qui peut tre considérée comme une gamme représentative des situations qu'un modèle utile doit être capable de prédéterminer. Pour la collecte des données, on procède comme indiqué à la figure 1. Plus explicitement, pour chaque vol de certification considéré, on demande à chaque pilote de l'avion, dans des phases et conditions du vol et dans des conditions de pilotage déterminées, constituant les paramètres de vol, de donner une cotation de sa charge de travail au moyen d'une échelle préétablie qui sera décrite plus loin, on consigne cette cotation en mime temps que les paramètres de vol correspondants et on mesure simultanément la valeur de données significatives desdits paramètres de vol, puis on établit finalement une corrélation statistique entre les valeurs cotées de la charge de travail d'une part et les données précitées et les paramètres de vol d'autre part. 2 6 1 5 6 4 Le modèle préférentiel de la présente invention. sélectionné en fonction de sa puissance statistique et de considérations pratiques, utilise les données avion et les informations de rythme cardiaque des deux membres d'équipage -(pilote et copilote). Un modèle n'utilisant que les données avion a également été élaboré car la mesure du rythme cardiaque en vol peut 9tre considérée, à tort ou à raison, comme une charge supplémentaire en vol. Ce modèle simplifié, que l'on décrira plus loin, constitue une alternative raisonnable dans des situations qui font obstacle à l'enregistrement ou à l'utilisation des données de rythme cardiaque. Toutefois! ses capacités de prédiction sont bien inférieures à celles d'un modèle utilisant des données avion et de rythme cardiaque. Dans le procédé illustré par les figures 1 et 9, on synchronise et fusionne quatre ensembles distincts de données pour créer les mesures finales avec lesquelles le modèle est construit. Trois de ces ensembles sont fournis sur bande magnétique. Ce sont les cotations des pilotes et les notations d'un observateur embarqué sur les phases de vol et les conditions de vol, les données de vol de l'avion et les données cardiaques des pilotes. Le quatrième ensemble de données comprend des notes de bord imprimées lors des différents vols et indiquant les phases de vol (montée, croisière, descente etc.) et les conditions de vol (normale! urgence etc.). On va maintenant décrire en détail chaque ensemble de données en se référant à chaque fois, en ce qui concerne les moyens de saisie de données, à la figure 2. NATURE ET SAISIE DES DONNES DE COTATION DE CHARGE DE TRAVAIL Pour évaluer leur charge de travail, les pilotes utilisent l'échelle suivante: 2 Charge de travail légère; capacité de réserve ample; très acceptable! 3 Charge de travail modérée; capacité de réserve plus que suffisante -(adéquate); tout à fait 6 1 5 6 4 1 acceptable. 4 Charge de travail normale i capacité de réserve suffisante; acceptable, Charge de travail importante: capacité de S réserve diminuée/réduite; élevée mais accep- table, 6 Charge de travail lourde; capacité de réserve faible; juste acceptable, 7 Charge de travail extrême; capacité de réserve à peu près nulle; non acceptable de façon continue. L'échelle précitée peut être précisée de la fagçon suivante: 2:toutes les tâches du pilote sont accomplies sans délai, -: certaines taches peuvent Etre interrompues, 4: interruptions de taches régulières auxquelles le pilote peut faire face sans effort déraisonnable, :allégement des taches souhaitables; interruptions de taches répétées et/ou niveau d'effort mental ou de stress prononceé, 6: charge de travail supérieure à ce que le pilote peut accepter continuellement dans des conditions normales; interruptions de tâches fréquentes et/ou niveau d'effort mental ou de stress accentué 7: tâches continuellement interrompues et/ou niveau d'effort mental ou de stress très accentué; risque élevé d'erreur ou d'omission. Pour coter, en vol, sa charge de travail au moyen de l'échelle précitée, chaque pilote dispose d'un boîtier 1 (voir figures 2, 4 et 5) placé en face de lui et à portée de main sur le tableau de bord 9 de la cabine de pilotage. Chaque boîtier pilote 1 présente en face avant six boutons poussoirs à impulsion la numérotés de 2 à 7 pour la cotation et un voyant lumineux lb qui s'allume lorsqu'une cotation est demandée au pilote. Chaque boîtier 1 est électriquement relié (voir figure 2) à un boîtier de notation observateur 2 261S641 dont la face avant est représentée à la figure 3A. Ce boîtier observateur 2 est par exemple installé dans la cabine de pilotage à portée de main de l'observateur embarqué. Ce boîtier 2 présente en face avant un bouton inverseur marche/arrêt 3, un afficheur numérique à deux digits 4 et seize boutons poussoirs 5 répérés respectivement E (effacement), CM1 (membre d'équipage i occupant la place de gauche). CM2 (membre d'équipage 2 occupant la place de droite), CS, DR, R, DEC, M,.C, DES, AP, AT, N, NC, PU, CU. Le bouton poussoir E est un bouton lumineux à voyant incorporé. L'inverseur 3 est utilisé pour la mise sous tension des boîtiers 1 et 2, la remise à zéro de l'afficheur 4 et la remise à zéro générale. Le bouton poussoir repéré E (effacement) s'allume lorsqu'on appuie sur les boutons CM1 ou CM2 ce qui signifie qu'une cotation a été demandée à l'un des deuxpilotes. Dans ce cas, le voyant lb du boîtier pilote i correspondant s'allume puis s'éteint une fois que le pilote a appuyé sur l'un des boutons de cotation la. Le bouton E du boîtier 2 s'allume également lorsqu'on appuie sur le bouton CS (condition de vol seulement). Il s'éteint: - 5 secondes après extinction du voyant pilote lb si on demande une cotation (ou 5 secondes après la dernière commande de l'observateur si elle est postérieure); - 5 secondes après la dernière commande de l'observateur lorsque l'on ne demande pas de cotation; - Dès qu'il est enfoncé. Les 5 secondes de temporisation peuvent tre ajustées (de 2s/minimum à lOs/maximum). L'afficheur numérique 4 est incrémenté lorsqu'on appuie sur les boutons CMI, CM2, ou CS et est remis à zéro par l'inverseur 3. Le bouton poussoir à impulsion CS, lorsqu'il est enfoncé, génère une information "pilote non testé" sans allumer les voyants lb des boîtiers 1 PILOTE ou COPILOTE et indique le changement de phase de vol et de condition de vol sans cotation du pilote ou copilote. La séquence se poursuit par la commande des boutons poussoirs phase de vol et condition de vol (voir figure 3). Les boutons poussoirs "phase de vol" commandent par tout moyen approprié connu en lui-mme la génération d'un code, par exemple binaire, de valeur 0 à 7 suivant la correspondance suivante: - 0 pour le bouton DR (remise des gaz ou descente rapide), - 1 pour le bouton R (roulage avant décollage>), - 2 pour le bouton DEC (décollage), - 3 pour le bouton M (montée), - 4 pour le bouton C (croisière)! - 5 pour le bouton DES (descente), - 6 pour le bouton AP (approche), - 7 pour le bouton AT (atterrissage). Les boutons poussoirs "condition de vol" commandent la génération d'un code de valeur 0 à 3 suivant la correspondance suivante: 0: normal (bouton N) 1: liste de vérification normale (bouton NC) 2: panne - urgence (bouton PU) 3: liste de vérification d'urgence (bouton CU) Le bouton N n'est pas connecté. Il est toutefois présent pour que l'observateur appuie, à chaque observation, sur un bouton poussoir 5 des trois catégories précitées. Comme on le voit à la figure 2, les bottiers pilote i et observateur 2 sont électriquement reliés à un boîtier de traitement logique 6 qui, en liaison avec un enregistreur 7 et une horloge 8 définissant une base de temps commune aux différents moyens de mesure de l'avion (heure GMT), permet la centralisation et l'enregistrement des cotations des pilotes et des notations de l'observateur. Le bottier 6 regoit l'alimentation 28V continu du réseau avion et assure l'alimentation des boîtiers 1 et 2. Il génère par ailleurs des signaux logiques représentatifs des cotations et notations précitées qui sont appliquées au boîtier d'acquisition de signaux 7a de l'enregistreur 7. Le boîtier 6 délivre ses signaux sous forme parallèle sur 23 sorties "tout ou rien" 6a. Comme on le décrira plus loin, re boîtier 6 est également relié à un système de télévision embarqué utilisé pour la consignation des paramètres de vol. Alternativement, l'observateur peut disposer d'un boîtier 2 dont la face avant est représentée à la figure 3B. Ce boîtier est plus particulièrement adapté à vérifier l'adéquation des cotations faites par les pilotes. Comme le boîtier de la figure 3A, il comprend un inverseur 3, trois boutons poussoirs 5 repérés respectivement E, CM1 et CM2 et un afficheur 4 qui ont des fonctions identiques aux éléments décrits plus haut. Le boîtier 2 de la figure 3B comprend également six boutons poussoirs 5 repérés de 2 à 7, trois boutons poussoirs 5b repérés M, I et S et un voyant lumineux 2a. Tout au long du vol, l'observateur enregistre au moyen du boîtier 2 précité son appréciation propre des charges de travail respectives des deux membres d'équipage en utilisant les boutons poussoirs 5 repérés 2 à 7. L'observateur procède ainsi à la fin de chaque phase ou sous- phase de vol ou chaque fois qu'il détecte un changement de la charge de travail de l'un ou l'autre des pilotes. Il réalisera également une cotation en réponse à une illumination du voyant 2a qui indique que l'un des deux membres d'équipage vient lui-même de réaliser une cotation. L'observateur peut également consigner les erreurs de pilotage en utilisant les boutons M, I et S pour: M: erreur mineure I: erreur importante telle que anomalie corrigée ou non corrigée ou erreur non relative à la sécurité du vol non corrigée S: erreur relative à la sécurité, corrigée ou non. Afin de remplir ses fonctions de manière aussi efficace que possible, l'observateur doit Etre lui-même un pilote parfaitement familiarisé avec l'échelle de cotation utilisée et les différences entre les tâches opérationnelles, de surveillance et de planification. VARIABLE DERIVEE DES DONNES DE COTATION Une étude des cotations faites par les pilotes a mis en évidence des inégalités importantes dans la répartition de ces cotations en fonction du pilote considéré. L'importance des différences de pilote à pilote dans l'utilisation de l'échelle de cotation conduit à la nécessité d'interpréter ces cotations d'une manière qui annule les variations extérieures dues à la façon dont les pilotes interprètent l'échelle de cotation. Les travaux de GOPHER et BRAUNE (1984) sur la psychophysique des charges de travail ont permis de retenir la relation de transformation suivante: P - Pmin - 1 TLog10 () dans laquelle Pmax - Pmin T est la cotation transformée, P est la cotation brute, Pmin est la cotation minimale donnée par un pilote au cours de tous les vols d'essai, Pmax est la cotation maximale donnée par ce pilote au cours des mîmes vols. La transformation des cotations brutes selon GOPHER et BRAUNE présente au moins deux avantages. D'une part l'adéquation du modèle développé aux données disponibles est améliorée. D'autre part les variations "parasites" dues aux différences de cotation de pilote à pilote sont éliminées. Par ailleurs, compte tenu de l'ensemble des cotations relevées pour tous les pilotes, la valeur de Pmin peut tre fixée égale à 2. Un troisième avantage de la cotation transformée est qu'elle peut Etre considérée comme un nouvel indice de performance du système pilote-avion, d'application plus large que l'échelle utilisée pour la cotation. Toutefois, l'effet de "lissage" de la transformée précitée conduit à une certaine ambiguité si on désire exprimer la charge de travail prédéterminée suivant l'échelle de 2 à 7. On s'est donc orienté vers la relation suivante = PRi = (10iY) x R + 1 dans laquelle: PRi est la cotation prédéterminée sur l'échelle de 2 à 7, Yi est la cotation prédéterminée transformée, c'est-à-dire la cotation calculée, O1 R est la plage de cotation du pilote considéré. Afin de disposer d'une relation applicable pour tous les pilotes, on retient pour R une valeur égale à 3, ce qui conduit à la relation: PRi = 3xlOt' + 1 TRAITEMENT INFORMATIQUE DES DONNES DE COTATION Les cotations des pilotes et les notations de l'observateur, classées par fichier (un fichier par vol), sont enregistrées sur bande magnétique sous. forme binaire sur l'enregistreur 7, à raison de 25 données par seconde. Ces données sont préalablement traitées par ordinateur (non représenté), selon l'organigramme ci-après, pour déterminer l'instant auquel une cotation a été demandée, l'instant à laquelle cette cotation a été donnée, la valeur de la cotation, le pilote donnant la cotation (CM1 ou CM2) et la condition de vol et la phase de vol correspondantes. On obtient ainsi un fichier "cotation" par vol. PROGRAMME DON N EE5 COTATION LIRE LES DONNEES DE COTATION A PARTIR DES BANDES AVION DETERMINER L'INSTANT ET LA VALEUR DE LA COTATION ET L'INTERVALLE DE TEMPS ENTRE LA DEMANDE DE COTATION ET LA REPONSE DETERMINER LA COND[TION OBSERVATEUR ET LES CODES DES PHASES DE VOL CREER UN FICHIER DESTINE A ETRE CHARGE COMME ENTREE DU PROGRAMME COTATION Chaque fichier précité est déchargé de l'unité centrale de l'ordinateur utilisé pour le prétraitement précité vers un microcalculateur (non représenté) pour un nouveau traitement selon l'organigramme ci-après. PROGRAMME COTATION LIRE LES ENREGISTREMENTS DE COTATION VOL PAR VOL QUI ONT ETE PREALABLEMENT FORMATES A PARTIR DE LA BANDE ORIGINALE VERIFIER QUE L'INFORMATION DE VOL EST CORRECTE AJOUTER LE NUMERO DU VOL i AJOUTER LE CODE DE CONDITION DE VOL FICHIER DE COTATIONS DESTINE AU PROGRAMME MODELE Le microcalculateur précité est programmé pour introduire les particularités du vol, c'est-à-dire dans l'exemple ici décrit: - la condition de vol: panne/urgence, liste de verification de panne, liste de vérification normale, vol normal, - la phase de vol: remise des gaz ou descente d'urgence, préroulage, roulage avant décollage, décollage, montée, croisière, descente, approche, atterrissage, roulage après atterrissage, - le r8le du pilote: pilote CM1 (c'est-à-dire occupant le siège de gauche dans la cabine de pilotage) aux commandes de l'avion, pilote CM1 non aux commandes, pilote CM2 (c'est-à-dire occupant le siège de droite dans la cabine de pilotage) aux commandes de l'avion, pilote CM2 non aux commandes; - le nombre de systèmes ECAM en fonctionnement, c'est-à-dire le nombre de systèmes électroniques de gestion centralisée de l'avion (ECAM) fonctionnant au moment de la cotation. - le nombre de systèmes FMS en fonctionnement. c'est-à-dire le nombre de systèmes de gestion de vol (FMSB en fonctionnement au moment de la cotation, - l'état de l'automanette (enclenchée ou non). c'est-à-dire de la manette permettant la régulation automatique de la puissance moteur. Ces paramètres sont fournis en partie par les notations de l'observateur en cabine de pilotage et en partie par des notations consignées par un ou plusieurs observateurs installés en cabine lors des vols d'essai. Le microcalculateur précité est également programmé pour afficher chaque cotation individuelle, l'instant o elle a été donnée, et l'identité du pilote qui l'a donnée. L'instant de cotation est ensuite comparé aux notations des observateurs en cabine dites données d'observàtions. Ces données d'observations sont de 2 6 1 5 6 4 1 préférence codées en utilisant un code à quatre chiffres. Ce code est donné par le tableau I ci-après. TABLEAU I Code à quatre chiffres Le premier chiffre sert à identifier quatre situations: - 0: panne - 1: liste de vérification relative à une panne 2: liste de vérification de routine - 4: vol normal Lorsque le premier chiffre est égal à 2, la séquence observée est une liste de vérification de routine. Les deux chiffres suivants permettent d'identifier le type de liste de vérification, ou une autre tâche du pilote: - 01: préparation du poste de pilotage - 02: briefing - 03: avant démarrage - 04: après démarrage - 05: roulage - 06: avant décollage - 07: après décollage - 08: descente - 09: approche - 10: final - 11: après atterrissage - 20: navigation - 30: communication radio Lorsque le premier chiffre est égal à O, la séquence observée correspond à un incident de vol ou à une panne; s'il est égal à 1, la séquence observée correspond à la liste de vérification, ou à la procédure relative à la panne ou à l'incident. Dans ces deux cas, les deux chiffres suivants permettent d'identifier le type et le sous-type d'incident ou de panne: 2615641 TABLEAU I (suite) - 1: Centrale électrique (APU) - i: feu moteur ou dommage sévère - 2: APU en feu - 3: deux moteurs en feu - 4: panne d'automanette - 5: pression huile moteur basse - 6: précaution moteur ou coupure obligatoire - 7: feu sur un moteur - 8: moteur en panne - 2: Carburant - 1: alimentation carburant sans pompe - 2: pompes internes ou externes en panne - 3: panne simple de la pompe centrale - 4: panne double de la pompe centrale - 5: panne de transfert automatique entre pompes - 3: Electricité - 1: fumée sur pupitre de vol -2: fumée sur tableau d'équipements principal (phase 1) - 3: fumée batteries - 4: fumée sur équipement de cabine - 5: panne de générateur simple - 6: panne de générateur double - 9: fumée sur tableau principal (phase 2) - 4: Hydraulique - 2: perte de fluide hydraulique - 3: autre panne - 5: Commandes de vol -1: détection de volets coincés ou assymétriques 26 2615641 TABLEAU I (suite) - 6: Train d'atterrissage - 1: train d'atterrissage normal - 2: autre panne - 7: Pneumatique - 1: pannes diverses sur moteur ou APU - 2: autre panne - 8: Air conditionné et pressurisation - 1: panne compartiment cargo - 2: autre panne - 9: Divers - 1: problème de structure - 2: panne de pilote automatique - 3: autre panne Le quatrième chiffre sert à identifier la phase de vol ou de préparation de l'avion. - O: préparation - i: préroulage - 2: roulage avant décollage - 3: décollage - 4: montée - 5: croisière - 6: descente - 7: approche - 8: atterrissage - 9: roulage après atterrissage TABLEAU I (suite) Un certain nombre de codes identifie des séquences de vol ou événements très spécifiques, par exemple: - 4030: atterrissage forcé - 4040: vol difficile - 4050: descente d'urgence - 4060: atterrissage en surcharge - 4070: remise des gaz - 4081: incapacité CM1 - 4082: incapacité CM2 Les données d'observation, enregistrées périodiquement pendant les vols, comprennent des informations relatives à la phase de vol, à la condition de vol et au lieu et à la nature de toute panne ou situation d'urgence. Les instants de cotation et les données d'observation précitées sont alors synchronisés manuellement et le code approprié est entré pour chaque cotation. Les fichiers de cotation sont alors rechargés dans l'unité centrale de l'ordinateur pour tre utilisées dans le programme de traitement principal qu'on décrira plus loin. SAISIE DES PARAMETRES DE VOL Pour l'enregistrement des données d'observation, les observateurs installés en cabine disposent d'un circuit de télévision embarqué qui est représenté schématiquement à la figure 2. Ce circuit comprend essentiellement un écran télévision couleur 10 relié d'une part à une première caméra 11il par l'intermédiaire d'un dispositif mélangeur vidéo 12 et d'une unité de contr8le de caméra lia et d'autre part à un dispositif incrustateur vidéo 13 lui-mime relié au boîtier de traitement logique 6 précédemment décrit. La caméra 11 est une caméra grand angle filmant les postes de pilotage des pilotes tandis que le dispositif incrustateur 13 permet de visualiser sur l'écran 10, en surimposition par rapport à l'image de la caméra 11, les cotations de charge de travail et/ou'les données d'observation émanant respectivement des boitiers pilote 1 et observateur 2 par l'intermédiaire du boîtier 6. Les observateurs installés en cabine peuvent ainsi relever les cotations des pilotes à titre de vérification ultérieure ainsi que les paramètres de vol fournis par l'observateur tout en les comparant aux informations fournies par la caméra 11 sur l'écran 10. Le mélangeur vidéo 12 est également relié à un magnétoscope vidéo 14 luimime relié à une seconde caméra 15 par l'intermédiaire d'un second dispositif incrustateur vidéo 16 et d'une unité de contrSle de caméra 15a. La camera prend par exemple des vues du panneau supérieur et de la console centrale de pilotage qui peuvent ëtre enregistrées par le magnétoscope 14 et/ou visualisées sur l'écran 10 en même temps que celles du poste de pilotage grâce au mélangeur 12. En outre. l'incrustateur 16 qui est relié également à la base de temps de référence 8 de l'avion permet l'enregistrement et/ou la visualisation en surimpression de l'heure GMT. L'ensemble du circuit de télévision embarqué est piloté à partir d'un poste de commande 17 installé en cabine et comprenant une unité de commande de magnétoscope 17a reliée au magnétoscope 14 et une unité de commande de sélection de caméra 17b reliée au mélangeur vidéo 12. NATURE ET SAISIE DES DONNES AVION-MESURES EN RESULTANT Les données avion enregistrées au cours des vols d'essais précités, par exemple sur bande magnétique, au moyen d'un enregistreur 18 en synchronisme avec la base de temps de référence 8. comprennent 22 données de vol et 22 données de commandes de vol automatiques. Pour ces dernières données, l'information utile est l'information "commande enclenchée" ou non au moment de la mesure des données de vol. Les données précitées sont répertoriées dans le tableau II ci-après. TABLEAU II N DNE DONNEE AERONEF UNITE PLAGE UTILE DONNEE 1 Commande de Directeur de Vol: Profondeur Degrés +/- 180 2 Commande de Directeur de Vol: Roulis Degrés +/- 180 3 Débattement de Gouverne de Profondeur Gauche Degrés. +/- 180 4 Débattement de Gouverne de Direction Degrés +/- 180 Position des Volets Degrés 0 à 41 6 Débattement d'Aileron Gauche Degrés +/- 180 7 Débattement d'Aileron Droit Degrés +/- 180 8 Angle de Profondeur Degrés +/- 180 9 Angle de Roulis Degrés +/- 180 Vitesse de Profondeur (Corps) Degrés/Sec +/- 128 il Vitesse de Roulis (Corps) Degrés/Sec +/- 128 12 Vitesse verticale (Barométrique) Pieds/Min + /15000 13 Vitesse Aérodynamique Noeuds 15 à 450 14 Débattement du Plan Fixe Horizontal Degrés +/- 180 Angle d'incidence Degrés +/- 60 16 Accélération Longitudinale (Terre) g +/- 0,5 17 Accélération Verticale (Terre) g -0,5 à 2,5 18 Altitude Barométrique (Std) Pieds -1000 à +50 19 Vitesse N1 du moteur 1 t/m 0 à 4400 Vitesse N1 du moteur 2 t/m O à 4400 21 Angle de Manette des Gaz (Moteur 1) Degrés 15 à 87,5 22 Angle de Manette des Gaz (Moteur 2) Degrés 15 à 87,5 TABLEAU II (Suite) N DE DONNEE DONNEE (COMMANDES DE VOL AUTOMATIQUES) 23 Altitude 24 Vitesse Verticale (Vz) Arrondi 26 Directeur de Décollage (remise des gaz) 27 Directeur de Décollage (Décollage) 28 Changement de Niveau 29 Radioalignement de Site (Glide Slope) Modes longitudinaux Suivi 31 Modes longitudinaux - Capture 32 Navigation 33 Tenue de Cap 34 Sélection de Cap - Radioalignement d'Azimuth (Localizer) 36 Roulage 37 Décollage (Lacet) 38 Roulage (Remise des Gaz) 39 Modes latéraux - Capture VOR 41 Modes latéraux - Suivi 42 Alignement 43 Conditionnement Moteur Pilotage Transparent 44 Commande Conditionnement Moteur Pour les données de vol, les unités et les plages utiles ont également été indiquées au tableau II. Pour l'élaboration du modèle conforme & l'invention, les données de vol suivantes ont finalement été retenues: l'angle de roulis, la vitesse aérodynamique, l'angle d'incidence. la vitesse N1 du moteur 1, le débattement de la gouverne de direction et la vitesse de roulis. Les données brutes ne sont cependant pas utilisées dans l'élaboration du modèle pour au moins deux raisons. D'une part, la valeur brute ou instantanée de la plupart des données au moment de la cotation d'un pilote risque de ne pas donner une image exacte de la donnée à prendre en compte et d'autre part l'expérience a montré que les pilotes avaient tendance à baser leur cotation instantanée sur un ensemble d'impressions intégrées sur une certaine période de temps. Ces observations ont ainsi conduit à décider d'un "lissage" des données et la technique de lissage retenue est la moyenne exponentielle. Dans la moyenne exponentielle, on accorde progressivement moins de poids aux valeurs des données à mesure qu'elles vieillissent. La pondération des valeurs antérieures est exponentielle, ce qui signifie que l'importance de chaque valeur antérieure successive diminue du mme rapport. La formule générale pour une moyenne exponentielle avec une constante de temps T (seconde) est: DT Ao = Ai +T (So - A1) ou AO est la valeur de la moyenne exponentielle calculée,. Al est la valeur précédente de la moyenne exponentielle, DT est l'intervalle de solution, c'est-à- dire la période entre les points de données successifs (qui est une seconde pour les données avion et l'intervalle RR réel pour les données de rythme cardiaque ainsi qu'on le décrira plus loin), T est la constante de temps ou le facteur de lissage de la moyenne en seconde, SO est la valeur de la donnée au temps O, c'est-à- dire la valeur correspondant à la cotation, La moyenne exponentielle en donnant plus de poids aux valeurs des données les plus récentes fournit une mesure qui est plus sensible aux événements récents tout en subissant encore les effets d'informations beaucoup plus anciennes. Ces moyennes sont faciles à calculer parce que chaque nouvelle valeur ne dépend que de la valeur présente de la mesure et de la valeur précédente de la moyenne. Elles se sont révélées comme d'excellents moyens de prédiction des cotations de charge de travail des pilotes dans l'étude préliminaire faite par la demanderesse. Selon les cas, les valeurs de T retenues sont 8B 32, 64 et 128 secondes. Par ailleurs, pour le calcul des moyennes exponentielles, la valeur absolue de certaines données est préférée. Par exemple, pour le débattement de la gouverne de direction qui peut prendre des valeurs positives et négatives, il est préférable de relier à la charge de travail l'amplitude absolue de l'écart de valeur par rapport à la valeur nominale. On applique également le lissage par moyenne exponentielle à des variables telles que le taux de passage par zéro, le taux d'inversions ou le logarithme décimal de certaines données. Ainsi, pour des données pouvant prendre à la fois des valeurs positives et négatives et pour lesquelles zéro constitue la valeur nominale (par exemple l'angle ou la vitesse de roulis), ce taux est calculé en comparant le signe de la valeur instantanée avec celui de la valeur qui 2 615641 la précède immédiatement. Si ces deux signes sont identiques. une valeur de zéro est introduite dans la moyenne exponentielle correspondante et s'ils sont différents une valeur de 60 est introduite dans ladite moyenne exponentielle pour indiquer une valeur de 60 transitions par minute. Le taux d'inversions est le nombre d'inversions par minute introduit dans une moyenne exponentielle. Une inversion est définie comme un changement dans l'orientation d'une valeur de donnée comparée aux deux valeurs précédentes de cette donnée. Par exemple les suites de valeurs 1-2-3 ou 7-4-1 ne présentent pas d'inversion tandis que la suite 20-10-11 en présente une. * On a également introduit dans l'élaboration du modèle le logarithme décimal du nombre de commandes de vol automatiques enclenchées au moment de la cotation (voir tableau II). L'hypothèse adoptée est que la charge de travail d'un pilote est inversement proportionnelle à l'aide apportée à ce pilote par le système de commandes de vol automatiques et que ce degré d'assistance peut tre évalué par le nombre de commandes enclenchées. NATURE ET SAISIE DES DONNES PHYSIOLOGIQUES Suivant une forme d'exécution préférée ou perfectionnée de l'invention, on prend en compte, pour l'élaboration du modèle, des données physiologiques relatives aux deux pilotes de l'avion testé. Les données physiologiques ou cardiaques sont enregistrées sur bande magnétique sous la forme de battements ou périodes cardiaques séquentiels (intervalles RR) enregistrés en synchronisme avec une horloge 480 Hz. On utilise à cet effet un enregistreur magnétique à pile et cassette 19 (voir figure 2) dont un canal sert à enregistrer l'électrocardiogramme (ECG) du pilote 20, c'est-à-dire la tension entre deux électrodes 19a et 19b fixées sur la peau de ce pilote selon une dérivation dite CM5 (dérivation orientée approximativement suivant le grand axe du coeur). Une électrode 19a est placée au niveau du manubrium central 21 (partie supérieure du sternum) du pilote 20 et l'autre électrode 19b est placée en regard du bord gauche de la cage thoracique, au niveau de la sixième c8te. Simultanément, on enregistre sur une autre voie un signal de fréquence 480 Hz produit par une horloge à quartz. Sur ce meme canal est enregistré un repère temporel, introduit par un observateur à l'aide d'un bouton approprié, du début de l'enregistrement de l'ECG. L'observateur aidé d'un chronométreur tient à jour un fichier d'événements qui constitue un découpage de l'activité du pilote en une suite de segments temporels ou séquences dont chacune constitue une étape ou un événement déterminé de son activité, désignée par le code préétabli décrit plus haut. L'électrocardiogramme se compose d'une suite de courbes correspondant chacune à un cycle cardiaque, et formée de segments successifs désignés conventionnellement par les lettres PQRST (voir figure 6). En particulier, le complexe QRS correspond à la contraction (systole) ventriculaire et comporte une pointe R facile à détecter. Un traitement informatique préliminaire des données enregistrées doit reconnaitre les pointes des ondes R successives et mesure les intervalles de temps qui les séparent conformément à l'organigramme ci-après. PROGRAMME DONNEES PHYSIO. LIRE LES DONNEES CARDIAQUES SUR BANDE MAGNETIQUE CONVERTIR LES TEMPS EN HEURES GMT ET SYNCHRONISATION IDENTIFIER ET ELIMINER LES VALEURS HORS NORME Il I _ CALCULER ET SORTIR LES INTERVALLES RR EN ms., LA DIFFERENCE PAR RAPPORT A LA MOYENNE SUR 16 BATTEMENTS ET L'HEURE GMT A LA SECONDE PRES I CALCULER LES STATISTIQUES PAR VOL ET PILOTE ET POUR LA MOYENNE DE CHAQUE VOL - I CREER UN FICHIER DES BATTEMENTS | CARDIAQUES TRAITES Les données lues sont donc découpées selon les séquences successives spécifiées dans le fichier événement. Et pour chaque séquence on a les informations suivantes: - Identification de la séquence par son code avec heure réelle du début de la séquence, heure relative au début de l'ensemble de l'observation, et durée de la séquence; - Courbe de la FC (fréquence cardiaque) calculée de 5 en 5 secondes (comme l'inverse de la durée moyenne des intervalles RR entrant dans cette période de 5 secondes); - Histogramme des intervalles RR: plus précisément cet histogramme est présenté suivant 24 classes de FC (de 10 à 240 battements/minute) et donne, dans chaque classe, le nombre des complexes QRS détectés, soit une estimation du temps pendant lequel le pilote s'y est tenu. Les figures 7 et 8 représentent deux exemples d'histogramme pour une population donnée de pilotes testés (pilotes aux commandes de l'avion) et permettent d'apprécier l'évolution du rythme cardiaque de ces pilotes en fonction des phases de vol et par là méme l'intérêt de l'introduction de ce facteur dans l'élaboration du modèle conforme à l'invention. Aux figures 7 et 8, les parties hachurées correspondent à un nombre de battements cardiaques par minute supérieur ou égal à 110 qui est une valeur typique en ce domaine. Sur les courbes considérées, on a porté en ordonnée Ie pourcentage de pilotes ayant les valeurs de battements/minute indiquées en abscisse. MESURES OU VARIABLES RESULTANT DES DONNEES PHYSIOLOGIQUES Pour l'élaboration du modèle, les variables suivantes dérivées de la mesure précitée de l'électrocardiogramme sont utilisées î - L'intervalle RR brut lissé par moyenne exponentielle comme décrit plus haut. - Valeur de différence: deux valeurs de différence calculées à partir de l'écart de l'intervalle RR 2 6 1 5641 par rapport à la moyenne linéaire rectiligne des 16 ou 32 intervalles RR précédents. Ceci fournit un type de mesures de mouvements sur 16 et 32 battements. Les moyennes exponentielles lissent ensuite ces mouvements pour en quelque sorte annuler tout artifice. Grâce à ce lissage. l'impact d'un seul cycle cardiaque anormalement haut ou bas, qui se dégagerait des premiers filtrages, sera quelque peu diminué en fonction de la constante de temps particulière utilisée dans la moyenne. - Valeur moyenne: en plus d'tre utilisées pour crier des valeurs de différence, les moyennes des 16 ou 32 battements présentent déjà de l'intérêt sans &tre introduites dans les moyennes exponentielles. Le fait d'effectuer la moyenne des 16 ou 32 derniers battements précédant une période de cotation déterminée est estimé comme un lissage suffisant. De plus, ces mesures en tant que moyennes linéaires pures sont considérées comme une mesure raisonnable de l'amplitude à court terme du rythme cardiaque précédant juste la période de cotation. - Différence par rapport à la moyenne de référence: un écart type standard et moyen pour les conditions de base de chaque vol (2 minutes après le décollage) est calculé pendant le traitement préliminaire conformément à la procédure de traitement de données schématisées par l'organigramme précité du programme "données physio". Cette mesure représente un ensemble de moyennes exponentielles basées sur la soustraction de l'intervalle RR de la moyenne de l'accélération du rythme cardiaque de chaque pilote lors d'événements survenus lors d'une période de vol relativement calme. Les intervalles RR individuels sont soustraits de la moyenne de référence de fagon à ce que les valeurs positives de la mesure puissent indiquer des rythmes cardiaques plus rapides que les rythmes de base (intervalles RR plus courts). * 35 - Poids de la référence: la différence par rapport à la moyenne de référence ne tient pas compte de l'importance de l'écart type standard par rapport à la moyenne de référence. Si l'écart type standard par rapport à la moyenne de référence est important, des écarts relativement importants par rapport à celle-ci peuvent ne pas se révéler significatifs. Par contre, si l'écart type standard par rapport à la moyenne de référence est faible, il faut accorder plus d'importance aux différences entre les intervalles RR et la moyenne de référence. Pour expliquer l'écart type standard par rapport à la base (SD) et la moyenne (M), les variables suivantes sont calculées de la fagon suivante avant d'être utilisées dans les moyennes exponentielles. M - SD x DBL B1 = poids de référence sigma 1 M + SD M - 2 SD x DBL B2 = poids de référence sigma 2 M + 2 SD o Bi est la vaieur de la mesure DBL est la différence par rapport à la moyenne de référence M SD est l'écart type standard par rapport à la moyenne de référence et 2 SD est 2 fois l'écart type standard par rapport à cette moyenne de référence. - Résidus de regression: une ligne de régression linéaire présente les meilleurs moindres carrés adaptés à un ensemble de données. Les résidus d'une ligne de régression représentent les écarts au carré des points réels issus de la ligne tracée. En tant que tels, leur somme est une mesure pure de la variabilité à court terme qui n'est influencée par aucun mouvement linéaire à court ou à long terme. Ces résidus doivent simplement mesurer l'étendue de la variabilité cardiaque à court terme précédant le temps d'une cotation. Les résidus issus de la ligne de régression des 64 battements qui précèdent immédiatement un instant de cotation déterminé et les 32 premiers battements des 64 TABLEAU IIIa Variable X j Logarithme décimal de l'angle de roulis (128)** Taux d'inversions de la vitesse aérodynamique (64) Taux d'inversions de l'angle d'incidence (64) Taux d'inversions de la vitesse N1 du moteur 1(128) Vitesse N1 du moteur N1 (128) Logarithme décimal du débattement de la gouverne de direction (32) Taux de passage par zéro de l'angle de roulis élevé au carré (128) *** Logarithme décimal du nombre de commandes de vol automatiques enclenchées (I) Taux de passage par zéro de la vitesse de roulis (128) Taux de passage par zéro de la vitesse de roulis élevé au carré (128) *** Pilote testé - Différence RR par rapport à la moyenne de référence (128) Pilote testé - Poids de référence Sigma RR2 (64) Pilote testé - Valeur RR brute (8) Pilote testé - Résidus de régression RR sur 32 battements (I) Pilote non-testé - Valeur RR brute (8) Pilote non-testé Différence RR par rapport à la moyenne sur 16 battements (128) Pilote non-testé - Différence RR par rapport à la moyenne de référence (64) Pilote non testé - Poids de référence Sigma RR1 (64) Pilote non-testé Poids de référence Sigma RR2 (64) Pilote non- testé - Différence par rapport à la moyenne sur 32 battements (128) Pilote non testé - Résidus de régression RR sur 32 battements (I) ** La valeur entre parenthèses est la constante de temps en secondes pour la moyenne mobile exponentielle. "I" indique plutôt une valeur instan- tanée qu'une moyenne mobile. *** Ces mesures sont utilisées en multipliant la mesure par elle-même - avant la multiplication par le coefficient. battements qui précèdent cet instant, sont utilisés comme variables. La figure 6A représente un ensemble de 64 intervalles RR et une ligne de régression tracée. La mesure. résiduelle des 64 battements telle qu'elle est utilisée dans la construction du modèle est la somme des distances au carré des points réels de données sur les intervalles RR par rapport à la ligne de régression tracée. Dans la figure 6A, la ligne de régression L tracée montre un mouvement descendant, c'est-à-dire un mouvement vers un rythme cardiaque accéléré (intervalle RR plus court). L'importance de ce mouvement est mesurée par le coefficient de l'équation de régression. Ce coefficient est -défini par la pente de la ligne de régression L. Puisque le mouvement au cours des 32 et 64 battements utilisés pour les calculs ou les résidus présente aussi de l'intért, les coefficients de régression sont également utilisés comme mesures. Il est important de remarquer que l'amplitude des résidus de régression n'est absolument pas liée au mouvement ou pente de la ligne de régression. Le tracé graphique représenté à la figure 6A peut tre soumis à une rotation de 360 0 entrainant ainsi un changement continuel de la pente de la ligne (coefficient de régression) sans aucune influence sur la somme des distances au carré des points individuels P par rapport à la ligne de régression. Les variables dérivées de données avion et de la mesure de l'électrocartiogramme des pilotes et utilisées pour l'élaboration du modèle sont énumérées dans le tableau IIIa ci-après. Dans ce tableau, pilote testé signifie que la mesure a été relevée sur le pilote donnant une cotation et pilote non testé signifie que la mesure a été relevée sur le pilote ne donnant pas de cotation. 2615641- TRAITEMENT PRINCIPAL DES DONNEES Comme on vient de le décrire, l'étape de traitement principal du procédé Visé par l'invention consiste essentiellement à combiner les données avion, de cotation de charge de travail et cardiaques et à calculer un jeu de mesures à partir de ces données préalablement synchronisées par rapport à une base de temps de référence. Cette étape de traitement principal est illustrée par les organigrammes simplifies ci-après. PROGRAMME DONNEES AVION LIRE LES DONNEES AVION ^ PARTIR DES BANDES MAGNETIQUES AVION ( 25 DONNEES PAR SECONDE) ELIMINER LES DONNEES HORS DES PLAGES UTILISEES i ETABLIR UNE MOYENNE POUR OBTENIR UNE DONNEE PAR SECONDE SORTIR LES DONNEES AVION CORRESPONDANTES PROGRAMME REGRESSION SOUS-PROGRAMME DE REGRESSION MULTIPLE DISTRIBUE PAR IMSL, Inc QUI CALCULE UN MODELE DE RE- REGRESSION ET TOUTES LES STATISTIQUES ASSOCIEES - IL EST UTILISE POUR CALCULER LES RESIDUS DE REGRESSION RR SUR 32 ET 64 BATTEMENTS CARDIAQUES [ PROGRAMME FICHIER SAS ! _ ENTRER TOUTES LES MESURES| AFFECTER UNE ETIQUETTE A CHAQUE MESURE CALCULER DES MESURES COMPOSITES CREER UN FICHIER DE SYSTEME SAS ! l FICHIER SAS POUR LE CALCUL DU MODELE PROGRAMME MODELE | LIRE LES DONNEES CARDIAQUES, LES COTATIONS ET LES DONNEES AVION ET SAISIR LES DONNEES MANQUANTES AFFICHER LES DONNEES AVIONS l POUR LES VALEURS VALIDES CALCULER 279 MESURES BASEES SUR LES DONNEES AVION POUR CHAQUE INSTANT DE COTATION CALCULER 54 MESURES DE VARIABILITE CARDIAQUE POUR CHAQUE PILOTE ET A CHAQUE INSTANT DE TATIO I SORTIR LES DONNEES I CREER UN FICHIER DE TOUTES LES MESURES POUR LE PROGRAMME SAS DE CALCUL DU MODELE 2 6 1 5 6 4 1 | PROGRAMME MODELE SAS ENTRER LE FICHIER SAS SUPPRIMER LES COTATIONS COMPORTANT DES DONNEES INCOMPLETES REPERTORIER LES VARIABLES DE CLASSIFICATION ESSENTIELLES __________ I ________ APPLIQUER LA TECHNIQUE DES MODELES LINEAIRES GENERAUX POUR SELECTIONNER LES MESURES SAUVEGARDER UN FICHIER DE COTATIONS REELLES ET PREDETERMINE CORRELER LES COTATIONS REELLES ET PREDETERMINEES I PRIMER LE M ODELE Le traitement commence par la lecture d'un fichier de cotation déterminé. L'enregistrement comprend & la fois l'instant auquel l'observateur a demandé la cotation et l'intervalle de temps entre la demande et la réponse. Seul l'instant auquel l'observateur a demandé la cotation est pris en compte. Une seconde lui est ajoutée pour prendre en compte le temps nécessaire au pilote interrogé pour appuyer sur un bouton la de son bottier 1. On obtient ainsi un temps cible. Une fois le temps cible établi, les données avion sont lues et les mesures correspondantes totalisées jusqu'à ce temps cible. Une fois que le temps cible a été atteint sur la bande avion, les données cardiaques sont entrées pour chaque pilote et les mesures correspondantes calculées. Au total, et dans le cadre d'un exemple de réalisation préférée, 279 mesures sont obtenues à partir des données avion et 54 mesures sont dérivées des données cardiaques de chaque pilote et pour chaque cotation. Afin de supprimer l'éventualité d'avoir trop de variables indépendantes, chaque mesure est soumise à un pré - filtrage à l'aide de techniques statistiques dites de régression multiple et par paliers. Un tiers des variables à la fois est introduit dans une régression par paliers contre une mesure appropriée de la cotation du pilote. La technique de régression procure les meilleurs moindres carrés linéaires pour toutes les mesures qui correspondent à ou dépassent un niveau de tolérance pré-défini. Chacune des opérations par paliers rapporte de 15 à 20 mesures qui offrent la meilleure prévision de la cotation. Celles-ci sont ensuite introduites dans une régression finale par paliers pour produire le jeu de mesures utilisé finalement dans l'élaboration du modèle qui en résulte. ELABORATION DU MODELE L'élaboration du modèle proprement dit utilise en particulier la procédure des Modèles Généraux Linéaires (GLM) telle qu'elle est mise en oeuvre dans l'ensemble des statistiques SAS (Institut SAS, CARY, CAROLINE DU NORD, USA. 1982). Cette technique permet la combinaison des mesures numériques continues issues de l'avion (données avion) et des données cardiaques avec des données discrètes, c'est-à- dire les paramètres de vol mentionnés plus haut. En utilisant la procédure précitée, on relie conformément à l'invention, la charge de travail prédéterminée PR* aux paramètres de vol par une relation du type: Yi = Ai Xli + A2 X2i +... + Ai Xii +.. (M) + An Xni + BI + B2 +... + Bk +. + Bm + C dans laquelle - Yi est une variable représentative de la charge de travail prédéterminée PRi et liée à celle-ci par la formule citée plus haut, - les Xij sont des variables représentatives des..DTD: données avion et cardiaques précitées. - les constantes Aij. Bk et C sont des coefficients multiplicatifs ou additifs caractéristiques du modèle. la valeur des coefficients Bk étant indicative des paramètres de vol précités. On calcule les coefficients Ai, Bk et C en appliquant, sur la base de la relation précitée, les méthodes statistiques dites de régression linéaire multiple et d'analyse de la variance à partir des valeurs cotées par le pilote de sa charge de travail et des données précitées et.des paramètres de vol mesurés ou relevés au cours des vols d'essai. La procédure "GLM" met en oeuvre une analyse du modèle de variation (ANOVA) lorsque des variables discrètes sont utilisées. La mathématique de l'analyse de l'approche de variation devient plus complexe que dans la condition du modèle de régression, et plus d'une solution devient possible (il n'y a probablement pas de groupe unique de coefficients). Ceci n'est pas considéré comme un désavantage dans le cadre de l'invention. La caractéristique la plus importante du modèle qui ressort de cette approche est qu'il est simple & calculer. Au lieu de ne trouver qu'un simple groupe de constantes multiplicatives comme dans le modèle de régression, l'analyse du modèle de variation conduit également à des coefficients additifs. Ils ne sont ajoutés à l'estimation du modèle que si la situation correspond à une catégorie particulière de variables discrètes (paramètres de vol). Un exemple peut aider à clarifier la structure d'une analyse de modèle linéaire général de variation tel qu'il est fourni par l'invention. Supposons que le modèle (M) comprenne une variable continue, la position des volets (F), et une variable discrète, par exemple que le pilote soit le CM1 ou le CM2 (appelée désignation du pilote ou D). Le modèle résultant pour l'évaluation des charges de travail comprend alors les termes suivants: - un coefficient qui devra tre multiplié par la position des volets et ajouté soit à: - un premier coefficient additif si le pilote aux commandes était le CM1 au moment de la mesure de F ou à - un second coefficient additif différent si F a été mesurée quand le pilote CM2 était aux commandes et à - une constante d'interception qui est ajoutée aux deux autres termes pour déterminer la prédiction définitive. L'équation du modèle s'écrit alors sous la forme: P = mF + (ai si D = CM1 ou a2 si D = CM2) + i oû P est la charge de travail prédéterminée, m est le coefficient multiplicateur pour la position des volets, les a sont les coefficients additifs pour le CM1 et le CM2 respectivement (les valeurs de D) et i est la constante d'interception. En plus de l'utilisation des variables discrètes uniques dans les modèles de ce type, il est également possible d'utiliser des interactions entre des variables de classification. Par exemple, la contribution d'une interaction de condition de vol et de phase de vol peut &tre introduite comme une variable intéressante pour les prédictions. Cette interaction peut apporter des coefficients additifs pour des associations de phase et de condition, par exemple une urgence ou une panne au décollage, liste de vérification pendant la descente. Comme cela sera montré dans la discussion qui va suivre, les variables de classification et leurs coefficients d'interaction ont joué un r8le important dans le modèle élaboré prouvant encore davantage la valeur de l'approche du modèle linéaire général. CONTROLE STATISTIQUE DE LA VALIDITE DU MODELE Plusieurs moyens d'analyse statistique décrivent la bonne adaptation du modèle linéaire aux données globales. La valeur la plus intéressante est la valeur "Multiple R". Cette statistique est le coefficient de corrélation entre les valeurs prédites de la variable dépendante calculée par le modèle, c'est-à-dire la charge de travail, et les valeurs réelles. Dans l'exemple, le Multiple R représente la corrélation entre les cotations des charges de travail réelles données par les pilotes et transformées en utilisant la formule décrite plus haut et les valeurs prédites de la cotation transformée. La valeur au carré du Multiple R peut être aussi directement interprétée comme la proportion de variation dans la variable dépendante (valeurs transformées) expliquée par le modèle. L'importance statistique du modèle tout entier peut Etre testée en utilisant une valeur F calculée. La valeur F représente le rapport de variation expliqué par le modèle, ses degrés de liberté par rapport à l'erreur ou au reste qu'il n'explique pas, et ses degrés de liberté. Lorsque l'importance d'une valeur F est testée, des niveaux de probabilité de 0,05 (5%. de chances pour qu'une adaptation aussi bonne que celle du modèle ait pu se produire) sont souvent acceptés comme significatifs au point de vue statistique. Des valeurs de probabilité de OiO0 au moins laissent peu de doute quant à l'importance statistique du modèle qui n'était sans doute pas dû au hasard. Un modèle doit &tre significatif au point de vue statistique avant de pourvoir être jugé comme valable. Toutefois, une valeur F significative au point de vue statistique n'est pas une garantie de validité. Elle signifie simplement que le modèle élaboré a su prédire mieux que le hasard les données dont il a été issu. Une validation sur un échantillon indépendant, c'est-à- dire un échantillon de mesures non utilisées dans le calcul du modèle, est nécessaire pour prouver la validité. La méthode des "moitiés partagées" dont les étapes sont schématisées sur la figure 10 est utilisable pour démontrer la validité du modèle. Selon cette méthode, la moitié des données recueillies pendant les vols d'essai est utilisée pour l'élaboration du modèle tandis que l'autre moitié est utilisée pour sa validation c'est-à-dire que pour cette seconde moitié. on calcule les valeurs prédéterminées de charge de travail au moyen de la relation précitée et on les compare aux valeurs réelles des cotations. On peut, en variante, réaliser des vols spéciaux pendant lesquels on relève des données utilisables pour valider le modèle. Le modèle préféré de l'invention est détaillé dans les tableaux IIIb et IV ci-après. TABLEAU IIIb Variable X. COEFFICIENT J MULTIPLICATIF A Logarithme décimal de l'angle de roulis (128)** 0,05771089 Taux d'inversions de la vitesse aérodynamique (64) 0,00114733 Taux d'inversions de l'angle d'incidence (64) -0,00403249 Taux d'inversions de la vitesse N1 du moteur 1(128) -0,00327516 Vitesse N1 du moteur N1 (128) - 0,00107951 Logarithme décimal du débattement de la gouverne de direction (32) 0,00967576 Taux de passage par zéro de l'angle de roulis élevé au carré (128) *** 0,00029009 Logarithme décimal du nombre de commandes de vol automatiques enclenchées (I) -0,13144402 Taux de passage par zéro de la vitesse de roulis (128) 0,00888219 Taux de passage par zéro de la vitesse de roulis élevé au carré (128) *** -0,00026902 Pilote testé - Différence RR par rapport à la moyenne de référence (128)- -0,00202818 Pilote testé Poids de référence Sigma RR2 (64) 0,00365653 Pilote testé - Valeur RR brute (8) 0,00031228 Pilote testé - Résidus de régression RR sur 32 battements (I) 5,3414072E-06 Pilote non-testé - Valeur RR brute (8) -0,00037668 Pilote non-testé - Différence RR par rapport à la moyenne sur 16 battements (128) 0,01038512 Pilote non-testé - Différence RR par rapport à la moyenne de référence (64) 0,11703052 Pilote non testé - Poids de référence Sigma RR1 (64) -0,26649027 Pilote non-testé - Poids de référence Sigma RR2 (64) 0,15181157 Pilote non-testé - Différence par rapport à la moyenne sur 32 battements (128) -0,01318688 Pilote non testé - Résidus de régression RR sur 32 battements (I) -0,00326222 ** La valeur entre parenthèses est la constante de temps en secondes pour la moyenne mobile exponentielle. "I" indique une valeur instantanée. *** Ces mesures sont utilisées en multipliant la mesure par elle-même nvant la miiltinliratinn nar lp coefficient. TABLEAU IV PARAMETRES DE VOL COEFFICIENT ADDITIF Bk Estimation par le pilote CM1 aux commandes -0,01128010 Estimaitiocn par le pilote CM1 non aux commandes -0,03695012 Estimation par le pilote CM2 aux commandes -0,02018463 Estimation par le pilote CM2 non aux commandes 0,00000000 Deux systèmes ECAM en panne (aucun système ECAM en fonctionnement) 0,03846630 Un système ECAM en panne (un système ECAM en fonctionnement) -0,00036457 Deux systèmes ECAM en fonctionnement, 0,00000000 Remise des gaz ou descente rapide avec les deux FMS (système de gestion de vol) en fonctionnement 0,00000000 Préroulage avec les deux FMS en panne -0,01271110 Préroulage avec un FMS et un FMS en fonctionnement -0,22427695 Préroulage avec les deux FMS en fonctionnement 0,00000000 Roulage avant décollage avec les deux FMS en panne 0,01688276 Roulage avant décollage avec un FMS en panne et un FMS en fonctionnement 0,07965894 Roulage avant décollage avec les deux FMS en fonctionnement 0, 00000000 Décollage avec les deux FMS en panne -0,01649307 Décollage avec un FMS en panne et un FMS en fonctionnement -0,02033955 Décollage avec les deux FMS en fonctionnement 0,00000000 Montée avec les deux FMS en panne -0,03001335 Montée avec un FMS en panne et un FMS en fonctionnement 0,01695141 Montée avec les deux FMS en fonctionnement 0,00000000 Croisière avec les deux FMS en panne 0,05614202 Croisière avec un FMS en panne et un FMS en fonctionnement -0,10278220 Croisière avec les deux FMS en fonctionnement 0,00000000 Descente avec les deux FMS en panne - 0, 03471076 Descente avec un FMS en panne et un FMS en fonctionnement 0, 01862299 Descente avec les deux FMS en fonctionnement 0,00000000 2 6 1 5 6 41 TABLEAU IV (suite) PARAMETRE DE VOL COEFFICIENT ADDITIF Bk Panne ou Urgence y compris Incapacité du pilote CM1 0,04277797 Liste de vérification pour Panne ou Urgence 0,06829986 Liste de vérification pendant un vol normal 0,03892804 Vol normal sans liste de vérification.0, 00000000 Remise des Gaz ou Descente d'Urgence 0,05458590 Pré-roulage - 0, 09502411 Roulage avant décollage - 0,07419961 Décollage 0,02750070 Montée 0,02952791 Croisière 0,01699097 Descente 0,00673707 Approche 0,03832675 Atterrissage 0,03482687 Roulage après atterrissage 0,00000000 Approche avec les deux FMS en panne -0,1405300 Approche avec un FMS en panne et un FMS en fonctionnement -0,04206698 Approche avec les deux FMS en fonctionnement 0,00000000 Atterrissage avec les deux FMS en panne -0, 07434618 Atterrissage avec les deux FMS en fonctionnement 0,00000000 Roulage après atterrissage avec les deux FMS en panne -0,11597344 Roulage après atterrissage avec un FMS en panne et un FMS en fonctionnement -0, 07514224 Roulage après atterrissage avec les deux FMS en fonctionnement 0, 00000000 Remise des gaz ou descente d'urgence avec automanette enclenchée 0,00000000 Préroulage avec automanette non enclenchée 0,18152502 Préroulage avec automanette enclenchée 0,00000000 Roulage avec automanette non enclenchée -0,10494458 Roulage avec automanette enclenchée 0,00000000 Décollage avec automanette non enclenchée -0, 3796761 Décollage avec automanette enclenchée 0,00000000 261564 1 TABLEAU IV (suite) Montée avec automanette non enclenchée -0,02778540 Montée avec automanette enclenchée 0,00000000 Croisière avec automanette non enclenchée 0,13777431 Croisière avec automanette enclenchée 0, 00000000 Descente avec automanette non enclenchée 0,07989134 Descente avec automanette enclenchée 0,00000000 Approche avec automanette non enclenchée 0,04462338 Approche avec automanette enclenchée 0,00000000 Atterrissage avec automanette non enclenchée -0,01322323 Atterrissage avec automanette enclenchée 0,00000000 Le tableau IIIb montre les coefficients multiplicatifs Aij et leurs variables ou mesures avion et cardiaques correspondantes. Ces coefficients sont utilisés en multipliant leur valeur par la valeur de la mesure correspondante et en sommant les résultats pour les 21 mesures considérées (10 relatives à l'avion et il de nature cardiaque). Par exemple, le logarithme décimal de la moyenne exponentielle à 128 secondes de l'angle de roulis est multiplié par 0,05771089 et le résultat ajouté au produit de la-moyenne exponentielle à 64 secondes du taux d'inversion de la vitesse aérodynamique par 0,00114733 etc jusqu'à ce que les 21 mesures multipliées par leur coefficient respectif aient été sommées. On notera que deux mesures (la moyenne exponentielle à 128 secondes du taux de passage par zéro de l'angle de roulis et la moyenne exponentielle à 128 secondes du taux de passage par zéro de la vitesse de roulis) interviennent au carré. Le tableau IV fournit le détail des coefficients additifs du modèle. Ces coefficients sont associés à des variables d'état de vol et à leurs interactions. Ces variables sont: - la condition de vol, - la phase de vol, - le r8le du pilote, - le nombre de systèmes ECAM en fonctionnement, c'est-à-dire le nombre de systèmes de gestion électronique centralisée d'avion en fonctionnement au moment de la cotation, - le nombre de systèmes FMS-en fonctionnement, c'est-à-dire le nombre de systèmes de gestion de vol en fonctionnement au moment de la cotation. La condition de vol, la phase de vol, le r8le du pilote et le nombre de systèmes ECAM en fonctionnement présentent un impact prépondérant dans le modèle. En effet, il y a un coefficient additif pour chaque niveau de ces 26 15641 variables d'état. Par exemple, il y a des coefficients additifs distincts pour la montée (une condition de vol), le pilote CMI non aux commandes de l'avion (un niveau de r8le de pilote) et le cas de vol avec un système ECAM en panne. Le nombre de systèmes FMS fonctionnant et l'état de l'automanette interviennent en interagissant avec la phase de vol. Ainsi, un coefficient additif est calculé pour chaque combinaison de phase de vol et du nombre de systèmes FMS ou de l'état de l'automanette qui sont présents dans les données. Les coefficients additifs sont utilisés en les ajoutant à la somme des produits des coefficients multiplicatifs par les mesures de vol associées si et seulement si la cotation des charges de travail correspond à la condition d'état de vol considérée. Par exemple, si la cotation correspond au pilote CMI aux commandes de l'appareil pendant une liste de vérification de routine durant la descente avec l'automanette non enclenchée et tous les systèmes FMS et ECAM en fonctionnement, les coefficients suivants devront Etre ajoutés: - CM1 pilote aux commandes de l'avion:-0,01128010. - liste de vérification durant un vol normal: -0,03892804 - descente: 0,00673707 - les deux systèmes ECAM en fonctionnement: 0,00000000 - descente avec les deux systèmes FMS en fonctionnement: 0,00000000 descente avec l'automanette non enclenchée: 0,07989134. Il existe également une constante dite d'interception ayant une valeur de 0,18425416 qui est ajoutée pour chaque cotation. Le multiple R du modèle en question est égal à 0,664 ce qui, compte tenu du nombre de variables dépendantes Xii utilisé relativement réduit, indique une très bonne corrélation entre la valeur de charge de travail prédéterminée et la valeur cotée. La valeur F pour ce modèle vaut 15,21 avec 63 degrés de liberté. ce qui est significatif avec une probabilité inférieure à 0,0001. C'est-à-dire qu'il y a moins d'une chance sur 10000 pour qu'un modèleaussi bien adapté aux données relevées puisse tre dû au hasard. VALIDATION CONCRETE DU MODELE La validation du modèle précitée sur la base de la méthode des moitiés partagées a également conduit à démontrer la fiabilité du modèle élaboré. Sur la base des propriétés statistiques précitées du modèle, il est possible de conclure en toute confiance qu'un modèle valable a été élaboré. Ce modèle présente également un caractère pratique dans la mesure o les données utilisées pour son élaboration peuvent Etre aisément relevées en vol avec une instrumentation réduite. Dans ces conditions. on imagine que le modèle pourra &tre utilisé pour produire des cotations de charges de travail instantanées en vol ou sur simulateur à l'aide des calculateurs embarqués existants ou d'un microcalculateur programmé à cet effet. Une autre façon d'évaluer le caractère pratique du modèle est d'examiner si les relations qu'il suggère sont conformes à ce qui est réellement vécu pendant les vols. On examine ainsi les coefficients pour voir si les relations suggérées sont raisonnables. Par exemple, le fait que le logarithme du nombre de commandes de vols automatiques enclenchées soit introduit dans le modèle, fournit un coefficient multiplicatif négatif. Ceci indique que la charge de travail est inversement liée au nombre de commandes de vols automatiques enclenchées ou à l'importance de l'aide reçue par les pilotes grace aux commandes de vol automatiques, ce qui est conforme à la réalité. On a également vérifié la validité du modèle en calculant une charge de travail prédéterminée pour chacune des cotations utilisées pour construire le modèle. Les figures 11 à 15 donnent des courbes comparatives de la charge de travail prédéterminée et de la charge de travail réelle ou subjective en fonction de différents paramètres de vol tels que la phase de vol, le r8le du pilote ou les scénarios de vol. Les scénarios de vol en question sont ceux énumérés dans le tableau V ci-après. Ces courbes permettent d'apprécier la précision du modèle élaboré. TABLEAU V DESCRIPTION DES SCENARIOS DE VOL "o DE EQUIPEMENT SCENARIO DE LA MEL (1) PANNE CONSEQUENCES - générateur moteur - commande manuelle de LA Igauche Aucune température poste - commande de temps automatique poste - CRT (2) gauche - commande manuelle affichage naviga- temp. sortie groupe 1B tion Aucune - commande auto temp. groupe 1 LC - CRT ECAM gauche {Aucune - navigation traditionnelle - les deux FMS -TCC (3) - en croisière-panne - activation manuelle 1D - FMS gauche radio altimètre 2 manette des gaz R- adio altimètre 1 - début de descente: - commande manuelle manette 2 - TRP (4) incapacité du pilote des gaz (à la discrétion de - conduite par un seul l'estimateur) pilote - FAC 1 (5) - début de la descente: - vol manuel panne moteur gauche - un moteur hors service 3A- FAC 1 - à 2000 pieds remise des gaz à l'arrivée - FAC 2 - pendant la croisière: - vol manuel panne FMS droit 3B - FCC 2 (6) - début de descente: - un moteur hors service panne moteur gauche - FMS gauche -navigation traditionnelle - 5mn avant la descen- - vol avec informations 4 Aucun te: fumée électroréduites nique (soute électro- - application procédure nique) fumée électronique - les 2 CRT ECAM - pendant la croisière - démonstration de lecture * Panne pompe carburant! et surveillance des sys- Panne réchauffage tèmes avec perturbations A pitot co-pilote de la gestion de vol nor- Panne réchauffage male glaces Panne système press Restitution des. sys- tèmes après identi- fication pannes et application C/L TABLEAU V (suite) - en croisière: panne - démonstration de lecture FWC1 et surveillance des sys- B - SDAC (7) - début de descente: tèmes avec perturbation panne de pressuriza- de la gestion de vol tion normale - FWC1 (8) - pendant la croisière: - démonstration de lecture panne pompe carburant et surveillance des sys- panne réchauffage pi- tèmes avec perturbation tot co-pilote de la gestion de vol panne réchauffage normale C glaces panne système de pre- ssurisation restitution des sys- tèzes après identifi- cation pannes et ap- plication C/L - PMC (9) - juste avantTOD (10) - descente rapide suivie 6 panne du système pressu- d'un atterrissage en - Système pressuri- risation 2 diversion/dégagement sation 1 (1) Liste d'équipement minimale (Minimum equipment list) (2) Tube cathodique (Cathodic ray tube) (3) Calculateur de commande de poussée (Trust control computer) (4) Afficheur de poussée (Trust rating panel) (5) Calculateur de protection du domaine de vol (Flight augmentation computer) (6) Calculateur de commande de vol (Flight control computer) (7) Concentrateur-convertisseur digital analogique (8) Calculateur d'alarme de vol (Flight warning computer) (9) Calculateur de gestion des performances (Performance management computer) (10) Point de descente (Top of descent) 2 6 1 5 6 4 1 Afin de fournir une estimation de l'importance des capacités de prévision qui pourraient &tre perdues si les mesures de rythme cardiaque n'étaient pas introduites dans l'élaboration du modèle, un deuxième modèle n'utilisant que les données avion et des paramètres d'état a été élaboré. Les variables utilisées pour ce modèle et les coefficients correspondants sont énumérés dans le tableau VI ci-après. Bien que ce modèle simplifié se soit révélé significatif au plan statistique, ses capacités de prédiction sont inférieures à celles du modèle préféré de l'invention. Le multiple R dudit modèle n'est que de 0!52794 et cette valeur de R doit Etre jugée trop faible pour pouvoir considérer que le modèle simplifié en question présente une fiabilité suffisante. 2 6 1 5 6 4 1 Tableau VI COEFFICIENTS DU MODELE N'LTILISANT QUE LES DONSEES AVION COEFFICIENTS MULTIPLICATIFS* M.ESULE COEFFICIENT MULTIPLICATIF Logarithme Décimal de l'Angle de Roulis (128)** 0.05834188 Taux d'Inversions de Valeurs anémométriques (64) 0.00093520 Logarithme Décimal d'Angle d'incidence (16) 0.01313701 Taux d'Inversions d'Angle d'Incidence (8) 0.00072914 Taux d'Inversions d'Angle d'Incidence (64) -0.00516258 Ni Moteur # 1 Taux d'Inversions (128) -0.00033918 Angle de la Manette des Gaz # 1 Inversions (128) 0.00113424 Différence Absolue entre Inversions d'Angles de Manette Poussée (64) -0.00089320 Logarithme Décimal du Nombre de Modes AFC (CADV) engagés (I) -0.04726963 Logarithme Décimal du Débattement de l'Empennage Horizontal (32) -0.03979384 Logarithme Décimal de l'Ordre de Tangage du Directeur de Vol (8) -0.01314656 Logarithme Décimal du Débattement de la Gouverne de Direction(32) 0.01824641 Taux d'Inversions de Débattement de la Gouverne de Direction (32) 0.00070730 Taux d'Inversions d'Angle de Tangage (8) -0.00030736 Logarithme Décimal d'Accélération Longitudinale (32) 0.04581542 Valeur jusqu'à zéro inclus de l'Accélération Longitudinale (16) 0.00238442 Taux d'Inversions d'Accélération Longitudinale (16) -0.00060730 Logarithme Décimal de la Position des Volets (I) 0.00325275 Logarithme Décimal de la Vitesse de Tangage (128) 0.00425131 Valeur jusqu'à zéro inclus de la Vitesse de Tangage (64) -8.9090386E-OS Valeur jusqu'à zéro inclus de la Vitesse de Roulis (128) -0.00207023 Taux d'Inversions d'Accélération Verticale (16) 0.00028602 * Ces coefficients sont utilisés en multipliant la valeur absolue de la mesure p2 le coefficient puis en totalisant toutes les résultantes. Les coefficients anne de "E" sont lus comme la valeur précédant "E" multipliée par 10 à la puissance suivant "E". ** La valeur entre parenthèses est la constante de temps en secondes pour la moyel mobile exponentielle. "(I)" indique une valeur instantanée plutôt qu'une moyeni mobile. Tableau VI (suite) COEFFICIENTS DU MODELE N'LTILISANT QUE LES DONNES AVION COEFFICIENTS ADDITIFS *** CONDITION D'ESTIMATION DE CHARGE DE TRAVAIL*** COEFFICIENT ADDIIF Panne ou Urgence y compris Incapacité du CMI 0.087,2674 Checklist pour Panne ou Urgence 0.03175946 Checklist en Vol Normal -0.03447097 Vol Normal sans Checklist 0. 00000000 Remise des Gaz ou Descente en Urgence 0.08589053 Pré-roulage -0.33573995 Roulage avant Décollage -0.04527509 Décollage 0.00167.79 Montée -0.04318535 Croisière -0.00204135 Descente 0. 05925916 Approche 0.04052923 Atterrissage 0.07862515 Roulage après Atterrissage 0.00000000 Panne/Urgence en Pré-Roulage 0.20938602 Panne/Urgence en Roulage avant Décollage 0.30883981 Panne/Urgence en Montée -0.10548680 Panne/Urgence en Croisière 0.02265858 Panne/Urgence en Descente -0.05014735 Panne/Urgence pendant l'Atterrissage 0.02748625 Panne/Urgence en Roulage après Atterrissage 0.0000000C Checklist d'Urgence en Croisière 0.18715066 Checklist d'Urgence en Descente 0. 0332066C Checklist d'Urgence en Approche 0.01275090 Checklist d'Urgence pendant l'Atterrissage 0.00000000 Checklist normale en Pré-Roulage 0. 09255541 Checklist normale en Roulage avant Décollage 0.06003368 Checklist normale en Montée 0.07289595 Checklist normale en Croisière 0. 12555574 Checklist normale en Descente 0.05839350 Checklist normale en Approche 0.03380438 Checklist normale pendant l'Atterrissage -0.02361849 Checklist normale en Roulage après Atterrissage 0.00000000 Remise des Gaz ou Descente d'Urgence pendant les Scénarios 3A ou 3B -0.06047917 Remise des Gaz ou Descente d'Urgence pendant les Scénarios 5A, 5B ou 5C 0. 00547621 Remise des Gaz ou Descente d'Urgence pendant le Scénario 6 0. 00000000 Pré-Roulage pendant les Scénarios lA, lB, 1C ou ID 0.28182377 Pré-Roulage pendant le Scénario 2 0.22478104 Pré-Roulage pendant les Scénarios 3A ou 3B 0.24351973 Pré-Roulage pendant les Scénarios 5A, 5B ou 5C 0.32000730 *** Les coefficients additifs sont totalisés avec la somme globale des produits des coefficients multiplicatifs et leurs mesures si l'estimatioZ des charges de travail se fait dans les circonstances décrites dans la première colonne. 26 15641 Tableau VI (suite) COEFFICIENTS DU MODELE N'UTILISANT QUE LES DONNES AVION COEFFICIENTS ADDITIFS *** CONDITION DE L'ESTIMATION DE LA CHARGE DE TRAVAIL*** COEFFICIENT ADDITIF Roulage avant Décollage pendant les Scénarios 1A, 1B, 1C ou 1D 0.01120371 Roulage avant Décollage pendant le Scénario 2 -0.09609632 Roulage avant Décollage pendant les Scénarios 3A ou 3B -0.01332814 Roulage avant Décollage pendant les Scénarios 5A, 5B ou 5C 0.02574656 Roulage avant Décollage pendant le Scénario 6 0.00000000 Décollage pendant les Scénarios LA, 1B, 1C ou ID -0.03798744 Décollage pendant le Scénario 2 0. 05228940 Décollage pendant les Scénarios 3A ou 3B -0.02070000 Décollage pendant les Scénarios 5A, 5B ou 5C -0.01730159 Décollage pendant le Scénario 6 0.00000000 Montée pendant les Scénarios 1A, 1B, 1C ou ID 0. 00563219 Montée pendant le Scénario 2 -0.03161124 Montée pendant les Scénarios 3A ou 3B 0.01951183 Montée pendant les Scénarios 5A, 5B ou 5C 0. 07368082 Montée pendant le Scénario 6 0.00000000 Croisière pendant les Scénarios LA, 1B, LC ou 1D 0.01689474 Croisière pendant le Scénario 2. 0. 05292911 Croisière pendant les Scénarios 3A ou 3B 0.02314432 Croisière pendant les Scénarios 5A, 5B ou 5C 0.00938397 Croisière pendant le Scénario 6 0.00000000 Descente pendant les Scénarios LA, 1B, 1C ou 1D -0. 03121730 Descente pendant le Scénario 2 0.08713032 Descente pendant les Scénarios 3A ou 3B -0.02917356 Descente pendant les Scénarios 5A, 5B ou 5C -0.06624218 Descente pendant le Scénario 6 0.00000000 Approche pendant les Scénarios 1A, lB, 1C ou 1D 0.00523105 Approche pendant le Scénario 2 0.07847470 Approche pendant les Scénarios 3A ou 3B -0.02288651 Approche pendant les Scénarios 5A, 5B ou 5C -0.00445375 Approche pendant le Scénario 6 0.00000000 Atterrissage pendant les Scénarios 1A, 1B, 1C ou 1D -0.01729109 Atterrissage pendant le Scénario 2 -0.02280544 Atterrissage pendant les Scénarios 3A ou 3B -0.06915574 Atterrissage pendant les Scénarios 5A, 5B ou 5C 0.02389426 Atterrissage pendant le Scénario 6 0. 00000000 Roulage après Atterrissage pendant les Scénarios LA, 1B, 1C ou 1D 0.03157858 Roulage après Atterrissage pendant les Scénarios 3A'ou 3B 0. 08575654 Roulage après Atterrissage pendant les Scénarios 5A, 5B ou 5C 0. 06807901 Roulage après Atterrissage pendant le Scénario 6 0.00000000 Remise des Gaz ou Descente d'Urgence et Quotation du PF -0.03006523 Remise des Gaz ou Descente d'Urgence et Quotation du PNF 0.00000000 Pré- Roulage et Quotation du PF 0.02320635 Pré-Roulage et Quotation du PNF 0. 00000000 Roulage avant Décollage et Quotation du PF -0.03853689 Roulage avant Décollage et Quotation du PNF 0.00000000 Tableau VI (suite) COEFFICIENTS DU MODELE N'LTILISANT QUE LES DONNEES AVION COEFFICIENTS ADDITIFS*** CONDITION DE L'ESTIMATION DE LA CHARGE DE TRAVAIL*** COEFFICIENT ADDITIF Décollage et Quotation du PF 0.03387488 Décollage et Quotation du PNF 0. 00000000 Montée et Quotation du PF -0.01280829 Montée et Quotation du PNF 0.00000005 Croisière et Quotation du PF -0.03603436 Croisière et Quotation du PNF 0.0000000C Descente et Quotation du PF 0.00672002 Descente et Quotation du PNF 0.0000002 Approche et Quotation du PF 0. 00495925 Approche et Quotation du PNF O.00000000 Atterrissage et Quotation du PF 0.04742134 Atterrissage et Quotation du PNF 0.00000000 Roulage après Atterrissage et Quotation du PF -0.04454874 Roulage après Atterrissage et Quotation du PNF 0.00000002 INTERCEPTION (constante ajoutée dans toutes les estimations) 0.11550363 PF Pilote aux commandes PNF Pilote non aux commandes La validité, la fiabilité et le réalisme du modèle préféré élaboré ne lui garantissent pas nécessairement un caractère d'universalité. Toutefois, aucune des variables utilisées pour l'élaboration du modèle ne présente un caractère susceptible d'empêcher ce modèle de s'appliquer à tout avion de ligne ayant des caractéristiques et des performances voisines de celle de l'AIRBUS A 310 utilisé pour l'élaboration du modèle préféré de la présente invention. En outre, la technique des modèles linéaires généraux utilisée peut permettre si nécessaire d'introduire le type d'avion comme variable dépendante. APPLICATIONS DU MODELE La fiabilité prévisionnelle du modèle et sa capacité à suivre les cotations des pilotes dans la gamme complète des conditions de vol normales et de secours et des phases de vol disponibles dans les données, fournissent trois catégories principales d'application. Le modèle peut être utilisé comme aide à la mise au point d'un nouvel avion. Un des avantages essentiels de tout modèle mathématique est sa capacité à faciliter les prises de décision. Cependant, il est important de se rendre compte que les modèles ne font qu'aider cette prise de décision et sont sujets à erreur. Ils manipulent parfaitement les interactions entre des variables relativement nombreuses à condition que les relations entre ces variables soient bien comprises. Aucun effort de mise au point unique du système complexe "machine-humain" représenté par un avion de transport moderne ne pourrait espérer identifier et quantifier toutes ces relations. Les statistiques relatives au modèle indiquent qu'il prévoit nettement mieux que le hasard. Des personnes possédant une connaissance approfondie des systèmes avion devraient donc pouvoir utiliser le modèle pour améliorer le rendement des procédures d'utilisation. La disponibilité d'un ou plusieurs modèles qui récapitulent les interactions complexes entre l'avion, les pilotes et le système de contr8le aérien et leur influence sur le jugement des pilotes quant & leur charge de travail peut permettre d'intervenir favorablement dans le processus de mise au point d'un nouvel avion. En fait plusieurs utilisations différentes sont prévues. D'abord un modèle peut prendre part à n'importe quelle simulation d'avion nouveau ou modifié. Ceci permet aux concepteurs de prendre leurs décisions d'étude de système en fonction de leur influence probable sur la perception que les pilotes ont de leur charge de travail. A cet effet, il est plus facile d'utiliser un modèle qui ne comprend pas de données de rythme cardiaque comme celui décrit dans le tableau VI annexé, puisqu'en principe on ne recueille pas de données de rythme et variabilité cardiaques en condition de vol réel. D'une façon similaire, ces modèles sont utilisables en tant que système d'investigation pour fournir des informations sur l'impact possible de la variation d'un ou plusieurs paramètres du modèle sur les charges de travail. Tous les paramètres de la base de données utilisée dans ce but sont alors maintenus aux valeurs servant à construire le modèle sauf quelques uns que l'on fait systématiquement varier pour déterminer leur effet sur la charge de travail prévue. Par exemple, les mesures liées à la commande moteur comme les inversions de l'angle de la manette des gaz et les inversions de la vitesse NI pourraient varier pour déterminer la réduction de charge de travail possible si un calculateur de commande de poussée (TCC) supplémentaire était inclus pour continuer à assister les pilotes après une panne du TCC. L'examen rapide d'un vaste éventail de possibilités est une faculté particulière des modèles mathématiques. Dans ce cas, il est mime facilité car le calcul du modèle et le stockage de toutes les données correspondantes peuvent s'effectuer sur un micro-calculateur de taille modeste. Le modèle peut également tre utilisé d'une façon inverse pour assister la définition d'objectifs d'étude. Cette troisième approche est similaire à la seconde décrite ci-dessus sauf que la valeur ou l'éventail des charges de travail représente une cible et qu'on fait varier un ou plusieurs paramètres jusqu'à ce que la cible soit atteinte dans les conditions de vol souhaitées. Par exemple, on peut se poser la question suivante: "Quelle est l'importance de la réduction des inversions d'angle de la manette des gaz nécessaire pour maintenir la valeur de la charge de travail prévue à 5 ou moins m&me après des pannes de moteur ou de TCC ?". Une estimation correcte de la réduction nécessaire pourrait!tre obtenue en faisant dérouler des valeurs d'inversion d'angle de la manette des gaz tout en maintenant fixes les autres mesures aux valeurs des scénarios comportant des pannes de moteur et de TCC. En ce qui concerne la formation et/ou le perfectionnement des pilotes! les considérations de charge de travail peuvent présenter au moins trois intér-ts. D'abord lorsque les programmes de formation sont conçus et les protocoles pour les vols mis au point, un éventail de situations de charge de travail est souvent souhaité. Ceci permet d'assurer que les performances des pilotes en formation peuvent &tre observées dans diverses conditions. La disponibilité du modèle et de l'information qu'il contient lui permettent de concevoir des scénarios de formation qui atteignent le niveau souhaité de charge de travail subjective sans avoir besoin de faire faire des vols aux pilotes expérimentés et d'évaluer les scénarios proposés. En second lieu le modèle dans le cadre de la formation peut servir à mesurer le moment o le palier de performance est atteint. On peut supposer que des pilotes inexpérimentés ont des rythmes cardiaques différents de ceux des pilotes expérimentés dont les données sont utilisées pour construire le modèle. En répétant le meme scénario de vol au cours des diverses étapes de la formation et en comparant les résultats avec la prédiction du modèle, il peut &tre possible de déterminer le moment o un pilote est arrivé à la compétence suffisante pour donner des charges de travail subjectives équivalentes à celles des pilotes expérimentés. Evidemment, une base de données normative du résultat du modèle en fonction du pilote devrait tre établie avant que cette utilisation dans le cadre de la formation puisse devenir totalement opérationnelle. Enfin dans le même contexte, le modèle peut servir d'instrument de sélection. Les pilotes peuvent avoir tendance au cours d'entraînement à sous-estimer l'importance de la charge de travail qu'ils sont en train de subir pour tenter d'impressionner leurs instructeurs. Le modèle par contre n'est pas affecté par cette déformation puisqu'il est basé sur des variables que le pilote ne peut pas contr8ler consciemment, en particulier le rythme cardiaque. Le modèle peut donc 8tre capable d'identifier les pilotes qui ressentent un niveau de charge de travail inacceptable même s'ils ne l'admettent pas. Cependant il est nécessaire d'accumuler une expérience opérationnelle considérable et une base de données significative pour la mise en pratique d'une telle utilisation. Ainsi le modèle selon l'invention peut Etre utilisé pour contraler les performances et le degré d'adaptation à un nouvel avion des pilotes. L'application la plus directe du modèle est celle qui remplace le long et coûteux processus qui permet d'obtenir les cotations de charges de travail subjectives au cours de la certification de la composition de l'équipage de nouveaux avions. Le modèle ayant été élaboré dans des conditions de vol de certification présente une validité indiscutable dans ce contexte. Si le modèle est utilisé à la place des vols de certification, les données peuvent #tre acquises pendant que l'avion effectue des essais en vol, par exemple pendant d'autres parties de la campagne de certification et la présentation aux compagnies et aux services officiels. A partir de ces données, un profil des cotations de charges de travail prédites peut Etre obtenu. Ce profil peut tre comparé aux cotations réelles ou prédites des avions précédemment certifiés. Il serait nécessaire d'effectuer des vols supplémentaires spécifiques pour mesurer la charge de travail uniquement afin d'obtenir des prédictions sur les situations, c'est-à- dire les urgences qui ne sont pas démontrées par les vols précités. Le modèle utilisé pendant les vols de certification peut aussi Etre utilisé de manière très profitable pour remplacer les cotations subjectives. Il est démontré que les cotations subjectives sont soumises à des différences individuelles appréciables entre les pilotes y compris à des réactions différentes face à des situations de vol différentes. Actuellement, il n'existe pas d'autre moyen d'expliquer ou de contrôler ces différences que de faire exécuter la mime manoeuvre par plusieurs pilotes et de faire la moyenne de leurs cotations respectives. Si on considère, par exemple, ce qui se passerait si un pilote sur un vol donné lors d'une manoeuvre donnait une cotation de charge de travail subjective de 7, les conséquences d'une erreur de cotation pourraient Etre très sérieuses. Il serait nécessaire de savoir immédiatement si cette cotation était réelle, c'est-à-dire si c'était une cotation à laquelle on pouvait s'attendre de la part d'autres pilotes - mis dans les mmes conditions. Le modèle indiquerait que cette cotation était erronée, ou bien il pourrait mettre en évidence la source de la charge de travail excessive. Le modèle selon l'invention peut encore Utre mis en oeuvre pour déterminer objectivement le nombre -et la composition du personnel navigant technique nécessaire sur un appareil déterminé. La présente invention démontre par conséquent qu'un modèle pouvant prédéterminer la charge de travail d'un pilote d'avion peut tre élaboré à partir de paramètres objectifs de vol et présenter les capacités de prédétermination suffisantes pour tre d'une utilisation Z615641 pratique. Le modèle obtenu est un moyen de prédétermination valable et fiable des cotations subjectives faites par les pilotes. Bien entendu, un tel modèle pourrait tre également utilisé pour prédéterminer la charge de travail d'un pilote d'un autre type d'aéronef tel qu'un hélicoptère ou même encore d'un pilote automobile, de navette spatiale, de fusée ou de navire. REVENDICATIONS 1. Procédé d'élaboration d'un modèle statistique pour prédéterminer la charge de travail d'un pilote d'aéronef, notamment d'un avion de type commercial., en fonction de paramètres objectifs de vol, dans lequel au cours d'au moins un vol, réel ou simulé, de l'aéronef, on demande au pilote de l'aéronef, dans des phases et conditions du vol et dans des conditions de pilotage déterminées, constituant les paramètres de vol, de donner une cotation de sa charge de travail au moyen d'une échelle préétablie, on consigne cette cotation en mLme temps que les paramètres de vol correspondants et on mesure simultanément la valeur de données significatives desdits paramètres de vol, puis on établit finalement une corrélation statistique entre les valeurs cotées de la charge de travail d'une part et les données précitées et les paramètres de vol d'autre part. 2. Procédé conforme à la revendication 1, caractérisé en ce qu'on relie la charge de travail prédéterminée aux paramètres de vol par une relation du type: Yi = AI Xil + A2 Xi2 +... + Ai Xii +.. + An Xin + Bi + B2 +.. + Bk +... + Bm + C (M) dans laquelle: Yi est une variable représentative de la charge de travail prédéterminée, À les Xii sont des variables représentatives des données significatives d'un premier ensemble de paramètres de vol, les constantes Aji, Bk et C sont des coefficients multiplicatifs ou additifs caractéristiques du modèle précité, la valeur des coefficients Bk étant indicative d'un second ensemble de paramètres de vol, et en ce qu'on calcule les coefficients Aji. Bk et C en appliquant, sur la base de la relation précitée, les méthodes statistiques dites de régression linéaire multiple et d'analyse de la variance à partir de valeurs cotées par le pilote de sa charge de travail et des données précitées et des paramètres de vol mesurés ou relevés au cours du vol précité. 3. Procédé conforme à la revendication 2, caractérisé en ce que les variables Xii sont indicatives de O10 données de vol relatives à l'aéronef proprement dit et en ce que les paramètres de vol dont dépend la valeur des coefficients additifs B.. sont les conditions de vol et de pilotage de l'aéronef et la phase de vol à l'instant correspondant à la mesure des données précitées. 4. Procédé conforme à la revendication 3, caractérisé en ce que les données de vol relatives à l'aéronef et les paramètres de vol déterminant les coefficients Bk utilisés pour l'élaboration du modèle précité sont successivement énumérés aux tableaux ci-après: 2 615641 DONNEE AERONEF UNITE PLAGE UTILE Débattement de la gouverne de direction Degrés +/- 180 Angle de roulis Degrés +/- 180 Vitesse de roulis (corps) Degrés/Sec +/- 128 Vitesse aérodynamique Noeuds 15 à 450 Angle d'incidence Degrés +/- 60 Vitesse N1 du moteur 1 t/m O à 4400 N DE DONNEE DONNEE (COMMANDES DE VOL AUTOMATIQUES) 23 Altitude 24 Vitesse Verticale (VZ) Arrondi 26 Directeur de Décollage (remise des gaz) 27 Directeur de Décollage (Décollage) 28 Changement de Niveau 29 Radioalignement de Site (Glide Slope) Modes longitudinaux Suivi 31 Modes longitudinaux - Capture 32 Navigation 33 Tenue de Cap 34 Sélection de Cap Radioalignement d'Azimuth (Localizer) 36 Roulage 37 Décollage (Lacet) 38 Roulage (Remise des Gaz) 39 Modes latéraux - Capture VOR 41 Modes latéraux - Suivi 42 Alignement 43 Conditionnement Moteur Pilotage Transparent 44 Commande Conditionnement Moteur PARAMETRES DE VOL Panne ou Urgence y compris Incapacité du pilote CM1 Liste de vérification pour Panne ou Urgence Liste de vérification pendant un vol normal Vol normal sans Liste de vérification Remise des gaz ou Descente d'Urgence Pré-roulage Roulage avant Décollage Décollage Montée Croisière Descente Approche Atterrissage Roulage après Atterrissage Estimation par le pilote CM1 aux commandes Estimation par le pilote CM1 non aux commandes Estimation par le pilote CM2 aux commandes Estimation par le pilote CM2 non aux commandes Deux systèmes ECAM (système électronique de gestion centralisée d'avion) en panne (aucun système ECAM en fonctionnement) Un système ECAM en panne (un système ECAM en fonctionnement) Deux systèmes ECAM en fonctionnement Remise des gaz ou descente rapide avec les deux FMS (système de gestion de vol) en fonctionnement Préroulage avec les deux FMS en panne Préroulage avec un FMS en panne et un FMS en fonctionnement Préroulage avec les deux FMS en fonctionnement Roulage avant décollage avec les deux FMS en panne Roulage avant décollage avec un FMS en panne et un FMS en fonctionnement Roulage avant décollage avec les deux FMS en fonctionnement Décollage avec les deux FMS en panne Décollage avec un FMS en panne et un FMS en fonctionnement Décollage avec les deux FMS en fonctionnement Montée avec les deux FMS en panne Montée avec un FMS en panne et un FMS en fonctionnement Montée avec les deux FMS en fonctionnement Croisière avec les deux FMS en panne Croisière avec un FMS en panne et un FMS en fonctionnement Croisière avec les deux FMS en fonctionnement Descente avec les deux FMS en panne Descente avec un FMS en panne et un FMS en fonctionnement Descente avec les deux FMS en fonctionnement Approche avec les deux FMS en panne Approche avec un FMS en panne et un FMS en fonctionnement Approche avec les deux FMS en fonctionnement Atterrissage avec les deux FMS en panne Atterrissage avec les deux FMS en fonctionnement Roulage après atterrissage avec les deux FMS en panne Roulage après atterrissage avec un FMS en panne et un FMS en fonctionnement Roulage après atterrissage avec les deux FMS en fonctionnement Remise des gaz ou descente d'urgence avec automanette enclenchée Préroulage avec automanette non enclenchée Préroulage avec automanette enclenchée Roulage avec automanette non enclenchée Roulage avec automanette enclenchée Décollage avec automanette non enclenchée Décollage avec automanette enclenchée Montée avec automanette non enclenchée Montée avec automanette enclenchée Croisière avec automanette non enclenchée Croisière avec automanette enclenchée Descente avec automanette non enclenchée Descente avec automanette enclenchée Approche avec automanette non enclenchée Approche avec automanette enclenchée Atterrissage avec automanette non enclenchée Atterrissage avec automanette enclenchée 2 6 1 5 6 41 5. Procedé conforme a l'une des revendications -5 ou 4, caractérisé en ce que les variables Xij sont également indicatives du rythme et de la variabilité cardiaques du pilote de l'aéronef que l'on mesure pendant ledit vol d'essai lors de la cotation faite par celui-ci. 6. Procédé conforme à l'une des revendications précédentes, caractérisé en ce que pour l'élaboration du modèle, on prend en compte la cotation de charge de travail faite par le pilote qui est aux commandes de l'aéronef et de son co-pilote ainsi que la position du pilote aux commandes de l'aéronef dans le poste de pilotage. 7. Procédé conforme à l'une des revendications précédentes, caractérisé en ce que l'échelle de cotation de charge de travail comprend six indices de cotation s'étendant de 2 à 7 et exprimant une charge de travail croissante. 8. Procédé conforme & la revendication 7, caractérisé en ce que la charge de travail prédéterminée sex-prime suivant l'échelle précitée en fonction de la -20 variable Yi par la relation: PRi = -3 x 1OY + 1 9. Procédé conforme à l'une des revendications 3 à 8, caractérisé en ce que les variables Xii comprennent les variables énumérées dans le tableau ci-après: Variable X.j Logarithme décimal de l'angle de roulis (128) Taux d'inversions de la vitesse aérodynamique (64) Taux d'inversions de l'angle d'incidence (64) Taux d'inversion de la vitesse Ni du moteur 1 Vitesse N1 du moteur 1 (128) Logarithme décimal du débattement de la gouverne de direction (32) Taux de passage par zéro de l'angle de roulis élevé au carré (128) Logarithme décimal du nombre de commandes de vol automatiques enclenchées Taux de passage par zéro de la vitesse de roulis (128) Taux de passage -par zéro de la vitesse de roulis élevé au carré (128) 10. Procédé conforme à l'une des revendications 5 à 9 caractérise en ce que les variables Xij comprennent les variables énumérées dans le tableau ci-après: Variable X 1j Logarithme décimal de l'angle de roulis (128) Taux d'inversions de la vitesse aérodynamique (64) Taux d'inversions de l'angle d'incidence (64) Taux d'inversion de la vitesse Nl du moteur 1 Vitesse N1 du moteur 1 (128) Logarithme décimal du débattement de la gouverne de direction (32) Taux de passage par zéro de l'angle de roulis élevé au carré (128) Logarithme décimal du nombre de commandes de vol automatiques enclenchées Taux de passage par zéro de la vitesse de roulis (128) Taux de passage par zéro de la vitesse de roulis élevé au carré (128) Pilote testé - Différence RR par rapport à la moyenne de référence (128) Pilote testé - Poids de référence Sigma RRR (64) Pilote testé - Valeur RR Brute (8) Pilote testé - Résidus de régression RR sur 32 Battements Pilote non testé - valeur RR brute (8) Pilote non-testé Différence par rapport à la moyenne sur 16 battements (128) Différence par rapport à la moyenne de référence (64) Pilote non testé - Poids de référence Sigma RR1 (64) Pilote non-testé - Poids de référence Sigma RR2 (64) Pilote non-testé - Différence par rapport à la moyenne sur 32 battements (128) Pilote non-testé - Résidus de régression RR sur 32 battements - X.. COEFFICIENT 1J MULTIPLICATIF A Logarithme décimal de l'angle de roulis (128) 0,05771089 Taux d'inversions de la vitesse aérodynamique f64) 0,00114733 Taux d'inversions de l'angle d'incidence (64) -0,00403?49 Taux d'inversions de la vitesse N1 du moteur 1 -0,00327516 Vitesse N1 du moteur N1 (128) -0, 00107951 Logarithme décimal du débattement de la gouverne de direction (32) 0,00967576 Taux de passage par zéro de l'angle de roulis élevé au carré (128) 0,00029009 Logarithme décimal du nombre de commandes de vol automatiques enclenchées -0,13144402 Taux de passage par zéro de la vitesse de roulis (128) 0,00888219 Taux de passage par zéro de la vitesse de roulis élevé au carré (128) -0,00026902 Pilote testé - Différence RR par rapport à la moyenne de référence (128) -0,00202818 Pilote testé - Poids de référence Sigma RRL2 (64) 0,00365653 Pilote testé - Valeur RR brute (8) 0,00031228 Pilote testé - Résidus de régression RR sur 32 battements 5,3414072E-06 Pilote non-testé - Valeur RR brute (8) -0, 00037668 Pilote non-testé - Différence RR par rapport à la moyenne sur 16 battements (128) 0,01038512 Pilote non-testé - Différence RR par rapport à la moyenne de référence (64) 0,11703052 Pilote non testé - Poids de référence Sigma RR1 (64) -0,26649027 Pilote non-testé - Poids de référence Sigma RR2 (64) 0,15181157 Pilote non-testé - Différence par rapport à la moyenne sur 32 battements (128) -0,01318688 Pilote non testé - Résidus de régression RR sur 32 battements -0,00326222 11. Modele statistique obtenu selon le procédé conforme à l'une des revendications précédentes caractérisé par la relation suivante: Yi = A1 Xil + A2 Xi2 +... + Ai Xij +.. + An Xin + B1 + B2 +... Bk +... + Bm + C (M) dans laquelle: Yi est une variable représentative de la charge de travail prédéterminée, les Xij sont des variables représentatives de données significatives d'un premier ensemble de paramètres de vol, les constantes Aij, Bk et C sont des coefficients multiplicatifs ou additifs caractéristiques du modèle, la valeur des coefficients Bk étant indicative d'un deuxième ensemble de paramètres de vol. 12. Modèle conforme à la revendication 11, caractérisé en ce que la charge de travail prédéterminée F'Ri s'exprime suivant une échelle de cotation comportant six indices de 2 à 7, en fonction de la variable Yi par la relation: PRi = 3 x 10 y + 1, en ce que les variables Xii sont indicatives de données de vol relatives à l'aéronef et au rythme et à la variabilité cardiaques des pilotes et en ce que les paramètres dont dépend la valeur des coefficients Bk sont les conditions de vol et de pilotage de l'aéronef ainsi que la phase de vol à l'instant correspondant à la relevée des données précitées. 1Z. Modèle conforme à l'une des revendications Il ou 12, caractérisé en ce que les coefficients Ai et Bk sont respectivement énumérés selon les variables Xij et les paramètres de vol dans les tableaux ci-après et en ce que le coefficient C vaut 0,18425416. 261564 1 PARAMETRESDE VOL COEFFICIENT ADDITIF Bk Estimation par le pilote CMI aux commandes -0,01128010 Estimation par le pilote CM1 non aux commandes -0,03695012 Estimation par le pilote CM2 aux commandes -0,02018463 Estimation par le pilote CM2 non aux commandes 0,00000000 Deux systèmes ECAM en panne (aucun système ECAM en fonctionnement) 0,03846630 Un système ECAM en panne (un système ECAM en fonctionnement) -0,00036457 Deux systèmes ECAM en fonctionnement, 0, 00000000 Remise des gaz ou descente rapide avec les deux FMS (système de gestion de vol) en fonctionnement 0,00000000 Préroulage avec les deux FMS en panne -0,01271110 Préroulage avec un FMS et un FMS en fonctionnement 0,22427695 Préroulage avec les deux FMS en fonctionnement 0,00000000 Roulage avant décollage avec les deux FMS en panne 0,01688276 Roulage avant décollage avec un FMS en panne et un FMS en fonctionnement 0, 07965894 Roulage avant décollage avec les deux FMS en fonctionnement 0, 00000000 Décollage avec les deux FMS en panne -0,01649307 Décollage avec un FMS en panne et un FMS en fonctionnement -0,02033955 Décollage avec les deux FMS en fonctionnement 0,00000000 Montée avec les deux FMS en panne -0,03001335 Montée avec un FMS en panne et un FMS en fonctionnement 0,01695141 Montée avec les deux FMS en fonctionnement 0,00000000 Croisière avec les deux FMS en panne 0,05614202 Croisière avec un FMS en panne et un FMS en fonctionnement -0,10278220 Croisière avec les deux FMS en fonctionnement 0,00000000 Descente avec les deux FMS en panne 0, 03471076 Descente avec un FMS en panne et un FMS en fonctionnement 0, 01862299 Descente avec les deux FMS en fonctionnement 0,00000000 261564 1 84. PARAMETRE DE VOL COEFFICIENT ADDITIF Bk Panne ou Urgence y compris Incapacité du pilote CM1 0,04277797 Liste de vérification pour Panne ou Urgence 0,06829986 Liste de vérification pendant un vol normal 0,03892804 Vol normal sans liste de vérification 0, 00000000 Remise des Gaz ou Descente d'Urgence 0,05458590 Pré-roulage - 0, 09502411 Roulage avant décollage - 0,07419961 Décollage 0,02750070 Montée 0,02952791 Croisière 0,01699097 Descente 0,00673707 Approche 0,03832675 Atterrissage 0,03482687 Roulage après atterrissage 0,00000000 Approche avec les deux FMS en panne -0,1405300 Approche avec un FMS en panne et un FMS en fonctionnement -0,04206698 Approche avec les deux FMS en fonctionnement 0,00000000 Atterrissage avec les deux FMS en panne -0, 07434618 Atterrissage avec les deux FMS en fonctionnement 0,00000000 Roulage après atterrissage avec les deux FMS en panne -0,11597344 Roulage après atterrissage avec un FMS en panne et un FMS en fonctionnement -0, 07514224 Roulage après atterrissage avec les deux FMS en fonctionnement 0, 00000000 Remise des gaz ou descente d'urgence avec automanette enclenchée 0,00000000 Préroulage avec automanette non enclenchée 0,18152502 Préroulage avec automanette enclenchée 0,00000000 Roulage avec automanette non enclenchée -0,10494458 Roulage avec automanette enclenchée 0,00000000 Décollage avec automanette non enclenchée -0, 3796761 Décollage avec automanette enclenchée 0,00000000 PARAMETRES DE VOL COEFFICIENT ADDITIF Bk Montée avec automanette non enclenchée -0,02778540 Montée avec automanette enclenchée 0,00000000 Croisière avec automanette non enclenchée 0,13777431 Croisière avec automanette enclenchée 0,00000000 Descente avec automanette non enclenchée 0,07989134 Descente avec automanette enclenchée 0,00000000 Approche avec automanette non enclenchée 0,04462338 Approche avec automanette enclenchée 0,0000000 Atterrissage avec automanette non enclenchée -0,01322323 Atterrissage avec automanette enclenchée 0,00000000 14. Dispositif pour la mise en oeuvre du procédé conforme à l'une des revendications 1 à 10, caractérisé en ce qu'il comprend des moyens (1) pour consigner les cotations de charge de travail des membres d'équipage. des smoens (2) pour interroger sélectivement les membres d'équipage quant à leur charge de travail et consigner simultanément la phase de vol et la condition de vol correspondantes, l'ensemble des moyens précités étant relié à des moyens (6, 7) de traitement et d'enregistrement des données fournies par ledit ensemble de moyens en synchronisme avec une base de temps de référence (8) et des moyens (18) de consignation et/ou d'enregistrement des données significatives des paramètres de vol en synchronisme avec ladite base de temps de référence. 15. Dispositif conforme à la revendication 14, caractérisé en ce que les moyens pour consigner les cotations de charge de travail des membres d'équipage comprennent un boîtier (1) distinct affecté à chaque pilote et comprenant des moyens pour visualiser une demande de cotation et des moyens de cotation suivant une échelle préétablie et en ce que les moyens pour interroger sélectivement les membres d'équipage comprennent un boîtier (2) comprenant des moyens (5) pour sélectivement adresser une demande de cotation vers l'un des membres d'équipage et des moyens (5) pour consigner les paramètres de vol correspondants. 16. Dispositif conforme à la revendication 14, caractérisé en ce qu'il comprend en outre des moyens (19, 19a. 19b) de mesure et d'enregistrement de l'électrocardiogramme des membres d'équipage (20), en synchronisme avec la base de temps de référence. 17. Dispositif conforme à l'une des revendications 14 ou 15, caractérisé en ce qu'il comprend en outre des moyens de prise de vues (11, 15) du poste de pilotage, des moyens de visualisation (10) et d'enregistrement (14) de ces vues et des moyens d'incrustation (13) sur les moyens de aguTwÀagp loa^ ap UoT;enqTs aun suep aquTwjaqpgàd ITeAed3 ap a6Aeqm es queInmlez ua a4olTd qTpnp semuewwo+Jad sael anle^9 uo allanbel suep 'eauomJe.p aqolTd un=p +auojee IaAnou un e uoT;e3depep qj6ap np TATns ne qa samuewjo+Jad sap alseuom ne 'luswauuoT3oa+Jad ne luoTiewJo+ el V eol e T suoitemTpua^aJ sap aunT1I aw5Jo+uom ippoad al uolas nuaqqo no i1 ç 11 suoTqemTpua^ai sap aunri V awJo+uom anbT;sT;eqs alpow np uoT;eTilddV oz -s;olTd S sas ap ITe^ejq ap aGieym el ins jauoJej[ ap anbTisTi-e4ei alla^nou aunp amuanluTil JaInWTS inod (pTX)+ = TA aqinom el asTIT3n uo allanbel suep e+auone Ia^nou unp quawaddolaA^p ne 401 V 1 suoT;emTpua^aa sap aun-, eq awJo+uon apgmod al UolaS nuaqqo no úI ç IT suoT;emTpua^a sacp aun; 0 awJo+uom anbTisTqe4s asepow np uoTeDTilddV '6- -gT+TzI4Àa quSwwpaD9lJd amuaiej}ej ap +auoige unp salTan ç saeedwoD as.insua que ITeAedà ap a5ieqp ap sànale^ sao eaquTwiaep aljIaq aun q quawaA;-elaJ IteAe; 3p a5beqm ap saquepuodsasoL Sl sJnale^ saI qTnpqp uo sallanbsap ÀiTàed e TA alqeTieA el ap saquepuodsasJon snale^A sl ainile5 uo QTX sSlqeTJe^ ssl ssAT;enTpuT quos;uop IoA ap sasJweied sel iuoTmemTTam ap sal6Qj sap sa^tqequasgida- Io^ ap saseqd;s suotqTpuoz saquaj++Tp suep aJnsaw uo 6+auoJge IaAnou np Tessap 01 IoA un suTow neap sJnoa nernb.a ua a sTJg;eJei 'auonJe IaAnou unp a6edTnb;! ap uoTiTsodwom el Sp UOT4emT+Tqjam el ? 601 t I suoT;entpua^aj sap aunI e awJo5uoz gppoJd al uolas nuaqqo no úI V TI suoT3emTpua^as sap aun.1 V awuonuo. anbTwsT;eqs asQpow np uoTeatIdd '8T S 6aOedTnbp saJqwaw sam ap ITe^ej4 ap a5jeqz ap suoiTeom sael aubTsuom inod (1) suaAow sap a as5edTnbpp sajqwaw sel quawa^.Toals jeaoiaqUT Ànod (Z) suaAow sap queue-w saguuop sap (0) uoTe5sTILensTi L 9S L9Z