*HEADING
   3D Single leg bend Delamination model 
**
**preprint,contact=yes,history=yes
**Restart, Write
**
*parameter
** User input data defining the interface properties
** Fracture toughness
 GIc   = 2.0
 GIIc  = 2.0
 GIIIc = 30.0
** B-K parameter:
 eta = 1.75
 modeMixLaw=3
** mixed mode parameter
 am = 1.0
 an=1.5
 ao=1.0
** Damage and tolerance parameters
damv=0
**5.0e-5
tol=0.4
width=12.7 
*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
20001, 180,100,100, 10, 1,1 
*ELGEN,ELSET=upper
20051, 180,100,100, 10, 1,1
**
**
*Shell Section, elset=upper, material=MAT_upper,Orientation=Or1,OFFSET=-0.5
0.066,
** 
*Shell Section, elset=lower, material=MAT_lower,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_upper
*ELASTIC,TYPE=ENGINEERING CONSTANTS
8000000.0,8000000.0,8000000.0, 0.0, 0.0, 0.0,4000000.0,4000000.0,
4000000.0, 
**
*MATERIAL,NAME=MAT_lower
*ELASTIC,TYPE=ENGINEERING CONSTANTS
8000000.0,8000000.0,8000000.0, 0.0, 0.0, 0.0,4000000.0,4000000.0,
4000000.0,
*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,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
*initial conditions, type=contact
slave, master,ntied
*Contact Pair, interaction=SURFS, small sliding
slave, master
*Surface Interaction, name=SURFS
1.,
*Friction, slip tolerance=0.005
0.,
*clearance,slave=slave,master=master,value=3.0E-9
*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
** STEP: Step-1
**
*STEP, INC=1000
**,nlgeom
**
*STATIC
**,stabilize=1.0E-8
0.035,0.07
**contact controls,slave=slave, master=master,stabilize=5.E-6
*debond, slave=slave, master=master
*FRACTURE CRITERION,TYPE=VCCT,MIXED MODE BEHAVIOR=REEDER,TOLERANCE=<tol>
<GIc>, <GIIc>, <GIIIc>, <eta>
*CONTROLS,PARAMETERS=TIME INCREMENTATION
 , , , , , , , 30,
*BOUNDARY
LP, 3,3, -0.07
**
**
*print,contact=yes
** OUTPUT REQUESTS
** 
**Output, field, variable=PRESELECT
**Output, history, variable=PRESELECT
**NODE PRINT,NSET=LP
** U3, RF3
**
**
**output, field,freq=2
**node output
**u
**contact output,slave=slave,master=master
**sdv
**cstress,
**cdisp
*Output, history,freq=1
*NODE output,NSET=LP
 U3, RF3
**energy output,var=all
*End Step
*STEP, INC=1000
**,nlgeom
**
*STATIC
**,stabilize=1.0E-8
0.0005, 0.01,1.0E-8,0.0005
*contact controls,slave=slave, master=master,stabilize=5.E-6
**debond, slave=slave, master=master
**FRACTURE CRITERION,TYPE=VCCT,MIXED MODE BEHAVIOR=REEDER,TOLERANCE=<tol>
**<GIc>, <GIIc>, <GIIIc>, <eta>
*CONTROLS,PARAMETERS=TIME INCREMENTATION
 , , , , , , , 30,
*BOUNDARY
LP, 3,3, -0.08
**
**
**print,contact=yes
** OUTPUT REQUESTS
** 
**Output, field, variable=PRESELECT
**Output, history, variable=PRESELECT
**NODE PRINT,NSET=LP
** U3, RF3
**
**
*output, field,freq=10
*node output
u
*element output
s
*contact output,slave=slave,master=master
dbt,dbsf,dbs,openbc,crsts,enrrt,efenrrtr,bdstat
**cstress,
**cdisp
*Output, history,freq=1
*NODE output,NSET=LP
 U3, RF3
**energy output,var=all
*End Step