IntroductionI N the analysis point of view, the composite rotor blade has typically been analyzed through the one-dimensional beam assumption since the spanwise length of rotor blades is generally much longer than their lateral dimensions. In developing the beam theory, there may be coupling among extension, bending, and torsional deformations. These couplings generally invalidate the Euler-Bernoulli assumption: plane sections remain plane and are perpendicular to the elastic axis. The assumption leads to underestimation of beam displacements, especially in case of bending, because of constant shear distribution across the beam section. Moreover, for a composite beam in bending, this distribution of shear is nearly parabolic (piecewise in general). 1 The use of the shear correction factor (SCF) may be the most economical one for the transverse shear behavior without largely sacrificing the required accuracy of solution. In addition, warping and warping inhibition effects are to be considered in the analysis. 2 Therefore, an appropriate analytical model capturing these behaviors is inevitable to get more enhanced results from the aeroelastic analysis of a composite rotor blade.Hong and Chopra 3 used the nonlinear kinematic model of Hodges and Dowell 4 and extended it to the case of composite rotor. They used a simplified beam model, in which the transverse shear flexibility was not included in the formulation. Smith and Chopra 5 modified this one to include the transverse shear effects and other secondary structural modeling effects. They focused on the behavior of an elastically tailored composite blade and presented various results for vibratory hub loads of box-beam having different ply configurations, but they did not go further to consider the sectional distribution of shear stresses. To improve the theoretical results, an alternative approach, which has been developed by Jung and Kirn 6 for the effects of transverse shear and structural damping on the aeroelastic response of composite rotor blades in hover, involved the usage of SCF to account for the sectional distribution of shear stresses. They showed that the effects of transverse shear and structural damping can have a key role on the flutter boundary of the rotor, but the lay-up structure is confined to symmetric configuration only. In the present work, the formulation of Ref. 6, which considers the effects of transverse shear flexibility, torsion warping, and two-dimensional in-plane behavior, is extended to analyze arbitrary lay-up geometry including antisymmetric configuration. Numerical simulations are performed for a specific antisymmetric configuration to identify the transverse shear behavior on the aeroelastic stability of composite rotor.
Problem FormulationThe rotor structure is idealized as a laminated composite boxbeam whose constituent laminae are characterized by different ply orientation angles and different material and thickness properties as depicted in Ref. 6. The deformation of the blade in space is described by u, v, w, and