BOUNDARY ELEMENT ANALYSIS OF ELASTOPLASTIC CONTACT PROBLEMS

In this paper a boundary element method for two-dimensional elasto-plastic stress analysis of frictional contact problems is presented. The bodies in contact are treated as separate regions. The contact conditions are used to join the different system of equations for different regions of contacted bodies and hence an overall system of equations would be obtained. An incremental and iterative procedure can be used to find the contact load, or the contact extent and the proper contact conditions. To include the plastic deformation in the analysis the initial strain algorithms is employed. Elastic-perfectly plastic or work-hardening material behavior can be assumed. For the numerical analysis, isoparametric three-noded line elements are used to represent the boundary and eight-noded quadrilateral and six-noded triangular elements are used for the interior of the domain. The displacement rates and traction rates are assumed to vary quadratically and the shape functions for the interior strain rates are also of quadratic type. As an example the behavior of an elastic and elasto-plastic body with a smooth, circular inclusion under the increasing load is presented.