SSPquadUP Element: Difference between revisions

From OpenSeesWiki
Jump to navigation Jump to search
(Created page with '{{CommandManualMenu}} This command is used to construct a SSPquadUP element object. {| | style="background:yellow; color:black; width:800px" | '''element SSPquadUP $eleTag $iN...')
 
No edit summary
 
(15 intermediate revisions by the same user not shown)
Line 11: Line 11:
|  style="width:150px" | '''$eleTag ''' || unique integer tag identifying element object
|  style="width:150px" | '''$eleTag ''' || unique integer tag identifying element object
|-
|-
|  '''$iNode $jNode $kNode $lNode''' || the four nodes defining the element, input in counterclockwise order
|  '''$iNode $jNode $kNode $lNode''' || the four nodes defining the element, input in counterclockwise order (-ndm 2 -ndf 3)
|-
|-
|  '''$matTag''' || unique integer tag associated with previously-defined nDMaterial object
|  '''$matTag''' || unique integer tag associated with previously-defined nDMaterial object
Line 21: Line 21:
|  '''$fDen''' || mass density of the pore fluid
|  '''$fDen''' || mass density of the pore fluid
|-
|-
|  '''$k1 $k2''' || permeability coefficients in horizontal and vertical directions, respectively
|  '''$k1 $k2''' || permeability coefficients in global x- and y-directions, respectively
|-
|-
|  '''$void''' || voids ratio
|  '''$void''' || voids ratio
|-
|-
|  '''$alpha''' || spatial pressure field stabilization parameter (see notes below for more information)
|  '''$alpha''' || spatial pressure field stabilization parameter (see discussion below for more information)
|-
|-
|  '''$b1 $b2''' || constant body forces in horizontal and vertical directions, respectively (optional, default = 0.0)
|  '''$b1 $b2''' || constant body forces in global x- and y-directions, respectively (optional, default = 0.0) - '''See Note 3'''
|}
|}


Line 34: Line 34:




The SSPquadUP element is an extension of the SSPquad element for use in dynamic plane strain analysis of fluid saturated porous media.  A mixed displacement-pressure (u-p) formulation is used, based upon the work of Biot as extended by Zienkiewicz and Shiomi (1984).
The SSPquadUP element is an extension of the [[SSPquad Element]] for use in dynamic plane strain analysis of fluid saturated porous media.  A mixed displacement-pressure (u-p) formulation is used, based upon the work of Biot as extended by Zienkiewicz and Shiomi (1984).


The physical stabilization necessary to allow for reduced integration incorporates an assumed strain field in which the volumetric dilation and the shear strain associated with the the hourglass modes are zero, resulting in an element which is free from volumetric and shear locking.  The elimination of shear locking results in greater coarse mesh accuracy in bending dominated problems, and the elimination of volumetric locking improves accuracy in nearly-incompressible problems.  Analysis times are generally faster than corresponding full integration elements.   
The physical stabilization necessary to allow for reduced integration incorporates an assumed strain field in which the volumetric dilation and the shear strain associated with the the hourglass modes are zero, resulting in an element which is free from volumetric and shear locking.  The elimination of shear locking results in greater coarse mesh accuracy in bending dominated problems, and the elimination of volumetric locking improves accuracy in nearly-incompressible problems.  Analysis times are generally faster than corresponding full integration elements.   


