For the complex dynamics modelling process,low computational efficiency accuracy of large scale multibody system,an efficient and precise dynamics modelling method is proposed based on spatial operator algebraic.The method combines transformation operators,which include recursive relations of system topology structure,motion and strength,with recursive dynamics of Newton-Euler through the dynamics variables expressed by screw.This method achieves forward backward recursion of generalized velocity,acceleration force,strength and mass,and avoids cross operation and unnecessary integrals.It is valuable in scientific and engineering aspects for its brief modeling process,clear physical meaning,and suitability in computer programming and calculation.An example is given to verify this method.