*HEADING 2D Draw Bending problem for Numisheet'93 Mild steel, 2.45kN blankholder force Shell elements; Static analysis using rigid surfaces; Both punching and the springback are included. ** *PREPRINT, ECHO=NO **------------------------------------------ blank ** *NODE 1, 351, 0.175 *NGEN, NSET=BLANK1 1,351,1 *NCOPY, OLD=BLANK1,NEW=BLANK2,CHANGE=1000,SHIFT 0., 0.005 *NSET, NSET=BLANK BLANK1,BLANK2 *ELEMENT, TYPE=S4R 1,1,3,1003,1001 *ELGEN, ELSET=BLANK 1,175,2,2 *SHELL SECTION, ELSET=BLANK, MATERIAL=STEEL 0.00078,5 *MATERIAL,NAME=STEEL *DENSITY 7800., *ELASTIC 206.E9,0.3 *PLASTIC .16700E+09,.00000E+00 .19410E+09,.10000E-01 .21913E+09,.20000E-01 .23803E+09,.30000E-01 .25348E+09,.40000E-01 .26668E+09,.50000E-01 .27826E+09,.60000E-01 .28864E+09,.70000E-01 .29807E+09,.80000E-01 .30674E+09,.90000E-01 .31477E+09,.10000E+00 .32226E+09,.11000E+00 .32930E+09,.12000E+00 .33594E+09,.13000E+00 .34224E+09,.14000E+00 .34822E+09,.15000E+00 .35393E+09,.16000E+00 .35940E+09,.17000E+00 .36464E+09,.18000E+00 .36968E+09,.19000E+00 .37454E+09,.20000E+00 .37922E+09,.21000E+00 .38375E+09,.22000E+00 .38814E+09,.23000E+00 .39239E+09,.24000E+00 .39652E+09,.25000E+00 .40053E+09,.26000E+00 .40443E+09,.27000E+00 .40822E+09,.28000E+00 .41193E+09,.29000E+00 .41554E+09,.30000E+00 .41906E+09,.31000E+00 .42251E+09,.32000E+00 .42588E+09,.33000E+00 .42917E+09,.34000E+00 .43240E+09,.35000E+00 .43556E+09,.36000E+00 .43866E+09,.37000E+00 .44169E+09,.38000E+00 .44468E+09,.39000E+00 .44760E+09,.40000E+00 .45047E+09,.41000E+00 .45330E+09,.42000E+00 .45607E+09,.43000E+00 .45880E+09,.44000E+00 .46148E+09,.45000E+00 .46413E+09,.46000E+00 .46672E+09,.47000E+00 .46928E+09,.48000E+00 .47181E+09,.49000E+00 .47429E+09,.50000E+00 *POTENTIAL 1., 1.0402,1.24897,1.07895,1.,1. ** **------------------------------------------ die ** *NODE 2000, 0.026, 0.0025, -0.080 *SURFACE, TYPE=CYLINDER, NAME=DIE 0.0, 0.0, 0.0, 1.0, 0.0, 0.0 0.0, 1.0, 0.0 START, 0.026, -0.08000 LINE, 0.026, -0.00539 CIRCL, 0.031, -0.00039, 0.031, -0.00539 LINE, 0.180, -0.00039 *RIGID BODY,REF NODE=2000,ANALYTICAL SURFACE=DIE ** **------------------------------------------ punch ** *NODE 3000, 0.000, 0.0025, 0.001 *SURFACE,TYPE=CYLINDER,NAME=PUNCH 0.0, 0.0, 0.0, 1.0, 0.0, 0.0 0.0, 1.0, 0.0 START, 0.025, 0.081 LINE, 0.025, 0.006 CIRCL, 0.020, 0.001, 0.020, 0.006 LINE, 0.000, 0.001 *RIGID BODY,REF NODE=3000,ANALYTICAL SURFACE=PUNCH ** **------------------------------------------ holder ** *NODE 4000, 0.0585, 0.0025, 0.010 *SURFACE,TYPE=CYLINDER,NAME=HOLDER 0.0, 0.0, 0.0, 1.0, 0.0, 0.0 0.0, 1.0, 0.0 START, 0.086, 0.01000 LINE, 0.086, 0.00539 CIRCL, 0.081, 0.00039, 0.081, 0.00539 LINE, 0.031, 0.00039 LINE, 0.031, 0.01000 *RIGID BODY,REF NODE=4000,ANALYTICAL SURFACE=HOLDER ** **------------------------------------------ boundary conditions ** *BOUNDARY BLANK, 2,2 BLANK, 4,4 BLANK, 6,6 1, 1,1 1, 5,5 1001, 1,1 1001, 5,5 2000, 1,6 3000, 1,2 3000, 3,3 3000, 4,6 4000, 1,2 4000, 4,6 ** ***RESTART,WRITE, FREQ=100 ** *AMPLITUDE,NAME=APUNCH 0.,0.,0.00466666667,1.,0.00933333333,0. *ELSET,ELSET=CHECK 57,151 *NSET,NSET=CHECKN 1,20,57,2000,3000,4000,1001,1020,1057 ** **------------------------------------------ contact definitions *SURFACE, NAME=TOP BLANK,SPOS *SURFACE, NAME=BOTTOM BLANK,SNEG ** *SURFACE INTERACTION, NAME=ALLCONT *FRICTION 0.144, *SURFACE BEHAVIOR, PRESSURE-OVERCLOSURE=EXPONENTIAL 0.00039, 5.0E6 ** *CONTACT PAIR, INTERACTION=ALLCONT TOP, PUNCH TOP, HOLDER BOTTOM, DIE **------------------------------------------ analysis ** **------------------------------------------ apply the holder force ** *STEP, NLGEOM, INC=1000,Convert Sdi=No move the holder down using displacement boundary condition *STATIC 1e-8 ,1.0e-7 *PRINT, CONTACT=YES *BOUNDARY 4000,3, 3, -140.0e-7 ** *EL FILE, ELSET=CHECK S, E, PE , *NODE FILE, NSET=CHECKN U, RF ** *ENERGY FILE ALLSE,ALLPD,ALLIE,ALLWK,ETOTAL,ALLFD *EL PRINT, FREQ=0 *NODE PRINT, FREQ=0 *CONTACT FILE, FREQ=100, MASTER=HOLDER, SLAVE=TOP CSTRESS, CDISP *CONTACT FILE, FREQ=100, MASTER=DIE, SLAVE=BOTTOM CSTRESS, CDISP *OUTPUT,HISTORY,VARIABLE=PRESELECT,FREQ=2 *OUTPUT,FIELD,FREQ=100 *NODE OUTPUT U,COORD *ELEMENT OUTPUT S,LE,PE,PEEQ ** *END STEP *STEP, NLGEOM, INC=1000,Convert Sdi=No replace the displacement boundary condition with the known force on the holder. *STATIC 1e-6, 1.0e-5 *BOUNDARY, OP=NEW BLANK, 2,2 BLANK, 4,4 BLANK, 6,6 1, 1,1 1, 5,5 1001, 1,1 1001, 5,5 2000, 1,6 3000, 1,2 3000, 3,3 3000, 4,6 4000, 1,2 4000, 4,6 *CLOAD 4000, 3, -175.0 ** *END STEP ** **------------------------------------------ move punch down ** *STEP, NLGEOM, INC=1000, UNSYMM=YES,Convert Sdi=No move the punch down *STATIC **** To run the full punching problem uncomment the data line **** below and comment out the data line after that. **0.00934e-3 , 0.00933333333 0.00934e-3 , 0.004 *BOUNDARY,AMP=APUNCH,TYPE=VELOCITY, OP=NEW 3000,3,3,-15.0 *BOUNDARY,OP=NEW BLANK, 2,2 BLANK, 4,4 BLANK, 6,6 1, 1,1 1, 5,5 1001, 1,1 1001, 5,5 2000, 1,6 3000, 1,2 3000, 4,6 4000, 1,2 4000, 4,6 ** *CONTACT FILE, FREQ=100, MASTER=PUNCH, SLAVE=TOP CSTRESS, CDISP *END STEP ** **------------------------------------------ spring back ** *STEP,NLGEOM,Convert Sdi=No remove the contact pairs to obtain the spring-back of the blank. *STATIC 0.1,1. *BOUNDARY, OP=NEW BLANK,2,2 BLANK, 4,4 BLANK, 6,6 1, 1,1 1, 5,5 1001, 1,1 1001, 5,5 *BOUNDARY, OP=NEW, TYPE=VELOCITY 1,3,3 1001,3,3 2000,1,6 3000,1,6 4000,1,6 *CLOAD, OP=NEW *MODEL CHANGE, REMOVE, TYPE=CONTACT PAIR TOP, PUNCH TOP, HOLDER BOTTOM, DIE *CONTACT FILE, FREQ=0 *END STEP