The debonding of particle/matrix interfaces has an important effect on the behavior of composite materials.During last decades,great efforts have been made to simulate the behavior of the interface.Boundary element method,as an alternative effective numerical method,has great advantages in simulating interface debonding over FEM for dimension reduction and high accuracy.In this paper,a scheme using BEM is presented to simulate the interface debonding process.