SUMMARYIn this paper, the novel concept of probabilistic yielding is used for 1-D cyclic simulation of the constitutive behavior of geomaterials. Fokker-Planck-Kolmogorov equation-based probabilistic elastic-plastic constitutive framework is applied for obtaining the complete probabilistic (probability density function) material response. Both perfectly plastic and hardening-type material models are considered. It is shown that when uncertainties in material parameters are taken into consideration, even the simple, elastic-perfectly plastic model captures some of the important features of geomaterial behavior, for example, modulus reduction with cyclic strain, which, deterministically, is only possible with more advanced constitutive models. Furthermore, it is also shown that the use of isotropic and kinematic hardening rules does not significantly improve the probabilistic material response.