Abstract: Based on the influence of the stochastic nature of the parameter in the perturbation stochastic FEM (finite element method) on the coupled process of the seepage field and the stress field, finite element iterative formulas are derived, where the permeability coeffcient and the elasticity modulus are set as the basic stochastic quantities, and accordingly, the stochastic calculating programs used for calculating the seepage field and the stress field separately are re-formulated. By utilizing the improved programs, the H-M problem of a dam is calculated, and the results are compared with those calculated by COMSOL Multiphysics. It is shown that these two methods are consistent, which indicates that for fluid-solid coupling problems in engineering practice, the stochastic characteristics can be obtained by the program, which is not available in the commercial software.