In this paper a strategy to perform elastoplastic analysis by using the Symmetric Boundary Element Method (SBEM) for multidomain type problems is shown. This formulation uses a self-stresses equation to evaluate the trial stress in the predictor phase, and to provide the elastoplastic solution in the corrector one. Since the solution is obtained through a return mapping involving simultaneously all the plastically active bem-elements, the proposed strategy does not depend on the path of the plastic strain process and it is characterized by computational advantages due the considerable decrease of the plastic iterations number. This procedure has been developed inside Karnak.sGbem code [1] by introducing an additional module.