Taylor series numerical method is a novel time integration method for structural dynamics. In comparison with the well-known ones, Taylor series method has high accuracy and good convergence characteristics and thus is a good alternative for solving structural dynamics problems. In the paper, the recursive formulae of Taylor series numerical method in linear and nonlinear structural dynamics are derived. The desired user routine in Fortran 90, which was compiled and linked into the general software ANSYS program, is designed and programmed. User-defined menus are constructed via the User Interface Design Parameter. Numerical examples illustrated the validity of the customized plug-in. The research in this paper promotes the algorithm application in complicated engineering problems.