An assessment of numerical techniques to find energy‐minimizing microstructures associated with nonconvex potentials

Abstract

Summary Microstructural patterns emerge ubiquitously during phase transformations, deformation twinning, or crystal plasticity. Challenges are the prediction of such microstructural patterns and the resulting effective material behavior. Mathematically, the experimentally observed patterns are energy‐minimizing sequences produced by an underlying non‐(quasi)convex strain energy. Therefore, identifying the microstructure and effective response is linked to finding the quasiconvex, relaxed energy. Due to its nonlocal nature, quasiconvexification has traditionally been limited to (semi‐)analytical techniques or has been dealt with by numerical techniques such as the finite element method (FEM). Both have been restricted to primarily simple material models. We here contrast three numerical techniques—FEM, a Fourier‐based spectral formulation, and a meshless maximum‐entropy (max‐ent) method. We demonstrate their performance by minimizing the energy of a representative volume element for both hyperelasticity and finite‐strain phase transformations. Unlike FEM, which fails to converge in most scenarios, the Fourier‐based spectral formulation (FFT) scheme captures microstructures of intriguingly high resolution, whereas max‐ent is superior at approximating the relaxed energy. None of the methods are capable of accurately predicting both microstructures and relaxed energy; yet, both FFT and max‐ent show significant advantages over FEM. Numerical errors are explained by the energy associated with microstructural interfaces in the numerical techniques compared here.

Data-driven topology optimization of spinodoid metamaterials with seamlessly tunable anisotropy

“…for spinodoids, ρ ∈ {0} ∪ [ρ min , 1] and θ 1 , θ 2 , θ 3 ∈ {0} ∪ [0, π/2]. Similar observations were reported and explained previously by, e.g., Shiye and Jiejiang (2016) and Kumar et al (2020b).…”

A review of nonlinear FFT-based computational homogenization methods

“…FFT-based computational techniques proved useful for studying the viscoelasticity of cement paste [310], mortar samples [311], cement [312] and concrete [313]. Furthermore, FFT-based techniques were used for studying explosive materials [314], secondary creep in a porous nuclear fuel [315], the thermal expansion of an energetic material [316], optical properties of deposit models for paint [317], dynamic recrystallization [318], to compute geodesics in two-dimensional media [182], fitting microstructure-property relationships [319], topology optimization [320] and for finding emerging microstructures associated to non-convex potentials [93,321].…”

Numerical approaches for investigating quasiconvexity in the context of Morrey's conjecture

“…One may conjecture that minimizing W + magic in a finite element space fails to find microstructure because each finite element coefficient influences only a very local part of the deformation function. Additionally, when the strain energy density lacks quasiconvexity (which we consider possible for W + magic ), FEM-methods generally fail [44] due to the non-uniqueness of the solution, i.e. the microstructures.…”