Equal-order interpolation is used for the displacement and pressure fields, thus, the SSPquadUP element does not inherently pass the ''inf-sup'' condition in the incompressible-impermeable limit.  A stabilizing parameter is employed to permit the use of equal-order interpolation for this element.  This parameter '''$alpha''' can be computed as
Equal-order interpolation is used for the displacement and pressure fields, thus, the SSPquadUP element does not inherently pass the ''inf-sup'' condition, and is not fully acceptable in the incompressible-impermeable limit (the [http://opensees.berkeley.edu/OpenSees/manuals/usermanual/1098.htm QuadUP Element] has the same issue).  A stabilizing parameter is employed to permit the use of equal-order interpolation for the SSPquadUP element.  This parameter '''$alpha''' can be computed as
 
alpha = 0.25*(h^2)/(den*c^2)
 
where '''h''' is the element size, '''c''' is the speed of elastic wave propagation in the solid phase, and '''den''' is the mass density of the solid phase.  The '''$alpha''' parameter should be a small number.  With a properly defined '''$alpha''' parameter, the SSPquadUP element can produce comparable results to a higher-order element such as the [[Nine Four Node Quad u-p Element | 9_4_QuadUP Element]] at a significantly lower computational cost and with a greater ease in mesh generation.
 
The full formulation for the SSPquadUP element can be found in McGann et al. (2012) along with several example applications.


:::<math> \alpha = \alpha_0^{}\frac{\Delta t_{\mathrm{crit}}^2}{\rho}
</math>


'''NOTES:'''
'''NOTES:'''
   
 
# Valid queries to the SSPquad element when creating an ElementalRecorder object correspond to those for the nDMaterial object assigned to the element (e.g., 'stress', 'strain').  Material response is recorded at the single integration point located in the center of the element.   
# The SSPquadUP element will only work in dynamic analysis.
# The SSPquad element was designed with intentions of duplicating the functionality of the Quad element.  If an example is found where the SSPquad element cannot do something that works for the Quad element, please contact the developers listed below so the bug can be fixed.
# For saturated soils, the mass density input into the associated nDMaterial object should be the saturated mass density. 
# When modeling soil, the body forces input into the SSPquadUP element should be the components of the gravitational vector, not the unit weight.
# Fixing the pore pressure degree-of-freedom (dof 3) at a node is a drainage boundary condition at which zero pore pressure will be maintained throughout the analysis. Leaving the third dof free allows pore pressures to build at that node.
# Valid queries to the SSPquadUP element when creating an ElementalRecorder object correspond to those for the nDMaterial object assigned to the element (e.g., 'stress', 'strain').  Material response is recorded at the single integration point located in the center of the element.   
# The SSPquadUP element was designed with intentions of duplicating the functionality of the [http://opensees.berkeley.edu/OpenSees/manuals/usermanual/1098.htm QuadUP Element].  If an example is found where the SSPquadUP element cannot do something that works for the [http://opensees.berkeley.edu/OpenSees/manuals/usermanual/1098.htm QuadUP Element], e.g., material updating, please contact the developers listed below so the bug can be fixed.


'''EXAMPLES:'''
'''EXAMPLES:'''


SSPquad element definition with element tag 1, nodes 1, 2, 3, and 4, material tag 1, plane strain conditions, unit thickness, horizontal body force of zero, and vertical body force of -10.0
SSPquadUP element definition with element tag 1, nodes 1, 2, 3, and 4, material tag 1, unit thickness, bulk modulus of water (kPa), mass density of water (Mg/m^3), horizontal and vertical permeabilities of 1e-3, voids ratio of 0.7, alpha parameter of 6e-5, horizontal body force of zero, and vertical body force of -9.81


  element SSPquad 1  1 2 3 4  1 "PlaneStrain" 1.0 0.0 -10.0
  element SSPquadUP 1  1 2 3 4  1 1.0 2.2e6 1.0 1.0e-3 1.0e-3 0.7 6.0e-5 0.0 -9.81


Elemental recorders for stress and strain when using the SSPquad element (note the difference from the Quad element)
Elemental recorders for stress and strain when using the SSPquadUP element (note the difference from the [http://opensees.berkeley.edu/OpenSees/manuals/usermanual/1098.htm QuadUP Element])


  recorder Element -eleRange 1 $numElem -time -file stress.out  stress
  recorder Element -eleRange 1 $numElem -time -file stress.out  stress
  recorder Element -eleRange 1 $numElem -time -file strain.out  strain
  recorder Element -eleRange 1 $numElem -time -file strain.out  strain
Pore pressure recorder for the SSPquadUP element (pore pressure is the third degree-of-freedom)
recorder Node -nodeRange 1 $numNode -time -file porePressure.out -dof 3 vel
'''REFERENCES:'''
McGann, C. R., Arduino, P., and Mackenzie-Helnwein, P. (2012).  “Stabilized single-point 4-node quadrilateral element for dynamic analysis of fluid saturated porous media.”  ''Acta Geotechnica'', 7(4), 297-311.
Zienkiewicz, O.C. and Shiomi, T. (1984). "Dynamic behavior of saturated porous media; the generalized Biot formulation and its numerical solution." ''International Journal for Numerical Methods in Geomechanics'', 8, 71-96.





Latest revision as of 22:03, 11 January 2017




This command is used to construct a SSPquadUP element object.

element SSPquadUP $eleTag $iNode $jNode $kNode $lNode $matTag $thick $fBulk $fDen $k1 $k2 $void $alpha <$b1 $b2>


$eleTag unique integer tag identifying element object
$iNode $jNode $kNode $lNode the four nodes defining the element, input in counterclockwise order (-ndm 2 -ndf 3)
$matTag unique integer tag associated with previously-defined nDMaterial object
$thick thickness of the element in out-of-plane direction
$fBulk bulk modulus of the pore fluid
$fDen mass density of the pore fluid
$k1 $k2 permeability coefficients in global x- and y-directions, respectively
$void voids ratio
$alpha spatial pressure field stabilization parameter (see discussion below for more information)
$b1 $b2 constant body forces in global x- and y-directions, respectively (optional, default = 0.0) - See Note 3




The SSPquadUP element is an extension of the SSPquad Element for use in dynamic plane strain analysis of fluid saturated porous media. A mixed displacement-pressure (u-p) formulation is used, based upon the work of Biot as extended by Zienkiewicz and Shiomi (1984).

The physical stabilization necessary to allow for reduced integration incorporates an assumed strain field in which the volumetric dilation and the shear strain associated with the the hourglass modes are zero, resulting in an element which is free from volumetric and shear locking. The elimination of shear locking results in greater coarse mesh accuracy in bending dominated problems, and the elimination of volumetric locking improves accuracy in nearly-incompressible problems. Analysis times are generally faster than corresponding full integration elements.

Equal-order interpolation is used for the displacement and pressure fields, thus, the SSPquadUP element does not inherently pass the inf-sup condition, and is not fully acceptable in the incompressible-impermeable limit (the QuadUP Element has the same issue). A stabilizing parameter is employed to permit the use of equal-order interpolation for the SSPquadUP element. This parameter $alpha can be computed as

alpha = 0.25*(h^2)/(den*c^2)

where h is the element size, c is the speed of elastic wave propagation in the solid phase, and den is the mass density of the solid phase. The $alpha parameter should be a small number. With a properly defined $alpha parameter, the SSPquadUP element can produce comparable results to a higher-order element such as the 9_4_QuadUP Element at a significantly lower computational cost and with a greater ease in mesh generation.

The full formulation for the SSPquadUP element can be found in McGann et al. (2012) along with several example applications.


NOTES:

  1. The SSPquadUP element will only work in dynamic analysis.
  2. For saturated soils, the mass density input into the associated nDMaterial object should be the saturated mass density.
  3. When modeling soil, the body forces input into the SSPquadUP element should be the components of the gravitational vector, not the unit weight.
  4. Fixing the pore pressure degree-of-freedom (dof 3) at a node is a drainage boundary condition at which zero pore pressure will be maintained throughout the analysis. Leaving the third dof free allows pore pressures to build at that node.
  5. Valid queries to the SSPquadUP element when creating an ElementalRecorder object correspond to those for the nDMaterial object assigned to the element (e.g., 'stress', 'strain'). Material response is recorded at the single integration point located in the center of the element.
  6. The SSPquadUP element was designed with intentions of duplicating the functionality of the QuadUP Element. If an example is found where the SSPquadUP element cannot do something that works for the QuadUP Element, e.g., material updating, please contact the developers listed below so the bug can be fixed.

EXAMPLES:

SSPquadUP element definition with element tag 1, nodes 1, 2, 3, and 4, material tag 1, unit thickness, bulk modulus of water (kPa), mass density of water (Mg/m^3), horizontal and vertical permeabilities of 1e-3, voids ratio of 0.7, alpha parameter of 6e-5, horizontal body force of zero, and vertical body force of -9.81

element SSPquadUP 1  1 2 3 4  1 1.0 2.2e6 1.0 1.0e-3 1.0e-3 0.7 6.0e-5 0.0 -9.81

Elemental recorders for stress and strain when using the SSPquadUP element (note the difference from the QuadUP Element)

recorder Element -eleRange 1 $numElem -time -file stress.out  stress
recorder Element -eleRange 1 $numElem -time -file strain.out  strain

Pore pressure recorder for the SSPquadUP element (pore pressure is the third degree-of-freedom)

recorder Node -nodeRange 1 $numNode -time -file porePressure.out -dof 3 vel

REFERENCES:

McGann, C. R., Arduino, P., and Mackenzie-Helnwein, P. (2012). “Stabilized single-point 4-node quadrilateral element for dynamic analysis of fluid saturated porous media.” Acta Geotechnica, 7(4), 297-311.

Zienkiewicz, O.C. and Shiomi, T. (1984). "Dynamic behavior of saturated porous media; the generalized Biot formulation and its numerical solution." International Journal for Numerical Methods in Geomechanics, 8, 71-96.



Code Developed by: Chris McGann, Pedro Arduino, & Peter Mackenzie-Helnwein, at the University of Washington