This study presents the numerical free vibration analysis of aircraft wings created using different airfoil cross sections such as NACA 0009, NACA 2424, and NACA 4415. Aircraft wings were made of different lengths. Numerical frequency analyses were conducted Taguchi L9 orthogonal array with two control factors including three levels and so nine numerical modal analyses were performed. Airfoil cross sections and lengths of aircraft wings were used as the first and the second control factors. To detect the control factors with optimal levels, analysis of signal-to-noise (S/N) ratio was employed. In addition, analysis of variance (ANOVA) at the 95 % confidence level was implemented to carry out percent contributions of airfoil cross sections and lengths of aircraft wings on free vibration. As can be summarized from this study, the maximum free vibration behavior was obtained by using NACA 2424 wing profiles with a length of 5 meters. Also, the most dominant control factors were found to be airfoil type with 85.21 % effect and wing length with 12.87 % effect, according to ANOVA.