Skip to contents

Builds one Monte Carlo realisation of a layered soil column from a total height Hs and a USCS classification, samples random geotechnical properties (void ratio, unit weight, plasticity, shear modulus parameters), computes the fundamental period Ts via Rayleigh's method (fitModel.Ts) and the Vs30-equivalent shear-wave velocity. getSiteProperties calls this function repeatedly to assemble the Monte Carlo ensemble; call directly when one realisation suffices.

Usage

geSiteTable(
  Hs,
  Water = 0,
  USCS,
  Group = NULL,
  h = 0.5,
  DrID = NULL,
  Vref = 760,
  UniformDistribution = TRUE,
  POP = 100,
  IgnoreModelIntervals = TRUE,
  getSiteLayers = FALSE
)

Arguments

Hs

Numeric scalar. Total height to hard ground (m).

Water

Numeric scalar. Water-table depth as fraction of Hs: zw = Hs - Water * Hs.

USCS

Character vector. Unified Soil Classification System codes (e.g. c("GC", "CL", "ML")).

Group

Character vector or NULL. Soil groups, c("Gravels","Fines","Sands"); if NULL, derived from USCS.

h

Numeric scalar. Layer thickness in the discretisation (m). Default 0.50.

DrID

Character or NULL. Relative-density label, c("Very Dense","Dense","Compact","Loose","Very Loose"); NULL samples relative density from the standard prior.

Vref

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

UniformDistribution

Logical. If TRUE (default), sample soil properties from uniform distributions; if FALSE, use normal distributions.

POP

Numeric. Pre-consolidation pressure (kPa); enters the over-consolidation ratio OCR = (pm + POP) / pm. Default 100.

IgnoreModelIntervals

Logical. If TRUE (default), ignore model-specific layer-thickness intervals and use h as the discretisation step.

getSiteLayers

Logical. If TRUE, return the per-layer detail in addition to the per-site summary.

Value

data.table with one row per (Hs, POP, Hw, sample) combination and columns Hs, Hw, NL, Z500, Z1000, SID, Go, mo, Ts, VSo, VS30, plus USCS-fraction columns.