*HEADING
 CANTILEVER BEAM SUBJECT TO SEISMIC LOAD        
*NODE
 1, 
 21,20.
 22,20.
*NGEN
 1,21
*NSET,NSET=ALL,GENERATE
 1,21,1
*NSET,NSET=NPRINT
 10,21
*NSET,NSET=GROUND
 1,22
*ELEMENT,TYPE=B21
 1,1,2
*ELGEN,ELSET=BEAM
 1,20,1
*ELEMENT,TYPE=ITSUNI,ELSET=ITS
 21,21,22
*BEAM GENERAL SECTION,ELSET=BEAM,DENSITY=.004
 1.,2.,0.,2.,4.
 0.,0.,-1.
 30.E6,11.54E6
*ITS,ELSET=ITS
 0.,0.,1.,0.,0.,0.,1.,0.
*SPRING,NONLINEAR

 -2.,-1.5E-5
 0.,-.5E-5
 0.,0.
 0.,.5E-5
 2.,1.5E-5
*BOUNDARY
 1,1
 1,2
 1,6
 22,2
 ALL,1
*AMPLITUDE,TIME=TOTAL TIME,VALUE=RELATIVE,NAME=DYN
 0.,0.,0.001,0.0,0.011,10.,0.016,-10.
*STEP
*DYNAMIC,DIRECT
 1.E-6,1.E-6
*DLOAD,AMPLITUDE=DYN
 BEAM,PY,-0.004
*EL PRINT,FREQUENCY=0
*NODE PRINT,FREQUENCY=0
*NODE FILE,NSET=NPRINT,FREQUENCY=1
 U, 
*OUTPUT,FIELD,FREQ=1
*NODE OUTPUT,NSET=NPRINT
U,
*NODE FILE,NSET=GROUND,FREQUENCY=1
 RF, 
*OUTPUT,FIELD,FREQ=1
*NODE OUTPUT,NSET=GROUND
RF,
*END STEP
*STEP
*DYNAMIC,DIRECT
 .001,.001
*DLOAD,AMPLITUDE=DYN
 BEAM,PY,-0.004
*END STEP
*STEP,INC=1000
*DYNAMIC,ALPHA=0.,NOHAF,DIRECT
 3.125E-5,.019
*DLOAD,AMPLITUDE=DYN
 BEAM,PY,-0.004
*PRINT,RESIDUAL=NO,FREQUENCY=10
*EL PRINT,ELSET=BEAM,FREQUENCY=1000
*NODE PRINT,NSET=ALL,FREQUENCY=100
 U, 
*EL FILE,ELSET=BEAM,FREQUENCY=1000
 SF, 
*OUTPUT,FIELD,FREQ=1000
*ELEMENT OUTPUT,ELSET=BEAM
SF,
*END STEP