The nonlinear dynamic response, Limit Cycle Oscillations (LCOs), of High-Aspect-Ratio (HAR) wings using novel indicial aerodynamics in subsonic ow was investigated. Using the nonlinear beam theory, the structural model was derived with in-plane and out-of-plane bending and torsional motions, all nonlinearities up to cubic order arising from large deformation, mass distribution, and cross-sectional mass imbalance. Based on new approximations of the indicial functions, a comprehensive unsteady aerodynamic model was used. Coupling such indicial aerodynamics to nonlinear structural equations can result in a uni ed nonlinear aeroelastic formulation in both incompressible and subsonic compressible ows. The e ects of ight conditions, wing tip initial disturbances, Sti ness Ratio (SR) between bending modes, and nonlinearity due to inertia and cross-sectional mass imbalance on the characteristics of LCO are discussed. The results showed that compressibility could a ect the LCO boundary up to 12 percent, which implied that appropriate Mach-dependent aerodynamics was required to achieve a more reasonable and realistic description of dynamic behavior of the system. It was observed that the presence of LCO before the linear utter speed depended on the initial disturbances as well as wing characteristics.