The use of the subregion technique is a common way to deal with potential and elastic zoned problems by the boundary element method. Although this procedure has been successfully applied in many practical engineering analyses, nonlinear solutions can be easily deteriorated in the presence of an increasing number of discretized interface lines. This paper presents an appropriate formulation to model this kind of problem, avoiding some interface approximations, and reducing, therefore, the possible errors in the final solution.