Cam Clay Soft Soil

The Cam Clay Soft Soil (CCSS) [Verneer et al., 1999; PLAXIS 2016] is an elasto-plastic rate independent model.

Theory

The strain increment \(d \boldsymbol{\varepsilon}\) is assumed to be additively decomposed into elastic \(d \boldsymbol{\varepsilon}^e\) and plastic \(d \boldsymbol{\varepsilon}^p\) parts. Dividing both sides by time increment \(d t\), the total strain rate reads:

\(\dot{\boldsymbol{\varepsilon}}=\dot{\boldsymbol{\varepsilon}}^e+\dot{\boldsymbol{\varepsilon}}^p\)

The elastic response, in rate form, reads:

\(\dot{\boldsymbol{\varepsilon}}^e=\dot{\boldsymbol{\varepsilon}}-\dot{\boldsymbol{\varepsilon}}^p = \mathbf{D}^{-1} \dot{\boldsymbol{\sigma}}\)

where \(\dot{\boldsymbol{\sigma}}\) is the rate of the internal stress \(\boldsymbol{\sigma}=\left(\sigma_x,\sigma_y,\sigma_z,\tau_{xy},\tau_{yz},\tau_{xz} \right)^T\) and \(\mathbf{D}\) is the elastic matrix. Nonlinear isotropic elasticity is assumed and so the elastic matrix \(\mathbf{D}\) is function of the Poisson coefficient \(\nu\) and the Young modulus E. The Poisson coefficient \(\nu\) is a constant input parameter while E is stress-dependent and updated according to the following law:

\(E = - \left( 1 - 2 \nu \right) \frac{\sigma_x+\sigma_y+\sigma_z}{k^*}\)

where \(k^*\) is the modified swelling parameter.

The SSCM has a yield function and a purely (nonlinear) elastic domain and the plastic strain is non-zero only if the \textit{yield criterion F} is satisfied:

\(F = \mathrm{p}_c - \mathrm{p}_{c,y}\)

where the stress \(\mathrm{p}_{c}\) is a representation of the stress state in the plane of the stress invariants (p,q):

\(\mathrm{p}_c = \mathrm{p} + \frac{c} {\tan \phi } + \frac{\mathrm{q}^2} {\mathrm{M}^2 \left(\mathrm{p} + \frac{c} {\tan \phi }\right)}\)

with:

\(\mathrm{p} = -\frac{\sigma_x+\sigma_y+\sigma_z}{3}\)

\(\mathrm{q} = \sqrt{\sigma_x \left( \sigma_x-\sigma_y \right) + \sigma_y \left( \sigma_y-\sigma_z \right) + \sigma_z \left( \sigma_z-\sigma_x \right) + 3 \left( \tau^2_{xy} + \tau^2_{yz} + \tau^2_{zx} \right) }\)

and where \(c\) is the cohesion, \(\phi\) is the friction angle, \(\mathrm{M}\) is the critical state line and the stress \(\mathrm{p}_{c,y}\) is a reference stress that describes the hardening behavior:

\(\mathrm{p}_{c,y}=\mathrm{p}_{c,y,0} \exp \left[ -\frac{\varepsilon^p_x+\varepsilon^p_y+\varepsilon^p_z}{\lambda^*-k^*} \right]\)

where \(\lambda^{*}\) is the modified compression index and \(\mathrm{p}_{c,y,0}\) is a parameter related to the plastic strain developed by the material before loading. The ratio \(\frac{\mathrm{p}_{c,y,0}}{\mathrm{p}_{c,0}}\) is the geotechnical initial Over-Consolidation-Ratio \(OCR_0\). The plastic response is governed by the flow rule:

\(\dot {\boldsymbol \varepsilon}^p = \dot \gamma \frac{\partial \mathrm{p}_c}{\partial \boldsymbol \sigma}\)

where \(\dot \gamma\) is the plastic multiplier rate. The plastic multiplier rate is rate independent and it is function of the yield criterion:

\(\dot \gamma = \begin{cases} \dot\gamma \quad F = 0 \\ 0 \quad F < 0 \end{cases}\)

The stress \(\boldsymbol \sigma\), the reference stress \(\mathrm{p}_{c,y}\) and the plastic multiplier \(\gamma\) describe the current material state in terms of stress and strain (elastic and plastic) and they are defined as the state variables of the model.

Initialization

The initialization can be done using the tab_file or the restart file.

Input Notes

In the input file materials must be set:

• The material type MAT_TYPE = 5.

• The integer parameters MAT_IPAR[NI_MAT_TYPE] with NI_MAT_TYPE = 2:

• The real parameters MAT_RPAR[NR_MAT_TYPE] with NI_MAT_TYPE = 13: