A FE code with large deformation contact capability and large elasto plastic analysis in Henky strain

Hi every one...



I finished my master in mechanical engineearing about eight month ago ... and in my Master career, I wrote a Fortran Finte element code with capabality of  mortar cantact analysis[1,2] and rate independent plasticity with kinematic and isotropic hardening and with large deformation capabilty for isotropic materials[3,4]. My code converges very well and I achieved quadratic convergence and I didn't use any comercial Fe packages like ansys and abaqus .....

there are two samples here for you,the example shows a cylindrical die is pressed into the slab and then slid over the surface:

(1) Mooney-Rivlin material with frictional contact

(2) Armstrong-type kinematic hardening model with henky strain and frictional contact

[1] Puso, M.A. and T.A. Laursen, A mortar segment-to-segment contact method for large deformation solid mechanics. Computer Methods in Applied Mechanics and Engineering, 2004. 193(6–8): p. 601-629.



[2] Puso, M.A. and T.A. Laursen, A mortar segment-to-segment frictional contact method for large deformations. Computer Methods in Applied Mechanics and Engineering, 2004. 193(45–47): p. 4891-4913.



[3] Montáns, F.J. and K.-J. Bathe, Computational issues in large strain elasto-plasticity: an algorithm for mixed hardening and plastic spin. International Journal for Numerical Methods in Engineering, 2005. 63(2): p. 159-196.



[4] Montáns, F.J., J.M. Benítez, and M.Á. Caminero, A large strain anisotropic elastoplastic continuum theory for nonlinear kinematic hardening and texture evolution. Mechanics Research Communications, 2012. 43(0): p. 50-56.