Current challenges in industrial multibody system simulation are often beyond the classical range of application of existing industrial simulation tools. The present paper describes an extension of a recursive order-n multibody system (MBS) formulation to nonlinear models of flexible deformation that are of particular interest in the dynamical simulation of wind turbines. The floating frame of reference representation of flexible bodies is generalized to nonlinear structural models by a straightforward transformation of the equations of motion (EoM). The approach is discussed in detail for the integration of a recently developed discrete Cosserat rod model representing beamlike flexible structures into a general purpose MBS software package. For an efficient static and dynamic simulation, the solvers of the MBS software are adapted to the resulting class of MBS models that are characterized by a large number of degrees of freedom, stiffness, and high frequency components. As a practical example, the run-up of a simplified three-bladed wind turbine is studied where the dynamic deformations of the three blades are calculated by the Cosserat rod model.