In this article, a new finite element model (FEM) of an EulerBernoulli beam, developed through an absolute nodal coordinate formulation (ANCF), is presented for simulation and analysis of the performance of surface-bonded piezoelectric actuators in suppressing non-linear transverse vibrations that are induced by very fast slewing. The elastic deformations experienced are an order of magnitude larger than cases considered to date, and the model employs a unique cubic spline approximation to the beam's deformed elastic line that is in terms of node positions and curvatures. To ensure relevant commentary on the vibration suppression properties of the distributed piezoelectric actuators, a material damping model was introduced in the continuum equations to capture the non-linear damping of the very slender beam that is observed in experiments. Following the ANCF methodology, the constitutive damping moment is formulated in terms of the absolute nodal coordinates with care taken to ensure the calculation is singularity free. Galerkin's method of weighted residuals is applied to discretize the revised equations of motion derived for the beam continuum. The FE beam model exploits a synergy between the twisted spline geometry and the lumped mass approximation to halve the size of the matrix equations that must be solved on each time step. However, this condensation of the matrix equations requires the use of interelement boundaries at the edges of the surface-bonded piezos. Using a single-link flexible manipulator as an example, a number of static and dynamic simulation examples that illustrate the validity of our FEM are presented, including comparisons to theoretical and other existing numerical solutions in literature. In addition, active vibration control examples are presented using proportional-and derivative-based hub motion and piezoelectric actuator controls in suppressing dramatic vibrations induced by fast slewing.