Abstract. Beam models are used for the aeroelastic time and frequency domain analysis of wind turbines due to their computational efficiency. Many current aeroelastic tools for the analysis of wind turbines rely on Timoshenko beam elements with classical crosssectional properties (EA, EI, etc.). Those cross-sectional properties do not reflect the various couplings arising from the anisotropic behaviour of the blade material. A twonoded, three-dimensional Timoshenko beam element was therefore extended to allow for