*HEADING
   DCB test model using surface based cohesive behavior in 2D   
*PREPRINT,MODEL=YES,HISTORY=YES
*parameter
** Ultimate strenth in tensile and mode II:
 ultI  =  8.0E2
 ultII =  8.0E2
** cohesive layer modulus:
 Emod = 8.51264378633914E10
** Fracture toughness:
 GIc   = 0.8
 GIIc  = 0.8
** B-K parameter:
 eta=1.75
** width in the plane strain direction
 width =1.0
**
*NODE
    1,  0.0, 0.0, 0.0
    5,  0.0, 0.2, 0.0
    6,  0.0, 0.202, 0.0
   10,  0.0, 0.4, 0.0 
**  
**
 4001,  4.03, 0.0, 0.0
 4005,  4.03, 0.2, 0.0
 4006,  4.03, 0.202, 0.0
 4010,  4.03, 0.4, 0.0  
**
 8001,  8.03, 0.0, 0.0
 8005,  8.03, 0.2, 0.0
 8006,  8.03, 0.2, 0.0
 8010,  8.03, 0.4, 0.0  
**
 9001,  9.0, 0.0, 0.0
 9005,  9.0, 0.2, 0.0
 9006,  9.0, 0.2, 0.0
 9010,  9.0, 0.4, 0.0  
**
*NGEN,NSET=L1
    1,   5,1
    6,  10,1
*NGEN,NSET=L2
 4001,4005,1
 4006,4010,1
*NGEN,NSET=L5
 8001,8005,1
 8006,8010,1
*NGEN,NSET=L6
 9001,9005,1
 9006,9010,1
**
*NFILL,NSET=ALLN
L1,L2,40,100
*NFILL,NSET=ALLN
L2,L5,40,100
*NFILL,NSET=ALLN
L5,L6,10,100
**
*NSET,NSET=BC,GENERATE
9001, 9004
9006, 9010
*nset,nset=nboun
7
*BOUNDARY
BC,1,3
**
**
*ELEMENT,TYPE=CPE4
1, 1, 101,102, 2
6, 6, 106,107, 7
**
*ELGEN,ELSET=EBASE
1, 90,100,100, 4, 1,1
6, 90,100,100, 4, 1,1
*ELSET,ELSET=BOT,gen
 4, 8904, 100
*ELSET,ELSET=TOP,gen
 6, 8906, 100
** 
**
*SOLID SECTION,ELSET=EBASE,MATERIAL=BMAT
**
*MATERIAL,NAME=BMAT
*ELASTIC, TYPE=ISO
      8.0E+6,         0.3
**
**
*NSET,NSET=BNODES,GEN
 4205, 9005, 100
*SURFACE,NAME=BOT
 BOT, S3
*SURFACE,NAME=TOP
 TOP, S1 
*NSET,NSET=AD,GEN
 5,9005,100
*CONTACT PAIR, INTERACTION=FRACT,ADJUST=BNODES
**CONTACT PAIR, INTERACTION=FRACT
 BOT, TOP
*SURFACE INTERACTION, NAME=FRACT
1.0
*cohesive behavior,eligibility=specified contacts
**<Emod>,<Emod>
*DAMAGE INITIATION,CRITERION=QUADS
<ultI>, <ultII>
*DAMAGE EVOLUTION,TYPE=ENERGY,MIXED MODE BEHAVIOR=BK,POWER=<eta>
<GIc>, <GIIc>
**DAMAGE STABILIZATION
** 1.e-5
**
*INITIAL CONDITIONS,TYPE =CONTACT
 BOT, TOP, BNODES
**
**step
***static
**1,1
***endstep
*STEP,NLGEOM, INC=5000,convert sdi=no
*STATIC
0.005, 1,, 0.015
*CONTROLS,PARAMETERS=TIME INCREMENTATION
 , , , , , , 50
*BOUNDARY
  4, 2,2, -0.16
  7, 2,2, 0.16
**contact controls,SLAVE=BOT,MASTER=TOP,stabilize=1.e-3
**print,contact=yes
*output,field, freq=10
*element output
s,e
*node output
u
*output, history
*node output,nset=nboun
u2
rf2
*END STEP