*HEADING
TWO-STAGE FORGING OF A STEEL PART
Units: N, mm, seconds
*NODE,NSET=MOLD, INPUT=ale_forging_steelpartnode1.inp
*NODE,NSET=DIE1, INPUT=ale_forging_steelpartnode2.inp 
*ELEMENT, TYPE=C3D8R, ELSET=MOLD, INPUT=ale_forging_steelpartelem1.inp
*NODE,NSET=PUNCH1
2000001,-1.,-1.,29.
200000,-1.,-1.,29.
200001,44.,-1.,29.
200002,44.,20.,29.
200003,-1.,20.,29.
*ELEMENT, TYPE=R3D4, ELSET=PUNCH1
200000,200000,200003,200002,200001
*NODE,NSET=BASE
3000001,-1.,-1.,.000001
300000,-1.,-1.,.000001
300001,44.,-1.,.000001
300002,44.,20.,.000001
300003,-1.,20.,.000001
*ELEMENT, TYPE=R3D4, ELSET=BASE
300000,300000,300001,300002,300003
*ELEMENT, TYPE=R3D4, ELSET=DIE1,INPUT=ale_forging_steelpartelem2.inp
*NODE,NSET=PUNCH2, INPUT=ale_forging_steelpartnode3.inp 
*ELEMENT, TYPE=R3D4, ELSET=PUNCH2, INPUT=ale_forging_steelpartelem3.inp 
*ELEMENT, TYPE=R3D3, ELSET=PUNCH2, INPUT=ale_forging_steelpartelem4.inp 
*NODE,NSET=DIE2, INPUT=ale_forging_steelpartnode4.inp 
*ELEMENT, TYPE=R3D4, ELSET=DIE2, INPUT=ale_forging_steelpartelem5.inp 
*INCLUDE,INPUT=ale_forging_steelpartsets.inp
*BOUNDARY
1000001,1,6
2000001,1,2
2000001,4,6
3000001,1,6
4000001,1,6
5000001,1,1
5000001,3,6
YSYMM,YSYMM
*SOLID SECTION,ELSET=MOLD,MATERIAL=C15,CONTROLS=SECT
*SECTION CONTROLS,NAME=SECT,HOURGLASS=STIFFNESS, KINEMATICS=ORTHOGONAL
**
*MATERIAL,NAME=C15
*ELASTIC
210000.,0.3
*PLASTIC
   200.000,     0.000
   246.934,     0.010
   434.092,     0.100
   514.439,     0.200
   568.167,     0.300
   609.658,     0.400
   643.916,     0.500
   673.331,     0.600
   699.247,     0.700
   722.501,     0.800
   743.654,     0.900
   763.100,     1.000
   781.129,     1.100
   797.960,     1.200
   813.762,     1.300
   828.672,     1.400
*DENSITY
7.85E-9, 
*RESTART,WRITE,NUM=10
*NSET,NSET=P1
2000001,
*NSET,NSET=P2
4000001,
*SURFACE,TYPE=ELEMENT,NAME=PUNCH2
PUNCH2,SPOS
*SURFACE,TYPE=ELEMENT,NAME=TOP,REGION TYPE=SLIDING
TOP,S2
*SURFACE,TYPE=ELEMENT,NAME=DIE1
DIE1,SPOS
*SURFACE,TYPE=ELEMENT,NAME=DIE2
DIE2,SPOS
*SURFACE,TYPE=ELEMENT,NAME=SIDE,REGION TYPE=SLIDING
SIDE2,S5
SIDE4,S4
*SURFACE,TYPE=ELEMENT,NAME=BOTTOM,REGION TYPE=SLIDING
BOTTOM,S1
*SURFACE,TYPE=ELEMENT,NAME=BASE
BASE,SPOS
*SURFACE,TYPE=ELEMENT,NAME=PUNCH1
PUNCH1,SPOS
*RIGID BODY,ELSET=DIE1,REF NODE=1000001
*RIGID BODY,ELSET=PUNCH1,REF NODE=2000001
*RIGID BODY,ELSET=BASE,REF NODE=3000001
*RIGID BODY,ELSET=PUNCH2,REF NODE=4000001
*RIGID BODY,ELSET=DIE2,REF NODE=5000001
*SURFACE INTERACTION, NAME=FRIC
*FRICTION
0.1,
*NSET, NSET=REFN
1000001, 2000001, 3000001, 4000001, 5000001
*ELSET, ELSET=QA_TEST
TOP,
*NSET, NSET=QA_TEST, ELSET=QA_TEST
*FILTER, NAME=SMOOTH, TYPE=BUTTERWORTH
15000, 2
**
** FIRST FORMING PROCESS
** 
*STEP
*DYNAMIC,EXPLICIT
,.000423333
*BOUNDARY,TYPE=VELOCITY
2000001,3,3,-30000.
5000001,2,2,1181.1
**
*CONTACT
*CONTACT INCLUSIONS, ALL ELEMENT BASED
*CONTACT EXCLUSIONS
PUNCH1,DIE1
 ,PUNCH2
 ,DIE2
*CONTACT PROPERTY ASSIGNMENT
,,FRIC
*SURFACE PROPERTY ASSIGNMENT, PROPERTY=FEATURE EDGE CRITERIA
 ,20.0
**
**Output for QA purposes:
*OUTPUT, FIELD, TIME MARKS=YES, NUMBER INTERVAL=20
*ELEMENT OUTPUT, ELSET=QA_TEST
PEEQ,
*NODE OUTPUT, NSET=QA_TEST
U,
*OUTPUT, HIST, TIME INTERVAL=4E-6
*ENERGY OUTPUT, VAR=PRESELECT
**
**Output for postprocessing:
*OUTPUT, FIELD, VAR=PRESELECT
*OUTPUT, HIST, TIME INTERVAL=4E-6
*NODE OUTPUT, NSET=REFN
RF,
*OUTPUT, HIST, TIME INTERVAL=4E-6, FILTER=SMOOTH
*NODE OUTPUT, NSET=REFN
RF,
**
*END STEP
**
** MOVE DIE2 INTO PLACE
**
*STEP
*DYNAMIC,EXPLICIT
,.00001
*BOUNDARY,TYPE=VELOCITY
5000001,2,2,-50000.0
**
*CONTACT, OP=NEW
*CONTACT INCLUSIONS, ALL ELEMENT BASED
*CONTACT EXCLUSIONS
PUNCH2,DIE2
,PUNCH1
,DIE1
*CONTACT PROPERTY ASSIGNMENT
,,FRIC
*SURFACE PROPERTY ASSIGNMENT, PROPERTY=FEATURE EDGE CRITERIA
 ,20.0
**
*END STEP
**
** SECOND FORMING PROCESS
** 
*STEP
*DYNAMIC,EXPLICIT
,.00035
*BOUNDARY,TYPE=VELOCITY
4000001,3,3,-30000.
5000001,2,2,0.0
**
*END STEP