FFT-based homogenization of hypoelastic plasticity at finite strains

• Published in: Computer methods in applied mechanics and engineering 2019-06, Vol.349, p.499-521
• Main Authors: Ma, Ran, Truster, Timothy J.
• Format: Article
• Language: English
• Subjects: Algorithms; Boundary conditions; Configurations; Constitutive models; Crystal plasticity; Crystal pulling; Fast Fourier transformations; FFT; Finite element method; Finite strain; Fourier transforms; Galerkin method; Homogenization; Hypoelastic; Mathematical analysis; Mathematical models; Mechanical analysis; Multiscale analysis; Newton methods; Nonlinear programming; Parallel processing; Plastic properties; Polynomials; Robustness (mathematics); Stiffness
• ISSN: 0045-7825; 1879-2138
• DOI: 10.1016/j.cma.2019.02.037

Fast Fourier Transform (FFT) based methods are becoming increasingly popular for modeling the texture evolution and local mechanical response of polycrystalline materials within the representative volume element (RVE). Originally, the FFT-based method was formulated through the Lippman–Schwinger (L–S) equation in terms of a homogeneous elastic reference medium. More recently, the Fourier–Galerkin method was proposed to discretize the RVE weak form independently of the reference medium using trigonometric polynomials. The Fourier–Galerkin method, albeit efficient and robust, has not been extended to solve crystal plasticity (CP) problems with spatial heterogeneity and fine resolution. Also, its algorithmic homogenization has not yet been investigated, which is essential for concurrent multiscale modeling. In this paper, a general framework is proposed that connects the FFT-based method and objective rate constitutive models, in particular a generalized implementation of crystal plasticity. Consistent linearization is achieved by pulling back the tangent stiffness from unrotated configuration to reference configuration, which improves the convergence behavior of the Newton–Krylov solver. Applying the inexact Newton method further improves the numerical efficiency. Also, the algorithmic homogenized tangent stiffness for the Fourier–Galerkin method is derived to enable mixed boundary conditions and concurrent multiscale modeling. Lastly, the Fourier–Galerkin method’s accuracy and efficiency for solving CP problems are studied in a parallelized environment, and significant speedup is observed versus the finite element method and L–S based FFT method. •general interface between FFT-based method and objective rate constitutive models.•homogenization of Fourier–Galerkin method.•efficient FFT finite strain crystal plasticity.•inexact Newton solver acceleration.•numerical efficiency compared with FE and other FFT methods.