FR 3106675 A1 20210730 FR 2000708 A 20200124 L’invention concerne un procédé de détermination de paramètres de grenaillage par apprentissage pour une pièce, comprenant : a) Fournir des spécifications de grenaillage de la pièce ; b) Fournir un modèle tridimensionnel de la pièce c) Pour une pluralité de zones : i. déterminer une ou plusieurs combinaisons de paramètres de grenaillage d’une grenailleuse répondant aux spécifications de grenaillage pour ladite zone considérée ; ii. simuler la projection d’un jet de grenaille en chaque point de ladite zone et en déduire en chaque point un état de surface comprenant une information d’intensité Almen et de taux de recouvrement ; d) combiner les états de surface des zones de la pièce afin d’obtenir une cartographie d’intensité Almen et de taux de recouvrement pour le modèle tridimensionnel ; e) déterminer un modèle tridimensionnel post-grenaillage de la pièce : f) comparer dimensionnellement le modèle tridimensionnel post-grenaillage au modèle tridimensionnel de la pièce. Figure à publier avec l’abrégé : Figure n°1 Procédé de détermination de paramètres de grenaillage par apprentissage Domaine technique de l’invention L’invention porte sur les opérations de grenaillage de pièces destinées au domaine de l’aéronautique, et en particulier du paramétrage intelligent des machines de grenaillage. Etat de la technique antérieure Le grenaillage de précontrainte est largement utilisé dans le domaine de l’aéronautique. Il s’agit d’un traitement mécanique de surface également considéré comme un procédé de renforcement mécanique. Cette opération consiste à déformer plastiquement une couche superficielle de la surface d’une pièce en projetant des billes à grande vitesse. Lors de l’impact, chaque bille forme un renfoncement local à la surface de la pièce. Les opérations de grenaillage ont les effets techniques suivants : -améliorer la durée de vie en fatigue des pièces industrielles traitées, - compenser des abattements de durée de vie induits par les traitements de surface des pièces, - homogénéiser l’état de surface usinées des pièces. Afin que l’opération de grenaillage induise des propriétés homogènes en surface de la pièce traitée, la surface doit être entièrement recouverte par les impacts de billes. Les meilleurs résultats sont obtenus avec une grenaille homogène (par exemple en acier dur trempé et revenu dont la dureté doit être supérieure à celle des pièces traitées) projetée le plus perpendiculairement possible à la surface à traiter. Il est connu d’utiliser également des billes en verre ou céramique, ces billes étant également appelées grenaille. Dans le domaine de l’aéronautique, le grenaillage de précontraintes s’applique couramment sur des pièces industrielles métalliques complexes telles que des aubes, des disques turbine alvéolés, des disques aubagés monoblocs, des arbres de turbine, des composants de trains d’atterrissage, des roues, pignons. L’intensité de grenaillage est quantifiée en Almen selon la méthode Almen [SAE J442, SAE J443]. La détermination des paramètres d’un procédé de grenaillage pour une pièce industrielle complexe, c’est-à-dire des paramètres entrés en machine, comporte plusieurs étapes. La première étape consiste à analyser la spécification au plan propre à la pièce à traiter par le procédé de grenaillage. Ce document, réalisé par des ingénieurs, comprend les informations relatives au grenaillage de la pièce, par exemple : - les zones de la pièce doivent être grenaillées (par conséquent quelles zones ne doivent pas être grenaillées), - le type de grenaille devant être utilisé, - l’intensité Almen à utiliser pour chacune des zones, et - le taux de recouvrement pour chacune des zones à grenailler. La deuxième étape consiste à choisir de premiers paramètres machines. Cette étape est généralement effectuée par des spécialistes, capable de déterminer, par expérience, les paramètres machines pouvant être utilisés pour respecter les spécifications au plan de la pièce. Une troisième étape consiste à mesurer, suite au grenaillage de la pièce, les grandeurs suivantes : - recouvrement et intensité du grenaillage pour chacune des zones ; - distorsions de la pièce induite. Ce procédé, bien que ne représentant pas de difficultés particulières, s’appuie uniquement sur l’expertise de spécialistes de grenaillage, de leurs essais sur des pièces et sur des mannequins, ce qui pose un problème évident de reproductibilité puisque l’opération dépend d’un opérateur. Par ailleurs, ceci peut se révéler problématique, notamment lorsque ces spécialistes quittent les entreprises. Egalement, certaines exigences, en particulier celles liées aux contraintes résiduelles, au taux de recouvrement, à la rugosité surfacique, ou encore aux tolérances dimensionnelles de la pièce, ne peuvent être contrôlées avant que la pièce ne soit soumise à un procédé de grenaillage. Le contrôle de ces exigences ne peut ainsi être fait que a posteriori , nécessitant donc plusieurs itérations pour trouver les paramètres machines optimaux au regard des exigences des spécifications au plan de la pièce. En outre, la complexité géométrique de certaines zones à grenailler (par exemple des alvéoles dans un disque, perçages, diamètre interne d’arbres de transmission de puissance…) peut rendre le contrôle du grenaillage (intensité Almen, recouvrement, rugosité, contraintes résiduelles) difficile, voire impossible, dans ces zones. Le présent document vise, dans un premier temps, à limiter les tests de grenaillage, permettant de déduire les bons paramètres de la grenailleuse permettant un grenaillage en accord avec spécifications au plan de la pièce. Un autre objectif de l’invention est de proposer une solution permettant de capitaliser le savoir faires relatif au choix de paramètres machines pour le grenaillage. Plus globalement, le présent document a pour but d’apporter une solution simple, efficace et économique aux inconvénients de la technique actuelle exposée ci-dessus. Présentation de l’invention A cette fin, il est proposé un procédé de détermination de paramètres de grenaillage par apprentissage pour une pièce, par exemple de trains d’atterissage, de transmission de puissance, le procédé comprenant les étapes suivantes : a) Fournir des spécifications de grenaillage de la pièce ; b) Fournir un modèle tridimensionnel de la pièce, réalisé par éléments finis à l’aide d’un maillage tridimensionnel de la pièce et subdiviser le modèle tridimensionnel en une pluralité de zones ; c) Pour chacune des zones : i. déterminer une ou plusieurs combinaisons de paramètres de grenaillage d’une grenailleuse répondant aux spécifications de grenaillage pour ladite zone considérée ; ii. simuler la projection d’un jet de grenaille en chaque point de ladite zone et en déduire en chaque point un état de surface comprenant une information d’intensité Almen et de taux de recouvrement ; d) combiner les états de surface des zones de la pièce afin d’obtenir une cartographie d’intensité Almen et de taux de recouvrement pour le modèle tridimensionnel ; e) déterminer un modèle tridimensionnel post-grenaillage de la pièce : i. Calculer un tenseur de contraintes à partir de la cartographie réalisée à l’étape d) ; ii. Appliquer le tenseur de contraintes au modèle tridimensionnel de la pièce pour obtenir le modèle tridimensionnel post-grenaillage ; f) comparer dimensionnellement le modèle tridimensionnel post-grenaillage au modèle tridimensionnel de la pièce. Ce procédé, pouvant être destiné à être mis en œuvre sur un ordinateur, propose une méthode d’apprentissage permettant par simulation de faciliter la détermination des gammes de paramétrage des machines de grenaillage, et de choisir les paramètres optimaux. Contrairement à la technique antérieure, il est possible par simulation de connaitre, les effets induits par les opérations de grenaillage sur la pièce traitée. Ainsi, il est possible de déterminer si les paramètres choisis sont conformes aux exigences explicités dans la spécification au plan de la pièce, notamment les exigences relatives à l’intensité Almen à laquelle a été soumise la pièce, aux distorsions, aux contraintes résiduelles, au taux de recouvrement, à la rugosité surfacique, ou encore aux tolérances dimensionnelles de la pièce. Ce procédé permet donc de valider par calculs et simulation que les exigences de la spécification au plan de la pièce sont respectées en tout point de la pièce, et ceci avant de faire un essai réel sur la pièce. Cela permet une économie de temps non négligeable, dès lors que les essais réels sur la pièce sont réalisés avec des paramètres de grenaillage de la grenailleuse validés par simulation. Ainsi, les essais réels sur la pièce consistent en une étape de validation des paramètres déterminés par le procédé. Il peut également être prévu, pour chaque pièce produite en série, de conserver les paramètres de grenaillage, permettant ainsi une capitalisation et une pérennisation du savoir-faire en termes d’opération de grenaillage. Ce procédé permet en outre de réduire l’utilisation de la machine de grenaillage et le nombre d’essais sur pièces pour déterminer les paramètres de grenaillage requis pour la mise en production d’une série de pièces. Bien entendu, ce procédé est applicable au grenaillage de tout type de pièces en alliage métallique. Par ailleurs, les spécifications de grenaillage utilisé dans le procédé peuvent comprendre un ou plusieurs des éléments suivants : i. des parties de la pièce à grenailler ; ii. des parties de la pièce interdites au grenaillage ; iii. des parties de la pièce pour lesquelles le grenaillage est admissible ; iv. le type de grenaille ; v. des caractéristiques mécaniques et dimensionnelles de la grenaille ; vi. une plage d’intensité de grenaillage pour chacune des parties à grenailler ; vii. une plage de taux de recouvrement minimum et/ maximum pour chacune des parties à grenailler ; viii. profil de contraintes résiduelles minimal ; ix. rugosité surfacique maximum admissible ; x. des tolérances dimensionnelles de la pièce. Ces spécifications sont propres à chaque type de pièces. En effet, en fonction de leur condition d’utilisation, les opérations de grenaillage des pièces diffèrent. Les exigences devant être remplies par les pièces obtenues suite au grenaillage sont des données d’entrées du procédé, de sorte que ces exigences sont prises en considération lors de la sélection de paramètres machines, puis suite à la simulation de ces paramètres machine sur un modèle tridimensionnel de la pièce. L’étape c) i) peut consister à rechercher les paramètres de grenaillage de la grenailleuse, dans une première base de données préalablement établie liant les spécifications de grenaillage, en particulier l’intensité Almen et le taux de recouvrement auxdits paramètres de grenaillage de la grenailleuse. Ainsi, le choix des paramètres de grenaillage de la grenailleuse sont choisis en se basant en particulier sur les deux paramètres des spécifications de grenaillage propre à la pièce à grenailler, à savoir l’intensité Almen et le taux de recouvrement. Le stockage de ces données dans une base de données, permet en particulier une accessibilité accélérée à ces données. Par ailleurs, l’étape c) ii) peut être suivie d’une étape de vérification consistant à vérifier si l’intensité Almen simulée et le taux de recouvrement sont conformes aux plages d’intensité de grenaillage et de taux de recouvrement compris dans les spécifications de grenaillage. Cette étape permet d’obtenir, pour chacun des points de surface grenaillés, une information concernant l’intensité Almen à laquelle ils ont été soumis, ainsi que le taux de recouvrement. Cela permet ainsi, d’obtenir et de superposer sur le modèle 3D de la pièce à grenailler, une cartographie surfacique de la variation de l’intensité Almen ainsi du taux de recouvrement. Cela permet, suite à la simulation du jet de grenaillage sur la pièce à grenailler, avec les paramètres de grenaillage sélectionnés dans la base de donnée, de confirmer ou non les paramètres machines par rapport aux spécifications de grenaillage, notamment vérifier que les spécifications concernant l’intensité Almen ainsi que le taux de recouvrement sont bien respectés. En outre, l’étape d) peut être suivie d’une étape de vérification consistant à vérifier si l’intensité Almen simulée et le taux de recouvrement sont conformes aux plages d’intensité de grenaillage et de taux de recouvrement compris dans les spécifications de grenaillage. La vérification de la conformité de l’intensité Almen et du taux de recouvrement peut ainsi se faire en tout point de la pièce, et ceci sans avoir à grenailler une pièce en condition réelle. Une telle simulation permet de confirmer si des paramètres de grenaillage respectent les exigences des spécifications de grenaillage de la pièce. Pour faciliter le traitement de l’information, il est possible de visualiser l’hétérogénéité de l’intensité Almen auxquelles les zones de la surface sont soumises, ainsi que le taux de recouvrement de ces différentes zones. Pour réaliser l’étape e) : i. on peut calculer une distance z de chaque point d’intégration du maillage à un élément de surface du maillage le plus proche ; ii. on peut calculer le tenseur de contraintes résiduelles initial à appliquer au modèle tridimensionnel aux fins d’obtention du modèle tridimensionnel post-grenaillage, le calcul dudit tenseur étant réalisé pour tout point d’intégration du maillage en fonction de la distance z à l’élément de surface associé et du profil de contraintes résiduelles de type massif semi-infini correspondant aux conditions de grenaillage obtenues (recouvrement, intensité Almen…). Le grenaillage impliquant des distorsions et des déformations de la pièce, l’étape e) consiste en des calculs afin de déterminer pour tout point du volume de la pièce les contraintes résiduelles résultantes de l’opération de grenaillage simulée. Ces calculs participent à l’obtention d’une modélisation de la pièce déformée suite au grenaillage. L’étape ii. peut utiliser une deuxième base de données liant des paramètres de la grenailleuse à des profils de contraintes résiduelles. Cette deuxième base de données sert également de capitalisation du savoir-faire de l’entreprise, notamment des professionnels de grenaillage. Ainsi, la détermination du modèle tridimensionnel post-grenaillage peut comprendre les étapes consistant à : - Calculer pour chaque point du maillage tridimensionnel par éléments finis de l’étape b), la distance z à la surface la plus proche ; - Calculer pour chaque point du maillage bidimensionnel le profil de contraintes résiduelles selon z à partir de la combinaison de paramètres de la grenailleuse ; - Calculer le tenseur des déformations plastiques correspondant ε p et de contraintes résiduelles initial pour chaque point de maillage, en fonction de la combinaison de paramètres de la grenailleuse; - Calculer l’équilibrage statique de la pièce par éléments finis à partir l’objet modélisé tridimensionnel de la pièce et de propriétés élastiques d’un matériau constitutif de la pièce ; - Calculer, pour chaque point du maillage tridimensionnel, la déformation géométrique induite du modèle tridimensionnel, après un équilibrage statique de la pièce, à partir du tenseur de contraintes résiduelles et de l’équilibrage statique. L’identification des profils de contraintes peut s’appuyer sur des abaques, dans lesquelles le profil de contraintes résiduelles est lié aux conditions de grenaillage. Le calcul du tenseur de contraintes résultant de l’opération de grenaillage simulé est effectué pour chaque point du maillage dans le volume de la pièce. La modélisation de la pièce post-grenaillage prend également en considération l’équilibrage statique de la pièce suite au grenaillage, et donc permet d’obtenir un modèle tridimensionnel, où les distorsions induites par le grenaillage sont ainsi prises en compte la modélisation de la pièce suite au grenaillage. Cette étape permet alors de visualiser une prédiction des déformations induites par le grenaillage avec les paramètres de grenaillage sélectionné et utilisés pour la simulation. Il est alors possible, sans réaliser de grenaillage réel, de visualiser dans quelles mesures les paramètres de grenaillage induisent des distorsions, et si celles-ci sont acceptables par rapport aux tolérances acceptables pour la pièce grenaillée. Le présent document concerne en outre un programme d’ordinateur comprenant des instructions pour mettre en œuvre le procédé tel que décrit précédemment. Un tel programme permet ainsi de simuler le grenaillage d’une pièce avec des paramètres de grenaillage choisis en adéquation avec les exigences de grenaillage propre à la pièce à grenailler, de visualiser les distorsions et les déformations géométriques de la pièce en prenant en considération le matériau constitutif de la pièce à grenailler, et de valider ou non les paramètres de grenaillage si les déformations géométriques sont comprises dans les tolérances acceptées. Ceci permet donc de simuler les paramètres de grenaillage et déterminer, avant de grenailler en condition réelles une pièce, si ces paramètres de grenaillage sont conformes aux spécifications de grenaillage. Ainsi, contrairement à l’art antérieur, le grenaillage en conditions réelle sur une pièce n’est effectué qu’avec des paramètres pouvant être, selon la simulation, en adéquation avec les spécifications de grenaillage. Cela réduit donc le nombre d’essais sur une pièce réelle, et donc réduit de fait le nombre d’itérations nécessaires pour déterminer les paramètres de grenaillages adéquats pour une pièce donnée par une personne sans expérience dans le domaine du grenaillage. Par conséquent, il est aussi possible par le biais de ce programme, de mutualiser les connaissances relatives au grenaillage de pièces, de sorte que les personnes en charge du grenaillage de pièces peuvent s’appuyer sur l’expérience d’autres, dont les simulations de paramètres peuvent être conservées en mémoire. Le présent document porte également sur un support de stockage de données sur lequel est enregistré le programme d’ordinateur tel que mentionné ci-dessus. En outre, la présente demande concerne un système de CAO comprenant un processeur couplé à une mémoire et une interface utilisateur graphique et apte à communiquer avec une grenailleuse, tel que le programme décrit ci-dessus est enregistré en mémoire. Brève description des figures représente un logigramme du procédé selon le présent document. représente un exemple de spécifications de grenaillage pour un exemple de pièce à grenailler. représente un logigramme détaillé du procédé selon le présent document. dont la partie A représente une modélisation tridimensionnelle de la pièce à grenailler et la partie B représente une modélisation tridimensionnelle d’une portion de la pièce à grenailler. dont la partie A est une représentation tridimensionnelle d’une simulation de grenaillage de la pièce à grenailler dont la partie A représente une cartographie de l’intensité Almen observé par la surface de la portion de la pièce pour laquelle le grenaillage est simulé pour les paramètres des grenaillages de la machine et dont la partie B représente une cartographie du taux de recouvrement observé par la surface de la portion de la pièce pour laquelle le grenaillage est simulé pour les paramètres des grenaillages de la machine. comprend deux parties, une partie A illustrant le champ de contraintes résiduelles après grenaillage (équilibrés dans la pièce) et une partie B illustrant le champ de déformation de la pièce après grenaillage (comparaison entre la géométrie initiale et déformée (facteur 20 sur les déformations), cette dernière partie B illustre le déplacement avant et après grenaillage représente un schéma des étapes pour calculer des contraintes résiduelles induites par les opérations de grenaillage. Description détaillée de l’invention La figure 1 illustre le principe général du procédé selon l’invention, permettant par apprentissage, de déterminer les paramètres de grenaillage de la grenailleuse. Les paramètres de grenaillage de la grenailleuse sont des paramètres dits machines, et par exemple : - distance de la buse par rapport à la pièce (d buse ), - angle du jet de grenaille, - distance parcourue par la buse (d travail ), - pression du jet, en particulier de la pression l’air comprimé qui permet de projeter le jet à travers la buse de grenaillage ; - débit de la grenaille en sortie de buse, - vitesse de déplacement de la buse, - nombre de cycles, c’est-à-dire du nombre d’aller-retour de la buse par rapport à la zone traitée. La pièce à grenailler peut être, par exemple, une pièce de turbomachine. Cette pièce peut présenter des symétries, dans ce cas, ces symétries de géométrie seront exploitées dans le déroulé du procédé. La première étape A comprend la fourniture des spécifications de grenaillage de la pièce. Les spécifications de grenaillage rassemblent l’ensemble des informations nécessaire au grenaillage d’une pièce afin qu’elle soit utilisable et conforme aux tolérances. En particulier, les spécifications assurent le grenaillage de la pièce n’induise pas de problème de comportement mécanique de cette dernière. Un exemple de spécifications de grenaillage 2 est illustré à la figure 2. Comme on peut le voir, les spécifications comprennent des informations telles que : i. des parties 4, 6 de la pièce à grenailler. ii. des parties 10 de la pièce interdites au grenaillage ; iii. le type de grenaille ; iv. des caractéristiques mécaniques et dimensionnelles de la grenaille ; v. une plage d’intensité de grenaillage pour chacune des parties à grenailler ; vi. une plage de taux de recouvrement minimum et/ maximum pour chacune des parties à grenailler ; vii. profil de contraintes résiduelles minimal ; viii. rugosité surfacique maximum admissible ; ix. des tolérances dimensionnelles de la pièce. Il s’agit d’un disque pour turbine, constitutif d’un labyrinthe d’étanchéité. Dans les parties de la pièce à grenailler, différentes parties sont à distinguer dans l’exemple de la figure 2. Un premier 4 et deuxième 6 type de partie désigne des zones à grenailler d’intensité Almen différentes. Un troisième type 8 de partie désigne les zones, à ne pas grenailler, mais pour lesquels des rebonds de grenaille sont autorisés. Les parties 10 de la pièce interdite au grenaillage, sont des parties où même des rebonds de grenailles ne sont pas tolérés. Lors du grenaillage de la pièce, des protections dites serrées autour de ces parties, sont utilisées. Le type de grenaille peut être des billes en acier, en céramique ou en verre, de diamètres plus ou moins variable. Le choix de la grenaille combiné à l’intensité de grenaillage impacte alors la rugosité surfacique obtenue suite au grenaillage. Le taux de recouvrement est également précisé : il s’agit du ratio de la zone grenaillée sur la zone à grenailler. Pour s’assurer d’un grenaillage complet de la pièce, le taux de recouvrement minimal est généralement choisi dans la plage de valeurs comprises entre [1,25 ; 1,5]. Pour une intensité Almen donnée, plusieurs combinaisons de paramètres de la machine de grenaillage peuvent aboutir à des résultats similaires. L’objet de ce procédé est de fournir une ou plusieurs combinaisons de paramètres de grenaillage permettant d’obtenir un grenaillage conforme aux spécifications, tout en permettant de réduire le temps de grenaillage. La deuxième étape B du procédé consiste à fournir un modèle tridimensionnel de la pièce, réalisé par éléments finis à l’aide d’un maillage tridimensionnel de la pièce et subdiviser le modèle tridimensionnel en une pluralité de zones. Cela peut également être réalisé en utilisant un modèle bidimensionnel de la pièce. Le modèle obtenu, sera alors utilisé pour effectuer la simulation de grenaillage pour des combinaisons données de paramètres. Cette deuxième étape consiste en particulier à fournir un modèle modélisé selon la méthode des éléments finis. L’utilisation de la méthode des éléments finis est nécessaire pour permettre de quantifier les déformations de la pièce suite à la simulation de grenaillage. Pour cela, un modèle tridimensionnel de la pièce 12 (pouvant être obtenu via un outil de conception assisté par ordinateur) est importé dans un logiciel exploitant la méthode des éléments finis, tel que par exemple, le logiciel ABAQUS ®. Le modèle tridimensionnel 12 du disque ici considéré comme la pièce à grenailler, est illustré à la figure 4A. Comme on peut le voir, le disque 12 est une pièce annulaire, présentant une symétrie centrale. Selon que la géométrie de la pièce s’y prête ou non, il est possible de ne modéliser qu’une portion de la pièce 12 pour simplifier la modélisation. Dans l’exemple, la pièce étant un disque, la symétrie de révolution de la pièce permet de simplifier le problème en considérant uniquement une portion 14 de celui-ci. Un modèle tridimensionnel 14 de la portion considérée est alors utilisé pour effectuer un maillage, tel qu’illustré à la figure 4B. Cette méthode s’applique alors toute pièce axisymétrique. Une fois le modèle tridimensionnel de cette portion 14 obtenu, le modèle tridimensionnel est enrichi de données définissant les propriétés élastiques du matériau dans lequel le disque est conçu ainsi que les conditions aux limites de la portion. Les données en lien avec l’élasticité du matériau peuvent être par exemple le module de Young et le coefficient de poisson du matériau utilisé. Quant aux conditions aux limites, la symétrie, le bridage et les efforts auxquels sont soumis la pièce sont pris en considération. Un maillage de la portion 14 de pièce est par la suite calculé sur la surface de la portion de la pièce, de manière à définir au moins 5 éléments, de préférence 10. Un ensemble, représentant le maillage, comprenant les nœuds ainsi que les éléments surfaciques constitutifs du maillage de la portion 14, est enregistré. En parallèle, ou par la suite, un fichier est généré par le logiciel exploitant la méthode des éléments finis, comprenant des calculs d’équilibrage. Dans le cas particulier où une portion 14 de la pièce 12 présentant des symétries a été modélisée, l’étape B comprend en outre une étape visant à établir un modèle bidimensionnel axisymétrique de la pièce 12. En particulier, dans cette étrape, un modèle bidimensionnel de la pièce est fourni en entrée du logiciel de modélisation selon la méthode des éléments finis. Comme pour la première partie de la deuxième étape B, le modèle bidimensionnel est alors utilisé pour réaliser un maillage de la pièce, divisant la surface de la pièce en sous-surfaces. L’ensemble constituant le maillage, c’est-à-dire les nœuds et les éléments surfaciques sont enregistrés. En outre, le modèle bidimensionnel est enrichi de données relatives aux propriétés élastiques des matériaux constitutifs de la pièce. Les données en lien avec l’élasticité du matériau peuvent être par exemple le module de Young et le coefficient de poisson du matériau utilisé. Quant aux conditions aux limites, la symétrie, le bridage et les efforts auxquels sont soumis la pièce sont pris en considération. En parallèle, ou par la suite, un fichier est généré par le logiciel exploitant la méthode des éléments finis, comprenant des calculs d’équilibrage. Egalement, à la fin de cette deuxième étape B, la surface de la pièce 12 est divisée en zones, appelées par la suite zones élémentaires. En particulier, les zones à grenailler sont considérées individuellement, et dans le cas où les zones à grenailler présentent une grande surface, elles sont divisées en zones élémentaires pour lesquelles des instructions pour la grenailleuse sont déterminées individuellement. Chacune de ces zones élémentaires est alors considérées pour chacune des étapes suivantes. Il est bien-sûr envisageable que la pièce 12 ne possède qu’une zone à grenailler avec une seule ou plusieurs zones élémentaires. La troisième étape C consiste à mettre au point une gamme de grenaillage. Par mettre au point, on entend, identifier au moins un ensemble de paramètres machine permettant de remplir les exigences de grenaillage explicité dans les spécifications au plan de la pièce. Cette troisième étape C consiste dans un premier temps de déterminer, pour chacune des zones élémentaires, une combinaison de paramètres de grenaillage d’une grenailleuse répondant aux spécifications de grenaillage pour ladite zone élémentaire considérée. Cela est réalisé en recherchant les paramètres de grenaillage de la grenailleuse, dans une première base de données préalablement établie liant les spécifications de grenaillage, en particulier l’intensité Almen et le taux de recouvrement auxdits paramètres de grenaillage de la grenailleuse. La première base de données sert également de capitalisation du savoir-faire de l’entreprise, notamment des expériences des professionnels métier sur la grenailleuse. Cette première base de données, pouvant être appelé abaque numérique, a été réalisé en interne. Cet abaque lie ainsi les paramètres machine influant : - les propriétés macroscopiques du jet de grenaille émis par la grenailleuse et les impacts de la grenaille sur les zones élémentaires (tel que par exemple la vitesse moyenne de projection, taille moyenne de cratère résultant de l’impact d’une grenaille, la fréquence des impacts sur une surface donnée), et - les exigences de grenaillage (tel que par exemple l’intensité Almen et recouvrement). Cet abaque a été réalisé pour une pièce plate de type massif semi-infini. Ainsi, pour une intensité Almen et un taux de recouvrement donné, les paramètres de la grenailleuse obtenus par le biais de l’abaque peuvent résulter en grenaillage sur une pièce avec une géométrie plus complexe (avec notamment des surfaces non planes) avec une intensité Almen et un taux recouvrement différents. Une interface utilisateur a été développée pour permettre l’utilisation de l’abaque. Elle permet ainsi d’obtenir de manière ergonomique, à partir des exigences de grenaillage, les paramètres de grenaillage de la grenailleuse. Cette interface permet ainsi d’identifier une plage restreinte de paramètres machine qui permettent d’atteindre des taux de recouvrement et des intensités Almen conformes aux exigences au plan. Cela est notamment réalisé à l’aide de la première base de données, construite à partir de données expérimentales de grenaillage et du savoir-faire de l’entreprise. La troisième étape C dans un deuxième temps comprend la simulation de la projection d’un jet de grenaille en chaque point de ladite zone et en déduire en chaque point un état de surface comprenant une information d’intensité Almen et de taux de recouvrement. Dans un premier temps, des volumes de simulation du jet 22 de grenaille sont définis. Ces volumes servent à délimiter l’espace dans lequel la grenaille (les billes de grenaillage) peut être simulée. Les dimensions et géométrie de ces volumes sont choisis de sorte à englober les zones élémentaires de la pièce 12 et le volume balayé par le jet 22 de grenaille suivant la trajectoire de la buse 24. Plusieurs volumes sont créés en fonction des stratégies de grenaillage et de la géométrie de la pièce 12, en d’autres termes du nombre et de la position des zones élémentaires définies dans l’étape A les unes par rapport aux autres. Les volumes sont créés à l’aide d’une copie du modèle tridimensionnel de la deuxième étape B. Une fois les volumes définis, un maillage surfacique est réalisé à l’aide de méthode des éléments finis. Les informations relatives au maillage (c’est-à-dire les sommets et des éléments surfaciques délimités par des segments reliant les nœuds) des différents volumes sont stockées dans différents fichiers. Dans un deuxième temps, la projection d’un jet de grenaille est simulée en chaque point de ladite zone et en déduire en chaque point un état de surface 24 comprenant une information d’intensité Almen et de taux de recouvrement. Pour cela, pour chacune des zones élémentaires définies à l’étape A, les paramètres 20 de la grenailleuse obtenue à l’aide de l’abaque sont utilisés comme données d’entrée de la simulation de l’instruction de la grenailleuse pour une zone donnée. Ainsi, dans cette étape, il est précisément simulé l’instruction de travail en cours pour la zone considérée avec les paramètres retenus en début d’étape. Une interface utilisateur 26 a été développée pour visualiser la simulation, comme illustrée en figure 5. On voit donc sur la figure 5 que la trajectoire de chacune des billes de la grenaille du jet 22 de la grenailleuse pour chaque instruction pour chacune des zones élémentaires de la pièce est ainsi simulée. Bien que non représenté, il possible d’obtenir une représentation du taux de recouvrement pour la zone considérée. La quatrième étape D du procédé comprend la combinaison des états de surface des zones élémentaires de la pièce afin d’obtenir une cartographie d’intensité Almen et de taux de recouvrement pour le modèle tridimensionnel. La simulation des instructions de grenaillage étant faite successivement ou parallèlement pour chacune des zones élémentaires, lorsque l’ensemble des simulations pour chacune des zones élémentaires est réalisée, l’ensemble des données sont consolidées pour analyser et visualiser les intensités Almen ainsi que les taux de recouvrement obtenus pour les paramètres de la grenailleuse sélectionnés. La figure 6 représente la recombinaison des résultats de simulation pour chacune des zones élémentaires. Pour faciliter l’analyse, une cartographie 30 de l’intensité Almen et une cartographie 32 du taux de recouvrement sur une coupe bidimensionnelle du modèle tridimensionnel de la pièce 12 est réalisé. Chaque couleur correspondant à une gamme d’intensité Almen ou de taux de recouvrement, il peut être facilement visualiser si l’intensité Almen et/ou le taux de recouvrement sont conformes aux exigences de spécifications. Des critères de conformité ou de non conformité peuvent également être affichés à l’utilisateur, à partir des résultats de simulation et des exigences au plan pour la zone en question. Une telle cartographie s’avère instructive, car elle permet d’obtenir des informations additionnelles notamment sur le recouvrement pour les jonctions séparant les zones élémentaires. Ainsi, il est vérifié localement vérifier si l’intensité Almen simulée et le taux de recouvrement sont conformes aux plages d’intensité de grenaillage et de taux de recouvrement compris dans les spécifications de grenaillage. En d’autres termes, on valide que les paramètres de la grenailleuse sélectionnés permettent un grenaillage de la pièce selon les exigences de grenaillage explicités dans les spécifications. Si localement, les intensités Almen et/ou le taux de recouvrement ne sont pas en adéquation avec les spécifications de grenaillage, les zones élémentaires dans lesquelles les valeurs d’intensité Almen et/ou le taux de recouvrement ne sont pas conformes sont identifiées. La troisième étape C et la quatrième étape D sont alors à nouveau exécutées pour les zones identifiées en question. Une fois que localement, toutes les intensités Almen et le taux de recouvrement sont conformes aux spécifications de grenaillage, la cinquième étape E du procédé est réalisée. La cinquième étape E, comprend la détermination du modèle tridimensionnel de la pièce 12 suite au grenaillage, en prenant en considération les déformations induites par le grenaillage. Pour cela dans un premier temps, le tenseur de contraintes est calculé à partir de la cartographie. Le tenseur de contraintes est ensuite appliqué à chacun des points constitutifs de la pièce 12. Le calcul du tenseur de contraintes est effectué sous l’hypothèse d’un massif semi infini, ce qui signifie que le tenseur de contraintes ne dépend que de la profondeur z. Ainsi, dans un premier temps, on calcule une distance z de chaque point d’intégration de GAUSS du maillage à un élément de surface du maillage le plus proche. Les éléments de surface sont notés E n pour n zones élémentaires. Les points d’intégration de GAUSS font partie intégrante de la méthode des éléments finis. En ces points d’intégration de GAUSS, plusieurs équations différentielles sont résolues (par intégration) de sorte à obtenir les contraintes résiduelles qui sont induites par le grenaillage et qui sont responsables des distorsions de pièces. La Figure 8 illustre la méthode de calcul de la distance à la surface des points d’intégration du maillage selon la méthode des éléments finis. Cela est réalisé pour chaque zone définie à la première étape A du procédé, et ceci en parallèle. La distance z, de chaque élément 34 de la pièce, correspond à la distance la plus courte à la surface parmi les éléments de surface E 1 , E 2 , E 3 du maillage situé à proximité. Dans l’exemple illustré, la distance z est déterminée à l’aide l’élément de surface E 2 , les distances z 1 et z 2 étant supérieure à la distance z 3 . La distance z est alors la distance entre le projeté orthogonal d’un point 36 de l’élément 34 dans la pièce et le point 36 de l’élément 34 dans la pièce comme cela est visible sur la figure 8. A partir de cette distance, un calcul est effectué pour chaque point du maillage bidimensionnel le profil de contraintes résiduelles selon z à partir de la combinaison de paramètres de la grenailleuse. Le profil de contraintes résiduelles peut par exemple, au lieu d’être calculé, être identifié dans un abaque reliant la combinaison de paramètres de la grenailleuse avec le profil de contraintes résiduelles . Un exemple de profil de contraintes résiduelles 38 est illustré à la figure 8. Avec ces éléments, pour déterminer dans quelle mesure le grenaillage simulé, impacte la géométrie de la pièce, deux méthodes peuvent être utilisées, pour réaliser les étapes suivantes du procédé : - Calcul de la déformation plastique induite par le grenaillage ; - Calcul de l’état initial à introduire dans la pièce ; -Calcul de l’équilibrage élastique. Ces deux méthodes ne sont que des méthodes parmi d’autres méthodes mathématiques connues pouvant également être appliqué au cas d’espèce. La première méthode est la méthode dite des « Eigen Strain », prenant en considération l’inhomogénéité de la dilatation thermique, suite aux opérations de grenaillage. La déformation plastique induite par le grenaillage est calculée pour chaque point de GAUSS du maillage. Ce calcul est réalisé à partir de l’équation suivante : Ensuite, le calcul de l’état initial à introduire dans la pièce consiste à déterminer le tenseur de contrainte équivalent. Le tenseur de contraintes résiduelles initial est calculé pour chaque point de maillage, correspondant aux points d’intégration de GAUSS, en fonction de la combinaison de paramètres de la grenailleuse, suivant l’équation suivante : Ce calcul du tenseur de contraintes résiduelles initial est effectué à l’aide d’une deuxième base de données, dans laquelle les paramètres de contrôle du grenaillage (intensité de grenaillage, intensité de grenaillage, etc.) sont liés au profil de contraintes résiduelles massif semi-infini (σ xx res (z)). Cette deuxième base de données sert également de capitalisation du savoir-faire de l’entreprise, notamment des professionnels de grenaillage. Une fois ces calculs réalisés, un fichier est généré, de sorte à contenir l’ensemble des informations nécessaires au calcul d’équilibrage statique de la pièce suite au grenaillage, c’est-à-dire la déformation plastique induite par le grenaillage ainsi que l’état initial. Cela permet donc, de calculer l’équilibrage élastique de la pièce, en d’autres termes de déterminer les contraintes résiduelles finales, selon l’équation suivante : La deuxième méthode pouvant être mise en œuvre est la méthode des champs thermiques fictifs. Cette méthode connue ne sera pas détaillée ici. Ensuite, un calcul mécanique d’équilibrage élastique est réalisé par la suite afin de déterminer les contraintes résiduelles finales, par la méthode des éléments finis, à l’aide de l’équation Math 5, la formule générique du tenseur des contraintes résiduelles finales. Au terme de ce calcul, on peut alors obtenir la modélisation de la déformation induite suite au grenaillage simulé, et alors le visualiser tel qu’illustré à la figure 7. La forme initiale 40 de la pièce avant grenaillage est illustrée et la forme de la pièce après grenaillage 42 est également illustrée. La sixième étape F consiste à comparer dimensionnellement le modèle tridimensionnel post-grenaillage au modèle tridimensionnel de la pièce. Un code couleur peut permettre d’identifier/visualiser dans quelle mesure ces déformations sont acceptables au regard des tolérances imposées pour la pièce en question. Si les tolérances de dimension ne sont pas conformes, le procédé est à nouveau réalisé à partir de l’étape B, pour les zones pour lesquelles les tolérances ne sont pas respectées. Au terme de ces étapes, le procédé est répété de sorte à générer l’ensemble de la gamme des paramètres pour lesquels les exigences de spécifications de grenaillage sont respectées. Suite au déroulement du détaillé dans le présent document, une dernière étape vise à valider la gamme des paramètres de grenaillage, obtenus suite à l’exécution du procédé, sur une pièce réelle. Ainsi, suite à la simulation, un grenaillage en conditions réelles est réalisé sur différentes pièces afin de valider les différentes gammes de paramètres de la grenailleuse obtenus. Cela permet de valider sur les résultats expérimentaux sont bien conformes avec les exigences de grenaillage spécifiées pour la pièce en question. Le procédé tel que décrit dans le document a vocation à être implémenté et mise en œuvre sur un ordinateur. En particulier, un programme d’ordinateur comprenant des instructions pour mettre en œuvre le procédé peut être ainsi développé. Ce programme peut ainsi être enregistré sur une mémoire d’un ordinateur ou un support de stockage de données lisible par ordinateur, permettant l’exécution du procédé sur ordinateur. Egalement, le présent document concerne un système de CAO comprenant un processeur couplé à une mémoire et une interface utilisateur graphique et apte à communiquer avec une grenailleuse, tel que le programme est enregistré en mémoire. Ainsi, le procédé est exécuté sur ordinateur et dès lors que des paramètres de la grenailleuse sont identifiés, ceux-ci peuvent être envoyés à la grenailleuse, de sorte à tester en condition réelles ces paramètres au regard des exigences de grenaillage de la pièce. La simulation permet ainsi de ne tester en condition réelles de grenaillage, que des paramètres de grenailleuse qui ont de bonnes chances de permettre un grenaillage conforme aux spécifications de grenaillage. Cela permet un gain de temps, et également un gain économique, grâce à la simulation. Procédé de détermination de paramètres de grenaillage par apprentissage pour une pièce le procédé comprenant les étapes suivantes : a) Fournir des spécifications de grenaillage de la pièce ; b) Fournir un modèle tridimensionnel de la pièce, réalisé par éléments finis à l’aide d’un maillage tridimensionnel de la pièce et subdiviser le modèle tridimensionnel en une pluralité de zones ; c) Pour chacune des zones : i. déterminer une ou plusieurs combinaisons de paramètres de grenaillage d’une grenailleuse répondant aux spécifications de grenaillage pour ladite zone considérée ; ii. simuler la projection d’un jet de grenaille en chaque point de ladite zone et en déduire en chaque point un état de surface comprenant une information d’intensité Almen et de taux de recouvrement ; d) combiner les états de surface des zones de la pièce afin d’obtenir une cartographie d’intensité Almen et de taux de recouvrement pour le modèle tridimensionnel ; e) déterminer un modèle tridimensionnel post-grenaillage de la pièce : i. Calculer un tenseur de contraintes à partir de la cartographie réalisée à l’étape d) ; ii. Appliquer le tenseur de contraintes au modèle tridimensionnel de la pièce pour obtenir le modèle tridimensionnel post-grenaillage ; f) comparer dimensionnellement le modèle tridimensionnel post-grenaillage au modèle tridimensionnel de la pièce. Procédé selon la revendication 1, dans lequel les spécifications de grenaillage comprennent un ou plusieurs des éléments suivants : i. des parties de la pièce à grenailler ; ii. des parties de la pièce interdites au grenaillage ; iii. des parties de la pièce pour lesquelles le grenaillage est admissible; iv. le type de grenaille ; v. des caractéristiques mécaniques et dimensionnelles de la grenaille ; vi. une plage d’intensité de grenaillage pour chacune des parties à grenailler ; vii. une plage de taux de recouvrement minimum et/ maximum pour chacune des parties à grenailler ; viii. profil de contraintes résiduelles minimal ; ix. rugosité surfacique maximum admissible ; x. des tolérances dimensionnelles de la pièce. Procédé selon la revendication 2, dans lequel l’étape c) i) consiste à rechercher les paramètres de grenaillage de la grenailleuse, dans une première base de données préalablement établie liant les spécifications de grenaillage, en particulier l’intensité Almen et le taux de recouvrement auxdits paramètres de grenaillage de la grenailleuse. Procédé selon l’une des revendications 1 à 3, dans lequel l’étape c) ii) est suivie d’une étape de vérification consistant à vérifier si l’intensité Almen simulée et le taux de recouvrement sont conformes aux plages d’intensité de grenaillage et de taux de recouvrement compris dans les spécifications de grenaillage. Procédé selon l’une des revendications 1 à 4, dans lequel l’étape d) est suivie d’une étape de vérification consistant à vérifier si l’intensité Almen simulée et le taux de recouvrement sont conformes aux plages d’intensité de grenaillage et de taux de recouvrement compris dans les spécifications de grenaillage. Procédé selon l’une des revendications 1 à 5, dans lequel pour réaliser l’étape e) : i. on calcule une distance z de chaque point d’intégration du maillage à un élément de surface du maillage le plus proche ; ii. on calcule le tenseur de contraintes résiduelles initial à appliquer au modèle tridimensionnel aux fins d’obtention du modèle tridimensionnel post-grenaillage, le calcul dudit tenseur étant réalisé pour tout point d’intégration du maillage en fonction de la distance z à l’élément de surface associé et du profil de contraintes résiduelles. Procédé selon la revendication 6, dans lequel l’étape i. utilise une deuxième base de données liant des paramètres de la grenailleuse à des profils de contraintes résiduelles. Procédé selon l’une des revendications 1 à 7, dans lequel pour la détermination du modèle tridimensionnel post-grenaillage comprend les étapes consistant à : - Calculer pour chaque point du maillage tridimensionnel par éléments finis de l’étape b), la distance z à la surface la plus proche ; - Calculer pour chaque point du maillage bidimensionnel le profil de contraintes résiduelles selon z à partir de la combinaison de paramètres de la grenailleuse ; - Calculer le tenseur des déformations plastiques correspondant ε p et de contraintes résiduelles initial pour chaque point de maillage, en fonction de la combinaison de paramètres de la grenailleuse; - Calculer l’équilibrage statique de la pièce par éléments finis à partir l’objet modélisé tridimensionnel de la pièce et de propriétés élastiques d’un matériau constitutif de la pièce ; - Calculer, pour chaque point du maillage tridimensionnel, la déformation géométrique induite du modèle tridimensionnel, après un équilibrage statique de la pièce, à partir du tenseur de contraintes résiduelles et de l’équilibrage statique. Programme d’ordinateur comprenant des instructions pour mettre en œuvre le procédé selon l’une des revendications 1 à 8 lorsque ledit programme est exécuté sur un ordinateur. Support de stockage de données sur lequel est enregistré le programme d’ordinateur selon la revendication 9. Système de CAO comprenant un processeur couplé à une mémoire et une interface utilisateur graphique et apte à communiquer avec une grenailleuse, tel que le programme selon la revendication 9 est enregistré en mémoire.