16.4 Computer Program for Free Vibration Analysis of Framed Structures
567
Table 16.3 (continued)
dimension freq(20)
dimension title(15)
character ans*2
real ixx,iyy,izz,len,load,length
open(1,file='frame.in')
open(2,file='frame.out')
150
format(15a4)
151
format(/,'Youngs Modulus.....= ',g12.4,/'Density.....= ',g12.4)
152
format(/'No. of elements.....= ',i5,/'No. of nodes.....= ',i5)
153
format(/'Cross-Sectional Properties',/'Element
Area
> Ixx
Iyy
Izz',/(i5,4g12.4))
154
format(/'Nodal Connectivity',/'Element
Nodal
>Connectivity',/(i5,5x,2i8))
155
format(/'Nodal Coordinates',/'Node
X-Cor
Y-Cor
>
Z-Cor',/(i5,3g15.6))
156
format(/'No. of boundary points......= ',i5)
157
format(/'Bounary Codes',/'Node
Boundary Codes',/(i3,7x,6i2))
158
format(/'No. of loads acting on the structure......= ',i5)
159
format(/'Nodes on which load is acting and its magnitude',/
>(i3,6g10.4))
160
format(/'Total d.o.f in the structure.....= ',i5,/'Half Bandwidth
>.....= ',i5)
161
format(/'Element Stiffness Matrix',/(6g12.4))
162
format(/'Element Mass Matrix',/(6g12.4))
163
format(/'Overall Stiffness Matrix',/(6g12.4))
164
format(/'Overall Load Matrix',/(6g12.4))
165
format(/'Overall Mass Matrix',/(6g12.4))
166
format(/'Nodal Displacements',/(i4,6g12.4))
167
format(/'Natural Frequencies',/(i3,5x,g12.4))
168
format(a2)
print *,'Do you want to carry Static Analysis? If yes enter : y
> or Y.'
print *,'Else Free Vibration Analysis is carried out'
read 168,ans
read (1,150)title
write(2,150)title
*
Specify the material properties
read(1,*)emod,dens
write(2,151)emod,dens
*
Specify the no. of elements and nodes in the structure
read(1,*)nelm,nnod
write(2,152)nelm,nnod
*
Specify the properties of the elements
read(1,*)(area(iel),ixx(iel),iyy(iel),izz(iel),iel=1,nelm)
write(2,153)(iel,area(iel),ixx(iel),iyy(iel),izz(iel),iel=1,nelm)
*
Specify the nodal connectivity of the elements
read(1,*)((nod(iel,ind),ind=1,2),iel=1,nelm)
write(2,154)(iel,(nod(iel,ind),ind=1,2),iel=1,nelm)
*
Specify the nodal coordinates
read(1,*)(xg(ind),yg(ind),zg(ind),ind=1,nnod)
write(2,155)(ind,xg(ind),yg(ind),zg(ind),ind=1,nnod)
(continued)