Abstract: A nonlinear complementarity approach is presented to solve elastoplastic problems by boundary element method, in which the equations are formulated by stress equations of internal points and complementarity function obtained from plasticity constitutive law. The domain integrals involved are transformed into boundary integrals by radial integration method. The algorithm’s precision is demonstrated by stress, displacement and load-displacement results of a thick-walled cylinder subjected to internal pressure.