A non-local approach for Reissner–Mindlin shell elements in dynamic simulations: Application with a Gurson model

Abstract

Numerical prediction of ductile tearing during a car crash is necessary to ensure the structural integrity of the vehicle. Recently, encouraging results were obtained by Davaze et al., (2021) with a non-local damage approach, compatible with dynamic explicit simulations and preserving parallel computation: an implicit gradient model enabled to reproduce experimental results with DP450 steel sheets with 3D solid elements. However, in the automotive industry, thin sheets of metals are mostly modeled with shell elements for computation time efficiency. Shell elements have a more complex formulation: their kinematic description is enriched accounting for displacements but also rotations, leading to two different strain contributions: membrane strain and bending strain. The transverse strain is computed independently in a post-processing step from local internal variables with a zero normal stress assumption. This formulation is known to lead to strain localization after the onset of necking, while this is not the case with 3D solid elements. The problem is even more challenging when considering shells in association with a softening material accounting for damage as this also leads to strong mesh dependency, as with 3D solid elements. This complex problem is solved using non-local formulations for the material behavior. The proposed solution accounts for both in-plane and in-thickness damage gradients. It consists of a two-step regularization procedure introducing nested wire elements in the shell thickness. Using non-local variables to evaluate the transverse strain enables to obtain a continuous spatial distribution for strains, in agreement with the 3D solid case. An updated Lagrangian formulation is used to account for large strains. The performance of the strategy is demonstrated in three test cases, with a comparison with an experimental result for one of them. The results obtained following the proposed strategy lead are equivalent to the ones obtained with the non-local method proposed by Davaze et al., (2021) and are thus consistent with the experimental data.