*HEADING
Analysis of the forming of a channel
SI units (N, kg, m, s)
**
** Nodal coordinates for blank
**
*NODE
  1, 0.,  0.
101, 0.1, 0.
203, 0.,  0.001
303, 0.1, 0.001
*NGEN, NSET=BOT
1,101,1
*NGEN, NSET=TOP
203,303,1
*NFILL
BOT,TOP,2,101
**
** Element connectivity for blank
**
*ELEMENT, TYPE=CPE4I
1, 1,2,103,102
*ELGEN, ELSET=BLANK
1, 100,1,1, 2,101,100
** 
** Node and element sets
**
*NSET, NSET=CENTER
1,102,203
*NSET, NSET=MIDRIGHT, GENERATE
155,202,1
*NSET, NSET=MIDLEFT, GENERATE
102,147,1
*ELSET, ELSET=BLANK_B, GENERATE
1,100,1
*ELSET, ELSET=BLANK_T, GENERATE
101,200,1
**
** Nodal coordinates and rigid body reference node for die
**
*NODE
9001, 0.051, -0.06
9006, 0.051, -0.005
9016, 0.056, 0.
9021, 0.11,  0.
*NGEN
9001,9006,1
9016,9021,1
*NGEN, LINE=C
9006,9016,1,,0.056, -0.005, 0.
*NODE, NSET=REFDIE
9000,0.1,-0.06
**
** Rigid elements for die
**
*ELEMENT, TYPE=R2D2, ELSET=DIE
9001, 9001,9002
*ELGEN, ELSET=DIE
9001, 20,1,1
*RIGID BODY, ELSET=DIE, REF NODE=9000
*NORMAL, TYPE=CONTACT SURFACE
DIE, 9006,-1,0,0
DIE, 9016, 0,1,0
**
** Rigid blankholder
**
*NODE, NSET=REFHOLD
8000, 0.1,0.06
*SURFACE, TYPE=SEGMENTS, NAME=HOLDER, FILLET RADIUS=0.001
START,0.110, 0.001
LINE, 0.056,0.001
CIRCL,0.051,0.006, 0.056,0.006
LINE, 0.051,0.060
*RIGID BODY, ANALYTICAL SURFACE=HOLDER, REF NODE=8000
**
** Rigid punch
**
*NODE, NSET=REFPUNCH
7000, 0.000, 0.061
*SURFACE, TYPE=SEGMENTS, NAME=PUNCH, FILLET RADIUS=0.001
START,0.050, 0.061
LINE, 0.050, 0.007
CIRCL,0.045, 0.002, 0.045,0.007
LINE,-0.010, 0.002
*RIGID BODY, ANALYTICAL SURFACE=PUNCH, REF NODE=7000
**
** Physical and material properties
**
*SOLID SECTION, MATERIAL=STEEL, ORIENTATION=LOCAL, ELSET=BLANK
*ORIENTATION, NAME=LOCAL
1.,0.,0., 0.,1.,0.
3,0
*MATERIAL, NAME=STEEL
*ELASTIC
2.1E11,0.3
*PLASTIC
400.E6, 0.0E-2
420.E6, 2.0E-2
500.E6,20.0E-2
600.E6,50.0E-2
**
** Contact definition
**
*SURFACE, NAME=BLANK_B
BLANK_B, S1
*SURFACE, NAME=BLANK_T
BLANK_T, S3
*SURFACE, NAME=DIE
DIE, SPOS
*CONTACT PAIR, INTERACTION=FRIC
BLANK_B, DIE
*CONTACT PAIR, INTERACTION=FRIC
BLANK_T, HOLDER
*CONTACT PAIR, INTERACTION=NOFRIC
BLANK_T, PUNCH
*SURFACE INTERACTION, NAME=FRIC
*FRICTION
0.1, 
*SURFACE INTERACTION, NAME=NOFRIC
*NSET,NSET=NOUT
 REFDIE, REFHOLD, REFPUNCH
**
** Step 1
**
*STEP, INC=10, NLGEOM, UNSYMM=YES
Push the blankholder and die together 
*STATIC
1.0,1.0
*BOUNDARY
MIDRIGHT  ,2,2
MIDLEFT, 2, 2
CENTER  ,XSYMM
REFDIE  ,1,1
REFDIE  ,2,2, 1.E-8
REFDIE  ,6,6
REFPUNCH,ENCASTRE
REFHOLD ,1,1
REFHOLD ,2,2,-1.E-8
REFHOLD ,6,6
*PRINT, CONTACT=YES
*NODE PRINT, FREQUENCY=0
*EL PRINT, FREQUENCY=0
*OUTPUT, FIELD, FREQ=0
*OUTPUT,HISTORY,FREQ=1
*NODE OUTPUT, NSET=NOUT
 RF, U 
*END STEP
**
** Step 2 
**
*STEP, INC=10, NLGEOM, UNSYMM=YES
Remove the middle constraint
*STATIC
1.0,1.0
*BOUNDARY, OP=NEW
MIDLEFT  ,2,2
CENTER  ,XSYMM
REFDIE  ,1,1
REFDIE  ,2,2, 1.E-8
REFDIE  ,6,6
REFPUNCH,ENCASTRE
REFHOLD ,1,1
REFHOLD ,2,2,-1.E-8
REFHOLD ,6,6
*END STEP
**
** Step 3
**
*STEP, INC=10, NLGEOM, UNSYMM=YES
Apply prescribed force on blankholder
*STATIC
1.0,1.0
*BOUNDARY,OP=NEW
MIDLEFT  ,2,2
CENTER  ,XSYMM
REFDIE  ,1,1
REFDIE  ,2,2, 1.E-8
REFDIE  ,6,6
REFPUNCH,ENCASTRE
REFHOLD ,1,1
REFHOLD ,6,6
*CLOAD
REFHOLD,2,-4.4E5
*END STEP
**
** Step 4 
**
*STEP, INC=20, NLGEOM, UNSYMM=YES
Move punch down a little whilst applying small pressure to
BLANK_T
*STATIC
0.1,1.0 
*BOUNDARY, OP=NEW
CENTER  ,XSYMM
REFDIE  ,1,1
REFDIE  ,2,2, 1.E-8
REFDIE  ,6,6
REFPUNCH,1,1
REFPUNCH,2,2,-1.00000E-3
REFPUNCH,6,6
REFHOLD ,1,1
REFHOLD ,6,6
*DLOAD
BLANK_T, P3, -1000.0E0
*END STEP    
**
** Step 5
**
*STEP, INC=1000, NLGEOM, UNSYMM=YES
Move the punch down to full extent
*STATIC
.0001,1.0,1.E-6
*BOUNDARY
CENTER  ,XSYMM
REFDIE  ,1,1
REFDIE  ,2,2, 1.E-8
REFDIE  ,6,6
REFPUNCH,1,1
REFPUNCH,2,2,-0.03 
REFPUNCH,6,6
REFHOLD ,1,1
REFHOLD ,6,6
*DLOAD, OP=NEW
*OUTPUT, FIELD, FREQ=1000, OP=ADD
*ELEMENT OUTPUT, ELSET=BLANK
 S, PE, PEEQ
*NODE OUTPUT
 U
*CONTACT OUTPUT,MASTER=PUNCH
 CSTRESS
*END STEP