*HEADING

CAP PLASTICITY, UNI-AXIAL COMPRESSIVE STRAIN TEST, CPE4, LOAD CONTROL
*NODE, NSET=ALLN
1,0.,0.
2,1.,0.
3,1.,1
4,0.,1.
*ELEMENT, TYPE=CPE4, ELSET=EL1
1, 1,2,3,4
*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, DEPENDENCIES=1
** COHESN, BETA,    BIGR,   EVOLPI,   ALPHA,   BIGK,  TEMP, FV
173.20508, 25.0, 0.61858957, 0.027, 0.69258232, 1.0,   0.0,  0.0
173.20508, 35.0, 0.61858957, 0.027, 0.69258232, 1.0, 200.0,  0.0
173.20508, 25.0, 0.61858957, 0.024, 0.69258232, 1.0,   0.0, 40.0
173.20508, 35.0, 0.61858957, 0.030, 0.69258232, 1.0, 200.0, 40.0
*CAP HARDENING
** PB, EVOLP, TEMP
 223.0, 0.00,  0.0
 232.0, 0.01,  0.0
 252.0, 0.02,  0.0
 292.0, 0.03,  0.0
 372.0, 0.04,  0.0
 532.0, 0.05,  0.0
 852.0, 0.06,  0.0
1492.0, 0.07,  0.0
2772.0, 0.08,  0.0
 203.0, 0.00, 200.0
 212.0, 0.01, 200.0
 232.0, 0.02, 200.0
 272.0, 0.03, 200.0
 352.0, 0.04, 200.0
 512.0, 0.05, 200.0
 832.0, 0.06, 200.0
1472.0, 0.07, 200.0
2752.0, 0.08, 200.0
**
*INITIAL CONDITIONS, TYPE=TEMPERATURE
ALLN, 100.0
*INITIAL CONDITIONS, TYPE=FIELD, VARIABLE=1
ALLN, 20.0
**
*STEP, INC=10, UNSYMM=YES
*STATIC, DIRECT
1., 10.
*BOUNDARY
ALLN, 2
1, 1
4, 1
*DLOAD
EL1, P2, 800.0
*EL PRINT, F=5
S, 
SINV, 
E, 
PE, 
PEQC, 
*EL PRINT, POS=NODES, F=5
PEQC4, 
*NODE PRINT, F=5
U, 
RF   , 
*OUTPUT,FIELD,FREQ=5
*ELEMENT OUTPUT
PE,S,E,PEQC
*NODE OUTPUT
U,RF
*END STEP