Abstract
In this paper, a two-stage approach is presented for analyzing flexible beams undergoing large rotations. In the first stage, the symbolic forms of equations of motion and the Jacobian matrix are generated by means of MATLAB and written into a MATLAB script file automatically, where the flexible beams are described by the unified formulation presented in our previous paper. In the second stage, the derived equations of motion are solved by means of implicit numerical methods. Several comparison computations are performed. The two-stage approach proves to be much more efficient than pure numerical one.