Nonlinear analysis methods are developed which will enable the reliable prediction of the dynamic behavior of the space shuttle main engine(SSME) turbopumps in the presence of bearing clearances and other local nonlinearities. A computationally efficient convolution method, based on discretized Duhamel and transition matrix integral formulations, is developed for the transient analysis. In the formulation, the coupling forces due to the onlinearities are treated as external forces acting on the coupled subsystems. Iteration is utilized to determine their magnitudes at each time increament. The method is applied to a nonlinear generic model of the high pressure oxygen turthods, the convolution approach proved to be more accurate and highly more efficient. For determining the nonlinear, steady-state periodic responses, an incremental harmonic balance(IHB) method was also developed. The method was successfully used to determine dominantly harmonic and subharmonic(subsynchronous) responses of the HPOTP generic model with bearing clearances. A reduction method similar to the impedance formulation utilized with linear systems is used to reduce the housing-totor models to their coordinates at the bearing clearances.