*HEADING

CAP PLASTICITY, TEST ADJUSTMENT OF INITIAL CAP POSITION, 2 C3D8s,
ALL DOFS FIXED
*NODE, NSET=BOT
1,0.,0.,0.
2,1.,0.,0.
3,1.,1.,0.
4,0.,1.,0.
*NODE, NSET=MID
5,0.,0.,1.
6,1.,0.,1.
7,1.,1.,1.
8,0.,1.,1.
*NODE, NSET=TOP
 9,0.,0.,2.
10,1.,0.,2.
11,1.,1.,2.
12,0.,1.,2.
*NSET, NSET=ALLN
BOT, MID, TOP
*ELEMENT, TYPE=C3D8R, ELSET=ELS
1,1,2,3,4,5,6,7,8
2,5,6,7,8,9,10,11,12
*SOLID SECTION, ELSET=ELS, 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
12857.1429, 0.28571429
*CAP PLASTICITY
** COHESN, BETA,    BIGR,   EVOLPI,   ALPHA,   BIGK
173.20508, 30.0, 0.61858957, 0.0,  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
**
*INITIAL CONDITIONS, TYPE=RATIO
ALLN, 1.0
*INITIAL CONDITIONS, TYPE=STRESS
** Element 1 : (P,T)=(119.0, 125.0*sqrt(3.0))
1, -263.33757, -46.83122, -46.83122,  0.0, 0.0, 0.0
** Element 2 : (P,T)=(121.0, 125.0*sqrt(3.0))
2, -121.0, -121.0, -121.0, 125.0, 0.0, 0.0
**
*STEP, NLGEOM, INC=1
*STATIC
1.0, 1.0, 1.0
*BOUNDARY
ALLN,1,3
*EL PRINT
S, 
SINV, 
E, 
EE, 
PE, 
PEQC, 
*EL PRINT, POS=NODES
PEQC4, 
*NODE PRINT
U, 
RF   , 
*OUTPUT,FIELD
*ELEMENT OUTPUT
S,SINV,E,EE,PE
*OUTPUT,FIELD
*NODE OUTPUT
U,RF
*END STEP