FR3059777A1 - Procede d'essai - Google Patents
Procede d'essai Download PDFInfo
- Publication number
- FR3059777A1 FR3059777A1 FR1662001A FR1662001A FR3059777A1 FR 3059777 A1 FR3059777 A1 FR 3059777A1 FR 1662001 A FR1662001 A FR 1662001A FR 1662001 A FR1662001 A FR 1662001A FR 3059777 A1 FR3059777 A1 FR 3059777A1
- Authority
- FR
- France
- Prior art keywords
- test method
- turbine
- linear
- kalman filter
- physical system
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000010998 test method Methods 0.000 title claims abstract description 30
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 21
- 238000005259 measurement Methods 0.000 claims abstract description 19
- 238000005070 sampling Methods 0.000 claims abstract description 18
- 230000006978 adaptation Effects 0.000 claims abstract description 8
- 238000012360 testing method Methods 0.000 claims description 27
- 239000011159 matrix material Substances 0.000 claims description 20
- 238000000034 method Methods 0.000 claims description 8
- 230000009965 odorless effect Effects 0.000 claims description 8
- 238000011144 upstream manufacturing Methods 0.000 claims description 4
- 238000009530 blood pressure measurement Methods 0.000 claims description 2
- 238000009529 body temperature measurement Methods 0.000 claims description 2
- 238000004590 computer program Methods 0.000 claims description 2
- 230000007704 transition Effects 0.000 claims description 2
- 230000037213 diet Effects 0.000 claims 1
- 235000005911 diet Nutrition 0.000 claims 1
- 239000012530 fluid Substances 0.000 description 9
- 238000012545 processing Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 4
- 238000001914 filtration Methods 0.000 description 4
- 230000015654 memory Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 239000007788 liquid Substances 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 239000003380 propellant Substances 0.000 description 3
- 244000124209 Crocus sativus Species 0.000 description 2
- 230000006399 behavior Effects 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000009916 joint effect Effects 0.000 description 1
- 230000007257 malfunction Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005312 nonlinear dynamic Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000135 prohibitive effect Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M15/00—Testing of engines
- G01M15/14—Testing gas-turbine engines or jet-propulsion engines
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/048—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators using a predictor
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B17/00—Systems involving the use of models or simulators of said systems
- G05B17/02—Systems involving the use of models or simulators of said systems electric
Landscapes
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Artificial Intelligence (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Combustion & Propulsion (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
L'invention concerne le domaine des procédés d'essai et spécifiquement un procédé d'essai d'un système physique non-linéaire (10) comprenant les étapes suivantes : effectuer une succession d'échantillonnages de valeurs de commande et de valeurs de mesure dudit système physique non-linéaire ; soumettre lesdites valeurs de commande et de mesure à un filtre de Kalman discret non-linéaire utilisant un modèle numérique dudit système non-linéaire ; et mettre en œuvre un algorithme d'adaptation paramétrique pour optimiser au moins un premier et un deuxième paramètre dudit modèle numérique sur base de valeurs issues dudit filtre de Kalman discret non-linéaire .
Description
® RÉPUBLIQUE FRANÇAISE
INSTITUT NATIONAL DE LA PROPRIÉTÉ INDUSTRIELLE © N° de publication :
(à n’utiliser que pour les commandes de reproduction)
©) N° d’enregistrement national
059 777
62001
COURBEVOIE ©IntCI8: G 01 M 15/14 (2017.01)
DEMANDE DE BREVET D'INVENTION
A1
©) Date de dépôt : 06.12.16. | © Demandeur(s) : SAFRAN AIRCRAFT ENGINES — |
(30) Priorité : | FR. |
@ Inventeur(s) : LE GONIDEC SERGE, DANIEL et | |
TREANGLE CLEMENT. | |
(43) Date de mise à la disposition du public de la | |
demande : 08.06.18 Bulletin 18/23. | |
©) Liste des documents cités dans le rapport de | |
recherche préliminaire : Se reporter à la fin du | |
présent fascicule | |
(© Références à d’autres documents nationaux | ® Titulaire(s) : SAFRAN AIRCRAFT ENGINES. |
apparentés : | |
©) Demande(s) d’extension : | © Mandataire(s) : CABINET BEAU DE LOMENIE. |
Pty PROCEDE D'ESSAI.
FR 3 059 777 - A1 (P/ L'invention concerne le domaine des procédés d'essai et spécifiquement un procédé d'essai d'un système physique non-linéaire (10) comprenant les étapes suivantes:
effectuer une succession d'échantillonnages de valeurs de commande et de valeurs de mesure dudit système physique non-linéaire;
soumettre lesdites valeurs de commande et de mesure à un filtre de Kalman discret non-linéaire utilisant un modèle numérique dudit système non-linéaire; et mettre en oeuvre un algorithme d'adaptation paramétrique pour optimiser au moins un premier et un deuxième paramètre dudit modèle numérique sur base de valeurs issues dudit filtre de Kalman discret non-linéaire .
Tub
Tdt
Tdb
Arrière-plan de l'invention
La présente invention concerne un procédé d'essai d'un système physique non-linéaire.
L'objectif principal des essais de systèmes physiques est celui de caractériser le comportement de ces systèmes physiques, et en particulier identifier des valeurs de certains paramètres liés à leur fonctionnement.
Toutefois, l'identification de ces valeurs devient plus difficile quand elles correspondent à plusieurs paramètres différents liés entre eux. Ainsi par exemple, dans un système physique comprenant une turbine et une vanne de contournement de la turbine, comme notamment dans un circuit d'alimentation de moteur-fusée à ergols liquides, on peut souhaiter identifier un rendement de la turbine, sous forme d'un rapport entre une performance réelle et une performance théorique pendant l'essai, et un décalage entre une consigne de position de la vanne et une position réelle. Avec les procédés d'essai de l'état de la technique, séparer les effets de l'un et de l'autre sur le fonctionnement global d'un système physique non-linéaire exige une multiplication des régimes de fonctionnement différents testés dans chaque essai, voire une multiplication des essais. Le bruit de mesure et de traitement peut rendre cette identification encore plus ardue.
Objet et résumé de l’invention
La présente divulgation concerne un procédé d'essai visant à remédier à ces inconvénients, et plus particulièrement un procédé d'essai d'un système physique non-linéaire qui permette d'obtenir de manière plus rapide et facile des valeurs pertinentes d'au moins deux paramètres liés au fonctionnement du système physique pendant l'essai, et notamment leur identification dans une seule démarche d'essai et de calcul.
Pour cela, ce procédé comprend les étapes d'effectuer une succession d'échantillonnages de valeurs de commande et de valeurs de mesure dudit système physique non-linéaire, soumettre lesdites valeurs de commande et de mesure à un filtre de Kalman discret non-linéaire utilisant un modèle numérique dudit système non-linéaire, et mettre en œuvre un algorithme d'adaptation paramétrique pour optimiser au moins un premier et un deuxième paramètre dudit modèle numérique sur base de valeurs issues dudit filtre de Kalman discret non-linéaire.
Les filtres de Kalman pouvant non seulement fournir des estimations de variables réelles à partir de données de mesure et de commande bruitées, mais aussi des valeurs indicatives de la précision de ces estimations, telles que des variances et covariances de ces estimations, il est ainsi possible d'utiliser ces valeurs issues du filtre de Kalman discret non-linéaire dans un algorithme d'adaptation paramétrique pour créer une boucle de retour supplémentaire permettant d'ajuster les paramètres du modèle numérique utilisé dans le filtre de Kalman discret non-linéaire pour améliorer la précision des estimations.
Le filtre de Kalman discret non-linéaire peut notamment être un filtre de Kalman non-linéaire du type dit « sans-odeur », qui permet d'obtenir des résultats particulièrement précis pour un grand nombre de systèmes non-linéaires. Toutefois, d'autres types de filtres de Kalman discrets nonlinéaires, comme par exemple le filtre de Kalman étendu, peuvent être utilisés alternativement à un filtre de Kalman « sans-odeur ».
L'algorithme d'adaptation paramétrique peut notamment être un algorithme d'optimisation desdits aux moins deux paramètres du modèle numérique par minimisation d'une valeur calculée sur base d'au moins une matrice de covariance issue dudit filtre de Kalman discret non linéaire. Les matrices de covariance issues de filtres de Kalman discrets sont une source utile d'information sur la précision des estimations issues du filtre et, en conséquence, sur la qualité du modèle numérique sous-jacent. L'au moins une matrice de covariance peut par exemple comprendre une matrice de covariance d'innovation réelle dudit filtre de Kalman discret non-linéaire, et la valeur à minimiser peut comprendre la somme de tous les éléments d'au moins une matrice de covariance issue dudit filtre de Kalman non linéaire, voire la somme de tous les éléments d'une succession de matrices de covariances issue dudit filtre de Kalman discret non linéaire pendant une séquence d'échantillonnage.
Le système physique non-linéaire soumis à l'essai peut notamment comprendre une turbine et une vanne de contournement de ladite turbine. Ainsi, ce procédé d'essai permettra l'identification simultanée de paramètres de fonctionnement de la turbine et de la vanne de contournement. Par exemple, le premier paramètre du modèle numérique du système physique non-linéaire pourrait correspondre à un rapport entre une performance théorique et une performance réelle de la turbine pendant l'essai, et notamment à un rapport entre un couple théorique et un couple réel généré par la turbine pendant l'essai, et le deuxième paramètre du modèle numérique du système physique non-linéaire pourrait correspondre à un décalage entre une consigne de position et une position réelle de la vanne pendant l'essai. Lesdites valeurs de commande peuvent notamment comprendre une consigne d'ouverture de la vanne de contournement, une vitesse mesurée de rotation de la turbine, et/ou un couple mesuré de la turbine, ce dernier pouvant être mesuré, par exemple, à travers de mesures de pression et débit d'une pompe actionnée par la turbine. Lesdites valeurs de mesure peuvent notamment comprendre des mesures de pression et de température en amont et en aval de la turbine et de la vanne de contournement.
Ce procédé d'essai est particulièrement efficace pour obtenir des informations utiles à partir de fonctionnements non-stationnaires du système physique objet de l'essai. En conséquence, au moins une partie desdits échantillonnages successifs peut être effectuée pendant une période de transition d'un premier régime de fonctionnement à un deuxième régime de fonctionnement dudit système physique non-linéaire différent du premier régime.
Ce procédé d'essai peut notamment être mis en oeuvre à l'aide d'un ordinateur. Ainsi, cette divulgation concerne aussi un programme d'ordinateur configuré pour la mise en œuvre de ce procédé d'essai et un support informatique lisible par un dispositif informatique et comprenant des instructions pour la mise en œuvre de ce procédé d'essai par un dispositif informatique lisant ledit support. Par « support informatique » on entend tout support permettant le stockage de données, de manière durable et/ou transitoire, et leur subséquente lecture par un système informatique. Ainsi, on entend par « support informatique », entre autres, les bandes magnétiques, les disques magnétiques et/ou optiques, ou les mémoires électroniques d'état solide, tant volatiles que non-volatiles.
Brève description des dessins
L’invention sera bien comprise et ses avantages apparaîtront mieux, à la lecture de la description détaillée qui suit, d’un mode de réalisation représenté à titre d’exemple non limitatif. La description se réfère aux dessins annexés sur lesquels :
- la figure 1 est un dessin schématique illustrant un ensemble de turbine et vanne de contournement formant un système physique non-linéaire à mettre à l'essai ;
- la figure 2 est un schéma illustrant le traitement des données dans ce procédé d'essai ; et
- les figures 3A à 3D représentent des étapes successives de traitement de données suivant le schéma de la figure 2.
Description détaillée de l'invention
La figure 1 illustre un exemple de système physique non-linéaire 10 pouvant être mis à l'essai suivant un mode de réalisation spécifique du procédé d'essai objet de la présente divulgation. Ce système physique non-linéaire 10 comprend une turbine 11 et une vanne 12 de contournement de la turbine dans un conduit 13' de contournement de la turbine dérivé d'un conduit principal 13 de manière à ce que la turbine 11 et la vanne 12 puissent être alimentées en fluide moteur à travers une entrée commune ΓΓ. Le régime de la turbine 11 peut ainsi être en partie réglé par le degré d'ouverture de la vanne 12, qui peut ainsi permettre le contournement de la turbine à une proportion plus ou moins grande de fluide moteur. En outre, dans l'exemple illustré, la turbine 11 est couplée en rotation à une pompe 14, de manière à ce que la rotation de la turbine 11, actionnée par la détente du débit de fluide moteur ne contournant pas la turbine 11 à travers la vanne 12, actionne à son tour la pompe 14. La pompe 14 présente sa propre entrée IP de fluide. La turbine 11 et la pompe 14 forment donc ensemble une turbopompe telle qu'utilisée, par exemple, dans un circuit d'alimentation d'un moteur-fusée en au moins un ergol liquide, et notamment un ergol cryogénique tel que l'hydrogène liquide.
Avant de soumettre un tel système physique non-iinéaire 10 à un usage réel, il est typiquement nécessaire d'effectuer des essais afin de caractériser ses composants, c'est-à-dire déterminer des paramètres associés à chacun de ces composants dans un modèle numérique permettant de simuler le fonctionnement du système physique nonlinéaire 10 dans chacune des différentes situations pouvant être rencontrées dans un fonctionnement réel, pour ainsi optimiser sa commande et éviter les dysfonctionnements. Dans le cas illustré, un premier paramètre rp à déterminer correspond à un rapport entre un couple Mt théorique et un couple réel Mr généré par la turbine 11, et un deuxième paramètre ov correspond à un décalage entre une consigne de position et une position réelle de la vanne 12. Comme ces deux paramètres ont un effet conjoint sur le fonctionnement de ce système physique non-linéaire 10, il est difficile de déterminer ces deux paramètres distincts, d'autant plus que les signaux évalués pour cette détermination sont bruités. Dans une approche appliquant un modèle stationnaire, il peut ainsi être nécessaire d'effectuer des essais couvrant un grand nombre de points de fonctionnement, en maintenant chaque point de fonctionnement pendant une durée de temps suffisant. En conséquence, la durée et le coût des essais peuvent augmenter rapidement et devenir rédhibitoires, surtout si ce système physique non-linéaire 10 n'est qu'une partie d'un système plus complexe comprenant un nombre encore plus grand de composants avec des paramètres à identifier.
Afin d'éviter ces inconvénients, cette divulgation propose un procédé d'essai permettant de déterminer ces paramètres dans une seule démarche de calcul, sans imposer des points de fonctionnement stationnaires, et cela à partir de signaux soumis à du bruit, bruit provenant des instruments de mesure, de la conversion des signaux analogues continus en des signaux numériques discrets et du système physique non-linéaire 10 lui-même.
Pour procéder à l'essai, le système physique non-linéaire 10 illustré est instrumenté avec un capteur de pression 15 et un capteur de température 16 en amont de la turbine 11, un capteur de pression 17 et un capteur de température 18 en amont de la vanne 12, un capteur de pression 19 et un capteur de température 20 en aval de la turbine 11, un capteur de pression 21 et un capteur de température 22 en aval de la vanne 12, un capteur de pression 23 et un capteur de débit 24 en aval de la pompe 14 et un capteur de vitesse 25 sur l'arbre rotatif de la turbine
11. Par ailleurs, la vanne 12 comprend un actionneur avec encodage de position. Les capteurs 15 à 25 et l'actionneur de la vanne 12 sont tous connectés à un calculateur numérique 26, pouvant notamment être un ordinateur programmé pour mettre en œuvre le procédé d'essai.
Pendant l'essai, les entrées de fluide ΓΤ et IP sont connectées à une source de fluide qui les alimente en fluide. Afin de faire varier le régime de fonctionnement du système physique non-linéaire 10 pendant l'essai, il est possible de varier les pressions et températures des débits de fluide alimentant chaque entrée de fluide ΓΤ, IP, ainsi qu'une consigne de position cv de la vanne 12. Pendant l'essai, la vanne 12 transmet un signal représentant cette consigne de position cv au dispositif informatique 26. Simultanément, les capteurs de pression 15, 17, 19, 21 et 23 transmettent, respectivement, des signaux put, pUb, Pat, Pdb, Pdp correspondant chacun à la pression captée par le capteur correspondant, les capteurs de température 16, 18, 20, 22 transmettent, respectivement, des signaux Tut, TUb, Tdt, Tdb correspondant chacun à la température captée par le capteur correspondant, le capteur de débit 24 transmet un signal Vp correspondant au débit capté par ce capteur 24, et le capteur de vitesse 25 transmet un signal vpt correspondant à la vitesse de rotation de l'arbre rotatif de la turbine 11.
Tous ces signaux sont reçus par le calculateur numérique 26, qui peut prendre la forme d'un ordinateur programmable avec une unité centrale de traitement de données, une unité de mémoire, pouvant comprendre une mémoire morte et une mémoire volatile, et des interfaces d'entrée et de sortie. Le calculateur numérique 26 ne pouvant opérer qu'un traitement discret de ces signaux, et non pas continu, il procède à un échantillonnage successif des signaux dans le temps, avec une fréquence d'échantillonnage déterminée et suffisamment élevée pour obtenir une approximation satisfaisante d'un traitement continu des signaux.
La figure 2 montre un schéma fonctionnel du traitement de ces signaux dans le calculateur numérique 26. Ainsi, les signaux de pression Put, Pub, Pdt, Pdb, PdP, et de température Tut, Tub, Tdt, Tdb sont échantillonnées à des moments t successifs par le convertisseur analogique-numérique 50. Ces données échantillonnées de pression et de température, réunies pour chaque moment t, forment ensemble un état mesuré yt. En même temps, les signaux de vitesse vpt, débit Vp et position de la vanne cv sont échantillonnés aux mêmes moments t successifs par le convertisseur analogique-numérique 50 et les données échantillonnées correspondantes forment ensemble une commande ut pour chaque moment t. Un signal correspondant à la mesure, directe ou indirecte, d'un couple généré par la turbine 11 peut éventuellement aussi être intégré dans cette commande ut. Ces états mesurés yt et commandes ut alimentent un observateur d'états 27.
Cet observateur d'états TJ est à filtre de Kalman « sans odeur » ou UKF (pour son acronyme anglais, « Unscented Kalman Filter »). Cet algorithme de filtrage, décrit dans « The Unscented Kalman Filter for Nonlinear Estimation », Proceedings of Symposium 2000 on Adaptive Systems for Signal Processing, Communication and Control (AS-SPCC), IEEE, Lake Louise, Alberta, Canada, Oct, 2000, peut être utilisé non seulement pour filtrer le bruit de signaux, mais même de manière prédictive pour prévoir l'évolution à court terme d'un système dynamique non-linéaire comme système physique non-linéaire 10.
Dans un tel système dynamique, on peut assumer qu'il existe une séquence de Markov d'états latents xt évoluant dans, le temps suivant une fonction de dynamique modèle A comprenant comme paramètres, pour le système physique non-linéaire 10, le premier paramètre rp et le deuxième paramètre ov.
Ces états latents sont observés indirectement par les capteurs indiquant des états mesurés yt à travers une fonction de mesure G. Ainsi, xt et yt peuvent être exprimés suivant les formules :
xt = F(xt_i,ut)+ ε y, =G(Xl)+v
Les valeurs ε [epsilon] et v [nu] représentent, respectivement, le bruit inhérent au système et le bruit de mesure, et suivent toutes les deux des distributions gaussiennes.
L'objet d'un algorithme de filtrage est celui d'inférer l'état du système dynamique à travers les valeurs bruitées mesurées par les capteurs. Un filtre de Kalman offre une inférence rapide et précise pour les systèmes linéaires. Il n'est toutefois pas directement applicable aux systèmes nonlinéaires, dans lesquels le présent cas d'application pourrait potentiellement être classé. Parmi différentes alternatives pour adapter l'algorithme de filtrage de Kalman aux systèmes non-linéaires, on connaît notamment le filtre de Kalman « sans odeur ». Cet algorithme propage plusieurs estimations de l'état latent xt à travers les fonctions A et G et reconstruit une distribution gaussienne comme si les valeurs propagées provenaient d'un système linéaire. Les positions des estimations de l'état latent xt sont nommées « points sigma », et sont calculées à partir d'une moyenne et une variance initiales avec un schéma d'approximation appelé transformation non-parfumée.
Ainsi, dans l'observateur d'état 27, après une initialisation 28 fournissant une estimation initiale d'état x0 et une matrice de covariance correspondante Po, des points sigma initiaux Xo°, Xo1, Xo2, Xo3, Xo4 illustrés sur la figure 3A sont calculés, par une telle transformation non-parfumée dans une étape de calcul de points sigma 29, à partir de ces données et de matrices de bruit Qi et Ri, correspondant, respectivement, au bruit de mesure et au bruit inhérent au système, matrices fournies par un module d'estimation de bruit à partir de l'état mesuré yx et de la commande Ui au premier échantillonnage t=l.
Dans l'étape suivante, illustrée sur la figure 3B, des positions estimées X'i, X'i1, X'2, X'i, X'i des points sigma correspondant à un premier échantillonnage (t=l) sont prédites en appliquant, aux points sigma initiaux Xo, Xo1, Xo2, Xo, Xo, l'étape de prédiction 30 de l'algorithme de filtrage de Kalman avec la commande ui pour ce premier échantillonnage t=l. Dans l'étape suivante de correction 31, les points sigma Χ°, Xi1, Xi, X3, Xi effectifs, illustrés sur la figure 3C, sont calculés à partir de l'état mesuré yi pour cet échantillonnage (t=l), et une nouvelle estimation d'état latent Xi et une matrice de covariance correspondante Pi sont calculées à partir en comparant ces nouveaux points sigma Xi, Xi1, Xi, X3, X4 aux positions estimées X'i, X'i, X'3, X'if X'i, comme illustré sur la figure 3D. Cet algorithme est récursif, et chaque étape à partir de celle de calcul de points sigma 29 va être répétée à chaque nouvel échantillonnage. Pour cela, le module 35 d'estimation de bruit fournit des nouvelles matrices de bruit Qt, Rt pour chaque nouvel échantillonnage à partir des états mesurés yt et des commandes ut observés pendant la séquence temporelle. Alternativement, les matrices de bruit Qt, Rt pour l'ensemble de la séquence d'échantillonnage pourrait toutefois être préprogrammé.
Les différences entre les positions X't-i, X't.i, X't-2, X't-3, X't-4 prédites sur base des points sigma précédents Xt-i, Xt-i1, Xt-i, Xt-i, Xt-i et celles X°, Xt, X2, X3, X4 effectivement calculées sur base de chaque nouvel échantillonnage permettent d'obtenir de l'information, sur la fonction F d'évolution de l'état latent xt dans le temps. Cette information, contenue dans la matrice de covariance Pt, peut être utilisée, dans une boucle parallèle d'optimisation 32, avec une récurrence égale ou inférieure à celle de l'observateur d'états 27, pour ajuster les paramètres de la fonction de dynamique modèle F. Ainsi, dans le mode de réalisation illustré, une succession de matrices de covariance Pt sur une période de temps sont transmises à un intégrateur 33, dans lequel la somme de l'ensemble des éléments de chaque matrice Pt dans cette succession est intégrée pour obtenir une valeur C à minimiser par un algorithme d'adaptation paramétrique 34 , qui peut être, par exemple, l'algorithme du gradient, l'algorithme de Newton, ou les algorithmes de DFP ou BFGS (méthodes de Quasi-Newton).
Dans ce mode de réalisation, cet algorithme d'adaptation paramétrique 34 peut ajuster itérativement les valeurs du premier paramètre rp et du deuxième paramètre ov sur deux itérations successives ou, alternativement, de façon simultanée lors d'une itération de la boucle d'optimisation 32, de manière à minimiser la valeur C, permettant ainsi d'identifier des valeurs pour ces paramètres correspondant très étroitement au comportement réel de la turbine 11 et de la vanne 12 et de les utiliser dans l'étape de calcul de points sigma 29.
Parmi les différents types de matrices de covariance pouvant être produites pour chaque échantillonnage successif par le filtre de Kalman « sans odeur », on peut notamment utiliser pour cette boucle d'optimisation 32 la matrice de covariance de l'innovation réelle, c'est-àdire la matrice résultant, pour chaque moment t, de la matrice des différences entre l'état mesuré yt et une estimation d'état mesuré yt calculée à partir de l'estimation d'état latent xt, multipliée par la transposée de celle-ci. Ce critère traduit la variation de l'innovation au cours du temps. Une réduction de cette variation implique donc une plus grande fidélité de la fonction de dynamique modèle F à la réalité.
La valeur C à minimiser peut toutefois être calculée de manière différente, sans passer par des matrices de covariance Pt. Elle pourrait notamment être calculée à partir des différences entre une sélection de valeurs mesurées et les estimations correspondantes fournies par le filtre de Kalman « sans odeur », et en particulier à partir de sommes d'erreurs quadratiques. Ainsi, par exemple, C pourrait résulter de l'intégration d'une somme pondérée du carré de la différence entre valeur mesurée et estimée de la vitesse de rotation de l'arbre rotatif de la turbine 11 et le carré de la différence entre valeur mesurée et estimée de la pression en aval de la pompe 14. Une autre alternative envisageable basée sur la somme d'erreurs quadratiques serait le calcul de la valeur C en intégrant de la trace de la matrice résultant de la multiplication par sa matrice transposée de la matrice des différences entre l'état mesuré yt et une estimation d'état mesuré yt, suivant l'expression suivante :
C = Σΐ=ο tr((yk - yk)(yk - ÿk)T)
Finalement, encore une autre alternative envisageable basée sur la somme d'erreurs quadratiques ne prendrait plus seulement la trace de la matrice en considération, mais toutes les valeurs de la matrice, suivant l'expression :
t η n c = ” yk)(./k - ykY^j k=o i=iy=i où n est la dimension du vecteur d'état mesuré yt.
Par ailleurs, même si le filtre de Kalman « sans odeur » est spécifiquement utilisé dans ce mode de réalisation, d'autres algorithmes dérivés du filtre de Kalman et destinés à être appliqués à des systèmes non-linéaires, comme par exemple le filtre de Kalman étendu, peuvent alternativement être utilisés de manière analogue.
En outre, quoique le système physique non-linéaire 10 dans ce mode de réalisation soit spécifiquement un système relativement simple, comportant uniquement deux éléments à caractériser (la turbine 11 et la vanne 12), le procédé d'essai suivant l'invention est également applicable à des systèmes sensiblement plus complexes, comprenant plus de deux éléments à caractériser et donc plus de deux paramètres dont les valeurs doivent être ajustées. Ainsi, ce procédé pourrait même être appliqué à des essais de moteurs fusées complets, avec des pluralités de turbines et de vannes à caractériser.
Quoique la présente invention ait été décrite en se référant à un exemple de réalisation spécifique, il est évident que des différentes modifications et changements peuvent être effectués sur cet exemple sans sortir de la portée générale de l’invention telle que définie par les revendications. Par conséquent, la description et les dessins doivent être considérés dans un sens illustratif plutôt que restrictif.
Claims (15)
- REVENDICATIONS1. Procédé d'essai d'un système physique non-linéaire (10), comprenant les étapes suivantes :effectuer une succession d'échantillonnages de valeurs de commande et de valeurs de mesure dudit système physique non-linéaire ;soumettre lesdites valeurs de commande et de mesure à un filtre de Kalman discret non-linéaire utilisant un modèle numérique dudit système non-linéaire ; et mettre en oeuvre un algorithme d'adaptation paramétrique pour optimiser au moins un premier et un deuxième paramètre dudit modèle numérique sur la base de valeurs issues dudit filtre de Kalman discret nonlinéaire.
- 2. Procédé d'essai suivant la revendication 1, dans lequel ledit filtre de Kalman discret non-linéaire est un filtre de Kalman sans-odeur.
- 3. Procédé d'essai suivant l'une quelconque des revendications 1 ou 2, dans lequel l'algorithme d'adaptation paramétrique est un algorithme d'optimisation desdits aux moins deux paramètres du modèle numérique par minimisation d'une valeur calculée sur base d'au moins une matrice de covariance issue dudit filtre de Kalman discret non linéaire.
- 4. Procédé d'essai suivant la revendication 3, dans lequel ladite au moins une matrice de covariance comprend une matrice de covariance d'innovation réelle dudit filtre de Kalman discret non-linéaire.
- 5. Procédé d'essai suivant l'une quelconque des revendications 3 à 4, dans lequel ladite valeur à minimiser comprend la somme de tous les éléments d'au moins une matrice de covariance issue dudit filtre de Kalman non linéaire.
- 6. Procédé d'essai suivant la revendication 5, dans lequel ladite valeur à minimiser comprend la somme de tous les éléments d'une succession de matrices de covariances issue dudit filtre de Kalman non linéaire pendant une séquence d'échantillonnage.
- 7. Procédé d'essai suivant l'une quelconque des revendications précédentes, dans lequel ledit système physique non-linéaire (10) comprend une turbine (11) et une vanne (12) de contournement de ladite turbine.
- 8. Procédé d'essai suivant la revendication 7, dans lequel le premier paramètre du modèle numérique du système physique nonlinéaire (10) correspond à un rapport entre une performance théorique et une performance réelle de la turbine pendant l'essai, et le deuxième paramètre du modèle numérique du système physique non-linéaire (10) correspond à un décalage entre une consigne de position et une position réelle de la vanne pendant l'essai.
- 9. Procédé d'essai suivant la revendication 8, dans lequel lesdites performances théorique et réelle sont, respectivement, un couple théorique et un couple réel généré par la turbine pendant l'essai.
- 10. Procédé d'essai suivant l'une quelconque des revendications 7 à 9, dans lequel lesdites valeurs de commande comprennent une consigne d'ouverture de la vanne de contournement, une vitesse mesurée de rotation de la turbine, et/ou un couple mesuré de la turbine.
- 11. Procédé d'essai suivant la revendication 10, dans lequel ledit couple mesuré de la turbine est mesuré au travers de mesures de pression et débit d'une pompe (14) actionnée par la turbine (11).
- 12. Procédé d'essai suivant l'une quelconque des revendications 7 à 11, dans lequel lesdites valeurs de mesure comprennent des mesures de pression et de température en amont et en aval de la turbine (11) et de la vanne (12) de contournement.
- 13. Procédé d'essai suivant l'une quelconque des revendications précédentes, dans lequel au moins une partie desdits échantillonnages successifs est effectuée pendant une période de transition d'un premier régime de fonctionnement à un deuxième régime de fonctionnement dudit système physique non-linéaire différent du premier régime.
- 14. Programme d'ordinateur configuré pour la mise en œuvre d'un 5 procédé d'essai suivant l'une quelconque des revendications précédentes.
- 15. Support informatique lisible par un dispositif informatique et comprenant des instructions pour la mise en œuvre d'un procédé d'essai suivant l'une quelconque des revendications 1 à 13 par un dispositif10 informatique lisant ledit support.1/3
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR1662001A FR3059777B1 (fr) | 2016-12-06 | 2016-12-06 | Procede d'essai |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR1662001A FR3059777B1 (fr) | 2016-12-06 | 2016-12-06 | Procede d'essai |
FR1662001 | 2016-12-06 |
Publications (2)
Publication Number | Publication Date |
---|---|
FR3059777A1 true FR3059777A1 (fr) | 2018-06-08 |
FR3059777B1 FR3059777B1 (fr) | 2019-08-09 |
Family
ID=58707623
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR1662001A Active FR3059777B1 (fr) | 2016-12-06 | 2016-12-06 | Procede d'essai |
Country Status (1)
Country | Link |
---|---|
FR (1) | FR3059777B1 (fr) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040102890A1 (en) * | 2002-11-27 | 2004-05-27 | Brunell Brent Jerome | Methods and apparatus for model predictive control of aircraft gas turbine engines |
US20070073525A1 (en) * | 2005-09-27 | 2007-03-29 | General Electric Company | Method and system for gas turbine engine simulation using adaptive Kalman filter |
US20090198350A1 (en) * | 2008-01-31 | 2009-08-06 | Fisher-Rosemount Systems, Inc | Robust adaptive model predictive controller with tuning to compensate for model mismatch |
US20090240480A1 (en) * | 2008-03-19 | 2009-09-24 | Honeywell International Inc. | Target trajectory generator for predictive control of nonlinear systems using extended kalman filter |
EP2562612A1 (fr) * | 2011-08-24 | 2013-02-27 | General Electric Company | Procédés et systèmes de modélisation de turbine à gaz utilisant un filtrage adaptatif de Kalman |
WO2015012832A1 (fr) * | 2013-07-25 | 2015-01-29 | General Electric Company | Surveillance dynamique, diagnostic et commande de systèmes de tour de refroidissement |
US20150134270A1 (en) * | 2013-11-12 | 2015-05-14 | Aero Systems Engineering, Inc. | Systems and methods for improved accuracy |
-
2016
- 2016-12-06 FR FR1662001A patent/FR3059777B1/fr active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040102890A1 (en) * | 2002-11-27 | 2004-05-27 | Brunell Brent Jerome | Methods and apparatus for model predictive control of aircraft gas turbine engines |
US20070073525A1 (en) * | 2005-09-27 | 2007-03-29 | General Electric Company | Method and system for gas turbine engine simulation using adaptive Kalman filter |
US20090198350A1 (en) * | 2008-01-31 | 2009-08-06 | Fisher-Rosemount Systems, Inc | Robust adaptive model predictive controller with tuning to compensate for model mismatch |
US20090240480A1 (en) * | 2008-03-19 | 2009-09-24 | Honeywell International Inc. | Target trajectory generator for predictive control of nonlinear systems using extended kalman filter |
EP2562612A1 (fr) * | 2011-08-24 | 2013-02-27 | General Electric Company | Procédés et systèmes de modélisation de turbine à gaz utilisant un filtrage adaptatif de Kalman |
WO2015012832A1 (fr) * | 2013-07-25 | 2015-01-29 | General Electric Company | Surveillance dynamique, diagnostic et commande de systèmes de tour de refroidissement |
US20150134270A1 (en) * | 2013-11-12 | 2015-05-14 | Aero Systems Engineering, Inc. | Systems and methods for improved accuracy |
Also Published As
Publication number | Publication date |
---|---|
FR3059777B1 (fr) | 2019-08-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2744977C (fr) | Detection d'anomalie dans un moteur d'aeronef | |
EP3152528B1 (fr) | Procédé et système d'évaluation d'un débit d'un fluide | |
EP3941826A1 (fr) | Système et procédé de surveillance d'un moteur d'aéronef | |
FR3069668A1 (fr) | Procede et dispositif de recherche d'un defaut susceptible d'affecter un dispositif mecanique tournant de transmission de puissance | |
FR3035232A1 (fr) | Systeme de surveillance de l'etat de sante d'un moteur et procede de configuration associe | |
CA2918215C (fr) | Procede d'estimation sur une courbe d'un point pertinent pour la detection d'anomalie d'un moteur et systeme de traitement de donnees pour sa mise en oeuvre | |
FR3059777B1 (fr) | Procede d'essai | |
WO2021111093A1 (fr) | Procede de surveillance d'une turbomachine, dispositif, systeme, aeronef et produit programme d'ordinateur | |
WO2021089936A1 (fr) | Procédé, dispositif et programme d'ordinateur de surveillance d'une machine tournante d'un aéronef | |
EP3721300B1 (fr) | Procede d'optimisation des performances d'un asservissement d'un systeme mecatronique, dispositif adapte | |
FR3140181A1 (fr) | Procédé de commande d’un système mécatronique | |
EP3458827A1 (fr) | Estimation d'une grandeur physique en sortie d'un circuit hydraulique | |
FR3060157A1 (fr) | Procede de simulation du comportement d'une structure par calcul par elements finis | |
FR3043432A1 (fr) | Procede et dispositif d'ajustement de valeurs de consigne d'une pluralite de geometries variables determinees d'une turbomachine | |
FR3012882A1 (fr) | Procede d'essai technique | |
CA2745977A1 (fr) | Procede et systeme d'estimation d'une temperature de veine dans un turboreacteur | |
FR3012897A1 (fr) | Procede et dispositif de caracterisation d'un signal | |
WO2022096354A1 (fr) | Procédé de détermination de défaut d'un rendement d'un module d'un turbomoteur d'un aéronef | |
WO2021180726A1 (fr) | Procédé et système de contrôle non destructif d'une pièce aéronautique | |
WO2010004098A1 (fr) | Estimation de la réponse impulsionnelle d'un système à partir d'observations binaires | |
WO2023209111A1 (fr) | Procede et systeme de detection d'anomalie(s) dans une structure | |
EP4200580A1 (fr) | Procede de determination de panne pour un ensemble de capteurs | |
WO2020182808A1 (fr) | Correction de la mesure de pression d'un capteur de pression de cylindre | |
WO2023187287A1 (fr) | Procédé de surveillance de l'état de santé de turbomachine d'aéronef | |
EP3969967A1 (fr) | Recalage de modèle dans une turbomachine |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PLFP | Fee payment |
Year of fee payment: 2 |
|
PLSC | Publication of the preliminary search report |
Effective date: 20180608 |
|
PLFP | Fee payment |
Year of fee payment: 3 |
|
PLFP | Fee payment |
Year of fee payment: 4 |
|
PLFP | Fee payment |
Year of fee payment: 5 |
|
PLFP | Fee payment |
Year of fee payment: 6 |
|
PLFP | Fee payment |
Year of fee payment: 7 |
|
PLFP | Fee payment |
Year of fee payment: 8 |