Abstract: This paper presents analysis methods and the corresponding computer program for the problem in question. First, the static finite element model for wing-type composite structures is formulated. A variable-linking technique is used to reduce the order of design-variable space. The transformation from the static model to the dynamic one is done by use of static compliance method. Several algorithms for eigenproblems are employed to analyze natural vibration of structures. The opitimal design with the frequency constraint is studied by means of analytical derivatives. 3-D unsteady aerodynamic analysis of harmonically oscillating surfaces as well as flutter analysis is done using a subsonic doublet-lattice method. Lastly, the paper presents derivatives of the unsteady aerodynamic load with respect to design variables and investigates the optimum design subjected to flutter constraint by means of a feasible-direction method.