Numerical investigation of hyperelastic solids with contact interaction

The paper is devoted to the construction of a computational algorithm for the investigation of hyperelastic solids with the contact interaction. In the framework of the previously developed algorithm for the investigation of large deformations of three-dimensional solids, the solutions of contact problems based on the equation of the principle of virtual work in velocity term have been considered. Contact interaction has been modeled on the basis of the “master-slave” approach. The closest point projection procedure has been used to find the contact area. A contact functional has been built on the basis of the principle of virtual work in velocity term within the penalty method. The linearization of the kinematic relations and contact functional are based on the capacity on the possible velocities of penetration. The solution of the nonlinear system of equations has been obtained using the method of step loading. The numerical implementation is based on the finite element method.