*HEADING
: NAFEMS FORCED VIBRATION; TEST 5 [B32]
*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
*BOUNDARY
1, 1,4
11, 2,3
**
** Modes
** -----
**
*STEP,PERT
*FREQUENCY, eigensolver=lanczos
16, , 
*NODE PRINT,F=999
*EL PRINT,F=999
*MODAL FILE,F=1
*output,field
*node output
u,
*ENDSTEP
**
** Harmonic Load
** -------------
**
*STEP,PERT
*STEADY STATE DYNAMICS,FREQUENCY SCALE=LINEAR
20.,70.,50
*DLOAD,LOAD CASE=1
ALLE,PY,1.E6
*MODAL DAMPING,MODAL=DIRECT
1,16,0.02
*SELECT EIGENMODES, GENERATE
1,16,1
*NODE PRINT,NSET=NMID
U2, 
*EL PRINT,POSITION=AVERAGED AT NODES,ELSET=EMID
S11, 
*NODE FILE,NSET=NMID
U, 
*EL FILE,FRE=50,POSITION=AVERAGED AT NODES,ELSET=EMID
S, 
*output,field,frequency=50
*node output
u,
*element output
s,
*END STEP
**
** Transient Load
** --------------
**
*STEP,PERT,AMP=STEP
*MODAL DYNAMIC
0.0001,0.02
*DLOAD
ALLE,PY,1.E6
*MODAL DAMPING,MODAL=DIRECT
1,16,0.02
*SELECT EIGENMODES, GENERATE
1,16,1
*NODE PRINT,NSET=NMID
U2, 
*EL PRINT,POSITION=AVERAGED AT NODES,ELSET=EMID
S11, 
*NODE FILE,NSET=NMID
U, 
*EL FILE,FRE=116,POSITION=AVERAGED AT NODES,ELSET=EMID
S, 
*output,field,frequency=107
*node output,nset=nmid
u,
*output,history,frequency=107
*node output,nset=nmid
 U2,
*ENDSTEP
**
** Increase the damping
** (Note: an intermediate step was inserted here to see
**   detail of the velocity continuation, just after change
**   of step).
**
*STEP,PERT,AMP=STEP
*MODAL DYNAMIC,CONTINUE=YES
0.0001,0.02
*DLOAD
ALLE,PY,1.E6
*MODAL DAMPING,MODAL=DIRECT
1,16,0.10
*SELECT EIGENMODES, GENERATE
1,16,1
*NODE PRINT,NSET=NMID,F=100
U2, 
*EL PRINT,POSITION=AVERAGED AT NODES,ELSET=EMID,F=999
S11, 
*NODE FILE,NSET=NMID,F=100
U, 
*EL FILE,POSITION=AVERAGED AT NODES,ELSET=EMID,F=999
S, 
*output,field,frequency=107
*node output
u,
*output,history,frequency=107
*node output,nset=nmid
 U2,
*ENDSTEP
**
** Increase the timestep to find steady state
**
**
*STEP,PERT,AMP=STEP
*MODAL DYNAMIC,CONTINUE=YES
0.001,1.0
*DLOAD
ALLE,PY,1.E6
*MODAL DAMPING,MODAL=DIRECT
1,16,0.10
*SELECT EIGENMODES, GENERATE
1,16,1
*NODE PRINT,NSET=NMID,F=100
U2, 
*EL PRINT,POSITION=AVERAGED AT NODES,ELSET=EMID,F=999
S11, 
*NODE FILE,NSET=NMID,F=100
U, 
*EL FILE,POSITION=AVERAGED AT NODES,ELSET=EMID,F=999
S, 
*output,field,frequency=107
*node output
u,
*output,history,frequency=107
*node output,nset=nmid
 U2,
*ENDSTEP