*heading
 CPE4I: Alfano Crisfield non-symmetric multi-delamination
 analysis
*preprint, model=yes, echo=yes, history=yes
*parameter
** Fracture toughness (N/m):
 GIc   = 0.33e3
 GIIc  = 0.80e3
 GIIIc = 0.0
** B-K parameter:
 modeMixLaw=1
 eta=2.284
** mixed mode parameter
 am = 0.0
 an=0.0
 ao=0.0
** Damage and tolerance parameters
damv=0.0
tol=0.2 
** 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=crack2_top,gen
3061,3080,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=bonded_upper,generate
5041,5200,1
*nset,nset=bonded_lower,generate
2002,2061,1
2081,2201,1
*initial conditions,type=contact
top_bot, mid_top, bonded_upper
*initial conditions,type=contact
bot_top, mid_bot, bonded_lower
*surface,name=bot_top
bottom,s3
*surface,name=mid_bot
middle,s1
*surface,name=mid_top
middle,s3
*surface,name=top_bot
top,s1
*surface,name=crack2_top
crack2_top,s1
*contact pair,interaction=lower,adjust=bonded_lower
bot_top,mid_bot
*contact pair,interaction=upper,adjust=bonded_upper
top_bot,mid_top
*surface interaction,name=lower
20.e-3,
*surface interaction,name=upper
20.e-3,
*boundary
 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=1000
*static, stabilize, factor=7.e5
 0.01, 1.0, 1.e-8
*debond, slave=top_bot, master=mid_top
*FRACTURE CRITERION,TYPE=VCCT,MIXED MODE BEHAVIOR=BK
<GIc>,<GIIc>,<GIIIc>,<eta>
*debond, slave=bot_top, master=mid_bot
*FRACTURE CRITERION,TYPE=VCCT,MIXED MODE BEHAVIOR=BK
<GIc>,<GIIc>,<GIIIc>,<eta>
*boundary
 dispbc-bot, 2,2, -20.0e-3
 dispbc-top, 2,2,  20.0e-3
*output,field, freq=10
*element output
s
**s,e
*node output
u
*contact output,SLAVE=bot_top,MASTER=mid_bot
dbt,dbsf,dbs,openbc,crsts,enrrt,efenrrtr,bdstat
*contact output,SLAVE=top_bot,MASTER=mid_top
dbt,dbsf,dbs,openbc,crsts,enrrt,efenrrtr,bdstat
*output,history
*node output,nset=dispbc-top
u2, rf2
*contact output,SLAVE=bot_top,MASTER=mid_bot
dbt
*contact output,SLAVE=top_bot,MASTER=mid_top
dbt
**energy output
**allse
**allvd
**allsd
*monitor, dof=2, node=6001
*controls, parameter=field
 , 1.0
***controls, parameter=line search
** 6, 5.0, 0.20
*controls, analysis=discontinuous
*controls, parameter=time incrementation
 , , , , , , , 10,
*end step