Vertical Stress Calculations

Share:

Description

These functions calculate vertical stress at a point (sigmaZ) and versus depth (sigmaZ.profile)

Usage

1
2
3
sigmaZ(gamma, thk = NA, depth = NA, zw, zout, gammaW = NA, metric)
sigmaZ.profile(gamma, thk = NA, depth = NA, zw, zout = NA,
               gammaW = NA, metric)

Arguments

gamma

vector of unit weights (pcf or kN/m^3)

thk

vector of layer thicknesses (ft or m)

depth

vector of layer bottom depths (ft or m)

zw

depth of groundwater table (ft or m)

zout

desired depth of output (ft or m): a single value

gammaW

unit weight of water (default = 62.4 pcf for English units; 9.81 kN/m^3 for metric units)

metric

logical variable: TRUE (for metric units) or FALSE (for English units)

Details

Either layer thicknesses or depths to layer bottoms must be specified.

Value

Function sigmaZ outputs a three-element list giving the stresses at a specified depth:

  • sigmaZ.eff = effective vertical stress

  • sigmaZ.total = total vertical stress

  • u = pore water pressure

Function sigmaZ.profile outputs a four-element list giving the stress variation with depth (four vectors):

  • depth = depth

  • sigmaZ.eff = effective vertical stress

  • sigmaZ.total = total vertical stress

  • u = pore water pressure

Author(s)

James Kaklamanos <kaklamanosj@merrimack.edu> and Kyle Elmy <ElmyK@merrimack.edu>

See Also

stressHorizontal, stressInducedArea, stressInducedPoint, stressPlot

Examples

1
2
3
4
5
## Example code for Vertical Stress at a Point
sigmaZ(gamma = c(108, 116), depth = c(15, 40), zout = 18, zw = 15, metric = FALSE)

##  Example code for Vertical Stress Profile
sigmaZ.profile(gamma = c(108, 116), depth = c(15, 40), zw = 15, metric = FALSE)