The logical deduction on viscous-plastic numerical algorithm, in this paper, is translated under stochastic mathematical coverage coupled with stable time step by which the transcendental period of failure criterion evolution is measured in viscous-plastic stress space, by which, the nonlinear characteristics of geo-material is described deeply under Mohr-Coulomb failure criterion. Furthermore, constitution model of viscous-plastic nonlinear stochastic finite element method on 3-dimention and plane strain status is setup here on the basis of Partial Differentiation Method. Thereby, the numerical algorithm formulation on viscous-plastic nonlinear stochastic finite element method is introduced based on total strain theory. Levee structure construction as an objective case is simulated during whole applying course under the foregoing stochastic mathematical coverage, by which, the corresponding random evolution mechanism and reliability on dike structure reclamation working phase is studied comprehensively.