Numerical Stabilization for Flutter Analysis Procedure

Abstract

Severe mode switching is often observed when the PK-method is used in the flutter analysis of complex aircraft configurations, in particular when nearly 100 vibrational modes are considered. In the commonly used commercial software NASTRAN, the resulting eigenroots are sorted in an ascending order of frequency. Therefore, the appearance of massive mode-switching instances cannot be avoided in the PK-method flutter analyses, especially for engineering applications with real-world complex configurations. In this study, as a post-processing procedure, an extensive sorting capability was developed in order to compensate for NASTRAN’s lack of a mode-tracking procedure in between the airspeed steps. The capability was developed based on both the complex eigenvalues and their corresponding eigenvectors. In addition, numerical techniques commonly used in computational fluid dynamics (CFD) were introduced to improve the convergence of the traditional PK-method. A hybrid approach was applied to the initial guess of the reduced frequency, followed by a deferred correction scheme for the PK-iteration process. Additionally, mode matching was specifically addressed when locking eigenroots onto the aerodynamics within the PK iterations. In addition to the PK iterations, a damping iteration or modified g-method was implemented by extending the PK-method solver. The combination of these special techniques effectively improved the numerical stability of the iterations in the stability eigensolution process and significantly reduced the appearance of the misleading mode switching, minimizing risks in aircraft flight.

Keywords

• aeroelasticity
• flutter analysis
• PK-method
• <i>g</i>-method
• doublet lattice method
• mode tracking