*HEADING
 B31 with *BEAM ADDED INERTIA 
 Compare with b31_circ_mass_ri (mass/rotatryi)
*preprint,model=yes,history=yes,echo=yes
*NODE,NSET=NALL
   1,  0.,   0.
  11,  0., 300.
  101, 0.,   0.
  111, 0., 300.
*NSET,NSET=ENDS
 1,11,101,111 
*NSET,NSET=TIP
 11,111
*ELEMENT,TYPE=B31,ELSET=BS
 1, 1,11
*ELEMENT,TYPE=B31,ELSET=BGS
 101, 101,111
*BEAM GENERAL SECTION,SECTION=CIRC,ELSET=BGS,DENSITY=80.
 3.,
 0.,0.,-1.
 30.E8,
*BEAM ADDED INERTIA
2.0E3, 0.,-60.,0.,666.667,666.667 
2.0E3, -60.,  0.
*BEAM SECTION,SECTION=CIRC,MATERIAL=A1,ELSET=BS,ROTARY INERTIA=EXACT
 3.,
 0.,0.,-1.
*BEAM ADDED INERTIA
2.0E3, 0.,-60.,0.,666.667,666.667 
2.0E3, -60.,  0.
*MATERIAL,NAME=A1
*ELASTIC
 30.E8,  
*DENSITY
 80.,
*BOUNDARY
 1, 1,6
 11,3,3
 101,1,6
 111,3,3
*INITIAL CONDITION,TYPE=VELOCITY
 11,1,10.
 11,6,0.1
 111,1,10.
 111,6,0.1
*ELSET,ELSET=EALL
 BS,BGS
*STEP,INC=200,NLGEOM,UNSYMM=YES
 Initial velocity
*DYNAMIC,alpha=-0.05,direct
 .2,40.0
*NODE PRINT,freq=50
 U,
 V,
 A,
*ENERGY FILE,FREQUENCY=1
*OUTPUT,FIELD, FREQUENCY=1
*NODE OUTPUT,NSET=NALL
 U,V,A
*ELEMENT OUTPUT,ELSET=EALL
*OUTPUT,HISTORY, FREQUENCY=1
*NODE OUTPUT,NSET=TIP
 U,V,A,
*ENERGY OUTPUT,VARIABLE=ALL
*END STEP