Civil Engineering Reference
In-Depth Information
All nodes set to
f
= 100
at time t = 0.
Nodes maintained at
f
= 0 for t>0
y
plane of symmetry
z
101
105
49
52
36
33
26
17
20
30
x
1
2
3
4
5
1
3
2
4
24
6
5
8
28
64
1.0
11
125
48
12
9
32
1.0
16
13
16
50
21
25
1.0
element
'hexahedron'
planes of symmetry
nels nn nip nod ndim np_types
64 125 8 8 3 1
prop(cx,cy,cz)
1.0 1.0 1.0
etype(not needed)
g_coord
0.00 0.00 0.00 0.25 0.00 0.00 0.50 0.00 0.00 0.75 0.00 0.00
g_coord data for nodes 5-120 omitted here
0.00 1.00 -1.00 0.25 1.00 -1.00 0.50 1.00 -1.00 0.75 1.00 -1.00
1.00 1.00 -1.00
g_num
6 1 2 7 31 26 27 32 7 2 3 8 32 27 28 33
g_num data for elements 3-62 omitted here
98 93 94 99 123 118 119 124 99 94 95 100 124 119 120 125
dtim nstep theta npri nres
0.01 150 0.5 10 21
loads(i),i=1,neq
Top, right and back faces set to zero. All other nodes set to 100.0
fixed_freedoms,(node(i),value(i),i=1,fixed_freedoms)
61
1 0.0 2 0.0 3 0.0 4 0.0 5 0.0 10 0.0 15 0.0 20 0.0
fixed freedom data for 48 nodes omitted here
121 0.0 122 0.0 123 0.0 124 0.0 125 0.0
Figure 8.19
Mesh and data for Program 8.6 example
) coordinates of the mesh read into
g coord
,
followed by the 64 groups of 8 node numbers attached to each element read into
g num
.
If dealing with a three-dimensional 8-node element for example, the order in which the
node numbers are read must follow the ordering described in Appendix B. Due to the
volume of data required in this example, a truncated version of the data is actually shown
in Figure 8.19.
The next data involves the 125 (
x,y,z
