*HEADING
fv21t_s3r_vf: NAFEMS FORCED VIBRATION; TEST 21T; [S3R]
(DIRECT INTEGRATION)
*RESTART,WRITE,NUM=2, TIMEMARKS=NO
*NODE
1, 0.
5, 5.
41,0.,5.
45, 5.,5.
6,0.625,0.625
9,4.375,0.625
36,0.625,4.375
39,4.375,4.375
*NGEN,NSET=BOT
1,5
*NGEN,NSET=SIDEY
41,45
*NFILL,NSET=NALL
BOT,SIDEY,4,10
*NGEN,NSET=B
6,9,1
*NGEN,NSET=T
36,39,1
*NFILL, NSET=NALL
B,T,3,10
*NSET,NSET=LEFT,GEN
1,41,10
*NSET,NSET=EDGE
BOT,LEFT
*NSET,NSET=SIDEX,GEN
5,45,10
*NSET,NSET=NMID
45, 
*ELEMENT,TYPE=S3R
1,1,6,11
2, 1,2,6
3, 2,12,6
4,12,11,6
*ELGEN,ELSET=ALLE
1, 4,1,4, 4,10,16
2, 4,1,4, 4,10,16
3, 4,1,4, 4,10,16
4, 4,1,4, 4,10,16
*ELSET,ELSET=EMID
63,64
*SHELL SECTION,MATERIAL=MAT,ELSET=ALLE
1.0, 
*MATERIAL,NAME=MAT
*DENSITY
8000.,
*ELASTIC
200.E9,0.3
*DAMPING,ALPHA=5.772,BETA=6.929E-5
*BOUNDARY
NALL,1,2
NALL,6
EDGE,1,3
BOT,5
LEFT,4
SIDEX,XSYMM
SIDEY,YSYMM
**
** Transient Load
** --------------
**
*NSET,NSET=QA_TEST_SIDEX
SIDEX,
**
*STEP
*DYNAMIC,EXPLICIT
**0.0001,0.05
 ,0.05
*DLOAD
ALLE,P,1.E6
***OUTPUT,HISTORY,TIME INTERVAL=0.5E-3
***NODE OUTPUT,NSET=NMID
**U,
***ELEMENT OUTPUT,ELSET=EMID
**S,
***ENERGY OUTPUT
**ALLWK,ALLIE,ALLKE,ALLAE,ALLVD,ETOTAL,DT
*FILE OUTPUT,NUM=4, TIMEMARKS=YES
*NODE FILE, NSET=SIDEX
U, 
*ENERGY FILE
*OUTPUT, FIELD, NUMBER INTERVAL=4, TIMEMARKS=YES
*NODE OUTPUT, NSET=QA_TEST_SIDEX
U, 
*OUTPUT,FIELD,VARIABLE=PRESELECT
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=0.0125
*ENDSTEP
**
** Increase the damping and timestep to find steady state
**
*STEP
*DYNAMIC,EXPLICIT
,1.0
** ,2.0
*DLOAD
ALLE,P,1.E6
*FILE OUTPUT,NUM=4, TIMEMARKS=YES
*NODE FILE, NSET=SIDEX
U, 
*ENERGY FILE
*OUTPUT, FIELD, NUMBER INTERVAL=4, TIMEMARKS=YES
*NODE OUTPUT, NSET=QA_TEST_SIDEX
U, 
*OUTPUT,FIELD,VARIABLE=PRESELECT
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=0.25
*ENDSTEP