Impact of an Elastoplastic Bar

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

	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

Simulation Setup

Material:

• isotropic: E = 117 GPa, ν
  $\nu$

   

  = 0.35, σy

  ${\sigma }_y$

   

  = 400 MPa, E_T = 500 MPa

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]

Last updated: January 29th, 2019

Did you find this article helpful?

How can we do better?

We appreciate and value your feedback.

Send Your Feedback

Contents

Cookie Policy

SimScale uses cookies to improve your user experience. By using this website you consent to our cookie policy. Don't worry, we are not storing any personal information.

Accept Data Privacy