*HEADING
 BASIC VERIFICATION - UNIVERSAL MPC TYPE
 ELEMENT: B33H, TEST TYPE:5 ,TEST NUMBER:6  --  MPCUNIVERSAL
*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
 3,1.0,0.0,1.0
 4,1.0,1.0,0.0 
*NODE,NSET=CUBE
 11, -1.,-1.,-1.
 12,  0.,-1.,-1.
 13,  0., 0.,-1.
 14, -1., 0.,-1.
 15, -1.,-1., 0.
 16,  0.,-1., 0.
 18, -1., 0., 0.
*ELEMENT,TYPE=C3D8P,ELSET=CUBE
 10,11,12,13,14,15,16,10,18
*SOLID SECTION,ELSET=CUBE,MATERIAL=MAT1
*MATERIAL,NAME=MAT1
*ELASTIC
 200.0E9,0.0
*PERMEABILITY,SPECIFIC=1.
 1.,
*INITIAL CONDITIONS,TYPE=RATIO
 CUBE,0.5
*ELEMENT,TYPE=B33H,ELSET=A
 1,10,1
 2,2,20
*MPC
 PIN,1,2
 UNIVERSAL,2,1,3,4
*BEAM SECTION,SECT=CIRC,MATERIAL=MAT,ELSET=A
 0.1, 
*MATERIAL,NAME=MAT
*ELASTIC
 200.0E9,0.0
*BOUNDARY
 3,6,6
 4,6,6
 10,1,6
 CUBE,1,3
*STEP, UNSYMM=YES
 PULL ON THE JOINT
*STATIC
 1.0,1.0
*CLOAD,FOLLOWER
 20,1,1.0E6
*MONITOR,NODE=20,DOF=2
*EL PRINT,ELSET=A,FREQ=9999
 SF, 
*NODE PRINT,FREQ=9999
 U,COORD
 RF, 
*EL FILE,ELSET=A,FREQ=9999
 SF, 
*NODE FILE,FREQ=9999
 U,COORD
*OUTPUT,FIELD
*NODE OUTPUT
 U,COORD
*END STEP
*STEP,NLGEOM,INC=20, UNSYMM=YES
  TWIST THE 3-4 JOINT 45 DEGREES
*STATIC
 0.25,1.0         
*BOUNDARY
 3,6,6, 0.785398164
*END STEP
*STEP,NLGEOM,INC=20, UNSYMM=YES
  TWIST THE 3-5 JOINT 45 DEGREES
*STATIC
 0.25,1.0            
*BOUNDARY
 4,6,6,-0.785398164
*END STEP