FR 3008507 A1 20150116 FR 1356732 A 20130709 DOMAINE TECHNIQUE GENERAL La présente invention concerne la conception assistée par ordinateur. Plus précisément, elle concerne un procédé de modélisation d'une surface non-axisymétrique d'une pièce telle qu'un moyeu ou un carter. ETAT DE L'ART La nécessité d'amélioration constante des performances des équipements, en particulier aéronautiques, par exemple des rotors de turboréacteurs (c'est-à-dire l'ensemble formé d'un moyeu sur lequel sont fixées des aubes s'étendant radialement, tel que visible sur la figure la), a aujourd'hui imposé l'utilisation d'outils informatiques de modélisation. Ces outils permettent d'aider à concevoir des pièces en optimisant de façon automatisée certaines de leurs caractéristiques par l'exécution d'un grand nombre de calculs de simulation. La recherche d'un optimum géométrique aéromécanique et/ou acoustique sur les rotors conduit aujourd'hui à l'obtention de moyeux présentant une paroi localement non-axisymétrique (c'est-à-dire qu'une coupe selon un plan perpendiculaire à l'axe de rotation n'est circulaire) au niveau de la veine, c'est-à-dire l'ensemble des canaux entre les aubes pour l'écoulement de fluide (en d'autres termes les sections inter-aubes), au vu des conditions particulières qui y règnent. La veine non-axisymétrique définit une surface globalement annulaire d'un espace tridimensionnel (une « tranche » du moyeu). De façon générale, pour optimiser une caractéristique d'une pièce, il est nécessaire de paramétriser une loi (en d'autres termes une ou plusieurs courbes décrivant la valeur d'une grandeur physique selon une courbe le long de la pièce, dans un environnement donné) que l'on cherche à optimiser, c'est-à-dire d'en faire une fonction de N paramètres d'entrée. L'optimisation consiste alors à faire varier (en général aléatoirement) ces différents paramètres sous contrainte, jusqu'à déterminer leurs valeurs optimales pour un critère prédéterminé. Une courbe « lissée » est ensuite obtenue par interpolation à partir des points de passage déterminés. Le nombre de calculs nécessaires est alors directement lié (linéairement voire exponentiellement) au nombre de paramètres d'entrée du problème. De nombreuses méthodes de paramétrisation d'une loi existent, et on peut en particulier distinguer deux grandes catégories : - Modèle discret : la loi est définie par la position d'une pluralité de points, déplacés un à un lors de l'optimisation ; - Modèle paramétrique : la loi est définie via des courbes mathématiques connues dans la littérature, telles que les courbes de Bézier ou les courbes NURBS (B-splines rationnelles non uniformes). Dans le cas précis des veines non-axisymétriques, des méthodologies spécifiques de paramétrisation ont été développées. On peut distinguer une famille de méthodes, dites à modification par points d'attraction au niveau de la paroi, i.e. définition de creux et de bosses, qui consiste à partir d'une surface axisymétrique, et à la faire apparaître des creux et des bosses. Alternativement, il a été proposé des méthodes dites à modification par caractérisation de « sections ». On dessine alors chaque section de passage du fluide entre deux aubes, et l'évolution de la paroi selon ces 25 sections. En d'autres termes on paramétrise des lois selon une pluralité de chemins d'écoulement du fluide entre les aubes (en reprenant la méthodologie générale expliquées ci-avant), de sorte à définir une pluralité de surfaces élémentaires correspondant à la veine sur tout le pourtour de la 30 paroi. Bien que plus intéressantes du point de vue aérodynamique, les méthodes par caractérisation de section sont bien plus complexes et peu utilisées. En effet, comme l'on voit sur la figure lc qui sera décrite plus loin, une des difficultés qui en résulte porte sur la continuité entre les surfaces élémentaires, en particulier à l'amont et à l'aval de l'aube. En effet, la section peut n'être pas la même à l'intrados et à l'extrados de l'aube. Ce qui implique une « marche » problématique à l'amont et à l'aval de l'aubage. Il n'existe aujourd'hui pas de méthode par caractérisation de sections 10 qui permette de conserver une géométrie aéromécaniquement saine (i.e. sans aucune singularité). Il serait souhaitable de trouver une nouvelle façon de modéliser une surface non-axisymétrique dans sa globalité, répondant aux exigences de continuité de la paroi, de respect des contraintes de fabrication, et ce tout 15 en laissant à l'utilisateur une importante marge de manoeuvre dans le design. PRESENTATION DE L'INVENTION 20 La présente invention propose selon un premier aspect un procédé de modélisation d'au moins une partie d'une surface non-axisymétrique d'une portion d'une pièce, la portion étant limitée par un premier et un deuxième plan extrémal, le procédé étant caractérisé en ce qu'il comprend la mise en oeuvre, par des moyens de traitement de données d'un 25 équipement, d'étapes de : (a) Modélisation d'au moins deux surfaces élémentaires nonaxisymétriques composant ladite surface, chacune des surfaces élémentaires s'étendant entre les deux plans extrémaux de la portion et étant définie par au moins une 30 courbe de construction s'étendant selon un plan sensiblement parallèle aux plans extrémaux, une première et une deuxième surface élémentaire étant juxtaposées de sorte à présenter une jonction au moins Cl au niveau des plans extrémaux et au moins une discontinuité ; (b) Construction d'au moins une courbe de raccord de classe Cl définie par un premier et un deuxième point extrémal respectivement disposés sur la première et la deuxième surface élémentaire, ladite courbe de raccord étant tangente à la première et la deuxième surface respectivement en le premier et le deuxième point extrémal ; (c) Modification locale des première et deuxième surfaces élémentaires de sorte à respecter ladite courbe de raccord au voisinage de ladite discontinuité ; (d) Restitution sur une interface dudit équipement de l'au moins une partie de la surface obtenue. Le présent procédé correspond à une méthode par caractérisation de sections dans laquelle on ne cherche pas éviter les discontinuités, bien au contraire. Ces dernières sont obtenues, puis corrigées a posteriori sans modifier la géométrie de la surface non-axisymétrique, de sorte à garder tous les avantages de ces méthodes (qualités aérodynamiques et importante marge de manoeuvre), sans les désavantages (singularités). En outre, le procédé est particulièrement optimisé pour améliorer au mieux la qualité des modèles tout en diminuant la puissance de calcul requise. Selon d'autres caractéristiques avantageuses et non limitatives : - une pluralité d'aubes est disposée sur la circonférence de la portion de pièce, chaque surface élémentaire s'étendant entre deux aubes, de sorte que les première et deuxième surfaces élémentaires sont disposées de part et d'autre d'une aube, respectivement sur l'intrados et l'extrados ; - la jonction entre les première et deuxième surfaces élémentaires présente une première discontinuité s'étendant du premier plan extrémal au bord d'attaque de l'aube et une deuxième discontinuité s'étendant du bord de fuite de l'aube au deuxième plan extrémal ; - une courbe de raccord est construite à l'étape (b) pour chacune des première et deuxième discontinuités, chaque courbe de raccord étant tangente à l'aube au niveau du bord de fuite ou du bord d'attaque ; - une courbe de raccord passe par un point milieu disposé sur le bord de fuite ou du bord d'attaque de l'aube de sorte que la distance entre le point milieu et la première surface élémentaire est égale à la distance entre le point milieu et la deuxième surface élémentaire ; - une courbe de raccord est une spline entièrement définie par les points extrémaux et le point milieu ; - l'étape (c) comprend pour chacune des première et deuxième discontinuité la construction d'une surface de raccord présentant en tout point de son bord une continuité de tangence avec la première ou la deuxième surface élémentaire ; - chaque surface élémentaire présente au moins une courbe de construction disposée entre le bord d'attaque et le bord de fuite de l'aube, chacune des deux surfaces de raccord s'étendant sur chacune des première et deuxième surfaces élémentaires d'un bord extrémal jusqu'à la première courbe de construction au-delà de la courbe de raccord ; - l'étape (a) comprend pour chaque courbe de construction de chaque surface élémentaire la mise en oeuvre de sous-étapes de : (ai) Paramétrisation de la courbe de construction en tant que courbe de classe Cl représentant la valeur d'une grandeur physique caractéristique de ladite surface en fonction d'une position le long d'une largeur de la surface élémentaire, la courbe étant définie par : - Deux points de contrôle extrémaux, respectivement sur chacune des deux aubes entre lesquelles ladite surface élémentaire s'étend, définissant l'étendue de ladite largeur de la surface élémentaire ; - Au moins une spline ; la paramétrisation étant mise en oeuvre selon un ou plusieurs paramètres définissant au moins un des points de contrôle extrémaux ; (a2) Détermination de valeurs optimisées desdits paramètres de ladite courbe. - chaque courbe de construction de chaque surface élémentaire est également définie par au moins un point de contrôle intermédiaire disposé entre les points extrémaux, et au moins deux splines connectées au niveau dudit point de contrôle intermédiaire ; la paramétrisation étant également mise en oeuvre selon un ou plusieurs paramètres définissant ledit point de contrôle intermédiaire ; - le ou les paramètres définissant un point de contrôle intermédiaire ou extrémal définissant une courbe de construction sont choisis parmi une abscisse du point, une ordonnée du point, une orientation de tangente à la courbe au niveau du point et au moins un coefficient de tension associé à une demi-tangente à la courbe au niveau du point ; - un coefficient de tension associé à une demi-tangente à la courbe au niveau d'un point intermédiaire ou extrémal est fonction d'une valeur de dérivée seconde de la courbe au voisinage du point ; - chaque spline d'une courbe de construction est une courbe de Bézier 20 entièrement déterminée par les points de contrôle définissant ses extrémités. Selon un deuxième puis un troisième aspect, l'invention concerne un procédé de fabrication d'une pièce dont une portion présente une surface 25 non-axisymétrique, le procédé comprenant des étapes de : - Mise en oeuvre du procédé selon le premier aspect de sorte à modéliser au moins une partie de ladite surface nonaxisymétrique de la portion de la pièce ; - Fabrication de ladite pièce conformément à la modélisation de 30 l'au moins une partie de la surface obtenue ; Ainsi qu'une pièce dont une portion présente une surface non-axisymétrique obtenue via le procédé selon le deuxième aspect. Selon un quatrième aspect, l'invention concerne un équipement de modélisation d'au moins une partie d'une surface non-axisymétrique d'une portion d'une pièce, la portion étant limitée par un premier et un deuxième plan extrémal, l'équipement étant caractérisé en ce qu'il comprend des moyens de traitement de données configurés pour mettre en oeuvre : - Un module de modélisation d'au moins deux surfaces élémentaires non-axisymétriques composant ladite surface, chacune des surfaces élémentaires s'étendant entre les deux plans extrémaux de la portion et étant définie par au moins une courbe de construction s'étendant selon un plan sensiblement parallèle aux plans extrémaux, une première et une deuxième surface élémentaire étant juxtaposées de sorte à présenter une jonction au moins Cl au niveau des plans extrémaux et au moins une discontinuité ; - Un module de construction d'au moins une courbe de raccord de classe Cl définie par un premier et un deuxième point extrémal respectivement disposés sur la première et la deuxième surface élémentaire, ladite courbe de raccord étant tangente à la première et la deuxième surface respectivement en le premier et le deuxième point extrémal ; - Un module de modification locale des première et deuxième surfaces élémentaires de sorte à respecter ladite courbe de raccord au voisinage de ladite discontinuité ; - Un module de restitution sur une interface dudit équipement de l'au moins une partie de la surface obtenue. Selon un cinquième et un sixième aspect, l'invention concerne respectivement un produit programme d'ordinateur comprenant des instructions de code pour l'exécution d'un procédé selon le premier aspect de l'invention de modélisation d'au moins une partie d'une surface nonaxisymétrique d'une portion d'une pièce ; et un moyen de stockage lisible par un équipement informatique sur lequel un produit programme d'ordinateur comprend des instructions de code pour l'exécution d'un procédé selon le premier aspect de l'invention de modélisation d'au moins une partie d'une surface non-axisymétrique d'une portion d'une pièce. PRESENTATION DES FIGURES D'autres caractéristiques et avantages de la présente invention apparaîtront à la lecture de la description qui va suivre d'un mode de réalisation préférentiel. Cette description sera donnée en référence aux dessins annexés dans lesquels : - la figure la précédemment décrite représente un exemple de réacteur sur des parois duquel le procédé selon l'invention est mis en oeuvre ; - les figures I b-lc illustrent un détail d'une veine non-axisymétrique ; - la figure 2 représente un système pour la mise en oeuvre du procédé selon l'invention ; - les figures 3a-3b représentent deux exemples de courbes de construction d'une surface élémentaire utilisées par le procédé selon l'invention ; - la figure 4a-4d illustrent la construction d'une surface de raccord au niveau d'un bord d'attaque ; - la figure 5 représente une courbe de raccord utilisée par le procédé selon l'invention. DESCRIPTION DETAILLEE Le présent procédé vise à la modélisation d'au moins une partie d'une surface non-axisymétrique d'une portion d'une pièce. Cette pièce peut être un moyeu ou un carter, et de façon générale toute pièce d'un réacteur présentant une paroi interne/externe contre laquelle l'air circule. En référence à la figure la précédemment introduite (qui représente un réacteur CFM-56), on pourra citer les zones suivantes (qui ont été entourées pour faciliter leur localisation) : - Fan - OGV (« Outlet Guide Vane », i.e. aubages redresseurs du flux secondaire) + bras secondaire - Booster + bras primaire - Compresseur Haute Pression - Turbine Haute Pression - Turbine Basse Pression. Dans la suite de la présente description, on prendra l'exemple d'une modélisation du moyeu au niveau de la zone Compresseur Haute Pression (troisième zone entourée en partant de la gauche ou de la droite), mais l'homme du métier saura transposer à n'importe quelle autre pièce du réacteur. On observe sur la figure lb ce moyeu 1 et en particulier une portion 2 non-axisymétrique, limitée en amont par un premier plan extrémal, le « Plan de séparation » PS et en aval par un deuxième plan extrémal, le « Plan de raccord » PR. Il est à noter que le moyeu peut présenter une pluralité de portions 2 (par exemple dans le cas de plusieurs rangs 20 d'aubes), éventuellement séparées par des zones axisymétriques. La portion 2 présente une surface S définissant une forme sensiblement torique, dont au moins une partie est modélisée par le présent procédé. Comme expliqué, ce dernier reprend les méthodes dites à 25 modification par caractérisation de « sections ». Ce procédé propose une solution innovante aux problèmes de continuité : au lieu de chercher à tout prix à éviter les discontinuités entre deux sections, le procédé va au contraire ne rien faire pour empêche leur apparition, puis va les supprimer a posteriori en utilisant des courbes de raccord, qui seront décrites plus loin. 30 La partie de la surface S modélisée comprend ainsi au moins deux surfaces élémentaires (sur la figure 1 a, on compte quatre surfaces élémentaires référencées Sk, k E 41), surface élémentaires qui sont modélisées indépendamment via la construction d'au moins une courbe de construction, appelée également « Plan de construction », lors d'une première étape (a) qui va être décrite plus loin (sur la figure 1 a, trois courbes de construction PC-k définissent la première des surfaces élémentaires Sk). Chacune de ces surfaces élémentaires présente une forme sensiblement rectangulaire et s'étend entre les deux plans extrémaux PS, PR de la portion 2. La surface S est ainsi reconstituée sous la forme d'une pluralité de bandes juxtaposées. Deux surfaces élémentaires voisines Sk et Sk+i présentent une frontière particulière. Comme l'on voit toujours sur les figures, une première et une deuxième surfaces élémentaires voisines sont juxtaposées de sorte à présenter une jonction au moins C1 au niveau des plans extrémaux PS, PR. En d'autres termes, l'intersection entre le moyeu 1 et chacun des plans PR et PS définit une courbe continue et de dérivée continue (cette courbe est fermée et forme une boucle), courbe que respectent les surfaces élémentaires au niveau de leurs extrémités, d'où la continuité de valeur et de dérivée au niveau des points de juxtaposition. Entre leurs deux extrémités, le bord « commun » à deux surfaces élémentaires juxtaposées présente au moins une discontinuité, c'est-à-dire un point tel que les surfaces élémentaires soient disjointes. De façon plus formelle, en définissant un référentiel cylindrique centré sur l'axe du moyeu 1 (tout point de la surface du moyeu est défini par un triplet (r, 61, z), avec r le rayon, c'est-à-dire la distance entre le point de la paroi et l'axe du moyeu, (9 l'angle sur la circonférence du moyeu, et z la position selon l'axe du moyeu), on a: La portion 2 peut être définie par {(r, 61, z), z E [zs, zR]} ; PS et PR peuvent être respectivement définis par {(r , 61, z), r = rs, z = zs} et {(r, 61, z), r = rR, z = zR} ; Une surface élémentaire Sk peut être définie par une équation {(rk(61, z), 6 1, z), z E [zs, zR], E [61k, 61k+1],rk(61, zs) = rs,rk(61, zR) = zR}, telle que 3z E ]zs,zR[, rk(t9 k, z) # rk+i(t 9 k, z), ce qui définit le point de discontinuité. Généralement, cette discontinuité se vérifie sur tout l'intervalle ouvert Ps,zR[ (en d'autres termes, la seule continuité s'observe sur PR et PS). Comme on voit en particulier sur la figure lb, la paroi du moyeu 1 est traversée par des aubes (en d'autres termes des pales). Ces aubes 3, disposées régulièrement ou non (dans ce second cas, on peut procéder par mise à l'échelle) sur la circonférence de la paroi 2, s'étendent radialement. De façon préférée, on obtient une alternance d'aubes et de surfaces élémentaires. En d'autres termes, chaque surface élémentaire s'étend entre deux aubes, et encore en d'autres termes, une surface élémentaire est disposée de chaque côté d'une aube. L'un de ces cotés est l'extrados de l'aube, et l'autre est l'intrados. Ainsi, pour une aube 3 donnée (telle que celle représentée par la figure lb), on peut définir un surface élémentaire d'intrados SI et une surface élémentaire d'extrados SE, qui sont deux surfaces élémentaires juxtaposées. On repère en outre les courbes de constructions PC-I et PC-E de chacune de ces deux surfaces élémentaires SI, SE. De même, une surface élémentaire donné est à la fois la surface élémentaire d'intrados SI d'une aube la bordant, et la surface élémentaire d'extrados SE de l'autre aube la bordant. On comprendra par ailleurs qu'il est possible de modéliser toute surface élémentaire Sk comme l'union de deux « demi-surfaces » chacune s'étendant sur une moitié de la largeur de la veine. La demi-surface « de droite » joue donc le rôle de surface d'extrados d'une aube la bordant, et la demi-surface « de gauche » le rôle de surface d'intrados de l'autre aube voisine. Ce n'est toutefois qu'une convention de représentation parmi d'autres, et ne change rien dans la mise en oeuvre du procédé. Dans la suite de la présente description, on prendra le cas simple d'une aube 3 entourée de deux surfaces élémentaires juxtaposées : la première est la surface élémentaire d'extrados SE de l'aube, et la deuxième est la surface élémentaire d'intrados SI de l'aube. L'aube 3 occupe la majeure part de la jonction, mais laisse une première discontinuité « d'attaque » s'étendant du premier plan extrémal PS au bord d'attaque BA de l'aube 3 et une deuxième discontinuité s'étendant du bord de fuite BF de l'aube 3 au deuxième plan extrémal PR. Le présent procédé permet la construction d'une surface de raccord SR pour chacune des première et deuxième discontinuités (surfaces SR qui seront représentées sur des figures qui seront décrites plus loin). De façon générale, la partie de la surface S est modélisée dans le présent procédé, lors de sa conception, via un équipement informatique 10 du type de celui représenté sur la figure 2. Il comprend des moyens de traitement de données 11 (un ou plusieurs processeur), des moyens de stockage de données 12 (par exemple un ou plusieurs disques dur), des moyens d'interface 13 (composés de moyens de saisie tels qu'un clavier et une souris ou une interface tactile, et de moyens de restitution tels qu'un écran pour affichage des résultats). De façon avantageuse, l'équipement 10 est un supercalculateur, mais on comprendra qu'une mise en oeuvre sur des plateformes variées est tout à fait possible. Modélisation des surfaces élémentaires La première partie (a) du présent procédé voit la modélisation des surfaces élémentaires Sk, SI, SE indépendamment les unes de autres. Il est à noter que dans l'hypothèse avantageuse (mais non limitative) d'une périodicité dans la circonférence du moyeu (c'est-à-dire si les aubes 3 sont identiques et réparties uniformément), toutes les surfaces élémentaires seront identiques et il suffit d'en modéliser une puis de la dupliquer Le principe est de déterminer un optimum géométrique aéromécanique et/ou acoustique, de lois d'une surface élémentaire, en d'autres termes d'une ou plusieurs courbes décrivant la valeur d'une grandeur physique (en particulier le rayon, c'est-à-dire la distance d'un point de la surface élémentaire au centre du moyeu) le long d'une largeur de la surface élémentaire (comme expliqué avant, une surface élémentaire Sk est ainsi définie par ses courbes de construction PC-k, chacune s'étend selon un plan sensiblement parallèle aux plans extrémaux PS, PR de la portion 2) dans un environnement donné, par l'exécution d'un grand nombre de calculs de simulation. De nombreux critères peuvent être choisis comme critères à optimiser lors de la modélisation d'une surface élémentaire. A titre d'exemple, on peut tenter de maximiser des propriétés mécaniques telles que la résistance aux contraintes mécaniques, les réponses fréquentielles, les déplacements des aubes 3, des propriétés aérodynamiques telles que le rendement, l'élévation de pression, la capacité de débit ou la marge au pompage, etc. Pour cela il est nécessaire de paramétriser la loi que l'on cherche à optimiser, c'est-à-dire d'en faire une fonction de N paramètres d'entrée. L'optimisation consiste alors à faire varier (en général aléatoirement) ces différents paramètres sous contrainte, jusqu'à déterminer leurs valeurs optimales pour un critère prédéterminé. Une courbe « lissée » est ensuite obtenue par interpolation à partir des points de passage déterminés. Le nombre de calculs nécessaires est alors directement lié (linéairement voire exponentiellement) au nombre de paramètres d'entrée du problème. De nombreuses méthodes sont connues, mais de façon préférée on mettra en oeuvre une méthode similaire à celle décrite dans la demande de brevet FR1353439, qui permet une excellente qualité de modélisation, sans consommation élevée de puissance de calcul, tout en limitant le phénomène de Runge (« ondulation » excessive de la surface). De façon plus générale, le procédé comprend avantageusement pour chaque courbe de construction PC-k, PC-I, PC-E de chaque surface élémentaire Sk, SI, SE la mise en oeuvre de sous-étapes de : (ai) paramétrisation de la courbe de construction PC en tant que courbe de classe C1 représentant la valeur d'une grandeur physique caractéristique de ladite surface S (typiquement le rayon) en fonction d'une position le long d'au moins le long d'une largeur de la surface élémentaire Sk, SI, SE, la courbe étant définie par : - Deux points extrémaux PCUI, PCUE (pour « Point de Contrôle Utilisateur ») respectivement sur l'intrados et l'extrados de ladite largeur de la surface élémentaire Sk, SI, SE; - Au moins une spline ; la paramétrisation étant mise en oeuvre selon un ou plusieurs paramètres définissant au moins un des points de contrôle extrémaux PCUI, PCUE ; (a2) Détermination de valeurs optimisées desdits paramètres de ladite courbe. La première sous-étape (ai) est mise en oeuvre par les moyens de traitement de données 11 sous contrôle d'un opérateur. La position le long de la largeur de la surface élémentaire Sk, SI, SE est exprimée en fonction de la longueur de corde (en abscisse), et plus précisément la longueur de corde « normalisée », c'est-à-dire exprimée entre 0 et 1 lorsque l'on traverse la surface élémentaire Sk, SI, SE, en d'autres termes lorsqu'on passe d'une aube 3 à une autre. Il est à noter que comme expliqué, la surface élémentaire Sk, SI, SE est située au niveau de l'intrados de l'une de ces deux aubes 3, et de l'extrados de l'autre, c'est pourquoi le chemin correspond à un passage intrados-extrados. On comprendra toutefois que l'invention n'est aucunement limitée à l'expression d'une courbe représentant la valeur d'une grandeur en fonction d'une largeur, et que d'autres repères sont possibles. Cette courbe représentant la valeur d'une grandeur physique doit être comprise comme la modélisation d'une loi de cette grandeur physique (elle sera à ce titre désignée sous le nom de « courbe de modélisation » pour la distinguer des courbes de Bézier en termes de terminologie). Ladite grandeur physique peut être n'importe quelle grandeur ayant une signification aéromécanique et/ou acoustique pour la conception de surface élémentaires Sk, SI, SE, mais il s'agit typiquement du rayon mentionné avant. La courbe de modélisation est de classe de régularité au moins Cl, c'est-à-dire qu'elle correspond à une fonction continue et au moins de lè' dérivée continue sur son espace de définition (la largeur de la surface élémentaire). On verra plus loin l'importance de cette condition. En pratique, la courbe obtenue est Cc° par morceaux (fonctions indéfiniment dérivables sur chaque intervalle), avec continuité de la courbe et de la dérivée au niveau d'éventuels raccords (les points de contrôle intermédiaires qui vont être définis). On comprendra qu'il s'agit de conditions minimales et que la courbe peut tout à fait être par exemple Cc° sur l'ensemble de son espace de définition. La courbe est définie grâce à ses points de contrôle. Comme l'on voit sur la figure 3a, qui représente un premier exemple de modélisation, deux points de contrôle utilisateur extrémaux PCUI et PCUE (qui sont respectivement les points de contrôle « Intrados » et « Extrados », selon le coté de l'aube 3 au niveau duquel ils se trouvent) sont fixés et définissent l'étendue de la surface élémentaire (i.e. le domaine de définition de la courbe). La courbe de modélisation comprend au moins une spline qui relie ces points extrémaux. Les splines sont des courbes polynomiales paramétriques, parmi lesquelles on peut citer préférentiellement les courbes de Bézier définies comme combinaisons de N+1 polynômes élémentaires dits Polynômes de Bernstein : on définit une courbe de Bézier par l'ensemble de points riv_oBiN(t) - Pi ,t E [0,1], les BiN(t) = (Ni)tN(1 - étant les N+1 polynômes de Bernstein de degré N. Les points {Po, P1... PN} sont appelés points de contrôle « implicites » de la courbe et constituent les variables grâce auxquelles une loi d'une 30 surface élémentaire peut être modélisée par une courbe de Bézier. Ces points sont appelés « implicites » car une courbe de Bézier peut être vue comme l'ensemble des barycentres des N+1 points de contrôle pondérés d'un poids égal à la valeur du polynôme de Bernstein associé à chaque point de contrôle. En d'autres termes, ces points agissent comme des poids localisés attirants la courbe généralement sans qu'elle n'y passe (hormis le premier et le dernier, correspondant respectivement à t=0 et t=1, et certains cas d'alignement de points). De façon générale, dans les techniques de modélisation connues d'une loi utilisant une courbe de Bézier, les points de contrôle extrémaux Po et PN de la courbe utilisée sont fixés (ils définissent l'étendue sur laquelle la modélisation va être mise en oeuvre), mais les autres points {P1... présentent des coordonnées mobiles constituant les paramètres d'entrée pour l'algorithme d'optimisation. Dans le présent procédé, un (voire plusieurs) point de contrôle utilisateur intermédiaire PCUm (pour « Milieu ») est avantageusement disposé entre ces deux points extrémaux PCUI et PCUE. Ce cas est représenté par la figure 3b. Le ou les points intermédiaires sont des points de contrôle « explicite » et non implicite, car la courbe y passe. En effet, cette dernière comprend alors au moins deux splines connectées au niveau dudit point intermédiaire. Dans un mode de réalisation, la courbe de modélisation peut être 20 ainsi constituée seulement d'une séquence de courbes de Bézier, chacune s'étendant entre un point de contrôle extrémal PCUI, PCUE et un point intermédiaire PCUm, ou entre deux points intermédiaires PCUm. En d'autres termes, tous les points de contrôle utilisateur PCUI, PCUN, PCUE (extrémaux ou intermédiaires) de la courbe sont des points de 25 contrôle extrémaux Po, PN d'une courbe de Bézier. Le fait que la courbe soit de classe C1 impose que chaque point intermédiaire PCUm assure une continuité y compris sur la dérivée (même tangente). 30 Le traitement paramétrise la courbe de modélisation en fonction de paramètres de points de contrôle extrémaux (et éventuellement intermédiaires) définissant des extrémités de splines constituant la courbe de modélisation d'un plan de construction PC-k, PC-I, PC-E. En particulier, Si ces splines sont des courbes de Bézier, chacune peut être entièrement déterminée par les points PCU définissant ses extrémités. En d'autres termes, les paramètres des points PCU (en termes de coordonnées, et de dérivées) servent de conditions aux limites pour le calcul par les moyens de traitement de données 11 des coordonnées des points de contrôle implicites des différentes courbes de Bézier, qui sont choisies du degré minimal suffisant pour satisfaire ces conditions limites. L'étape (ai) comprend alors la définition des points implicites des courbes de Bézier en fonction des paramètres des points PCU constituant leurs extrémités. Le ou les paramètres définissant un point de contrôle PCU sont ainsi choisis parmi une abscisse du point, une ordonnée du point, une orientation de tangente à la courbe au niveau du point et un (dans le cas d'un point de contrôle extrémal PCUI, PCUE : on ne peut prendre en compte que la demi-tangente dans le domaine de définition de la courbe, celle à gauche pour PCUE et celle à droite pour PCU1) ou deux (dans le cas d'un point de contrôle intermédiaire PCUm) coefficients de tension chacun associé à une demi-tangente à la courbe au niveau du point. Sachant que la courbe est de classe C1 (dérivée continue), l'orientation de la tangente doit être la même de part et d'autre d'un point intermédiaire PCUm. En revanche, la « longueur » des deux demi-tangentes peut être différente de part et d'autre du point, longueur qui traduit la propension de chaque courbe de Bézier de part et d'autre du point à « coller » à la tangente. C'est ce que modélisent les « coefficients de tension » évoqués précédemment. De façon concrète, chaque coefficient de tension associé à une demi-tangente à la courbe au niveau d'un point intermédiaire est fonction d'une valeur de dérivée seconde de la courbe au voisinage du point. En effet, la valeur de la dérivée seconde au voisinage du point de contrôle traduit la « vitesse » avec laquelle la courbe s'écarte de la tangente. Et le fait que la courbe de modélisation ne soit pas forcément de classe C2 autorise des discontinuités de dérivée seconde au niveau des points intermédiaires. Il est à noter que tous les paramètres évoqués précédemment ne sont pas nécessairement effectivement utilisés. En effet, dans la majorité des cas, un ou plusieurs de ces paramètres (une abscisse du point, une ordonnée du point, une orientation de tangente à la courbe au niveau du point et deux coefficients de tension chacun associé à une demi-tangente à la courbe au niveau du point) présentent une valeur prédéterminée fixée par l'utilisateur, et ne sont donc pas utilisés comme « variables » pour la suite du procédé. Selon une deuxième sous-étape (a2), le procédé comprend une 15 étape de détermination par les moyens de traitement de données 11 de valeurs optimisées (et si possible optimales) desdits paramètres de ladite courbe. Il s'agit d'une étape d'optimisation. De nombreuses techniques pour la mise en oeuvre de cette étape sont connues de l'homme du métier, et on pourra par exemple simplement 20 faire varier pseudo-aléatoirement les paramètres choisis variables tout en effectuant des simulations pour déterminer ces valeurs optimisées (c'est-à-dire pour lesquelles le critère choisi est maximalisé) des paramètres des points de contrôle PCU. L'invention n'est toutefois pas limitée à cette possibilité. 25 Une fois les valeurs des paramètres déterminées, les moyens de traitement 11 obtiennent l'équation des courbes de constructions PC-k, PC-1, PC-E, et reconstruisent la totalité de la surface élémentaire Sk, SI, SE (en utilisant les équations des plans extrémaux PS, PR). L'étape est éventuellement répétée pour modéliser diverses surfaces 30 élémentaires Sk, SE, SI (alternativement, il est possible de considérer que toutes les surfaces élémentaires sont identiques et de simplement transposer la modélisation obtenue). Courbe de raccord Une fois les première surface élémentaire SE et deuxième surface élémentaire SI modélisées (étape (a)), on obtient l'état représenté par la figure 4a. On constate la présence de la première discontinuité « d'attaque ». Dans la suite de la présente description, on verra uniquement l'exemple du traitement de cette première discontinuité, mais on comprendra que l'homme du métier pourra mettre en oeuvre les mêmes étapes au niveau de la deuxième discontinuité « de fuite ». Il est à noter que les bords d'attaque BA et de fuite BF, sont ici des bords « réels » (c'est-à-dire les points qui sont réellement les plus en avant et en arrière de l'aube 3), quelques fois différents des bords théoriques (qui suivent le squelette de l'aube). Par exemple, la courbure de l'aube 3 fait que généralement le bord d'attaque BA est décalé vers l'intrados. Les bords BA, BF réels apparaissent naturellement par construction lors de la mise en oeuvre du présent procédé. Par exemple, lorsqu'on observe la figure 4a, on constate pour chacune des surfaces élémentaires SE, SI un point d'intersection avec l'avant de l'aube 3 au niveau de l'extrémité de la discontinuité. Ces deux points définissent le bord d'attaque réel BA, matérialisé sur la figure 4b. Les courbes des surfaces élémentaires SE, SI concernées sont références PC-BA sur cette figure. Le bord BA « réel » correspond en effet à une ligne de reflet, c'est-à-dire une courbe dont la normale à la surface (de l'aube 3) présente, en chaque point, le même angle par rapport à une direction déterminée, en l'espèce celle de la jonction entre les deux surfaces élémentaires SE, SI. Cette prise en compte des bords réels BA, BF permet une correction automatique des erreurs dues à l'extrapolation du squelette de l'aube 3, et aide la mise en oeuvre de l'étape suivante du présent procédé. Dans cette étape (b), est construite au moins une courbe de raccord CR de classe C1 définie par un premier et un deuxième point extrémal PCRI, PCRE (« point de contrôle de raccord ») respectivement disposés sur la première et la deuxième surface élémentaire Sk, SI, SE. Typiquement, chaque point extrémal est sur une ligne médiane de la surface élémentaire SE, SI, les deux points passant par un plan parallèle aux plans PR, PS, de sorte à ce que la courbe de raccord soit tangente à l'aube 3 au niveau du bord de d'attaque BA (respectivement du bord de fuite BF). La courbe de raccord CR est en outre tangente à la première et la deuxième surface SE, SI respectivement en le premier et le deuxième point extrémal PCRI, PCRE. On obtient la courbe CR telle que représentée sur la figure 4b. De façon particulièrement préférée, la courbe de raccord CR passe par un point milieu PCRm disposé sur le bord d'attaque BA (respectivement le bord de fuite BF) de l'aube 3, de sorte que la distance entre le point milieu PCRm et la première surface élémentaire SE est égale à la distance entre le point milieu PCRm et la deuxième surface élémentaire SI. La figure 5 permet de visualiser facilement cela : le point milieu PCRm est en d'autres termes au milieu de la marche existante entre les première et deuxième surfaces élémentaires SE, SI. On remarque 20 également le respect des tangentes. De nombreuses courbes de raccord CR satisfaisant ces conditions sont possible, mais de façon préférée cette dernière est une ou plusieurs splines (voire courbes de Béziers) entièrement définies par les points extrémaux PCRI, PCRE et le point milieu PCRm, de façon similaire à ce qui 25 est utilisé pour la modélisation des plans de constructions PC-I, PC-E (voir ci-avant). Surface de raccord 30 La courbe de raccord CR est utilisée par les moyens de traitement de données 11 pour définir une surface plus étendue, dite surface de raccord SR, qui recouvre la discontinuité. Ainsi, dans une étape (c), le procédé comprend la modification locale des première et deuxième surfaces élémentaires Sk, SI, SE de sorte à respecter ladite courbe de raccord au voisinage de ladite discontinuité. Dans le cas des première et deuxième discontinuités d'attaque et de fuite, cette étape consiste en la construction de la surface de raccord SR présentant en tout point de son bord une continuité de tangence avec la première ou la deuxième surface élémentaire Sk, SI, SE. La surface de raccord SR remplace localement les surfaces élémentaires. Comme représenté sur la figure 4c, cela comprend d'abord la construction d'une première moitié de la surface de raccord SR, qui s'étend d'un bord extrémal PR, PS de la portion 2 jusqu'à la courbe de construction CR. L'équation de chacun des bords de cette demi-surface est connue (dans le cas d'une forme sensiblement rectangulaire telle que visible sur la figure 4c, les quatre cotés sont respectivement définis par l'équation de SI, l'équation de PR, l'équation de SE, et l'équation de CR). En ce qui concerne l'intérieur de cette surface, n'importe quelle géométrie « saine » (continue, dérivable, et à la forme simple) peut faire l'affaire. Arrivé à ce point, il existe toujours une discontinuité entre la surface de raccord SR et les surfaces élémentaires SE, SI au niveau de la courbe 20 de raccord. Il suffit alors d'étendre la surface de raccord SR le long de l'aube 3, en particulier jusqu'à atteindre une courbe de construction PC-I, PC-E de chacune des surfaces élémentaires SE, SI (il suffit que chaque surface élémentaire SI, SE présente au moins une courbe de construction PC-k, 25 PC-I, PC-E disposée entre le bord d'attaque BA et le bord de fuite BF de l'aube 3 (en d'autres termes ailleurs qu'au niveau des discontinuités), de sorte que la surface de raccord SR s'étende jusqu'à la première courbe de construction PC-I, PC-E au-delà de la courbe de raccord CR), comme l'on voit sur la figure 4d. La surface de raccord SR continue à respecter en tout 30 point de son bord une continuité de tangence. La surface de raccord SR prend alors la forme d'un U autour de bord d'attaque BA (ou de fuite BF) de l'aube 3. L'ensemble formé par les surfaces élémentaires SE, SI et les surfaces de raccord SR est alors entièrement Cl et respecte le squelette de l'aube 3. Il définit la partie de la surface S non-axisymétrique modélisée. Dans une dernière étape (d), le modèle obtenu pour la surface S est restitué par les moyens d'interface 13 de l'équipement 10 pour exploitation, par exemple par l'affichage des surfaces en trois dimensions ou de leurs équations. Procédé de fabrication et pièce Une fois sa surface S modélisée, la pièce 1 (qui comme expliquée avant peut être notamment un moyeu ou un carter de réacteur) peut être fabriquée. Est ainsi proposé un procédé de fabrication d'une pièce dont une portion présente une surface non-axisymétrique, le procédé comprenant des étapes de - - Mise en oeuvre du procédé selon le premier aspect de sorte à modéliser au moins une partie de ladite surface nonaxisymétrique de la portion de la pièce ; - Fabrication de ladite pièce conformément à la modélisation de l'au moins une partie de la surface obtenue. La pièce ainsi obtenue (dont une portion présente une surface non- axisymétrique modélisée comme expliqué précédemment), présente une géométrie plus saine et de meilleure qualité que les pièces connues. Elle peut avoir une solidité augmentée, de meilleures propriétés aérodynamiques, des niveaux de bruit réduits, etc. Equipement L'équipement 10 (représenté sur la figure 2) pour la mise en oeuvre du procédé de modélisation d'au moins une partie d'une surface S non- axisymétrique d'une portion 2 d'une pièce 1 comprend des moyens de traitement de données 11 configurés pour mettre en oeuvre : - Un module de modélisation des au moins deux surfaces élémentaires Sk, SI, SE non-axisymétriques composant ladite surface S, chacune des surfaces élémentaires Sk, SI, SE s'étendant entre les deux plans extrémaux PS, PR de la portion 2 et étant définie par au moins une courbe de construction PC-k, PC-I, PC-E s'étendant selon un plan sensiblement parallèle aux plans extrémaux PS, PR, une première et une deuxième surface élémentaire Sk, SI, SE étant juxtaposées de sorte à présenter une jonction au moins C1 au niveau des plans extrémaux PS, PR et au moins une discontinuité ; - Un module de construction d'au moins une courbe de raccord CR de classe C1 définie par un premier et un deuxième point extrémal PCRI, PCRE respectivement disposés sur la première et la deuxième surface élémentaire Sk, SI, SE, ladite courbe de raccord CR étant tangente à la première et la deuxième surface respectivement en le premier et le deuxième point extrémal PCRI, PCRE ; - Un module de modification locale des première et deuxième surfaces élémentaires Sk, SI, SE de sorte à respecter ladite courbe de raccord CR au voisinage de ladite discontinuité ; - Un module de restitution sur une interface 13 dudit équipement 10 de l'au moins une partie de la surface S obtenue. Produit programme d'ordinateur Selon d'autres aspects, l'invention concerne un produit programme d'ordinateur comprenant des instructions de code pour l'exécution (sur des moyens de traitement de donnés 11, en particulier ceux de l'équipement 10) d'un procédé selon le premier aspect de l'invention de modélisation d'au moins une partie d'une surface S non-axisymétrique d'une portion 2 d'une pièce 1, ainsi que des moyens de stockage lisibles par un équipement informatique (par exemple une mémoire 12 de cet équipement 10) sur lequel on trouve ce produit programme d'ordinateur. REVENDICATIONS1. Procédé de modélisation d'au moins une partie d'une surface (S) non-axisymétrique d'une portion (2) d'une pièce (1), la portion (2) étant limitée par un premier et un deuxième plan extrémal (PS, PR), le procédé étant caractérisé en ce qu'il comprend la mise en oeuvre, par des moyens de traitement de données (11) d'un équipement (10), d'étapes de : (a) Modélisation d'au moins deux surfaces élémentaires (Sk, SI, SE) non-axisymétriques composant ladite surface (S), chacune des surfaces élémentaires (Sk, SI, SE) s'étendant entre les deux plans extrémaux (PS, PR) de la portion (2) et étant définie par au moins une courbe de construction (PC-k, PC-I, PC-E) s'étendant selon un plan sensiblement parallèle aux plans extrémaux (PS, PR), une première et une deuxième surface élémentaire (Sk, SI, SE) étant juxtaposées de sorte à présenter une jonction au moins C1 au niveau des plans extrémaux (PS, PR) et au moins une discontinuité ; (b) Construction d'au moins une courbe de raccord (CR) de classe C1 définie par un premier et un deuxième point extrémal (PCRI, PCRE) respectivement disposés sur la première et la deuxième surface élémentaire (Sk, SI, SE), ladite courbe de raccord (CR) étant tangente à la première et la deuxième surface respectivement en le premier et le deuxième point extrémal (PCRI, PCRE) ; (c) Modification locale des première et deuxième surfaces élémentaires (Sk, SI, SE) de sorte à respecter ladite courbe de raccord (CR) au voisinage de ladite discontinuité ; (d) Restitution sur une interface (13) dudit équipement (10) de l'au moins une partie de la surface (S) obtenue. 2. Procédé selon la revendication 1, dans lequel une pluralité d'aubes (3) est disposée sur la circonférence de la portion (2) de la pièce (1), chaque surface élémentaire (SI, SE) s'étendant entre deux aubes (3), de sorte que les première et deuxième surfaces élémentaires (SI, SE) sont disposées de part et d'autre d'une aube (3), respectivement sur l'intrados et l'extrados. 3. Procédé selon la revendication 2, dans lequel la jonction entre les première et deuxième surfaces élémentaires (SI, SE) présente une première discontinuité s'étendant du premier plan extrémal (PS) au bord d'attaque (BA) de l'aube (3) et une deuxième discontinuité s'étendant du bord de fuite (BF) de l'aube (3) au deuxième plan extrémal (PR). 4. Procédé selon la revendication 3, dans lequel une courbe de raccord (CR) est construite à l'étape (b) pour chacune des première et deuxième discontinuités, chaque courbe de raccord (CR) étant tangente à l'aube (3) au niveau du bord de fuite (BF) ou du bord d'attaque (BA). 5. Procédé selon la revendication 4, dans lequel une courbe de raccord (CR) passe par un point milieu (PCRm) disposé sur le bord de fuite (BF) ou du bord d'attaque (BA) de l'aube (3) de sorte que la distance entre le point milieu (PCRm) et la première surface élémentaire (SE) est égale à la distance entre le point milieu (PCRm) et la deuxième surface élémentaire (SI). 6. Procédé selon la revendication 5, dans lequel une courbe de raccord (CR) est une spline entièrement définie par les points extrémaux (PCRI, PCRE) et le point milieu (PCRm). 7. Procédé selon l'une des revendications 3 à 6, dans lequel l'étape (c) comprend pour chacune des première et deuxième discontinuité la construction d'une surface de raccord (SR) présentant en tout point de son bord une continuité de tangence avec la première ou la deuxième surface élémentaire (Sk, SI, SE). 8. Procédé selon la revendication 7, dans lequel chaque surface élémentaire (Sk, SI, SE) présente au moins une courbe de construction (PC-k, PC-I, PC-E) disposée entre le bord d'attaque (BA) et le bord de fuite (BF) de l'aube (3), chacune des deux surfaces de raccord (SR) s'étendant sur chacune des première et deuxième surfaces élémentaires (Sk, SI, SE) d'un bord extrémal (PR, PS) jusqu'à la première courbe de construction (PC-k, PC-I, PC-E) au-delà de la courbe de raccord (CR). 9. Procédé selon l'une des revendications 2 à 8, dans lequel l'étape (a) comprend pour chaque courbe de construction (PC-k, PC-1, PC-E) de chaque surface élémentaire (Sk, SI, SE) la mise en oeuvre de sous-étapes de : (ai) Paramétrisation de la courbe de construction (PC-k, PC-I, PC-E) en tant que courbe de classe C1 représentant la valeur d'une grandeur physique caractéristique de ladite surface (S) en fonction d'une position le long d'une largeur de la surface élémentaire (Sk, SI, SE), la courbe étant définie par : - Deux points de contrôle extrémaux (PC1.11, PCUE), respectivement sur chacune des deux aubes (3) entre lesquelles ladite surface élémentaire (Sk, SI, SE) s'étend, définissant l'étendue de ladite largeur de la surface élémentaire (Sk, SI, SE) ; - Au moins une spline ; la paramétrisation étant mise en oeuvre selon un ou plusieurs paramètres définissant au moins un des points de contrôle extrémaux (PC1.11, PCUE) ;(a2) Détermination de valeurs optimisées desdits paramètres de ladite courbe. 10. Procédé selon la revendication 9, dans lequel chaque courbe de construction (PC-k, PC-I, PC-E) de chaque surface élémentaire (Sk, SI, SE) est également définie par au moins un point de contrôle intermédiaire (PCUm) disposé entre les points extrémaux (PCUI, PCUE), et au moins deux splines connectées au niveau dudit point de contrôle intermédiaire (PCUm) ; la paramétrisation étant également mise en oeuvre selon un ou plusieurs paramètres définissant ledit point de contrôle intermédiaire (PCUm). 11. Procédé selon l'une des revendications 9 et 10, dans lequel le ou les paramètres définissant un point de contrôle intermédiaire ou extrémal (PCUI, PCUm, PCUE) définissant une courbe de construction (PC- k, PC-I, PC-E) sont choisis parmi une abscisse du point, une ordonnée du point, une orientation de tangente à la courbe au niveau du point et au moins un coefficient de tension associé à une demi-tangente à la courbe au niveau du point. 12. Procédé selon la revendication 11, dans lequel un coefficient de tension associé à une demi-tangente à la courbe au niveau d'un point intermédiaire ou extrémal est fonction d'une valeur de dérivée seconde de la courbe au voisinage du point. 13. Procédé selon l'une des revendications 9 à 12, dans laquelle chaque spline d'une courbe de construction (PC-k, PC-I, PC-E) est une courbe de Béziers entièrement déterminée par les points de contrôle (PCUI, PCUm, PCUE) définissant ses extrémités.30 14. Procédé de fabrication d'une pièce (1) dont une portion (2) présente une surface (S) non-axisymétrique, le procédé comprenant des étapes de : - Mise en oeuvre du procédé selon l'une des revendications 1 à 13 de sorte à modéliser au moins une partie de ladite surface (S) non-axisymétrique de la portion (2) de la pièce (1) ; - Fabrication de ladite pièce (1) conformément à la modélisation de l'au moins une partie de la surface (S) obtenue. 15. Pièce (1) dont une portion (2) présente une surface (S) non-axisymétrique obtenue via le procédé selon la revendication 14. 16. Equipement (10) de modélisation d'au moins une partie d'une surface (S) non-axisymétrique d'une portion (2) d'une pièce(1), la portion (2) étant limitée par un premier et un deuxième plan extrémal (PS, PR), l'équipement (10) étant caractérisé en ce qu'il comprend des moyens de traitement de données (11) configurés pour mettre en oeuvre : - Un module de modélisation d'au moins deux surfaces élémentaires (Sk, SI, SE) non-axisymétriques composant ladite surface (S), chacune des surfaces élémentaires (Sk, SI, SE) s'étendant entre les deux plans extrémaux (PS, PR) de la portion (2) et étant définie par au moins une courbe de construction (PC-k, PC-I, PC-E) s'étendant selon un plan sensiblement parallèle aux plans extrémaux (PS, PR), une première et une deuxième surface élémentaire (Sk, SI, SE) étant juxtaposées de sorte à présenter une jonction au moins C1 au niveau des plans extrémaux (PS, PR) et au moins une discontinuité ; - Un module de construction d'au moins une courbe de raccord (CR) de classe C1 définie par un premier et un deuxième point extrémal (PCRI, PCRE) respectivement disposés sur la première et la deuxième surface élémentaire(Sk, SI, SE), ladite courbe de raccord (CR) étant tangente à la première et la deuxième surface respectivement en le premier et le deuxième point extrémal (PCRI, PCRE) ; - Un module de modification locale des première et deuxième surfaces élémentaires (Sk, SI, SE) de sorte à respecter ladite courbe de raccord (CR) au voisinage de ladite discontinuité ; - Un module de restitution sur une interface (13) dudit équipement (10) de l'au moins une partie de la surface (S) obtenue. 17. Produit programme d'ordinateur comprenant des instructions de code pour l'exécution d'un procédé selon l'une des revendications 1 à 13 de modélisation d'au moins une partie d'une surface (S) non-axisymétrique d'une portion (2) d'une pièce (1). 18. Moyen de stockage lisible par un équipement informatique sur lequel un produit programme d'ordinateur comprend des instructions de code pour l'exécution d'un procédé selon l'une des revendications 1 à 13 de modélisation d'au moins une partie d'une surface (S) non-axisymétrique d'une portion (2) d'une pièce (1).