On the finite element modeling of the asymmetric cracked rotor

Abstract

The advanced phase of the breathing crack in the heavy duty horizontal rotor system is expected to be dominated by the open crack state rather than the breathing state after a short period of operation. The reason for this scenario is the expected plastic deformation in crack location due to a large compression stress field appears during the continuous shaft rotation. Based on that, the finite element modeling of a cracked rotor system with a transverse open crack is addressed here. The cracked rotor with the open crack model behaves as an asymmetric shaft due to the presence of the transverse edge crack. Hence, the time-varying area moments of inertia of the cracked section are employed in formulating the periodic finite element stiffness matrix which yields a linear time-periodic system. The harmonic balance method (HB) is used for solving the finite element (FE) equations of motion for studying the dynamic behavior of the system. The behavior of the whirl orbits during the passage through the subcritical rotational speeds of the open crack model is compared to that for the breathing crack model. The presence of the open crack with the unbalance force was found only to excite the 1/2 and 1/3 of the backward critical whirling speed. The whirl orbits in the neighborhood of these subcritical speeds were found to have nearly similar behavior for both open and breathing crack models. While unlike the breathing crack model, the subcritical forward whirling speeds have not been observed for the open crack model in the response to the unbalance force. As a result, the behavior of the whirl orbits during the passage through the forward subcritical rotational speeds is found to be enough to distinguish the breathing crack from the open crack model. These whirl orbits with inner loops that appear in the neighborhood of the forward subcritical speeds are then a unique property for the breathing crack model.