Absolute nodal coordinate finite element transfer matrix method is an effective strategy for dynamic simulation of flexible beam undergoing large overall motion. Its remarkable characteristic is that no global system dynamics equation is needed formally. However the method also depends on numerical integration method of ordinary differential equation, which is same with conventional dynamics method. In this investigation, by using eight different integration strategies in absolute nodal coordinate finite element transfer matrix method, such as classical four order Runge-Kutta method, explicit four order Adams method, two improved Newmark methods and four internally installed methods of MATLAB, the dynamic simulation of flexible beam undergoing large overall motion has been implemented. As to computational stability, numerical integration precision and calculated efficiency, several useful conclusions can be found to improve the numerical characteristics of absolute nodal coordinate finite element transfer matrix method.