Civil Engineering Reference
In-Depth Information
IF(numpe==1)THEN
write(11,'(A,E12.4)')"The displacement is ",totd_pp(1)
write(11,'(A)') " sigma z sigma x sigma y "
write(11,'(3E12.4)')tensor_pp(3,1,1),tensor_pp(1,1,1),tensor_pp(2,1,1)
write(11,'(A,I12)')"The total number of cj iterations was ",cjtot
write(11,'(A,I12)')"The number of plastic iterations was ",plasiters
write(11,'(A,F11.2)')"cj iterations per plastic iteration were
",
&
& REAL(cjtot)/REAL(plasiters)
END IF
IF (plasiters==plasits) EXIT
END DO load_increments
IF(numpe==1) WRITE(11,*) "This analysis took : ", elap_time( ) - timest(1)
CALL shutdown( )
END PROGRAM p122
New scalar integers:
cjiters
conjugate gradient iteration counter
conjugate gradient iteration ceiling
cjits
total number of conjugate gradient iterations
cjtot
number of load increments
incs
simple counter
iy
plastic iteration counter
plasiters
plastic iteration ceiling
plasits
New scalar reals:
big
largest component of a vector
cohesion
c
conjugate gradient iteration tolerance
cjtol
consolidation pressure
cons
Mohr-Coulomb plastic potential derivative
dq1
Mohr-Coulomb plastic potential derivative
dq2
Mohr-Coulomb plastic potential derivative
dq3
shear stress invariant
dsbar
dt
viscoplastic “time” step
f
current stress state
lode theta
Lode angle
phi
angle of internal friction
plastol
plastic iteration tolerance
presc
prescribed value of load/displacement
angle of dilation
psi
mean stress invariant
sigm
sine of angle of internal friction
sinph
New scalar logicals:
cj converged
set to
.TRUE.
if conjugate gradient iterations converged
plastic converged
set to
.TRUE.
if plastic iterations converged
