*HEADING
 CANTILEVER BEAM SUBJECT TO SEISMIC LOAD USING ITSCYL ELEMENT
*NODE
 1, 
 21,20.
 22,20.
*NGEN
 1,21
*NSET,NSET=ALL,GENERATE
 1,21,1
*NSET,NSET=NPRINT
 10,21,22
*NSET,NSET=GROUND
 1,22
*ELEMENT,TYPE=B21
 1,1,2
*ELGEN,ELSET=BEAM
 1,20,1
*ELEMENT,TYPE=ITSCYL,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.
*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
*FREQUENCY, eigensolver=lanczos
 6, , ,
*EL PRINT
*MODAL FILE
*END STEP
*STEP
*DYNAMIC
 1.E-6,1.E-6
*DLOAD,AMPLITUDE=DYN
 BEAM,PY,-0.004
*EL PRINT
*NODE PRINT
*EL FILE
*NODE FILE,NSET=NPRINT
 U, 
*NODE FILE,NSET=GROUND
 RF, 
*END STEP
*STEP
*DYNAMIC
 .001,.001
*DLOAD,AMPLITUDE=DYN
 BEAM,PY,-0.004
*END STEP
*STEP,INC=1000
*DYNAMIC,SUBSPACE
 3.125E-5,.019
*SELECT EIGENMODES,GENERATE
1,6,1
*DLOAD,AMPLITUDE=DYN
 BEAM,PY,-0.004
*EL PRINT,ELSET=BEAM,FREQUENCY=1000
SF, 
SE, 
*EL PRINT,ELSET=ITS,FREQUENCY=1000
S11, S12, S13, SS1, SS2,
PE12, PE13, E11, E12, E13,
*NODE PRINT,NSET=ALL,FREQUENCY=1000
 U,
*EL FILE,ELSET=BEAM,FREQUENCY=1000
SF, 
SE, 
*EL FILE,ELSET=ITS,FREQUENCY=1000
S,
PE,
E, 
*NODE FILE,NSET=ALL,FREQUENCY=1000
 U,
*END STEP