*HEADING
UNIAXIAL COMPRESSION TEST, 2 JOINTS, ALFA=20 THETA=0   
WITH STATIC LINEAR PERTURBATIONS
*NODE,NSET=ALLN
1,0.,0.,0.
2,1.,0.,0.
3,1.,1.,0.
4,0.,1.,0.
5,0.,0.,1.
6,1.,0.,1.
7,1.,1.,1.
8,0.,1.,1.
*ELEMENT,TYPE=C3D8,ELSET=ALLE
1,1,2,3,4,5,6,7,8
*SOLID SECTION,ELSET=ALLE,MATERIAL=ALLE
*MATERIAL,NAME=ALLE
*ELASTIC
300.E3,.3
*JOINTED MATERIAL
45.,45.,8000.
*JOINTED MATERIAL,JOINT DIRECTION=JOINT1
45.,45.,1000.
*JOINTED MATERIAL,JOINT DIRECTION=JOINT2
45.,45.,1000.
*ORIENTATION,NAME=JOINT1
-.9397,.342,0.,-.342,-.9397,0.
*ORIENTATION,NAME=JOINT2
.9397,.342,0.,-.342,.9397,0.
*STEP,INC=10
STEP 1 - NONLINEAR
*STATIC,DIRECT
1.,10.
*BOUNDARY
1,PINNED
2,2
5,2
6,2
4,1
5,1
8,1
2,3
3,3
4,3
7,2,,-.1
3,2,,-.1
4,2,,-.1
8,2,,-.1
*EL PRINT
S, 
SINV, 
E, 
PE, 
PEQC, 
EE, 
*NODE PRINT
U,RF
*EL FILE,FREQUENCY=2
S,E,PE
*END STEP
*STEP,PERTURBATION,INC=1
STEP 2 - LINEAR PERTURBATION
*STATIC,DIRECT
.1,.1
*BOUNDARY
1,PINNED
2,2
5,2
6,2
4,1
5,1
8,1
2,3
3,3
4,3
7,2,,-.01
3,2,,-.01
4,2,,-.01
8,2,,-.01
*EL PRINT
S, 
SINV, 
E, 
*NODE PRINT
U,RF
*EL FILE,FREQUENCY=2
S,E
*END STEP
*STEP,INC=10
STEP 3 - REST OF NONLINEAR
*STATIC,DIRECT
1.,10.
*BOUNDARY
7,2,,-.2
3,2,,-.2
4,2,,-.2
8,2,,-.2
*END STEP
*STEP,PERTURBATION,INC=1
STEP 4 - LINEAR PERTURBATION
*STATIC,DIRECT
.1,.1
*BOUNDARY
1,PINNED
2,2
5,2
6,2
4,1
5,1
8,1
2,3
3,3
4,3
7,2,,-.02
3,2,,-.02
4,2,,-.02
8,2,,-.02
*EL PRINT
S, 
SINV, 
E, 
*NODE PRINT
U,RF
*EL FILE,FREQUENCY=2
S,E
*END STEP