Abstract
This paper investigates the local stability analysis of periodic solutions corresponding to the nonlinear vibration response of an industrial compressor blade, NASA rotor 37, on which are applied different types of nonlinearities. These solutions are obtained using a harmonic balance method-based approach presented in a previous paper. It accounts for unilateral contact and dry friction of the rotating blade against a rigid casing through a regularized penalty law. A Lanczos filtering technique is also employed to mitigate spurious oscillations related to the Gibbs phenomenon thus enhancing the robustness of the solver. In addition, a component mode synthesis technique is used to reduce the dimension of the numerical model. Stability assessment of the computed solutions relies on Floquet theory. It is performed through the computation of the monodromy matrix as well as Hill's method. Both methodologies are applied and thoroughly compared as the severity of the nonlinearity is gradually increased from a cubic spring to three-dimensional contact conditions on a deformed casing. While the presented results underline the applicability of both stability assessment methodologies for all types of nonlinearities, they also put forward the much higher computational effort required when computing the monodromy matrix. Indeed, it is shown that Hill's method yields converged results for significantly lower values of both the number of retained harmonics and the considered number of time steps thus making it a far more efficient method when dealing with industrial models. It is also underlined that the presented results are in excellent agreement with reference solution points obtained with time domain solution methods. Specific implementation tweaks that were found to be of critical importance in order to efficiently assess the stability of computed solutions are also detailed in order to provide a comprehensive view of the challenges inherent to such numerical developments.
