The efficient nonlinear vibration analysis of a rotating elastic structure with a crack is examined. In particular, the solution of the forced vibration response of a cracked turbine engine blade is investigated. Starting with a finite element model of the cracked system, the Craig–Bampton method of component mode synthesis is used to generate a reduced-order model that retains the nodes of the crack surfaces as physical degrees of freedom. The nonlinearity due to the intermittent contact of the crack surfaces, which is caused by the opening and closing of the crack during each vibration cycle, is modeled with a piecewise linear term in the equations of motion. Then, the efficient solution procedure for solving the resulting nonlinear equations of motion is presented. The approach employed in this study is a multiharmonic hybrid frequency∕time-domain technique, which is an extension of the traditional harmonic balance method. First, a simple beam model is used to perform a numerical validation by comparing the results of the new method to those from transient finite element analysis (FEA) with contact elements. It is found that the new method retains good accuracy relative to FEA while reducing the computational costs by several orders of magnitude. Second, a representative blade model is used to examine the effects of crack length and rotation speed on the resonant frequency response. Several issues related to the rotation are investigated, including geometry changes of the crack, shifts in resonant frequencies, and the existence of multiple solutions. For the cases considered, it is found that the nonlinear vibration response exhibits the jump phenomenon only when rotation is included in the model.