*HEADING
:NAFEMS FORCED VIBRATION; TEST 13H;  [S8R]     
(DIRECT INTEGRATION) 
 GAUSS INTEGRATION AT 5 SECTION POINTS
*RESTART, WRITE
*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=S8R
1, 1,3,23,21,2,13,22,11
*ELGEN,ELSET=ALLE
1, 4,2,1, 4,20,10
*ELSET,ELSET=EMID
12,13,22,23
*SHELL SECTION,MATERIAL=MAT,ELSET=ALLE, SECTION=GAUSS
0.05,5
*MATERIAL,NAME=MAT
*DENSITY
8000.,
*ELASTIC
200.E9,0.3
*DAMPING,ALPHA=0.299,BETA=1.339E-3
*BOUNDARY
NALL,1,2
NALL,6
EDGE,1,3
BOT,5
TOP,5
LEFT,4
RIGHT,4
**
** Harmonic Load (direct)
** ----------------------
**
*STEP
*STEADY STATE DYNAMICS,DIRECT,FREQUENCY SCALE=LINEAR
0.1,4.16,200,,
*DLOAD,LOAD CASE=1
ALLE,P,100.
*NODE PRINT,NSET=NMID
U3, 
*EL PRINT,POSITION=AVERAGED AT NODES,ELSET=EMID
S11,
SF,
*NODE FILE,NSET=NMID
U, 
**
** REMOVE OUTPUT RESTRICTION IN FUTURE QA
**
*OUTPUT,FIELD
*NODE OUTPUT,NSET=NMID
U,
**
**REMOVEOUTPUTRESTRICTIONINFUTUREQA
**
*EL FILE,FREQ=116,POSITION=AVERAGED AT NODES,ELSET=EMID
S,
SF,
*ENDSTEP