*heading
 Alfano Crisfield non-symmetric multi-delamination analysis
 3-dimensional model in ABAQUS/Explicit 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
** ultI   =  66.0E6
** ultII  =  160.0E6
** cohesive layer modulus (N/m**2):
** Emod = 8.5E9
 Emod = 6.42E13
** Density
** density1 = 1.e3
** density2 = 0.75e3
 density1 = 2.0388e6
 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
*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
** 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,
*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
*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
*contact clearance,name=clearbot,search nset=bondbot
*contact clearance,name=cleartop,search nset=bondtop
*surface interaction,name=lower
*contact damping, definition=critical damping fraction
.1,
*cohesive behavior,eligibility=original contacts
<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>
*surface interaction,name=upper
*contact damping, definition=critical damping fraction
.1,
*cohesive behavior,eligibility=original contacts
<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>
*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, 1.0, 1.0
*filter,type=butterworth, name=one
20,
*step, nlgeom=yes
*dynamic, explicit, fixed time incrementation, scale factor=1.
 ,1.0
*variable mass scaling,type=below min,dt=7.e-7,elset=all_solid,freq=100
*contact
*contact inclusions
bot_top,mid_bot
top_bot,mid_top
*contact property assignment
bot_top,mid_bot,lower
top_bot,mid_top,upper
*contact clearance assignment
bot_top,mid_bot,clearbot
top_bot,mid_top,cleartop
**bulk viscosity
**0.09,1.8
*boundary,amplitude=load
 dispbc-bot, 2,2, -20.0e-3
 dispbc-top, 2,2,  20.0e-3
*nset,nset=qa_test_nset
dispbc,
*output,field,NUMBER INTERVAL=60,TIMEMARKS=YES
*node output
u,
*element output
s,
*node output,nset=qa_test_nset
 u, rf, 
*contact output
cforce,csdmg,csquadscrt
***element output,elset=qa_test_elset
** dmicrt,sdeg
************************
*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