Indentation with different materials

This example illustrates the use of adaptive meshing in indentation problems with different materials.

This page discusses:

Problem description
Adaptive meshing
Results and discussion
Input files
Figures

ProductsAbaqus/Explicit

Problem description

The finite element model for the problem is axisymmetric, as shown in Figure 1. The model consists of a rigid punch and a deformable blank. The blank, which is meshed with CAX4R elements, has a radius of 600 mm and a height of 300 mm. The punch is modeled as an analytical rigid surface. Coulomb friction is modeled between the punch and the blank with a friction coefficient of 0.1. Symmetry boundary conditions are defined at r=0 for the blank. The bottom of the blank is also fully constrained.

Several analyses are performed using the following material models for the blank: hyperelasticity, hyperelasticity with viscoelasticity, hyperfoam, Hill plasticity, Mises plasticity, rate-dependent Mises plasticity, Drucker-Prager plasticity, Drucker-Prager cap plasticity, crushable foam plasticity, and porous metal plasticity. The parameters and constants used for each material model can be found in the input files that are included with the Abaqus release. The punch is fully constrained except in the vertical direction, in which motion is prescribed such that the maximum indentation depth is 250 mm. The blank is indented dynamically when it is modeled with a rate-dependent material (hyperelasticity with viscoelasticity or rate-dependent Mises plasticity). A ramped velocity profile is prescribed such that the maximum velocity is 2000 mm/sec. The blank is indented quasi-statically when it is modeled with the remaining (rate-independent) material models. A smooth-step amplitude curve is used to specify the displacement of the punch and promote a quasi-static solution.

Adaptive meshing

A single adaptive mesh domain that incorporates the entire blank is used. The symmetry boundary conditions are defined as Lagrangian boundary regions (the default), and contact surfaces are defined as sliding boundary regions (the default). For some materials a large amount of localized deformation occurs immediately under the punch; therefore, nondefault values are used to specify the frequency of adaptive meshing and the number of mesh sweeps to be performed in each adaptive mesh increment. These values are chosen to solve the problem economically yet retain a good mesh throughout the simulation.

Results and discussion

With the exception of the hyperfoam material, an indentation of this depth cannot be simulated using a pure Lagrangian simulation. The crushable foam material can be indented a majority of this depth using a pure Lagrangian approach. For each of the materials continuous adaptive meshing maintains the quality of the mesh throughout the analysis, and contours of key variables appear to be reasonable. For the hyperfoam material comparisons of the results from the adaptive mesh simulation and a pure Lagrangian simulation can be made. Figure 2 and Figure 3 show contours of the y-component of nominal strain in the blank after punching using the adaptive mesh and pure Lagrangian approaches, respectively. The results are in good agreement.

Input files

ale_matverif_indentmises.inp: Mises plasticity.
ale_matverif_indenthyper.inp: Hyperelasticity.
ale_matverif_indentvishyper.inp: Hyperelasticity with viscoelasticity.
ale_matverif_indenthyperfoam.inp: Hyperfoam.
ale_matverif_indentplfoam.inp: Crushable foam plasticity with volumetric hardening.
ale_matverif_indentcrushfoam.inp: Crushable foam plasticity with isotropic hardening.
ale_matverif_indentporous.inp: Porous plasticity.
ale_matverif_indenthill.inp: Hill plasticity.
ale_matverif_indentdrucker.inp: Drucker-Prager plasticity.
ale_matverif_indentcap.inp: Drucker-Prager cap plasticity.
ale_matverif_indentratedep.inp: Mises plasticity with rate dependence.
lag_matverif_indenthyperfoam.inp: Pure Lagrangian analysis of the hyperfoam model.

Figures