*HEADING
 DYNAMIC ANALYSIS COMPARISON  B31                  
 (exact rotary inertia formulation versus isotropic:
  b31_moddyn_iso.inp))                 
*NODE
1 ,, 300.
11, 
*NSET,NSET=ENDS
1,11
*ELSET,ELSET=PR
1, 
*ELEMENT,TYPE=B31,ELSET=EALL
1,1,11
*MATERIAL,NAME=A1
*ELASTIC
30.E6, 
*DENSITY
 .000728, 
*BEAM SECTION,SECTION=RECT,MATERIAL=A1,ELSET=EALL,ROTARY INERTIA=EXACT
1. , 2.
*BOUNDARY
11, 1, 6
*AMPLITUDE,TIME=STEP TIME,VALUE=RELATIVE,NAME=EQ
0.0,0.0,0.25,1.0,0.5,0.0,0.75,-1.0,
1.0,0.0
**
*STEP
*FREQUENCY,EIGENSOLVER=LANCZOS, SIM=NO
7,
*MODAL FILE
*END STEP
*STEP
 SINE TYPE INPUT
*MODAL DYNAMIC
.025,1.0
*BASE MOTION,DOF=1,AMPLITUDE=EQ,SCALE= 386.09
*PRINT,RESIDUAL=NO,FREQUENCY=50
*EL PRINT,ELSET=PR,FREQUENCY=50
S,E
*NODE PRINT,FREQUENCY=50
U,
TU,
*NODE FILE,NSET=ENDS
TU,TV,TA
*ENERGY FILE,FREQUENCY=1
*OUTPUT,FIELD, FREQUENCY=1,OP=NEW
*NODE OUTPUT,NSET=ENDS
U,V,A
*OUTPUT,HISTORY, FREQUENCY=1,OP=NEW
*NODE OUTPUT,NSET=ENDS
U,V,A,
TU,TV,TA,
*ENERGY OUTPUT,VARIABLE=ALL
*END STEP