The analysis in this paper considers the effects of in-plane bending, i.e. about the vertical axis, on torsional and plunging bending as well as flutter instabilities in any one or more bending modes (DOFs) present. The inclusion of the in-plane bending displacements and in particular their associated velocities changes the problem from a linear one to a nonlinear one for slow flying vehicles such as UAVs and MAVs. In both instances the spatial dependences of the five governing relations are eliminated by Galerkin's method. The resulting nonlinear elastic ODEs are solved by using the Runge-Kutta numerical approach. The equivalent viscoelastic relations are nonlinear integral ordinary differential equations (IODEs) with variable coefficients, which unfortunately also can only be solved numerically thus making general parametric solutions, studies and conclusions unreachable. A number of illustrative problem examples with results are presented and discussed. The character and stability (flutter) of the solutions are examined. NOMENCLATURE a = distance as fraction of b from center of twist to aerodynamic center a ij , b ijkl = relaxation modulus coefficientsexternal vibratory force g, g w , g θ = structural damping coefficient k = reduced frequency N = summation limit in modulus, compliance Prony series O[n] = order n P , Q = viscoelastic differential operators t = time u = generalized displacement V, V (t) = free stream air velocity w = bending deflection x = {x 1 , x 2 , x 3 } = Cartesian coordinates α = trim angle of attack α = angle of zero lift α = Eq. (10) α r = rigid body angle of attack ∆p = aerodynamic pressure kl = strain components λ = exponent ρ air = air density σ kl = stress components θ = torsional deflection τ klmn = relaxation time ω = frequency 2