Abstract. Due to the uncertainty of the initial geometrical imperfection, it is difficult to evaluate its influence on bucking behavior of structures. In this study, the double trigonometric series is adopted to describe the initial geometric imperfection and the deformed deflection, and then the relation between them is deduced from the bending equilibrium equation. Therefore, the total displacement can be obtained by superposing these two deflections. Based on this analytical model, the range of the maximum deformation is determined by Monte Carlo simulation. The analysis provides a new analytical model for the buckling problem of plate with uncertain initial geometrical imperfection.