Industrial machines are sometimes exposed to the danger of earthquake. In the design of a mechanical system, this factor should be accounted for from the viewpoint of reliability to analyze a complex nonlinear structure system under random excitation is proposed. First, the actual random excitation, such as earthquake, is approximated to the corresponding Gaussian process for the statistical analysis. The modal equations of overall system are expanded sequentially. Then, the perturbed equations are synthesized into the overall system and solved in probabilistic way. Several statistical properties of a random process that are of interest in random vibration are evaluated in each substructure. Comparing with the results of the numerical simulation proved the efficiency of the proposed method.