*HEADING
 EARTHQUAKE ANALYSIS  S8R                 
*NODE
1, 0.,300.
21, 
101,0.,300.,-.5
121,0.,0.,-.5
201,0.,300.,-1.
221,0.,0.,-1.
*NSET,NSET=ENDS
1,21,201,221
*NGEN,NSET=NALL
 1,21
 101,121
 201,221
*ELSET,ELSET=PR
10, 
*ELEMENT,TYPE=S8R
1,1,3,203,201,2,103,202,101
*ELGEN,ELSET=EALL
1,10,2
*MATERIAL,NAME=A1
*ELASTIC
30.E6, 
*DENSITY
 .000728, 
*SHELL SECTION,MATERIAL=A1,ELSET=EALL
2., 
*BOUNDARY
21, 1, 6
121, 1, 6
221, 1, 6
NALL,3,5
** AMPLITUDE CURVE DATA FROM CANTILEVER_QUAKEDATA
** THIS INPUT ASSUMES THAT CANTILEVER_QUAKEDATA.INP
** HAS BEEN COPIED TO QUAKE.AMP
*AMPLITUDE,VALUE=RELATIVE,INPUT=QUAKE.AMP,NAME=EQ
*RESTART,WRITE,FREQUENCY=25
*STEP
*FREQUENCY,EIGENSOLVER=LANCZOS, SIM=NO
10,
*EL PRINT,FREQUENCY=0
*NODE PRINT
U,
*MODAL FILE
*output,field
*node output
u,
*END STEP
*STEP
 EARTHQUAKE INPUT
*MODAL DYNAMIC
.01,2.0
*SELECT EIGENMODES,GENERATE
1,6,1
*MODAL DAMPING
1,6
*BASE MOTION,DOF=1,AMPLITUDE=EQ,SCALE=386.09
*PRINT,FREQUENCY=50
*EL PRINT,ELSET=PR,FREQUENCY=50
S,E
*NODE PRINT,FREQUENCY=100
U,
TU,
V,
A,
*output,field,frequency=100
*node output
u,tu,v,a
*output,history
*node output,nset=ends
U,V,A,
*modal output
GU,GV,GA,BM
*NODE FILE,NSET=ENDS
U,V,A,
*MODAL FILE
GU,GV,GA,BM
*END STEP