Civil Engineering Reference
In-Depth Information
L
1/3 kN
1/6 kN
z
5
6
1/6 kN
1/3 kN
unit cube
y
2
E = 100 kN/m
2
u
= 0.3
x
1
7
8
3
element nod
'tetrahedron' 4
nels nn nip nodof nst ndim np_types
6 8 1 3 6 3 1
prop(e,v,
g
)
100.0 0.3 0.0
etype(not needed)
g_coord
0.0 0.0 0.0 1.0 0.0 0.0 0.0 0.0 -1.0 1.0 0.0 -1.0
0.0 1.0 0.0 1.0 1.0 0.0 0.0 1.0 -1.0 1.0 1.0 -1.0
g_num
1 3 4 7 1 4 2 7 1 2 5 7
6 4 8 7 6 2 4 7 6 5 2 7
nr,(k,nf(:,k),i=1,nr)
1 0 0 1 2 1 0 1 3 0 0 0 4 1 0 0
5 0 1 1 7 0 1 0 8 1 1 0
loaded_nodes,(k,loads(nf(:,k)),i=1,loaded_nodes)
1 0.0 0.0 -0.1667 2 0.0 0.0 -0.3333
5 0.0 0.0 -0.3333 6 0.0 0.0 -0.1667
fixed_freedoms
0
4
Figure 5.30 Mesh and data for second Program 5.4 example
There are 12 equations and the skyline storage is 69
Node x-disp y-disp z-disp
1 0.0000E+00 0.0000E+00 -0.1000E-01
2 0.3000E-02 0.0000E+00 -0.1000E-01
3 0.0000E+00 0.0000E+00 0.0000E+00
4 0.3000E-02 0.0000E+00 0.0000E+00
5 0.0000E+00 0.3000E-02 -0.9999E-02
6 0.3000E-02 0.3000E-02 -0.1000E-01
7 0.0000E+00 0.3000E-02 0.0000E+00
8 0.3000E-02 0.3000E-02 0.0000E+00
The integration point (nip= 1) stresses are:
Element x-coord y-coord z-coord
sig_x sig_y sig_z tau_xy tau_yz tau_zx
1 0.2500E+00 0.2500E+00 -0.7500E+00
-0.1965E-04 -0.2100E-04 -0.1000E+01 0.0000E+00 0.0000E+00 0.0000E+00
2 0.5000E+00 0.2500E+00 -0.5000E+00
0.2302E-05 0.1389E-04 -0.1000E+01 -0.6475E-05 0.2974E-04 0.2326E-04
3 0.2500E+00 0.5000E+00 -0.2500E+00
0.1230E-04 0.2075E-05 -0.9999E+00 0.0000E+00 0.3641E-04 0.2974E-04
4 0.7500E+00 0.7500E+00 -0.7500E+00
-0.1087E-04 -0.9023E-05 -0.1000E+01 0.1556E-05 -0.7467E-05 -0.9316E-05
5 0.7500E+00 0.5000E+00 -0.5000E+00
0.6102E-05 -0.1297E-05 -0.1000E+01 -0.8752E-05 -0.2541E-04 -0.3189E-04
6 0.5000E+00 0.7500E+00 -0.2500E+00
0.9809E-05 0.1536E-04 -0.9999E+00 0.2158E-05 -0.3632E-04 -0.4299E-04
Figure 5.31 Results from second Program 5.4 example
