QzSimple1 Material

• Command_Manual
  • Tcl Commands
  • Modeling_Commands
    • model
    • uniaxialMaterial
    • ndMaterial
    • frictionModel
    • section
    • geometricTransf
    • element
    • node
    • sp commands
      • fix
      • fixX
      • fixY
      • fixZ
    • mp commands
      • equalDOF
      • rigidDiaphragm
      • rigidLink
    • timeSeries
    • pattern
    • mass
    • block commands
      • block2D
      • block3D
    • region
    • rayleigh
  • Analysis Commands
    • constraints Command
    • numberer Command
    • system Command
    • test Command
    • algorithm Command
    • integrator Command
    • analysis Command
    • eigen Command
    • analyze Command
  • Output Commands
    • recorder
    • record
    • print
    • printA
    • logFile
    • RealTime_Output_Commands
  • Misc Commands
    • eleResponse Command
    • getTime Command
    • getEleTags Command
    • getNodeTags
    • loadConst Command
    • nodeCoord Command
    • nodeDisp Command
    • nodeVel Command
    • nodeAccel Command
    • nodeEigenvector Command
    • nodeBounds
    • print Command
    • remove Command
    • reset Command
    • testIter
    • testNorms
    • wipe Command
    • exit Command
  • DataBase Commands
    • Template:DataBaseCommandsMenu

---

This command is used to construct a QzSimple1 uniaxial material object:

uniaxialMaterial QzSimple1 $matTag $qzType $qult $Z50 <$suction $c>

---

$matTag	integer tag identifying material
$qzType	qzType = 1 Backbone of q-z curve approximates Reese and O'Neill's (1987) relation for drilled shafts in clay. qzType = 2 Backbone of q-z curve approximates Vijayvergiya's (1977) relation for piles in sand.
$qult	Ultimate capacity of the q-z material. SEE NOTE 1.
$Z50	Displacement at which 50% of qult is mobilized in monotonic loading. SEE NOTE 2.
$suction	Uplift resistance is equal to suction*qult. Default = 0.0. The value of suction must be 0.0 to 0.1.*
$c	The viscous damping term (dashpot) on the far-field (elastic) component of the displacement rate (velocity). Default = 0.0. Nonzero c values are used to represent radiation damping effects.*

NOTES:

1. $qult: Ultimate capacity of the q-z material. Note that "q" or "qult" are stresses [force per unit area of pile tip] in common design equations, but are both loads for this uniaxialMaterial [i.e., stress times tip area].
2. $Y50: Displacement at which 50% of pult is mobilized in monotonic loading. Note that Vijayvergiya's relation (qzType=2) refers to a "critical" displacement (zcrit) at which qult is fully mobilized, and that the corresponding z50 would be 0. 125zcrit.
3. optional args $suction and $c must either both be omitted or both provided.

THEORY:

The equations describing QzSimple1 behavior are similar to those for p-y materials by Boulanger, R. W., Curras, C. J., Kutter, B. L., Wilson, D. W., and Abghari, A. (1999). "Seismic soil-pile-structure interaction experiments and analyses." Journal of Geotechnical and Geoenvironmental Engineering, ASCE, 125(9): 750-759. Modifications were required for representing the different responses of a <math>q-z</math> material in compression versus uplift. The nonlinear q-z behavior is conceptualized as consisting of elastic (q-<math>z^e</math>), plastic (q-<math>z^p</math>), and gap (<math>q-z^g</math>) components in series. Radiation damping is modeled by a dashpot on the “far-field” elastic component (<math>q-z^e</math>) of the displacement rate. The gap component consists of a bilinear closure spring (qc-zg) in parallel with a nonlinear drag spring (qd-zg). Note that <math>z = z^e + z^p + z^g</math>, and that <math>q = q^d + q^c</math>.

The plastic component has an initial range of rigid behavior between <math>-Cr q_\text{ult} < q < Cr q_\text{ult}</math> with Cr = the ratio of <math>\frac{q}{q_{\text{ult}}}</math> when plastic yielding first occurs in virgin loading. The rigid range of q, which is initially <math>2 Cr q_\text{ult}</math>, translates and grows with plastic yielding. The rigid range of q is constrained to a maximum size of <math>0.7q_{\text{ult}}</math>. Beyond the rigid range, loading of the plastic (<math>q-z^p</math>) component is described by:

<math>q = q_{\text{ult}} - (q_{\text{ult}} - q_0) \left [\frac{c * z_{50}}{c * z_{50} + | z_p - z^p_0|} \right ] </math>

where <math>q_ult</math> = the ultimate resistance of the <math>q-z</math> material in the current loading direction, <math>q_o = q</math> at the start of the current plastic loading cycle, p <math>z^p_o = z^p</math> at the start of the current plastic loading cycle, and c and n are constants that control the shape of <math>q-z^p</math> curve.

The closure (<math>q^c-z^g</math>) component is simply a bilinear elastic spring, which is relatively rigid in compression and extremely flexible in tension (uplift).

The nonlinear drag (<math>q^d-z^g</math>) component is used to allow thethe specification of some minimum “suction” on the pile tip during uplift. It is described by:

<math>q^d = C_d q_\text{ult} - (C_d q_{\text{ult}} - q^d_0) \left [\frac{z_{50}}{z_{50} + 2| z^g - z^g_0|} \right ] </math>

where <math>C_d</math> = ratio of the maximum drag (suction) force to the ultimate resistance of the <math>q-z</math> material, <math>q^d_o = q^d</math> at the start of the current loading cycle, and <math>z^g_o = z^g</math> at the start of the current loading cycle.

The flexibility of the above equations can be used to approximate different q-z backbone relations. Reese and O’Neill’s (1987) recommended backbone for drilled shafts in clay is closely approximated using <math>c = 0.35</math>, <math>n = 1.2</math>, and <math>C_r = 0.2</math>. Vijayvergiya’s (1977) recommended backbone for piles in sand is closely approximated using <math>c = 12.3</math>, <math>n = 5.5</math>, and <math>C_r = 0.3</math>.

QzSimple1 is currently implemented to allow use of these two default sets of values. Values of <math>q_\text{ult}</math>, <math>z_50</math>, and suction (i.e., Cd) must then be specified to define the <math>q-z</math> material behavior.

Viscous damping on the far-field (elastic) component of the <math>q-z</math> material is included for approximating radiation damping. For implementation in OpenSees the viscous damper is placed across the entire material, but the viscous force is calculated as proportional to the component of velocity (or displacement) that developed in the far-field elastic component of the material. For example, this correctly causes the damper force to become zero during load increments across a fully formed gap in uplift. In addition, the total force across the <math>q-z</math> material is restricted to <math>q_\text{ult}</math> in magnitude so that the viscous damper cannot cause the total force to exceed the near-field soil capacity. Users should also be familiar with numerical oscillations that can develop in viscous damper forces under transient loading with certain solution algorithms and damping ratios. In general, an HHT algorithm is preferred over a Newmark algorithm for reducing such oscillations in materials like QzSimple1.

Examples of the monotonic backbones and cyclic loading response of QzSimple1 are given in the following plots.

EXAMPLE:

REFERENCES:

"Seismic Soil-pile-strcture interaction experiments and analysis", Boulanger, R.W., Curras, C.J., Kutter, B.L., Wilson, D.W., and Abghari, A. (1990). Jornal of Geotechnical and Geoenvironmental Engineering, ASCS, 125(9):750-759.

---

Code Developed by: Ross Boulanger, UC Davis This command is used to construct a PySimple1 uniaxial material object: