*Heading
 Alfano Crisfield non-symmetric multi-delamination analysis
 3-dimensional model in ABAQUS/Explicit using VCCT
**
*Parameter
** User input data defining the interface properties
 time = 0.1
 maxdisp = 20.0e-3
 negmaxdisp = -maxdisp 
** Cohesive layer modulus (N/m**2):
 Emod = 6.42E13
** Density
 density1 = 2.0388e4
 density2 = 1.5291e6
** 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
 width = 1.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>
**
 11001,       0.,   0., <width>
 11201, <length>,   0., <width>
 12001,       0., <z1>, <width>
 12201, <length>, <z1>, <width>
 13001,       0., <z2>, <width>
 13201, <length>, <z2>, <width>
 14001,       0., <z3>, <width>
 14201, <length>, <z3>, <width>
 15001,       0., <z4>, <width>
 15201, <length>, <z4>, <width>
 16001,       0., <z5>, <width>
 16201, <length>, <z5>, <width>
*Ngen, nset=NFRONT
 1001, 1201
 2001, 2201
 3001, 3201
 4001, 4201
 5001, 5201
 6001, 6201
*Ngen, nset=NBACK
 11001, 11201
 12001, 12201
 13001, 13201
 14001, 14201
 15001, 15201
 16001, 16201
*Nset, nset=NALL
 NFRONT, NBACK
*Nset, nset=RIGHT
 1201,
 2201,
 3201,
 4201,
 5201,
 6201,
 11201,
 12201,
 13201,
 14201,
 15201,
 16201, 
**
*Element, type=C3D8R
 1001, 2001,2002,12002,12001, 1001,1002,11002,11001,
 3001, 4001,4002,14002,14001, 3001,3002,13002,13001, 
 5001, 6001,6002,16002,16001, 5001,5002,15002,15001, 
*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, controls=EHC
*Section controls, name=EHC, hourglass=ENHANCED
*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>
*Elset, elset=MIDDLEBONDTOP
 3041, 3042, 3043, 3044, 3045, 3046, 3047, 3048, 
 3049, 3050, 3051, 3052, 3053, 3054, 3055, 3056, 
 3057, 3058, 3059, 3060, 3061, 3062, 3063, 3064, 
 3065, 3066, 3067, 3068, 3069, 3070, 3071, 3072, 
 3073, 3074, 3075, 3076, 3077, 3078, 3079, 3080, 
 3081, 3082, 3083, 3084, 3085, 3086, 3087, 3088, 
 3089, 3090, 3091, 3092, 3093, 3094, 3095, 3096, 
 3097, 3098, 3099, 3100, 3101, 3102, 3103, 3104, 
 3105, 3106, 3107, 3108, 3109, 3110, 3111, 3112, 
 3113, 3114, 3115, 3116, 3117, 3118, 3119, 3120, 
 3121, 3122, 3123, 3124, 3125, 3126, 3127, 3128, 
 3129, 3130, 3131, 3132, 3133, 3134, 3135, 3136, 
 3137, 3138, 3139, 3140, 3141, 3142, 3143, 3144, 
 3145, 3146, 3147, 3148, 3149, 3150, 3151, 3152, 
 3153, 3154, 3155, 3156, 3157, 3158, 3159, 3160, 
 3161, 3162, 3163, 3164, 3165, 3166, 3167, 3168, 
 3169, 3170, 3171, 3172, 3173, 3174, 3175, 3176, 
 3177, 3178, 3179, 3180, 3181, 3182, 3183, 3184, 
 3185, 3186, 3187, 3188, 3189, 3190, 3191, 3192, 
 3193, 3194, 3195, 3196, 3197, 3198, 3199, 3200,
 3040
*Surface, name=BOT_TOP
 BOTTOM, S1
*Surface, name=MID_BOT
 MIDDLE, S2
*Surface, name=MID_TOP
 MIDDLEBONDTOP, S1
*Surface, name=TOP_BOT
 TOP, S2
*Nset, nset=BONDBOT, gen
  2001,  2061
  2081,  2201
 12001, 12061
 12081, 12201
*Nset, nset=BONDTOP, gen
  5041,  5201
 15041, 15201
**
******************************************************
** Initial crack configuration definition
*Contact clearance, name=CLEARBOT, search nset=BONDBOT
*Contact clearance, name=CLEARTOP, search nset=BONDTOP
*Surface interaction, name=LOWER
*Cohesive behavior, eligibility=ORIGINAL CONTACTS
 <Emod>, <Emod>, <Emod>
*Fracture criterion, type=VCCT, mixed mode behavior=BK
 <GIc>, <GIIc>, <GIIc>, <eta>
*Surface interaction, name=UPPER
*Cohesive behavior, eligibility=ORIGINAL CONTACTS
 <Emod>, <Emod>, <Emod>
*Fracture criterion, type=VCCT, mixed mode behavior=BK
 <GIc>, <GIIc>, <GIIc>, <eta>
******************************************************
*Boundary
 RIGHT, 2,2
 DISPBC-BOT, 1,1
 DISPBC-TOP, 1,1
 NFRONT,3,3,
 NBACK,3,3,
*Nset, nset=DISPBC-BOT
 1001, 11001
*nset, nset=DISPBC-TOP
 6001, 16001
*Nset, nset=DISPBC
 DISPBC-BOT, DISPBC-TOP
*Amplitude, definition=SMOOTH STEP, name=LOAD
 0.0, 0.0, <time>, 1.0
*Filter, type=BUTTERWORTH, name=ONE
 200,
**
**
*Step, nlgeom=YES
*Dynamic, Explicit
 , <time>
**
*Contact
*Contact inclusions
 BOT_TOP, MID_BOT
 TOP_BOT, MID_TOP
**
******************************************************
** Initial crack configuration activation
*Contact clearance assignment
 BOT_TOP, MID_BOT, CLEARBOT
 TOP_BOT, MID_TOP, CLEARTOP
*Contact property assignment
 BOT_TOP, MID_BOT, LOWER
 TOP_BOT, MID_TOP, UPPER
******************************************************
*Boundary, amplitude=LOAD
 DISPBC-BOT, 2,2, <negmaxdisp>
 DISPBC-TOP, 2,2, <maxdisp>
*Nset, nset=QA_TEST_NSET
 DISPBC,
**
*Output, field, number interval=200
*Node output
 U,
*Element output
 S,
*Node output, nset=QA_TEST_NSET
 U, RF,
*Contact output
 DBT, DBSF, DBS, OPENBC, CRSTS, ENRRT, EFENRRTR, BDSTAT,
 CFORCE
**
*Output, history, time interval=0.00001, filter=ONE
*Node output, nset=QA_TEST_NSET
 RF2
**
*Output, history, time interval=0.00001
*Node output, nset=DISPBC
 U, RF
*Energy output, variable=PRESELECT
*Monitor, dof=2, node=6001
*End step