*HEADING
:NAFEMS FORCED VIBRATION; TEST 5T;  [B32]
(DIRECT INTEGRATION) 
*RESTART,WRITE,F=10
*NODE
1, 0.
11, 10.
*NGEN
1,11
*NSET,NSET=NMID
6, 
*ELEMENT,TYPE=B32
1, 1,2,3
*ELGEN,ELSET=ALLE
1, 5,2,1
*ELSET,ELSET=EMID
3, 
*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
11, 2,3
**
** Transient Load
** --------------
**
*STEP,AMP=STEP,INC=1000
*DYNAMIC,NOHAF
0.0001,0.02
*DLOAD
ALLE,PY,1.E6
*NODE PRINT,NSET=NMID
U2, 
*EL PRINT,FRE=118,POSITION=AVERAGED AT NODES,ELSET=EMID
S11, 
*NODE FILE,NSET=NMID
U, 
*EL FILE,FRE=118,POSITION=AVERAGED AT NODES,ELSET=EMID
S, 
*ENDSTEP
**
** Increase the timestep to find steady state
**
*STEP,AMP=STEP,INC=1000
*DYNAMIC,NOHAF
0.01,1.0
*DLOAD
ALLE,PY,1.E6
*NODE PRINT,NSET=NMID,F=1000
U2, 
*EL PRINT,POSITION=AVERAGED AT NODES,ELSET=EMID,F=1000
S11, 
*NODE FILE,NSET=NMID,F=1000
U, 
*EL FILE,POSITION=AVERAGED AT NODES,ELSET=EMID,F=1000
S, 
*ENDSTEP