AbstractA novel approach, referred to as the homotopy stochastic finite element method, is proposed to solve the eigenvalue problem of a structure associated with some amount of uncertainty based on the homotopy analysis method. For this approach, an infinite multivariate series of the involved random variables is proposed to express the random eigenvalue or even a random eigenvector. The coefficients of the multivariate series are determined using the homotopy analysis method. The convergence domain of the derived series is greatly expanded compared with the Taylor series due to the use of an approach function of the parameter h. Therefore, the proposed method is not limited to random parameters with small fluctuation. However, in practice, only single‐variable and double‐variable approximations are employed to simplify the calculation. The numerical examples show that with a suitable choice of the auxiliary parameter h, the suggested approximations can produce very accurate results and require reduced or similar computational efforts compared with the existing methods.