FR 3059799 A1 20180608 FR 1662048 A 20161207 1 DOMAINE GENERAL L'invention concerne un procédé de simulation de répartition d'aubes autour d'un disque de façon à former une roue aubagée pour turbomachine. ETAT DE L'ART Des aubes destinées à être montées autour d'un disque de turbomachine sont sujettes à de faibles dispersions au cours de leur fabrication. Ces dispersions peuvent être mesurables à froid sur chaque aube fabriquée, c'est-à-dire avant leur montage sur le disque et la mise en rotation de la roue aubagée obtenue à l'issue de ce montage. Ces dispersions peuvent également être mesurées à chaud sur les différentes aubes, c'est-à-dire au cours de la rotation de la roue aubagée. De telles dispersions sont susceptibles de rendre la roue aubagée déséquilibrée. Les conséquences d'un tel déséquilibre sont multiples: vibrations de l'ensemble du moteur, et/ou génération d'un bruit acoustique important pendant la mise en rotation de la roue aubagée. Pour répartir des aubes de façon équilibrée autour d'un disque de 20 turbomachine, il a été proposé de mettre en oeuvre un algorithme configuré pour rechercher, dans un ensemble de répartitions d'aubes possibles, une répartition optimisant une fonction de coût prédéterminée. Pour n aubes à monter autour d'un disque, une telle recherche est effectuée dans un ensemble de (n-1)! /2 répartitions d'aubes, en 25 considérant la roue identique par rotation et symétrie. Une telle recherche requiert donc un grand nombre de calculs si le nombre n est élevé. A titre d'exemple, pour n=18, qui est un nombre classique d'aubes dans une roue aubagée pour turbomachine, l'ensemble des répartitions sélectionnables comprend plus de 2.1014 répartitions. 30 Afin de limiter une telle charge de calcul, le document « POPMUSIC - Partial Optimization Metaheuristic under Special Intensification Conditions 3059799 2 », par Eric Taillard et publié le ler janvier 2002, propose de rechercher une répartition d'aubes optimisant une fonction de coût prédéterminée au moyen d'un algorithme de type Tabou. L'utilisation d'un tel algorithme de type Tabou permet d'explorer 5 l'espace défini par l'ensemble des répartitions possibles de façon relativement rapide. L'algorithme proposé dans ce document passe d'une répartition à l'autre par une permutation de deux aubes. Toutefois, cette manière d'explorer l'ensemble des répartitions possibles 10 peut requérir l'examen d'un nombre substantiel de répartitions lorsque ta fonction de coût utilisée dépend de plusieurs paramètres d'équilibrage d'aube. En d'autres termes, cet algorithme converge vers une solution optimale vérifiant la fonction coût de façon relativement lente. 15 PRESENTATION DE L'INVENTION Un but de l'invention est de répartir des aubes autour d'un disque de turbomachine de façon équilibrée, moyennant une charge de calcul limitée. Afin d'atteindre ce but l'invention propose un procédé de simulation 20 de répartition d'aubes sur un disque de turbomachine de sorte à former une roue aubagée, le procédé comprenant les étapes de : - fourniture d'une pluralité de configurations d'aube, chaque configuration d'aube étant associée à une aube et comprenant au moins une mesure de paramètre d'équilibrage mesurée sur l'aube 25 associée, - recherche et sélection d'une répartition de roue aubagée associant les configurations d'aube fournies avec des positions sur le disque , la répartition de roue aubagée favorisant l'atteinte d'au moins un critère défini d'après une fonction de coût prédéterminée dépendant 30 des mesures de paramètres d'équilibrage, la recherche et la sélection étant conduites par itérations successives, te procédé étant 3059799 3 caractérisé par le fait qu'une itération courante de la recherche et de la sélection comprend les étapes de: o génération d'une pluralité de répartitions voisines à partir d'une répartition de référence, ladite pluralité comprenant au 5 moins une première répartition voisine générée selon une stratégie favorisant l'atteinte d'un premier ensemble de critères de ressemblance prédéterminé avec la répartition de référence et au moins une deuxième répartition voisine générée selon une autre stratégie favorisant l'atteinte d'un 10 deuxième ensemble de critères de ressemblance prédéterminé avec la répartition de référence, o sélection, parmi les répartitions voisines, d'une répartition voisine favorisant une optimisation de la fonction de coût comme nouvelle répartition de référence. 15 Au cours de l'étape de génération de répartitions voisines, des critères de ressemblance différents sont utilisés pour générer ces voisins. Ces critères de ressemblance différents définissent des voisinages différents de la répartition de référence. Le fait d'utiliser différents critères de ressemblance (et donc de 20 définir différents voisinages pour une répartition de référence) permet de définir plusieurs stratégies d'exploration de l'espace des répartitions d'aube. Une répartition voisine peut ainsi être relativement équilibrée du point de vue d'un premier type de paramètre d'équilibrage, mais 25 relativement déséquilibrée d'un point de vue d'un autre type de paramètre d'équilibrage. Ainsi, les voisins définis par les différents ensembles de critères de ressemblance ne permettent pas de se rapprocher dans les mêmes proportions d'une répartition solution hypothétique vérifiant la fonction de coût. L'étape de sélection permet de réaliser un arbitrage entre 30 les différentes stratégies offertes par les différents ensembles de critères 3059799 4 de ressemblance utilisés, et donc de converger plus rapidement vers une solution optimale. L'invention peut également être complétée par les caractéristiques suivantes, prises seules ou en une quelconque de leurs combinaisons 5 techniquement possibles. La recherche et la sélection peut être mise en oeuvre au moyen d'un algorithme Tabou, dans lequel au moins une répartition de roue aubagée ayant participé à au moins une sélection au cours d'une itération précédente n'est pas traitée comme une répartition voisine au cours de 10 l'itération courante. L'algorithme de type Tabou permet d'explorer l'espace défini par l'ensemble des répartitions possibles de façon rapide. En effet, au cours de chaque itération de la recherche et de la sélection, certaines répartitions ayant déjà été testées précédemment ne sont pas traitées à nouveau. 15 Comme on le verra dans la suite, un algorithme de type Tabou utilisé au cours de la recherche fait en sorte de diminuer le nombre de répartitions d'aube testées au cours de la recherche. Au moins une répartition voisine peut résulter d'un ou deux échanges de positions d'aube entre configurations d'aubes de la répartition de 20 référence. Une des positions échangées dans le cadre d'un des deux échanges et une des positions échangées dans le cadre de l'autre des deux échanges peuvent être diamétralement opposées. L'au moins une répartition voisine peut résulter d'une ou de deux 25 permutations circulaires de configurations d'aubes appliquées à la répartition de référence, chaque permutation circulaire étant opérée sur un sous-ensemble de configurations d'aube associées à des positions consécutives autour du disque dans la répartition de référence. Chacune des deux permutations circulaires appliquées peut être 30 représentative d'un déplacement de configuration d'aube depuis une position initiale vers une position de destination, les deux positions initiales 3059799 5 étant diamétralement opposées dans la répartition de référence, et les deux positions finales étant diamétralement opposées dans la répartition voisine obtenue. 5 Il peut être sélectivement généré au cours de l'itération courante un nombre de répartitions voisines inférieur au nombre total de premières répartitions voisines de la répartition de référence. Le procédé peut comprendre, pour la répartition de référence et pour chaque répartition voisine, le calcul d'au moins un balourd à partir de 10 mesures de paramètre d'équilibrage d'un même type contenues dans les configurations d'aube de la répartition, et à partir des positions associées aux configurations d'aube dans la répartition, la fonction de coût dépendant du balourd. La fonction de coût peut dépendre de plusieurs balourds, chaque 15 balourd se rapportant à un type de mesure de paramètre d'équilibrage respectif, et dans lequel au moins un des balourds fait l'objet d'un seuillage par application de la fonction de coût. Les mesures de paramètre d'équilibrage sont d'au moins un des types suivants: angle de calage d'aube, moment axial d'aube, moment radial 20 d'aube, et moment tangentiel d'aube. Il est également proposé un procédé de fabrication d'une roue aubagée pour turbomachine à partir d'aubes et d'un disque, le procédé comprenant les étapes de mise en oeuvre du procédé de simulation qui précède, et d'assemblage des aubes sur le disque d'après la répartition de 25 roue aubagée sélectionnée par le procédé de simulation. Il est en outre proposé un produit programme d'ordinateur comprenant des instructions de code de programme pour exécuter les étapes du procédé de simulation qui précède, lorsque ce programme est exécuté par une unité de calcul. 30 3059799 6 DESCRIPTION DES FIGURES D'autres caractéristiques, buts et avantages de l'invention ressortiront de la description qui suit, qui est purement illustrative et non limitative, et qui doit être lue en regard des dessins annexés sur lesquels : 5 La figure 1 représente une roue aubagée pour turbomachine. Les figures 2a, 2b et 2c détaillent certains paramètres mesurables sur une aube de turbomachine. La figure 3 représente schématiquement un dispositif de répartition d'aubes selon un mode de réalisation de l'invention. 10 La figure 4 est un logigramme d'étapes d'un procédé de répartition d'aubes selon un mode de réalisation de l'invention. La figure 5 est un logigramme détaillant une étape représentée en figure 4. La figure 6 représente de façon schématique une répartition de roue 15 aubagée de référence et des répartitions voisines de cette répartition de références, générées au cours d'un procédé de répartition d'aube selon un mode de réalisation de l'invention. La figure 7 représente une fonction de coût selon un mode de réalisation de l'invention. Sur l'ensemble des figures, les éléments similaires portent des références identiques. DESCRIPTION DETAILLEE DE L'INVENTION En référence à la figure 1, une roue aubagée pour turbomachine 25 comprend un disque D et une pluralité d'aubes P réparties sur le disque autour de son axe. On a illustré sur les figures 2a, 2b, 2c une aube P destinée à être assemblée sur le disque D. L'aube P présente notamment les caractéristiques intrinsèques 30 suivantes : un moment axial A, un moment tangentiel T, un moment radial R. L'aube P présente également un angle de calage. 3059799 7 En référence à la figure 2a, l'angle de calage a peut être décrit comme l'angle entre la corde de l'aube P et un axe de rotation Ux du disque, également appelé « axe moteur » dans ce qui suit. En référence aux figures 2b et 2c, l'aube P présente un centre de 5 gravité G. Un moment statique est le produit de la masse de l'aube P et d'une des coordonnées (Dxcdg, DYcdg, Dzcdg) du centre de gravité G de l'aube P pris dans un repère direct orthonormé (0, Ux, Uy, Uz) avec Ux étant l'axe moteur, et O étant le centre du repère sur l'axe moteur. Pour le moment radial R, la coordonnée multipliée à la masse de 10 l'aube est celle du centre de gravité de l'aube prise sur l'axe U. Pour le moment axial A, la coordonnée multipliée à la masse de l'aube est celle du centre de gravité de l'aube prise sur l'axe U. Pour le moment tangentiel T, la coordonnée multipliée à la masse de l'aube est celle du centre de gravité de l'aube prise sur l'axe Uy. 15 En référence ta figure 3, un dispositif 1 pour répartir une pluralité d'aube P sur le disque D pour former une roue aubagée de turbomachine, comprend une mémoire de stockage 2, une unité de calcul 4, une mémoire tampon 6 et une interface homme-machine 8. La mémoire de stockage 2 est configurée pour stocker une base de 20 données contenant des données représentatives d'aubes (par exemple les moments et calage à froid, la masse de l'aube, l'angle de calage de l'aube, la hauteur de l'aube, de manière non limitative) et de positions de ces aubes sur un disque D. L'unité de calcul 4 est configurée pour mettre en oeuvre un 25 algorithme de recherche de type Tabou, qui sera détaillé ci-après, à partir des données stockées dans la mémoire 2. L'unité de calcul 4 a accès en lecture et en écriture à la mémoire de stockage 2. La mémoire tampon 6 est également accessible en lecture et en écriture par l'unité de calcul 4. Cette mémoire tampon 6 est configurée 30 pour stocker des données temporaires calculées par l'unité de calcul 4. En variante, les mémoires 2 et 6 forment une seule et même mémoire. 3059799 8 L'interface homme-machine 8 comprend par exemple un écran d'affichage et des moyens de saisie telle qu'un clavier manipulable par un opérateur. L'interface homme-machine 8 offre la possibilité à cet opérateur de modifier des paramètres d'exécution de l'algorithme mis en oeuvre par 5 l'unité de calcul 4. On va maintenant décrire en relation avec la figure 4 les étapes d'un procédé de simulation de répartition sollicitant le dispositif 1. On considère que n aubes P individuelles ont précédemment été fabriquées, ainsi que le disque D, comprenant n emplacements libres pour 10 aube, répartis autour de son axe. Dans une étape préliminaire 100, une ou plusieurs mesures de paramètres d'équilibrage sont effectuées sur chacune des n aubes. Dans la suite, on considère qu'un paramètre d'équilibrage est une grandeur physique caractérisant une aube individuelle. 15 Il peut être envisagé comme mesure de paramètres d'équilibrage au moins l'un des types de mesures suivants: mesure d'angle de calage, mesure de moment axial, mesure de moment tangentiel, mesure de moment radial, au moyen de méthodes de mesures connues de l'homme du métier (par exemple à l'aide de balances). 20 Les valeurs de mesure sont mémorisées par la mémoire 2 dans une base de données, au cours d'une étape de mémorisation 102. La base de données est par exemple organisée sous la forme d'une table à plusieurs entrées, chaque entrée se rapportant à une des n aubes considérées. Chaque entrée comprend une ou plusieurs mesures réalisées 25 sur une même aube, l'ensemble formé par ces mesures étant dans la suite appelé « configuration d'aube ». En complément de ces mesures, une entrée peut comprendre un identifiant d'aubes unique et/ou un numéro de série. SN Angle (°) Radial (inch Axial (inch Tangentiel (inch g) Posi_Instrum Fixe ? g) g) KLJLKAER -0,22 12345 - 507 1 Oui 2000 3059799 9 FHFGHFGHDFH -0,31 12500 - 420 5 Non 2500 FHGFH 0,250 13000 - 560 12 Non 1500 HGHGGH -0,02 12000 - 670 6 Non 1800 GHFG 0,101 12700 - 450 2 Oui 1702 FHGF 0,200 11800 - 590 10 Non 1504 FHGF 0,500 13500 - 360 17 Non 1203 ... ... ... Dans la table ci-dessus, les lignes représentent des entrées respectives de la base de données. Les colonnes du tableau contiennent, de gauche à droite : un numéro de série SN unique permettant d'identifier une 5 aube, une valeur d'angle de calage (par rapport à un angle de référence) mesuré sur l'aube, une valeur de moment radial mesurée sur l'aube, une valeur de moment axial mesurée sur l'aube, une valeur de moment tangentiel mesurée sur l'aube, et une information indiquant si l'aube doit être positionnée à une position imposée ou fixe, ou que cette aube peut au 10 contraire être librement placée à n'importe laquelle des n positions prédéterminées. Sont également initiées au cours de l'étape 100 dans la base de données, n informations de position, chaque information définissant la position autour de l'axe du disque D d'un emplacement pour aube respectif. 15 Dans la suite, ces informations de positions seront simplement dénommées « positions d'aube », ou plus simplement positions. Une position d'aube peut être exprimée comme un angle, ou bien comme un indice i variant de 1 à n dans un repère lié à la roue aubagée. Dans la suite, on considèrera que les positions sont des indices à partir 20 desquels des angles correspondants dans le repère de la roue aubagée, peuvent être déterminés. A ce stade du procédé, les n positions d'aube mémorisées et les n configurations d'aube ne sont pas encore associées. 3059799 10 Dans une étape 106 de recherche et de sélection, l'unité de calcul 4 produit une pluralité de répartitions des aubes sur la roue aubagée, et recherche et sélectionne une de ces répartitions de roue aubagée qui optimise une fonction de coût prédéterminée. Cette fonction de coût est 5 représentative d'un critère d'équilibrage et dépend des paramètres d'équilibrage considérés. On définit une répartition de roue aubagée comme un n-uplet d'associations, chaque association associant une des n configurations d'aube et une des n position d'aube. (n-1)!/2 répartitions d'aubes 10 différentes peuvent donc être théoriquement générées au cours de l'étape 106. Comme on le verra dans la suite, un algorithme de type Tabou utilisé au cours de la recherche fait en sorte de diminuer le nombre de répartitions d'aube testées au cours de la recherche. Dans une étape d'assemblage 108, les n aubes individuelles sont 15 assemblées sur le disque D, en ses n emplacements, d'après la répartition sélectionnée par le dispositif 1 au cours de l'étape précédente 106 de recherche et de sélection. Recherche et sélection au moyen d'un algorithme de type Tabou 20 On va maintenant détailler en relation avec la figure 5 l'étape de recherche et de sélection 106 de la répartition d'aube à utiliser pour l'assemblage. A titre préliminaire, deux listes sont allouées dans la mémoire tampon 6 : une liste de candidats, et une liste dite « tabou ». La recherche est mise en oeuvre au moyen d'un algorithme de type Tabou, comprenant des itérations successives. Une itération courante comprend les sous-étapes suivantes. L'unité de calcul détermine une répartition de roue aubagée de référence. Dans le cas particulier de la première itération, la répartition de référence peut par exemple être générée aléatoirement et mémorisée dans 3059799 11 la mémoire tampon 6 et/ou dans la base de données de la mémoire 2. En variante, la première répartition de référence est définie par l'opérateur au cours d'une étape de paramétrage 104, via l'interface homme/machine 8. 5 Dans toute itération excepté la première, la répartition de référence a été générée au cours d'une itération précédente. L'unité de calcul 4 génère à partir de la répartition de référence plusieurs répartitions voisines de la répartition de référence, d'après plusieurs critères de ressemblance prédéterminés, chaque répartition 10 voisine se rapportant à l'un des critères de ressemblance prédéterminés. En d'autres termes, si q critères de ressemblance différents sont utilisés au cours de cette étape de génération, au moins q voisins définis selon des critères de ressemblance différents sont générés. Chaque critère de ressemblance prédéterminé correspond à une 15 fonction mathématique bijective de permutation (composition de transpositions) entre deux répartitions. En d'autres termes, une répartition qui respecte un critère de ressemblance prédéterminé est une répartition voisine au sens de ce critère. Pour chacun des critères de ressemblance utilisés, il existe un 20 nombre total de répartitions voisines satisfaisant ce critère de ressemblance susceptibles d'être générées. Toutefois, pour limiter la charge de calcul du procédé, seule une partie des répartitions voisines possibles au sens d'un critère de ressemblance donné est générée, leur sélection étant aléatoire. 25 Les répartitions voisines sont mémorisées dans la liste de candidats, par la mémoire tampon 6, après vidage de cette liste de candidats. L'unité de calcul 4 calcule le coût de chaque configuration voisine, par application de la fonction de coût aux mesures de paramètres d'équilibrages contenus dans la configuration voisine. 30 Plus précisément, l'unité de calcul 4 calcule, pour chaque répartition voisine, au moins un balourd se rapportant à un ou plusieurs types de mesure 3059799 12 de paramètre d'équilibrage ; la fonction de coût est appliquée au balourd de chaque répartition voisine. Un balourd correspond à l'excentrement du centre de gravité de l'ensemble de la roue. Ce centre de gravité est calculé par sommation 5 vectoriel de tous les moments de chacune des aubes P. L'unité de calcul 4 procède ensuite à une sélection, dans la liste de candidats, d'une répartition qui optimise le plus la fonction de coût (c'est-à-dire la répartition qui offre le meilleur équilibrage parmi les répartitions voisines). En fonction du type de fonction considéré, la répartition 10 sélectionnée sera la répartition présentant un coût minimal ou maximal. La répartition voisine sélectionnée est alors considérée comme nouvelle répartition de référence. Les répartitions voisines générées sont en outre mémorisées dans la liste « tabou », par la mémoire tampon 6. 15 Ces étapes sont ensuite répétées dans une itération suivante, sur la base de cette nouvelle répartition de référence, et ainsi de suite dans des itérations ultérieures, chaque nouvelle itération considérant une nouvelle répartition de référence. L'algorithme se termine lorsqu'un critère de sortie prédéterminé est 20 vérifié, par exemple, une valeur de la fonction de coût est atteinte par une répartition de référence et/ou un nombre d'itérations est atteint. Dans les deux cas, la répartition optimale choisie est la dernière répartition de référence considérée par l'algorithme de type Tabou mis en oeuvre. 25 Une particularité de l'algorithme de type Tabou utilisé est que, au cours de toute itération exceptée la première, il est fait en sorte que tes répartitions de roue aubagée ayant participées à au moins une élection au cours d'une itération précédente ne sont pas traitées comme des répartitions voisines au cours d'une itération courante. 30 Ce principe d'exclusion peut être mis en oeuvre dans une itération courante par comparaison, après l'étape de génération des répartitions 3059799 13 voisines, entre le contenu de la liste « tabou » (mise à jour au cours de l'itération immédiatement précédente) et le contenu de la liste de candidats. Toute répartition trouvée dans les deux listes est alors effacée de la liste de candidats. 5 L'exclusion peut être limitée aux seuls voisinages générés au cours de l'itération immédiatement précédente, afin de limiter la mémoire consommée par la liste tabou. En variante, l'exclusion peut porter sur les répartitions voisines traitées au cours des K dernières itérations précédente l'itération courante. 10 Cette variante permet d'explorer plus rapidement l'espace des répartitions de roue aubagée possibles. Cette opération d'exclusion peut en outre conduire à augmenter la valeur de la fonction (dans un problème de minimisation) : c'est le cas lorsque toutes les répartitions voisines ont un coût plus élevé que la 15 répartition de référence, qui forme alors un minimum local. Le risque cependant est qu'à l'étape suivante, on retombe dans le minimum local auquel on vient d'échapper. C'est pourquoi il faut que l'heuristique ait de la mémoire : le mécanisme consiste à interdire (d'où le nom de tabou) de revenir sur les dernières positions explorées. 20 La taille de la liste Tabou peut être un paramètre d'exécution de l'algorithme, susceptible d'être changé par un opérateur via l'interface homme/machine 8. Génération de répartitions voisines 25 Dans un mode de réalisation, quatre critères de ressemblances sont utilisés pour déterminer/générer les répartitions voisines. Une répartition voisine définie d'après un premier critère de ressemblance est une répartition de roue aubagée obtenue après échange des positions associées à deux configurations d'aube dans la répartition de 30 référence. 3059799 14 Est représenté en figure 6 un exemple de répartition de référence REF et une répartition voisine V1 générée selon ce premier critère de ressemblance. Les positions 2 et 3 de la répartition REF sont échangées pour obtenir la répartition V1. 5 Une répartition voisine définie d'après un deuxième critère de ressemblance est une répartition de roue aubagée obtenue après deux échanges des positions associées à deux configurations d'aube dans la répartition de référence. Deux paires de positions sont échangées, par exemple distinctes. 10 Par ailleurs, une des positions échangées dans le cadre d'un des deux échanges, et une des positions échangées dans le cadre de l'autre des deux échanges peuvent être diamétralement opposées. A titre d'exemple, les positions 2 et 3 de la répartition REF sont échangées, de même tes positions 12 et 11, pour obtenir la répartition V2, pour un nombre d'aubes égal à 18. 15 Le nombre d'échanges maximal autorisé pour respecter le deuxième critère peut être augmenté à une valeur supérieure à 2. Les positions définissant un échange peuvent être adjacentes (dans l'application du premier critère ou du deuxième critère). Une répartition voisine définie d'après un troisième critère de 20 ressemblance est une répartition de roue aubagée résultant d'une permutation circulaire de configurations d'aube appliquée à un sous-ensemble de configurations d'aubes associées à des positions consécutives autour du disque dans la répartition de référence. Une telle permutation représente virtuellement le déplacement 25 d'une aube d'une roue aubagée définie d'après la répartition de référence depuis une première position vers une deuxième position, et un décalage des aubes intercalées entre la première position et la deuxième position. A titre d'exemple, la répartition V3 illustrée en figure 6 est obtenue par une permutation circulaire appliquée à l'ensemble [2, ... 13] de la 30 configuration de référence REF. Cette permutation circulaire représente le 3059799 15 déplacement d'une aube virtuellement positionnée à la position 13, à la position 2. Une répartition voisine définie d'après un quatrième critère de ressemblance est une répartition de roue aubagée résultant de deux 5 permutations circulaires conformes au troisième critère. Chacune des deux permutations circulaires appliquées étant représentative d'un déplacement de configuration d'aube depuis une position initiale vers une position de destination, il peut être prévu que les deux positions initiales soient diamétralement opposées dans la répartition 10 de référence, et les deux positions finales diamétralement opposées dans la répartition voisine obtenue. A titre d'exemple, la répartition V4 illustrée en figure 6 est obtenue d'après le quatrième critère. Le nombre de permutations circulaires maximal autorisé pour 15 respecter le quatrième critère peut être augmenté à une valeur supérieure à 2. Au cours de l'étape de génération des voisins, plusieurs des quatre critères décrits ci-dessus sont utilisés pour générer au moins une répartition voisine de la répartition de référence. En d'autres termes, si les quatre 20 critères sont utilisés alors au minimum quatre répartitions voisines sont générées au cours d'une même itération. Les quatre critères de ressemblance différents exposés ci-dessus définissent des voisinages différents de la répartition de référence. Le fait d'utiliser différents critères de ressemblance (et donc de 25 définir différents voisinages pour une répartition de référence) est avantageux dans le cas où la fonction de coût dépend de plusieurs types de paramètre d'équilibrage. En effet, l'utilisation d'un seul voisinage, défini à l'aide d'un unique critère de ressemblance pourrait aboutir à la sélection d'une répartition 30 relativement équilibrée du point de vue d'un premier type de paramètre 3059799 16 d'équilibrage, mais relativement déséquilibré d'un point de vue d'un autre type de paramètre d'équilibrage. En conséquence, le fait d'utiliser ces différents critères permet d'explorer l'espace des répartitions possibles selon différentes stratégies 5 possibles. L'algorithme peut être mis en oeuvre avec des contraintes prédéterminées que la répartition d'aubes à choisir doit respecter. Le caractère fixe ou bloqué de telle ou telle aube, déjà abordé précédemment, constitue un exemple de contrainte. 10 Par ailleurs, il peut être choisi qu'au moins un des paramètres de contrainte suivants doive être inférieur à un seuil prédéterminé dans La répartition : - Balourd « Radial » résiduel : résultante des moments radiaux sur la roue ; 15 - Balourd « Radial+Tangentiel » résiduel : résultante des moments radiaux et tangentiels sur la roue ; - Balourd « Radial généré par le couple Axial dans le plan du palier 1 résiduel : résultante des composantes radiales générées par le couple axial dans le plan du palier 1 sur la roue ; 20 - Balourd « Radial+Tangentiel+Radial généré par le couple Axial dans un plan du palier 1 » résiduel : résultante des moments radiaux, tangentiels et radiaux Liés au couple axial sur l'ensemble de la roue ; - Delta moment Radial max aubes face à face : le maximum des écarts de moment entre 2 aubes face à face relevé sur la roue ; 25 - Delta moment Axial max aubes face à face : le maximum des écarts de moment entre 2 aubes face à face relevé sur la roue ; - Delta moment Tangentiel max aubes face à face : le maximum des écarts de moment entre 2 aubes face à face relevé sur la roue ; - Balourd « Calage » résiduel : en assimilant les calages des aubes de 30 la roue à des vecteurs radiaux, on calcule la résultante radiale des calages sur la roue (calage mesuré à une certaine hauteur de coupe) ; 3059799 17 - Delta Calage max aube à aube : le maximum des écarts en calage relevé entre 2 aubes successives sur l'ensemble de la roue (calage mesuré à une certaine hauteur de coupe). En définitive, quatre familles de contraintes peuvent être examinées : 5 - Les balourds : sur chaque position de la roue, une caractéristique d'aube (angle, moment, masse, ...) est représentée par un vecteur. Le balourd est la longueur du déplacement engendré par l'addition de l'ensemble des vecteurs associés à une ou plusieurs caractéristiques d'aube. 10 - Les deltas de caractéristique d'aubes face à face : sur certaines configurations de roue, en face de chaque aube se trouve une aube. Dans ce document, une paire d'aubes, désignera 2 aubes qui sont montées l'une en face de l'autre. Se donner une contrainte sur tes paires consiste à regarder la différence des caractéristiques des 15 aubes appariées. - Les deltas de caractéristique d'aubes adjacentes : la différence maximale de caractéristiques d'une aube avec ses plus proches voisines est une quantité d'intérêt. - Les positions bloquées : l'algorithme d'optimisation ne testera pas 20 de roue avec des aubes bloquées qui ne seraient pas à leur place. Par exemple, « déplacer » une aube, au sens du troisième critère, conduit à déplacer beaucoup d'aubes de leur position initiale dans la configuration de référence, cassant l'appairage initial des aubes dans la 25 configuration de référence. Ceci a pour conséquence d'être incompatible avec la contrainte d'avoir des aubes fixes (c'est-à-dire des aubes indéplaçables) ou, dans une moindre mesure, des contraintes sur des paires d'aubes diamétralement opposées sur la roue. Si un certain nombre d'aubes perdent leur position initiale (généralement déplacée d'une case), leur 30 position vis-à-vis de leur voisine est peu modifiée, du coup ce voisinage privilégie une amélioration sur les critères de type « delta adjacent ». 3059799 18 Le tableau ci-dessous synthétise les avantages et inconvénients des différents voisinages en fonction du critère d'équilibrage recherché. Voisinage Aubes fixes Delta critère adjacent Delta critère paire Balourd déplacement - - + - + 1 aube permutation ++ - - + 2 aubes déplacement - - + + + 1 paire permutation ++ - + + 2 paires 5 Sélection d'une répartition voisine « optimale » Le balourd d'une configuration voisine donnée est calculé à partir des n valeurs de mesure de paramètre d'équilibrage de même type (par exemple n angles de calage) contenues dans chacune des n configurations d'aube, et à partir des n positions associées aux n configurations d'aube de 10 la répartition voisine. Par exemple, le balourd d'une configuration voisine se rapportant aux angles de calage peut être calculé comme la somme de n vecteurs élémentaires ; chaque vecteur élémentaire d'indice i allant de 1 à n est défini à l'aide de la mesure d'angle de calage et est consignée dans la 15 configuration d'aube d'indice i, et la position d'indice j. L'orientation du vecteur somme calculé, par rapport à l'axe virtuel du disque D est représentative d'un degré d'équilibre de la configuration voisine, du point de vue des angles de calage. Un vecteur somme colinéaire à cet axe est représentatif d'une roue parfaitement équilibrée. 20 Des balourds se rapportant à un moment (tangentiel, axial, ou radial) peuvent être calculés de la même manière. Lorsque les configurations d'aube contiennent plus d'un type de paramètre d'équilibrage, plusieurs balourds peuvent être calculés pour chacune des répartitions voisines, chaque balourd se rapportant à un type 25 de paramètre d'équilibrage. Par exemple, pour chacune des répartitions 3059799 19 voisines, un balourd relatif aux angles de calages, un balourd relatif aux moments axiaux, un balourd relatif aux moments tangentiels, etc., peuvent être calculés. D'autres balourds plus complexes, représentatifs d'une combinaison 5 de mesures de paramètre d'équilibrage, peuvent également être calculés. Fonction de coût Dans un mode de réalisation, la fonction de coût dépend des valeurs prédéterminées suivantes: 10 - C_ba : une référence de balourd axial - C_br : une référence de balourd radial - C_bt : une référence de balourd tangentiel - C_brt : une référence de balourd Radial+Tangentiel - C_brta : une référence de balourd Radial+Tangentiel+Axial 15 - C_bang : une référence de balourd angulaire - C_adjacent_ang : une référence d'écart d'angles adjacents - C_paire_a : une référence de critère paire moments axial - C_paire_t : une référence de critère paire moments tangentiel - C_paire_r : une référence de critère paire moments radial 20 - C_paire_ang : une référence de critère paire angle de calage Ces valeurs sont déterminées par un expert selon les seuils de vibration et d'acoustique qu'il souhaite atteindre. La fonction de coût dépend en outre des paramètres suivants : - BA(Roue) : balourd axial calculé pour la répartition de roue aubagée 25 Roue - BR(Roue) : balourd radial calculé pour la répartition de roue aubagée Roue - BT(Roue) : balourd tangentiel calculé pour la répartition de roue aubagée Roue 30 - BRT(Roue) : balourd radial+tangentiel calculé pour la répartition de roue aubagée Roue 3059799 20 - BRTA(Roue) : balourd radial+tangentiel+axial calculé pour la répartition de roue aubagée Roue - BAng(Roue) : balourd angulaire calculé pour la répartition de roue aubagée Roue 5 - Dadja_Ang(Roue) : écart maximum d'angles de calage adjacents trouvé dans la répartition Roue. La fonction de coût est construite à l'aide des termes suivants, appelés « éloignements » : - E_ba(Roue) = max(1, BA(Roue)/C_ba) 10 - E_br(Roue) = max(1, BR(Roue)/C_br) - E_bt(Roue) = max(1, BT(Roue)/C_bt) - E_brt(Roue) = max(1, BRT(Roue)/C_brt) - E_brta(Roue) = max(1, BRTA(Roue)/C_brta) - E_bang(Roue) = max(1, Bang(Roue)/C_bang) 15 - E_delta_ang(Roue) = max(1, Dadja_Ang (Roue)/C_adjacent_ang) Ces éloignements ont la forme illustrée en figure 7. On comprend que chaque référence de l'axe d'abscisse définit un critère d'équilibrage. Quel que soit l'éloignement considéré, lorsque la roue respecte le critère associé, l'éloignement stagne à 1. L'intérêt de cette borne inférieure, ou 20 seuillage, réside dans l'idée que tant qu'un critère est atteint l'algorithme se concentre sur les autres critères. La fonction de coût dépend également d'un autre paramètre Nb_paire(Roue), appelé « nombre de delta paire moment tangentiel hors critère en radial ou en tangentiel », et est défini comme suit : 3059799 21 Nb_paire(Roue) n 2 )(1{1 M_Ri-M_Rj +C paire_r} M_Ti-M_Ti+nl REVENDICATIONS1. Procédé de simulation de répartition d'aubes (P) sur un disque (D) de turbomachine, le procédé comprenant les étapes de : - fourniture d'une pluralité de configurations d'aube, chaque configuration d'aube étant associée à une aube et comprenant au moins une mesure de paramètre d'équilibrage (R, T, A) mesurée sur l'aube (P) associée, - recherche et sélection (106) d'une répartition de roue aubagée associant les configurations d'aubes fournies avec des positions sur le disque, la répartition de roue aubagée favorisant l'atteinte d'au moins un critère défini d'après une fonction de coût prédéterminée dépendant des mesures de paramètres d'équilibrage, la recherche et la sélection étant conduites par itérations successives, le procédé étant caractérisé par le fait qu'une itération courante de la recherche et de la sélection (106) comprend les étapes de: o génération d'une pluralité de répartitions voisines (V1, V2, V3, V4) à partir d'une répartition de référence (REF), ladite pluralité comprenant au moins une première répartition voisine générée selon une stratégie favorisant l'atteinte d'un premier ensemble de critères de ressemblance prédéterminé avec la répartition de référence (REF) et au moins une deuxième répartition générée selon une autre stratégie favorisant l'atteinte d'un deuxième ensemble de critères de ressemblance prédéterminé avec la répartition de référence, o sélection, parmi les répartitions voisines, d'une répartition voisine favorisant une optimisation de la fonction de coût comme nouvelle répartition de référence. 2. Procédé selon la revendication 1, dans lequel la recherche et la sélection sont conduites au moyen d'un algorithme de type Tabou, et dans lequel au 3059799 24 moins une répartition de roue aubagée ayant participée à au moins une sélection au cours d'une itération précédente n'est pas traitée comme une répartition voisine au cours de l'itération courante. 5 3. Procédé selon la revendication 2, dans lequel au moins une répartition voisine (V1, V2) résulte d'un ou deux échanges de positions d'aube entre configurations d'aubes de la répartition de référence (REF). 4. Procédé selon la revendication 3, dans lequel, une des positions 10 échangées dans le cadre d'un des deux échanges et une des positions échangées dans le cadre de l'autre des deux échanges sont diamétralement opposées. 5. Procédé selon l'une des revendications 2 à 4, dans lequel au moins une 15 répartition voisine (V3, V4) résulte d'une ou de deux permutations circulaires de configuration d'aube appliquées à la répartition de référence, chaque permutation circulaire étant opérée sur un sous-ensemble de configurations d'aube associées à des positions consécutives autour du disque dans la répartition de référence. 20 6. Procédé selon la revendication 5, dans lequel chacune des deux permutations circulaires appliquées est représentative d'un déplacement de configuration d'aube depuis une position initiale vers une position de destination, les deux positions initiales étant diamétralement opposées 25 dans la répartition de référence, et les deux positions finales étant diamétralement opposées dans la répartition voisine obtenue. 7. Procédé selon l'une des revendications 2 à 6, dans lequel est sélectivement généré au cours de l'itération courante un nombre de 30 répartitions voisines inférieur au nombre total de premières voisines de la répartition de référence. 3059799 25 8. Procédé selon l'une des revendications 2 à 7, comprenant, pour la répartition de référence et pour chaque répartition voisine, le calcul d'au moins un balourd à partir de mesures de paramètre d'équilibrage d'un 5 même type contenues dans les configurations d'aube de la répartition, et à partir des positions associées aux configurations d'aube dans la répartition, la fonction de coût dépendant du balourd. 9. Procédé selon l'une des revendications 2 à 8, dans lequel la fonction de 10 coût dépend de plusieurs balourds, chaque balourd se rapportant à un type de mesure de paramètre d'équilibrage respectif, et dans lequel au moins un des balourds fait l'objet d'un seuillage par application de la fonction de coût. 15 10. Procédé selon l'une des revendications 2 à 9, dans lequel les mesures de paramètre d'équilibrage sont d'au moins un des types suivants: angle de calage d'aube, moment axial d'aube, moment radial d'aube, et moment tangentiel d'aube. 20 11. Procédé de fabrication d'une roue aubagée pour turbomachine à partir d'aubes et d'un disque, le procédé comprenant les étapes de - mise en oeuvre d'un procédé de simulation selon l'une des revendications précédentes, et - assemblage des aubes sur le disque d'après la répartition de roue 25 aubagée sélectionnée par le procédé de simulation. 12. Produit programme d'ordinateur comprenant des instructions de code de programme pour exécuter les étapes du procédé de simulation selon l'une des revendications 1 à 10, lorsque ce programme est exécuté par une 30 unité de calcul (4).