DBFFT: A displacement based FFT approach for non-linear homogenization of the mechanical behavior

Abstract Most of the FFT methods available for homogenization of the mechanical response use the strain/deformation gradient as unknown, imposing their compatibility using Green’s functions or projection operators. This implies the allocation of redundant information and, when the method is based in solving a linear equation, the rank-deficiency of the resulting system. In this work we propose a fast, robust and memory-efficient FFT homogenization framework in which the displacement field on the Fourier space is the unknown: the displacement based FFT (DBFFT) algorithm. The framework allows any general non-linear constitutive behavior for the phases and direct strain, stress and mixed control of the macroscopic load. In the linear case, the method results in a linear system defined in terms of linear operators in the Fourier space and that does not require a reference medium. The system has an associated full rank Hermitian matrix and can be solved using iterative Krylov solvers and allows the use of preconditioners. A preconditioner is proposed to improve the efficiency of the system resolution. Finally, some numerical examples including elastic, hyperelastic and viscoplastic materials are solved to check the accuracy and efficiency of the method. The computational cost reduction respect the Galerkin-FFT was around 30%.

[1]  Luc Dormieux,et al.  FFT-based methods for the mechanics of composites: A general variational framework , 2010 .

[2]  J. Segurado,et al.  Computational Homogenization of Polycrystals , 2018, 1804.02538.

[3]  M. Schneider,et al.  Efficient fixed point and Newton–Krylov solvers for FFT-based homogenization of elasticity at large deformations , 2014 .

[4]  L. Dormieux,et al.  Combining Galerkin approximation techniques with the principle of Hashin and Shtrikman to derive a new FFT-based numerical method for the homogenization of composites , 2012 .

[5]  V. G. Kouznetsova,et al.  Multi-scale computational homogenization: Trends and challenges , 2010, J. Comput. Appl. Math..

[6]  Hervé Moulinec,et al.  A numerical method for computing the overall response of nonlinear composites with complex microstructure , 1998, ArXiv.

[7]  M. Geers,et al.  A finite element perspective on nonlinear FFT‐based micromechanical simulations , 2016, 1601.05970.

[8]  Javier Segurado,et al.  On the accuracy of spectral solvers for micromechanics based fatigue modeling , 2018, Computational Mechanics.

[9]  J. Segurado,et al.  An algorithm for stress and mixed control in Galerkin‐based FFT homogenization , 2019, International Journal for Numerical Methods in Engineering.

[10]  Jaroslav Vondrejc,et al.  An FFT-based Galerkin method for homogenization of periodic media , 2013, Comput. Math. Appl..

[11]  F. Willot,et al.  Fourier‐based schemes with modified Green operator for computing the electrical response of heterogeneous media with accurate local fields , 2013, 1307.1015.

[12]  Hervé Moulinec,et al.  Comparison of three accelerated FFT‐based schemes for computing the mechanical response of composite materials , 2014 .

[13]  Jan Novák,et al.  Accelerating a FFT-based solver for numerical homogenization of periodic media by conjugate gradients , 2010, J. Comput. Phys..

[14]  M. Schneider,et al.  Computational homogenization of elasticity on a staggered grid , 2016 .

[15]  Graeme W. Milton,et al.  A fast numerical scheme for computing the response of composites using grid refinement , 1999 .

[16]  R. Lebensohn N-site modeling of a 3D viscoplastic polycrystal using Fast Fourier Transform , 2001 .

[17]  G. Bonnet,et al.  A polarization‐based FFT iterative scheme for computing the effective properties of elastic composites with arbitrary contrast , 2012 .

[18]  H. Moulinec,et al.  A fast numerical method for computing the linear and nonlinear mechanical properties of composites , 1994 .

[19]  M. Schneider,et al.  Mixed boundary conditions for FFT-based homogenization at finite strains , 2016 .

[20]  M. Geers,et al.  Finite strain FFT-based non-linear solvers made simple , 2016, 1603.08893.

[21]  A. Cruzado,et al.  Multiscale modeling of the mechanical behavior of IN718 superalloy based on micropillar compression and computational homogenization , 2015 .

[22]  Hervé Moulinec,et al.  A Computational Method Based on Augmented Lagrangians and Fast Fourier Transforms for Composites with High Contrast , 2000 .