*HEADING
 BASIC VERIFICATION - *COUPLING WITH *KINEMATIC SUBOPTION
 EQUIVALENT TO REVOLUTE MPC TYPE
*PREPRINT, MODEL=YES
*NODE
10,0.0,0.0,0.0
20,2.0,0.0,0.0
1,1.0,0.0,0.0
2,1.0,0.0,0.0
105,1.0,0.0,1.0 
*ELEMENT,TYPE=B33H,ELSET=A
2,2,20
*ORIENTATION, NAME=RELEASE_AXIS
1.,0.,1.,-1.,0.,1.
**
** *KINEMATIC COUPLING,REF NODE=10, ORIENT=RELEASE_AXIS
** 2,1,3
** 2,4,5
**
*SURFACE, NAME=NODE_BASED_SURFACE, TYPE=NODE
 2, 1.0
*COUPLING, CONSTRAINT NAME=REVOLUTE, REF NODE=10, SURFACE=NODE_BASED_SURFACE, ORIENTATION=RELEASE_AXIS
*KINEMATIC
 1,3
 4,5
*BEAM SECTION,SECT=CIRC,MATERIAL=MAT,ELSET=A
0.1, 
*MATERIAL,NAME=MAT
*ELASTIC
200.0E9,0.0
*BOUNDARY
10,1,6
*STEP, UNSYMM=YES
 PULL ON THE JOINT
*STATIC
1.0,1.0
*CLOAD
20,1,1.0E6
*MONITOR,NODE=20,DOF=2
*EL PRINT,FREQ=9999
SF, 
*NODE PRINT,FREQ=9999
U,COORD
RF, 
*EL FILE,FREQ=9999
SF, 
*NODE FILE,FREQ=9999
U,COORD
*OUTPUT, FIELD, FREQ=9999
*NODE OUTPUT
 U,COORD
*ELEMENT OUTPUT
 SF,
*END STEP
*STEP,NLGEOM, UNSYMM=YES
  TWIST THE REFERENCE NODE AND THE REVOLUTE
  AXIS 90 DEGREES ABOUT X
*STATIC
0.25,1.0
*BOUNDARY
10,4,4, 1.5707963
*END STEP
*STEP,NLGEOM, UNSYMM=YES
  TWIST THE REFERENCE NODE 45 DEGREES ABOUT Z
*STATIC
0.25,1.0            
*BOUNDARY
10,6,6, 0.785398164
*BOUNDARY,TYPE=VELOCITY
10,4,4,0.0
*END STEP