FR 3030037 A1 20160617 FR 1462368 A 20141212 1 Arrière-plan de l'invention La présente invention concerne le domaine des essais techniques et plus particulièrement un procédé d'essai de machine. Typiquement, dans l'essai d'une machine, on commande le fonctionnement de la machine suivant une séquence prédéterminée de points de fonctionnement, de manière à permettre l'évaluation du fonctionnement de la machine dans chacun de ces points de fonctionnement. Chaque point de fonctionnement est défini par une durée et un ensemble de paramètres de fonctionnement de la machine. Ces essais peuvent être très coûteux, notamment quand ils concernent une machine de grande complexité technique, telle que par exemple un moteur-fusée à ergols liquides. Si l'essai doit être interrompu avant d'avoir atteint tous les points de fonctionnement que l'on souhaitait évaluer, soit par dysfonctionnement de la machine essayée, soit par épuisement de ressources telles que les ergols dans un essai de moteur-fusée, il faudra effectuer des nouveaux essais coûteux en temps et ressources, jusqu'à compléter l'ensemble des points de fonctionnement. Objet et résumé de l'invention La présente invention vise à remédier à ces inconvénients. L'invention vise à proposer un procédé d'essai d'une machine qui permette d'effectuer chaque essai suivant une séquence optimale de points de fonctionnement, de manière à maximiser le nombre et/ou l'importance des points de fonctionnement atteints à la fin de l'essai. Ce but est atteint grâce au fait qu'il comprend au moins les étapes suivantes : détermination d'une pluralité de points de fonctionnement de ladite machine, chaque point de fonctionnement étant défini par une durée minimale et une valeur spécifique d'au moins un paramètre de fonctionnement de la machine ; calcul d'un ensemble de distances entre paires de points de fonctionnement ; sélection d'une séquence optimale de points de fonctionnement en application d'un algorithme de résolution du problème du voyageur de 3030037 2 commerce audit ensemble de distances entre paires de points de fonctionnement ; et commande du fonctionnement de ladite machine suivant ladite séquence optimale de points de fonctionnement. 5 Grâce à ces dispositions, et notamment grâce à la sélection d'une séquence optimale en application de l'algorithme de résolution du problème du voyageur de commerce, il est possible de maximiser le nombre et/ou l'importance des points de fonctionnement atteints pendant l'essai. 10 La machine peut être une machine thermique, et notamment comprenant une ou plusieurs turbomachines. En particulier, la machine peut être un moteur-fusée à ergols liquides. Dans ce domaine, les essais sont particulièrement coûteux et soumis à d'importantes contraintes en temps et autres ressources, telles que la quantité d'ergols disponibles. En 15 outre, le risque de défaillance d'un tel moteur-fusée à ergols liquides pendant certains points de fonctionnement, entraînant ainsi une interruption prématurée de l'essai, peut être particulièrement élevé. En conséquence, les avantages de l'invention sont particulièrement marqués dans son application à une telle machine. Néanmoins, elle est aussi 20 applicable aux essais d'autres types de machines. Quand la machine à essayer est un moteur-fusée à ergols liquides, chaque point de fonctionnement peut être caractérisé par des valeurs spécifiques d'au moins une pression en chambre de combustion et un rapport entre débits d'ergols. Toutefois, d'autres paramètres de 25 fonctionnement peuvent aussi être utilisés alternativement ou en complément à ceux-là, comme par exemple des pressions des ergols en entrée de pompes d'alimentation. Pour calculer une distance d'un point de fonctionnement dans une paire de points de fonctionnement à l'autre point de fonctionnement dans 30 la paire, on peut notamment estimer un lapse de temps pour passer de l'un point de fonctionnement dans la paire à l'autre point de fonctionnement dans la paire. Ainsi, la séquence optimale sélectionnée suivant l'algorithme de résolution du problème du voyageur de commerce permettra de minimiser le temps nécessaire pour atteindre l'ensemble des 35 points de fonctionnement et/ou maximiser le nombre de points de fonctionnement atteints dans un lapse de temps limité. En outre, la 3030037 3 sélection de la séquence optimale pourra ainsi tenir compte de lapses de temps différents pour les passages entre points de fonctionnement d'une paire dans un sens et dans l'autre. Néanmoins, la distance entre points de fonctionnement peut aussi être basée sur d'autres critères : par exemple, 5 en complément ou alternativement à l'estimation du lapse de temps, on pourrait estimer la quantité d'ergol à consommer pour passer de l'un point de fonctionnement de la paire à l'autre point de fonctionnement de la paire. Dans certaines circonstances, on peut souhaiter ou être contraint à 10 faire succéder dans l'essai au moins deux points de fonctionnement suivant une sous-séquence prédéterminée. Afin de forcer cette succession dans la séquence optimale sélectionnée suivant l'algorithme de résolution du problème du voyageur de commerce, les distances d'un point de fonctionnement désigné précédent à tous les autres points de 15 fonctionnement sauf un point de fonctionnement désigné suivant peuvent être majorées de manière à forcer une succession directe du point de fonctionnement désigné précédent par le point de fonctionnement désigné suivant. En outre, certains points de fonctionnement peuvent être plus 20 prioritaires que d'autres. Pour prendre en compte cette priorité dans la sélection de la séquence optimale, au moins une distance d'un premier point de fonctionnement à un deuxième point de fonctionnement peut être minorée ou majorée suivant un coefficient de priorité associé audit deuxième point de fonctionnement. Ainsi, on peut minorer les distances 25 pour arriver aux points de fonctionnement les plus prioritaires, depuis tous les autres points de fonctionnement, facilitant ainsi que les points de fonctionnement les plus prioritaires soient arrangés en début de la séquence optimale. D'autre part, d'autres points de fonctionnement peuvent présenter 30 un risque particulièrement élevé d'interruption prématurée de l'essai. Pour prendre en compte ce risque dans la sélection de la séquence optimale, une distance d'un premier point de fonctionnement à un deuxième point de fonctionnement peut être majorée ou minorée suivant un coefficient de risque associé audit deuxième point de fonctionnement. Ainsi, on peut 35 majorer les distances pour arriver aux points de fonctionnement les plus risqués, depuis tous les autres points de fonctionnement, facilitant ainsi 3030037 4 que les points de fonctionnement les plus risqués soient arrangés en queue de la séquence optimale. Les algorithmes de résolution du problème du voyageur de commerce sont normalement adaptés pour donner une solution en forme 5 de circuit, avec le point final dans la séquence à proximité du point de départ. Toutefois, dans cette application, il peut ne pas être souhaité d'obtenir une séquence en forme de circuit. Pour éviter cela, il est possible de majorer toutes les distances d'un point de fonctionnement de départ à tous les autres points de fonctionnement. En rendant ainsi le point de 10 fonctionnement de départ identiquement distant, ou presque, de tous les autres points de fonctionnement, on minimise la tendance des algorithmes de résolution du voyageur de commerce à sélectionner une séquence optimale en forme de circuit, avec un point de fonctionnement final proche du point de fonctionnement de départ. 15 Dans certains cas, il peut être souhaitable de réarranger la séquence optimale en cours d'exécution d'essai. Ainsi, par exemple, après que tous les points de fonctionnement avec une priorité maximale aient été atteints, on peut neutraliser le coefficient de priorité de chacun des points de fonctionnement restants, de manière à ne plus prendre en compte la 20 priorité des points de fonctionnement restants dans une mise à jour subséquente de la séquence optimale. On peut aussi, après que tous les points de fonctionnement avec un risque minimum aient été atteints, neutraliser le coefficient de risque de chacun des points de fonctionnement restants, de manière à ne plus prendre en compte le 25 risque des points de fonctionnement restants dans une mise à jour subséquente de la séquence optimale. Ainsi, après qu'au moins un point de fonctionnement ait été atteint par la machine, on peut procéder à modifier des coefficients de majoration et/ou minoration desdites distances, calculer à nouveau ledit ensemble de distances avec les 30 coefficients modifiés, et mettre à jour ladite séquence optimale en appliquant ledit algorithme de résolution du problème du voyageur de commerce audit ensemble de distances après leur nouveau calcul. Le procédé d'essai suivant cette divulgation peut notamment être mis en oeuvre dans un banc d'essai à commande numérique. En 35 conséquence, cette invention concerne également un programme informatique comprenant un ensemble d'instructions destinées à mettre 3030037 5 en oeuvre le procédé d'essai susmentionné, ainsi qu'un dispositif informatique programmé pour mettre en oeuvre ce procédé d'essai, et un support informatique comprenant des instructions exécutables par un dispositif informatique pour mettre en oeuvre ce procédé d'essai. Dans le 5 présent contexte, on entend par « support informatique » tout dispositif de stockage de données pouvant être lu par un dispositif informatique. Un tel support informatique peut notamment être un dispositif de stockage magnétique de données, tel qu'un disque ou bande magnétique, un dispositif de stockage optique de données, tel qu'un disque optique, ou un 10 dispositif de stockage électronique de données, tel qu'une mémoire électronique volatile ou non-volatile. Brève description des dessins 15 L'invention sera bien comprise et ses avantages apparaîtront mieux, à la lecture de la description détaillée qui suit, de modes de réalisation représenté à titre d'exemples non limitatifs. La description se réfère aux dessins annexés sur lesquels : - la figure 1 est une vue schématique d'un moteur-fusée à ergols 20 liquides ; - la figure 2 est un schéma des étapes d'un procédé d'essai suivant un mode de réalisation ; - la figure 3 est un graphique représentant l'évolution de coefficients de risque pendant l'exécution de l'essai suivant un 25 mode de réalisation du procédé. Description détaillée de l'invention La figure 1 montre une vue schématique d'un moteur-fusée 1 à 30 ergols liquides. Ce moteur-fusée 1 comprend une chambre de combustion 2, prolongée par une tuyère convergente-divergente 3, ainsi qu'une première turbopompe 4, une deuxième turbopompe 5 et un générateur de gaz 6 destinés à assurer l'alimentation en ergols de la chambre de combustion 2 à partir de réservoirs correspondants (non illustrés). 35 Chacune des turbopompes comprend au moins une pompe et une turbine 3030037 6 reliées par un arbre rotatif commun, de manière à ce que la pompe puisse être actionnée par la turbine à travers l'arbre rotatif. En fonctionnement, un premier ergol liquide, qui peut par exemple être un ergol cryogénique comme l'hydrogène liquide, est pompé par la 5 première turbopompe 4 vers la chambre de combustion 2 et vers le générateur de gaz 6. Un deuxième ergol, qui peut aussi être un ergol cryogénique comme l'oxygène liquide, est pompé par la deuxième turbopompe 5 aussi vers la chambre de combustion 2 et vers le générateur de gaz 6. Dans la chambre de combustion 2, une réaction 10 chimique exothermique entre les deux ergols génère des gaz de combustion à haute température, accélérés ensuite à vitesse supersonique par leur détente dans la tuyère 3, de manière à générer une poussée par réaction dans le sens opposé. Comme illustré, les parois de la tuyère 3 peuvent être refroidies par les ergols avant leur injection dans la chambre 15 de combustion 2. D'autre part, la fraction des ergols alimentant le générateur de gaz 6 y réagit aussi de manière exothermique pour générer des gaz de combustion à haute température, dont la détente dans les turbines des deux turbopompes 4, 5 en assure l'actionnement. Des moteurs-fusées suivant des configurations différentes sont 20 également envisageables. Ainsi, par exemple, dans les moteurs-fusées dits à cycle expandeur, les turbopompes ne sont pas actionnées par des gaz de combustion provenant d'un générateur de gaz, mais par au moins un des ergols après avoir été chauffé par de la chaleur générée dans la chambre de combustion. La présente invention est également applicable 25 dans ces configurations alternatives, ainsi qu'aux essais d'autres types de machines, et en particulier de machines thermiques et de turbomachines. La figure 2 illustre schématiquement les étapes d'un procédé d'essai d'un tel moteur-fusée suivant un mode de réalisation. Dans une étape d'initialisation S101, on détermine l'ensemble de points de fonctionnement 30 du moteur-fusée qu'on souhaite tester. Chaque point de fonctionnement peut être caractérisé par une durée minimale et un ensemble de paramètres de fonctionnement, comme par exemple une pression dans la chambre de combustion 2 et un rapport entre un débit (massique ou volumique) du premier ergol en entrée de la première turbopompe 4 et un 35 débit (massique ou volumique) du deuxième ergol en entrée de la deuxième turbopompe 5. Alternativement ou en complément à ces 3030037 7 paramètres de fonctionnement, on peut utiliser d'autres paramètres de fonctionnement, et notamment ceux connus comme « conditions d'interface », c'est-à-dire la pression du premier ergol en entrée de la première turbopompe 4 et la pression du deuxième ergol en entrée de la 5 deuxième turbopompe. Dans l'étape suivante S102, on détermine un premier ensemble de distances entre points de fonctionnement. Dans ce premier ensemble, chaque distance c(A,B) d'un premier point de fonctionnement A à un deuxième point de fonctionnement B dans ce premier ensemble peut être 10 définie comme correspondant au temps nécessaire pour faire passer le moteur-fusée du point de fonctionnement A au point de fonctionnement B. Alternativement au temps, cette distance peut-être un volume et/ou masse d'ergol consommé pendant la transition, ou bien une fonction d'au moins un du temps et de la quantité d'ergol. Il faut noter que ces 15 distances ne sont pas nécessairement symétriques, c'est-à-dire que la distance c(A,B) peut être différente de la distance c(B,A), puisqu'il peut être plus ou moins difficile de passer du point de fonctionnement A au point de fonctionnement B que vice versa. Cet ensemble de distances peut être exprimé sous la forme d'une matrice carrée. 20 Ensuite, dans l'étape suivante S103, les distances du point de fonctionnement de départ à tous les autres points de fonctionnement sont majorées. Idéalement, ces distances pourraient être rendues infinies, mais à cause des ressources numériques limitées disponibles, on peut se contenter de majorer ces distances de plusieurs ordres de grandeur. 25 L'objet de cette majoration est d'éviter que, lors de l'application subséquente d'un algorithme de résolution du problème du voyageur de commerce, soit calculée une séquence de points de fonctionnement qui finisse nécessairement à proximité du point de fonctionnement de départ. En effet, les algorithmes de résolution du problème du voyageur de 30 commerce résultent normalement en une solution en forme de circuit, avec un point d'arrivée proche du point de départ. Une telle solution n'étant pas nécessairement adaptée au présent problème, on éloigne donc artificiellement le point de fonctionnement de départ de tous les autres points de fonctionnement afin de l'éviter. 35 La majoration des distances peut aussi servir à forcer la succession de certains points par d'autres dans la séquence résultante de l'application 3030037 8 subséquente de l'algorithme de résolution du problème du voyageur de commerce. Ainsi, dans l'étape suivante S104, quand on a désigné deux points de fonctionnement comme devant nécessairement se succéder dans l'essai, les distances du point de fonctionnement désigné précédent à 5 tous les autres points de fonctionnement sauf le point de fonctionnement désigné suivant sont majorées de manière à forcer une succession directe, dans la séquence qui sera calculée par l'algorithme de résolution du problème du voyageur de commerce, du point de fonctionnement désigné précédent par le point de fonctionnement désigné suivant. Ceci peut être 10 naturellement mis en oeuvre pour plus de deux points de fonctionnement. Idéalement, ces distances majorées pourraient aussi être rendues infinies, mais à cause des ressources numériques limitées disponibles, on peut également se contenter de majorer ces distances de plusieurs ordres de grandeur. 15 Dans l'étape suivante S105, un coefficient de priorité kp et un coefficient de priorité kr sont appliqués aux distances c(A,B) de ce premier ensemble de distances, pour obtenir un deuxième ensemble de distances c'(A,B) suivant la formule c'(A,B)=c(A,B)-kp-kr. Ce deuxième ensemble de distances peut également être exprimé en forme de matrice carrée, 20 comme le premier ensemble. Le coefficient de priorité kp appliqué à chaque distance c(A,B) d'un premier point de fonctionnement A à un deuxième point de fonctionnement B dépend d'une priorité préalablement accordée à ce deuxième point de fonctionnement B. Plus précisément, le coefficient de 25 priorité kp diminue avec une priorité croissante du deuxième point de fonctionnement B. Ainsi, toutes les distances vers les points de fonctionnement les plus prioritaires peuvent être minorées, tandis que toutes celles vers les points de fonctionnement les moins prioritaires peuvent être majorées, de manière à privilégier le passage par les points 30 de fonctionnement les plus prioritaires lors de la résolution du problème du voyageur de commerce. Par exemple, avec trois niveaux de priorité, le coefficient de priorité kp appliqué aux distances vers les points de fonctionnement les plus prioritaires peut recevoir une valeur de 0,8, et celui appliqué aux distances vers les points de fonctionnement les moins 35 prioritaires peut recevoir une valeur de 1,2, tandis que les distances vers les points de fonctionnement à priorité intermédiaire reçoivent une valeur 3030037 9 neutre de 1. Toutefois, d'autres valeurs, voire un nombre différent de niveaux de priorité peuvent être appliqués suivant les circonstances. Le coefficient de risque kr appliqué à chaque distance c(A,B) d'un premier point de fonctionnement A à un deuxième point de 5 fonctionnement B dépend d'un niveau de risque associé à ce deuxième point de fonctionnement B. Plus précisément, au moins initialement, le coefficient de risque augmente avec un niveau de risque croissant associé au deuxième point de fonctionnement. Ainsi, toutes les distances vers les points de fonctionnement les plus risqués peuvent être initialement 10 majorées, de manière à privilégier, ou même forcer le passage par les points de fonctionnement les moins risqués d'abord. Toutefois, comme il sera expliqué ci-dessous, les valeurs du coefficient de risque kr peuvent évoluer pendant l'exécution de l'essai. Dans l'étape suivante S106, un algorithme de résolution du problème 15 du voyageur de commerce, tel que par exemple l'algorithme heuristique de Lin-Kernighan (algorithme LKH), est appliqué au deuxième ensemble de distances c'(A,B), de manière à sélectionner une séquence optimale de points de fonctionnement, qui commence à être appliquée dans l'étape S107, dans laquelle on exécute l'essai en commandant le fonctionnement 20 du moteur-fusée 1 suivant d'abord au moins les points de fonctionnement initiaux suivant la séquence optimale. Pendant l'exécution de l'essai, on peut néanmoins procéder à une ou plusieurs mises à jour de cette séquence optimale, en actualisant notamment les données concernant les points de fonctionnement déjà 25 couverts, et les ressources consommées pendant ces points de fonctionnement, en particulier le temps et/ou la quantité d'ergols. En outre, comme brièvement mentionné ci-dessus, on peut modifier dans l'étape S108 les valeurs du coefficient de risque kr pendant l'essai. Un exemple d'une telle évolution du coefficient de risque kr en fonction du 30 temps t est illustré sur la Figure 3. Suivant cet exemple, on attribue trois niveaux de risque R1, R2 et R3 aux différents points de fonctionnement. Initialement, le coefficient de risque kr appliqué aux distances vers les points de fonctionnement avec le niveau de risque moins élevé, R1, reçoit une valeur neutre de 1, tandis que ceux appliqués aux distances vers les 35 points de fonctionnement avec les niveaux de risque plus élevés R2 et R3 reçoivent des valeurs très élevées de manière à forcer l'algorithme de 3030037 10 résolution du problème du voyageur de commerce à commencer la séquence optimale par les points de fonctionnement moins risqués. Quand on approche la fin de l'essai, la valeur du coefficient de risque kr appliqué aux distances vers les points de fonctionnement avec le niveau de risque 5 le moins élevé, R1, peut être augmenté sensiblement, de manière à bloquer les points de fonctionnement à bas risque qui n'aient pas encore été atteints, tandis que les valeurs de ceux appliqués aux distances vers les points de fonctionnement avec les niveaux de risque plus élevés R2 et R3 sont diminuées, de manière à permettre leur sélection. Dans un 10 premier temps, la valeur du coefficient de risque kr appliqué aux distances vers les points de fonctionnement à risque intermédiaire, niveau de risque R2, peuvent diminuer plus que les distances vers les points de fonctionnement au niveau de risque le plus élevé, R3, de manière à favoriser le passage par les premiers plutôt que les seconds. Vers la fin de 15 l'essai, cette relation peut encore être inversé, de manière à favoriser plutôt le passage par les points de fonctionnement au risque le plus élevé. Après l'actualisation de ces données, on répète les étapes S105 et S106 pour recalculer le deuxième ensemble de distances c'(A,B) avec les données actualisées et la séquence optimale avec ce deuxième ensemble 20 de distances, et on poursuit l'essai avec l'étape S107. Ce procédé d'essai peut être mis en oeuvre par le moyen d'un banc d'essai à commande numérique. Un tel banc d'essai peut comprendre une unité de commande apte à être relié à la machine à tester afin de commander son fonctionnement. Cette unité de commande peut aussi 25 être programmée pour effectuer la sélection de la séquence optimale avant l'essai et/ou pour sa mise à jour en cours d'essai, en fonction des données actualisées grâce à des mesures obtenues par des capteurs intégrés ou connectés audit banc d'essai. Quoique la présente invention ait été décrite en se référant à des 30 exemples de réalisation spécifiques, il est évident que des différentes modifications et changements peuvent être effectués sur ces exemples sans sortir de la portée générale de l'invention telle que définie par les revendications. En outre, des caractéristiques individuelles des différents modes de réalisation évoqués peuvent être combinées dans des modes de 35 réalisation additionnels. Par conséquent, la description et les dessins doivent être considérés dans un sens illustratif plutôt que restrictif. REVENDICATIONS1. Procédé d'essai d'une machine, comprenant les étapes suivantes : détermination (S101) d'une pluralité de points de fonctionnement de ladite machine, chaque point de fonctionnement étant défini par une durée minimale et une valeur spécifique d'au moins un paramètre de fonctionnement de la machine ; calcul d'un ensemble de distances entre paires de points de 10 fonctionnement ; sélection (S106) d'une séquence optimale de points de fonctionnement en application d'un algorithme de résolution du problème du voyageur de commerce audit ensemble de distances entre paires de points de fonctionnement ; et 15 commande (S107) du fonctionnement de ladite machine suivant ladite séquence optimale de points de fonctionnement. 2. Procédé d'essai suivant la revendication 1, dans lequel ladite machine est un moteur-fusée (1) à ergols liquides. 3. Procédé d'essai suivant la revendication 2, dans lequel chaque 20 point de fonctionnement est caractérisé par des valeurs spécifiques d'au moins une pression en chambre de combustion (2) et un rapport entre débits d'ergols. 4. Procédé d'essai suivant l'une quelconque des revendications précédentes, dans lequel, pour calculer une distance d'un point de 25 fonctionnement dans une paire de points de fonctionnement à l'autre point de fonctionnement dans la paire, on estime un lapse de temps pour passer de l'un point de fonctionnement dans la paire à l'autre point de fonctionnement dans la paire. 5. Procédé d'essai suivant l'une quelconque des revendications 30 précédentes, dans lequel les distances d'un point de fonctionnement désigné précédent à tous les autres points de fonctionnement sauf un point de fonctionnement désigné suivant sont majorées de manière à forcer une succession directe, dans la séquence optimale, du point de fonctionnement désigné précédent par le point de fonctionnement désigné 35 suivant. 3030037 12 6. Procédé d'essai suivant l'une quelconque des revendications précédentes, dans lequel une distance d'un premier point de fonctionnement à un deuxième point de fonctionnement est minorée ou majorée suivant un coefficient de priorité associé audit deuxième point de 5 fonctionnement. 7. Procédé d'essai suivant l'une quelconque des revendications précédentes, dans lequel une distance d'un premier point de fonctionnement à un deuxième point de fonctionnement est majorée ou minorée suivant un coefficient de risque associé audit deuxième point de fonctionnement. 8. Procédé d'essai suivant l'une quelconque des revendications précédentes, dans lequel sont majorées toutes les distances d'un point de fonctionnement de départ à tous les autres points de fonctionnement. 9. Procédé d'essai suivant l'une quelconque des revendications précédentes, comprenant en outre au moins les étapes supplémentaires suivantes, après qu'au moins un point de fonctionnement ait été atteint par la machine : modification de coefficients de majoration et/ou minoration desdites distances ; nouveau calcul dudit ensemble de distances avec les coefficients modifiés ; mise à jour de ladite séquence optimale en appliquant ledit algorithme de résolution du problème du voyageur de commerce audit ensemble de distances après leur nouveau calcul. 10. Programme informatique comprenant un ensemble d'instructions destinées à mettre en oeuvre un procédé d'essai suivant l'une quelconque des revendications précédentes.