*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
** START, 0.026, -0.07999997
** LINE,  0.026, -0.00538997
** CIRCL, 0.031, -0.00038997, 0.031,-0.00538997 
** LINE,  0.180, -0.00038997
*RIGID BODY,REF NODE=2000,ANALYTICAL SURFACE=DIE
**
**------------------------------------------  punch
**
*NODE
 3000, 0.000, 0.0025, 0.001
*NSET,NSET=PUNCH-REF
 3000
*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, 
**
*CONTACT PAIR,INTERACTION=ALLCONT,TYPE=SURFACE
TOP, PUNCH
TOP, HOLDER
BOTTOM, DIE
**------------------------------------------  analysis
**
**------------------------------------------  apply the holder force 
**
*STEP, NLGEOM, INC=1000, unsymm=yes, convert sdi=yes
 move the holder down using displacement boundary condition
*STATIC
1e-8 ,1.0e-7
*PRINT, CONTACT=YES 
*contact controls,stif=0.01,stabilize
*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, MAIN=HOLDER, SECONDARY=TOP
CSTRESS, CDISP 
*CONTACT FILE, FREQ=100, MAIN=DIE, SECONDARY=BOTTOM
CSTRESS, CDISP 
*OUTPUT,HISTORY,VARIABLE=PRESELECT,FREQ=2
*NODE OUTPUT,NSET=PUNCH-REF
 RF,U
*OUTPUT,FIELD,FREQ=100
*NODE OUTPUT
 U,COORD
*ELEMENT OUTPUT
S,LE,PE,PEEQ
**
*END STEP
*STEP, NLGEOM, INC=1000
 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 
 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, MAIN=PUNCH, SECONDARY=TOP
CSTRESS, CDISP 
*END STEP
**
**------------------------------------------  spring back  
**
*STEP,NLGEOM
  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