rudnicki86: Rudnicki's response to point and step injection

Description Usage Arguments References Examples

View source: R/rudnicki86.R

Description

Rudnicki's response to point and step injection

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
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 Feb. 15, 2022, 6:24 p.m.