16.4 Computer Program for Free Vibration Analysis of Framed Structures
573
Table 16.3 (continued)
enddo
enddo
emas(1,1) = mass
emas(2,2) = mass
emas(3,3) = mass
emas(7,7) = mass
emas(8,8) = mass
emas(9,9) = mass
return
end
*
Subroutine for assembling stiffness coefficients
subroutine assem(ind,ielm,estif)
parameter (ielt=50,inod=50,ntdf=inod*6,inb=50)
parameter (ibh=25,itdf=ntdf*ibh)
common /a/ ntv,nbh,gstif(itdf),gmass(itdf),gload(ntdf),kind(ntdf),
>mind(ntdf)
common /b/ nod(ielt,2),xg(inod),yg(inod),zg(inod)
dimension estif(12,12)
ndf =6
nnpe=2
do 60 i=1,nnpe
ig=(nod(ielm,i)-1)*ndf
do 60 j=1,ndf
mmi=ig+j
mmm=mmi-1
imi=(i-1)*ndf+j
do 60 k=1,nnpe
jg=(nod(ielm,k)-1)*ndf
do 60 l = 1,ndf
nnj=jg+l
jnj=(k-1)*ndf+l
if(nnj.gt.mmi) go to 55
if (mmi.le.nbh) then
kjj=mmm*(mmm+1)/2+nnj
kkj=kjj
else
kkj=nbh*(nbh+1)/2
lrow=(mmm-nbh+1)*nbh+nnj-mmi
kkj=kkj+lrow
endif
if(ind.eq.1) gstif(kkj)=gstif(kkj)+estif(imi,jnj)
if(ind.eq.2) gmass(kkj)=gmass(kkj)+estif(imi,jnj)
55
continue
60
continue
return
end
*
Subroutine crunch
subroutine crunch(ind)
parameter (ielt=50,inod=50,ntdf=inod*6,inb=50)
parameter (ibh=25,itdf=ntdf*ibh)
(continued)