As one of the approaches to solve the problem, aeroelastic tailoring techniques has been applied to the optimal design of composite wing structures [3, 4].
For aeroelastic tailoring, this becomes more concerned since such calculation is usually conducted in the gradient calculation of each design variable and each iterative cycle in a traditional optimization process.
However, when more detailed structure optimization such as aeroelastic tailoring is required, dynamic response result is required in each iteration step or more frequently.