*HEADING

CAP PLASTICITY, HYDROSTATIC CYCLIC TEST, C3D8, DISPLACEMENT CONTROL
*NODE, NSET=BOT
1,0.,0.,0.
2,1.,0.,0.
3,1.,1.,0.
4,0.,1.,0.
*NODE, NSET=TOP
5,0.,0.,1.
6,1.,0.,1.
7,1.,1.,1.
8,0.,1.,1.
*NSET, NSET=ALLN
BOT, TOP
*NSET, NSET=BACK
1,4,5,8
*NSET, NSET=FRONT
2,3,6,7
*NSET, NSET=LHS
1,2,5,6
*NSET, NSET=RHS
3,4,7,8
*ELEMENT, TYPE=C3D8, ELSET=EL1
1,1,2,3,4,5,6,7,8
*Distribution Table, name=distTableElastic
 modulus, ratio
*Distribution, name=distElastic, table=distTableElastic, location=element
 , 200.e3, 0.25
 EL1, 12857.1429, 0.28571429
*SOLID SECTION, ELSET=EL1, MATERIAL=CAPL
*MATERIAL,NAME=CAPL
** The material data is chosen to give G = 5000.0, BULK = 10000.0,
** first yield in pure shear        at P =   0.0,  S12 = 100.0, 
** first yield in hydro compression at P = 270.0,  S12 =   0.0,
** first yield in hydro tension     at P =-300.0,  S12 =   0.0,
** first yield (P = PA)             at P = 120.0,  S12 = 125.0
*ELASTIC
 distElastic,
*CAP PLASTICITY
** COHESN, BETA,    BIGR,   EVOLPI,   ALPHA,   BIGK
173.20508, 30.0, 0.61858957, 0.027, 0.69258232, 1.0
*CAP HARDENING
** PB, EVOLP
 213.0, 0.00
 222.0, 0.01
 242.0, 0.02
 282.0, 0.03
 362.0, 0.04
 522.0, 0.05
 842.0, 0.06
1482.0, 0.07                                              
2762.0, 0.08
**
*STEP, INC=8, UNSYMM=YES
Step 1: Load to yield in compression
*STATIC, DIRECT
1.0, 8.0
*BOUNDARY
BACK,1
LHS, 2
BOT, 3
FRONT,1,1, -0.024
RHS,  2,2, -0.024
TOP,  3,3, -0.024
*EL PRINT, F=999
S, 
SINV, 
E, 
PE, 
PEQC, 
*EL PRINT, POS=NODES, F=999
PEQC4, 
*NODE PRINT, F=999
U, 
RF   , 
*OUTPUT,FIELD,FREQ=999
*ELEMENT OUTPUT
PE,S,E,PEQC
*NODE OUTPUT
U,RF
*END STEP
**
*STEP, INC=2, UNSYMM=YES
Step 2: Unload, still in compression
*STATIC, DIRECT
1.0, 2.0
*BOUNDARY
FRONT,1,1, -0.018
RHS,  2,2, -0.018
TOP,  3,3, -0.018
*END STEP
**
*STEP, INC=4, UNSYMM=YES
Step 3: Reload in compression
*STATIC, DIRECT
1.0, 4.0
*BOUNDARY
FRONT,1,1, -0.030
RHS,  2,2, -0.030
TOP,  3,3, -0.030
*END STEP
**
*STEP, INC=20, UNSYMM=YES
Step 4: Unload in compression and load to yield in tension
*STATIC, DIRECT
1.0, 20.0
*BOUNDARY
FRONT,1,1,  0.024
RHS,  2,2,  0.024
TOP,  3,3,  0.024
*END STEP
**
*STEP, INC=20, UNSYMM=YES
Step 5: Unload in tension and load to yield in compression
*STATIC, DIRECT
1.0, 20.0
*BOUNDARY
FRONT,1,1, -0.024
RHS,  2,2, -0.024
TOP,  3,3, -0.024
*END STEP
**
*STEP, INC=2, UNSYMM=YES
Step 6: Unload
*STATIC, DIRECT
1.0, 2.0
*BOUNDARY
FRONT,1,1, -0.018
RHS,  2,2, -0.018
TOP,  3,3, -0.018
*END STEP