calc_cr_fbmw: Calculate root reinforcement at current strain according to...

In GJMeijer/FBMcw: Generic fibre bundle model for soil reinforcement by plant roots

Description

This function calculates the current root reinforcement according to the FBMw, assuming discrete root classes, given the current level of (normalised) strain in the reference root

Usage

calc_cr_fbmw(
  epsr0rel,
  dr,
  phir,
  tru,
  betaF,
  tru0,
  kappa,
  k = 1.2,
  dr0 = 1,
  sumoutput = T
)

Arguments

epsr0rel	tensile strain in reference root, normalised by the tensile strain to peak in this root (numeric array, size n)
dr	representative diameter of root class (numeric arra, size m)
phir	root area ratio of roots in class (numeric arra, size m)
tru	tensile strength of roots in class (numeric array, size m)
betaF	load sharing parameter (numeric scalar)
tru0	Tensile strength of root with reference diameter (numeric array)
kappa	Weibull shape parameter required for the root survival function (scalar)
k	Wu/Waldron factor accounting for root orientations at failure (numeric scalar, default 1.2)
dr0	Reference diameter (numeric scalar, default 1.0)
sumoutput	If TRUE, a total reinforcement is outputted for every strain step. If FALSE, a reinforcement matrix is outputted (rows: strain step, columns: classes)

Value

numeric array with root reinforcement preductions 'cr_fbmw'

Examples

calc_cr_fbmw(seq(0,2.5,l=101), seq(1,4), rep(0.01,4), 10e3*seq(4)^-0.5, 1, 10e3, 4)

GJMeijer/FBMcw documentation built on Dec. 17, 2021, 9:23 p.m.