A robust and scalable unfitted adaptive finite element framework for nonlinear solid mechanics

Santiago Badia,Manuel Caicedo,Alberto F. Martín,Javier Principe

DOI: https://doi.org/10.1016/j.cma.2021.114093

2021-07-25

Abstract:In this work, we bridge standard adaptive mesh refinement and coarsening on scalable octree background meshes and robust unfitted finite element formulations for the automatic and efficient solution of large-scale nonlinear solid mechanics problems posed on complex geometries, as an alternative to standard body-fitted formulations, unstructured mesh generation and graph partitioning strategies. We pay special attention to those aspects requiring a specialized treatment in the extension of the unfitted h-adaptive aggregated finite element method on parallel tree-based adaptive meshes, recently developed for linear scalar elliptic problems, to handle nonlinear problems in solid mechanics. In order to accurately and efficiently capture localized phenomena that frequently occur in nonlinear solid mechanics problems, we perform pseudo time-stepping in combination with h-adaptive dynamic mesh refinement and rebalancing driven by a-posteriori error estimators. The method is implemented considering both irreducible and mixed (u/p) formulations and thus it is able to robustly face problems involving incompressible materials. In the numerical experiments, both formulations are used to model the inelastic behavior of a wide range of compressible and incompressible materials. First, a selected set of benchmarks are reproduced as a verification step. Second, a set of experiments is presented with problems involving complex geometries. Among them, we model a cantilever beam problem with spherical hollows distributed in a Simple Cubic array. This test involves a discrete domain with up to 11.7M Degrees Of Freedom solved in less than two hours on 3072 cores of a parallel supercomputer.

Numerical Analysis,Mathematical Software

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is to provide a robust and scalable mesh - free adaptive finite - element framework for large - scale nonlinear solid mechanics problems. Specifically, it aims to address the following aspects:

1. **Computation on complex geometries**: Traditional methods require the generation of unstructured meshes conforming to geometric boundaries when dealing with complex geometries, which is time - consuming and difficult to automate. The method proposed in this paper embeds the computational domain in an easily - generated background mesh by using the unfitted finite - element method, thereby greatly reducing the requirements for geometric constraints.

2. **Capturing local phenomena**: Nonlinear solid mechanics problems often involve local phenomena (such as strain localization and fracture), which require high - precision meshes to capture. This paper introduces pseudo - time stepping combined with h - adaptive dynamic mesh refinement and re - balancing techniques, driven by a posteriori error estimates, to accurately and efficiently capture these local phenomena.

3. **Parallel computing and load balancing**: To achieve large - scale parallel computing, the method proposed in this paper utilizes a tree - based scalable background mesh (such as an octree) and solves the load - imbalance problem caused by localization through a dynamic load - balancing strategy.

4. **Ill - conditioned small - cut - element problem**: Small - cut elements in the embedded method can lead to discretization problems with poor condition numbers. For this reason, this paper adopts the Aggregated Finite - Element Method (AgFEM) and stabilizes the discretization system by adding additional constraint conditions.

5. **Treatment of history variables**: For nonlinear material behaviors that depend on strain history, this paper proposes a method of representing history variables as fields and selects a discontinuous Lagrange finite - element space to approximate these variables, thereby simplifying the variable transfer between different meshes and reducing memory requirements.

In summary, this paper is committed to developing a framework that can efficiently solve nonlinear solid mechanics problems on large - scale and complex geometries, with particular attention to how to improve computational efficiency and accuracy through advanced numerical techniques and parallel computing strategies.