*HEADING
:  Rigid body rotation about end point for 3-D beam.
*NODE,NSET=NALL
1,0.,0.
2,.707107,.707107
21,0.,0.
22,.707107,.707107
23,.3535535,.3535535
*ELEMENT,TYPE=B31,ELSET=EALL
1,1,2
*ELEMENT,TYPE=B32,ELSET=EALL2
21,21, 23, 22
*ELSET, ELSET=EALL_ALL
EALL, EALL2
*BEAM SECTION,SECTION=CIRC,MATERIAL=STEEL,ELSET=EALL,POISSON=.3
.02, 
0.,.5,-.866025
*BEAM SECTION,SECTION=CIRC,MATERIAL=STEEL,ELSET=EALL2,POISSON=.3
.02, 
0.,.5,-.866025
*MATERIAL,NAME=STEEL
*ELASTIC
2.1E7,.3
*DENSITY
7800., 
*BOUNDARY
1,1,4
21,1,4
***the following BC should not be 
** needed; no twist should take place
**1,6,6
**21,6,6
*TRANSFORM,NSET=NALL
.707107,.707107,0.,-.707107,.707107,0.
** Note that the initial velocity must be specified
** in global system even when transform is defined
** at the nodes.
*INITIAL CONDITIONS,TYPE=VELOCITY
NALL,4,-3.668997
NALL,5, 3.668997
2,3,-5.188745
22,3,-5.188745
23,3,-2.5943725
*INITIAL CONDITIONS,TYPE=STRESS
1,1.0524E5
21,1.0524E5
**
*NSET, NSET=QA_TEST_ALLNODES
NALL, 
**
*STEP
*DYNAMIC,EXPLICIT
,6.0545
*BULK VISCOSITY
0.,0.
*RESTART,WRITE,NUMBER INTERVAL=1
***OUTPUT,HISTORY,TIME INTERVAL=0.0
***ELEMENT OUTPUT,ELSET=EALL_ALL
**SE,S11
***NODE OUTPUT,NSET=NALL
**U,V,UR1,UR2,UR3,VR1,VR2,VR3
***ENERGY OUTPUT
**ALLKE,ALLSE,ETOTAL
*FILE OUTPUT,NUMBER=4, TIMEMARKS=YES
*NODE FILE
U,
*OUTPUT, FIELD, NUMBER INTERVAL=4, TIMEMARKS=YES
*NODE OUTPUT, NSET=QA_TEST_ALLNODES
U,
*OUTPUT,HISTORY,VARIABLE=PRESELECT,TIMEINT=6.0545e-1
*END STEP
*STEP
*DYNAMIC,EXPLICIT
,6.0545
*BULK VISCOSITY
0.06,1.2
*END STEP