*HEADING
 Simple test of cohesive behavior with damamge in standard [3D]
*PREPRINT, MODEL=YES,HISTORY=YES,ECHO=YES
**
** CASE1:
**
*NODE, NSET=BLOCK1A
 1, 0.0,0.0,0.0
 2, 0.0,1.0,0.0
 3,-1.0,1.0,0.0
 4,-1.0,0.0,0.0
11, 0.0,0.0,1.0
12, 0.0,1.0,1.0
13,-1.0,1.0,1.0
14,-1.0,0.0,1.0
**
*ELEMENT, TYPE=C3D8R, ELSET=BLOCK1A
 1, 1,2,3,4,11,12,13,14
**
*NODE, NSET=BLOCK1B
101, 1.0,0.0,0.0
102, 1.0,1.0,0.0
103, 0.0,1.0,0.0
104, 0.0,0.0,0.0
111, 1.0,0.0,1.0
112, 1.0,1.0,1.0
113, 0.0,1.0,1.0
114, 0.0,0.0,1.0
*ELEMENT, TYPE=C3D8R, ELSET=BLOCK1B
101, 101, 102, 103, 104, 111, 112, 113, 114
**
*ELSET,ELSET=BLOCK1
BLOCK1A,BLOCK1B
**
*SOLID SECTION,ELSET=BLOCK1,MAT=MAT1
1., 
*MATERIAL,NAME=MAT1
*ELASTIC
1.5E11,.3
*NSET, NSET=BLOCK1
BLOCK1A,BLOCK1B
*SURFACE, NAME=SECONDARY1
BLOCK1A,S3
*SURFACE, NAME=MAIN1
BLOCK1B,S5
**
*SURFACE INTERACTION, NAME=INTPROP1
*COHESIVE BEHAVIOR, Eligibility = Current Contacts 
*DAMAGE INITIATION, CRITERION=MAXS
1.E+10,2.E+12,2.E+12
*DAMAGE EVOLUTION, TYPE=DISPLACEMENT, SOFTENING=LINEAR
.002,
*DAMAGE STABILIZATION
1.E-6
**
*CONTACT PAIR, INTERACTION=INTPROP1
SECONDARY1, MAIN1
**
** CASE3:
**
*NODE, NSET=BLOCK3A
1001, 0.0,2.0,0.0
1002, 0.0,3.0,0.0
1003,-1.0,3.0,0.0
1004,-1.0,2.0,0.0
1011, 0.0,2.0,1.0
1012, 0.0,3.0,1.0
1013,-1.0,3.0,1.0
1014,-1.0,2.0,1.0
**
*ELEMENT, TYPE=C3D8R, ELSET=BLOCK3A
 1001, 1001,1002,1003,1004,1011,1012,1013,1014
**
*NODE, NSET=BLOCK3B
1101, 1.0,2.0,0.0
1102, 1.0,3.0,0.0
1103, 0.0,3.0,0.0
1104, 0.0,2.0,0.0
1111, 1.0,2.0,1.0
1112, 1.0,3.0,1.0
1113, 0.0,3.0,1.0
1114, 0.0,2.0,1.0
*ELEMENT, TYPE=C3D8R, ELSET=BLOCK3B
1101, 1101, 1102, 1103, 1104, 1111, 1112, 1113, 1114
**
*ELSET,ELSET=BLOCK3
BLOCK3A,BLOCK3B
**
*SOLID SECTION,ELSET=BLOCK3,MAT=MAT1
1., 
*NSET, NSET=BLOCK3
BLOCK3A,BLOCK3B
*SURFACE, NAME=SECONDARY3
BLOCK3A,S3
*SURFACE, NAME=MAIN3
BLOCK3B,S5
**
*SURFACE INTERACTION, NAME=INTPROP3
*COHESIVE BEHAVIOR, REPEATED CONTACTS, Eligibility = Current Contacts 
*DAMAGE INITIATION, CRITERION=MAXS
1.E+10,2.E+12,2.E+12
*DAMAGE EVOLUTION, TYPE=DISPLACEMENT, SOFTENING=LINEAR
.002,
*DAMAGE STABILIZATION
1.E-6
**
*CONTACT PAIR, INTERACTION=INTPROP3
SECONDARY3, MAIN3
**
** CASE2:
**
*NODE, NSET=BLOCK2A
2001, 0.0,4.5,0.0
2002,-0.01,5.5,0.0
2003,-1.0,5.5,0.0
2004,-1.0,4.5,0.0
2011, 0.0,4.5,1.0
2012,-0.01,5.5,1.0
2013,-1.0,5.5,1.0
2014,-1.0,4.5,1.0
**
*ELEMENT, TYPE=C3D8R, ELSET=BLOCK2A
 2001, 2001,2002,2003,2004,2011,2012,2013,2014
