Docs

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

ImpactOfAnElastoplasticBar-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
ImpactOfAnElastoplasticBar-mesh

Mesh used for the SimScale

Simulation Setup

Material:

  • isotropic: E = 117 GPa, \(\nu\) = 0.35, \(\sigma_y\) = 400 MPa, ET = 500 MPa
ImpactOfAnElastoplasticBar-material

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]