*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