L'optimisation topologique est une méthode mathématique (et logicielle) qui permet de trouver la répartition de matière optimale dans un volume donné soumis à des contraintes. Elle se distingue notamment de l'optimisation de forme qui ne fait varier que la frontière de la pièce, qu'il faut avoir dessinée au départ[1].
Les méthodes d'optimisation topologiques connues à ce jour permettent d'optimiser la résistance élastico-mécanique, la conductivité thermique ou certains problèmes d'écoulement fluide.
Histoire
Cette méthode issue des mathématiques a été clairement définie, expliquée et rendue exploitable pour la mécanique dans les années 2000, notamment avec l'article fondateur de Ole Sigmund[2].
Des logiciels d’optimisation topologique de plus en plus sophistiqués permettent aux ingénieurs d'économiser la matière possible d’un objet tout en conservant ou améliorant sa solidité ou sa souplesse (au besoin) et en tenant compte des contraintes qui s’exerceront sur lui, travail autrefois basé sur l’intuition, la méthode des essais et des erreurs ou le génie des créateurs ou des ingénieurs de fabrication[3].
Un exemple très simple est la réduction optimisée du nombre de rayons d’une roue de bicyclette[3]. Jusqu’ici seules des formes simples étaient concernées, car ces logiciels sont très gourmands en calcul ou étaient rapidement limités par la complexité du travail demandé[3].
En , dans la revue Nature[4], des chercheurs d’une université danoise présentent une méthode permettant de faire ce travail pour des objets de grande taille, en améliorant la résolution possible. Une image 2D est composée de pixels alors qu’une image 3D est composée de voxels. Jusqu’à il y a peu[Quand ?], la résolution des modèles 3D optimisés était limitée à 5 millions de voxels, mais un nouveau programme optimise des objets jusqu'à 1 milliard de voxels. Ceci permet par exemple de modéliser et concevoir une aile de Boeing 777 en l’optimisant : elle se trouve allégée de 5 % et renforcée de l’intérieur par des longerons courbes et des nervures diagonales, et non plus par une structure en forme de grille. L'économie attendue est de 200 tonnes de kérosène/an[4]. 5 jours de calculs par un superordinateur ont été nécessaires, et ce design (qui évoque l’intérieur de certains os ou l'intérieur de parties d’exosquelettes d’insectes) est actuellement « infabricable ». Cependant, les progrès de l'impression 3D pourraient bientôt le mettre à notre portée[3],[4].
Formalisme mathématique
Énoncé
Le problème général d'optimisation topologique est posé comme
où
La variable de conception (design) est la densité du matériau. 1 signifie présence de matériau et 0 signifie absence de matériau. La variable d'état est la solution à une équation différentielle qui modélise le phénomène en jeu, selon le design (généralement les équations d'élasticité);
est la fonctionnelle coût. Elle représente une quantité d'intérêt qu'on veut minimiser ou maximiser (généralement la flexibilité ou la rigidité);
est l'espace de conception, c'est-à-dire l'ensemble des points qui sont sujets à l'optimisation;
Les contraintes permettent d'imposer des particularités à la forme finale. Par exemple, impose une contrainte sur la quantité maximale de matériau à avoir.
Interpolation des propriétés matérielles
Résoudre ce problème vise à déterminer une répartition optimale du matériau dans l'espace de conception . Le module de Young est la propriété qui quantifie la capacité d'un matériau à se déformer. Cette caractéristique doit être interpolée selon la densité (les endroits où il y a présence de matériau sont rigides, contrairement à ceux où il y a du vide). Par exemple, l'approche SIMP (Solid Isotropic Material with Penalization) utilise une loi de puissance pour déterminer le module de Young selon la densité en un point:
Le paramètre (typiquement[2]) pénalise les valeurs entre 0 et 1 pour obtenir un résultat binaire. est pris très près de 0 pour représenter le vide, mais non nul pour éviter les termes nuls dans la matrice de résolution de la méthode des éléments finis.
D'autres approches sont utilisées[5], comme la méthode RAMP, qui passe par une interpolation rationnelle.
Caractérisation de l'optimalité
L'annulation du gradient en un point est une condition nécessaire à l'optimalité dans un problème sans contrainte. La dérivée au sens de Gateaux permet de calculer le gradient d'une fonctionnelle dans une direction :
Une solution analytique à ce problème est difficile à obtenir. Les méthodes numériques d'optimisation comme la descente du gradient utilisent pour obtenir des designs presque optimaux. Les contraintes peuvent être prises en compte dans ces algorithmes en utilisant diverses techniques de pénalisation.
Mise en œuvre
En mécanique, la résolution d'un problème d'optimisation topologique passe par la modélisation de la pièce, ou de l'ensemble de pièces, à optimiser par la méthode des éléments finis. Une méthode classique d'optimisation topologique consiste alors à considérer en tout point du volume d'optimisation une densité de matière variant entre 0 et 1. D'autres méthodes considèrent l'orientation locale du matériau (pour les matériaux non isotropes) ou même d'autres caractéristiques. Optimiser consiste généralement dans ces méthodes à minimiser l'énergie de déformation de la structure, ce qui revient à peu près à trouver la structure la plus rigide possible. On peut ainsi soit fixer la quantité de matière mise en œuvre pour mettre en avant des formes optimales, pour orienter un design et une optimisation faite par ailleurs, soit chercher directement à définir une forme minimisant la matière à mettre en œuvre pour alléger au maximum la structure, en respectant une contrainte à ne pas dépasser. En pratique, pour un résultat exploitable (car c'est un nuage de densité de matière qui est généré a priori et non pas un solide avec une frontière bien définie) il est ajouté des méthodes de pénalisation, de filtrage[6], et de seuillage [7], en particulier pour imposer les contraintes géométriques spécifiques liées au processus de fabrication (symétries, autorisation de volume creux ou non, …, plan de joint)[8].
Les principales étapes et difficultés à franchir sont généralement les suivantes :
Espace réellement disponible : il est souvent bien plus grand que la pièce éventuellement existante, et peut être encore agrandi en remettant à plat la fonction réellement à remplir et les contraintes environnantes à cette pièce, ou de l'ensemble de pièces à reconcevoir. Il ne faut pas oublier les zones où la matière est imposée ou interdite (pour des raisons fonctionnelles ou esthétiques).
Liaisons mécaniques avec l'environnement : il faut bien remettre à plat les liaisons possibles avec les pièces voisines, car il y a souvent bien plus de liberté pour les zones de fixations que celles envisagées a priori. Il n'est parfois pas évident de distinguer quelles zones bloquer ou quelles zones sont chargées par des forces, le plus pragmatique étant alors d'imaginer comment pourrait être testée la pièce sur un banc d'essai, avec des liaisons fixes et des vérins, par exemple.
Efforts mécaniques subis : il faut bien prendre en compte tous les chargements mécaniques vus par la pièce, au-delà de la fonction principale, c.-à-d. les efforts liés à des étapes de fabrication (notamment l'usinage), les efforts liés à la manipulation de la pièce (montage/démontage de la pièce, transport), les efforts accidentels (chocs), par exemple.
Symétries et conditions de fabrication (cela étant de mieux en mieux pris en compte par les logiciels de calcul).
Lancer le calcul d'optimisation topologique : il faut adapter la finesse du maillage à la précision spatiale recherchée et aux moyens informatiques dont on dispose ; les calculs peuvent être longs, on cherche donc à faire des premiers calculs à l'échelle de quelques minutes, pour les affiner ensuite. Il faut aussi bien vérifier comment sont pris en compte les différents cas de chargement par l'algorithme. En effet si l'on cherche juste la structure la plus rigide possible pour une masse donnée, les énergies des différents chargements sont simplement sommées, il convient alors de les pondérer entre eux, éventuellement. Par contre, si l'objectif est d'obtenir la pièce la plus légère possible qui ne casse pas, nul besoin de pondération.
Analyse du résultat : Pour montrer une pièce facilement compréhensible (avec du vide et du plein bien définis), le résultat est généralement filtré par les logiciels pour l'affichage (par ex. le plein correspond aux zones de densité matière supérieures à 50%, sinon c'est du vide). Il faut donc bien prendre en compte qu'en général c'est une matière plus ou moins dense/poreuse qui est vraiment considérée par l'algorithme, et que les zones éventuelles de matière non reliée au reste sont tout à fait possibles, à l'affichage, car elles sont liées au reste par de la matière à faible densité, non affichée. Le résultat est donc à intreprêter pour définir une pièce faite de vide et de plein, au plus près de ce que l'algorithme propose. Il existe des paramètres (parfois cachés) pour explorer en détail ces subtilités : seuil d'affichage de la matière (par défaut 50% en général), pénalisation (paramètre limitant les zones de densités autour de 50%, mais pouvant dégrader la convergence des algorithmes), filtrage/lissage (filtre permettant d'éliminer des détails jugés trop petits), et bien sûr la finesse du maillage (qui permet de faire apparaître des détails plus ou moins fins). On se rend fréquemment compte à ce stade que la forme obtenue est absurde, généralement à la suite de l'oubli d'une contrainte majeure, ou car le problème était mal posé (par exemple s'il n'y a pas assez de liaisons au bâti pour maintenir la pièce, ou parce que des blocages ou des efforts ont été appliqués sur une zone où la matière est interdite).
Dessin et vérification : une fois l'interprétation des résultats consolidée, la pièce peut être dessinée, aussi proche que possible de la topologie obtenue (nombre de barres/plaques, orientation, épaisseurs relatives), mais éventuellement plus agréable à l'œil, car les formes dites «organiques» obtenues par optimisation topologique ne conviennent pas toujours. C'est pourquoi on impose parfois une peau extérieure à la pièce (la partie visible) en limitant l'optimisation topologique uniquement à l'intérieur de la pièce à alléger (partie invisible). Il est optimum si possible d'utiliser des structures lattices (c.-à-d. un réseau serré de poutres ou de parois, comme des mousses), afin de mettre de la matière de densité intermédiaire là où le calcul la fait apparaître (c'est ainsi que la nature exploite des densités de bois différentes dans l'arbre, ou des densités d'os variables chez les vertébrés.
↑ a et b(en) Sigmund O., « A 99 line topology optimization code written in Matlab », Struct Multidisc Optim 21, 120–127, Springer-Verlag 2001, (lire en ligne)
↑(en) Joshua D. Deaton et Ramana V. Grandhi, « A survey of structural and multidisciplinary continuum topology optimization: post 2000 », Structural and Multidisciplinary Optimization, vol. 49, no 1, , p. 1–38 (ISSN1615-1488, DOI10.1007/s00158-013-0956-z, lire en ligne, consulté le )
↑Bourdin B (2001) Filters in topology optimization. Int. J. Numer. Methods Eng. 50, 2143–2158
Chin T & Kennedy G (2016) Large-scale compliance-minimization and buckling topology optimization of the undeformed Common Research Model wing. In 57th AIAA/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference 0939 (AIAA)|résumé
Evgrafov, A., Rupp, C. J., Maute, K. & Dunn, M. L. (2008) Large-scale parallel topology optimization using a dual-primal substructuring solver. Struct. Multidiscipl. Optim. 36, 329–345 |résumé
Groen J.P & Sigmund O. (2017) Homogenization-based topology optimization for high-resolution manufacturable microstructures. Int. J. Num. Methods Eng. https://doi.org/10.1002/nme.5575%7Crésumé
Nomura T, Sato K, Taguchi K, Kashiwa T & Nishiwaki S (2007) Structural topology optimization for the design of broadband dielectric resonator antennas using the finite difference time domain technique. Int. J. Numer. Methods Eng. 71, 1261–1296|résumé
Rao J, Kiran S, Kamesh J, Padmanabhan M & Chandra S (2009) Topology optimization of aircraft wing. J. Aerosp. Sci. Technol. 61, 402–414
Saad Y (2003) Iterative Methods for Sparse Linear Systems Ch. 9.4.1 (SIAM)
Questa voce sull'argomento università del Cile è solo un abbozzo. Contribuisci a migliorarla secondo le convenzioni di Wikipedia. Universidad de los Andes UbicazioneStato Cile CittàSantiago del Cile Dati generaliNome latinoVniversitas Stvdiorvm Andinensis Fondazione1989 Tipoprivata, cattolica RettoreJosé Antonio Guzmán Cruzat Mappa di localizzazione Sito web Modifica dati su Wikidata · Manuale L'Universidad de los Andes (con acronimo UANDES) è un'università privat...
Ector County, TexasLokasi di negara bagian TexasLokasi negara bagian Texas di Amerika SerikatDidirikan1893SeatOdessaWilayah • Keseluruhan902 sq mi (2.336 km2) • Daratan901 sq mi (2.334 km2) • Perairan1 sq mi (3 km2), 0.07%Populasi • (2000)121.123 • Kepadatan135/sq mi (52/km²)Situs webwww.co.ector.tx.us County Ector adalah county yang terletak di negara bagian Texas, Amerika Serikat. ...
Snap mor di pantai Samber, Biak Suasana Snap mor di pantai Samber-Binyeri, Distrik Yendidori, Kabupaten Biak Numfor Ikan yang ditangkap dengan tombak atau kalawai saat Snap mor Snap mor adalah tradisi menyebar jaring atau menangkap ikan di Biak, Papua. Kegiatan snap mor bisa dilakukan pada saat air surut, yang ditandai oleh angin timur dan curah hujan yang dominan, serta dilakukan di daerah yang dangkal. Snap mor berasal dari bahasa Biak, yaitu Snap dan mor. Snap adalah koral atau batu kecil ...
قرية ليونز الإحداثيات 43°03′48″N 76°59′27″W / 43.0633°N 76.9908°W / 43.0633; -76.9908 [1] تقسيم إداري البلد الولايات المتحدة[2] التقسيم الأعلى لايونز خصائص جغرافية المساحة 12.049666 كيلومتر مربع10.750437 كيلومتر مربع (1 أبريل 2010) ارتفاع 125 متر عدد السكا...
American baseball player (born 1982) Baseball player Jered WeaverWeaver with the Angels in 2012PitcherBorn: (1982-10-04) October 4, 1982 (age 41)Northridge, California, U.S.Batted: RightThrew: RightMLB debutMay 27, 2006, for the Los Angeles Angels of AnaheimLast MLB appearanceMay 19, 2017, for the San Diego PadresMLB statisticsWin–loss record150–98Earned run average3.63Strikeouts1,621 Teams Los Angeles Angels of Anaheim / Los Angeles Angels (2006–2016)...
I Want to Know What Love IsLagu oleh Foreignerdari album Agent ProvocateurSisi-BStreet Thunder (Marathon Theme)Dirilis13 November 1984Format7GenrePop rock, soft rockDurasi5:00 (album/7 single) 6:23 (12 single)4:15 (radio edit/berbagai versi sunting)LabelAtlanticPenciptaMick JonesProduserMick Jones, Alex SadkinSampel Musik I Want to Know What Love Is noiconnoicon I Want to Know What Love Is adalah sebuah power ballad karya band rock Inggris-Amerika Foreigner. Lagu tersebut dirilis pada Novembe...
British-bred Thoroughbred racehorse La SagesseSireWisdomGrandsireBlinkhoolieDamSt MaryDamsireHermitSexMareFoaled1892[1]CountryUnited KingdomColourBayBreederGeorge Alexander BairdOwnerSir James Percy Miller, 2nd BaronetTrainerMartin GurryRecord27: 8-1-2Major winsNewmarket Breeders' Plate (1894)Fitzwilliam Stakes (1894)Oaks Stakes (1895)Derby Cup (1896) La Sagesse (1892 – 1909) was a British Thoroughbred racehorse and broodmare. She was highly tried as a juvenile in 1894, winning ...
Cet article est une ébauche concernant la république démocratique du Congo. Vous pouvez partager vos connaissances en l’améliorant (comment ?) selon les recommandations des projets correspondants. Pour les articles homonymes, voir Lwano. Lwano Caractéristiques Bassin collecteur Congo Cours Source au nord-est du centre-ville de Lubumbashi Confluence Rwashi Géographie Pays traversés République démocratique du Congo Régions traversées Haut-Katanga modifier La Lwano, auss...
GondangrejoKecamatanPeta lokasi Kecamatan GondangrejoNegara IndonesiaProvinsiJawa TengahKabupatenKaranganyarPemerintahan • Camat-Populasi • Total60,251 (2.003) jiwaKode Kemendagri33.13.13 Kode BPS3313130 Luas56,8 km²Desa/kelurahan13 Gondangrejo (bahasa Jawa: ꦒꦺꦴꦤ꧀ꦝꦤꦔꦿꦺꦗ, translit. Gondhangreja) adalah sebuah kecamatan di Kabupaten Karanganyar, Jawa Tengah. Batas Wilayah Utara Kecamatan Kalijambe, Kabupaten Sragen dan Kecamatan Plupu...
South Korean badminton player In this Korean name, the family name is Kim. Badminton playerKim Min-jungPersonal informationCountrySouth KoreaBorn (1986-07-29) 29 July 1986 (age 37)Jeju, South KoreaHeight1.68 m (5 ft 6 in)Weight63 kg (139 lb)HandednessRightWomen's & mixed doublesHighest ranking2 (WD)11 (XD) Medal record Women's badminton Representing South Korea World Cup 2005 Yiyang Women's doubles Sudirman Cup 2009 Guangzhou Mixed team 2011 Qingdao Mix...
Foto rumah Streisand di Malibu yang menjadi asal usul nama efek ini Efek Streisand adalah fenomena ketika upaya untuk menyembunyikan, menghapus, atau menyensor informasi malah membuat informasi tersebut tersebar lebih luas, biasanya dibantu oleh Internet. Efek Streisand merupakan contoh reaktansi psikologis: ketika masyarakat sadar bahwa ada informasi yang disembunyikan, mereka akan berusaha mengaksesnya dan menyebarkannya.[1] Nama efek ini diambil dari nama penyanyi asal Amerika Seri...
Stimmzettel beim Referendum Am 22. Januar 2012 fand in Kroatien ein obligatorisches Referendum über die Frage des anstehenden EU-Beitritts Kroatiens statt. Die große Mehrheit der Abstimmenden votierten für den EU-Beitritt. Die Wahlbeteiligung lag bei 43,51 %. Inhaltsverzeichnis 1 Verhältnis Kroatiens zur Europäischen Union 2 Das Referendum 3 Ergebnisse 4 Folgen 5 Einzelnachweise Verhältnis Kroatiens zur Europäischen Union Kroatien war nach dem Zerfall Jugoslawiens 1991 ein unabhä...
New Zealand mountaineer (1919–2008) Sir Edmund HillaryKG ONZ KBEHillary in c. 1953BornEdmund Percival Hillary(1919-07-20)20 July 1919Auckland, New ZealandDied11 January 2008(2008-01-11) (aged 88)Auckland, New ZealandKnown forWith Tenzing Norgay, first to reach summit of Mount EverestSpouses Louise Mary Rose (m. 1953; died 1975) June Mulgrew (m. 1989) Children3, including PeterMilitary ...
Claire CarletonCarleton di film tahun 1946, A Close Call for Boston BlackieLahir(1913-09-28)28 September 1913New York City, A.S.Meninggal11 Desember 1979(1979-12-11) (umur 66)Northridge, California, A.S.MakamForest Lawn Memorial Park (Hollywood Hills)PekerjaanAktrisTahun aktif1932–69Suami/istriFred E. Sherman(m. 19??; wafat 1969) Claire Carleton (28 September 1913 – 11 Desember 1979) adalah seorang aktris Amerika yang karirnya berlangsung selama empat dekade dari ...
Bài này viết về sân vận động được khánh thành vào năm 1923. Đối với sân vận động mới được thay thế cho sân này, xem Sân vận động Wembley. Đối với đấu trường trong nhà, xem Wembley Arena. Sân vận động WembleyTòa Tháp đôi của Sân vận động WembleyTên cũSân vận động Đế quốcSân vận động Triển lãm Đế quốc AnhVị tríWembley, Luân Đôn, AnhTọa độ51°33′20″B 0°16′47″T / ...
Friedrich SchleiermacherFriedrich SchleiermacherSinh(1768-11-21)21 tháng 11, 1768Breslau, Silesia, PhổMất12 tháng 2, 1834(1834-02-12) (65 tuổi)Berlin, Brandenburg, PhổTrường lớpĐại học Halle (1787–90)[1]Thời kỳThế kỷ 18/19VùngTriết học phương TâyTrường pháiGerman Idealism[2]Jena Romanticism[3]Berlin Romanticism[4]Romantic hermeneutics[5]Tổ chứcĐại học Halle (1804–07)Đại học Berlin (1810–18...