Skip to content

Private Interface

Index

Interfaces

# JuEQ.HomogeneousElasticProperties — Type.

Okada's dc3d only applies on isotropic materials, therefore, elastic modulus are constrained to be scalars.

source

# JuEQ.ODEStateVariable — Type.

Intermediate variable in solving ODEs aimed to avoid allocation overheads.

# JuEQ.applied_unit_dislocation — Method.

For noraml fault, it should of course be [0., -1., 0.]. However, in term of force balance, it is quivalent to thrust fault if dip angle are constrained within [0, π/2] in fact.

The unit of unit dislocation below is the same of v * t at set by user so to avoid normalization step.

# JuEQ.shear_traction — Method.

1
shear_traction(::Type{<:PlaneFault}, u, λ, μ, dip)

Calculate the shear traction on the fault plane w.r.t. fault types.

Arguments

  • u::AbstractArray{<:Number, 1}: the output from dc3d_okada
  • λ::Number: Lamé's first parameter
  • μ::Number: shear modulus
  • dip::Number: plane dip angle

Reference

# JuEQ.stiffness_periodic_boundary_condition! — Method.

Periodic boundary condition for 2D faults.

Arguments

  • u::AbstractVector: In-place output which is a 12-elements vector (exactly the output of dc3d_okada). No assertion here imposed.
  • same as dc3d_okada, see dc3d for details.
  • nrept::Integer: (half) number of repetition, as denoted by -npret: nrept
  • lrept::Number: length of repetition interval, see Note below

Note

  • The buffer block length is (buffer_ratio - 1) multipled by along-strike length.