Civil Engineering Reference
In-Depth Information
There are 125 equations and the skyline storage is 3225
Node Total Head Flow rate
1 0.1543E-19 -0.1543E+01
2 0.2732E-19 -0.2732E+01
3 0.2028E-19 -0.2028E+01
4 0.1032E-19 -0.1032E+01
5 0.3514E-20 -0.3514E+00
6 0.3319E+01 -0.1332E-14
7 0.3183E+01 -0.7772E-15
8 0.2152E+01 0.6661E-15
9 0.1106E+01 -0.4441E-15
10 0.1026E-19 -0.1026E+01
11 0.9727E+01 0.3553E-14
12 0.6905E+01 0.8882E-15
13 0.5206E+01 0.4441E-15
14 0.2141E+01 0.8882E-15
15 0.2004E-19 -0.2004E+01
16 0.1121E+02 0.3553E-14
17 0.2098E+02 0.5995E-14
18 0.6877E+01 0.4885E-14
19 0.3154E+01 -0.1998E-14
20 0.2684E-19 -0.2684E+01
21 0.1689E+03 0.1000E+03
.
.
.
120 0.2361E-20 -0.2361E+00
121 0.5872E-20 -0.5872E+00
122 0.8971E-20 -0.8971E+00
123 0.7082E-20 -0.7082E+00
124 0.3505E-20 -0.3505E+00
125 0.1255E-20 -0.1255E+00
Inflow Outflow
0.1000E+03 -0.1000E+03
Figure 7.21 Results from Program 7.4 example
A truncated version of the output from the program is shown in Figure 7.21. The total
head is greatest at the central node, and equals about 169. Outflow occurs at all the outside
nodes of the mesh where the total head was fixed to zero. For example, the outflow at node
number 5 equals
−
0.3514.
Program 7.5 General two- (plane) or three-dimensional analysis of steady seepage.
No global conductivity matrix assembly. Diagonally preconditioned conjugate gradient
solver.
PROGRAM p75
!-------------------------------------------------------------------------
! Program 7.5 General two- (plane) or three-dimensional analysis of steady
! seepage. No global conductivity 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,fixed_freedoms,i,iel,k,loaded_nodes,nci,ndim, &
nels,neq,nip,nod,nn,np_types
REAL(iwp)::alpha,beta,cg_tol,det,one=1.0_iwp,penalty=1.0e20_iwp,up,
&
zero=0.0_iwp
