In this paper, the high-order perturbation techniques and Karhunen-Loéve series expansion are combined to set up a series of deterministic recursive equations first, and then the coefficients of polynomial expansion can be obtained, so the coefficients can be modified by using the Galerkin projection scheme, and more accurate expansion of structural response is definitely determined. Furthermore, an example is investigated in order to demonstrate the efficiency and precision of the proposed method, the results show that this method has high calculation accuracy and much higher efficiency than traditional direct Monte Carlo simulation.