Fourier solvers have become efficient tools to establish structure–property relations in heterogeneous materials. Introduced as an alternative to the finite element (FE) method, they are based on fixed-point solutions of the Lippmann–Schwinger type integral equation. Their computational efficiency results from handling the kernel of this equation by the fast Fourier transform (FFT). However, the kernel is derived from an auxiliary homogeneous linear problem, which renders the extension of FFT-based schemes to nonlinear problems conceptually difficult. This paper aims to establish a link between FE-based and FFT-based methods in order to develop a solver applicable to general history-dependent and time-dependent material models. For this purpose, we follow the standard steps of the FE method, starting from the weak form, proceeding to the Galerkin discretization and the numerical quadrature, up to the solution of nonlinear equilibrium equations by an iterative Newton–Krylov solver. No auxiliary linear problem is thus needed. By analyzing a two-phase laminate with nonlinear elastic, elastoplastic, and viscoplastic phases and by elastoplastic simulations of a dual-phase steel microstructure, we demonstrate that the solver exhibits robust convergence. These results are achieved by re-using the nonlinear FE technology, with the potential of further extensions beyond small-strain inelasticity considered in this paper.
|Tijdschrift||International Journal for Numerical Methods in Engineering|
|Nummer van het tijdschrift||10|
|Status||Gepubliceerd - 7 sep 2017|