rudnicki86: Rudnicki's response to point and step injection

In abarbour/deform: Crustal Deformation Modeling Tools in R

Rudnicki's response to point and step injection

Description

Rudnicki's response to point and step injection

Usage

rudnicki86(
  r,
  z,
  t,
  zinj = 0,
  mu = 1,
  diffusiv = 1,
  nu = 0.25,
  nuu = 0.33,
  B = 1,
  response = NULL
)

## S3 method for class 'rudnicki.pt'
plot(x, response = FALSE, ...)

Arguments

r, z	numeric; receiver position (radial, depth-positive) [m]
t	numeric; time [s]
zinj	numeric; injection depth [m]
mu	numeric; the isotropic elastic shear modulus (keep at unity for 'normalized')
diffusiv	numeric; the isotropic hydraulic diffusivity (keep at unity for 'normalized')
nu	numeric; the drained Poisson's ratio [0,0.5]
nuu	numeric; the undrained Poisson's ratio [0,0.5] (nu < nuu)
B	numeric; Skempton's coefficient [0,1]
response	character; the type of response to output: either the points-source impulsive ('impulse') or Heaviside ('step') response
x	an object with class 'rudnicki.pt', e.g. from rudnicki86
...	additional arguments

References

Rudnicki, J. W. (1986), Fluid mass sources and point forces in linear elastic diffusive solids, Mechanics of Materials, 5(4), 383-393, https://doi.org/10.1016/0167-6636(86)90042-6

Examples

r <- 1
zi <- 1
t <- 10^seq(-4, 1, length.out=301)

# Impulse response
rt <- rudnicki86(r,1,t,zi)
rtdeep <- rudnicki86(r,10,t,zi)

# Step response
rt_s <- rudnicki86(r,1,t,zi, response='step')
rtdeep_s <- rudnicki86(r,10,t,zi, response='step')

layout(matrix(1:10, 2, byrow=TRUE))
plot(rt, col=1)
plot(rt, response=TRUE, col=1)

plot(rtdeep, col=2)
plot(rtdeep, response=TRUE, col=2)

plot(rt_s, col=3)
plot(rt_s, response=TRUE, col=3)

plot(rtdeep_s, col=4)
plot(rtdeep_s, response=TRUE, col=4)

abarbour/deform documentation built on Dec. 5, 2024, 1:55 p.m.