*HEADING
Forced Motion of a Rigid Body in STD
using connector elements type BEAM  and EULER
*PARAMETER
theta=20.0
L=1.0
pi = 3.141592654
thetaRad = theta*pi/180. 
z2=L*cos(thetaRad)
y2=-L*sin(thetaRad)
y2p=-y2
load=-20.0
*NODE,NSET=fxpnt
1,0.0,0.0,0.0
*node,nset=rbodycm
2,0.0,<y2>,<z2>
****************************************
** BEAM CONNECTOR TO MODEL THE TOP    **
****************************************
*ELEMENT, TYPE=CONN3D2, ELSET=BEAM
2000,1,2
*CONNECTOR SECTION, ELSET=BEAM
BEAM,
*ELEMENT,TYPE=ROTARYI,ELSET=ROTARY_INERTIA1
10,1
*ELEMENT,TYPE=ROTARYI,ELSET=ROTARY_INERTIA2
20,2
*ELEMENT,TYPE=MASS,ELSET=TRANSLATIONAL_INERTIA
100,1
200,2
*MASS,ELSET=TRANSLATIONAL_INERTIA
1e-8, 
*ROTARY INERTIA,ELSET=ROTARY_INERTIA1,ORIENTATION=BDYAXES
5.,5.,1.
*ROTARY INERTIA,ELSET=ROTARY_INERTIA2
1e-12,1e-12,1e-12
*ORIENTATION,NAME=BDYAXES
1,0.0,0.0,0.0,<z2>,<y2p>
*INITIAL CONDITIONS,TYPE=VELOCITY
1,4,0.5
1,5,-17.101
1,6,46.9846
*BOUNDARY
1,1,3,0.0
*nset,nset=nall
 rbodycm, fxpnt
*******************************************
**ADD DUMMY EULER CONNECTOR TO THE TOP   **
** TO MEASURE THE EULER's ANGLES DIRECTLY**
*******************************************
*NODE,NSET=FIXED_REF
5000,0.0,0.0,0.0
*ORIENTATION, NAME=GLOBAL
1.,0.,0., 0.,1.,0.
*ELEMENT, TYPE=CONN3D2, ELSET=EULER
1000,5000,2
*CONNECTOR SECTION, ELSET=EULER 
euler,
GLOBAL,BDYAXES
*BOUNDARY
5000,1,6
*STEP,INC=3000,NLGEOM, UNSYMM=YES
*DYNAMIC,ALPHA=0.,INITIAL=NO,DIRECT
1e-3,3.0
*** 1e-3,15.0
*cload
2,3,<load>
*controls,parameter=field
 1.e-4,
**
**
***output,field,variable=all,freq=50
*output,history,freq=20
*energy output,variable=all
*node output,nset=nall
u,v,a
*node file,nset=nall,freq=50
u,v,a
*el file,elset=euler,freq=50
CP
*print,frequency=100
*OUTPUT,HISTORY, FREQ=20
*ELEMENT OUTPUT, ELSET=EULER
CP
*ELEMENT OUTPUT, ELSET=BEAM
CP,CU
*end step