*HEADING
 B31 with mass and rotary inertia elements
 Compare with b31_circ_bai.inp
*preprint,model=yes,history=yes,echo=yes
*NODE,NSET=NALL
   1,  0.,   0.
  11,  0., 300.
 101,  0., 0.,60.
 102, 60., 0.
 111,  0., 300.,60.
 112, 60., 300.
*NSET,NSET=ENDS
 1,11 
*NSET,NSET=TIP
 11,
*ELEMENT,TYPE=B31,ELSET=BEAM
 1, 1,11
*BEAM SECTION,SECTION=CIRC,MATERIAL=A1,ELSET=BEAM,ROTARY INERTIA=EXACT
 3.,
 0.,0.,-1.
*MATERIAL,NAME=A1
*ELASTIC
 30.E8,  
*DENSITY
 80.,
*ELEMENT,TYPE=MASS,ELSET=MASS
 111,111
 112,112
 101,101
 102,102
*MASS,ELSET=MASS
 3.E5
*ELEMENT,TYPE=ROTARYI,ELSET=ROTARYI
 212,112
 202,102
*ROTARY INERTIA,ELSET=ROTARYI
100000.,200000.,100000.
*BOUNDARY
 1, 1, 6
*INITIAL CONDITION,TYPE=VELOCITY
 11,1,10.
 11,6,0.1
*MPC
 BEAM,101,1
 BEAM,102,1
 BEAM,111,11
 BEAM,112,11
*ELSET,ELSET=EALL
 BEAM, MASS
**STEP,inc=200,NLGEOM,UNSYMM=YES
*STEP,inc=200,nlgeom
 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