*HEADING
FORCED MOTION OF A RIGID BODY IN EXPLICIT
USING CONNECTOR ELEMENTS TYPE SLOT+LINK AND EULER
*RESTART,WRITE
*NODE,NSET=RBODY
*NODE,NSET=FXPNT
****************************************
** BEAM CONNECTOR TO MODEL THE TOP    **
****************************************
*ELEMENT, TYPE=CONN3D2, ELSET=SLOT
*CONNECTOR SECTION, ELSET=SLOT
SLOT,
SLOTORIENT,
*ELEMENT, TYPE=CONN3D2, ELSET=LINK
*CONNECTOR SECTION, ELSET=LINK
LINK,
SLOTORIENT,
*ORIENTATION, NAME=SLOTORIENT
*ELEMENT,TYPE=ROTARYI,ELSET=ROTARY_INERTIA1
*ELEMENT,TYPE=ROTARYI,ELSET=ROTARY_INERTIA2
*ROTARY INERTIA,ELSET=ROTARY_INERTIA1,ORIENTATION=BDYAXES
5.,5.,1.
*ROTARY INERTIA,ELSET=ROTARY_INERTIA2
1E-12,1E-12,1E-12
*ELEMENT,TYPE=MASS,ELSET=TRANSLATIONAL_INERTIA
***MASS,ELSET=TRANSLATIONAL_INERTIA, TYPE=ANISOTROPIC, ORIENTATION=BDYAXES
**2e-8,2e-8,2e-8
*MASS,ELSET=TRANSLATIONAL_INERTIA
2e-8,
*ORIENTATION,NAME=BDYAXES
*INITIAL CONDITIONS, TYPE=VELOCITY
100,4,0.5
100,5,-17.101
100,6,46.9846
*BOUNDARY,TYPE=DISP
FXPNT,1,3
*******************************************
**ADD DUMMY EULER CONNECTOR TO THE TOP   **
** TO MEASURE THE EULER'S ANGLES DIRECTLY**
*******************************************
*NODE,NSET=FIXED_REF
*ORIENTATION, NAME=GLOBAL
*ELEMENT, TYPE=CONN3D2, ELSET=EULER
*CONNECTOR SECTION, ELSET=EULER
EULER,
GLOBAL,BDYAXES
*BOUNDARY,TYPE=DISP
5000,1,6
*******************************************
*NSET,NSET=NALL
RBODY, FXPNT
****************
*STEP
*DYNAMIC,EXPLICIT,DIRECT USER CONTROL
***1E-3,15.0
1E-3,3.0
*BULK VISCOSITY
*CLOAD
 1, 3, -20.0
*OUTPUT,FIELD,VARIABLE=ALL,NUMBER INTERVAL=60
*OUTPUT,HISTORY,TIME INTERVAL=2E-2
*NODE OUTPUT,NSET=NALL
U,V,A
*ELEMENT OUTPUT, ELSET=EULER
CP
*ELEMENT OUTPUT, ELSET=SLOT
CP,CU
**
** Output requests for qa testing
**
*NSET,NSET=QA_TEST
NALL,
*OUTPUT,FIELD,NUMBER INTERVAL=4,TIME MARKS=YES
*NODE OUTPUT, NSET=QA_TEST
U,
*OUTPUT,HISTORY,VAR=PRESELECT,TIME INTERVAL=0.75
*FILE OUTPUT,NUMBER INTERVAL=4, TIME MARKS=YES
*NODE FILE, NSET=NALL
U,
*END STEP