A viscoplastic constitutive model for plastic silts and clays for static slope stability applications

A viscoplastic model for representing plastic silts and clays in geotechnical static slope stability applications is presented. The PM4SiltR model builds on the stress ratio-controlled, critical state-based, bounding surface plasticity model PM4Silt and is coded as a dynamic link library for use in the finite difference program FLAC 8.1. PM4SiltR incorporates strain rate-dependent shear strength, stress relaxation, and creep using a consistency approach combined with an internal strain rate and auto-decay process. The model does not include a cap, and as such cannot simulate strain rate-dependent consolidation under increasing overburden stress. Six parameters control the viscous response for PM4SiltR while the parameters controlling the nonviscous components of the response are the same as for PM4Silt. Single element simulations are presented to illustrate the influence of viscoplasticity on the constitutive response in direct simple shear loading and undrained creep. Single element responses are shown to be consistent with observed experimental results. Simulations of a hypothetical tailings dam constructed using the upstream method are performed to illustrate use of PM4SiltR at field scale. Results of field scale simulations show PM4SiltR can model undrained creep and progressive failure leading to delayed slope instability after relatively minor changes in loading conditions at field scale.