We propose a method to treat the strain hardening effect in the hyperelastic model. Coupled with the level set interface capturing method and a real ghost fluid method, it is formulated in a fully Eulerian framework. The HLLD approximate Riemann solver is used to obtain the interface state. We use the elastic inverse deformation gradient tensor to track the material deformation state. Once the solid starts to yield, along the normal direction in stress space, an iterative algorithm is used to relax the elastic inverse deformation gradient tensor to target value where the equivalent stress meets the requirement of consistency condition in plasticity. This approach is intuitive and can avoid the use of the physical delay time, which is difficult to calibrate from experiments. We use the Wilkins flying plate impact, cylindrical shell collapse, Taylor impact and the conical-nosed projectile penetration test cases to validate the method. The results show good agreement with the theoretical analysis and experiments. This approach can also apply without further extension when thermal softening effect or damage models are considered.
Read full abstract