A dual-time scheme is used to solve the unsteady compressible Navier-Stokes equations to obtain the needed airload.The model structural equations of motions are computed simultaneously using the Runge-Kutta method to calculate the aeroelastic static deformations and responses of wings in the time-domain.The O-H type grid around the deforming wings is generated using the improved transfinite interpolation method so that grid generation in each time-step is not too time-consuming.The numerical results show the transonic flutter speed drop and nonlinear characteristics.