Table of Contents
Fetching ...

Efficient Fine-Scale Simulation of Nonlinear Hyperelastic Lattice Structures

Clément Guillet, Thibaut Hirschler, Pierre Jolivet, Pablo Antolin, Robin Bouclier

Abstract

With the growing maturity of additive manufacturing, the fabrication of architected or lattice-based metamaterials has become a reality for industrial applications. These materials combine lightweight design with tailored mechanical properties, most of which exhibit pronounced nonlinear, especially large-deformation, behaviors. The main numerical challenge therefore lies in performing nonlinear simulations of such lattice structures, which may contain thousands of geometrically intricate unit cells, while lacking sufficient scale separation for multiscale homogenization schemes to be applicable straightforwardly. In this work, we propose a dedicated solver for the full volumetric fine-scale simulation of nonlinear hyperelastic lattice structures that drastically reduces both memory and computational costs. The key idea is to exploit the intrinsic self-similarity of the cells through a reduced-order modeling strategy applied within a domain-decomposition framework. At each Newton iteration, a limited set of principal cells is identified through a dedicated, weakly intrusive, EIM-like approach, allowing all local tangent operators to be expressed as linear combinations of a few principal ones. This enables fast and memory-efficient operator assembly, and then feeds an efficient inexact FETI-DP based preconditioner at the solution stage, resulting in a quasi matrix-free algorithm for the nonlinear analysis. Numerical experiments in two and three dimensions demonstrate significant computational gains, with runtime reductions from several hours to a few tens of minutes and memory savings by factors of about three, while maintaining full fine-scale accuracy. Notably, the proposed strategy enables the computation of problems involving thousands of cells (i.e., millions of degrees of freedom) within a few minutes on an off-the-shelf laptop.

Efficient Fine-Scale Simulation of Nonlinear Hyperelastic Lattice Structures

Abstract

With the growing maturity of additive manufacturing, the fabrication of architected or lattice-based metamaterials has become a reality for industrial applications. These materials combine lightweight design with tailored mechanical properties, most of which exhibit pronounced nonlinear, especially large-deformation, behaviors. The main numerical challenge therefore lies in performing nonlinear simulations of such lattice structures, which may contain thousands of geometrically intricate unit cells, while lacking sufficient scale separation for multiscale homogenization schemes to be applicable straightforwardly. In this work, we propose a dedicated solver for the full volumetric fine-scale simulation of nonlinear hyperelastic lattice structures that drastically reduces both memory and computational costs. The key idea is to exploit the intrinsic self-similarity of the cells through a reduced-order modeling strategy applied within a domain-decomposition framework. At each Newton iteration, a limited set of principal cells is identified through a dedicated, weakly intrusive, EIM-like approach, allowing all local tangent operators to be expressed as linear combinations of a few principal ones. This enables fast and memory-efficient operator assembly, and then feeds an efficient inexact FETI-DP based preconditioner at the solution stage, resulting in a quasi matrix-free algorithm for the nonlinear analysis. Numerical experiments in two and three dimensions demonstrate significant computational gains, with runtime reductions from several hours to a few tens of minutes and memory savings by factors of about three, while maintaining full fine-scale accuracy. Notably, the proposed strategy enables the computation of problems involving thousands of cells (i.e., millions of degrees of freedom) within a few minutes on an off-the-shelf laptop.
Paper Structure (29 sections, 60 equations, 12 figures, 3 tables, 1 algorithm)

This paper contains 29 sections, 60 equations, 12 figures, 3 tables, 1 algorithm.

Figures (12)

  • Figure 1: Brake pedal lattice structure, inspired from ntopo and defined by spline composition of a reference microstructure unit-cell and a macro-geometry mapping
  • Figure 2: RB strategy to compute local tangent operators.
  • Figure 3: Examples of identified principal cells for (top) a rectangular lattice and (bottom) a brake-pedal lattice, both based on a cross-truss unit cell. The reduced-basis tolerance is set to $\epsilon = 3 \times 10^{-4}$. Results are shown for both the linear (small-deformation) and nonlinear (hyperelastic large-deformation) cases. As expected, increasing the complexity of the macroscopic mapping or the deformation level leads to a larger number of principal cells. Conversely, as the number of cells increases, the likelihood of encountering mechanically similar ones also increases, thereby enhancing the computational savings achieved by the proposed strategy.
  • Figure 4: Non definiteness of a local remaining matrix for a deformation with buckling of bottom-hand unit-cell edge under the action of a vertical downward force applied to the right face of the structure. All matrices become definite when adding one primal DOF on each edge of the unit cell.
  • Figure 5: Different analysis-suitable models (obtained by applying $k$-refinement) of the unit cells (UC) considered. Characteristic mesh size is defined as $h = 1/n_{e}$, where $n_e$ is approximately the number of element per direction per patch of the underlying model.
  • ...and 7 more figures

Theorems & Definitions (7)

  • Remark 2.1
  • Definition 3.1: Finite-dimensional approximation space
  • Remark 3.1
  • Remark 3.2
  • Remark 3.3
  • Remark 3.4
  • Remark 4.1