In order to improve the accuracy of concrete faced rockfill dam (CFRD) settlement, the stochastic finite element method (SFEM) was introduced into the face rockfill dam and the applicability of the E-B constitutive model is enhanced by considering the spatial variability of the parameters of the E-B constitutive model, and a noninvasive stochastic finite element method (NSFEM) for the CFRD based on the Cholesky decomposition was established. First, the random field of parameters of rockfill materials was simulated; Second, the secondary development was achieved based on the existing commercial finite element software MSC.Marc. By calling Marc directly in DOS environment through UltraEdit, each unit was assigned different material parameters, thus the “noninvasive” stochastic finite element analysis method was established. Then, combined with the settlement calculation of Gongboxia dam, the in1uence law of the coefficient of variation of dam material parameters on settlement was analyzed. The calculation results showed that: considering the spatial variability of material parameters, the maximum settlement was obviously larger than the conventional finite element calculation value, which was closer to the actual observed value and more in line with the actual situation of the project. We conclude that the settlement is ranked from large to small: φ0、K、Kb，according to the sensitivity of the material parameters and the sensitivity of the material parameters of the 3BII area is larger than that of the 3C area. Using NSFEM and considering the spatial variability of dam material parameters is of great significance for accurately grasping CFRD settlement and ensuring engineering safety.