Civil Engineering Reference
In-Depth Information
n
f
for the stream contour map, can then be computed from,
Qn
d
kH
n
f
=
(7.2)
where
Q
is the flow rate computed from the potential problem,
k
is the (isotropic) perme-
ability of the soil,
n
d
is the number of equipotential drops chosen for the total head contour
3 m
q
Fixed potential,
f
= 100
z
1
2
3
4
5
6
7
8
Q
10
=25
(sink)
Fixed
potential,
f
= 100
Fixed
potential,
f
=0
9
10
11
12
5 m
13
14
15
16
17
18
19
20
type_2d dir
'axisymmetric' 'r'
21
22
23
24
r
nxe nye np_types
3 5 1
fixed potential,
f
= 0
prop(k
x
,k
y
)
1.0 1.0
etype(not needed)
x_coords, y_coords
0.0 1.0 2.0 3.0
0.0 -1.0 -2.0 -3.0 -4.0 -5.0
loaded_nodes,(k,loads(k),i=1,loaded_nodes)
1
10 -25.0
fixed_freedoms,(node(i),value(i),i=1,fixed_freedoms)
16
1 100.0 2 100.0 3 100.0 4 0.0 5 100.0 8 0.0 9 100.0
12 0.0 13 100.0 16 0.0 17 100.0 20 0.0 21 0.0 22 0.0
23 0.0 24 0.0
20
Figure 7.9 Mesh and data for second Program 7.2 example
