QzSimple1 Material: Difference between revisions

From OpenSeesWiki
Jump to navigation Jump to search
(Created page with '{{CommandManualMenu}} This command is used to construct a QzSimple1 uniaxial material object: {| | style="background:yellow; color:black; width:800px" | '''uniaxialMaterial Qz...')
 
No edit summary
 
(One intermediate revision by the same user not shown)
Line 4: Line 4:


{|  
{|  
| style="background:yellow; color:black; width:800px" | '''uniaxialMaterial QzSimple1 $matTag $qzType $qult $Y50 <$suction $c>'''
| style="background:yellow; color:black; width:800px" | '''uniaxialMaterial QzSimple1 $matTag $qzType $qult $Z50 <$suction $c>'''
|}
|}


Line 18: Line 18:
|  '''$qult ''' || Ultimate capacity of the q-z material. SEE NOTE 1.  
|  '''$qult ''' || Ultimate capacity of the q-z material. SEE NOTE 1.  
|-
|-
|  '''$Y50 ''' || Displacement at which 50% of qult is mobilized in monotonic loading. SEE NOTE 2.
|  '''$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.*
| '''$suction''' || Uplift resistance is equal to suction*qult. Default = 0.0. The value of suction must be 0.0 to 0.1.*
Line 32: Line 32:




THEORY:
EQUATIONS and EXAMPLE RESPONSES:


The equations describing QzSimple1 behavior are similar to those for p-y materials by
The equations describing QzSimple1 behavior are similar to those for p-y materials by
Line 38: Line 38:
soil-pile-structure interaction experiments and analyses." Journal of Geotechnical and
soil-pile-structure interaction experiments and analyses." Journal of Geotechnical and
Geoenvironmental Engineering, ASCE, 125(9): 750-759. Modifications were required for
Geoenvironmental Engineering, ASCE, 125(9): 750-759. Modifications were required for
representing the different responses of a q-z material in compression versus uplift.
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
The nonlinear q-z behavior is conceptualized as consisting of elastic (q-<math>z^e</math>), plastic (q-<math>z^p</math>), and
gap (q-zg) components in series. Radiation damping is modeled by a dashpot on the “far-field”
gap (<math>q-z^g</math>) components in series. Radiation damping is modeled by a dashpot on the “far-field”
elastic component (q-ze) of the displacement rate. The gap component consists of a bilinear
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 z = ze + zp + zg,
closure spring (<math>q^c-z^g</math>) in parallel with a nonlinear drag spring (<math>q^d-z^g</math>). Note that <math>z = z^e + z^p + z^g</math>, and that <math>q = q^d + q^c</math>.
and that q = qd + qc.
 
The plastic component has an initial range of rigid behavior between -Cr qult < q < Cr qult with
The plastic component has an initial range of rigid behavior between <math>-C_r q_\text{ult} < q < C_r q_\text{ult}</math> with
Cr = the ratio of q/qult when plastic yielding first occurs in virgin loading. The rigid range of q,
<math>C_r</math> = 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 2 Cr qult, translates and grows with plastic yielding. The rigid range of q is
which is initially <math>2 C_r q_\text{ult}</math>, translates and grows with plastic yielding. The rigid range of q is
constrained to a maximum size of 0.7qult. Beyond the rigid range, loading of the plastic (q-zp)
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:
component is described by:


[[Image:TzSimple1.gif]]
:<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., <math>C_d</math>) 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.
 
[[Image:QzSimple1A.gif]]
 
[[Image:QzSimple1B.gif]]


EXAMPLE:
EXAMPLE:

Latest revision as of 21:59, 22 March 2010




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.


EQUATIONS and EXAMPLE RESPONSES:

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 (<math>q^c-z^g</math>) in parallel with a nonlinear drag spring (<math>q^d-z^g</math>). 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>-C_r q_\text{ult} < q < C_r q_\text{ult}</math> with <math>C_r</math> = 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 C_r 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., <math>C_d</math>) 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: