The symmetric boundary element method, based on the Galerkin hypotheses, has found an application in the nonlinear analysis of plasticity and in contact-detachment problems, but both dealt with separately. In this paper, we want to treat these complex phenomena together as a linear complementarity problem.A mixed variable multidomain approach is utilized in which the substructures are distinguished into macroelements, where elastic behavior is assumed, and bem-elements, where it is possible that plastic strains may occur. Elasticity equations are written for all the substructures, and regularity conditions in weighted (weak) form on the boundary sides and in the nodes (strong) between contiguous substructures have to be introduced, in order to attain the solving equation system governing the elastoplastic-contact/detachment problem. The elastoplasticity is solved by incremental analysis, called for active macro-zones, and uses the well-known concept of self-equilibrium stress field here shown in a discrete form through the introduction of the influence matrix (self-stress matrix). The solution of the frictionless contact/detachment problem was performed using a strategy based on the consistent formulation of the classical Signorini equations rewritten in discrete form by utilizing boundary nodal quantities as check elements in the zones of potential contact or detachment. the integral equations satisfying the boundary and domain conditions, the Signorini equations regarding the contact/detachment conditions rewritten as LCP, and the constitutive relations for rate-independent plasticity.