Civil Engineering Reference
In-Depth Information
0
0
-5
y
-10
10
0
20
10
30
x
20
40
z
30
50
40
60
Figure 6.48
Deformed mesh at failure from Program 6.10 example
deformed mesh corresponding to this unconverged solution. The slumping of the slope
towards its centreline is clearly seen.
Program 6.11 Three-dimensional strain analysis of an elastic-plastic (Mohr-Cou-
lomb) slope using 20-node hexahedra. Viscoplastic strain method. No global stiffness
matrix assembly. Diagonally preconditioned conjugate gradient solver.
PROGRAM p611
!-------------------------------------------------------------------------
! Program 6.11 Three-dimensional strain analysis of an elastic-plastic
! (Mohr-Coulomb) slope using 20-node hexahedra. Viscoplastic
! strain method. No global stiffness matrix assembly.
! Diagonally preconditioned conjugate gradient solver.
!-------------------------------------------------------------------------
USE main; USE geom; IMPLICIT NONE
INTEGER,PARAMETER::iwp=SELECTED_REAL_KIND(15)
INTEGER::cg_iters,cg_limit,cg_tot,i,iel,ifix,iters,iy,k,limit,ndim=3, &
ndof=60,nels,neq,nip=8,nn,nod=20,nodof=3,nprops=6,np_types,nsrf,nst=6, &
nx1,nx2,ny1,ny2,nze
REAL(iwp)::alpha,beta,cf,cg_tol,ddt,det,dq1,dq2,dq3,dsbar,dt=1.0e15_iwp, &
d1,d4=4.0_iwp,d180=180.0_iwp,e,f,fmax,h1,h2,lode_theta,one=1.0_iwp,phi,&
phif,pi,psi,psif,sigm,snph,start_dt=1.e15_iwp,s1,tnph,tnps,tol,
&
two=2.0_iwp,up,v,w1,w2,zero=0.0_iwp
CHARACTER(LEN=80)::element='hexahedron'; LOGICAL::converged,cg_converged
!-----------------------dynamic arrays------------------------------------
INTEGER,ALLOCATABLE::etype(:),g(:),g_g(:,:),g_num(:,:),nf(:,:),num(:)
REAL(iwp),ALLOCATABLE::bdylds(:),bee(:,:),bload(:),coord(:,:),d(:),
&
dee(:,:),der(:,:),deriv(:,:),devp(:),diag_precon(:),eld(:),eload(:),
&
