DuncanSelig

Plane Strain Duncan-Selig Soil Model

References

1. 10.3141/2511-07

Syntax

material DuncanSelig (1) (2) (3) (4) (5) (6) (7) (8) (9) (10) [11]
# (1) int, unique material tag
# (2) double, reference atmospheric pressure, postive, p_a
# (3) double, reference elatic modulus, E_r
# (4) double, elastic modulus exponent, n
# (5) double, reference bulk modulus, B_r
# (6) double, bulk modulus exponent, m
# (7) double, initial friction angle, \phi_i
# (8) double, friction angle slope, \Delta\phi
# (9) double, ratio of actual failure to ultimate failure, r_f
# (10) double, cohesion, c
# [11] double, density, default: 0.0

Theory

The constitutive relationship can be expressed as

\[ \dot\sigma=\dfrac{3B}{9B-E}\begin{bmatrix} 3B+E&3B-E&0\\ 3B-E&3B+E&0\\ 0&0&E \end{bmatrix}\dot\varepsilon. \]

Note it is an incremental form of the constitutive relationship. Symbols \(B\) and \(E\) denote bulk and elastic modulus, respectively.

The elastic modulus \(E\) is a function of stress.

\[ E=E_i\left(1-\dfrac{\sigma_d}{\sigma_{d,max}}\right)^2, \]

with

\[ E_i=E_r\left(\dfrac{\sigma_3}{p_a}\right)^n,\quad \sigma_{d,max}=\dfrac{2}{r_f}\dfrac{c\cos\phi+\sigma_3\sin\phi}{1-\sin\phi}. \]

The friction angle \(\phi\) decreases with increasing \(\sigma_3\).

\[ \phi=\phi_i-\Delta\psi\log_{10}\left(\dfrac{\sigma_3}{p_a}\right). \]

The deviatoric stress \(\sigma_d\) is the difference between the major and minor principal stresses.

\[ \sigma_d=\sigma_1-\sigma_3. \]

The bulk modulus \(B\) is a function of \(\sigma_3\).

\[ B=B_r\left(\dfrac{\sigma_3}{p_a}\right)^m. \]