*HEADING
fv13t_s3r_c: NAFEMS FORCED VIBRATION; TEST 13T; [S3RS]
(DIRECT INTEGRATION)
*RESTART,WRITE,NUM=2,TIMEMARKS=NO
*NODE
1, 0.
5, 5.
41, 0., 5.
45, 5.,5.
*NGEN,NSET=BOT
1,5
*NGEN,NSET=SIDEY
41,45
*NFILL,NSET=NALL
BOT,SIDEY,4,10
*NSET,NSET=LEFT,GEN
1,41,10
*NSET,NSET=EDGE
BOT,LEFT
*NSET,NSET=SIDEX
5,25,45
*NSET,NSET=NMID
45, 
*ELEMENT,TYPE=S3RS
1, 1,12,21
2,1,3,12
3,3,23,12
4,23,21,12
*ELGEN,ELSET=ALLE
1, 2,2,4, 2,20,8
2, 2,2,4, 2,20,8
3, 2,2,4, 2,20,8
4, 2,2,4, 2,20,8
*ELSET,ELSET=EMID
15,16
*SHELL SECTION,MATERIAL=MAT,ELSET=ALLE, POISSON=0.3
0.05, 
*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
LEFT,4
SIDEX,XSYMM
SIDEY,YSYMM
**
** Transient Load
** --------------
**
*STEP
*DYNAMIC,EXPLICIT
 ,0.5
*DLOAD
ALLE,P,100.
***OUTPUT,HISTORY,TIME INTERVAL=0.005
***NODE OUTPUT,NSET=NMID
**U3,
***ELEMENT OUTPUT,ELSET=EMID
**S11,
***ENERGY OUTPUT
**ALLWK,ALLIE,ALLKE,ALLAE,ALLVD,ETOTAL,DT
*OUTPUT,FIELD,VAR=PRESELECT,NUMBER INTERVAL=5
**
** Output requests for qa testing
**
*NSET,NSET=QA_TEST
SIDEX,
*OUTPUT,FIELD,NUMBER INTERVAL=4,TIMEMARKS=YES
*NODE OUTPUT, NSET=QA_TEST
U,
*OUTPUT,HISTORY,VAR=PRESELECT,TIME INTERVAL=0.125
*FILE OUTPUT, NUMBER INTER=4, TIMEMARKS=YES
*NODE FILE, NSET=SIDEX
U, 
*ENERGY FILE
*ENDSTEP
**
** Increase the time period to find steady state
**
*STEP
*DYNAMIC,EXPLICIT
** ,25.0
,5.0
*DLOAD
ALLE,P,100.
*ENDSTEP