stressPlot: Plot of Stress Profile

In geotech: Geotechnical Engineering

Description

These functions produce plots of vertical stresses versus depth (plot.sigmaZ) or horizontal stresses versus depth (plot.sigmaX).

Usage

sigmaX.plot(depth, sigmaX.eff, sigmaX.total = NA, u = NA, metric)
sigmaZ.plot(depth, sigmaZ.eff, sigmaZ.total = NA, u = NA, metric)

Arguments

depth	vector of depths (ft or m)
sigmaX.eff	vector of effective horizontal stresses (psf or kPa)
sigmaX.total	vector of total horizontal stresses (psf or kPa)
sigmaZ.eff	vector of effective vertical stresses (psf or kPa)
sigmaZ.total	vector of total vertical stresses (psf or kPa)
u	vector of effective vertical stresses (psf or kPa)
metric	logical variable: TRUE (for metric units) or FALSE (for English units)vector of total vertical stresses (psf or kPa)

Details

• Arguments sigmaX.eff and sigmaX.total are used in conjunction with sigmaX.plot, and arguments sigmaZ.eff and sigmaZ.total are used in conjunction with sigmaZ.plot.

• If total stresses and pore water pressure are left blank, the plot is only constructed for effective stress.

• Once constructed, additional profiles may be added by the user to this plot (for example, for induced stress or maximum past pressure).

Value

This function creates a plot of stresses; no numerical output is obtained. To perform numerical stress calculations, see stressHorizontal or stressVertical.

Author(s)

James Kaklamanos <[e-mail]> and Kyle Elmy <[e-mail]>

See Also

stressHorizontal, stressVertical

Examples

##  Example code for vertical stress plot
temp <- sigmaZ.profile(gamma = rep(100, 3), depth = c(10, 20, 30),
                       zw = 10, metric = FALSE)
depth <- temp$depth
sigmaTotal <- temp$sigmaZ.total
u <- temp$u
sigmaEff <- temp$sigmaZ.eff
sigmaZ.plot(depth = depth, sigmaZ.eff = sigmaEff, metric = FALSE,
            sigmaZ.total = sigmaTotal, u = u)

##  Example code for horizontal stress plot
##  Site with constant unit weight = 100 pcf, GWT at 10 ft depth
temp <- sigmaX.profile(gamma = rep(100, 3), depth = c(10, 20, 30),
                       K = c(0.35, 0.30, 0.28), zw = 10, metric = FALSE)
depth <- temp$depth
sigmaTotal <- temp$sigmaX.total
u <- temp$u
sigmaEff <- temp$sigmaX.eff
sigmaX.plot(depth = depth, sigmaX.eff = sigmaEff, metric = FALSE,
            sigmaX.total = sigmaTotal, u = u)

geotech documentation built on May 1, 2019, 9:09 p.m.