Brevet FR2848307A1 — publié le 20040611

L'invention concerne un procédé et un dispositif pour le contrôle en pilotage d'un engin propulsé qui est muni de moyens de commande en orientation, notamment un lanceur muni de tuyères orientables. Rappels Pour le contrôle de vol des lanceurs, on distingue le guidage et le pilotage. Grossièrement on peut dire que: - le guidage donne les instructions permettant d'accomplir la mission (consignes d'attitudes). Ces instructions permettent de définir le mouvement du centre de gravité du lanceur pour que celui-ci atteigne une cible 15 prédéterminée (orbite de transfert généralement); - le pilotage suit les instructions fournies par le guidage en stabilisant le mouvement du lanceur autour de son centre de gravité et en s'adaptant aux variations de l'environnement extérieur (perturbations telles que vent, rafale, ...) et intérieur (dispersions et méconnaissances sur les 20 paramètres du lanceur, .. .). On peut affiner les deux propositions précédentes en disant que: A partir de la connaissance du point cible et de la position courante du lanceur mesurée par des capteurs inertiels, le guidage calcule la trajectoire que doit suivre le repère de consigne d'attitude du lanceur (son origine est fixée 25 sur le centre de gravité du lanceur et un de ses axes donne la direction que doit suivre l'axe de l'engin; les mesures inertielles du lanceur fournissent les attitudes du lanceur définies dans un repère supposé galiléen choisi fixe au sol). Lorsque le lanceur n'est soumis à aucune perturbation et que la 30 boucle de pilotage répond instantanément, l'axe de l'engin est confondu avec un des axes du repère de consigne. Le pilotage est une boucle d'asservissement ayant pour but premier de stabiliser les mouvements du lanceur autour de son centre de gravité (notion de stabilité), typiquement en ajustant l'orientation de la poussée. D'autre part, le pilotage réalise dans la mesure du possible la consigne d'attitude envoyée par 5 le guidage (le pilotage est un système asservi "suiveur") et assure une certaine immunité du lanceur face aux perturbations qu'il est susceptible de rencontrer (le pilotage est une boucle de régulation). Pour cela, la fonction pilotage, à partir des mesures inertielles d'attitudes (les mêmes que celles utilisées par le guidage), détermine l'orientation convenable de la poussée. Une perturbation telle que le vent perturbe le vol du lanceur en générant un mouvement autour de son centre de gravité. En réalité, le vent génère aussi un mouvement d'ensemble du lanceur (mouvement en translation ou dérive). Comme le guidage tient compte du point courant du lanceur (à l'aide des mesures inertielles) pour établir une nouvelle trajectoire jusqu'au point 15 cible, il n'est pas faux de considérer uniquement les mouvements autour du centre de gravité pour l'établissement des performances de la fonction de pilotage. La fonction pilotage doit tenir compte de certaines particularités dont 20 le fait, notamment, que par nature, le lanceur est un système non stationnaire. Ses caractéristiques physiques (position du centre de gravité, masse, inertie) et aérodynamiques évoluent en fonction du vol. Des phases transitoires vont venir interagir avec la fonction pilotage, parmi lesquelles la séparation d'étage, la présence de perturbations externes (vent) et internes (excentrement de la 25 poussée par exemple). Enfin, comme tout système physique, les paramètres du lanceur sont affectés d'incertitudes. La fonction pilotage est établie à partir d'un cahier des charges tenant compte explicitement des éléments précédents. Ce cahier des charges 30 est typiquement le suivant: - stabiliser le lanceur - Rattraper le lanceur et minimiser l'incidence a en phase de vol atmosphérique (voir figure 1). - Suivre la consigne d'attitude générée par le guidage - Rejeter les perturbations externes et internes - Tenir compte des performances de la servo-gouverne L'algorithme de pilotage est généralement intégré sous la forme d'équations récurrentes dans un calculateur numérique. L'algorithme de guidage, à partir des informations fournies par les détecteurs (attitudes, vitesses d'attitudes, accélérations) et de la navigation, calcule les consignes 10 d'attitude que doit réaliser la boucle de pilotage, et qui servent donc à l'algorithme de pilotage. En général, les lois de pilotage sont déterminées en supposant un découplage des axes tangage, lacet et roulis. Cette hypothèse permet de synthétiser la loi en utilisant un modèle un-axe valable aux petits angles. L'effet 15 du couplage des trois lois (une pour chaque axe) est vérifié a posteriori en simulation trois-axes. La figure 2 établit le schéma fonctionnel un-axe de la loi de pilotage. L'algorithme de pilotage est vu comme une boîte "noire" disposant d'informations provenant des détecteurs du lanceur (centrale inertielle pour la 20 mesure d'attitude et gyromètre pour la mesure de la vitesse d'attitude) et définissant une commande pour l'actionneur un-axe, ici une tuyère. La représentation un-axe suppose que le mouvement du lanceur est un mouvement plan. D'autre part, l'actionneur est asservi en position. Il rallie une position équivalente à la consigne de braquage traduite en élongation. 25 Etat de la technique La principale difficulté que rencontre le concepteur de l'ensemble de pilotage est la retranscription des objectifs physiques que doit satisfaire la loi de pilotage pour déduire un critère mathématique qui peut ensuite être optimisé. 30 Cette retranscription nécessite de contourner la nature physique des objectifs pour trouver une représentation mathématique adéquate. A titre d'exemple, la stabilité (d'un point de vue physique) est constatée par la qualité d'amortissement de la réponse temporelle de la boucle de pilotage soumise à une perturbation ou à une consigne. Mathématiquement, on peut aussi dire que le degré de stabilité (valeur de l'amortissement) est lié au 5 coefficient de surtension de la boucle de pilotage (résonance de la boucle de pilotage à une fréquence donnée). La relation explicite entre l'amortissement et le coefficient de surtension n'existe que pour des systèmes dits académiques (systèmes du second ordre) . Malgré ces restrictions, les synthèses actuelles de lois de pilotage 10 utilisent la notion de surtension plutôt que d'amortissement pour prendre en compte un objectif de stabilité car celle-ci s'adapte très bien aux critères établis à partir de la minimisation d'une norme mathématique. Comme la relation entre surtension et amortissement n'est pas directe, ce n'est qu'a posteriori par simulation temporelle, qu'on constate l'amortissement. Les lois de pilotage connues utilisent les algorithmes de synthèse suivants: - la méthode LQG (A4,M45,A5), - la méthode Hoo (A5). Chaque algorithme propose des paramètres de synthèse permettant 20 d'optimiser la loi pour que celle-ci réponde aux exigences et aux contraintes requises pour le pilotage (stabilité, performance en rattrapage, robustesse en stabilité et en performance). Les méthodes de synthèse précédentes cherchent à globaliser les différents objectifs en utilisant la notion de système standard. Un grand nombre 25 de problèmes d'analyse aussi bien que de synthèse de commande peuvent être écrits par un seul et même schéma standard sur lequel il est possible de formuler les objectifs de la lois sous la forme d'un critère mathématique dont on cherche à minimiser la norme (norme H2 ou Hco). La notion de problème standard est donc extrêmement puissante pour globaliser les objectifs sous la 30 forme d'un critère mathématique unique. Cependant, cette globalisation se paie en s'éloignant du sens physique que l'on peut donner au critère résultant. D'autre part, la mise sous forme standard engendre généralement une complexité de l'algorithme. Les paramètres de synthèse n'étant pas forcément adaptés à la résolution du problème posé, on est obligé de se donner des degrés de liberté supplémentaires (d'o une augmentation de l'ordre) pour retranscrire correctement les objectifs visés. Les méthodes modernes actuellement utilisées apportent un plus incontestable dans l'analyse et la synthèse des objectifs d'une loi de commande. Par contre, le critère résultant à partir duquel la synthèse de la loi sera effectuée, conduit à mettre en place une procédure "essaierreur" pour 10 converger vers un résultat satisfaisant. En effet, les solutions obtenues avec ces algorithmes ne permettent pas d'obtenir une relation entre le comportement en fréquence et le comportement temporel autrement que par simulation, en particulier le fort niveau de surtension de la boucle fermée corrélé avec un amortissement 15 convenable de cette dernière. On ne peut donc rechercher la loi optimale qu'en itérations successives "synthèse-simulation- vérification". Problème technique et présentation de l'invention L'invention a pour objet un procédé et un dispositif de pilotage d'un 20 engin orientable, asservi en attitude sur une trajectoire et soumis à des perturbations extérieures, mettant en oeuvre un algorithme de pilotage dont les paramètres peuvent être corrélés à des objectifs physiques fixés pour le pilotage. Elle propose à cet effet un procédé de pilotage d'un engin orientable 25 au moyen d'un actionneur, asservi en attitude sur une trajectoire et soumis à des perturbations extérieures selon lequel, au moins en tangage: - on définit au préalable une loi de pilotage comportant un terme correcteur du premier ordre, un terme de transfert comportant des paramètres caractéristiques de l'engin et des perturbations extérieures, un gain en attitude 30 et un gain en vitesse d'attitude; - on définit au préalable des relations de corrélation entre le gain en vitesse d'attitude et les paramètres du terme correcteur en fonction du gain en attitude, des paramètres du terme de transfert de l'engin et de deux termes représentant des consignes de rapidité et d'amortissement de la boucle de pilotage; et - pour au moins une phase de vol de l'engin, on détermine des 5 valeurs des paramètres du terme du transfert de l'engin et on choisit le gain en attitude et les termes de rapidité et d'amortissement, on en déduit les paramètres du terme correcteur et le gain en vitesse, et on applique à l'actionneur la loi de pilotage avec ces valeurs des paramètres. Elle propose de même un dispositif de pilotage adapté à la mise en 10 oeuvre du procédé tel que décrit ci-dessus, comportant une entrée recevant au moins une mesure en attitude de tangage, une sortie destinée à appliquer une commande d'un actionneur, au moins une boucle d'asservissement en tangage connectée entre cette entrée et cette sortie et comportant une filtre correcteur du premier ordre, un terme de transfert comportant des paramètres 15 caractéristiques de l'engin et des perturbations extérieures, un gain en attitude et un gain en vitesse d'attitude. On peut noter que le choix d'une loi de pilotage aussi simple va a priori à 'encontre de l'analyse actuelle de l'homme de métier selon laquelle un algorithme simple, donc avec peu de paramètres, est peu adapté à un bon 20 pilotage. Pourtant, contre toute attente, le choix d'une loi aussi simple est apprau permettre l'élaboration de lois permettant de corréler, compte-tenu de certaines simplifications de calcul par exemple par développement limité de la boucle fermée, les paramètres d'une telle loi à des paramètres physiques tels que l'amortissement (donc, la stabilité) et la rapidité de l'asservissement (donc 25 des performances temporelles), au gain en attitude autour de l'axe considéré, ainsi qu'à des paramètres physiques de la modélisation. Liste des figures Les objets, caractéristiques et avantages de l'invention seront décrits 30 en référence aux figures sur lesquelles: - la figure 1 est un schéma classique montrant le régime des contraintes sur un engin soumis à une pression dynamique, - la figure 2 est un schéma classique d'un système de pilotage unaxe, - la figure 3 est un schéma classique définissant divers angles, pour un axe donné, - la figure 4 est un schéma d'une architecture générique d'une loi de pilotage analytique utilisée selon l'invention, pour un axe donné, - la figure 5 en est une variante, dans le cas o la vitesse angulaire autour de l'axe considéré n'est pas mesurée, - la figure 6 est une courbe corrélant le gain et la phase de la boucle 10 ouverte dans le cas nominal et dans les cas dits BF et HF, - les figures 7 et 8 sont des courbes donnant la réponse en attitude, braquage et incidence de la boucle fermée de pilotage pour le BF (A6mn, KI max) et HF (A6max, K1min). Description détaillée Selon l'invention, on utilise l'algorithme de pilotage strictement nécessaire, en terme de complexité, pour assurer le pilotage d'un engin. L'algorithme repose sur une structure générique dont les paramètres sont calculés de manière formelle en fonction de ceux de l'engin à piloter et des 20 caractéristiques désirées du pilotage (suivi de trajectoire, minimisation de l'incidence, phases transitoires, ... ). La loi tient compte explicitement des incertitudes des paramètres de l'engin et d'un retard de chaîne globalisé (retard induit par la dynamique de l'actionneur + retard des capteurs + retard d à l'échantillonnage). Cette loi est applicable dans un cadre non stationnaire (évolution des paramètres du lanceur). La loi non stationnaire est obtenue rapidement car les paramètres de celle-ci sont formels. Compte tenu de la structure générique, il n'y a pas de problèmes associés aux techniques d'interpolation largement utilisées dans les pilotages des lanceurs. L'apport indéniable de cette méthode est que les paramètres de synthèse sont des paramètres physiques de haut niveau (amortissement de la boucle fermée, bande passante,). On sait associer à ces paramètres une grandeur physique et un objectif correspondants. Le réglage est donc obtenu directement à partir d'exigences physiques exprimées par le concepteur. Cela ressort des explications suivantes qui, par mesure de simplifications, se limitent au pilotage du mode rigide d'un engin. La présence 5 éventuelle des modes souples peut être prise en compte en choisissant une politique "prudente", c'est-à-dire en évitant d'exciter les modes souples par le gain de boucle de la loi de pilotage (par ajout d'un filtre atténuateur par exemple). L'établissement d'une loi de pilotage nécessite l'introduction d'un 10 modèle du lanceur appelé modèle de synthèse. La dynamique de boucle ouverte non corrigée peut être représentée par deux équations couplées décrivant les mouvements d'attitude et de dérive du lanceur (cf. fig. 3). L'établissement de la loi de pilotage nécessite le choix d'un modèle approprié de l'engin appelé modèle de synthèse. La dynamique de boucle ouverte non corrigée peut être représentée, en référence aux grandeurs définies à la figure 3, par deux équations couplées décrivant les mouvements d'attitude et de dérive de l'engin. Tout d'abord, en se limitant aux petits angles, le moment en tangage peut être défini par les équations (1) de l'annexe tandis que les forces 20 transverses par rapport à l'axe de la vitesse de l'air peuvent être définies par les équations (2). Par ailleurs, les divers angles de la figure 3 sont liés par l'équation (3) qui, après dérivation par rapport au temps, donne les équations (4). Après élimination de y entre les équations (2) et (4) on obtient le 25 système d'équations (5) o s est la variable de Laplace. Si on suppose que les termes identifiés dans la condition (C) sont petits pour les valeurs de ce ces contenus dans la bande passante de la boucle de pilotage, on peut approcher l'expression (5) par l'expression (6). Cette expression montre que le système est d'ordre 2 avec une 30 dynamique du premier ordre, instable à la pulsation A6.. Le système d'équations (6) constitue le modèle de synthèse à partir duquel vont être établies les bases de la synthèse analytique. Il en découle que le transfert de commande, c'est-à-dire le transfert entre la variable contrôlée (0) et la commande (braquage Pf) s'exprime par le terme K1/(s2-A6). Les comportements de la servo-gouverne et de la chaîne de mesure peuvent être modélisés par un retard pur global de chaîne, T, en série avec le transfert de commande, d'o l'expression (7), o G' représente, du point de vue automatique, une représentation mathématique du système physique que l'on 10 veut boucler. Le retard pur a, dans l'expression précitée, été approché par une approximation de PADE d'ordre 1. Le bouclage est réalisé en sorte que la boucle fermée correspondante répondre aux objectifs de stabilité et de performance attendus. 15 Exemples d'architectures génériques de la loi de pilotage analytique. La figure 4 montre l'architecture générique de la loi de pilotage analytique. Elle est présentée en supposant que les mesures en attitude mais aussi en vitesse d'attitude sont accessibles mais, comme cela apparaîtra plus 20 loin, cette hypothèse ne constitue pas une limitation de la méthode. Dans cette architecture, le retard T somme le retard de chaîne, le retard induit par l'échantillonage-blocage ainsi que le retard équivalent au déphasage de la gouverne. Ce retard peut aussi inclure le déphasage d'un filtre destiné à atténuer localement le gain de boucle de la loi de pilotage pour 25 assurer une robustesse en gain d'un mode souple lanceur. On impose un correcteur du premier ordre dans la chaîne directe et un gain en vitesse d'attitude (Kv) et un gain en attitude (Kp). L'ordre de s+ (02 la loi de pilotage est donc égal à 1. On évalue la première dynamique du procédé en boucle fermée pour 30 fixer cette dernière sur un modèle du 2ème ordre de paramètres (4, )nJ Cette politique permet de condenser le compromis performance-stabilité d'une manière élégante. L'amortissement induit un certain degré de stabilité tandis que C n (rapidité de l'asservissement) engendre des performances temporelles. Les paramètres (4,Cn) constituent donc les paramètres de réglage de " haut niveau ". On effectue ensuite un développement limité de la boucle fermée jusqu'à l'ordre 3. On obtient un système de trois équations à partir desquelles trois des quatre paramètres de réglage de la loi de pilotage (Kv,c"1,co2) sont déduits en fonction de (Kp, Con, 0) par les équations (8). Afin d'éviter une certaine lourdeur de l'exposé, les solutions 10 précédentes ont été établies sans tenir compte du retard global de la boucle z. En présence de celui-ci, ces dernières sont nettement plus complexes. Cependant, les solutions existent toujours et peuvent être aisément implantées au coeur d'un logiciel de synthèse de loi. Les remarques suivantes peuvent être formulées: - Le pôle du correcteur de la chaîne directe est instable (oe2 - Le zéro du correcteur compense la dynamique stable de l'engin (.= AX). - Pour A6 = 0, le pôle et le zéro du correcteur de la chaîne directe tendent vers 0. La stabilité du système est assurée par une avance de phase bien connue des techniques classiques de commande et définie par Kp + sK,. La détermination des paramètres de la loi montre que ceux-ci ont été calculés à partir de paramètres de haut niveau (amortissement et rapidité de la 30 boucle fermée). Cela nous permet d'établir le lien entre cet amortissement et la marge de stabilité de l'asservissement exprimée sous la forme d'une marge de il gain Mo (rappelons que la marge de gain est la valeur du gain de la boucle ouverte lorsque l'entrée et la sortie de cette dernière sont déphasées de 1800, la juste instabilité étant obtenue pour une valeur de la marge de gain = 1) donnée par l'équation (9). Nous observons que la marge de gain est affine en (amortissement de la boucle fermée). Notamment, elle est minimale pour un gain en attitude nul (Kp = 0). Sa valeur est donnée par l'expression (10). Si on choisit une dynamique de boucle fermée définie par con = A6 (dans ce cas, le lanceur est stabilisé), la marge de gain s'exprime par 10 l'expression (11). Cette expression précédente établit le lien direct entre l'amortissement et la marge de stabilité. La méthode précédente est applicable même si la mesure de vitesse d'attitude n'est pas accessible. Dans ce cas, on estime cette dernière par un 15 filtre estimateur de type passe-haut propre (existence d'un terme de transmission directe aux hautes fréquences) ou strictement propre (coupure du gain de filtre aux hautes fréquences avec une pente de 20 dB par décade au moins). Choisissons un filtre passe-haut avec une coupure de 20 dB par décade. Son expression est définie par l'équation (12). La loi de pilotage devient alors une loi mono bouclage dont l'architecture est définie par la figure 5. Comme pour le cas avec retour de vitesse d'attitude, les expressions de la loi de pilotage (Coic2,Kv) peuvent être entièrement déterminées par développement de la boucle fermée à l'ordre 3 et en fonction des paramètres 25 KP et T. Choix des paramètres de réglage Le suivi d'attitude engendre forcément une prise d'incidence pouvant être gênante du point de vue des efforts généraux appliqués sur le lanceur. En 30 conséquence, le choix du gain en attitude est fait de la manière suivante: Une valeur non nulle de Kp devra être fixée si on désire un bon suivi d'attitude (décollage, phase de vol hors atmosphère). Une valeur faible ou nulle de Kp devra être choisie pour minimiser l'incidence dans des phases de vol critiques comme celles correspondant à la pression dynamique maximale. La valeur de con détermine la rapidité de contre réaction de la boucle de pilotage. La pulsation oe,, doit être convenablement sélectionnée pour que le rattrapage du lanceur à la pression dynamique maximale soit compatible avec l'incidence admissible et les performances de la servogouverne. Le paramètre 10 de réglage own est donc étroitement lié avec l'objectif de performance de la loi de pilotage. Le choix de C n tiendra compte des incertitudes sur les paramètres du lanceur (objectif de robustesse en performance). Les objectifs de stabilité et de robustesse en stabilité sont fixés par le paramètre 4 de manière à ce que tous les lanceurs définis par les variations 15 des paramètres A6 et K1 soient stabilisés avec un amortissement au moins égal à l'amortissement spécifié par le concepteur. Evolution non stationnaire de la loi Compte tenu des différentes contraintes intervenant au cours du vol 20 (décollage, suivi de trajectoire, stabilité de l'engin aux phases critiques, précision d'injection en orbite (lanceur civil) ou de tir (missile), la loi la plus intéressante résulte d'une interpolation de plusieurs lois linéaires calculées à différents instants du vol. Les paramètres de la loi deviennent alors des fonctions du temps. La structure du correcteur étant générique, l'interpolation 25 de ces paramètres peut s'effectuer de manière directe, ceux-ci gardant le même sens physique d'un instant de vol à un autre. A un instant de vol t, on définit: - les paramètres du lanceur A6(t) et K, (t), - un gain en attitude Kp(t), - une rapidité de la boucle de pilotage co(t), - un amortissement 4 (t) qui peut être choisi constant pour tout le vol. Les paramètres de la loi correspondant à l'instant t sont évalués à partir des formules analytiques fournies par le développement limité de la boucle fermée. Le processus de synthèse est réitéré pour n points de vol de l'étage considéré. La loi globale non stationnaire est obtenue en interpolant les n lois stationnaires. Elargissement de la structure de la loi La loi analytique a deux formes génériques de base selon que l'on dispose ou non d'un retour de vitesse d'attitude, données par les équations (13) et (14). Ces formes peuvent être élargies bien que la nature de la solution 15 soit définie par l'une des deux formes de base. Une première proposition consiste à complexifier l'une ou l'autre des formes de bases en multipliant ces dernières par un filtre d'ordre quelconque. Ce filtre aura pour objectifs, soit de réduire localement en fréquences le gain pour éviter d'exciter des modes souples du lanceur à des fréquences particulières, soit d'apporter une 20 atténuation en gain de la boucle à partir d'une certaine fréquence définie par le concepteur. Bien entendu, il est tout à fait possible d'adjoindre un autre type de filtre n'ayant aucune fonction particulière si ce n'est de complexifier apparemment l'une des deux formes de base. 25 Sens physique et valeurs typiques des paramètres Compte tenu de l'expression analytique des paramètres KV,eo1, eteo2, il est possible d'associer pour certains d'entre eux un sens physique donné: La pulsation ew, est égale à Agi, dynamique propre du lanceur. On 30 réalise une compensation de cette dernière. Plus précisément, on peut écrire que la pulsation w1 est proportionnelle à la pression dynamique de l'instant de vol considéré et inversement proportionnel à l'inertie du lanceur. Sa valeur est a priori difficile à quantifier, puisque celle-ci dépend des caractéristiques du lanceur. L'expression de la pulsation od2 est fonction de la rapidité optée pour la boucle de pilotage. De la même manière que précédemment, il est difficile de 5 donner une plage de valeurs numériques. Cependant qualitativement, on peut préciser que plus la valeur de oe2 sera grande, plus la boucle de pilotage sera rapide. En pratique, la rapidité de la boucle est choisie telle que cette dernière ne force pas le procédé (lanceur) au-delà de sa dynamique propre (XA6). L'expression du gain Kv est relativement complexe. Outre les termes 10 Kl, A6, Kp et cor, il apparaît le terme d'amortissement , paramètre fixant la stabilité et la robustesse en stabilité de la boucle de pilotage. Qualitativement, le degré de stabilité de l'asservissement est d'autant plus important que la valeur numérique du gain KV est élevée, celui-ci étant proportionnel à . Exemple L'exemple ci-après montre l'application de la synthèse analytique d'une loi de pilotage stationnaire pour un lanceur X. Les valeurs des paramètres A6 et Ki à la pression dynamique maximale (correspondant à une trajectoire fixée) sont données par: A6=35 15s-2 KI = 30 1os-2 Le cas majorant pour la stabilité du lanceur correspond à A6=A6max=50 s-2 Ki =KI min=20 s-2 La loi de pilotage est synthétisée à la pression dynamique maximale avec les paramètres du lanceur ci-dessus. On fixe l'amortissement de la dynamique de boucle fermée à: 4=0.3 (c'est l'amortissement requis correspondant au cas majorant) Dans le but de minimiser l'incidence à la pression dynamique 30 maximale, on fixe un gain en position nul soit: Kp=0 Le choix du paramètre con (paramètre fixant la rapidité de l'asservissement) est déterminé en vérifiant que sa valeur permet de ne pas dépasser la vitesse et l'excursion de braquage atteignables par la gouverne. On suppose que la valeur de on est connue et est donnée par: con=3.5rad/s Si tel n'est pas le cas on réitère l'étape de synthèse du correcteur jusqu'à trouver unevaleur de con compatible avec les performances de la gouverne vérifiées par simulation temporelle. Les paramètres du correcteur sont fournis par les équations à la base de l'invention. Nous avons (pour Kp=0) col = A6max = 7 rad/s ______-___ V (3. 5)2 Co - m 6m = - 5) - 1.74 rad / s A6max + KIKp 50 K = 2onA46max + A6max (2n +A6max) 2 * 0.3 * 3.5 * (352 50) .2 Kv = 2 *.33550Jb(.22.23 sKlnmin on2 20*3.52 L'application des formules suppose que la chaîne de pilotage n'induit pas de retard. Si cette hypothèse n'est pas vérifiée, il est nécessaire de réévaluer les expressions analytiques des paramètres du correcteur en tenant compte de la présence du retard de chaîne. La marge de gain basse fréquence est donnée par l'expression 20 analytique M0 =1+ 2 1+ 2 * 0.- 3 1.25soit 2 dBavecx A6x - =x+l/x 2+1/2 o 3.5 Une simulation temporelle montre que l'on obtient un braquage maximal de la gouverne de 5.6 et une vitesse maximale de braquage de 30 . Ces chiffres sont obtenus pour le cas majorant du point de vue de la stabilité 25 (cas A6max,Kl min dit cas BF). Les performances fréquentielles sont données par le figure 6. En particulier, on retrouve, par simulation, la marge de gain calculée dans le cas majorant (A6max,Kl min). Les performances temporelles obtenues avec le réglage calculé sont données par les courbes des figures 7 et 8 " La loi de pilotage analytique permet d'atteindre rapidement les objectifs de stabilité et de performance pour l'ensemble des lanceurs décrits par 5 les variations des paramètres A6 et KM. Contrairement aux méthodes du type LQG ou Hoo, l'itération fréquentielle, qui consiste à vérifier a posteriori les marges de stabilité, a été éliminée. L'itération temporelle, qui permet de choisir la rapidité de la boucle de pilotage, reste nécessaire car l'objectif de performance est directement lié à la nature des perturbations agissant sur le 10 lanceur (vent, rafale, perturbations internes, ...). L'extension non stationnaire est directe et ne pose pas de problèmes particuliers du point de vue de la synthèse. La loi de pilotage analytique peut éventuellement être étendue pour des engins autres que les lanceurs. L'étude de base montre qu'elle est élaborée 15 à partir d'un modèle de synthèse déduit des équations de la mécanique du vol, l'hypothèse forte étant l'instabilité de l'engin. Il n'y a donc aucune limitation pour que ce type de loi puisse être applicable pour des avions instables sur un ou plusieurs axes de vol (tangage, lacet et roulis). Pour cela, il suffit d'adapter les équations de la mécanique du vol aux avions instables. Cette adaptation 20 permettra d'obtenir un nouveau modèle de synthèse duquel on déduira les paramètres de la loi de pilotage analytique. De manière préférée: - on détermine pour la phase de vol des valeurs (A6max, Kl min) des paramètres du terme de transfert qui majorent les perturbations extérieures. - on applique un test de validation par simulation aux choix du gain en attitude et des termes de rapidité et d'amortissement pendant toute la durée du vol et on applique ces choix pendant tout le vol si ces choix sont validés, sinon on décompose le temps de vol en plusieurs phases. - le terme de transfert est défini pour des perturbations extérieures 30 essentiellement causées par l'air atmosphérique. - on applique en roulis et en lacet des lois de pilotage définies selon la même méthode que la loi de pilotage en tangage. ANNEXE: FORMULES fô=q = =q :I = Lf (QSC( a) + LtPcfJ q = A6a + K,8 (1) QSC,, + Pc Pc mVpQSC,,.a+Pcsin(f +a)z=Ä a+ -f (2) mV mV avec m: Q: Cza I: S: Pc: Lf: Lt: masse du lanceur pression dynamique gradient de portance inertie en tangage surface de réference de l'engin pour le calcul de la force de portance Poussée de l'engin distance centre de gravité - foyer aérodynamique distance centre de gravité - centre de rotation tuyère (9=r +a + arctan. r -+a+ +W Vir aVir (3) 0?;+: à+g (4) t 1 dw 1 dw dz Vzdw ov, g.... (gradient de vent) V t dt V dz dtl V dz vria (La vitesse Vz est la projection de la vitesse du lanceur sur l'axe vertical) e a i_ I__ _ _ _ aO (2 + QSCza, + PC mV 1 QSCa.K, + (KI -A6)Pc + s mV K -Pc K ---s mV (5) Pc et QSCzaKt + (K - A6)Pc KmV retv K1mV M Vc (C) a) ((s2 - A6 V V (6) O, (1-:s/2) KI f (1+Ts/2) s2-A6 (7) K_=:2 2ó0 (2 6 + KlK 'p) + 6((O n + A6) + A6 KIlKp" V 2 r- -lK con (8) C0o = \A6 4A6 Co2 (0 2 -- A6 + K1Kp 2 e( con + 'z KIK p) -A6 >1 Cn +A6 +KKp MAo =: 1 + (9) 214on A60 __2_ _A60 nw 1+ __Mo(Kp =OO=A)=1+ A6o +(Op= n X + 1 /X O, (10) (11) Af Oz=zO S (12) (S2 /wc2 +s/wc+ 1) f, = (KpO + KvO) (- ) (loi avec retours de vitesse à et d'attitude 8) (s + 02,) (13) (s2 /w0 (Kpc+ slv) (s+ 'I) (loi sans retour de vitesse O) (s5 /WC' + s/wC + 1) (S+ 2) (14) REVENDICATIONS 1. Procédé de pilotage d'un engin orientable au moyen d'un actionneur, asservi en attitude sur une trajectoire et soumis à des perturbations extérieures selon lequel, au moins en tangage: - on définit au préalable une loi de pilotage comportant un terme correcteur du premier ordre, un terme de transfert comportant des paramètres caractéristiques de l'engin et des perturbations extérieures (A6, K1), un gain en attitude (Kp) et un gain en vitesse d'attitude (Kv); - on définit au préalable des relations de corrélation entre le gain en vitesse d'attitude et les paramètres du terme correcteur en fonction du gain en attitude, des paramètres du terme de transfert de l'engin et de deux termes représentant des consignes de rapidité (O)n) et d'amortissement ( l), et - pour au moins une phase de vol de l'engin, on détermine des 15 valeurs des paramètres du terme de transfert et on choisit le gain en attitude et les termes de rapidité et d'amortissement, on en déduit les paramètres du terme correcteur et le gain en vitesse, et on applique à l'actionneur la loi de pilotage avec ces valeurs des paramètres. 2. Procédé selon la revendication 1, caractérisé en ce qu'on définit 20 les lois de corrélation à partir d'un développement limité de l'équation de la boucle fermée constituée par l'engin et son actionneur. 3. Procédé selon la revendication 1 ou 2, caractérisé en ce qu'on détermine pour la phase de vol des valeurs (A6max, Kl min) des paramètres du terme de transfert qui majorent les perturbations extérieures. 4. Procédé selon l'une quelconque des revendications 1 à 3, caractérisé en ce qu'on prévoit en outre dans la loi de pilotage un filtre. 5. Procédé selon l'une quelconque des revendications 1 à 4, caractérisé en ce qu'on applique la loi de pilotage en utilisant des mesures de l'attitude et de la vitesse en attitude. 6. Procédé selon l'une quelconque des revendications 1 à 4, caractérisé en ce qu'on applique la loi de pilotage en utilisant un filtre estimateur de vitesse d'attitude dans la loi de pilotage. 7. Procédé selon l'une quelconque des revendications 1 à 6, caractérisé en ce qu'on applique un test de validation par simulation aux choix du gain en attitude et des termes de rapidité et d'amortissement pendant toute la durée du vol et on applique ces choix pendant tout le vol si ces choix sont validés, sinon on décompose le temps de vol en plusieurs phases. 8. Procédé selon l'une quelconque des revendications 1 à 7, caractérisé en ce que le terme de transfert est défini pour des perturbations extérieures essentiellement causées par l'air atmosphérique. 9. Procédé selon l'une quelconque des revendications 1 à 8, 10 caractérisé en ce que la loi de pilotage comporte en outre un terme de retard. 10. Procédé selon l'une quelconque des revendications 1 à 9, caractérisé en ce qu'on applique en roulis et en lacet des lois de pilotage définies selon la même méthode que la loi de pilotage en tangage. 11. Dispositif de pilotage adapté à la mise en oeuvre du procédé 15 selon l'une quelconque des revendications 1 à 10 comportant une entrée recevant au moins une mesure en attitude de tangage, une sortie destinée à appliquer une commande d'un actionneur, au moins une boucle d'asservissement en tangage connectée entre cette entrée et cette sortie et comportant une filtre correcteur du premier ordre, un terme de transfert 20 comportant des paramètres caractéristiques de l'engin et des perturbations extérieures, un gain en attitude et un gain en vitesse d'attitude.