FR 3090460 A1 20200626 FR 1873368 A 20181219 Description Titre de l'invention : PRÉDICTION DE LA FORME D'UN OBJET TRIDIMENSIONNEL SOUMIS À UN PROCESSUS DE DIFFUSION Domaine technique 100011 La présente invention concerne de manière générale le domaine de la fabrication d'objets tri-dimensionnels à partir d'une forme initiale destinée à subir un processus de diffusion, et plus particulièrement une méthode de prédiction de la forme de tels objets au terme dudit processus. Technique antérieure [0002] De nombreuses techniques de fabrication additive d'objets tri-dimensionnels ont vu le jour ces dernières années. Ces techniques font généralement appel à un dépôt de couches successives subissant une transformation physique telle qu'une polymérisation par exemple. Ainsi, la demande EP-A-820855 décrit une méthode de fabrication d'objets tri-dimensionnels à partir de couches successives obtenues par stéréolithographie et subissant une polymérisation sous l'effet d'un rayonnement énergétique. [0003] Quelle que soit la technique utilisée, il est connu de procéder à une refusion (reflow) de l'objet formé par ces couches successives. Toutefois, la forme finale de l'objet après refusion est difficile à déterminer à l'avance, a fortiori lorsque le coefficient de diffusion du matériau et le temps de refusion sont élevés. 11 est donc nécessaire de procéder à des campagnes d'essais longues et coûteuses, avec des couches de formes différentes et de positions relatives différentes avant d'atteindre la forme finale souhaitée. [0004] Un but de la présente invention est par conséquent de proposer une méthode de prédiction de la forme d'un objet tri-dimensionnel après que celui-ci a subi un processus de diffusion. Un but subsidiaire de la présente invention est de proposer une méthode de fabrication d'un objet tri-dimensionnel de forme de consigne prédéterminée et de matériau donné, à partir d'une forme initiale donnée. Exposé de l'invention [0005] La présente invention est définie par une méthode de prédiction de la forme d'un objet tri-dimensionnel soumis à un processus de diffusion pendant une durée prédéterminée, ladite méthode comprenant : - une phase de calibration dans laquelle une pluralité d'échantillons de formes initiales différentes subissent, pendant ladite durée déterminée, le processus de diffusion, au terme duquel les échantillons adoptent des formes finales, et dans laquelle, pour chaque échantillon : 2 a) on mesure la hauteur initiale de l'échantillon et l'on acquiert les premiers contours d'une pluralité tranches horizontales de l'échantillon dans sa forme initiale, à des hauteurs normalisées par la hauteur initiale ; h) on mesure la hauteur finale de l'échantillon et l'on acquiert des seconds contours d'une même pluralité de tranches horizontales de l'échantillon dans sa forme finale, à des hauteurs normalisées par la hauteur finale, égales aux dites hauteurs normalisées par la hauteur finale ; c) on estime les paramètres d'une loi de morphing verticale permettant de passer de la hauteur initiale à la hauteur finale et on estime les paramètres d'une loi de morphing latérale permettant de passer d'une représentation, dans un espace de Fourier, des premiers contours à une représentation, dans ce même espace, des seconds contours ; - une phase de prédiction dans laquelle : d) on estime la hauteur finale dudit objet à partir de sa hauteur initiale et des paramètres de la loi de morphing verticale ; e) on obtient une représentation dans ledit espace de Fourier d'une pluralité de contours des tranches de l'objet dans sa forme initiale, aux dites hauteurs normalisées ; f) on obtient une représentation dans ledit espace de Fourier d'une pluralité de contours des tranches de l'objet dans sa forme finale, aux hauteurs normalisées, à partir de ladite représentation des tranches de l'objet dans sa forme initiale et des paramètres de la loi de morphing latérale ; g) on détermine les contours de la pluralité de tranches dudit objet dans sa forme finale, aux hauteurs normalisées, à partir de la représentation des contours de ces tranches dans ledit espace de Fourier obtenue à l'étape (f) ; h) on obtient la forme finale de l'objet tridimensionnel à partir de la hauteur finale dudit objet et des contours des tranches dudit objet dans sa forme finale, précédemment déterminées aux hauteurs normalisées. [0006] Selon un premier mode de réalisation la représentation dans ledit espace de Fourier des premiers et seconds contours à l'étape (c) ainsi que la représentation, dans ce même espace, des contours des tranches de l'objet dans sa forme initiale à l'étape (f), est obtenue par une transformée de Fourier curviligne de ces contours. [0007] Réciproquement, à l'étape (g), la détermination des contours de la pluralité de tranches dudit objet dans sa forme finale est obtenue par une transformation de Fourier curviligne inverse de la représentation des contours de ces tranches dans ledit espace de Fourier obtenue à l'étape (f). [0008] Selon un second mode de réalisation la représentation dans ledit espace de Fourier des premiers et seconds contours à l'étape (c) ainsi que la représentation dans cet espace des contours des tranches de l'objet dans sa forme initiale à l'étape (f), est obtenue en effectuant les étapes suivantes : 3 (f1) on sélectionne une pluralité de points de ce contour, et on détermine les normales de ce contour en ces points ; (f2) on calcule les coefficients d'une transformée de Fourier spatiale bidimensionnelle d'une fonction de courbe de niveau prenant une valeur nulle aux points ainsi sélectionnés et dont les normales en ces points sont respectivement égales aux normales du contour en ces mêmes points. [0009] Réciproquement, à l'étape (g), la détermination des contours de la pluralité de tranches dudit objet dans sa forme finale est obtenue en effectuant une transformée de Fourier spatiale bidimensionnelle inverse des coefficients calculés à l'étape (f2). [0010] L'invention est également définie par une méthode de fabrication d'un objet tri- dimensionnel de matériau donné et ayant une forme de consigne donnée, la méthode de fabrication comprenant une phase de calibration dans laquelle on construit à partir de mesures d'échantillons de cc matériau, avant et après un processus de diffusion de durée T, une base de données de paramètres d'une loi de morphing verticale et d'une loi de morphing latérale pour une pluralité de hauteurs normalisées desdits échantillons, puis on applique une méthode de prédiction telle que définie ci-dessus à une pluralité de formes initiales pour prédire leurs formes finales respectives au terme du processus de diffusion, et l'on sélectionne la forme initiale dont la forme finale correspondante est la plus proche de ladite forme de consigne au sens d'une distance prédéterminée, on réalise un objet tridimensionnel ayant la forme initiale ainsi sélectionnée et l'on applique enfin le processus de diffusion pendant la durée T pour obtenir un objet tridimensionnel ayant ladite forme de consigne. Brève description des dessins [0011] D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture d'un mode de réalisation préférentiel de l'invention, décrit en référence aux figures jointes parmi lesquelles : [0012] [fig.1] représente, de manière schématique, la forme initiale d'un objet tridi- mensionnel ; [0013] [fig.2] représente, de manière schématique, la forme de l'objet tri-dimensionnel de la Fig. 1 après que celui-ci a subi un processus de diffusion ; [0014] [fig.3A] et [0015] [fig.3B] représentent l'ordinogramme d'une méthode de prédiction, selon un premier mode de réalisation, de la forme d'un objet tri-dimensionnel ayant subi un processus de diffusion ; [0016] [fig.4A] et [0017] [fig.4B] représentent l'ordinogramme d'une méthode de prédiction de la forme d'un objet tri-dimensionnel ayant subi un processus de diffusion, selon une variante du premier mode de réalisation ; [0018] [fig.5À] et [0019] [fig.SB] représentent l'ordinogramme d'une méthode de prédiction, selon un second mode de réalisation, de la forme d'un objet tri-dimensionnel ayant subi un processus de diffusion ; [0020] [fig.6À] et [0021] [fig.6B] représentent l'ordinogramme d'une méthode de prédiction de la forme d'un objet tri-dimensionnel ayant subi un processus de diffusion, selon une variante du second mode de réalisation ; [0022] [fig.7] représente l'ordinogramme d'une méthode de fabrication d'un objet tri- dimensionnel de matériau donné et de forme de consigne prédéterminée. [0023] EXPOSÉ DÉTAILLÉ DE MODES DE RÉALISATION PARTICULIERS [0024] Nous considérerons dans la suite un objet tri-dimensionnel de forme initiale donnée et de matériau donné. Cet objet peut par exemple avoir été réalisé au moyen d'une technique de fabrication additive ou par photolithographie 3D. Le matériau peut être par exemple un polymère. [0025] On a représenté en Fig. 1 une vue en perspective d'un tel objet, de hauteur totale hm«,. Cet objet peut être divisé en tranches horizontales successives désignées par e chaque tranche étant caractérisée par son propre contour dans le plan ho- rizontal. [0026] On suppose que cet objet subit une transformation au moyen d'un processus de diffusion d'origine physique voire physico-chimique. Ce processus de diffusion peut être par exemple une refusion (rellow), telle qu'une refusion thermique du matériau constitutif de l'objet ou encore une érosion de ce matériau en présence d'un solvant en phase gazeuse. [0027] La Fig. 2 représente l'objet tri-dimensionnel de la Fig. 1 après que celui-ci a subi un processus de diffusion pendant un temps T. En raison du processus de diffusion et de la gravité, la forme de l'objet s'est modifiée : la forme s'est globalement affaissée et les zones correspondant à des fréquences spatiales élevées ont été lissées. [0028] L'idée à la base de la présente invention est de modéliser la transformation de forme (encore dénommée morphing) de l'objet tri-dimensionnel. Il est important de comprendre que cette modélisation n'a pas pour objectif de modéliser le processus physique sous-jacent mais simplement le changement de forme résultant de l'application du processus. [0029] De manière originale, l'évolution de la forme de l'objet est modélisée, d'une part, par une première loi donnant la hauteur totale de l'objet dans le temps et, d'autre part, par une seconde loi décrivant l'évolution du contour d'une tranche horizontale de l'objet dans le temps. [0030] La hauteur de l'objet tri-dimensionnel dépend du temps et est notée pour cette raison h(t), h(0) = han, désignant la hauteur initiale de l'objet et h(1) la hauteur de ce même objet après qu'il a subi un processus de diffusion pendant un temps T. [0031] On a pu montrer que, lorsque l'objet présente un axe de symétrie vertical, la hauteur knde l'objet au terme du processus de diffusion variait en fonction de la hauteur initiale de l'objet et du rayon ro du cercle équivalent à la base de cet objet. Par cercle équivalent à la hase de l'objet, on entend un cercle dont l'aire est égale à la surface de la base de cet objet. [0032] [Math .1] avec : [0033] [Math.2j hr .r (t- où a(i), b(7') sont des coefficients strictement positifs dépendant du matériau utilisé et de la durée du processus de diffusion T notés plus simplement dans la suite a4./5 L'expression fl (ro T) est dénommée dans la suite facteur de morphing vertical relatif à la durée T [0034] On considère ensuite une tranche horizontale de l'objet, située à une hauteur Z par rapport à sa base. On note 2 = z / h la hauteur normalisée de la tranche. Cette hauteur normalisée reste inchangée au cours du processus de diffusion. On note 0(t) le contour de la tranche, celui-ci évoluant pendant l'application du processus de diffusion. Le contour 16 (t) représenté comme une fonction de son abscisse curviligne u et peut être développé en série de Fourier curviligne, à savoir : [0035] [Math.3] (1)5t (3-1) [0036] [Math.4] (3-2) où (x0 4. (1, y(1_0) sont les coordonnées cartésiennes d'un point P(a,/)du contour ei) d'abscisse curviligne et ak bk(0, ck(t), d[ (t) sont les coefficients de Fourier 6 correspondant à la fréquence spatiale k. On notera que ces coefficients de Fourier dépendent du temps dans la mesure où le contour 0(t) de la tranche évolue dans le temps en raison de l'application du processus de diffusion. [0037] Les équations (3-1) et (3-2) peuvent être cxpri niées de manière plus compacte au moyen de coordonnées complexes : [0038] [Math.5] 0 OÙ OU = À-0 (I) ± ..020 (I) représente l'évolution dans le temps du centre de gravité du contour, 01Q) représente l'évolution dans le temps de l'ellipse équivalente, c'est-à-dire approximant au mieux le contour de la tranche. [0039] Lorsque la tranche est isolée, le centre de gravité du contour n'a pas de raison de se déplacer et l'on peut supposer, sans perte de généralité, que celui-ci est confondu avec l'origine, autrement dit 4(X)=-0. [0040] Lc processus de diffusion a tendance à rendre le contour de la tranche circulaire. Autrement dit, l'ellipse équivalente tend vers un cercle et : [0041] [Math.6] où lifi(i) est le rayon du contour de la tranche considérée à la hauteur normalisée 2 au bout d'un temps infini. [0042] De manière plus générale, on suppose qu'au bout d'un temps infini la tranche se contracte latéralement d'un coefficient s(1): [0043] [Math.7] [0044] Le coefficient 42) sera dénommé dans la suite coefficient de contraction. Le co- efficient de contraction dépend de la hauteur normalisée à laquelle se trouve la tranche dans l'objet tri-dimensionnel. En effet, en raison du processus de diffusion la forme de l'objet au bout d'un temps infini est celle d'une sphère tronquée. Il en résulte que s(1) = 0 , puisque le sommet de la sphère peut. être considéré comme une tranche de taille nulle. En revanche, le coefficient 40) peut être supérieur à 1, autrement dit la base peut se dilater dans le temps. [0045] Le contour de la tranche tendant. vers un cercle, la puissance des harmoniques de rang 7 > 1 tend vers zéro lorsque le temps tend vers l'infini, c'est-à-dire : [0046] [Math.8] [0047] On suppose que la puissance des harmoniques décroit de manière exponentielle dans le temps, la rapidité de la décroissance dépendant du rang k de l'harmonique, à savoir : [0048] [Math.9] où a et 7 sont des réels strictement positifs. [0049] En résumé, nous pouvons donc écrire de manière condensée : [0050] [Math.10] (9) avec kb1(1,ic)1- s(2)-10,(0)let OAK = (1 - (t - 1)0, )+ (1+ sa +1» ,(7) où est le symbole de Dirac. [0051] De manière équivalente, on a : [0052] [Math.11] où m, (t) = exp (-a H t) est dénommé facteur de morphing latéral pour l'harmonique de rang k 100531 Si l'on considère maintenant non plus une tranche isolée mais l'ensemble des tranches de l'objet, H faut tenir compte du fait que les centres de gravité de ces différentes tranches, notés e(!), où f représente l'indice de la tranche, ne sont pas alignés à l'instant initial/ = 0. De surcroît, les tranches peuvent interagir entre elles et les centres de gravité des différentes tranches peuvent dériver dans le temps. 100541 Dans des cas de formes simples, on peut légitimement supposer que les tranches sont indépendantes et leurs centres de gravité restent fixes dans le temps. Il est alors possible de décrire l'évolution de ces tranches dans le temps par : 100551 IMath.121 où A I; (t) = exp( a" Ikr 1) est le facteur de morphing relatif à la tranche et où a if ) 8 et sont les coefficients a et relatifs à la couche f',10[(,6)1= s ')-19K(0)let - - - 1)0i' (x)) + -8(k + I) O` (x:) où g(e) est le facteur de contraction relatif à la couche e. [0056] Cependant, dans d'autres cas, notamment lorsque les contours des tranches sont concaves, par exemple se présentent sous la forme de U, l'hypothèse d'indépendance des tranches n'est pas réaliste. Dans des cas extrêmes, cette hypothèse d'indépendance conduit à des résultats aberrants comme celui par exemple où une tranche supérieure ne serait pas supportée entièrement par une tranche inférieure. [0057] De manière originale, la dépendance entre tranches est traitée comme un problème de minimisation d'énergie, l'énergie d'interaction entre deux tranches successives étant élevée quand celles-ci ne se chevauchent pas et faible lorsqu'elles se chevauchent. [0058] On note que, lorsque les tranches sont indépendantes, l'évolution dans le temps des harmoniques curvilignes des différents contours peut être obtenue à partir de l'équation (8) [0059] [Math.13] (12) [0060] En revanche, lorsque deux tranches successives fi et e2 sont couplées entre elles, l'évolution des contours de ces tranches est régie par les équations différentielles couplées : [0061] [Math.14] [0062] [Math.l5] [0063] Ces équations peuvent s'écrire de manière plus compacte sous la forme matricielle : [0064] [Math.16] (14) OÙ eitek (t) est un vecteur de taille L dont les éléments sont les harmoniques curvilignes 9 de rang k: Ak est une matrice symétrique tri-diagonale de taille L x L dont les termes de la diagonale principale sont égaux à cx()Ikle et dont les termes de la diagonale inférieure ainsi que de la diagonale supérieure sont égaux à ek, où ace) et yce) sont des paramètres strictement positifs et C' i. est un coefficient de couplage pour l'harmonique de rang k entre deux tranches successives de l'échantillon. [0065] En absence de couplage, c'est-à-dire lorsque les tranches sont indépendantes, la matrice Ak se réduit à une simple matrice diagonale. A l'inverse, le couplage peut s'étendre au-delà de la tranche immédiatement inférieure et de la tranche immédiatement supérieure, auquel cas la matrice comporte des sous-diagonales supérieures et inférieures composées de termes non nuls. En pratique, toutefois, un couplage avec la tranche immédiatement inférieure et de la tranche immédiatement supérieure est généralement suffisant pour décrire l'évolution des contours des différentes tranches de l'objet. [0066] Il résulte de l'équation (14) que : [0067] [Math.17] >ffik (k), ('S) où e avec matrice identité de taille L x L. exp(Akt)= (A,) 7 [0068] De manière similaire, avec les mêmes conventions que précédemment : [0069] [Math.18] (1) (1) = e )41>110)-+-:( (16) 100701 Le modèle ci-dessus permet de prédire la forme d'un objet tri-dimensionnel, après avoir subi un processus de diffusion pendant un temps prédéterminé T [0071] Les Fig. 3A-3B représentent l'ordinogramme d'une méthode de prédiction, selon un premier mode de réalisation, de la forme d'un objet tri-dimensionnel ayant subi un processus de diffusion pendant un temps prédéterminé. 100721 La méthode de prédiction comporte tout d'abord une phase de calibration comprenant les étapes 310-345 en Fig. 3A. 100731 A l'étape 310, on dispose d'une pluralité d'échantillons d'objets tri-dimensionnels de hauteurs différentes. Ces échantillons présentent également avantageusement des profils différents. 100741 A l'étape 315, on mesure la forme initiale de chacun des échantillons, par exemple au moyen d'un microscope à force atomique (AFM). La hauteur de chaque échantillon est mesurée et divisée en L intervalles, chaque intervalle définissant une tranche horizontale C de l'objet située à une hauteur normalisée, 5"(. On obtient à partir de la mesure de la forme initiale, le contour de chaque tranche e = 1 , . A l'étape 320, on effectue une décomposition de Fourier curviligne du contour de chacune de ces tranches. [0075] A l'étape 325, on applique aux différents échantillons le processus de diffusion pendant une durée T. [0076] A l'étape 330, on mesure la forme finale de chacun des échantillons, de préférence par les mêmes moyens que précédemment. On en déduit les nouvelles hauteurs des échantillons et les contours des tranches situées aux mêmes hauteurs normalisées que précédemment. [0077] A l'étape 335, on estime les coefficients de la loi de morphing verticale. Plus pré- cisément, à partir des hauteurs mesurées à l'étape 330 et des rayons équivalents des tranches formant la base des échantillons, obtenues à l'étape 315, on détermine les coefficients a, h de l'expression (2), à l'aide d'une technique d'ajustement de courbe. [0078] On estime ensuite la loi de morphing latérale. [0079] Pour ce faire, on effectue à l'étape 340 une décomposition de Fourier curviligne de chacun des contours des tranches C. =1,.., L, tels qu'obtenus à l'étape 330, cc pour chacun des échantillons. [0080] A l'étape 345, on détermine par une technique d'ajustement de courbe (eurve fating) pour chaque tranche t les paramètres sm, am, y( qui minimisent l'erreur quadratique entre les valeurs 0, (T) données par le modèle (expression (11) et les valeurs expérimentales correspondantes. On obtient ainsi une estimation de ces paramètres [0081] Avantageusement, on pourra utiliser une régression polynomiale pour estimer ces pa- ramètres. Plus précisément, on a observé de manière heuristique que les paramètres s, a, 7 pouvaient s'exprimer en fonction de la hauteur normalisée ±7= z /km, où I est la hauteur de la tranche considérée et =h(0) est la hauteur de l'échantillon initial. [0082] [Math.19] (17-1) [0083] [Math.20] 100841 11 [Math.21] (17-3) où P, (g), Pt 0) sont respectivement des polynômes de degré 2, 3 et 4 de la hauteur normalisée 2.On remarque que, dans l'expression (17-3), le coefficient de contraction s'annule pour 2 =1. Autrement dit, dès lors que la durée T est suffisamment longue, la tranche au sommet de l'objet tri-dimensionnel se réduit à un point. [0085] En pratique, les valeurs des paramètres s,',7 sont estimées pour une pluralité L de hauteurs normalisées (le cas échéant plusieurs estimations relatives à une même hauteur normalisée pouvant être issues de différents échantillons) et les coefficients des polynômes sont ensuite calculés pour minimiser l'écart quadratique moyen aux valeurs des paramètres ainsi estimées [0086] La phase de calibration est suivie d'une phase de prédiction proprement dite, comprenant les étapes 350-370 en Fig. 3B. [0087] La phase de prédiction vise à prédire la forme d'un objet tri-dimensionnel (de forme initiale donnée) lorsque celui-ci aura subi le processus de diffusion pendant une durée T. La forme initiale de l'objet est caractérisée par sa hauteur totale, hrtx, d'une part, et par les contours des tranches horizontales de l'objet situées aux valeurs de hauteurs normalisées g e, d'autre part. [0088] A l'étape 350, on calcule le rayon équivalent, c;, de la tranche (e =1)à la base de l'objet et l'on détermine, à partir de la loi de morphing verticale, la hauteur de l'objet au terme du processus de diffusion : [0089] [Math.22] [0090] A l'étape 355, on considère les tranches ; aux hauteurs normalisées ge, A cette étape, on effectue une transformée de Fourier curviligne des contours de ces tranches et l'on obtient les coefficients complexes 0' (0),k = N , . 100911 A l'étape 360, on détermine à partir de la loi de morphing latérale, les contours des tranches L au terme du processus de diffusion. Plus précisément, on calcule les coefficients 0 (T) des harmoniques de rang IkI >1 au moyen de : 100921 IMath.231 `; avec ild;(T)= exp (libration. 19) ) où 6/1', firf) ont été obtenus dans la phase de ca- ' Ikl 12 [0093] De manière similaire, on calcule les coefficients 01" ,(.0 pour la fréquence fon- damentale =1: [0094] [Math.24] 0" (20) où les coefficients de contraction r ), (2= 1,..,L, ont été obtenus dans la phase de calibration. [0095] A l'étape 365, on effectue une transformée curviligne inverse pour calculer le contour de chaque tranche : [0096] [Ma th.25] [0097] Enfin, à l'étape 370, on reconstruit la forme de l'objet-tridimensionnel à partir des contours des tranches C obtenus à l'étape 365, ces tranches étant désormais situées respectivement aux hauteurs h (T)f =1,.., L.Les Fies. 4A-4B représentent de manière schématique l'ordinogramme d'une méthode de prédiction de la forme d'un objet tri-dimensionnel ayant subi un processus de diffusion, selon une variante du premier mode de réalisation. [0098] Dans cette variante, on ne suppose plus que les tranches d'un l'objet tri-dimensionnel sont indépendantes et l'on effectue la calibration conjointement sur toutes les tranches de chaque échantillon. [0099] Comme précédemment, la méthode de prédiction comporte une phase de calibration (Fig. 4A) avant la phase de prédiction proprement dite (Fig. 4B). [0100] Les étapes 410-440 sont respectivement identiques aux étapes 310-340 et leur des- cription ne sera donc pas reprise ici. [0101] A l'étape 445, on estime les coefficients =1,..,L et C'A. vérifiant pour chacun des échantillons : [0102] [Math.26] exp(Az T).(11 (22) où la matrice tri-diagonale A, a pour termes diagonaux a, o _ Â et pour éléments de la diagonale inférieure et de la diagonale supérieure, le terme constant C., e k (in est un vecteur de taille L dont les éléments sont les harmoniques curvilignes de rang k des contours des différentes tranches de l'échantillon au terme de l'application du processus de diffusion et e k (0) est un vecteur de taille L dont les éléments sont les harmoniques curvilignes de rang k des contours des différentes tranches de 13 l'échantillon initial. [0103] En pratique, la résolution du système d'équations (22) se fait par ajustement de courbe en ne considérant que les premiers termes du développement en série de exp(AkT). [0104] La phase de prédiction proprement dite comprend les étapes 450-470. Elle vise comme précédemment à prédire la forme d'un objet tri-dimensionnel (de forme initiale donnée) lorsque celui-ci aura subi le processus de diffusion pendant une durée T. [0105] L'étape 450 de détermination de la hauteur de l'objet tri-dimensionnel, au moyen de la loi de morphing verticale dont les coefficients ont été estimés dans la phase de calibration, est identique à l'étape 350. [0106] A l'étape 455, on obtient les coefficients complexes 0s,(0), de la transformée de Fourier curviligne des contours des tranches CE, f =1,..,L,de l'objet tri-dimensionnel comme à l'étape 355, puis on forme des vecteurs 0,(o), Ikl>1, dont les éléments sont les coefficients O',:(0), À° =1,..,L. [0107] A l'étape 460, on détermine à partir de la loi de morphing latérale, les contours des tranches ;, À° =1,..,L au ternie du processus de diffusion. Plus précisément, on calcule les vecteurs e,'. (4'), 141, 1 au moyen de : [0108] [Math.27] (23) où la matrice Ak, a été déterminée dans la phase de calibration. [0109] Les éléments de chaque vecteur 4 (n, k I > i, donnent les coefficients des har- moniques de rang k,01,(n, dès différentes tranches C [0110] Par ailleurs, on calcule les coefficients 95';,(T),des différentes tranches (I: =1,.., L, pour la fréquence fondamentale = I à partir de l'expression (20) comme à l'étape 360. [0111] Les étapes 465 et 470 sont respectivement identiques aux étapes 365 et 370. Autrement dit, on calcule en 465 la transformée de Fourier curviligne inverse de 0"),:(T) ,k =-N,.'N, pour obtenir les contours des différentes tranches et on reconstruit enfin la forme (le l'objet en 470. [0112] Selon un second mode de réalisation de l'invention, la représentation des contours des tranches de l'objet tri-dimensionnel ne fait pas appel à une décomposition en série de Fourier curviligne mais à une méthode de surfaces de niveau, développées en coefficients de Fourier. De manière générale, les méthodes par surfaces de niveau permettent de représenter l'évolution d'une interface dans le temps. Une interface F à une dimension ou ligne de niveau, bornant à l'instant t une région 0dans un plan est définie par une équation du type : [0113] [Math.28] 14 (24) où la fonction P est une fonction lisse, positive à l'intérieur de la courbe 17 (donc sur Q) et négative à l'extérieur de celle-ci. [0114] On considère maintenant un instant to donné et l'on note f(x,y) la fonction dé- finissant la courbe f à cet instant. Connaissant un certain nombre de points appartenant à cette courbe et les normales à la courbe en ce point, il est possible de calculer les coefficients de Fourier (spatiaux) de la fonction f(xy). [0115] En effet, soit un ensemble de N points =1,..,N, de coordonnées (x' y") ap- partenant à la courbe F et soient N = (V, pl les vecteurs normaux à cette courbe en ces points, la fonction ,fix,y) doit vérifier les contraintes : 101161 IMath.291 101171 1Math.301 es:e: [0118] La fonction f(x,y) peut être développée en série de Fourier spatiale : [0119] [Math.31] (25-2) (25) avec A->: = 2.7z7 L, et k = Ly où L. et sont les dimensions d'un rectanele incluant la courbe, autrement dit k,. et k,' les fréquences spatiales les plus basses selon l'axe des -k" et des Y pour décrire la courbe r 101201 .11 résulte de l'expression (26) que : 101211 IMath.321 ('27-1 [0122] [Math.33] (:27-2) [0123] Les contraintes (25-1) et (25-2) peuvent alors se réécrire sous la forme de 3N équations que doivent vérifier les coefficients de Fourier .f'_' : [0124] [Math.34] [0125] [Math.35] (28-2 101261 IMath.361 -3) [0127] Si le nombre N de points est suffisant, on peut déterminer les coefficients pour un nombre prédéterminé d'harmoniques. Plus précisément, si l'on souhaite aller jusqu'aux harmoniques de rang ce,,P,, on aura besoin de N > 4 a p points sur la courbe F. - 3 [0128] On supposera dans la suite que a. = P, et l'on notera Q = = 21jy. [0129] Les équations précédentes (28-1), (28-2), (28-3) peuvent être représentées sous une forme matricielle plus compacte : [0130] [Math.371 FE, FM/ (29) où F est une matrice de taille QxQ dont les éléments sont les coefficients de Fourier, , 4. est une matrice de taille Q80 dont les éléments sont les valeurs exp(IkvErnx, + fkn-v,) Wà est une matrice de taille QxQ dont les éléments sont les valeurs jk-inexp(iknix, ik,ny W ( est une matrice de taille QxQ dont les éléments sont les valeurs jkynexp (Pym + ç? est un vecteur de taille Q dont les éléments sont nuls, v 7 est un vecteur de taille Q dont les éléments sont tous égaux à v7, est un vecteur de taille Q dont les éléments sont tous égaux à V,. [0131] Lorsque la courbe F évolue dans le temps, les coefficients de Fourier L. sont également fonction du temps et sont alors notés ft',,(0. 16 [0132] Si la courbe F est le contour d'une tranche horizontale d'un objet tri-dimensionnel, le phénomène de diffusion induit un lissage de la courbe de sorte que la puissance des harmoniques f', '(Otend vers zéro, autrement dit, dès lors que 7m2 + n2 > i f ,'(0 = 0. 101331 Comme dans le premier mode de réalisation, on suppose que la puissance des har- moniques détroit de manière exponentielle dans le temps, la rapidité de la décroissance dépendant du rang de l'harmonique, à savoir : [0134] [Math.38] (30) [0135] Ceci conduit à une évolution des haimoniqucs du type de celle de l'équation (10) : [0136] [Math.39] (31) [0137] f',.,,(00)1= s(f).1.f',,,(0)151 Jin et ' (30) 0 si \/', où s(i) est ic facteur de contraction de la tranche comme dans le premier mode de réalisation, et où le facteur de morphing latéral M7,,,,,, (t) est donné par : [0138] [Math.40] [0139] Lorsque l'on peut considérer les tranches e de l'objet tri-dimensionnel comme indépendantes l'évolution dans le temps des harmoniques des différents contours peut être obtenue à partir de l'équation (30) : [0140] [Math.41] (33) où f;:, (I) est le coefficient de Fourier (t) du contour d'une tranche f de l'échantillon. [0141] En revanche, lorsque les tranches sont dépendantes, on peut utiliser un modèle de couplage entre harmoniques de même rang de tranches successives, de manière similaire à (13-1) et (13-2) : [0142] 17 [Math.42] t [0143] [Math.43] tiet- ? qui peut également s'écrire sous forme matricielle : 101441 IMath.441 (35) où F', n(t) est un vecteur dc taille L dont les éléments sont les harmoniques spatiales d'indices in, n des différentes couches : f:'(t), est une matrice symétrique tri- diagonale dc taille L x L dont les termes de la diagonale principale sont égaux à kr'° _cm, et dont les termes de la diagonale inférieure ainsi que de la diagonale supérieure sont égaux à (7m.,. 101451 11 en résulte que le vecteur Fm, ()peut s'écrire : 101461 IMath.451 [0147] De manière similaire, avec les mêmes conventions que précédemment : [0148] [Math.46] ) [0149] [Math.47] ) e'p(A, F e\p(A: -2) [0150] Ce modèle, comme celui à la hase du premier mode de réalisation, permet de prédire la forme d'un objet tri-dimensionnel, après avoir subi un processus de diffusion pendant un temps prédéterminé T. [0151] Les Figs. 5A-5B représentent l'ordinogramme d'une méthode de prédiction, selon le second mode de réalisation, de la forme d'un objet tri-dimensionnel ayant subi un 18 processus de diffusion pendant un temps prédéterminé. [0152] La méthode de prédiction comporte tout d'abord une phase de calibration comprenant les étapes 510-545 en Fig. SA. [0153] A l'étape 510, on dispose d'une pluralité d'échantillons tri-dimensionnels de hauteurs différentes. Ces échantillons présentent également avantageusement des profils différents. [0154] A l'étape 515, on mesure la forme initiale de chacun des échantillons, par exemple au moyen d'un microscope à force atomique (AFM). La hauteur de chaque échantillon est mesurée et divisée en L intervalles, chaque intervalle définissant une tranche horizontale (7, de l'objet située à une hauteur normalisée, 5"(. On choisit sur le contour F de chaque tranche Ce, P =1,..,L, une pluralité N de points F', i , de co- ordonnées (x;, y;) appartenant à la courbe F, et on détermine les vecteur normaux N' rm,^ ) les vecteurs normaux à la courbe F en ces points. On choisira avanta- , geusement les points P,' équidistribués le long du contour. [0155] A l'étape 517, on détermine les courbes de niveau ,f (x, y), £ =1,..,L, vérifiant les contraintes : 101561 IMath.481 [0157] [Math.49] -2) [0158] A l'étape 520, on calcule, pour chacun (les échantillons, les coefficients de Fourier spatiaux f '(yA 2(0) de .f (.r, y), où la valeur D indique que l'objet tri-dimensionnel est considéré à l'instant initial avant l'application du processus de diffusion. Pour le calcul des coefficients de Fourier, on choisit les fréquences spatiales fondamentales k =2z /4,, et k,. = 2r t y où L, et 1/, sont les dimensions d'un rectangle incluant la projection de l'objet sur un plan horizontal. A l'étape 525, on applique aux différents échantillons le processus de diffusion pendant une durée T. [0159] A l'étape 530, on mesure la forme finale de chacun des échantillons par les mêmes moyens que précédemment. On en déduit les nouvelles hauteurs des échantillons et les contours des tranches situées aux mêmes hauteurs normalisées que précédemment. [0160] A l'étape 535, on estime les coefficients de la loi de morphing verticale. Plus pré- cisément, à partir des hauteurs mesurées à l'étape 530 et des rayons équivalents des tranches fondant la base des échantillons, obtenues à l'étape 515, on détermine les co- 19 efficients (-1,1) de l'expression (2), à l'aide d'une technique d'ajustement de courbe. [0161] On estime ensuite la loi de morphing latérale. [0162] Pour ce faire, à l'étape 540, on choisit à nouveau sur le contour r t de chaque tranche C', e =1,.., L, une pluralité N de points r/,i=1,..,N de coordonnées (2c',y,) appartenant à courbe F g et on détermine les vecteur normaux N =k'',v'-`l les vecteurs normaux à la courber f en ces points. On comprendra, bien que nous ayons ici adopté pour des raisons de simplification les mêmes notations qu'à l'étape 515, qu'il s'agit ici des contours après diffusion, à l'instant t = T.On détermine les courbes de niveau (x, y), vérifiant les contraintes (38-1) et (38-2), cette fois ci pour les points et les vecteurs normaux N obtenus à cette même étape, c'est à dire ceux relatifs aux contours des tranches à l'instant t = T. [0163] A l'étape 543, on calcule les coefficients de Fourier spatiaux je1:1 (r) des fonctions (x, y), où la valeur T indique que l'objet tri-dimensionnel est considéré après l'application du processus de diffusion. On utilise à cette fin les mêmes fréquences spatiales fondamentales k,., qu'à l'étape 520. [0164] A l'étape 545, on détermine par une technique d'ajustement de courbe (curve pour chaque tranchet les paramètres s(f) a() ,rçf) qui minimisent l'erreur quadratique entre les valeurs1;;;;,(T) données par le modèle (expression (31) et les valeurs expérimentales correspondantes. On obtient ainsi une estimation de ces paramètres : j;cet 10165] Avantageusement, on pourra utiliser une régression polynomiale comme dans le premier mode de réalisation pour estimer ces paramètres. Cette variante ne sera pas à nouveau décrite ici. 10166] La phase de cal ibration est suivie d'une phase de prédiction proprement dite, comprenant les étapes 550-570 en Fig. 5B. 10167] La phase de prédiction vise à prédire la forme d'un objet tri-dimensionnel (de forme initiale donnée) lorsque celui-ci aura subi le processus de diffusion pendant une durée T. 10168] A l'étape 550, on calcule le rayon équivalent, ro', de la tranche (P, =1)à la hase de l'objet et l'on détermine, à partir de la loi de morphing verticale, la hauteur de l'objet au terme du processus de diffusion, comme dans l'expression (18). 10169] A l'étape 553, on considère les tranches C ; aux hauteurs normalisées ge, e A cette étape, pour chaque tranche Cc, on sélectionne une pluralité N de points , de coordonnées (x', y';') sur le contour F , de chaque tranche et l'on détermine les vecteurs normaux à la courbe, ^ en ces points. On choisira avantageusement les points r'; angulairement équidistribués sur [0,2z]. [0170] A l'étape 555, on détermine les courbes de niveau f (x, = L, vérifiant les contraintes : [0171] [Math.50] (39-1 [0172] [Math.51] ]0l73] A l'étape 557, on calcule les coefficients de Fourier spatiaux f (0) des fonctions f (x, y), où la valeur I) indique que l'objet tri-dimensionnel est considéré à l'instant initial avant l'application du processus de diffusion. Les fréquences spatiales fondamentales utilisées pour le calcul des coefficients de Fourier sont les mêmes que la phase de calibration. ]0174] A l'étape 560, on détermine à partir de la loi de morphing latérale, les contours des tranches e=1,..,L au terme du processus de diffusion. Plus précisément, on calcule les coefficients l':,'(T) des harmoniques de rang vim2+3.7'- > j au moyen de (32) c'est à dire : 101751 IMath.521 (40) avec mp,,,(T) ( exp_â T = où ) ont été obtenus dans la phase de ca- libration. 101761 De manière similaire, on calcule les coefficients f tro,, Cr) et (T ) par : 101771 IMath.531 0 (41-1) ]0178] ]Math.54] 41-2 où les coefficients de contraction ,i7»:), e=1,..,L, ont été obtenus dans la phase de calibration. [0179] A l'étape 565, on effectue une transformée de Fourier spatiale inverse des coef- ficients J '':''(r)pour calculer la courbe de niveau f (2ç y): [0180] 21 [Math.55] 442) [0181] Enfin, à l'étape 570, on reconstruit la forme de l'objet-tridimensionnel à partir de courbes de niveau des tranches C ;, e = L, obtenues à l'étape 565, ces tranches étant désormais situées respectivement aux hauteurs (T)2, «=1,.., L. [0182] Les Figs. 6A-6B représentent l'ordinogramme d'une méthode de prédiction de la forme d'un objet tri-dimensionnel ayant subi un processus de diffusion, selon une variante du second mode de réalisation. [0183] Dans cette variante, on ne suppose plus que les tranches de l'objet tri-dimensionnel sont indépendantes et l'on effectue la calibration sur l'ensemble des tranches de chaque échantillon. Comme précédemment, la phase de calibration (Fig. 6A) est suivie de la phase de prédiction proprement dite (Fig. 6B). [0184] Les étapes 610-643 sont respectivement identiques aux étapes 541-543 du second mode de réalisation et leur description sera donc omise [0185] A l'étape 645, on estime les paramètres in, y(') vérifiant pour chacun des échantillons: [0186] [Math.561 ) (0) [0187] [Math.57] t) 17,,,,(C 3-2) [0188] [Math.58] ) - e.m)(A, )11 (43-3) où rn~a (0) (resp. F',,' (n) est le vecteur de taille L dont les éléments sont .f;:;1 (0) (resp. .f;;,':`,), (in) obtenus dans la phase de calibration et A.',' est la matrice symétrique tri- diagonale de taille L x L dont les termes de la diagonale principale sont égaux à umm, :4( = L, et dont les termes de la diagonale inférieure ainsi que de la diagonale supérieure sont égaux à C',,', où C'.' est un coefficient de couplage pour l'harmonique spatiale d'indices mm. entre deux tranches successives de l'échantillon. En pratique, la résolution du système d'équations (43-1), (43-2), (43-3) se fait par ajustement de courbes en ne considérant que les premiers du développent en série des exponentielles de matrices. [0189] On obtient ainsi une estimation de ces paramètres :,r), fiv), y [0190] La phase de calibration est suivie d'une phase de prédiction proprement dite, comprenant les étapes 650-670 en Fig. 6B. [0191] La phase de prédiction consiste à prédire, à partir de la forme d'un objet tri- dimensionnel à un instant t = 0, la forme de ce dernier après application d'un processus de diffusion pendant une durée T [0192] Les étapes 650-657 sont respectivement identiques aux étapes 550-557 du second mode de réalisation. [0193] En bref, on calcule le rayon équivalent, ru, de la tranche ( e =1)à la hase de l'objet et l'on détermine, à partir de la loi de morphing verticale, la hauteur de l'objet au terme du processus de diffusion selon (18). [0194] On calcule ensuite les coefficients de Fourier spatiaux J '',2' (0) des tranches C., aux hauteurs normalisées fe, L =1,.., L. On obtient ainsi une pluralité Ci vecteurs de taille L F ''',' (0), où Q est le nombre d'harmoniques pris en considération dans la transformée de Fourier. [0195] A l'étape 660, on détermine, à partir de la loi de morphing latérale, les contours des tranches L =1,.., L au ternie du processus de diffusion. Plus précisément, on calcule les coefficients ) des harmoniques de rang i',2+ > 1 au moyen de (36) c'est à dire : [0196] [Math.59] ( ) )1 (44) où na) est le vecteur des coefficients de Fourier spatiaux des tranches C; aux hauteurs normalisées fe, ; [0197] et pour les fondamentaux , à partir de (37-1) et (37-2) : [0198] [Math.60] (45-i) [0199] [Math.61] (45-2) [0200] A l'étape 665, on effectue une transformée de Fourier spatiale inverse des coef- ficients .1" (r) obtenus à l'étape précédente. [0201] Enfin, à l'étape 670, on reconstruit la forme de l'objet-tridimensionnel à partir de courbes de niveau des tranches C obtenues à l'étape 665, ces tranches étant désormais situées respectivement aux hauteurs h (T)2 = [0202] La méthode de prédiction de la forme d'un objet tri-dimensionnel, selon les premier et second modes de réalisation ainsi que leurs variantes exposées ci-dessus, permet de réduire considérablement les campagnes de mesure car la modélisation des lois de 23 morphing verticale et latérale ne nécessite qu'un nombre d'échantillons relativement faible. [0203] La Fig. 7 représente l'ordinogramme d'une méthode de fabrication d'un objet tri- dimensionnel de matériau donné et de forme de consigne prédéterminée. [0204] La méthode de fabrication comporte une phase de calibration, une phase de prédiction, une phase d'optimisation et une phase de fabrication proprement dite. [0205] Dans la phase de calibration, 710, on construit à partir de mesures d'échantillons avant et après un processus de diffusion de durée T, une base de données contenant les paramètres :i.''),â('),P") pour une pluralité de hauteurs normalisées =1,.., L. Cette phase de calibration peut être celle du premier mode de réalisation ou du second mode de réalisation, ou encore de leurs variantes. Le cas échéant, la base de données peut contenir ces paramètres pour une pluralité de valeurs de durée. [0206] Une fois la hase construite, on teste en 720 une pluralité de formes initiales d'un objet tri-dimensionnel et l'on applique à chacune de ces formes, la méthode de prédiction décrite dans les étapes 350-370 ou 450-470 ou 550-570 ou 650-670 selon le mode de réalisation choisi, pour prédire les formes au bout du processus de diffusion (le cas échéant en prenant en compte différentes durées de diffusion), dites formes finales. [0207] Dans une phase d'optimisation, 730, on détermine la forme finale, dite optimale, qui se rapproche le plus de la forme de consigne au sens d'une distance prédéterminée. [0208] Le cas échéant, on peut sélectionner une pluralité de formes finales se rapprochant le plus de la forme de consigne et interpoler les formes initiales correspondantes pour obtenir une nouvelle forme initiale à tester en 720. On peut ainsi, par itérations successives, obtenir une forme initiale optimale telle que la distance entre la forme finale correspondante et la forme de consigne soit alors inférieure à un seuil prédéterminé. [0209] Quelle que soit la manière dont la forme finale optimale est obtenue, clans la phase de fabrication, 740, on réalise un objet tri-dimensionnel ayant la forme initiale correspondante, par exemple par lithographie 3D ou fabrication additive, et l'on applique le processus de diffusion pendant la durée stockée dans la base de données. Au terme du processus, on obtient un objet tri-dimensionnel dont la forme est identique à ou proche de la forme de consigne. 24 [Revendication 1]