*HEADING 2D Single leg bend Delamination model tf=0.8 ply ts=24 ply (4% zpins not) use vcct based low cycle fatigue *parameter ** Fracture toughness: GIc = 5.63 GIIc = 9.305 GIIIc = 0.0 ** Spatially varying properties GIc1=1.505 GIc2=1.6686 GIc3=1.8279 GIc4=2.14222 GIc5=2.4436 GIc6=2.732 GIc7=3.011 GIc8=3.2744 GIc9=3.524 GIc10=3.7607 GIc11=3.98443 GIc12=4.1952 GIc13=4.39308 GIc14=4.57799 GIc15=4.74997 GIc16=4.909 GIc17=5.0551 GIc18=5.188 GIc19=5.30848 GIc20=5.4157 GIc21=5.51 GIc22=5.59 GIc23=5.63 ** B-K parameter: modeMixLaw=1 eta=1.75 ** mixed mode parameter am = 0.0 an=0.0 ao=0.0 damv=0 **5.0E-5 tol=0.2 ** width in the plane strain direction width =1.0 ** All Nodes *NODE 100000, 0.00, 0.0000, 0.0 100024, 0.00, 0.1992, 0.0 200024, 0.00, 0.1992, 0.0 200048, 0.00, 0.3984, 0.0 ** 160000, 6.00, 0.0000, 0.0 160024, 6.00, 0.1992, 0.0 260024, 6.00, 0.1992, 0.0 260048, 6.00, 0.3984, 0.0 ** *NGEN,NSET=FL 100000, 100024, 1 *NGEN,NSET=FR 160000, 160024, 1 *NFILL,NSET=NFLANGE FL,FR,200,300 ** *NGEN,NSET=SL 200024, 200048, 1 *NGEN,NSET=SR 260024, 260048, 1 *NFILL,NSET=NSKIN SL,SR,200,300 ** ** FLANGE ELEMENTS *ELEMENT, TYPE=CPE4I 700017,100016, 100316, 100317, 100017, *ELGEN 700017,8,1,1 ** *ELGEN, ELSET=EFLANGE_0 **700017,160,300,300 700018,160,300,300 **700019,160,300,300 700020,160,300,300 700021,160,300,300 **700022,160,300,300 700023,160,300,300 **700024,160,300,300 ** *ELGEN, ELSET=EFLANGE_45 700017,160,300,300 **700018,160,300,300 700019,160,300,300 **700020,160,300,300 **700021,160,300,300 700022,160,300,300 **700023,160,300,300 700024,160,300,300 ** ** SKIN ELEMENTS *ELEMENT, TYPE=CPE4I 800025,200024, 200324, 200325, 200025, *ELGEN 800025, 24,1,1 ** *ELGEN, ELSET=ESKIN_0 **800025,200,300,300 800026,200,300,300 **800027,200,300,300 800028,200,300,300 **800029,200,300,300 800030,200,300,300 **800031,200,300,300 800032,200,300,300 **800033,200,300,300 800034,200,300,300 **800035,200,300,300 800036,200,300,300 800037,200,300,300 **800038,200,300,300 800039,200,300,300 **800040,200,300,300 800041,200,300,300 **800042,200,300,300 800043,200,300,300 **800044,200,300,300 800045,200,300,300 **800046,200,300,300 800047,200,300,300 **800048,200,300,300 ** *ELGEN, ELSET=ESKIN_45 800025,200,300,300 **800026,200,300,300 800027,200,300,300 **800028,200,300,300 800029,200,300,300 **800030,200,300,300 800031,200,300,300 **800032,200,300,300 800033,200,300,300 **800034,200,300,300 800035,200,300,300 **800036,200,300,300 **800037,200,300,300 800038,200,300,300 **800039,200,300,300 800040,200,300,300 **800041,200,300,300 800042,200,300,300 **800043,200,300,300 800044,200,300,300 **800045,200,300,300 800046,200,300,300 **800047,200,300,300 800048,200,300,300 ** *BOUNDARY 200024, 1 200024, 2 260024, 2 ** *ORIENTATION, NAME=OID1 1,0,0, 0,1,0 *SOLID SECTION,ELSET=ESKIN_0, MATERIAL=PW8552_0, ORIENTATION=OID1 *SOLID SECTION,ELSET=ESKIN_45, MATERIAL=PW8552_0, ORIENTATION=OID1 *SOLID SECTION,ELSET=EFLANGE_0, MATERIAL=PW8552_0, ORIENTATION=OID1 *SOLID SECTION,ELSET=EFLANGE_45, MATERIAL=PW8552_0, ORIENTATION=OID1 ** ** TYPE=ENGINEERING CONSTANTS ** E1, E2, E3, Nu12, Nu13, Nu23, G12, G13 ** G23 ** *MATERIAL,NAME=PW8552_0 *ELASTIC 9100000.0,0.33 ** ** Fracture elements (Rcurve) and MPC TIE included *MPC TIE, 100024, 200024 TIE, 100324, 200324 TIE, 100624, 200624 TIE, 100924, 200924 TIE, 101224, 201224 TIE, 101524, 201524 TIE, 101824, 201824 TIE, 102124, 202124 TIE, 102424, 202424 TIE, 102724, 202724 TIE, 103024, 203024 TIE, 103324, 203324 TIE, 103624, 203624 TIE, 103924, 203924 TIE, 104224, 204224 TIE, 104524, 204524 TIE, 104824, 204824 TIE, 105124, 205124 TIE, 105424, 205424 TIE, 105724, 205724 TIE, 106024, 206024 TIE, 106324, 206324 TIE, 106624, 206624 TIE, 106924, 206924 TIE, 107224, 207224 TIE, 107524, 207524 TIE, 107824, 207824 TIE, 108124, 208124 TIE, 108424, 208424 TIE, 108724, 208724 TIE, 109024, 209024 TIE, 109324, 209324 TIE, 109624, 209624 TIE, 109924, 209924 TIE, 110224, 210224 TIE, 110524, 210524 TIE, 110824, 210824 TIE, 111124, 211124 TIE, 111424, 211424 TIE, 111724, 211724 TIE, 112024, 212024 TIE, 112324, 212324 TIE, 112624, 212624 TIE, 112924, 212924 TIE, 113224, 213224 TIE, 113524, 213524 TIE, 113824, 213824 TIE, 114124, 214124 TIE, 114424, 214424 TIE, 114724, 214724 TIE, 115024, 215024 TIE, 115324, 215324 TIE, 115624, 215624 TIE, 115924, 215924 TIE, 116224, 216224 TIE, 116524, 216524 TIE, 116824, 216824 TIE, 117124, 217124 TIE, 117424, 217424 TIE, 117724, 217724 TIE, 118024, 218024 TIE, 118324, 218324 TIE, 118624, 218624 TIE, 118924, 218924 TIE, 119224, 219224 TIE, 119524, 219524 TIE, 119824, 219824 ** *nset, nset=bond, generate 120124, 145324, 300 *elset, elset=top, generate 820125, 859725, 300 *elset, elset=bot, generate 720124, 747724, 300 *surface, type=element, name=topsurf top, s1 *surface, type=element, name=botsurf bot, s3 *contact pair, interaction=crack, adjust=bond,small sliding botsurf, topsurf *surface interaction, name=crack 1.0 *initial conditions, type=contact botsurf, topsurf, bond ** *NSET, NSET=LP 230024 *nset,nset=L145324 145324 *NODAL ENERGY RATE bond,, , 145324,,, 145024,,, 144724,,, 144424,,, 144124,,, 143824,,, 143524,,, 143224,,, 142924,,, 142624,,, 142324,,, 142024,,, 141724,,, 141424,,, 141124,,, 140824,,, 140524,,, 140224,,, 139924,,, 139624,,, 139324,,, 139024,,, 138724,,, ****** Step #1 ************* *amplitude,name=test 0,0,0.5,1,1,0 *amplitude,name=test2 0,0,0.5,1,1,0.0001 *STEP, NLGEOM, INC=5000 *STATIC,direct 0.01, 1,1.e-12 ** *controls, analysis=discontinuous *controls, parameter=time incrementation , , , , , , , 30, *CONTROLS, PARAMETER=FIELD,FIELD=DISPLACEMENT ,1.0 **contact controls,automatic tolerances *BOUNDARY,amp=test2 230024,2,2,-0.12 *debond, slave=botsurf, master=topsurf *FRACTURE CRITERION,TYPE=VCCT,NODAL ENERGY RATE , , , *output, field,freq=1 *element output s *node output u, *contact output,SLAVE=botsurf, MASTER=topsurf cstress cdisp *output, history,freq=5 *node output, nset=LP u2, rf2 *END STEP *STEP,INC=5000 *direct cyclic,fatigue 0.01,1,,,25,25,,20 3,6,20000,1.1 *debond, slave=botsurf, master=topsurf *FRACTURE CRITERION,TYPE=fatigue,NODAL ENERGY RATE 0.5,-1.44,4.8768E-6,1.15,,0.95,,, ,1.75 *controls,type=direct cyclic ,100,1E-3,100,1E-3 **contact controls,automatic tolerances *BOUNDARY,amp=test 230024,2,2,-0.12 *output,field, freq=50 *element output s *node output u *contact output,SLAVE=botsurf,MASTER=topsurf bdstat *output, history,freq=5 *node output,nset=LP u2 rf2 *END STEP