4.15.1.13. PM4Silt¶
- nDMaterial('PM4Silt', matTag, S_u, Su_Rat, G_o, h_po, Den <Su_factor, Patm, nu, nG, h0, eInit, lambda, phicv, nb_wet, nb_dry, nd, Ado, ru_max, zmax, cz, ce, Cgd, ckaf, m_m, CG_consol>)
Code Developed by: Long Chen and pedro at U.Washington.
This command is used to construct a 2-dimensional PM4Silt material.
matTag(int)integer tag identifying material
S_u(float)Undrained shear strength
Su_Rat(float)Undrained shear strength ratio.
G_o(float)Shear modulus constant
h_po(float)Contraction rate parameter
Den(float)Mass density of the material
Su_factor(float)Optional: Undrained shear strength reduction factor
P_atm(float)Optional: Atmospheric pressure
nu(float)Optional: Poisson’s ratio. Default value is 0.3.
nG(float)Optional: Shear modulus exponent. Default value is 0.75.
h0(float)Optional: Variable that adjusts the ratio of plastic modulus to elastic modulus. Default value is 0.5.
eInit(float)Optional: Initial void ratios. Default value is 0.90.
lambda(float)Optional: The slope of critical state line in e-ln(p) space. Default value is 0.060.
phicv(float)Optional: Critical state effective friction angle. Default value is 32 degrees.
nb_wet(float)Optional: Bounding surface parameter for loose of critical state conditions \(1.0 \geq nb_wet \geq 0.01\). Default value is 0.8. in cyclic loading
nb_dry(float)Optional: Bounding surface parameter for dense of critical state conditions \(nb_dry \geq 0\). Default value is 0.5.
nd(float)Optional: Dilatancy surface parameter \(nd \geq 0\). Default value is 0.3.
Ado(float)Optional: Dilatancy parameter. Default value is 0.8. with accumulation of fabric
ru_max(float)Optional: Maximum pore pressure ratio based on p’.
z_max(float)Optional: Fabric-dilatancy tensor parameter
cz(float)Optional: Fabric-dilatancy tensor parameter. Default value is 100.0.
ce(float)Optional: Variable that adjusts the rate of strain accumulation in cyclic loading
cgd(float)Optional: Variable that adjusts degradation of elastic modulus with accumulation of fabric. Default value is 3.0. in stress ratio space)
ckaf(float)Optional: Variable that controls the effect that sustained static shear stresses have on plastic modulus. Default value is 4.0.
m_m(float)Optional: Yield surface constant (radius of yield surface in stress ratio space). Default value is 0.01.
CG_consol(float)Optional: Reduction factor of elastic modulus for reconsolidation. \(CG_consol \geq 1\). Default value is 2.0.
Note
The material formulation for the PM4Silt object is PlaneStrain
If both S_u and Su_Rat values are specified the value of S_u is used.
Valid Element Recorder queries are stress, strain, alpha (or backstressratio) for \(\mathbf{\alpha}\), fabric for \(\mathbf{z}\), and alpha_in (or alphain) for \(\mathbf{\alpha_{in}}\)
Elastic or response could be enforced by
updateMaterialStage('-material', matTag, '-stage', 0)
Elastoplastic by
updateMaterialStage('-material', matTag, '-stage', 1)
The program will use the default value of a secondary parameter if a negative input is assigned to that parameter, e.g. Ado = -1. However, FirstCall is mandatory when switching from elastic to elastoplastic if negative inputs are assigned to stress-dependent secondary parameters, e.g. Ado and zmax. FirstCall can be set as,
setParameter('-value', 0, '-ele', elementTag, 'FirstCall', matTag)
Post-shake reconsolidation can be activated by
setParameter('-value', 1, '-ele', elementTag, 'Postshake', matTag)
The user should check that the results are not sensitive to time step size.
Refer to PM4Sand for examples.
R.W.Boulanger, K.Ziotopoulou. “PM4Silt(Version 1): A Silt Plasticity Model for Earthquake Engineering Applications”. Report No. UCD/CGM-18/01 2018