Civil Engineering Reference
In-Depth Information
type_2d
'axisymmetric'
element nod dir
'quadrilateral' 4 'z'
nxe nye nip np_types
3 2 9 2
prop(e,v)
100.0 0.3
1000.0 0.45
etype
1 2 1 2 1 2
x_coords, y_coords
0.0 4.0 10.0 30.0
0.0 -4.0 -10.0
nr,(k,nf(:,k),i=1,nr)
8
1 0 1 2 0 1 3 0 0 6 0 0 9 0 0 10 0 1 11 0 1 12 0 0
loaded_nodes,(k,loads(nf(:,k)),i=1,loaded_nodes)
3
1 0.0 -2.6667 4 0.0 -23.3333 7 0.0 -24.0
fixed_freedoms
0
Figure 5.17
(
Continued from page 182
)
There are 12 equations and the skyline storage is 58
Node r-disp z-disp
1 0.0000E+00 -0.3176E-01
2 0.0000E+00 -0.3231E-02
3 0.0000E+00 0.0000E+00
4 0.1395E-02 -0.3990E-01
5 0.1165E-02 -0.2498E-02
6 0.0000E+00 0.0000E+00
7 0.1704E-02 -0.6046E-02
8 0.1330E-02 -0.4421E-03
9 0.0000E+00 0.0000E+00
10 0.0000E+00 0.2588E-02
11 0.0000E+00 0.3091E-03
12 0.0000E+00 0.0000E+00
The integration point (nip= 1) stresses are:
Element r-coord z-coord sig_r sig_z tau_rz sig_t
1 0.2000E+01 -0.2000E+01 -0.4140E+00 -0.1073E+01 -0.3452E-01 -0.4140E+00
2 0.2000E+01 -0.7000E+01 -0.4776E+00 -0.9072E+00 0.6508E-01 -0.4776E+00
3 0.7000E+01 -0.2000E+01 -0.2933E+00 -0.7099E+00 0.1180E+00 -0.2810E+00
4 0.7000E+01 -0.7000E+01 -0.4316E+00 -0.6101E+00 0.1308E+00 -0.3796E+00
5 0.2000E+02 -0.2000E+01 -0.3200E-01 -0.5814E-01 0.1082E-01 -0.2325E-01
6 0.2000E+02 -0.7000E+01 -0.1090E+00 -0.9367E-01 0.4470E-01 -0.7455E-01
Figure 5.18 Results from fifth Program 5.1 example
