The method of component mode synthesis is developed to determine the forced response of nonlinear, multishaft, rotor-bearing systems. The formulation allows for simulation of system response due to blade loss, distributed unbalance, base shock, maneuver loads, and specified fixed frame forces. The motion of each rotating component of the system is described by superposing constraint modes associated with boundary coordinates and constrained precessional modes associated with internal coordinates. The precessional modes are truncated for each component and the reduced component equations are assembled with the nonlinear supports and interconnections to form a set of nonlinear system equations of reduced order. These equations are then numerically integrated to obtain the system response. A computer program, which is presently restricted to single shaft systems has been written and results are presented for transient system response associated with blade loss dynamics, with squeeze film dampers, and with interference rubs.