Multiscale design of nonlinear materials using a Eulerian shape optimization scheme

• Published in: International journal for numerical methods in engineering Vol. 122; no. 12; pp. 2981 - 3014
• Main Authors: Najafi, Ahmad R., Safdari, Masoud, Tortorelli, Daniel A., Geubelle, Philippe H.
• Format: Journal Article
• Language: English
• Published: Hoboken, USA John Wiley & Sons, Inc 30.06.2021; Wiley Subscription Services, Inc; Wiley
• Subjects: Analytical sensitivity analysis; Cost function; Derivatives; ENGINEERING; Fixed grid; Gradient-based shape optimization; Mechanical and civil engineering; Microstructure; Nonlinear GFEM; NURBS; Optimization; Particulate composites; Sensitivity analysis; Shape functions; Shape optimization
• ISSN: 0029-5981; 1097-0207
• DOI: 10.1002/nme.6650

Motivated by recent advances in manufacturing, the design of materials is the focal point of interest in the material research community. One of the critical challenges in this field is finding optimal material microstructure for a desired macroscopic response. This work presents a computational method for the mesoscale‐level design of particulate composites for an optimal macroscale‐level response. The method relies on a custom shape optimization scheme to find the extrema of a nonlinear cost function subject to a set of constraints. Three key “modules” constitute the method: multiscale modeling, sensitivity analysis, and optimization. Multiscale modeling relies on a classical homogenization method and a nonlinear NURBS‐based generalized finite element scheme to efficiently and accurately compute the structural response of particulate composites using a nonconformal discretization. A three‐parameter isotropic damage law is used to model microstructure‐level failure. An analytical sensitivity method is developed to compute the derivatives of the cost/constraint functions with respect to the design variables that control the microstructure's geometry. The derivation uncovers subtle but essential new terms contributing to the sensitivity of finite element shape functions and their spatial derivatives. Several structural problems are solved to demonstrate the applicability, performance, and accuracy of the method for the design of particulate composites with a desired macroscopic nonlinear stress‐strain response.