*HEADING DEEP DRAWING OF CYLINDRICAL CUP WITH S4/S3R DEEPDRAWCUP_S4 *RESTART,WRITE,FREQUENCY=25 *PREPRINT,CONTACT=YES ** Define end nodes in x-z plane *NODE,NSET=ENDS 101, 181,9.951847E-2,-9.801714E-3,0. ** Generate end nodes rotated by theta=11.25 degrees *NCOPY, OLD SET=ENDS, CHANGE NUMBER=1000, SHIFT, NEW SET=ENDS2 0.,0.,0. 0.,0.,0.,0.,0.,1.,11.25 **Generate nodes for blank *NGEN,NSET=BOT 101,181,2 1101,1181,2 *NSET,NSET=PLANE1,GENERATE 103,181,2 *NSET,NSET=PLANE2,GENERATE 1101,1181,2 *NSET,NSET=PLANE PLANE1,PLANE2 *NSET,NSET=WRKPC BOT, **Define rigid body nodes for the die, punch and holder *NODE,NSET=DIE 100,0.,0.,-0.05 *NODE,NSET=PUNCH 200,0.,0.,.05 *NODE,NSET=HOLDER 300,0.,0.,0.05 *NSET,NSET=TOOLS PUNCH,DIE,HOLDER *NSET,NSET=CENTER 101, **Put a triangle in the centre of the wedge *ELEMENT,TYPE=S3R,ELSET=TRIA 201,101,121,1121 **Create quads for the rest of the blank *ELEMENT,TYPE=S4 202,121,123,1123,1121 *ELGEN,ELSET=QUADS 202,30,2,2 *ELSET,ELSET=BLANK TRIA,QUADS **Define contact surfaces between the blank and the holder *ELSET,ELSET=HOLDCON,GENERATE 222,260,2 **Define interface elements between each node of the blank and the die *ELSET,ELSET=DIECON,GENERATE 212,260,2 **Define interface elements between each node of the blank and the punch *ELSET,ELSET=PUNCHCON,GENERATE 201,201,1 202,230,2 **Define the thickness and material properties for the blank *SHELL SECTION,MATERIAL=STEEL,ELSET=BLANK,POISSON=0.5 .00082,5 *MATERIAL,NAME=STEEL *ELASTIC 2.1E11,0.3 *PLASTIC,HARDENING=ISOTROPIC 0.91294E+08, 0.00000E+00 0.10129E+09, 0.21052E-03 0.11129E+09, 0.52686E-03 0.12129E+09, 0.97685E-03 0.13129E+09, 0.15923E-02 0.14129E+09, 0.24090E-02 0.15129E+09, 0.34674E-02 0.16129E+09, 0.48120E-02 0.17129E+09, 0.64921E-02 0.18129E+09, 0.85618E-02 0.19129E+09, 0.11080E-01 0.20129E+09, 0.14110E-01 0.21129E+09, 0.17723E-01 0.22129E+09, 0.21991E-01 0.23129E+09, 0.26994E-01 0.24129E+09, 0.32819E-01 0.25129E+09, 0.39556E-01 0.26129E+09, 0.47301E-01 0.27129E+09, 0.56159E-01 0.28129E+09, 0.66236E-01 0.29129E+09, 0.77648E-01 0.30129E+09, 0.90516E-01 0.31129E+09, 0.10497E+00 0.32129E+09, 0.12114E+00 0.33129E+09, 0.13916E+00 0.34129E+09, 0.15919E+00 0.35129E+09, 0.18138E+00 0.36129E+09, 0.20588E+00 0.37129E+09, 0.23287E+00 0.38129E+09, 0.26252E+00 0.39129E+09, 0.29502E+00 0.40129E+09, 0.33054E+00 0.41129E+09, 0.36929E+00 0.42129E+09, 0.41147E+00 0.43129E+09, 0.45729E+00 0.44129E+09, 0.50696E+00 0.45129E+09, 0.56073E+00 0.46129E+09, 0.61881E+00 0.47129E+09, 0.68145E+00 0.48129E+09, 0.74890E+00 0.49129E+09, 0.82142E+00 0.50129E+09, 0.89928E+00 0.51129E+09, 0.98274E+00 0.52129E+09, 0.10721E+01 **Define rigid surfaces (main surfaces) for the die, holder and punch *rigid body,analytical surface=RDIE,REFNODE=100 *surface,TYPE=REVOLUTION,NAME=RDIE,FILLET RADIUS=0.00541 0.,0.,0.,0.,0.,1. START,0.05125,-0.06041 LINE,0.05125,-0.00041 LINE,0.12,-0.00041 *rigid body,analytical surface=RHOLDER,REFNODE=300 *surface,TYPE=REVOLUTION,NAME=RHOLDER,FILLET RADIUS=0.001 0.,0.,0.,0.,0.,1. START,0.12,0.0004100175 LINE,0.05625,0.0004100175 LINE,0.05625,.06041 *rigid body,analytical surface=RPUNCH,REFNODE=200 *surface,TYPE=REVOLUTION,NAME=RPUNCH,FILLET RADIUS=.01341 0.,0.,0.,0.,0.,1. START,0.05,0.06041 LINE,0.05,2.91782E-3 CIRCL,0.0,.00051,0.0,1.00051 ** define secondary surfaces *surface,NAME=DIE DIECON,S2 *surface,NAME=HOLDER HOLDCON,S1 *surface,NAME=PUNCH PUNCHCON,S1 ** define contact pairs *CONTACT *CONTACT INCLUSIONS DIE,RDIE HOLDER,RHOLDER PUNCH,RPUNCH *CONTACT PROPERTY ASSIGNMENT DIE,RDIE,DIESURF HOLDER,RHOLDER,HOLDSURF PUNCH,RPUNCH,PUNSURF ** *SURFACE INTERACTION,NAME=DIESURF *FRICTION 0.1, *SURFACE INTERACTION,NAME=HOLDSURF *FRICTION 0.1, *SURFACE INTERACTION,NAME=PUNSURF *FRICTION 0.25, **Apply a transform to all nodes at theta=(+/-)5.625 *TRANSFORM,TYPE=C,NSET=PLANE 0.,0.,0.,0.,0.,1. ****************************************************** *STEP, NLGEOM, UNSYMM=YES STEP 1: PUSH THE BLANKHOLDER DOWN BY A PRESCRIBED DISPLACEMENT *STATIC 1.,1. *BOUNDARY CENTER,1,2 CENTER,4,6 WRKPC,2,2 WRKPC,4,4 WRKPC,6,6 DIE,1,6 PUNCH,1,6 HOLDER,1,2 HOLDER,3,3,-1.77E-8 HOLDER,4,6 *MONITOR,NODE=200,DOF=3 *PRINT,CONTACT=YES *NODE PRINT,NSET=TOOLS,FREQUENCY=100 COORD,U RF, *EL PRINT,ELSET=BLANK,FREQUENCY=500 S,E STH, *CONTACT PRINT,SECONDARY=HOLDER,MAIN=RHOLDER,FREQUENCY=500 *CONTACT PRINT,SECONDARY=DIE,MAIN=RDIE,FREQUENCY=500 *CONTACT PRINT,SECONDARY=PUNCH,MAIN=RPUNCH,FREQUENCY=500 *NODE FILE,NSET=TOOLS,FREQUENCY=10 U,RF *OUTPUT, FIELD, FREQ=0 *END STEP ****************************************************** *STEP,NLGEOM, UNSYMM=YES STEP 2: APPLY THE PRESCRIBED FORCE ON THE BLANKHOLDER AND RELEASE THE DISPLACEMENT *STATIC 1.,1. *BOUNDARY,OP=NEW CENTER,1,2 CENTER,4,6 WRKPC,2,2 WRKPC,4,4 WRKPC,6,6 DIE,1,6 PUNCH,1,6 HOLDER,1,2 HOLDER,4,6 *CLOAD HOLDER,3,-3125. *END STEP ****************************************************** *STEP,INC=500,NLGEOM,UNSYMM=YES STEP 3: MOVE THE PUNCH DOWN *STATIC .01,1.,1.E-6 *CONTACT CONTROLS,STIF=1.0 *controls,parameter=time ,,,,,,20 *BOUNDARY,OP=NEW CENTER,1,2 CENTER,4,6 WRKPC,2,2 WRKPC,4,4 WRKPC,6,6 DIE,1,6 PUNCH,1,2 PUNCH,3,3,-.06 PUNCH,4,6 HOLDER,1,2 HOLDER,4,6 *CLOAD HOLDER,3,-3125. *output,field,freq=20 *node output U, *END STEP ****************************************************** *STEP,INC=1,NLGEOM, UNSYMM=YES STEP 4: FIX ALL NODES AND REMOVE THE CONTACT SURFACES *STATIC 1.,1.,1.,1. *BOUNDARY,FIXED WRKPC,1,6 ** *MODEL CHANGE,TYPE=CONTACT PAIR,REMOVE *CONTACT *CONTACT EXCLUSIONS DIE,RDIE HOLDER,RHOLDER PUNCH,RPUNCH *CLOAD,OP=NEW HOLDER,3,0. *output,field,freq=20,var=preselect *END STEP ****************************************************** *STEP,NLGEOM, UNSYMM=YES STEP 5: REPLACE THE BOUNDARY CONDITIONS BY THE REGULAR SET *STATIC .1,1.,1.E-6 *BOUNDARY,OP=NEW CENTER,1,2 CENTER,4,6 WRKPC,2,2 WRKPC,4,4 WRKPC,6,6 181,3,3 *output,field,freq=20 *node output u,rf *element output,var=preselect *contact output,var=preselect *END STEP