# Impact of an elastoplastic bar¶

## Overview¶

The aim of this test case is to validate the following functions:

• Dynamic analysis
• Frictional contact
• Elastoplasticity
• Large deformations

The simulation results of SimScale were compared to the numerical results presented in case B of [SDNV103]. In addition to the quadratic mesh consisting of prisms and hexahedrals described in [SDNV103], the linear hexahedral mesh was also used.

Import validation project into workspace

## Geometry¶

Geometry of the elastoplastic bar and rigid plate

A B C D E F G H I J K L M N
x [m] 0 0.016 0.016 0 0 0.016 0.016 0 0 0.0032 0 0 0.0032 0
y [m] 0 0 0.016 0.016 0 0 0.016 0.016 0 0 0.0032 0 0 0.0032
z [m] 0 0 0 0 0.001 0.001 0.001 0.001 0.00327 0.00327 0.00327 0.03567 0.03567 0.03567

## Analysis type and Domain¶

Tool Type : Code-Aster

Analysis Type : Dynamic (nonlinear elastoplastic)

Mesh and Element types :

Case Mesh type Number of nodes Number of prisms Number of hexahedrals Element type
(A) linear prisms and hexahedrals 764 120 481 3D isoparametric
(B) quadratic prisms and hexahedral 2861 120 481 3D isoparametric reduced integration
(C) quadratic prisms and hexahedral 2861 120 481 3D isoparametric standard

Mesh used for the SimScale

## Simulation Setup¶

Material:

• isotropic: E = 117 GPa, $$\nu$$ = 0.35, $$\sigma_y$$ = 400 MPa, ET = 500 MPa

Dependency of the Young’s modulus

Initial Boundary Condition:

• Velocity of 227 m/s to volume IJKLM (bar)

Constraints:

• Volume ABCDEFGH (plate) zero x-y-z-displacement
• Face IKLN zero x-displacement
• Face IJLM zero y-displacement

## Results¶

Comparison of the displacement Dx and Dy computed on points J and L respectively at time 9e-5 obtained with SimScale and the results presented in case B of [SDNV103].

Comparison of Dx and Dy in [m] at point J and L
Case Point [SDNV103] ¹ [SDNV103] ² SimScale Error (%) - [SDNV103] ¹ Error (%) - [SDNV103] ²
(A) J (Dx) 0.00387 0.00291 0.00131 66.202 55.07
L (Dz) -0.01346 -0.01269 -0.01386 2.998 9.247
(B) J (Dx) 0.00387 0.00291 0.00287 25.889 1.48
L (Dz) -0.01346 -0.01269 -0.01266 5.983 0.278
(C) J (Dx) 0.00387 0.00291 0.00277 28.47 4.911
L (Dz) -0.01346 -0.01269 -0.01269 5.738 0.018

## References¶

 [SDNV103] (1, 2, 3, 4, 5, 6, 7, 8, 9) SDNV103 - Impact d’une barre de Taylor élastoplastique

¹ Results by Stainier et al. as mentioned in [SDNV103]: L. Stainier and J.Ph. Ponthot. An improved one-point integration method for large strain elastoplastic analysis. Computer Methods in Applied Mechanics and Engineering, 118(1–2):163 – 177, 1994

² Code_Aster results as mentioned in case B of [SDNV103]