feap * * infinitesimal deformation model for rotation of an elastic beam
0,0,1,2,2,4

NOPRint ! suppress full printing
! create the mesh
BLOCk
4,2,10,1,1,1
1,0,1
2,0,-1
3,10,-1
4,10,1

! define a phantom node for plotting the full range of the deformation
COORdinates
34,,-10,

! define a phantom element for plotting the full range of the deformation
ELEMents
21,,1,34,34,34,34

! set the boundary conditions
BOUNdary conditions
2,,1,1
34,,1,1

! set the load
FORCes
32,0,0,-1

! specify the type of element and the material properties
MATErial
solid
plane strain
small deformation 
displacement method
elastic isotropic 100000,0
density mass 0.1
quad 2 2 2




END
! specify batch execution
BATCh mode
! set up the plot contents
PLOT,defo
PLOT,mesh
PLOT,load
PLOT,axis
PRINt,off
NOPRint
! specify Newmark method
TRANsient,newm,0.25,0.5
! set up time step
DT,,.1
! advance ten time steps using the given load
LOOP,,10
TIME
! perform Newton-Raphson iterations 
LOOP,,20
TANGent,newton,1
NEXT
PLOT,wipe
PLOT,mesh
PLOT,load
NEXT
PROPoprional,1
! continue the analysis without the load
LOOP,,1990
TIME
! perform Newton-Raphson iterations 
LOOP,,20
TANGent,newton,1
NEXT
PLOT,wipe
PLOT,mesh
NEXT
END
1,0,0,0,1000000000,0,0




INTEractive mode
STOP