*HEADING
 BASIC VERIFICATION - *COUPLING WITH *KINEMATIC SUB OPTION
 EQUIVALENT TO UNIVERSAL 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
0.,1.,0.,0.,0.,1.
**
** *KINEMATIC COUPLING,REF NODE=10,ORIENT=RELEASE_AXIS
** 2,1,3
** 2,6,6
**
*NSET, NSET=NREFNODE
 10
*SURFACE, NAME=SURF, TYPE=NODE
 2,
*COUPLING, CONSTRAINT NAME=UNIVERSAL, REF NODE=NREFNODE, SURFACE=SURF, ORIENTATION=RELEASE_AXIS
*KINEMATIC
 1, 3
 6, 6
*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, 
*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