16.6 Computer Program for Finite Element Free Vibration Analysis of Plates
595
Table 16.9 (continued)
k=i-4
15
nod(ii,i)= (nx21+nx1)*(nxi+1)+2*(ii-nx*nxi)+2-k
10
continue
write(8,*) '
node
numbers'
do 50 ii=1,nelm
write (8,*) '
The element number.', ii
write (8,101) (nod(ii,i),i=1,nn)
50
continue
c****************************************************************************
c****************************************************************************
c
automatic mesh generation over
c****************************************************************************
C
generation of rgd matrix starts
c****************************************************************************
do 23 i=1,5
do 23 j=1,5
23
rgd(i,j)=0.0
rgd(1,1) = dx
rgd (1,2)= d1
rgd(2,1) = d1
rgd(2,2) =dy
rgd(3,3)= dxy
rgd(4,4) =sx
rgd(5,5) = sy
c****************************************************************************
c
generation of the rgd matrix over
c****************************************************************************
c****************************************************************************
do 12 i=1,nenp
do 12 j=1,nenp
12
stiff(i,j)=0.0
c****************************************************************************
c
element coordinate information is prescribed
c****************************************************************************
divx=nx
divy=ny
exl=a/divx
eyl=b/divy
xg(1) = eyl*cos(phai)
xg(2) = exl/2.0+xg(1)
yg(1) = eyl*sin(phai)
yg(2) = yg(1)
xg(3) = exl+xg(1)
yg(3) = yg(1)
xg(4) = exl+xg(1)/2.0
yg(4) = yg(1)/2.0
xg(5) = exl
yg(5) = 0.0
xg(6) = exl/2.0
yg(6) = 0.0
xg(7) = 0.0
yg(7) = 0.0
xg(8) = xg(1)/2.0
yg(8) = yg(1)/2.0
zx(1)=-1.0
zy(1)=1.0
zx(2)=0.0
zy(2)=1.0
zx(3)=1.0
zy(3)=1.0
zx(4)=1.0
zy(4)=0.0
(continued)