*heading CPE4I: Alfano Crisfield non-symmetric multi-delamination low cycle fatigue analysis using direct cyclic approach *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, , 0. 2001, 0., 2201, , 3001, 0., 3201, , 4001, 0., 4201, , 5001, 0., 5201, , 6001, 0., 6201, , *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,small sliding bot_top,mid_bot *contact pair,interaction=upper,adjust=bonded_upper,small sliding 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 *amplitude,name=test 0,0,0.5,1,1,0 *amplitude,name=test2 0,0,0.5,1,1,0.0001 *step, nlgeom, inc=1000 *static,direct 0.01, 1.0, 1.e-8 **contact controls,automatic tolerances *debond, secondary=top_bot, main=mid_top *FRACTURE CRITERION,TYPE=VCCT,MIXED MODE BEHAVIOR=BK ,,, *debond, secondary=bot_top, main=mid_bot *FRACTURE CRITERION,TYPE=VCCT,MIXED MODE BEHAVIOR=BK ,,, *boundary,amp=test2 dispbc-bot, 2,2, -1.0e-3 dispbc-top, 2,2, 1.0e-3 *output,field, freq=50 *element output s *node output u *contact output,SECONDARY=bot_top,MAIN=mid_bot dbt,dbsf,dbs,openbc,crsts,enrrt,efenrrtr,bdstat *contact output,SECONDARY=top_bot,MAIN=mid_top dbt,dbsf,dbs,openbc,crsts,enrrt,efenrrtr,bdstat *output,history,freq=50 *node output,nset=dispbc-top u2, rf2 *controls, parameter=field , 1.0 *controls, analysis=discontinuous *controls, parameter=time incrementation , , , , , , , 10, *end step *step, inc=10000 *direct cyclic,fatigue 0.01,1,,,25,25,,20 3,6,91,1.1 **contact controls,automatic tolerances *controls,type=direct cyclic ,100,1E-3,100,1E-3 *debond, secondary=top_bot, main=mid_top *FRACTURE CRITERION,TYPE=FATIGUE 0.5,-1.44,4.8768E-6,1.15,,,0.33e3,0.80e3, 0,2.284 *debond, secondary=bot_top, main=mid_bot *FRACTURE CRITERION,TYPE=FATIGUE 0.5,-1.44,4.8768E-6,1.15,,,0.33e3,0.80e3, 0,2.284 *boundary,amp=test dispbc-bot, 2,2, -1.0e-3 dispbc-top, 2,2, 1.0e-3 *output,field, freq=50 *element output s *node output u *contact output,SECONDARY=bot_top,MAIN=mid_bot dbt,dbsf,dbs,openbc,crsts,enrrt,efenrrtr,bdstat *contact output,SECONDARY=top_bot,MAIN=mid_top dbt,dbsf,dbs,openbc,crsts,enrrt,efenrrtr,bdstat *output,freq=50,history *node output,nset=dispbc-top u2, rf2 *controls, parameter=field , 1.0 *controls, analysis=discontinuous *controls, parameter=time incrementation , , , , , , , 10, *end step