*HEADING
 REVOLUTE MPC VERIFICATION: FORM AND LOAD A CRANK     
*NODE
1,0.0
2,400.0
3,400.0
4,800.0
5,800.0
6,1200.0
12,400.0,1.0
14,800.0,1.0
*ELEMENT,TYPE=B31H,ELSET=EALL
1,1,2
2,3,4
3,5,6
*BEAM SECTION,SECTION=RECT,MATERIAL=MAT,ELSET=EALL
20.0,20.0
*MATERIAL,NAME=MAT
*ELASTIC
200000.0, 
*MPC
PIN,2,3
PIN,4,5
REVOLUTE,2,3,12
REVOLUTE,4,5,14
*BOUNDARY
1,1,6
12,6
14,6
*STEP,PERTURBATION
 PERTURBATION STEP WITH ZERO RELATIVE ROTATION IN THE REVOLUTE JOINT
*STATIC
*CLOAD
6,3,-500.0
*NODE FILE,FREQUENCY=99
U,COORD
CF, 
*OUTPUT,FIELD,FREQ=99
*NODE OUTPUT
U,COORD
CF,
*NODE PRINT,FREQUENCY=99
U,COORD
CF, 
RF, 
*EL PRINT,FREQUENCY=0
*END STEP
*STEP,NLGEOM
*STATIC
*CLOAD
6,3,-500.0
*NODE FILE,FREQUENCY=99
U,COORD
CF, 
*OUTPUT,FIELD,FREQ=99
*NODE OUTPUT
U,COORD
CF,
*NODE PRINT,FREQUENCY=99
U,COORD
CF, 
RF, 
*EL PRINT,FREQUENCY=0
*END STEP
*STEP,NLGEOM
 PRESCRIBED ROTATIONS 90 DEG. IN THE JOINTS
*STATIC
0.25,1.0
*BOUNDARY,OP=NEW
1,1,6
12,6,6,-1.5707963
14,6,6, 1.5707963
*END STEP
*STEP,NLGEOM
 PRESCRIBED ROTATIONS 90 DEG. IN THE JOINTS + 90 DEG. ABOUT THE Z-AXIS
*STATIC
0.25,1.0
*BOUNDARY,OP=NEW
1,1,5
1,6,6,  1.5707963
12,6,6,-1.5707963
14,6,6, 1.5707963
*END STEP