NUMERICAL ANALYSIS OF ELASTIC-PLASTIC CONTACT PROBLEMS

The algorithm for numerical realization of boundary conditions of contact interaction of deformable bodies taking into account a friction is worked out. Contact conditions are formulated in the form of quasivariational inequality with one-sided constraint. This constraint corresponds to the condition of nonpenetration of deformable bodies into each other. Dynamic interaction of elastic-plastic materials with beforehand unknown, time-dependent zone of contact is described by the model taking into account small strains and finite rotations. This model consists of the system of equations of motion, the Hooke law for elastic constituents of the strain tensor, the equation for the rotation angle and the principle of maximum of the energy dissipation rate describing the process of plastic flow. The transition of material from elastic state to plastic one is determined by the Mises yield condition. Parallel shock-capturing algorithm is proposed for implementation of the model on multiprocessor computer systems. It is based on a combination of splitting methods with respect to physical processes and spatial variables. The results of computations of an oblique impact of two deformable plates are represented.