Civil Engineering Reference
In-Depth Information
Program 11.1
Warburton(1964)
t
0
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6
1.8
Figure 11.3
Computed tip displacement-time history from Program 11.1 example
the analytical solution (e.g. Warburton, 1964) for this problem during the loading phase
when 0
≤
t
≤
T
is given by:
441 sin
πt
T
216 sin
2
πt
T
v(t)
=
0
.
−
0
.
(11.2)
and for the free vibration phase when
t>T
by,
v(t)
=−
0
.
432 sin
{
6
.
284
(t
−
T)
}
(11.3)
It is clear from Figure 11.3 that this analytical solution is virtually indistinguishable from
the computed result.
Program 11.2 Forced vibration analysis of an elastic solid in plane strain using 4- or
8-node rectangular quadrilaterals. Lumped mass. Mesh numbered in
x
-or
y
-direction.
Modal superposition.
PROGRAM p112
!-------------------------------------------------------------------------
! Program 11.2 Forced vibration of an elastic solid in plane strain
! using 4- or 8-node rectangular quadrilaterals. Lumped mass.
! Mesh numbered in x- or y-direction. Modal superposition.
!-------------------------------------------------------------------------
USE main; USE geom; IMPLICIT NONE
INTEGER,PARAMETER::iwp=SELECTED_REAL_KIND(15)
INTEGER::i,idiag,iel,ifail,j,jj,k,nband,ndim=2,ndof,nels,neq,nmodes,nn, &
nod,nodof=2,npri,nprops=3,np_types,nr,nres,nstep,nxe,nye
