For nonlinear system,solving its dynamic ordinary differential equations (ODEs) is always difficult.Traditional time-domain integral method,such as central difference method and classical fourth-order Runge-Kutta (RK4),cannot accurately identify the inflection point between loading and unloading states of a nonlinear frictional system,further possibly lead to inaccurate results and even unacceptable errors.