*HEADING
   3D Single leg bend Delamination model 
**
**
*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   = 2.0
 GIIc  = 2.0
 GIIIc = 30.0
** B-K parameter:
 eta = 1.75
** Analysis time 
 time_step = 0.15
**
*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, offset=-0.5
0.066,
** 
*Shell Section, elset=LOWER_ELEMENTS, material=MAT, orientation=OR1, offset=0.5
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
20001,37901,100
20002,37902,100
20003,37903,100
20004,37904,100
20005,37905,100
20006,37906,100
20007,37907,100
20008,37908,100
20009,37909,100
20010,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=SURFS 
*Cohesive Behavior, eligibility=ORIGINAL CONTACTS
*Fracture Criterion,type=VCCT, mixed mode behavior=REEDER
 <GIc>, <GIIc>, <GIIIc>, <eta>
******************************************************
*MPC
BEAM, 12051, 12056
BEAM, 12052, 12056
BEAM, 12053, 12056
BEAM, 12054, 12056
BEAM, 12055, 12056
BEAM, 12057, 12056
BEAM, 12058, 12056
BEAM, 12059, 12056
BEAM, 12060, 12056
BEAM, 12061, 12056
**
** Boundary constraints
**
*Nset, nset=RL, generate
 6051, 6061, 1
*Nset, nset=RR, generate
18051, 18061, 1
*Nset, nset=NOSLIDE
 6051, 18051
*Nset, nset=NOSLIP
 6051
*Nset, nset=LP
 12056
**
*Boundary
 RL, 3
 RR, 3
 NOSLIDE, 2
 NOSLIP, 1
**
*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, SURFS
*Surface Property Assignment, property=THICKNESS
SLAVE, 0.0
MASTER, 0.0
******************************************************
*Boundary, amplitude=LOAD
LP, 3,3, -0.15
**
** OUTPUT REQUESTS 
*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