FR 3039677 A1 20170203 FR 1557301 A 20150730 1 PROCEDE DE CONCEPTION DE PIECES MECANIQUES, NOTAMMENT D'AUBES DE TURBOMACHINE DOMAINE TECHNIQUE GÉNÉRAL ET ART ANTÉRIEUR Dans le cadre de la conception de pièces mécaniques tels que des aubes de turbomachines, on utilise classiquement des techniques d'optimisation paramétrique afin d'identifier le meilleur des compromis entre d'une part un objectif, par exemple un objectif de performance aérodynamique des aubes et d'autre part d'autres contraintes métiers, telles que celles liées aux coûts de fabrication ou de certification (thermique, acoustique, mécanique, etc..), qui sont souvent contradictoires avec l'objectif à atteindre. Pour ce faire, des chaînes de calcul multi-métiers automatisées mettent en oeuvre des modèles de calcul qui s'appuient sur des codes de simulation variés (aérodynamiques, acoustiques, mécaniques, ...) et permettent une optimisation 15 multidisciplinaire. En raison du nombre important de variables géométriques et de contraintes, le problème d'optimisation qui en découle est particulièrement complexe. Compte tenu du coût de calcul unitaire induit par les codes de simulations chargés du calcul des réponses des modèles, le nombre de calculs à réaliser pour explorer largement 20 l'espace de conception rend illusoire l'emploi de stratégies d'optimisation paramétrique globale et directe dans un délai compatible des cycles de conception moteur (de type algorithme génétique par exemple). A titre illustratif sur un problème type de conception d'aubes de turbomachines, les coûts de calcul des simulations aérodynamiques et mécaniques, 25 peuvent représenter 3 à 5h de calcul pour une géométrie donnée. Pour résoudre cette problématique, il est classique, dans des stratégies d'optimisation globale indirecte automatisée, d'utiliser des surfaces de réponse (« estimateurs » ou « méta-modèles » dans la suite du présent texte) qui sont des représentations mathématiques simplifiées des réponses des calculs à réaliser 30 (aussi appelées « sorties ») dont le coût unitaire d'évaluation est négligeable. Des exemples en ce sens, utilisant des méta-modèles de Krigeage ou des méta-modèles de type RBF (« Radial Basis Function » selon la terminologie anglo-saxonne généralement utilisée), sont par exemple décrits dans les publications suivantes : 3039677 2 [1] M.J. Sasena. Flexibilty and efficiency enhancements for constrained global design optimization with Kriging approximations.PhD thesis, University of Michigan, 2002. [2] Donald R. Jones, Matthias Schonlau and William J. Welch. Efficient 5 Global Optimization of expensive black-box functions.Journal of Global Optimization, 13(4):455-492, December 1998 [3] Rommel G. Regis. Constrained optimization by radial basis function interpolation for high-dimensional expensive black-box problems with infeasible initial points. Engineering Optimization, 00(0) : 1-28,2012. 10 [4] "Multidisciplinary Design Optimization with Minamo" : http://www.cenaero.be/Page_Generale.asp?Docl D=15336&1a=1&langue=EN Les méta-modèles de réponse sont généralement construits de façon itérative par apprentissage en s'appuyant sur une sélection automatique et intelligente des points de l'espace de conception qui vont permettre de maximiser la 15 connaissance des réponses des calculs et la qualité prédictive des méta-modèles. L'optimisation s'arrête lorsque l'algorithme a convergé vers un optimum jugé global. Cette approche présente plusieurs inconvénients. Le coût de construction d'un méta-modèle peut ne pas être négligeable, 20 notamment si la base de données d'apprentissage et le nombre de variables sont importants : le calcul des méta-modèles de dizaines d'objectifs et contraintes peut donc s'avérer très coûteux à chaque itération de l'optimisation ; c'est le cas en particulier des modèles de Krigeage qui peuvent présenter des limitations importantes lorsqu'il s'agit d'inverser des matrices de covariance de taille 25 significative. Par ailleurs, la qualité des méta-modèles (c'est-à-dire leur capacité prédictive), pour être acceptable, peut nécessiter un grand nombre de calculs de simulation. En outre, cette approche est très mal adaptée au cas de problèmes très 30 contraints présentant des domaines de solutions faisables de taille très restreinte par rapport à l'espace de recherche, entraînant un coût calcul excessif en phase exploratoire, parfois sans que cela ne débouche sur une solution - on entend par solution faisable une solution respectant toutes les contraintes métiers à respecter. PRÉSENTATION GÉNÉRALE DE L'INVENTION 3039677 3 Un but général de l'invention est de pallier ces inconvénients et de proposer une solution fiable, qui converge rapidement et qui soit particulièrement peu coûteuse en temps de calcul. Un but de l'invention est en particulier de proposer un procédé de conception 5 d'une pièce mécanique permettant de déterminer de manière rapide et peu coûteuse en temps de calcul un dimensionnement d'une pièce représentant un optimum par rapport à un objectif à atteindre, tout en respectant un ensemble de contraintes. 10 A cet égard, l'invention a pour objet un procédé de conception d'une pièce mécanique définie par un ensemble de variables de définition prenant des valeurs dans un espace de conception, pour déterminer un dimensionnement de la pièce optimisé du point de vue d'un objectif à atteindre, tout en respectant un ensemble de contraintes, le procédé comprenant la mise en oeuvre d'au moins une itération de la succession des étapes: a) à partir d'un plan d'expérience comprenant au moins un dimensionnement d'une pièce, chaque dimensionnement comprenant des valeurs pour l'ensemble des variables de définition, et, pour chaque dimensionnement, une valeur de fonction objectif correspondante, construire un méta-modèle de krigeage pour chaque contrainte et pour la fonction objectif, b) évaluer les méta-modèles de Krigeage sur l'espace de conception, et déterminer un dimensionnement de prédiction optimale d'atteinte de l'objectif et de respect des contraintes, c) évaluer le respect des contraintes et un degré d'atteinte de l'objectif par ledit dimensionnement, et enrichir le plan d'expérience précédent avec ledit dimensionnement et les valeurs obtenues par l'évaluation, le procédé étant caractérisé en ce que la construction des méta-modèles de Krigeage comprend la mise en oeuvre d'une régression par moindres carrés partiels 30 pour diminuer le nombre de paramètres dans la construction de chaque méta-modèle. Avantageusement, mais facultativement, le procédé selon l'invention peut en outre comprendre au moins l'une des caractéristiques suivantes : l'étape de construction des méta-modèles de Krigeage peut comprendre l'utilisation, 35 pour l'expression de prédicteurs de Krigeage, d'une matrice de corrélation 3039677 4 calculée à partir d'un nombre de paramètres de pondération strictement inférieur au nombre de variables de définition. - l'étape de construction des méta-modèles de Krigeage peut comprendre l'utilisation d'une matrice de corrélation calculée à partir d'une fonction de 5 corrélation définie par : d = exp 2 lui 02) =1 i=t où u et y sont deux vecteurs de coordonnées u,, y, avec i entier compris entre 1 et d, d est le nombre de variables de définition, W est une matrice contenant des coefficients de pondération des variables initiales dans la construction des variables latentes , 6, est un paramètre de pondération de la variable de définition j, et q est le 10 nombre réduit de dimensions. - le procédé peut en outre comprendre, lors de la première itération, avant la mise en oeuvre de l'étape b), la mise en oeuvre d'une étape de recherche d'au moins un dimensionnement satisfaisant à l'ensemble des contraintes. - l'étape de recherche d'un dimensionnement satisfaisant aux contraintes peut 15 comprendre : d) l'évaluation des méta-modèles de Krigeage sur l'espace de conception pour déterminer un dimensionnement présentant une probabilité maximale de respect des contraintes, e) l'évaluation du respect des contraintes et du degré d'atteinte de 20 l'objectif par ledit dimensionnement, et l'enrichissement du plan d'expérience précédent avec ledit dimensionnements et les valeurs obtenues par l'évaluation, et f) la mise à jour des méta-modèles de Krigeage à partir du plan d'expérience enrichi. 25 - les étapes d) à f) peuvent être mises en oeuvre de manière itérative jusqu'à déterminer un dimensionnement respectant l'ensemble des contraintes, ou jusqu'à un nombre maximal d'itérations pré-déterminé. - l'étape b) peut être mise en oeuvre en évaluant les méta-modèles de Krigeage générés sur la base du dernier plan d'expérience enrichi. - lors de la première itération, le plan d'expérience initial peut être un plan d'expérience pré-calculé ou un plan de type Latin Hypercube. 3039677 5 les étapes a) à c) peuvent être itérées jusqu'à atteindre un nombre d'itérations pré-déterminé, et un dimensionnement optimal de la pièce est le dernier dimensionnement ayant enrichi le plan d'expérience à l'étape c). 5 L'invention a également pour objet un produit programme d'ordinateur, comprenant des instructions de code pour la mise en oeuvre du procédé selon la description qui précède, lorsqu'il est exécuté par un processeur. L'invention a enfin pour autre objet une unité de traitement, comprenant un 10 processeur configuré pour la mise en oeuvre du procédé selon la description qui précède. Le procédé de conception selon l'invention permet de réduire de façon importante le temps de calcul pour générer les méta-modèles de Krigeage. En effet, 15 le couplage d'une méthode de Krigeage avec une méthode de régression par moindres carrés partiels permet de diminuer de façon importante le nombre de paramètres de réglage du méta-modèle de Krigeage. En particulier, une étape de recherche de dimensionnement satisfaisant à l'ensemble des contraintes peut être mise en oeuvre, en recalculant à chaque 20 nouveau dimensionnement candidat les méta-modèles de Krigeage, sans nécessiter un temps de calcul supplémentaire important. Cette étape est donc rapide à mettre en oeuvre et accélère ensuite la convergence des étapes de recherche du dimensionnement optimal. 25 Le gain de temps ainsi dégagé, sans pour autant dégrader la qualité des méta-modèles de Krigeage, permet de mieux respecter les cycles de conception et/ou d'affiner la recherche du concept optimal. La solution proposée permet en outre de lever les difficultés techniques liées à la construction d'un méta-modèle de Krigeage lorsque la base d'apprentissage et 30 la dimension du problème sont de taille importante. PRÉSENTATION DES FIGURES D'autres caractéristiques et avantages de l'invention ressortiront encore de la description qui suit, laquelle est purement illustrative et non limitative, et doit être lue en regard des figures annexées sur lesquelles : 3039677 6 - la figure 1 est une représentation schématique d'un profil d'aubage de turbomachine sur lequel on a porté différentes variables géométriques intervenant dans sa définition ; - la figure 2 est un organigramme illustrant différentes étapes d'un traitement 5 d'optimisation paramétrique conforme à un mode de mise en oeuvre de l'invention. - les figures 3a et 3b illustrent schématiquement la mise en oeuvre d'une analyse par composantes principales. DESCRIPTION D'UN OU PLUSIEURS MODES DE MISE EN OEUVRE ET DE RÉALISATION 10 En référence à la figure 1, on a représenté un exemple de section d'une pièce - en l'occurrence la pièce est une aube de turbomachine -susceptible d'être optimisée au moyen du procédé proposé. Sur la figure 1 sont également représentées différentes variables géométriques définissant le dimensionnement de la pièce, ci-après appelées variables de définition, caractéristiques du profil de cette 15 pièce et destinées à constituer des paramètres de simulation pour celles-ci : positions du bord d'attaque BA et du bord de fuite BF, position du centre de gravité CG, épaisseur maximale EPMax, longueur de corde, angle y de la ligne de corde, anglef31 d'entrée d'air au niveau du bord d'attaque BA, angle f32 de sortie d'air au niveau du bord de fuite BF. 20 Le nombre de variables de définition peut être très important. Dans l'exemple qui précède, les variables de définition sont définies sur plusieurs niveaux de hauteur d'une aube par rapport au pied de l'aube. Le nombre de variables total augmente donc presque linéairement avec le nombre de sections d'aube considérées. 25 L'optimisation d'une telle aube se fait par exemple avec un objectif général consistant par exemple à maximiser la performance aérodynamique de la pale, sous un certain nombre de contraintes aérodynamiques et mécaniques. Par exemple, l'objectif général peut consister à optimiser le rendement isentropique 11DP des points DP de conception correspondant aux différentes 30 variables entrant dans la définition de l'aubage dans l'espace de conception DS, sous différentes contraintes aérodynamiques (contraintes relatives par exemple aux débits massiques d'air ThDp, au rendement isentropique à la marge de poussée Tism, à la pression 1157,4 à la marge de poussée, etc..), ainsi que sous différentes 3039677 7 contraintes mécaniques (contraintes relatives aux déplacements de l'aubage au niveau de son épaulement de contact, aux déplacements de l'aubage au niveau de ses bords de fuite BF et d'attaque BA, aux valeurs de contraintes de Von Mises, aux différentiels entre les harmoniques du rotor et les fréquences propres des 5 différentes aubages, etc.). Pour une même aube, le nombre de variables de définition peut être de l'ordre d'une ou plusieurs dizaines, voire de l'ordre de la centaine. Par exemple, une aube peut comporter environ 30 variables de définition. Dans la suite, on note d le nombre de variables de définition. 10 Le nombre de contraintes est également élevé : plusieurs dizaines, voire plusieurs centaines. Dans la suite, on note m le nombre de contraintes. On note en outre mathématiquement x un dimensionnement donné pour une pièce mécanique, x étant donc un vecteur de taille d : (x1,...,xd), où les x, sont les valeurs de chacune des variables de définition. 15 On note également g, une fonction mathématique représentative du respect d'une contrainte i (i entre 1 et m). Le respect de chaque contrainte par un dimensionnement x donné peut être ramené à une expression d'inégalité du type gi(x) O. En effet, pour illustrer, une contrainte g, peut être le respect d'une valeur de 20 masse maximale de la pièce. La fonction contrainte correspondante est dans ce cas formulée comme « masse (x) - masse maximale ». On observe que la fonction évaluée en x est négative si la masse de la pièce pour le dimensionnement x est effectivement inférieure à la masse maximale. Enfin, on note y une valeur correspondant à un degré d'atteinte de l'objectif, 25 et f la fonction qui associe à un dimensionnement x donné la valeur y (un scalaire) représentative du degré d'atteinte d'un objectif, on note f(x)=y. Pour reprendre l'exemple qui précède, y peut correspondre à une valeur rendement isentropique II Dp pour un dimensionnement donné. 30 Le traitement d'optimisation illustré sur la figure 2 met en oeuvre une stratégie d'optimisation adaptative automatisée assistée par méta-modèles de Krigeage. En reprenant les notations qui précèdent, le problème à résoudre est un problème d'optimisation mono-objectif paramétrique continue et sous contrainte et il 35 est formulé comme suit : 3039677 8 /min y = min f (x) x x avec x E Rd [gi(x) Ohi = 1.. m} Ce procédé est mis en oeuvre au moyen d'un ou plusieurs ordinateurs, par exemple par un processeur (non représenté) configure pour exécuter un programme approprié. 5 Dans une première étape (étape 1), on renseigne le processeur avec un plan d'expérience comprenant n dimensionnements x(i) d'une pièce, i=1..n, avec n entier, l'ensemble des dimensionnements étant noté (x(1), ...,x(n)), et pour chaque dimensionnement la valeur y(i)=f(x(i)) correspondante de degré d'atteinte de l'objectif et les valeurs gi(e) de respect de contraintes préalablement calculées pour ces 10 dimensionnements. Le calcul des gi(e) et des y(i) est coûteux en temps de calcul puisqu'il fait appel à des modèles exacts de simulation. Chaque point x(i) est un vecteur de d valeurs (où d est un nombre entier) correspondant aux différentes variables entrant dans la définition de la pièce. Ce plan d'expérience est par exemple déterminé à partir d'un 15 échantillonnage de type « Latin Hypercube ». Alternativement, le plan d'expérience peut être un plan ayant été préalablement calculé et chargé par le processeur pour la mise en oeuvre du procédé. A partir du plan d'expérience ainsi renseigné, le processeur construit par 20 traitement de Krigeage un méta-modèle de réponse (étape 2) pour chacune des sorties g, (fonctions contraintes) et f (fonction objectif). On note 3î la prédiction d'un méta-modèle de Krigeage pour la valeur de y, et i=1..m les prédictions des méta-modèles de Krigeage respectifs pour les résultats des fonctions g,. On note également yx=f(x) la valeur réelle de la fonction f 25 en x. La construction des méta-modèles met en oeuvre une régression par moindre carrés partiels (« Partial Least Square » ou « PLS » selon la terminologie anglo-saxonne généralement utilisée). Dans le cadre d'un traitement de Krigeage, une réalisation d'un processus 30 gaussien (le résultat d'un tirage selon une loi de probabilité normale), notée par exemple A(x) est décomposée en tout point x en la somme d'un prédicteur de Krigeage p.(x) et d'une erreur Z(x) (confère paragraphe 3.2 de la publication [1] de 3039677 9 Sasena citée ci-avant) distribuée selon une loi normale .7V'(0, 0-2) où la variance 0-2 est un paramètre pouvant être estimé : A(x) = u(x) + Z(x) ,u (x) =1,g jhj(x) j=i Avec un nombre s de fonctions de base h,.(x). Le Krigeage repose sur la supposition si deux points u et y sont proches, les 5 erreurs respectives Z(u) et Z(v) seront proches également, c'est-à-dire que l'erreur Z(x) n'est pas indépendante. Une corrélation ;1,, entre deux points u et y est exprimée par l'équation : d ruv=flexp(-Bjluj-vj 1=1 où 01 est un paramètre mesurant l'effet de la direction j du domaine de recherche sur la variation de la sortie (plus 01 est important, et plus le terme 10 exponentiel correspondant est petit), et pi mesure un caractère plus ou moins lisse de l'exponentielle dans la direction]. On a une matrice de corrélation : R = (rx(1),..,rx(n)) avec rx = (rx,x(1),..rxx(n))t Alors le meilleur prédicteur de Krigeage non biaisé de la valeur y (x'), où x' 15 est un dimensionnement pour lequel la valeur réelle yx'est inconnue, est : = h(xnv)t/ - HP) Avec h le vecteur contenant l'ensemble des fonctions de régression, fi l'estimation des paramètres de régression et y le vecteur contenant les réponses, définis comme suit : h(x) = (hi(x), hs(x))t h(x(1))t\ H= Vi(x(n))t F'1 /3 = (1) = (HtR-1H)-1HtR 1y fis PI) 3039677 10 Y = 1) 37x( YX (n) La difficulté de construire un tel méta-modèle de Krigeage lorsque le nombre de dimensions d des variables est élevé provient de la difficulté d'estimer les paramètres O dans l'expression de la corrélation. C'est typiquement le cas lors de 5 l'application des méta-modèles de Krigeage à la conception de pièces telles que des aubes, qui comprennent plusieurs dizaines de variables de définition (correspondant au nombre d de dimensions). Pour remédier à cette problématique, la technique de Krigeage est couplée avec la méthode des moindres carrés partiels PLS (pour l'anglais Partial Least 10 Squares). La méthode des moindres carrés partielles permet de projeter un ensemble de variables initiales sur un ensemble plus réduit de variables appelées variables latentes (ou composantes principales) qui sont obtenues par combinaison linéaire des variables initiales : 15 T=XW W est une matrice contenant le vecteur de poids, c'est-à-dire des coefficients de pondération des variables initiales dans la construction des variables latentes (en anglais, « loading vectors ») de d lignes et q colonnes, ou q est très inférieur à d ; typiquement inférieur ou égal à 4 (par exemple égal à 2, 3 ou 4). q correspond au 20 nombre de variables latentes retenues. Les variables latentes sont des vecteurs formant les colonnes de la matrice T de n lignes et q colonnes, où n est toujours le nombre de dimensionnement du plan d'expérience initial. X est la matrice contenant les dimensionnements du plan d'expérience 25 initial. En référence à la figure 3a, après avoir centré et normalisé la matrice X et le vecteur y de degré d'atteinte d'objectif correspondant, la première variable latente est calculée en maximisant la covariance au carré entre la première variable latente XW( -0 et le vecteur y- où W( désigne la emecolonne de la matrice W: (max wi Xt yyt Xwi W(: 1) = t. q. wfwi = 1 3039677 11 La maximisation est obtenue lorsque w1 est le vecteur propre de la matrice XtyytX correspondant à la plus grande valeur propre. La plus grande valeur propre de ce système peut être estimée par la méthode itérative décrite dans : [5] C. Lanczos, « An iteration method for the solution of the eigenvalue 5 problem of linear differential and integral operators », Journal of Research of the National Bureau of Standards, 45(4) :255-282, Octobre 1950. Une fois W( nobtenu les matrices résiduelles de X et de y sont calculées, notées Xil:' et 37,1-'respectivement : fXres = X - X We Yrés = Y - clXVIT(,t) Où piest le vecteur des coefficients de régression de la régression locale de 10 X sur la première variable latente et c1 est le coefficient de régression de la régression de y sur la première variable latente. Puis la deuxième variable latente est orthogonale à la première et peut être séquentiellement calculée en résolvant le même système que précédemment en remplaçant respectivement X et y par X,1.-es et 37,1-es. 15 Ceci est mis en oeuvre itérativement jusqu'à obtenir le nombre q de variables latentes correspondant au nombre de dimensions réduits à obtenir. La combinaison du Krigeage avec la méthode des moindres carrés partiels consiste à substituer la fonction de corrélation spatiale classique décrite ci-avant par une fonction de corrélation dans laquelle les variables sont de dimension réduite 20 grâce à la méthode PLS. L'expression de la nouvelle fonction de corrélation est, en reprenant les mêmes notations que précédemment : d 1 ru,' = exp (-0iI1W1,11 lui 02) =1 i=i Dans cette équation 1W01 correspond à l'important de la ième variable initiale (i=1,...,d) pour la construction de la jème variable latente (j=1,...,q) et Bi traduit 25 l'importance de la jème variable latente sur la variable de sortie. Par conséquent on constate que les variables latentes transmettent l'influence des variables initiales sur la variable de sortie. De plus, on constate que le nombre de paramètres Bipasse de d à q et qu'il est donc considérablement réduit, ce qui réduit le nombre de calculs à mettre en 30 oeuvre et donc augmente la vitesse de construction des méta-modèles de Krigeage. Typiquement, q peut être limité à 2, 3 ou 4. 3039677 12 Enfin les propriétés d'une matrice de corrélation classique Gaussienne sont conservées avec cette nouvelle expression de la fonction de corrélation. Une fois les méta-modèles des réponses ainsi construits, le procédé peut 5 comprendre une étape facultative (étape 3) de recherche d'au moins un point de faisabilité. Cette étape est en effet mise en oeuvre si et seulement si le plan d'expérience initial ne comporte aucun point faisable. Par point faisable, on entend un dimensionnement x satisfaisant toutes les contraintes, c'est-à-dire un dimensionnement pour lequel gi(x) 0 pour tout i de 1 10 à m. A cet effet, en utilisant les méta-modèles de Krigeage obtenus à l'étape précédente, et dont le coût unitaire d'évaluation en un point x est négligeable, le processeur recherche au cours d'une sous-étape 3a, parmi tout le domaine de 15 conception - et donc pour tous les dimensionnements x possibles de ce domaine et pas seulement ceux du plan d'expérience initial - le dimensionnement x maximisant la probabilité de respect de l'ensemble des contraintes : max nnei(x) 200 expériences). En outre, l'utilisation du critère d'enrichissement LWB2 au lieu du critère El 15 habituel, permet une convergence plus rapide tout en assurant une exploration efficace de l'espace de conception (le critère LWB2 étant un critère plus local que le critère El). Une telle méthode est particulièrement avantageuse pour la conception d'aubage amélioré de turbomachines mais elle peut également s'appliquer à 20 d'autres types de pièces que des aubages. Elle trouve en particulier avantageusement application dans le cas de conception de pièces de grande taille, très contraintes en termes de certification de leur dimension. 25 REVENDICATIONS1. Procédé de conception d'une pièce mécanique définie par un ensemble de variables de définition prenant des valeurs dans un espace de conception, pour déterminer un dimensionnement de la pièce optimisé du point de vue d'un objectif à atteindre, tout en respectant un ensemble de contraintes, le procédé étant mis en oeuvre par une unité de traitement et comprenant la mise en oeuvre d'au moins une itération de la succession des étapes: a) à partir d'un plan d'expérience comprenant au moins un dimensionnement d'une pièce, chaque dimensionnement comprenant des valeurs pour l'ensemble des variables de définition, et, pour chaque dimensionnement, une valeur de fonction objectif correspondante, construire (2, 4c) un méta-modèle de krigeage pour chaque contrainte et pour la fonction objectif, b) évaluer les méta-modèles de Krigeage sur l'espace de conception, et déterminer (4a) un dimensionnement de prédiction optimale d'atteinte de l'objectif et de respect des contraintes, c) évaluer le respect des contraintes et un degré d'atteinte de l'objectif par ledit dimensionnement, et enrichir (4b) le plan d'expérience précédent avec ledit dimensionnement et les valeurs obtenues par l'évaluation, le procédé étant caractérisé en ce que la construction des méta-modèles de Krigeage comprend la mise en oeuvre d'une régression par moindres carrés partiels pour diminuer le nombre de paramètres dans la construction de chaque méta-modèle. 2. Procédé selon la revendication 1, dans lequel l'étape de construction (2, 3c, 4c) des méta-modèles de Krigeage comprend l'utilisation, pour l'expression de prédicteurs de Krigeage, d'une matrice de corrélation calculée à partir d'un nombre de paramètres de pondération strictement inférieur au nombre de variables de définition. 3039677 18 3. Procédé selon la revendication 2, dans lequel l'étape de construction (2, 3c, 4c) des méta-modèles de Krigeage comprend l'utilisation d'une matrice de corrélation calculée à partir d'une fonction de corrélation définie par : d ,2. lllj=1111111 i=1 où u et y sont deux vecteurs de coordonnées u' v' avec i entier compris entre 1 et 5 d, d est le nombre de variables de définition, W est une matrice contenant des coefficients de pondération des variables initiales dans la construction des variables latentes , 0, est un paramètre de pondération de la variable de définition j, et q est le nombre réduit de dimensions. 10 4. Procédé selon l'une des revendications 1 à 3, comprenant en outre, lors de la première itération, avant la mise en oeuvre de l'étape b), la mise en oeuvre d'une étape (3) de recherche d'au moins un dimensionnement satisfaisant à l'ensemble des contraintes. 15 5. Procédé selon la revendication 4, dans lequel l'étape de recherche d'un dimensionnement satisfaisant aux contraintes comprend : d) l'évaluation (3a) des méta-modèles de Krigeage sur l'espace de conception pour déterminer un dimensionnement présentant une probabilité maximale de respect des contraintes, 20 e) l'évaluation (3b) du respect des contraintes et du degré d'atteinte de l'objectif par ledit dimensionnement, et l'enrichissement du plan d'expérience précédent avec ledit dimensionnements et les valeurs obtenues par l'évaluation, et f) la mise à jour (3c) des méta-modèles de Krigeage à partir du plan 25 d'expérience enrichi. 6. Procédé selon la revendication 5, dans lequel les étapes d) à f) sont mises en oeuvre de manière itérative jusqu'à déterminer un dimensionnement respectant l'ensemble des contraintes, ou jusqu'à un nombre maximal d'itérations pré-30 déterminé. 3039677 19 7. Procédé selon l'une des revendications 4 ou 5, dans lequel l'étape b) est mise en oeuvre en évaluant les méta-modèles de Krigeage générés sur la base du dernier plan d'expérience enrichi. 5 8. Procédé selon l'une des revendications précédentes, dans lequel, lors de la première itération, le plan d'expérience initial est un plan d'expérience pré-calculé ou un plan de type Latin Hypercube. 9. Procédé selon l'une des revendications précédentes, dans lequel les étapes a) à 10 c) sont itérées jusqu'à atteindre un nombre d'itérations pré-déterminé, et un dimensionnement optimal de la pièce est le dernier dimensionnement ayant enrichi le plan d'expérience à l'étape c). 10. Produit programme d'ordinateur, comprenant des instructions de code pour la 15 mise en oeuvre du procédé selon l'une des revendications précédentes, lorsqu'il est exécuté par un processeur. 11. Unité de traitement, comprenant un processeur configuré pour la mise en oeuvre du procédé selon l'une des revendications 1 à 9. 20