*HEADING
Forced Motion of a Rigid Body - Standard  
*PARAMETER
theta=20.0
L=1.0
pi = 3.14159
thetaRad = theta*pi/180. 
z2=L*cos(thetaRad)
y2=-L*sin(thetaRad)
y2p=-y2
load=-20.0
*NODE, NSET=NREF
998, 0.0, 0.0, 0.0
*NODE,NSET=fxpnt
1,0.0,0.0,0.0
*****Parametric node definition*****
*node,nset=rbodycm
2,0.0,<y2>,<z2>
*ELEMENT,TYPE=RB3D2,ELSET=LINK
1,1,2
*ELEMENT,TYPE=ROTARYI,ELSET=ROTARY_INERTIA
10,998
*ELEMENT,TYPE=MASS,ELSET=TRANSLATIONAL_INERTIA
100,998
*MASS,ELSET=TRANSLATIONAL_INERTIA
1e-8, 
*ROTARY INERTIA,ELSET=ROTARY_INERTIA,ORIENTATION=BDYAXES
5.,5.,1.
********************************************
** Dummy element - Postprocessing Purposes *
********************************************
*NODE
3,0.2,<y2>,<z2>
4,-0.2,<y2>,<z2>
*element,type=RB3D2,elset=link
101,2,3
102,2,4
*RIGID BODY,ELSET=link,REF NODE=998
********************************************
*ORIENTATION,NAME=BDYAXES
1,0.0,0.0,0.0,<z2>,<y2p>
*INITIAL CONDITIONS,TYPE=VELOCITY
998,4,0.5
998,5,-17.101
998,6,46.9846
*BOUNDARY
998,1,3,0.0
*nset,nset=nall
 rbodycm, fxpnt
*STEP,INC=3000,NLGEOM, UNSYMM=YES
*DYNAMIC,ALPHA=0.,INITIAL=NO,DIRECT
1e-3,3.0
*** 1e-3,15.0
*cload
2,3,<load>
*output,field,variable=all,freq=50
*output,history,freq=20
*energy output,variable=all
*node output,nset=nall
 u,v,a
*el print,freq=100
 hc, ho,
 vc,
*energy print,freq=100
*node print,freq=100
 u,
 v,
 a,
 rf,
*node file,nset=nall,freq=50
u,v,a
*print,frequency=100
*controls,parameter=field
 1.e-4,
*end step