*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,6
YSYMM,YSYMM
RIG1,YSYMM
RIG2,YSYMM
RIG3,YSYMM
RIG4,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
***
*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
*STEP
*DYNAMIC,EXPLICIT
,.000423333
*BOUNDARY,TYPE=VELOCITY
2000001,3,3,-30000.
**
*CONTACT PAIR,INTERACTION=FRIC
BOTTOM,BASE
BOTTOM,DIE1
SIDE,DIE1
TOP,PUNCH1
**
*SURFACE INTERACTION, NAME=FRIC
*FRICTION
0.1, 
**
**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
**
**
*STEP
*DYNAMIC,EXPLICIT
,.00035
*BOUNDARY,TYPE=VELOCITY
4000001,3,3,-30000.
**
*CONTACT PAIR,INTERACTION=FRIC,OP=DELETE
BOTTOM,DIE1
SIDE,DIE1
TOP,PUNCH1
**
*CONTACT PAIR,INTERACTION=FRIC1
BOTTOM,DIE2
SIDE,DIE2
TOP,PUNCH2
**
*SURFACE INTERACTION, NAME=FRIC1
*FRICTION
0.1, 
**
*END STEP