Civil Engineering Reference
In-Depth Information
type_2d dir
'plane' 'x'
nxe nye cg_tol cg_limit np_types
5 5 0.0001 100 1
prop(cx,cy)
1.0 1.0
etype(not needed)
x_coords, y_coords
0.0 0.2 0.4 0.6 0.8 1.0
0.0 -0.2 -0.4 -0.6 -0.8 -1.0
dtim nstep theta npri nres ntime
0.01 150 0.5 10 31 100
loads(i),i=1,neq
0.0 0.0 0.0 0.0 0.0 0.0
100.0 100.0 100.0 100.0 100.0 0.0
100.0 100.0 100.0 100.0 100.0 0.0
100.0 100.0 100.0 100.0 100.0 0.0
100.0 100.0 100.0 100.0 100.0 0.0
100.0 100.0 100.0 100.0 100.0 0.0
fixed_freedoms,(node(i),value(i),i=1,fixed_freedoms)
11
1 0.0 2 0.0 3 0.0 4 0.0 5 0.0 6 0.0
12 0.0 18 0.0 24 0.0 30 0.0 36 0.0
nci
10
Figure 8.12
Data for Program 8.3 example
There are 36 equations
Time Pressure (node 31) cg iters
0.0000E+00 0.1000E+03
0.1000E+00 0.9009E+02 3
0.2000E+00 0.5845E+02 3
0.3000E+00 0.3580E+02 2
0.4000E+00 0.2178E+02 2
0.5000E+00 0.1324E+02 2
0.6000E+00 0.8051E+01 2
0.7000E+00 0.4895E+01 2
0.8000E+00 0.2976E+01 2
0.9000E+00 0.1809E+01 2
0.1000E+01 0.1100E+01 2
0.1100E+01 0.6687E+00 2
0.1200E+01 0.4065E+00 2
0.1300E+01 0.2472E+00 2
0.1400E+01 0.1503E+00 2
0.1500E+01 0.9135E-01 2
Figure 8.13
Results from Program 8.3 example
