*HEADING
fv5t_b22_c: NAFEMS FORCED VIBRATION; TEST 5T;  [B22]
(DIRECT INTEGRATION) 
*RESTART,WRITE,NUM=2,TIMEMARKS=NO
*NODE
1, 0.
7, 10.
*NGEN
1,7,1
*NSET,NSET=NMID
4, 
*ELEMENT,TYPE=B22
1, 1,2,3
*ELGEN,ELSET=ALLE
1, 3,2,1
*ELSET,ELSET=EMID
2,
*BEAM SECTION,SECT=RECT,MATERIAL=MAT,ELSET=ALLE
2.0,2.0
0.,0.,-1.
*MATERIAL,NAME=MAT
*DENSITY
8000.,
*ELASTIC
200.E9,0.3
*DAMPING,ALPHA=5.36,BETA=7.46E-5
*BOUNDARY
1, 1,4
7, 2,3
**
** Transient Load
** --------------
**
*NSET, NSET=CHECK
4, 7 
*NSET,NSET=QA_TEST_CHECK
CHECK,
**
*STEP
*DYNAMIC, EXPLICIT, DIRECT
0.0001,0.02
** ,0.02
*BULK VISCOSITY
1.0E-5,0.02
*DLOAD
ALLE,PY,1.E6
***OUTPUT,HISTORY,TIME INTERVAL=1.0E-3
***NODE OUTPUT,NSET=NMID
**U2,
***ELEMENT OUTPUT,ELSET=EMID
**S11,
***ENERGY OUTPUT
**ALLWK,ALLKE,ALLIE,ALLVD,ETOTAL
*FILE OUTPUT, NUMBER INTER=4,TIMEMARKS=YES
*NODE FILE, NSET=CHECK
U, 
*ENERGY FILE
*OUTPUT, FIELD, NUMBER INTERVAL=4, TIMEMARKS=YES
*NODE OUTPUT, NSET=QA_TEST_CHECK
U, 
*OUTPUT,FIELD,VARIABLE=PRESELECT
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=0.005
*ENDSTEP
**
** Increase the time period to reach steady state
**
*STEP
*DYNAMIC,EXPLICIT
**0.01,1.0
 ,1.0
*DLOAD
ALLE,PY,1.E6
*FILE OUTPUT, NUMBER INTER=4,TIMEMARKS=YES
*NODE FILE, NSET=CHECK
U, 
*ENERGY FILE
*OUTPUT, FIELD, NUMBER INTERVAL=4, TIMEMARKS=YES
*NODE OUTPUT, NSET=QA_TEST_CHECK
U, 
*OUTPUT,FIELD,VARIABLE=PRESELECT
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=0.25
*ENDSTEP