Modified p-k method for flutter solution with damping iteration

The g method equation is reformulated and solved by a modified p-k method with damping iteration and implications of the second order of damping term are also drawn. Linear unsteady aerodynamics computed by panel method in the frequency domain or by linearized frequency-domain CFD results are assumed. To adapt the standard p-k method to the dependence of matrix A on damping, it is suggested to include the damping iteration in addition to the reduced-frequency iteration in the solution algorithm of the standard p-k method. The generalized mass, stiffness, and aerodynamic matrices are obtained from ZAERO output file. The comparison of the flutter characteristics between the modified p-k method and the g method show a good agreement between these methods. The V-g and V-f diagrams seem to be almost identical between both approaches.