Civil Engineering Reference
In-Depth Information
END IF
IF(f>=zero.OR.(converged.OR.iters==limit))THEN
eload=MATMUL(devp,bee); bload=bload+eload*det*weights(i)
END IF
!-----------------------update the Gauss Point stresses-------------------
IF(converged.OR.iters==limit)tensor(:,i,iel)=stress
END DO gauss_pts_2
!-----------------------compute the total bodyloads vector----------------
bdylds(g)=bdylds(g)+bload; bdylds(0)=zero
END DO elements_3; IF(converged.OR.iters==limit)EXIT
END DO its; totd=totd+loads
WRITE(11,'(I5,2E12.4,I5)')iy,ptot,totd(nf(2,node(1))),iters
IF(iters==limit)EXIT
END DO load_incs
CALL dismsh(loads,nf,0.05_iwp,g_coord,g_num,13)
CALL vecmsh(loads,nf,0.05_iwp,0.1_iwp,g_coord,g_num,14)
STOP
END PROGRAM p61
Scalar integers:
i
simple counter
iel
simple counter
incs
number of load increments
iters
counts plastic iterations
iwp SELECTED REAL KIND(15)
iy
counts load increments
k
node number
limit
plastic iteration ceiling
loaded nodes
number of loaded nodes
ndim
number of dimensions
number of degrees of freedom per element
ndof
number of elements
nels
number of degrees of freedom in the mesh
neq
number of integrating points per element
nip
number of nodes in the mesh
nn
number of nodes per element
nod
number of degrees of freedom per node
nodof
number of material properties
nprops
number of different property types
np types
number of restrained nodes
nr
nst
number of stress (strain) terms
nxe
number of elements in
x
-direction
nye
number of elements in
y
-direction
Scalar reals:
ddt
used to find the critical time step
determinant of the Jacobian matrix
det
plastic potential derivative,
∂Q/∂σ
m
dq1
plastic potential derivative,
∂Q/∂J
2
dq2
