Newmark displacement curve Dn(ky) for one scenario

Evaluates ensemble and per-submodel Dn over a vector of ky values. The scenario is sampled once and the same epsilon is reused across all ky so that each realisation s represents a coherent physical curve Dn_s(ky).

Usage

fitDnCurve(uhs, ky, Ts, Mw = 6.5, NS = 30L, Rrup = 100, weights, NFC = "D100")

Arguments

uhs

data.table with columns Tn, p, Sa (one scenario)

ky

numeric vector, yield accelerations (g)

Ts

numeric scalar, fundamental period (s)

Mw

numeric scalar, moment magnitude

NS

integer, Monte Carlo samples for curve quantiles and per-realisation draws.

Rrup

numeric scalar, rupture distance (km)

weights

named numeric vector, ensemble weights by IDn

NFC

character, near-fault component selector for model BM19: "D100" (maximum component, default) or "D50". Ignored by all other submodels.

Value

list(curve, draws). curve = data.table(ky, p, Dn, IDn, w). draws = data.table(ky, s, Dn, IDn) — always populated (possibly empty when no model produces output). Dn in cm.

Examples

if (FALSE) { # \dontrun{
out <- fitDnCurve(uhs, ky = getDnKy(uhs, Ts = 0.3), Ts = 0.3, Mw = 7.5,
                  NS = 100, weights = c(AM88=1, BT07=1, SR08=1, BM17=1))
} # }