*HEADING
 B21 with *BAI (compare with b31_circ_bai)
*preprint,model=yes,history=yes,echo=yes
*NODE,NSET=NALL
   1,  70.710678, 70.710678
  11,  282.84271247, 282.84271247
  101,  70.710678, 70.710678
  111,  282.84271247, 282.84271247
*NSET,NSET=ENDS
 1,11,101,111 
*NSET,NSET=TIP
 11,111
*ELEMENT,TYPE=B21,ELSET=BS
 1, 1,11
*ELEMENT,TYPE=B21,ELSET=BGS
 2, 101,111
*ELSET,ELSET=EALL
 BS,BGS
*BEAM GENERAL SECTION,SECTION=CIRC,ELSET=BGS,ROTARY INERTIA=EXACT,DENSITY=80.
 3.,
0.,0.,-1.
30.E8
*BEAM ADDED INERTIA
 2.0E3, 0.0, 60.  
*BEAM SECTION,SECTION=CIRC,MATERIAL=A1,ELSET=BS,ROTARY INERTIA=EXACT
 3.,
*BEAM ADDED INERTIA
 2.0E3, 0.0, 60.  
*MATERIAL,NAME=A1
*ELASTIC
 30.E8,  
*DENSITY
 80.,
*BOUNDARY
 1, 1, 6
 101,1,6
*transform,nset=tip
0.70710678,-0.70710678,0., 0.70710678,0.70710678,0.
*INITIAL CONDITION,TYPE=VELOCITY
***velocity of 10 in the local 2 direction
 11,1,7.0710678
 11,2,-7.0710678
 111,1,7.0710678
 111,2,-7.0710678
*STEP,inc=200,nlgeom
 Initial velocity
*DYNAMIC,alpha=-0.05
 .2,40.0
*NODE PRINT,FREQ=50
 U,
 V,
 A,
*NODE PRINT,global=yes,freq=50
 U,
 V,
 A,
*ENERGY FILE,FREQUENCY=50
*OUTPUT,FIELD, FREQUENCY=50
*NODE OUTPUT,NSET=NALL
 U,V,A
*ELEMENT OUTPUT,ELSET=EALL
*OUTPUT,HISTORY, FREQUENCY=50
*NODE OUTPUT,NSET=TIP
 U,V,A
*ENERGY OUTPUT,VARIABLE=ALL
*END STEP