Civil Engineering Reference
In-Depth Information
nxe nye fftol ltol np_types
8 4 -1.0e-6 5.0e-5 1
prop(c
u
,e,v)
100.0 1.0e5 0.3
etype(not needed)
x_coords, y_coords
0.0 1.0 2.0 3.0 4.0 5.5 7.0 9.0 12.0
0.0 -1.25 -2.5 -3.75 -5.0
nr,(k,nf(:,k),i=1,nr)
33
1 0 1 2 0 1 3 0 1 4 0 1 5 0 1 6 0 1
7 0 1 8 0 1 9 0 0 14 0 0 23 0 0 28 0 0
37 0 0 42 0 0 51 0 0 56 0 0 65 0 0 70 0 0
79 0 0 84 0 0 93 0 0 98 0 0 107 0 0 112 0 0
113 0 0 114 0 0 115 0 0 116 0 0 117 0 0 118 0 0
119 0 0 120 0 0 121 0 0
loaded_nodes,(node(i),val(i,:),i=1,loaded_nodes)
5
1 0.0 -0.166667 10 0.0 -0.666667 15 0.0 -0.333333
24 0.0 -0.666667 29 0.0 -0.166667
tol limit
0.001 50
incs,(qinc(i),i=1,incs)
10
200.0 100.0 50.0 50.0 50.0 30.0 20.0 10.0 5.0 5.0
Figure 6.25
Data for Program 6.5 example
There are 184 equations and the skyline storage is 4130
step load disp iters
1 0.2000E+03 -0.6592E-02 2
2 0.3000E+03 -0.1154E-01 4
3 0.3500E+03 -0.1614E-01 4
4 0.4000E+03 -0.2285E-01 4
5 0.4500E+03 -0.3278E-01 4
6 0.4800E+03 -0.4184E-01 4
7 0.5000E+03 -0.5031E-01 4
8 0.5100E+03 -0.5613E-01 4
9 0.5150E+03 -0.6058E-01 4
10 0.5200E+03 -0.3307E+00 28
Figure 6.26
Results from Program 6.5 example
Program 6.6 Plane strain bearing capacity analysis of an elastic-plastic (von Mises)
material using 8-node rectangular quadrilaterals. Initial stress method. Tangent stiff-
ness. Consistent return algorithm. No global stiffness matrix assembly. Diagonally
preconditioned conjugate gradient solver.
PROGRAM p66
!-------------------------------------------------------------------------
! Program 6.6 Plane strain bearing capacity analysis of an elastic-plastic
! (von Mises) material using 8-node rectangular
! quadrilaterals. Initial stress method. Tangent stiffness.
! 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,incs,iters,iy,k,limit, &
loaded_nodes,ndim=2,ndof=16,nels,neq,nip=4,nn,nod=8,nodof=2,nprops=3, &
