Civil Engineering Reference
In-Depth Information
thus 1
5 to each mid-side
node in each direction. This weighting is assigned to the 8-node element by subroutine
elmat
.
The mesh and data shown in Figure 10.6 are virtually identical to those used for the
analysis with 4-node elements. The only differences lies in the data for
nod
, which is now
read as 8, and the boundary condition data at the built-in end of the cantilever involves 3
fixed nodes.
/
20 of the total mass is assigned to each corner node and 1
/
There are 20 equations, the half-bandwidth is 7
and the skyline storage is 114
The eigenvalues are:
0.4595E-02 0.1053E+00 0.1529E+00 0.5169E+00 0.1226E+01 0.1288E+01
0.2058E+01 0.2352E+01 0.2415E+01 0.2710E+01 0.2719E+01 0.3114E+01
0.3125E+01 0.3188E+01 0.3381E+01 0.3665E+01 0.3960E+01 0.4211E+01
0.4962E+01 0.6097E+01
Eigenvector number 1 is:
0.7195E-01 -0.6932E-01 -0.7195E-01 -0.6932E-01 0.1224E+00 -0.2362E+00
-0.1224E+00 -0.2362E+00 0.1526E+00 -0.4662E+00 -0.1526E+00 -0.4662E+00
0.1662E+00 -0.7285E+00 -0.1662E+00 -0.7285E+00 0.1695E+00 -0.1000E+01
-0.1695E+00 -0.1000E+01
Eigenvector number 2 is:
0.2487E+00 -0.4240E+00 -0.2487E+00 -0.4240E+00 0.1022E+00 -0.8727E+00
-0.1022E+00 -0.8727E+00 -0.2584E+00 -0.7980E+00 0.2584E+00 -0.7980E+00
-0.5732E+00 -0.7556E-01 0.5732E+00 -0.7556E-01 -0.6870E+00 0.1000E+01
0.6870E+00 0.1000E+01
Eigenvector number 3 is:
0.3090E+00 0.1555E-14 0.3090E+00 0.1572E-14 0.5878E+00 0.3059E-14
0.5878E+00 0.2986E-14 0.8090E+00 0.2674E-14 0.8090E+00 0.2728E-14
0.9511E+00 -0.5347E-16 0.9511E+00 0.0000E+00 0.1000E+01 -0.3174E-14
0.1000E+01 -0.3478E-14
Eigenvector number 4 is:
0.2413E+00 -0.9151E+00 -0.2413E+00 -0.9151E+00 -0.4612E+00 -0.8557E+00
0.4612E+00 -0.8557E+00 -0.5141E+00 0.3742E+00 0.5141E+00 0.3742E+00
0.4120E+00 0.6579E+00 -0.4120E+00 0.6579E+00 0.1000E+01 -0.7050E+00
-0.1000E+01 -0.7050E+00
Eigenvector number 5 is:
0.2575E+00 0.9693E+00 -0.2575E+00 0.9693E+00 0.8105E+00 -0.2304E+00
-0.8105E+00 -0.2304E+00 -0.3806E+00 -0.5583E+00 0.3806E+00 -0.5583E+00
-0.6787E-01 0.6923E+00 0.6787E-01 0.6923E+00 0.1000E+01 -0.2476E+00
-0.1000E+01 -0.2476E+00
Figure 10.5 Results from first Program 10.2 example
1
4
6
E= 1.0
u
1.0
=0.0
r
= 1.0
2
3
5
28
4.0
Figure 10.6 Mesh and data for second Program 10.2 example (
Continued on page 451
)
