*HEADING
 B21H with *BEAM ADDED INERTIA 
 1. Frequency extration
 2. Static rigid body rotation
 3. Frequency extraction
*preprint,model=yes,history=yes,echo=yes
*NODE,NSET=NALL
   1,  0.,   0.
  11,  0., 300.
*NSET,NSET=TIP
 11,
*ELEMENT,TYPE=B21H,ELSET=BEAM
 1, 1,11
*BEAM SECTION,SECTION=CIRC,MAT=MAT,
ELSET=BEAM,ROTARY INERTIA=EXACT
 3.,
 0.,0.,-1.
*BEAM ADDED INERTIA
 6.0E3, 0.,-20., 0., 2.4E6
*MATERIAL,NAME=MAT
*ELASTIC
 30.E8, 
*DENSITY
 80. 
*BOUNDARY
 1, 1,2
*ELSET,ELSET=EALL
 BEAM,
*BOUNDARY
 1, 6,6
*STEP,NLGEOM
 FREQ 1
*FREQUENCY, SIM=NO
 3,
*END STEP
*STEP,inc=2000,nlgeom
 rotate
*STATIC
 .1,1.0
*NODE PRINT
 U,
*EL PRINT
 SF,
*BOUNDARY
 1, 6,6, -1.57079633
***********************
*OUTPUT,FIELD, FREQUENCY=0
*NODE OUTPUT,NSET=NALL
 U,
*ELEMENT OUTPUT,ELSET=EALL
*OUTPUT,HISTORY, FREQUENCY=1
*NODE OUTPUT,NSET=TIP
 U,
**************************
*CONTROLS,PARAMETER=FIELD
 1.E-8,
*END STEP
*STEP,NLGEOM
 FREQ
*FREQUENCY, SIM=NO
 3,,,-1.
*BOUNDARY
 1, 6,6
*END STEP