Civil Engineering Reference
In-Depth Information
nxe nye cg_tol cg_limit fftol ltol np_types
8 4 0.0001 100 -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)
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.27
Data for Program 6.6 example
There are 184 equations
step load disp iters cg iters/plastic iter
1 0.2000E+03 -0.6593E-02 1 46.00
2 0.3000E+03 -0.1154E-01 4 51.50
3 0.3500E+03 -0.1614E-01 4 56.00
4 0.4000E+03 -0.2284E-01 4 63.00
5 0.4500E+03 -0.3277E-01 4 69.50
6 0.4800E+03 -0.4183E-01 4 78.75
7 0.5000E+03 -0.5031E-01 4 78.25
8 0.5100E+03 -0.5614E-01 4 88.50
9 0.5150E+03 -0.6062E-01 4 84.75
10 0.5200E+03 -0.2117E+00 17 99.59
Figure 6.28
Results from Program 6.6 example
as compared with the data for Program 6.5, is limited to the conjugate gradient tolerance
and iteration limit, set respectively to
cg
tol=0.0001
and
cg limit=100
.
The output is listed as Figure 6.28 which can be compared with Figure 6.26. Between 46
and 100 conjugate gradient iterations per plastic iteration were required, but the program
ran faster than Program 6.1 for the solution of this problem, even in scalar mode. In
parallel implementations, there is a trade-off between constant stiffness and tangent stiffness
methods because for the latter, all yielded elements are different, although their geometries
and elastic properties may be identical.
