*HEADING
 EARTHQUAKE ANALYSIS (S8R) VIA MULTIPLE BASE MOTIONS  
 MODAL DYNAMICS MODEL WITH PRIMARY BASE AND SECONDARY BASE (NODE21)
 A TIME SHIFT IS APPLIED TO THE AMP RECORD FOR THE SECONDARY BASE
*NODE
1,300.
21, 
101,300.,.5
121,0.,.5
201,300.,1.
221,0.,1.
*NSET,NSET=ENDS
1,21,201,221
*NSET,NSET=NODE11
11,
*NSET,NSET=NOUT
1,11,21,101,111,211,201,211,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
NALL,YSYMM
  1, 1
  1, 3
  1, 5
101, 1
101, 3
101, 5
201, 1
201, 3
201, 5
 21, 1
 21, 5
121, 1
121, 5
221, 1
221, 5
** THIS INPUT ASSUMES THAT bmk/multibasemotion_quake.inp 
** HAS BEEN COPIED TO QUAKE01.AMP 
** AMPLITUDE CURVE DATA FROM FIRST 2 SEC OF 3-1-2-I
*AMPLITUDE,VALUE=RELATIVE,NAME=EQ,INPUT=QUAKE01.AMP
** THIS INPUT ASSUMES THAT bmk/multibasemotion_quake_shift.inp 
** HAS BEEN COPIED TO QUAKE21.AMP 
** AMPLITUDE CURVE DATA FROM FIRST 2 SEC OF 3-1-2-I WITH TIME SHIFT = .25 SEC
*AMPLITUDE,VALUE=RELATIVE,NAME=EQ2,INPUT=QUAKE21.AMP
*RESTART,WRITE,FREQUENCY=25
*STEP
*FREQUENCY,EIGENSOLVER=LANCZOS, SIM=NO
10,,,-1000.
*BOUNDARY,BASE NAME=NODE21
 21,3
121,3
221,3
*EL PRINT,FREQUENCY=0
*NODE PRINT
U,
*MODAL FILE
*OUTPUT,FIELD,FREQ=0
*END STEP
*STEP
 EARTHQUAKE INPUT
*MODAL DYNAMIC
.01,2.0
*SELECT EIGENMODES,GENERATE
1,6
*BASE MOTION,DOF=3,AMPLITUDE=EQ,SCALE=386.09
*BASE MOTION,DOF=3,AMPLITUDE=EQ2,SCALE=386.09,BASE NAME=NODE21
*PRINT,FREQUENCY=50
*EL PRINT,ELSET=PR,FREQUENCY=50
S,E
*NODE PRINT,FREQUENCY=100
U,
TU,
V,
A,
*MODAL PRINT,FREQUENCY=50
GU,GV,GA,
BM, 
*NODE FILE,NSET=NODE11
U,TU,
*NODE FILE,NSET=NOUT,FREQUENCY=5
U,V,A,TU
*MODAL FILE,FREQUENCY=25
GU,GV,GA,BM
*OUTPUT,HIST,FREQ=25
*MODAL OUTPUT
GU,GV,GA,BM
*OUTPUT,HIST,FREQ=5
*NODE OUTPUT,NSET=NOUT
U,V,A,TU
*OUTPUT,HIST,FREQ=1
*NODE OUTPUT,NSET=NODE11
U,TU
*OUTPUT,FIELD,FREQ=5
*NODE OUTPUT
U, 
*END STEP