Civil Engineering Reference
In-Depth Information
!-----------------------extract the eigenvectors--------------------------
READ(10,*)nmodes
DO i=1,nmodes
kv=kh;kv(kdiag)=kv(kdiag)-diag(i); kv(1)=kv(1)+penalty
udiag=zero; udiag(1)=kv(1)
CALL sparin_gauss(kv,kdiag); CALL spabac_gauss(kv,udiag,kdiag)
udiag=rrmass*udiag; WRITE(11,'(A,I3,A)')" Eigenvector number",i," is:"
WRITE(11,'(6E12.4)')udiag(1:)/MAXVAL(ABS(udiag(1:)))
END DO
STOP
END PROGRAM p101
Scalar integers:
i
simple counter
iel
simple counter
idiag
skyline bandwidth
iwp SELECTED REAL KIND(15)
ifail
warning flag from
bisect
subroutine
j
simple counters
k
simple counters
nband
bandwidth of upper triangle
ndof
number of degrees of freedom per element
nels
number of elements
neq
number of degrees of freedom in the mesh
nmodes
number of eigenvectors required
nn
number of nodes in the mesh
nod
number of nodes per elements
nodof
number of degrees of freedom per node
nprops
number of material properties
np types
number of different property types
nr
number of restrained nodes
Scalar reals:
d12
set to 12
.
0
set to 1
.
0
one
10
20
penalty
set to 1
×
.
pt5
set to 0
5
10
−
30
eigenvalue tolerance set to 1
×
etol
set to 0
.
0
zero
Dynamic integer arrays:
etype
element property type vector
element steering vector
g
g
g
global element steering matrix
kdiag
diagonal term location vector
nf
nodal freedom matrix
num
element node number vector