vertical_stress_profile: Calculate verticals stresses using layers

In GJMeijer/soilmech: R functions and plots for Soil Mechanics teaching

Calculate verticals stresses using layers

Description

Function takes a soil profile and returns the vertical stresses at all layer interfaces

Usage

vertical_stress_profile(
  z_soil = 0,
  z_watertable = 0,
  z_max = 10,
  z_interval = NULL,
  gamma_b = 19,
  gamma_w = 10,
  q = 0,
  ...
)

Arguments

z_soil	depth of top of each layer
z_watertable	depth of the water table
z_max	maximum depth of interest
z_interval	additional interval at which to calculate stresses
gamma_b	unit weight of soil, in kN/m3. Should be a scalar or a vector with the same length as 'z_soil'
gamma_w	unit weight of water: 10 kN/m3
q	surcharge at soil surface
...	additional fields to be added to the returned data

Value

a tibble with depth 'z', vertical total stress 'sigma_v', vertical effective stress ‘sigma’_v' and pore pressure 'u'. All fields in '...' are added too

Examples

#only calculate stresses at layer interfaces and water table
vertical_stress_profile(
  c(0, 4, 6),
  z_watertable = 2,
  gamma_b = c(18, 20, 22),
  q = 10
)

#add stresses at regular intervals
vertical_stress_profile(
  c(0, 4, 6),
  z_watertable = 2,
  gamma_b = c(18, 20, 22),
  q = 10,
  z_interval = 1
)

GJMeijer/soilmech documentation built on May 22, 2022, 10:39 a.m.