*heading
 3-dimensional model in ABAQUS/Standard using COH3D8 elements
*preprint, model=yes, echo=yes, history=yes
*parameter
** User input data defining the interface properties
** Ultimate strenth in tensile and mode II:
 ultI   =  3.30E6
 ultII  =  7.0E6
** cohesive layer modulus (N/m**2):
 Emod = 8.5E8
** Density
 density1 = 2.0388e6
 density2 = 1.5291e6
** Fracture toughness (N/m):
 GIc   = 0.33e3
 GIIc  = 0.80e3
** B-K parameter:
  eta = 2.284
** cohesive layer thickness
 thick = 0.01325e-3
** width in the plane strain direction
 w1=200
 width = w1*0.001
**delamination details
 dlx=90
 dly=50
 sdlx=5
 sdly=5
*parameter
** dl1&dl2 incdicate the bottom and upper  layeres. delamination is present bitween these layers
 dl1=2
 dl2=6
 thickLayer = 0.1325e-3
 t1 = dl1*thickLayer
 t2 = dl2*thickLayer
 l1=250
 length = l1*0.001
 z1 = t1
 z2 = z1 + thick
 z3 = z2 + t2
 *node
 100001,       0.,    0.
 100251, <length>,    0.
 150201,       0.,   <width>,   0.   
 150451, <length>,   <width>,   0.
 200001,       0.,        0.,  <z1>
 200251, <length>,        0.,  <z1>
 250201,       0.,   <width>,  <z1>
 250451, <length>,   <width>,  <z1>
 300001,       0.,        0.,  <z2>
 300251, <length>,        0.,  <z2>
 350201,       0.,   <width>,  <z2>
 350451, <length>,   <width>,  <z2>
 400001,       0.,        0.,  <z3>
 400251, <length>,        0.,  <z3>
 450201,       0.,   <width>,  <z3>
 450451, <length>,   <width>,  <z3>
*****************************************
 *ngen, nset=b1_front
 100001,100251
 *ngen, nset=b1_back
 150201,150451
 *nfill, nset=b1,bias=1
 b1_front,b1_back, 200,251
******************************************
 *ngen, nset=b2_front
 200001,200251
 *ngen, nset=b2_back
 250201,250451
 *nfill, nset=b2,bias=1
 b2_front,b2_back, 200,251
********************************************* 
 *ngen, nset=b3_front
 300001,300251
 *ngen, nset=b3_back
 350201,350451
 *nfill, nset=b3,bias=1
 b3_front,b3_back, 200,251
********************************************
 *ngen, nset=b4_front
 400001,400251
 *ngen, nset=b4_back
 450201,450451
 *nfill, nset=b4,bias=1
 b4_front,b4_back, 200,251
******************************************
 *nset,nset=nall
b1,b2,b3,b4
*ngen,nset=right_1
100001,150201,251
*ngen,nset=right_2
200001,250201,251
*ngen,nset=right_3
300001,350201,251
*ngen,nset=right_4
400001,450201,251
*nset,nset=right
right_1,right_2,right_3,right_4
********************************************
**manipulations for delamination location
*parameter
** zone 2 parameters
re2=dlx+1+200000
rn21=200001+dlx
rn22=200001+dlx+1
rn23=200001+dlx+1+251
rn24=200001+dlx+251
rn25=300001+dlx
rn26=300001+dlx+1
rn27=300001+dlx+251+1
rn28=300001+dlx+251
****************************************
**forzone 3 parameters
re3= 200001+dlx+sdlx
rn31=200001+dlx+sdlx
rn32=200001+dlx+sdlx+1
rn33=200001+dlx+sdlx+251+1
rn34=200001+dlx+sdlx+251
rn35=300001+dlx+sdlx
rn36=300001+dlx+sdlx+1
rn37=300001+dlx+sdlx+1+251
rn38=300001+dlx+sdlx+251
***********************************************
**for zone 4 parameters
re4  =200001+250*dly+250*sdly+90
rn41=200001+251*dly+251*sdly+90
rn42=200001+251*dly+251*sdly+90+1
rn43=200001+251*dly+251*sdly+1+90+251
rn44=200001+251*dly+251*sdly+90+251
rn45=300001+251*dly+251*sdly+90
rn46=300001+251*dly+251*sdly+90+1
rn47=300001+251*dly+251*sdly+1+90+251
rn48=300001+251*dly+251*sdly+90+251
************************************************
dv1=250-dlx-sdlx
dv2=dly+sdly
dv3=250-dlx
dv4=200-dly-sdly
****************************************************
*ELEMENT,TYPE=COH3D8, elset=cohesion
200001,  200001, 200002, 200253, 200252,  300001, 300002, 300253, 300252
<re2>,  <rn21>,<rn22>,<rn23>,<rn24>, <rn25>,<rn26>,<rn27>,<rn28>
<re3>,  <rn31>,<rn32>,<rn33>,<rn34>, <rn35>,<rn36>,<rn37>,<rn38>
<re4>,  <rn41>,<rn42>,<rn43>,<rn44>, <rn45>,<rn46>,<rn47>,<rn48>
**
*elgen, elset=cohesion
**zone1
200001,<dlx>,   1,   1,   200,  251, 250
**zone2
 <re2>,<sdlx>,   1,   1,     <dly>,  251, 250
**zone3
<re3>, <dv1>,   1,   1,  <dv2>,  251, 250
**zone4
<re4>, <dv3>,   1,   1,   <dv4>,  251, 250
*******************************************************
*COHESIVE SECTION,elset=cohesion,material=adhesive,THICKNESS=GEOMETRY,response=traction separation
*MATERIAL,NAME=adhesive
*Density
2000.,
*ELASTIC,TYPE=TRACTION
<Emod>,<Emod>,<Emod>
*DAMAGE INITIATION,CRITERION=QUADS
<ultI>, <ultII>, <ultII>
*DAMAGE EVOLUTION,TYPE=ENERGY,MIXED MODE BEHAVIOR=POWER LAW,POWER=1.0
<GIc>, <GIIc>, <GIIc>
***********************************************************
*element, type=C3D8I
100001,  100001,100002,100253,100252,200001,200002,200253,200252,
300001,  300001,300002,300253,300252,400001,400002,400253,400252, 
*elgen,elset=bottom
100001,250,1,1,200,251,250
*elgen,elset=top
300001,250,1,1,200,251,250
*elset,elset=all_solid
 bottom,top
*solid section, elset=all_solid, mat=mat1, orient=global
*orientation, name=global
 1,0,0, 0,1,0
*material, name=mat1
*elastic, type=engineering
**    E1,    E2,    E3, nu12, nu13, nu23,   G12,   G13,
**   G23
 115.0E9, 8.5E9, 8.5E9, 0.29, 0.29,  0.3, 4.5E9, 3.3E9
 4.5E9,
*density
<density1>
*surface,name=bot_top
bottom,s2
*surface,name=top_bot
top,s1
*contact pair,interaction=top&bottom_pair
bot_top,top_bot
*surface interaction,name=top&bottom_pair
*Boundary
 right, 1, 1
 right, 2, 2
 right, 3, 3
*STEP
*FREQUENCY, EIGENSOLVER=Lanczos, NORM=Mass 
,,300.0
**
*EL PRINT,FREQUENCY=0
*MODAL FILE
*PRINT, SOLVE=YES
*END STEP