*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.
405, 0.,  0.001
505, 0.1, 0.001
*NGEN, NSET=BOT
1,101,1
*NGEN, NSET=TOP
405,505,1
*NFILL
BOT,TOP,4,101
**
** Element connectivity for blank
**
*ELEMENT, TYPE=CPE4R
1, 1,2,103,102
*ELGEN, ELSET=BLANK
1, 100,1,1,4,101,100
** 
** Node and element sets
**
*NSET, NSET=CENTER
1,102,203,304,405
*NSET, NSET=MIDRIGHT
303,
*NSET, NSET=MIDLEFT
203,
*ELSET, ELSET=BLANK_B, GENERATE
1,100,1
*ELSET, ELSET=BLANK_T, GENERATE
301,400,1
**
** Nodal coordinates and rigid body reference node for die
**
*NODE, NSET=REFDIE
9000, 0.1, -0.06
*SURFACE, TYPE=SEGMENTS, NAME=DIE, FILLET RADIUS=0.001
START,0.051,-0.060
LINE, 0.051,-0.005
CIRCL,0.056,0.,0.056,-0.005
LINE, 0.11, 0.
*RIGID BODY, ANALYTICAL SURFACE=DIE, REF NODE=9000
**
** 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
**
*NSET,NSET=NOUT
REFDIE, REFHOLD, REFPUNCH
**
** Physical and material properties
**
*ORIENTATION, NAME=LOCAL
1.,0.,0., 0.,1.,0.
1, 0
*SOLID SECTION, MATERIAL=STEEL, ORIENTATION=LOCAL, ELSET=BLANK, CONTROL=EC-1
*SECTION CONTROLS, NAME=EC-1, HOURGLASS=ENHANCED
*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
*CONTACT PAIR, INTERACTION=FRIC
BLANK_B, DIE
BLANK_T, HOLDER
*CONTACT PAIR, INTERACTION=NOFRIC
BLANK_T, PUNCH
*SURFACE INTERACTION, NAME=FRIC
*FRICTION
0.1, 
*SURFACE INTERACTION, NAME=NOFRIC
**
** Step 1
**
*STEP, NLGEOM=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
*OUTPUT, FIELD, FREQ=0
*OUTPUT,HISTORY,FREQ=1
*NODE OUTPUT, NSET=NOUT
 RF, U 
*END STEP
**
** Step 2 
**
*STEP, NLGEOM=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, NLGEOM=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, NLGEOM=YES
Move punch down a little while 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.E-3
REFPUNCH,6,6
REFHOLD ,1,1
REFHOLD ,6,6
*DLOAD
BLANK_T, P3, -1000.0
*END STEP    
**
** Step 5
**
*STEP, NLGEOM=YES, INC=1000
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.031 
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