*HEADING Predict postbuckling response of the delaminated composite plate using Abaqus/Explicit VCCT ** ** Units: in, lbf, psi ** *Parameter ** Density density1 = 0.191e0 ** User input data defining the interface properties ** Cohesive elastic modulus (lbf/in**2) Emod = 1.0E6 ** Fracture toughness: GIc = 0.46863 GIIc = 3.171825 GIIIc = 3.171825 ** B-K parameter: eta = 1.75 ** *NODE ** Center points for upper plate 1, 0.0, 0.0, 0.0005 65, 0.0, 0.0, 0.0005 ** Delamination front for upper plate 801, 1.25, 0.0 865, 1.25, 0.0 ** ** Center points for lower plate 100001, 0.0, 0.0, -0.0005 100065, 0.0, 0.0, -0.0005 ** ** Plate - center 2401, 2.25, 0.0 2409, 2.25, 2.25 2425, -2.25, 2.25 2441, -2.25, -2.25 2457, 2.25, -2.25 2465, 2.25, 0.0 ** Plate - upper sticking out piece 10409, 2.25, 4.5, -0.001 10425, -2.25, 4.5, -0.001 ** Lower 10441, -2.25, -4.5, -0.001 10457, 2.25, -4.5, -0.001 ** 999998, 0.0, 0.0, -0.0005 999999, 0.0, 0.0, 0.0005 ** *Ngen, nset=CENTER_UP 1, 65, 1 *Ngen, nset=CENTER_BOT 100001, 100065, 1 *Ngen, nset=DELAM_FRONT_UP, line=C 801, 865, 1, 999998, , , , 0., 0., 1. ** *Ncopy, change number=100000, old set=DELAM_FRONT_UP, new set=DELAM_FRONT_BOT, shift 0.0, 0.0, 0.0 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0 *Ngen, nset=PLATE_EDGE 2401, 2409, 1 2409, 2425, 1 2425, 2441, 1 2441, 2457, 1 2457, 2465, 1 ** **For tie *Nset, nset=STARTING_EDGE_UP, generate 101, 2401, 100 *Nset, nset=ENDING_EDGE_UP, generate 165, 2465, 100 *Nset, nset=STARTING_EDGE_BOT, generate 100101, 102401, 100 *Nset, nset=ENDING_EDGE_BOT, generate 100165, 102465, 100 ** *Ngen, nset=TOP_EDGE 10409, 10425, 1 *Ngen, nset=BOTTOM_EDGE 10441, 10457, 1 *Nset, nset=EDGE1, generate 2409, 2425, 1 *Nset, nset=EDGE2, generate 2441, 2457, 1 *Nset, nset=EDGE3-UP, generate 2402, 2408, 1 2458, 2465, 1 *Nset, nset=EDGE3-BOT, generate 102402, 102408, 1 102458, 102465, 1 *Nset, nset=EDGE4-UP, generate 2426, 2440, 1 *Nset, nset=EDGE4-BOT, generate 102426, 102440, 1 ** *Nfill, nset=DELAM_FACE_UP, bias=1.2 CENTER_UP, DELAM_FRONT_UP, 8,100 *Nfill, nset=DELAM_FACE_BOT, bias=1.2 CENTER_BOT, DELAM_FRONT_BOT, 8,100 *Nfill, nset=UPPER_PLATE DELAM_FRONT_UP, PLATE_EDGE, 16,100 *Nfill, nset=UPPER_PLATE_OUT EDGE1, TOP_EDGE, 8, 1000 EDGE2, BOTTOM_EDGE, 8, 1000 ** *Ncopy, change number=100000, old set=UPPER_PLATE, new set=LOWER_PLATE, shift 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0 *Ncopy, change number=100000, old set=UPPER_PLATE_OUT, new set=LOWER_PLATE_OUT, shift 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0 ** *MPC TIE, CENTER_UP, 999999 TIE, STARTING_EDGE_UP, ENDING_EDGE_UP TIE, UPPER_PLATE_OUT, LOWER_PLATE_OUT TIE, EDGE3-UP, EDGE3-BOT TIE, EDGE4-UP, EDGE4-BOT ** *Nset, nset=FIXED_EDGE_BOT, generate 110441, 110457, 1 30441, 40441, 10000 30457, 40457, 10000 ** *Nset, nset=CLAMPED 110449 *Nset, nset=LOADED_EDGE_BOT, generate 110409, 110425, 1 30409, 40409, 10000 30425, 40425, 10000 *Nset, nset=STABLIZED_EDGE, generate 102402, 102408, 1 102409, 110409, 1000 102425, 110425, 1000 102426, 102440, 1 102441, 110441, 1000 102457, 110457, 1000 102458, 102465, 1 ** ** *Nset, nset=LEFT_1, generate 2425, 10425, 1000 2425, 2441, 1 2441, 10441, 1000 *Ncopy, change number=20000, old set=LEFT_1, shift, new set=LEFT_2 -0.125, 0.0, 0.0 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0 *Ncopy, change number=30000, old set=LEFT_1, shift, new set=LEFT_3 -.25, 0.0, 0.0 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0 ** *Element, type=S4 100001, 2425, 3425, 23425,22425 110001, 22425,23425, 33425,32425 100011, 3441, 2441, 22441,23441 110011, 23441,22441, 32441,33441 100021, 2426, 2425, 22425,22426 110021, 22426,22425, 32425,32426 ** *Elgen, elset=FIX_UPPER 100001, 8,1000,1000 110001, 8,1000,1000 100011, 8,1000,1000 110011, 8,1000,1000 100021, 16, 1, 1 110021, 16, 1, 1 ** ** *Nset, nset=RIGHT_1, generate 2409, 10409, 1000 2401, 2409, 1 2457, 2465, 1 2457, 10457, 1000 ** *MPC TIE, 22401, 22465 TIE, 32401, 32465 ** *Ncopy, change number=20000, old set=RIGHT_1, shift 0.125, 0.0, 0.0 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0 *Ncopy, change number=30000, old set=RIGHT_1, shift 0.25, 0.0, 0.0 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0 ** *Element, type=S4 100101, 2409,22409, 23409, 3409 110101, 22409,32409, 33409,23409 100111, 3457,23457, 22457, 2457 110111, 22457,23457, 33457,32457 100121, 2401,22401, 22402, 2402 110121, 22401,32401, 32402,22402 100131, 2457,22457, 22458, 2458 110131, 22457,32457, 32458,22458 ** *Elgen, elset=FIX_UPPER 100101, 8,1000,1000 110101, 8,1000,1000 100111, 8,1000,1000 110111, 8,1000,1000 100121, 8, 1, 1 110121, 8, 1, 1 100131, 8, 1, 1 110131, 8, 1, 1 ** *Elcopy, element shift=100000, old set=FIX_UPPER, new set=FIX_LOWER, shift nodes=0 ** *Element, type=S3 1000001, 1,101,102 2000001, 100001,100101,100102 *Element, type=S4 1000101, 101,201,202,102 1000801, 801,901,902,802 2000101, 100101,100201,100202,100102 2000801, 100801,100901,100902,100802 1100001, 2409,3409,3410,2410 1200001, 2441,3441,3442,2442 ** 2100001, 102409,103409,103410,102410 2200001, 102441,103441,103442,102442 *Elgen, elset=UPPER_DELAM 1000001, 64,1,1 1000101, 64,1,1,7,100,100 *Elgen, elset=UPPER_IN 1000801, 64,1,1,16,100,100 *Elgen, elset=UPPER_OUT 1100001, 16,1,1,8,1000,1000 1200001, 16,1,1,8,1000,1000 *Elgen, elset=LOWER_DELAM 2000001, 64,1,1 2000101, 64,1,1,7,100,100 *Elgen, elset=LOWER_IN 2000801, 64,1,1,16,100,100 *Elgen, elset=LOWER_OUT 2100001, 16,1,1,8,1000,1000 2200001, 16,1,1,8,1000,1000 *Element, type=s4 2000101,100165,100265,100202,100102 2000201,100265,100365,100302,100202 2000301,100365,100465,100402,100302 2000401,100465,100565,100502,100402 2000501,100565,100665,100602,100502 2000601,100665,100765,100702,100602 2000701,100765,100865,100802,100702 2000801,100865,100965,100902,100802 2000901,100965,101065,101002,100902 2001001,101065,101165,101102,101002 2001101,101165,101265,101202,101102 2001201,101265,101365,101302,101202 2001301,101365,101465,101402,101302 2001401,101465,101565,101502,101402 2001501,101565,101665,101602,101502 2001601,101665,101765,101702,101602 2001701,101765,101865,101802,101702 2001801,101865,101965,101902,101802 2001901,101965,102065,102002,101902 2002001,102065,102165,102102,102002 2002101,102165,102265,102202,102102 2002201,102265,102365,102302,102202 2002301,102365,102465,102402,102302 *Element, type=s3 2000001,100001, 100101, 100102 2000002,100001, 100102, 100103 2000003,100001, 100103, 100104 2000004,100001, 100104, 100105 2000005,100001, 100105, 100106 2000006,100001, 100106, 100107 2000007,100001, 100107, 100108 2000008,100001, 100108, 100109 2000009,100001, 100109, 100110 2000010,100001, 100110, 100111 2000011,100001, 100111, 100112 2000012,100001, 100112, 100113 2000013,100001, 100113, 100114 2000014,100001, 100114, 100115 2000015,100001, 100115, 100116 2000016,100001, 100116, 100117 2000017,100001, 100117, 100118 2000018,100001, 100118, 100119 2000019,100001, 100119, 100120 2000020,100001, 100120, 100121 2000021,100001, 100121, 100122 2000022,100001, 100122, 100123 2000023,100001, 100123, 100124 2000024,100001, 100124, 100125 2000025,100001, 100125, 100126 2000026,100001, 100126, 100127 2000027,100001, 100127, 100128 2000028,100001, 100128, 100129 2000029,100001, 100129, 100130 2000030,100001, 100130, 100131 2000031,100001, 100131, 100132 2000032,100001, 100132, 100133 2000033,100001, 100133, 100134 2000034,100001, 100134, 100135 2000035,100001, 100135, 100136 2000036,100001, 100136, 100137 2000037,100001, 100137, 100138 2000038,100001, 100138, 100139 2000039,100001, 100139, 100140 2000040,100001, 100140, 100141 2000041,100001, 100141, 100142 2000042,100001, 100142, 100143 2000043,100001, 100143, 100144 2000044,100001, 100144, 100145 2000045,100001, 100145, 100146 2000046,100001, 100146, 100147 2000047,100001, 100147, 100148 2000048,100001, 100148, 100149 2000049,100001, 100149, 100150 2000050,100001, 100150, 100151 2000051,100001, 100151, 100152 2000052,100001, 100152, 100153 2000053,100001, 100153, 100154 2000054,100001, 100154, 100155 2000055,100001, 100155, 100156 2000056,100001, 100156, 100157 2000057,100001, 100157, 100158 2000058,100001, 100158, 100159 2000059,100001, 100159, 100160 2000060,100001, 100160, 100161 2000061,100001, 100161, 100162 2000062,100001, 100162, 100163 2000063,100001, 100163, 100164 2000064,100001, 100164, 100165 *Elset, elset=LOWER_PLATE LOWER_DELAM, LOWER_IN, LOWER_OUT, FIX_LOWER *Elset, elset=LOWER_IN2, generate 2000801,2000864 2000901,2000964 2001001,2001064 2001101,2001164 2001201,2001264 2001301,2001364 2001401,2001464 2001501,2001564 2001601,2001664 2001701,2001764 2001801,2001864 2001901,2001964 2002001,2002064 2002101,2002164 2002201,2002264 *Elset, elset=UPPER_PLATE UPPER_DELAM, UPPER_IN, UPPER_OUT, FIX_UPPER *Elset, elset=ESLAVE LOWER_DELAM, LOWER_IN2, ** *Surface, name=SLAVE ESLAVE, SPOS *Surface, name=MASTER UPPER_PLATE, SNEG *Elset, elset=ETIED, generate 2000801,2000864 2000901,2000964 2001001,2001064 2001101,2001164 2001201,2001264 2001301,2001364 2001401,2001464 2001501,2001564 2001601,2001664 2001701,2001764 2001801,2001864 2001901,2001964 2002001,2002064 2002101,2002164 2002201,2002264 *Nset, nset=NTIED, elset=ETIED ******************************************************** *Contact clearance, name=CLEAR56, search nset=NTIED *Surface interaction, name=SURFS *Cohesive behavior, eligibility=ORIGINAL CONTACTS , , *Fracture criterion, type=VCCT, mixed mode behavior=BK , , , ******************************************************** *Orientation, system=R, name=OR_0 1.0, 0.0, 0.0, 0.0, 1.0, 0.0 3, 90. *Orientation, system=R, name=or_90 1.0, 0.0, 0.0, 0.0, 1.0, 0.0 3, 180. *Orientation, system=R, name=or_pos_45 1.0, 0.0, 0.0, 0.0, 1.0, 0.0 3, 135. *Orientation, system=R, name=or_neg_45 1.0, 0.0, 0.0, 0.0, 1.0, 0.0 3, 45. *Orientation, system=R, name=or_pos_60 1.0, 0.0, 0.0, 0.0, 1.0, 0.0 3, 150. *Orientation, system=R, name=or_neg_60 1.0, 0.0, 0.0, 0.0, 1.0, 0.0 3, 30. *Orientation, system=R, name=or_pos_15 1.0, 0.0, 0.0, 0.0, 1.0, 0.0 3, 105. *Orientation, system=R, name=or_neg_15 1.0, 0.0, 0.0, 0.0, 1.0, 0.0 3, 75. ** *Shell section, elset=LOWER_PLATE, offset=SPOS, composite 0.005315, 3, AS4_3501_6, OR_POS_45 0.005315, 3, AS4_3501_6, OR_NEG_45 0.005315, 3, AS4_3501_6, OR_90 0.005315, 3, AS4_3501_6, OR_0 0.005315, 3, AS4_3501_6, OR_POS_45 0.005315, 3, AS4_3501_6, OR_NEG_45 0.005315, 3, AS4_3501_6, OR_90 0.005315, 3, AS4_3501_6, OR_0 0.005315, 3, AS4_3501_6, OR_POS_60 0.005315, 3, AS4_3501_6, OR_NEG_60 0.005315, 3, AS4_3501_6, OR_POS_15 0.005315, 3, AS4_3501_6, OR_NEG_15 0.005315, 3, AS4_3501_6, OR_NEG_15 0.005315, 3, AS4_3501_6, OR_POS_15 0.005315, 3, AS4_3501_6, OR_NEG_60 0.005315, 3, AS4_3501_6, OR_POS_60 0.005315, 3, AS4_3501_6, OR_0 0.005315, 3, AS4_3501_6, OR_90 0.005315, 3, AS4_3501_6, OR_NEG_45 ** *Shell section, elset=UPPER_PLATE, offset=SNEG, composite 0.005315, 3, AS4_3501_6, OR_POS_45 0.005315, 3, AS4_3501_6, OR_0 0.005315, 3, AS4_3501_6, OR_90 0.005315, 3, AS4_3501_6, OR_NEG_45 0.005315, 3, AS4_3501_6, OR_POS_45 ** ** *Material, name=AS4_3501_6 *Elastic, type=LAMINA 18.50E+6, 1.640E+6, 0.30, 0.871E+6, 0.871E+6, 5.220E+5 *Density ** ** *Shell section, elset=UPPER_DELAM, offset=SNEG, composite 0.005 , 3, TEFLON , OR_0 0.005315, 3, AS4_3501_6, OR_POS_45 0.005315, 3, AS4_3501_6, OR_0 0.005315, 3, AS4_3501_6, OR_90 0.005315, 3, AS4_3501_6, OR_NEG_45 0.005315, 3, AS4_3501_6, OR_POS_45 ** *Material, name=TEFLON *Elastic, type=LAMINA 18.50E+4, 1.640E+4, 0.30, 0.871E+4, 0.871E+4, 5.220E+3 *Density ** ** *Elset, elset=UPPER_CON, generate 1000001,1000064,1 1000101,1000164,1 1000201,1000264,1 1000301,1000364,1 1000401,1000464,1 1000501,1000564,1 *Elset, elset=LOWER_CON, generate 2000001,2000064,1 2000101,2000164,1 2000201,2000264,1 2000301,2000364,1 2000401,2000464,1 2000501,2000564,1 ** ** *Elset, elset=PRINT 207001,217001,207101,217101, 2107001,2107002,2107003,2107004,2107005,2107006,2107007,2107008, 2107009,2107010,2107011,2107012,2107013,2107014,2107015,2107016 *Amplitude, definition=SMOOTH STEP, name=LOAD 0.0, 0.0, 0.75, 1.0 **------------------------------------------------ ** *NSET,NSET=QA_TEST DELAM_FRONT_BOT,DELAM_FRONT_UP, DELAM_FACE_BOT,DELAM_FACE_UP *Step, nlgeom=YES *Dynamic, Explicit , 0.75 *Variable mass scaling, type=BELOW MIN, dt=5.e-6, freq=10 **------------------------------------------------ *Contact *Contact inclusions SLAVE, MASTER *Contact clearance assignment SLAVE, MASTER, CLEAR56 *Contact property assignment SLAVE, MASTER, SURFS **------------------------------------------------ *Boundary, amplitude=LOAD CLAMPED, 1, 1 STABLIZED_EDGE, 3, 3 FIXED_EDGE_BOT, 2, 6 LOADED_EDGE_BOT, 2, 2, -0.03 LOADED_EDGE_BOT, 3, 6, 0.0 **------------------------------------------------ *Output, field, number interval=200 *Node output U *Contact output DBT, DBSF, DBS, OPENBC, CRSTS, ENRRT, EFENRRTR, BDSTAT, CFORCE ** *Output, history, freq=1 *Node output, nset=LOADED_EDGE_BOT U2, RF2 ***QA purpose*** *Output, field, num =4 *Node output, nset=QA_TEST UT *Output, history *Energy output, variable=PRESELECT *End step