seismic_elasticity: Calculate elastic properties from a velocity model

In abarbour/stress: Tools for Working with Stresses and Focal-mechanism Inversions

Description

Calculate elastic properties from a velocity model

Usage

seismic_elasticity(vp, vs, rho)

Arguments

vp	numeric; the P-wave speed [m/s]
vs	numeric; the S-wave speed [m/s]
rho	numeric; the density [kg/m^3]
...	additional parameters (e.g., passing nu)

Value

bulk modulus, shear modulus, Lame's first constant; in Pa if SI given for inputs

Examples

# inputs in m/s and kg/m^3 give results in Pascal
Pa <- seismic_elasticity(7200, 3500, 2650)

# use km/s and g/cm^3 instead to get giga-Pascal
GPa <- seismic_elasticity(7.2, 3.5, 2.65)
all.equal(Pa/1e9, GPa)

abarbour/stress documentation built on Oct. 5, 2019, 11:20 a.m.