Civil Engineering Reference
In-Depth Information
nxe nye nod
5 1 4
np_types
1
prop(e,v,rho)
1.0 0.0 1.0
etype (not needed)
x_coords, y_coords
0.0 0.8 1.6 2.4 3.2 4.0
0.0 -1.0
nr,(k,nf(:,k),i=1,nr)
2
1 0 0 2 0 0
nmodes
5
el er lalfa leig lx lz acc
0.0 5.0 500 20 80 500 1.0e-6
Figure 10.10 Data for Program 10.3 and 10.4 examples
There are 20 equations and the half-bandwidth is 7
It took 22 iterations
Eigenvalues in the range 0.0000E+00 and 0.5000E+01 are:
0.4931E-02 0.1499E+00 0.1555E+00 0.9573E+00 0.1493E+01 0.3043E+01
0.4688E+01
Eigenvector number 1 is:
-0.7275E-01 0.7082E-01 0.7275E-01 0.7082E-01 -0.1229E+00 0.2394E+00
0.1229E+00 0.2394E+00 -0.1522E+00 0.4700E+00 0.1522E+00 0.4700E+00
-0.1648E+00 0.7311E+00 0.1648E+00 0.7311E+00 -0.1677E+00 0.1000E+01
0.1677E+00 0.1000E+01
Eigenvector number 2 is:
-0.2156E+00 0.4086E+00 0.2156E+00 0.4086E+00 -0.6107E-01 0.7885E+00
0.6107E-01 0.7885E+00 0.2617E+00 0.6519E+00 -0.2617E+00 0.6519E+00
0.5123E+00 -0.5054E-01 -0.5123E+00 -0.5054E-01 0.5925E+00 -0.1000E+01
-0.5925E+00 -0.1000E+01
Eigenvector number 3 is:
0.3090E+00 -0.2022E-06 0.3090E+00 -0.2015E-06 0.5878E+00 -0.3972E-06
0.5878E+00 -0.3961E-06 0.8090E+00 -0.3516E-06 0.8090E+00 -0.3508E-06
0.9511E+00 -0.3366E-07 0.9511E+00 -0.3273E-07 0.1000E+01 0.4045E-06
0.1000E+01 0.4053E-06
Eigenvector number 4 is:
-0.1193E+00 0.9255E+00 0.1193E+00 0.9255E+00 0.4733E+00 0.5949E+00
-0.4733E+00 0.5949E+00 0.3215E+00 -0.6339E+00 -0.3215E+00 -0.6339E+00
-0.5413E+00 -0.5437E+00 0.5413E+00 -0.5437E+00 -0.1000E+01 0.9777E+00
0.1000E+01 0.9777E+00
Eigenvector number 5 is:
-0.8090E+00 -0.4830E-07 -0.8090E+00 -0.5335E-07 -0.9511E+00 -0.4265E-07
-0.9511E+00 -0.4216E-07 -0.3090E+00 0.3175E-07 -0.3090E+00 0.3655E-07
0.5878E+00 0.3762E-07 0.5878E+00 0.4084E-07 0.1000E+01 -0.7482E-07
0.1000E+01 -0.7508E-07
Figure 10.11 Results from Program 10.3 example
The output from the program is listed as Figure 10.11. It consists of the number of
Lanczos iterations (22 in this case) and the computed eigenvalues in the range 0
2
.
0
<ω
<
5
0. The first five eigenvectors are also printed.
Since the consistent mass assumption was made, the eigenvalues differ somewhat from
those computed by Program 10.2. For example the first eigenvalue is computed as
.
ω
=
