Download FEAP - - A Finite Element Analysis Program

Transcript
CHAPTER 5. ADDING ELEMENTS
character*4
common /bdata/
o,head
o,head(20)
integer
common /cdata/
numnp,numel,nummat,nen,neq,ipr
numnp,numel,nummat,nen,neq,ipr
31
integer
nstep,niter,naugm, titer,taugm,tform
common /counts/ nstep,niter,naugm, titer,taugm,tform
integer
iaugm, iform,intvc,iautl, nstepa
common /counts/ iaugm, iform,intvc,iautl, nstepa
real*8
dm
integer
n,ma,mct,iel,nel
common /eldata/ dm,n,ma,mct,iel,nel
real*8
tt
common /elplot/ tt(200)
real*8
bpr,
ctan,
psil
common /eltran/ bpr(3),ctan(3),psil
real*8
ut
common /eluser/ ut(200)
integer
common /hdata/
nh1,nh2,nh3,nlm
nh1,nh2,nh3,nlm
integer
ior,iow,ilg
common /iofile/ ior,iow,ilg
integer
nph
real*8
erav,j_int
,jshft
common /prstrs/ nph,ner,erav,j_int(3),jshft
integer
common /sdata/
ndf,ndm,nen1,nst,nneq,ndl,nnlm,nadd
ndf,ndm,nen1,nst,nneq,ndl,nnlm,nadd
real*8
common /tdata/
ttim,dt,c1,c2,c3,c4,c5, chi
ttim,dt,c1,c2,c3,c4,c5, chi
real*8
integer
common /
hr
mr
/ hr(1),mr(1000)
Figure 5.2: FEAP Element Common Blocks.
expressed in terms of nodal displacements, velocities, and accelerations given by ui (t),
u̇i (t), and üi (t), respectively. We denote the differential equation for node-i as the