CoulombStress: Calculate Coulomb stresses on a plane from a given stress...

In abarbour/strain: Tools for working with strainmeter and NOTA borehole data

Calculate Coulomb stresses on a plane from a given stress tensor

Description

Calculate Coulomb stresses on a plane from a given stress tensor

Usage

CoulombStress(S, n = NULL, mu = 0.6, B = 0.85, in.bar = FALSE, unit.scale = 1)

projection_matrix(n = c(1, 1, 1))

Arguments

S	symmetric stress tensor, in Pascals
n	an optional vector normal to the plane on which to project the stresses... not used
mu	coefficient of internal friction; Byerlee puts this at 0.6 for deep rock
B	Skempton's coefficient for an undrained material; 0.85 is a reasonable number for bulk Westerly granite
in.bar	logical; should the stresses be in units of bar rather than Pa?
unit.scale	numeric; stresses will be multiplied by this value. Useful to return, say, gigaPascals (use unit.scale=1e-9 then); but, note that because Pore pressure is derived from volume stress, that too will be scaled by this value.

See Also

stresses

abarbour/strain documentation built on Dec. 10, 2024, 10:58 p.m.