Fundamental periods of the soil profile and truncated shear beam for every realization

Applies solveWedgePeriod() and solveColumnPeriod() to each layer table of a list, with the paired height, truncation ratio, exponent and basal velocity of the same realization, so that the numerical periods and their controls stay paired with the realization before any summary.

Usage

solveProfilePeriods(layers, Hs, lo, mo, VSo, method, elements)

Arguments

layers: A non-empty list of layer tables, such as the layers column of simulateSoilProfiles(); see solveColumnPeriod().
Hs: Slope heights in metres: one for every table, or one per table.
lo: Truncation ratios in \([0, 1)\): one for every table, or one per table.
mo: Inhomogeneity exponents of the analytic route, one per table.
VSo: Basal shear-wave velocities of the analytic route in m/s, one per table.
method: The solution method of the root model: "exact", or the surrogate "lm" or "nlm"; see buildRootModel().
elements: One positive whole number of linear elements per layer of the declared mesh; the control doubles it.

Value

A data.table with one row per layer table, in the order of the list, and the columns of solveWedgePeriod() followed by those of solveColumnPeriod(). An invalid table stops the call with its position in the list.

Examples

set.seed(20260915)
Profiles <- simulateSoilProfiles(Hs = 10, USCS = c("ML", "MH"), NS = 2, h = 1,
                                 water = 0, POP = 100, Vref = 760)
solveProfilePeriods(Profiles$layers, Hs = 10, lo = 0.351, mo = Profiles$mo,
                    VSo = Profiles$VSo, method = "exact", elements = 2)
#>           Ts TsMeshChange TsAnalytic TsAnalyticDifference       an anExact
#>        <num>        <num>      <num>                <num>    <num>  <lgcl>
#> 1: 0.2045104 0.0002263781  0.1922096          -0.06014783 3.323958    TRUE
#> 2: 0.1951046 0.0003668317  0.1808622          -0.07299876 3.768778    TRUE
#>    anWithinSupport        Tr TrMeshChange TrRayleighWeighted
#>             <lgcl>     <num>        <num>              <num>
#> 1:            TRUE 0.2470676 0.0002935533          0.2443280
#> 2:            TRUE 0.2355099 0.0004297911          0.2327509
#>    TrRayleighWeightedDifference
#>                           <num>
#> 1:                  -0.01108829
#> 2:                  -0.01171491