A surrogate model for computational homogenization of elastostatics at finite strain using high‐dimensional model representation‐based neural network

• Published in: International journal for numerical methods in engineering 2020-11, Vol.121 (21), p.4811-4842
• Main Authors: Minh Nguyen‐Thanh, Vien, Trong Khiem Nguyen, Lu, Rabczuk, Timon, Zhuang, Xiaoying
• Format: Article
• Language: English
• Subjects: Boundary value problems; Collocation methods; computational homogenization; Computer architecture; data‐driven; Elasticity; Elastostatics; Equilibrium equations; FFT‐based methods; Finite element method; Flux density; Homogenization; Iterative methods; Mathematical analysis; Mathematical models; Neural networks; nonlinear elasticity; Numerical stability; Representations; Solvers; Strain
• ISSN: 0029-5981; 1097-0207
• DOI: 10.1002/nme.6493

Summary We propose a surrogate model for two‐scale computational homogenization of elastostatics at finite strains. The macroscopic constitutive law is made numerically available via an explicit formulation of the associated macroenergy density. This energy density is constructed by using a neural network architecture that mimics a high‐dimensional model representation. The database for training this network is assembled through solving a set of microscopic boundary value problems with the prescribed macroscopic deformation gradients (input data) and subsequently retrieving the corresponding averaged energies (output data). Therefore, the two‐scale computational procedure for nonlinear elasticity can be broken down into two solvers for microscopic and macroscopic equilibrium equations that work separately in two stages, called the offline and online stages. The finite element method is employed to solve the equilibrium equation at the macroscale. As for microscopic problems, an FFT‐based collocation method is applied in tandem with the Newton‐Raphson iteration and the conjugate‐gradient method. Particularly, we solve the microscopic equilibrium equation in the Lippmann‐Schwinger form without resorting to the reference medium. In this manner, the fixed‐point iteration that might require quite strict numerical stability conditions in the nonlinear regime is avoided. In addition, we derive the projection operator used in the FFT‐based method for homogenization of elasticity at finite strain.