*HEADING
   3D Double cantilever beam test model 
**
** Units: in, lbf, psi 
**
*Parameter
**
** Higher density is used here to reduce analysis time. To
** reproduce results as shown in the Example Problem Manual,
** lower the density by an order of 2.
**
** Density (lbf/in**3)
 density1 = 1.91E-1
** density1 = 1.91E-3
** User input data defining the interface properties
** Fracture toughness (lbf/in):
 GIc   = 1.2
 GIIc  = 6.0
 GIIIc = 6.0
** B-K parameter:
 eta = 1.75
** Analysis time 
 time_step = 0.3
**
*Node
  1,   0.0, 0.0
 11,   0.0, 1.0
 51,   0.0, 0.0
 61,   0.0, 1.0
** 
**
 9001,   4.5, 0.0
 9011,   4.5, 1.0
 9051,   4.5, 0.0
 9061,   4.5, 1.0
**
18001,   9.0, 0.0
18011,   9.0, 1.0
18051,   9.0, 0.0
18061,   9.0, 1.0
**
*Ngen, nset=NS1
    1 , 9001, 100
   51 , 9051, 100
 9001 ,18001, 100
 9051 ,18051, 100
*Ngen, nset=NS2
   11 , 9011, 100
   61 , 9061, 100
 9011 ,18011, 100
 9061 ,18061, 100
**
*Nfill
NS1,NS2,10,1
**
**
*Element, type=S4
20001, 1, 101, 102, 2
20051,51, 151, 152,52
**
*Elgen, elset=LOWER_ELEMENTS
20001, 180,100,100, 10, 1,1 
*Elgen, elset=UPPER_ELEMENTS
20051, 180,100,100, 10, 1,1
**
*Shell Section, elset=UPPER_ELEMENTS, material=MAT, orientation=OR1
0.2,
** 
*Shell Section, elset=LOWER_ELEMENTS, material=MAT, orientation=OR1
0.2,
** 
*Orientation, name=OR1
1.0, 0.0, 0.0, 0.0, 1.0, 0.0
3, 0.0 
**
*Material, name=MAT
*Elastic, type=ENGINEERING CONSTANTS
8.0E+6, 8.0E+6, 8.0E+6, 0.0, 0.0, 0.0, 4.0E+6, 4.0E+6,
4.0E+6, 
*Density
<density1>
**
*Elset, elset=SLAVE, generate
20101,37901,100
20102,37902,100
20103,37903,100
20104,37904,100
20105,37905,100
20106,37906,100
20107,37907,100
20108,37908,100
20109,37909,100
20110,37910,100
*Surface, type=element, name=SLAVE
SLAVE, SPOS
*Surface, type=element,name=MASTER
UPPER_ELEMENTS, SNEG
*Elset, elset=ETIED, generate
29101,37901,100
29102,37902,100
29103,37903,100
29104,37904,100
29105,37905,100
29106,37906,100
29107,37907,100
29108,37908,100
29109,37909,100
29110,37910,100
**
*Nset, nset=NTIED, elset=ETIED
**
********************************************************
*Contact Clearance, name=CLEARMID, search nset=NTIED
*Surface Interaction, name=MIDDLE
*Cohesive Behavior, eligibility=ORIGINAL CONTACTS
*Fracture Criterion, type=VCCT, mixed mode behavior=REEDER
 <GIc>, <GIIc>, <GIIIc>, <eta>
********************************************************
*Nset, nset=FIX, generate
18001, 18011
18051, 18061
**
*Boundary
FIX,1,6
**
*MPC
BEAM, 1, 6
BEAM, 2, 6
BEAM, 3, 6
BEAM, 4, 6
BEAM, 5, 6
BEAM, 7, 6
BEAM, 8, 6
BEAM, 9, 6
BEAM,10, 6
BEAM,11, 6
**
BEAM, 51, 56
BEAM, 52, 56
BEAM, 53, 56
BEAM, 54, 56
BEAM, 55, 56
BEAM, 57, 56
BEAM, 58, 56
BEAM, 59, 56
BEAM, 60, 56
BEAM, 61, 56
**
**
*Nset, nset=LP
 6, 56
**
*Amplitude, definition=SMOOTH STEP, name=LOAD
 0.0, 0.0, <time_step>, 1.0
**
*Step
*Dynamic, Explicit
 , <time_step>
*Bulkviscosity
0.06, 1.2
**
*Contact
*Contact Inclusions
SLAVE, MASTER
**
******************************************************
*Contact Clearance Assignment
SLAVE, MASTER, CLEARMID
*Contact Property Assignment
SLAVE, MASTER, MIDDLE
*Surface Property Assignment, property=THICKNESS
SLAVE, 0.0
MASTER, 0.0
******************************************************
*Boundary, amplitude=LOAD
**
 6,3,3, -0.15
56,3,3,  0.15
**
*Output, field, number interval=100
*Node Output
 U,
*Element Output
 S,
*Contact Output
 DBT, DBSF, DBS, OPENBC, CRSTS, ENRRT, EFENRRTR, BDSTAT,
** CFORCE, CSTRESS
**
*Output, history, time interval=0.0001
***Node Output, nset=LP
** U3, RF3
*Energy Output, variable=PRESELECT
*End Step