*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 = 6.42E13
** 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=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,
*nset,nset=bondbot,generate
2001,2061
2081,2201
*nset,nset=bondtop,generate
5041,5201
*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,adjust=bondbot
bot_top,mid_bot
*contact pair,interaction=upper,adjust=bondtop
top_bot,mid_top
*surface interaction,name=lower
20.e-3,
*cohesive behavior,eligibility=original contacts
<Emod>,<Emod>,<Emod>
*DAMAGE INITIATION,CRITERION=QUADS
<ultI>, <ultII>
*DAMAGE EVOLUTION,TYPE=ENERGY,MIXED MODE BEHAVIOR=POWER LAW,POWER=1.0
<GIc>, <GIIc>
*DAMAGE STABILIZATION
 1.e-7
*surface interaction,name=upper
20.e-3,
*cohesive behavior,eligibility=original contacts
<Emod>,<Emod>,<Emod>
*DAMAGE INITIATION,CRITERION=QUADS
<ultI>, <ultII>
*DAMAGE EVOLUTION,TYPE=ENERGY,MIXED MODE BEHAVIOR=POWER LAW,POWER=1.0
<GIc>, <GIIc>
*DAMAGE STABILIZATION
 1.e-7
*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
*nset,nset=qa_test_nset
dispbc-top,
*output,field, freq=9999
*node output
u,
*element output
s,
*output,field,freq=9999
*node output
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