This paper presents a numerical implementation of the geometrically exact beam theory based on the Legendre-spectral-finite-element (LSFE) method. The displacement-based geometrically exact beam theory is presented, and the special treatment of three-dimensional rotation parameters is reviewed. An LSFE is a high-order finite element with nodes located at the Gauss-Legendre-Lobatto points. These elements can be an order of magnitude more computationally efficient than low-order finite elements for a given accuracy level. The new module, BeamDyn, is implemented in the FAST modularization framework for dynamic simulation of highly flexible composite-material wind turbine blades within the FAST aeroelastic engineering model. The framework allows for fully interactive simulations of turbine blades in operating conditions. Numerical examples are provided to validate BeamDyn and examine the LSFE performance as well as the coupling algorithm in the FAST modularization framework. BeamDyn can also be used as a stand-alone high-fidelity beam tool.Wind power is becoming one of the most important renewable energy sources in the USA. In recent years, the size of wind turbines has been increasing immensely to lower the cost of energy, which, because of weight restrictions, requires highly flexible turbine blades. This huge electromechanical system poses a significant challenge for engineering design and analysis. Although possible with modern super computers, direct three-dimensional (3-D) structural analysis is so computationally expensive that engineers are always seeking for more efficient, highly accurate models, especially in the context of coupled aeroelastics.Beam models are widely used to represent and analyse engineering structures that have one dimension that is much larger than the other two. Many engineering components can be idealized as beams: structural members of buildings and bridges in civil engineering, joists and lever arms in heavy-machine industries and helicopter rotor blades. The blades, tower and shaft in a wind turbine system can be analysed as beams. In the weight-critical applications of beam structures, like high-aspect-ratio wings in aerospace and wind energy applications, composite materials are attractive because of their superior strength-to-weight and stiffness-to-weight ratios. However, analysis of composite-material structures is more difficult than their isotropic counterparts due to elastic-coupling effects. Furthermore, wind turbine blades are further complicated by their high flexibility and initial twist/curvatures, which must be treated in the underlying analysis. The geometrically exact beam theory (GEBT), first proposed by Reissner, 1 is a method that has proven powerful for analysis of highly flexible composite beams in the helicopter engineering community. During the past several decades, much effort has been invested in GEBT. Simo 2 and Simo and Vu-Quoc 3 extended Reissner's work to deal with 3-D dynamic problems. Jelenić and Crisfield 4 implemented GEBT using the finite-element me...