*heading
 Alfano Crisfield non-symmetric multi-delamination analysis
 2-dimensional model in ABAQUS/Standard using COH2D4 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
** ultI   =  66.0E6
** ultII  =  160.0E6
** cohesive layer modulus (N/m**2):
** Emod = 8.5E9
 Emod = 8.5E8
** Emod = 2.18625E11
** Fracture toughness (N/m):
 GIc   = 0.33e3
 GIIc  = 0.80e3
** B-K parameter:
 eta = 2.284
** thickness
 thick = 0.01325e-3
** width in the plane strain direction
 width = 20.0e-3
*parameter
 thickLayer = 0.1325e-3
 t1 = 10.0*thickLayer
 t2 =  2.0*thickLayer
 t3 = 12.0*thickLayer
 length = 200.0e-3
 z1 = t1
 z2 = z1 + thick
 z3 = z2 + t2
 z4 = z3 + thick
 z5 = z4 + t3
*node
 1001,       0.,   0.
 1201, <length>,   0.
 2001,       0., <z1>
 2201, <length>, <z1>
 3001,       0., <z2>
 3201, <length>, <z2>
 4001,       0., <z3>
 4201, <length>, <z3>
 5001,       0., <z4>
 5201, <length>, <z4>
 6001,       0., <z5>
 6201, <length>, <z5>
*ngen, nset=nfront
 1001, 1201
 2001, 2201
 3001, 3201
 4001, 4201
 5001, 5201
 6001, 6201
*nset, nset=nall
 nfront,
*nset, nset=right
 1201,
 2201,
 3201,
 4201,
 5201,
 6201,
*****************************************
*ELEMENT,TYPE=COH2D4, elset=cohesion
 2001, 2001,2002, 3002,3001
 2081, 2081,2082, 3082,3081
 4041, 4041,4042, 5042,5041
*elgen, elset=cohesion
 2001,  60, 1, 1
 2081, 120, 1, 1
 4041, 160, 1, 1
*COHESIVE SECTION,elset=cohesion,material=adhesive,THICKNESS=GEOMETRY,response=traction separation
,<width>
*MATERIAL,NAME=adhesive
*ELASTIC,TYPE=TRACTION
<Emod>,<Emod>,<Emod>
*DAMAGE INITIATION,CRITERION=QUADS
<ultI>, <ultII>
*DAMAGE EVOLUTION,TYPE=ENERGY,MIXED MODE BEHAVIOR=POWER LAW,POWER=1.0
<GIc>, <GIIc>
*****************************************
*element, type=CPE4I
 1001, 1001,1002, 2002,2001
 3001, 3001,3002, 4002,4001
 5001, 5001,5002, 6002,6001
*elgen,elset=bottom
 1001, 200, 1, 1
*elgen,elset=middle
 3001, 200, 1, 1
*elgen,elset=top
 5001, 200, 1, 1
*elset,elset=all_solid
 bottom, middle, top
*solid section, elset=all_solid, mat=mat1, orient=global
 20.e-3,
*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,
*surface,name=bot_top
bottom,s3
*surface,name=mid_bot
middle,s1
*surface,name=mid_top
middle,s3
*surface,name=top_bot
top,s1
*contact pair,interaction=lower
bot_top,mid_bot
*contact pair,interaction=upper
top_bot,mid_top
*surface interaction,name=lower
20.e-3,
*surface interaction,name=upper
20.e-3,
*boundary
**  right, 1,2
 dispbc-bot, 1,1
 dispbc-top, 1,1
*nset, nset=dispbc-bot
 1001,
*nset, nset=dispbc-top
 6001,
*nset, nset=dispbc
 dispbc-bot, dispbc-top
*step, nlgeom, inc=5000
*static, stabilize, factor=7.e5
 0.01, 1.0, 1.e-8
*boundary
 dispbc-bot, 2,2, -20.0e-3
 dispbc-top, 2,2,  20.0e-3
*elset,elset=qa_test_elset
cohesion,
*nset,nset=qa_test_nset
dispbc-top,
*output,field, freq=9999
*node output
u,
*element output
s,
*output,field,freq=9999
*node output,nset=qa_test_nset
u,
**element output,elset=qa_test_elset
** dmicrt,sdeg
*output,history, freq=1
*node output, nset=qa_test_nset
 u2,rf2
*monitor, dof=2, node=6001
*controls, parameter=field
 ,1.0
*controls, parameter=line search
6,5.0,0.2,0.25,0.15
*controls, analysis=discontinuous
*controls, parameter=time incrementation
 , , , , , , , 10,
*end step