A computational procedure for the analysis of loads due to steady and oscillatory control surface deflections is presented. The numerical algorithm is based on a time accurate solution of the transonic full potential equation in a body-fitted coordinate system. Control surface deflection and motion are modeled using an equivalent body velocity that circumvents the need for generating time-dependent grids and interpolation between planes and discontinuity. Viscous effects, including mild separation, are modeled using an interactive inverse boundary layer and the transpiration velocity approach. Numerical results are presented for a three-dimensional transport aircraft wing with supercritical sections, for control surfaces located at the trailing edge and the leading edge, for steady as well as oscillatory deflections. Results are compared with experimental data as well as with linear theory.