geophys (version 1.2-1)

Maxstress: Maximum Stress

Description

Maximum Normal and Shear stress in a plane

Usage

Maxstress(NN, Stensor)

Arguments

NN

Normal Vector

Stensor

Stress tensor

Value

• NN
  Normal vector to plane is returned
• sigNORMmax
  maximum normal stress
• tauSHEARmax
  maximum shear stress in the plane

Details

Stress tensor can be entered as a vector of 3 or six values or as a 3 by 3 matrix. If three values are provided provided they are assumed to be the principal stresses. If six values are provided, they are c(x11, x12, x13, x22, x33). In either case eigen values and eigen vectors are calculated and used in the following analysis.

See Also

stress

Examples

S = stressSETUP( )

Nvec = NORMvec(S$PPs, S$xscale, S$Rview, S$aglyph , add = FALSE)

Stensor = matrix(c(
15, 0, 0,
0, 8, 0,
0,  0, 5), ncol=3)

 

Mstress  = Maxstress(Nvec, Stensor)

DoMohr(Stensor)
 axis(1)
axis(2)


points(Mstress$sigNORMmax , Mstress$tauSHEARmax, pch=21, col='blue'  , bg='gold' )

u=par('usr')

segments(0, Mstress$tauSHEARmax, Mstress$sigNORMmax ,
Mstress$tauSHEARmax, lty=2, col='green'  )

text(mean(c(0, Mstress$tauSHEARmax)),  Mstress$tauSHEARmax,
"MaxShear in Plane", pos=3)


segments(Mstress$sigNORMmax , u[3] , Mstress$sigNORMmax ,
Mstress$tauSHEARmax, lty=2, col='purple'  )

text(Mstress$sigNORMmax , u[3], "MaxNormal stress", adj=c(0,-1) )