Monte Carlo site properties from synthetic soil profiles

Generates NR realisations of synthetic soil profiles for a given total height Hs and USCS classification, then summarises the resulting Ts, mo, Go, VSo and VS30 by user-requested quantile levels (and/or the mean). Each realisation is built by geSiteTable().

Usage

getSiteProperties(
  Hs,
  USCS,
  POP = 100,
  Water = 0,
  NR = 1,
  h = 1,
  levels = c(0.05, 0.5, "mean", 0.95),
  Vref = 760
)

Arguments

Hs

Numeric vector. Total height(s) to hard ground (m).

USCS

Character vector. USCS soil classification codes (e.g. c("GC", "CL", "ML")).

POP

Numeric. Pre-consolidation pressure (kPa). Default 100.

Water

Numeric. Water-table depth as fraction of Hs. Default 0 (no water).

NR

Integer. Number of Monte Carlo realisations. Default 1.

h

Numeric. Layer thickness (m) in the discretisation. Default 1.00.

levels

Vector. Quantile levels to report; entries can be numeric in (0, 1) and/or the literal "mean". Default c(0.05, 0.5, "mean", 0.95).

Vref

Numeric. Reference shear-wave velocity (m/s) for site-class assignment. Default 760.

Value

data.table with one row per (Hs, POP, Hw, level) combination and columns USCS, Go, mo, Ts, VSo, VS30, Z500, Z1000, level. Numeric properties are rounded for readability.

Examples

# \donttest{
getSiteProperties(
  Hs   = 30,
  USCS = c("GC", "CL", "ML"),
  NR   = 50,
  levels = c(0.16, "mean", 0.84)
)
#>       Hs   POP    Hw     USCS    Go    mo    Ts   VSo  VS30  Z500  Z1000  level
#>    <num> <num> <num>   <char> <num> <num> <num> <num> <num> <num>  <num> <char>
#> 1:    30   100     0 GC CL ML 145.0 0.419  0.64 268.5 203.9   2.1  9.400   0.16
#> 2:    30   100     0 GC CL ML 188.3 0.575  0.71 305.5 221.2   3.3 21.300   0.84
#> 3:    30   100     0 GC CL ML 165.6 0.505  0.68 287.7 212.0   2.8 16.238   mean
# }