*HEADING
: NAFEMS FORCED VIBRATION; TEST 21 [S4R5]
*NODE
1, 0.
9, 10.
81, 0., 10.
89, 10.,10.
*NGEN,NSET=BOT
1,9
*NGEN,NSET=TOP
81,89
*NFILL,NSET=NALL
BOT,TOP,8,10
*NSET,NSET=LEFT,GEN
1,81,10
*NSET,NSET=RIGHT,GEN
9,89,10
*NSET,NSET=EDGE
TOP,BOT,RIGHT,LEFT
*NSET,NSET=NMID
45, 
*ELEMENT,TYPE=S4R5
1, 1,2,12,11
*ELGEN,ELSET=ALLE
1, 8,1,1, 8,10,10
*ELSET,ELSET=EMID
34,35,44,45
*SHELL SECTION,MATERIAL=MAT,ELSET=ALLE
1.0, 
*HOURGLASS STIFFNESS
,1.28192E9,2.88462E8,10
*MATERIAL,NAME=MAT
*DENSITY
8000.,
*ELASTIC
200.E9,0.3
*BOUNDARY
NALL,1,2
NALL,6
EDGE,1,3
BOT,5
TOP,5
LEFT,4
RIGHT,4
**
** Modes
** -----
**
*STEP,PERT
*FREQUENCY, eigensolver=lanczos
17, , 
*MODAL FILE
*NODE PRINT,F=0
*EL PRINT,F=0
*OUTPUT,FIELD,MODE LIST
1,4,9
*NODE OUTPUT
U
*ENDSTEP
**
** Harmonic Load (modal)
** ---------------------
**
*STEP,PERT
*STEADY STATE DYNAMICS,FREQUENCY SCALE=LINEAR
1.,78.17,50
*DLOAD,LOAD CASE=1
ALLE,P,1.E6
*MODAL DAMPING,MODAL=DIRECT
1,17,0.02
*SELECT EIGENMODES, GENERATE
1,17,1
*NODE PRINT,NSET=NMID,FREQ=50
U3, 
*EL PRINT,FRE=50,POSITION=AVERAGED AT NODES,ELSET=EMID
S11, 
*NODE FILE,NSET=NMID,FREQ=50
U, 
*EL FILE,FRE=50,ELSET=EMID
S,
*OUTPUT,FIELD,FREQ=50
*NODE OUTPUT,NSET=NMID
U 
*ENDSTEP
**
** Transient Load
** --------------
**
*STEP,PERT,AMP=STEP
*MODAL DYNAMIC
0.0001,0.05
*DLOAD
ALLE,P,1.E6
*MODAL DAMPING,MODAL=DIRECT
1,17,0.02
*SELECT EIGENMODES, GENERATE
1,17
*NODE PRINT,NSET=NMID,FREQ=109
U3, 
*EL PRINT,POSITION=AVERAGED AT NODES,ELSET=EMID,FREQ=109
S11, 
*NODE FILE,NSET=NMID,FREQ=109
U, 
*EL FILE,FRE=109,ELSET=EMID
S,
*OUTPUT,FIELD,FREQ=109
*NODE OUTPUT,NSET=NMID
U
*ENDSTEP
**
** Increase the damping and timestep to find steady state
**
*STEP,PERT,AMP=STEP
*MODAL DYNAMIC,CONTINUE=YES
0.002,1.0
*DLOAD
ALLE,P,1.E6
*MODAL DAMPING,MODAL=DIRECT
1,17,0.10
*SELECT EIGENMODES, GENERATE
1,17,1
*NODE PRINT,NSET=NMID,F=100
U3, 
*EL PRINT,POSITION=AVERAGED AT NODES,ELSET=EMID,F=100
S11, 
*NODE FILE,NSET=NMID,F=999
U, 
*EL FILE,ELSET=EMID,F=999
S,
*OUTPUT,FIELD,FREQ=100
*NODE OUTPUT,NSET=NMID
U 
*ENDSTEP