*HEADING
 EARTHQUAKE ANALYSIS (B23) VIA MULTIPLE BASE MOTION  
 MODAL DYNAMICS MODEL WITH TWO SECONDARY BASES (NODE1 and NODE21)
 NO TIME SHIFT
*NODE
1 ,300.
21, 
*NSET,NSET=ENDS
1,21
*NSET,NSET=NOUT
1,11,21
*NGEN
 1,21
*ELSET,ELSET=PR
10, 
*ELEMENT,TYPE=B23
1,1,2
*ELGEN,ELSET=EALL
1,20
*MATERIAL,NAME=A1
*ELASTIC
30.E6, 
*DENSITY
 .000728, 
*BEAM SECTION,SECTION=RECT,MATERIAL=A1,ELSET=EALL
1. , 2.
** THIS INPUT ASSUMES THAT MULTIBASEMOTION_QUAKE.INP HAS BEEN COPIED 
** TO QUAKE01.AMP 
** AMPLITUDE CURVE DATA FROM FIRST 2 SEC OF cantilever_quakedata.inp
*AMPLITUDE,VALUE=RELATIVE,NAME=EQ,INPUT=QUAKE01.AMP
*RESTART,WRITE,FREQUENCY=100
*BOUNDARY
 1, 1
 1, 6
21, 1 
21, 6
*STEP
*FREQUENCY,EIGENSOLVER=LANCZOS
10,,,-10000.
*EL PRINT,FREQUENCY=0
*BOUNDARY,BASE NAME=NODE1
 1,2
*BOUNDARY,BASE NAME=NODE21
21,2
*NODE PRINT
U,
*MODAL FILE
*END STEP
*STEP
 EARTHQUAKE INPUT 
*MODAL DYNAMIC
 .01 , 2.0
*SELECT EIGENMODES,GENERATE
1,7
*BASE MOTION,DOF=2,AMPLITUDE=EQ,SCALE=386.09,BASE NAME=NODE1
*BASE MOTION,DOF=2,AMPLITUDE=EQ,SCALE=386.09,BASE NAME=NODE21
*PRINT,FREQUENCY=50
*EL PRINT,ELSET=PR,FREQUENCY=500
S,E
*NODE PRINT,FREQUENCY=100
U,TU,
V,TV,
A,TA,
*MODAL PRINT,FREQUENCY=100
GU,GV,GA,
BM, 
*NODE FILE,NSET=NOUT
U,V,A,TU,TV,TA
*MODAL FILE,FREQUENCY=10
BM,
*OUTPUT,FIELD
*NODE OUTPUT,NSET=NOUT
U,V,A,TU,TV,TA
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=NOUT
U,V,A,TU,TV,TA
*OUTPUT,HISTORY,FREQUENCY=10
*MODAL OUTPUT
BM,
*END STEP