**
*NODE, NSET=BLOCK2B
2101, 1.0,4.0,0.0
2102, 1.0,5.0,0.0
2103, 0.0,5.0,0.0
2104, 0.0,4.0,0.0
2111, 1.0,4.0,1.0
2112, 1.0,5.0,1.0
2113, 0.0,5.0,1.0
2114, 0.0,4.0,1.0
2115, 1.0,6.0,0.0
2116, 1.0,6.0,1.0
2117, 0.0,6.0,0.0
2118, 0.0,6.0,1.0
*ELEMENT, TYPE=C3D8R, ELSET=BLOCK2B
2101, 2101, 2102, 2103, 2104, 2111, 2112, 2113, 2114
2102, 2102, 2115, 2117, 2103, 2112, 2116, 2118, 2113
**
*ELSET,ELSET=BLOCK2
BLOCK2A,BLOCK2B
**
*SOLID SECTION,ELSET=BLOCK2,MAT=MAT1
1., 
*NSET, NSET=BLOCK2
BLOCK2A,BLOCK2B
*SURFACE, NAME=SECONDARY2
BLOCK2A,S3
*SURFACE, NAME=MAIN2
BLOCK2B,S5
**
*SURFACE INTERACTION, NAME=INTPROP2
*COHESIVE BEHAVIOR, ELIGIBILITY=ORIGINAL 
5.E+13,5.E+13,5.E+13
**
*CONTACT PAIR, INTERACTION=INTPROP2
SECONDARY2, MAIN2
**
** CASE4:
**
*NODE, NSET=BLOCK4A
3001, 0.0,7.0,0.0
3002, 0.0,8.0,0.0
3003,-1.0,8.0,0.0
3004,-1.0,7.0,0.0
3011, 0.0,7.0,1.0
3012, 0.0,8.0,1.0
3013,-1.0,8.0,1.0
3014,-1.0,7.0,1.0
**
*ELEMENT, TYPE=C3D8R, ELSET=BLOCK4A
 3001, 3001,3002,3003,3004,3011,3012,3013,3014
**
*NODE, NSET=BLOCK4B
3101, 1.0,7.0,0.0
3102, 1.0,8.0,0.0
3103, 0.0,8.0,0.0
3104, 0.0,7.0,0.0
3111, 1.0,7.0,1.0
3112, 1.0,8.0,1.0
3113, 0.0,8.0,1.0
3114, 0.0,7.0,1.0
*ELEMENT, TYPE=C3D8R, ELSET=BLOCK4B
3101, 3101, 3102, 3103, 3104, 3111, 3112, 3113, 3114
**
*ELSET,ELSET=BLOCK4
BLOCK4A,BLOCK4B
**
*SOLID SECTION,ELSET=BLOCK4,MAT=MAT1
1., 
*NSET, NSET=BLOCK4
BLOCK4A,BLOCK4B
*SURFACE, NAME=SECONDARY4
BLOCK4A,S3
*SURFACE, NAME=MAIN4
BLOCK4B,S5
**
*SURFACE INTERACTION, NAME=INTPROP4
*COHESIVE BEHAVIOR, TYPE=COUPLED, Eligibility = Current Contacts 
4.E5, 2.E5, 2.E5, 1.E5, 1.E5, 1.E5
*DAMAGE INITIATION, CRITERION=MAXS
8000., 2000., 2000.
*DAMAGE EVOLUTION, TYPE=DISPLACEMENT, SOFTENING=LINEAR
0.06
**
*CONTACT PAIR, INTERACTION=INTPROP4
SECONDARY4, MAIN4
**
*NSET, NSET=TOP
3,4,13,14,1003,1004,1013,1014
*NSET, NSET=TOP2
2003,2004,2013,2014
*NSET, NSET=TOP4
3003,3004,3013,3014
*NSET, NSET=BOT
101,102,111,112,1101,1102,1111,1112,2101,2102,2111,2112,
2115,2116,3101,3102,3111,3112
*BOUNDARY
BOT,1,3,0.0
BLOCK1,2,3
BLOCK3,2,3
**
*AMPLITUDE, NAME=AMP13, DEFINITION=SMOOTH STEP
0.0,0.0,7.E-4,1.E-4,14.E-4,0.,7.E-3,-3.E-3,
20.E-3,1.E-4,30.E-3,-2.E-3
*AMPLITUDE, NAME=AMP2, DEFINITION=SMOOTH STEP
0.0,0.0,7.E-4,1.E-4,14.E-4,0.,7.E-3,-3.E-3,
20.E-3,1.E-4,30.E-3,-1.E-3
*AMPLITUDE, NAME=AMP4, DEFINITION=SMOOTH STEP
0., 0., 30.E-3, 1.E-3
**
*STEP, inc=5000
*STATIC
1.E-4,30.E-3,,1.E-4
*BOUNDARY, AMPLITUDE=AMP13
TOP,1,1,50.0
*BOUNDARY, AMPLITUDE=AMP2
TOP2,1,1,50.0
*BOUNDARY, AMPLITUDE=AMP4
TOP4,1,3,50.0
**Output requests for post-processing
*OUTPUT,FIELD
*NODE OUTPUT
U
***CONTACT OUTPUT
**CFORCE, CSTRESS
*OUTPUT, HISTORY
*CONTACT OUTPUT, SECONDARY=SECONDARY1,MAIN=MAIN1
CFT1,CSDMG,CSMAXSCRT
*CONTACT OUTPUT, SECONDARY=SECONDARY2,MAIN=MAIN2
CFT1
*CONTACT OUTPUT, SECONDARY=SECONDARY3,MAIN=MAIN3
CFT1,CSDMG,CSMAXSCRT
*CONTACT OUTPUT, SECONDARY=SECONDARY4,MAIN=MAIN4
CFT,CSDMG,CSMAXSCRT
**
**
*END STEP