*HEADING
EXTRUSION OF A BAR WITH FRICTIONAL HEAT GENERATION
EXPLICIT [CAX4RT]
PENALTY CONTACT
*NODE
1,0.,0.
61,0.,.3
2001,.1,0.
2061,.1,.3
*NGEN,NSET=AXIS
1,61,1
*NGEN,NSET=OUTSIDE
2001,2061,1
*NFILL,NSET=ALL
AXIS,OUTSIDE,20,100
*NSET, NSET=TEMP
2025, 2027, 2029, 2031
*ELEMENT, TYPE=CAX4RT, ELSET=WORK
1,1,201,203,3
*ELGEN,ELSET=WORK
1,30,2,1,10,200,100
*ELSET,ELSET=BOT,GENERATE
1,901,100
*ELSET,ELSET=SIDE,GENERATE
901,930,1
*ELSET,ELSET=TOP,GENERATE
30,930,100
*SOLID SECTION,ELSET=WORK,MATERIAL=METAL,
CONTROLS=HGLASS
*SECTION CONTROLS,NAME=HGLASS,HOURGLASS=COMBINED
*MATERIAL,NAME=METAL
*ELASTIC
6.9E10,.33
*PLASTIC
** STRAIN RATE APPX .1 **
60.E6,0.0   ,20.
90.E6,.125  ,20.
113.E6,.25  ,20.
124.E6,.375 ,20.
133.E6,0.5  ,20.
165.E6,1.0  ,20.
166.E6,2.0  ,20.
60.E6,0.    ,50.
80.E6,.125  ,50.
97.E6,.25   ,50.
110.E6,.375 ,50.
120.E6,0.5  ,50.
150.E6,1.0  ,50.
151.E6,2.0  ,50.
50.E6,0.0   ,100.
65.E6,.125  ,100,
81.5E6,.25  ,100.
91.E6,.375  ,100.
100.E6,0.5  ,100.
125.E6,1.0  ,100.
126.E6,2.0  ,100.
45.E6,0.0   ,150.
63.E6,.125  ,150.
75.E6,.25   ,150.
89.E6,.5    ,150.
110.E6,1.   ,150.
111.E6,2.   ,150.
*SPECIFIC HEAT
880., 
*DENSITY
2700., 
*CONDUCTIVITY
204.,0.
225.,300.
*EXPANSION,ZERO=20.0 
8.42E-5, 
*INELASTIC HEAT FRACTION
.9, 
*NODE, NSET=REFNODE
9999, 0.2, 0.0, 0.0
*ELEMENT, TYPE=HEATCAP, ELSET=CAP
99001,9999
*HEATCAP, ELSET=CAP
1.,
**
*NSET,NSET=TOP,GENERATE
61,2061,100
*NSET,NSET=ALL
1,2061,1
*INITIAL CONDITIONS,TYPE=TEMPERATURE
ALL,20. 
*NSET,NSET=NFILEOUT
AXIS,OUTSIDE,TOP
*ELSET,ELSET=EFILEOUT
BOT,SIDE,TOP
*SURFACE, TYPE=S, NAME=DIE, FILLET=0.075
START,.25,-.18
LINE,.0866,-.18
LINE,.0666,-.18
LINE,.0666,-.17
LINE,.0666,-.15
LINE,.0666,-.1
LINE,.099999,-.05
LINE,.099999,0.0
LINE,.099999,.3
LINE,.099999,.31
LINE,.2,.41
*SURFACE,TYPE=ELEMENT, NAME=BAR
BOT, S1
SIDE, S2
*RIGID BODY, REFNODE=9999, ISOTHERMAL=YES, 
ANALYTICAL SURFACE =DIE
**
**
*Filter, name=filter_hfl
500,4
*STEP
STABILIZE WORKPIECE INSIDE DIE
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT
 ,1.
*FIXED MASS SCALING, ELSET=WORK, FACTOR=1.E5
*CONTACT PAIR, INTERACTION=CONTACT,
MECHANICAL CONSTRAINT=PENALTY
BAR, DIE
*SURFACE INTERACTION, NAME=CONTACT
*FRICTION
0.1 , 
*BOUNDARY, TYPE=VELOCITY
REFNODE,1,6,0.0
AXIS,1,1,0.0
TOP,2,2,-.000125
2061,1,1,0.0
*BOUNDARY 
ALL,11,11,20.0
REFNODE,11,11,20.0
*OUTPUT,FIELD,NUM=1,VAR=PRESELECT
*OUTPUT,HISTORY,FREQ=50
*NODE OUTPUT, NSET=TEMP
NT, 
*Output, field, num=10, filter=filter_hfl
*ELEMENT OUTPUT,ELSET=WORK
HFL
*END STEP
**
**
*STEP
EXTRUSION
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT
 ,10.
*FIXED MASS SCALING, ELSET=WORK, FACTOR=1.E5
*BOUNDARY,OP=NEW, TYPE=VELOCITY 
REFNODE,1,6,0.0
AXIS,1,1,0.0
TOP,2,2,-.0249875
2061,1,1,0.0
*BOUNDARY,OP=NEW
REFNODE,11,11,20.0
*ELSET, ELSET=QA_TEST
EFILEOUT,
*NSET, NSET=QA_TEST
NFILEOUT,
*OUTPUT, FIELD, TIME MARKS=YES, NUMBER INTER=1
*ELEMENT OUTPUT, ELSET=QA_TEST

PEEQ,
*NODE OUTPUT, NSET=QA_TEST
U, 
*FILE OUTPUT, NUM=1
*EL FILE,ELSET=EFILEOUT
PE,
*NODE FILE,NSET=NFILEOUT
NT,
*OUTPUT,HISTORY,FREQ=50
*NODE OUTPUT, NSET=TEMP
NT, 
*OUTPUT,FIELD,NUM=1
*ELEMENT OUTPUT,ELSET=WORK
S,PEEQ, TEMP, HFL
*NODE OUTPUT
NT, U
*Output, field, num=10, filter=filter_hfl
*ELEMENT OUTPUT,ELSET=WORK
HFL
*END STEP
**
**
*STEP
REMOVE CONTACT 
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT
 ,.1
*FIXED MASS SCALING, ELSET=WORK, FACTOR=1.E5
*CONTACT PAIR, OP=DELETE
BAR,DIE
*BOUNDARY,OP=NEW, TYPE=VELOCITY 
REFNODE,1,6,0.0
AXIS,1,1,0.0
TOP,2,2,0.0
*END STEP
**
**
*STEP
LET WORKPIECE COOL DOWN--I (ADD VISCOUS PRESSURE)
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT
 ,10.
*FIXED MASS SCALING,ELSET=WORK,FACTOR=1.E6 
*AMPLITUDE, NAME=RAMP,TIME=TOTAL TIME
0.0,0.0,11.1,0.0,10011.1,1.0
*FILM, AMP=RAMP
BOT,F1,20.,10.
TOP,F3,20.,10.
SIDE,F2,20.,10.
**
*DLOAD
WORK,VP1,1.E8
WORK,VP2,1.E8
WORK,VP3,1.E8
WORK,VP4,1.E8
**
*OUTPUT, FIELD, TIME MARKS=YES, NUMBER INTER=1
*ELEMENT OUTPUT, ELSET=QA_TEST

PEEQ,
*NODE OUTPUT, NSET=QA_TEST
U, 
*FILE OUTPUT, NUM=1
*EL FILE, ELSET=EFILEOUT
PE, 
*NODE FILE, NSET=NFILEOUT
NT, 
*OUTPUT,HISTORY,FREQ=50
*NODE OUTPUT, NSET=TEMP
NT, 
*OUTPUT, FIELD, NUM=1
*ELEMENT OUTPUT,ELSET=WORK
S,PEEQ, TEMP, HFL
*NODE OUTPUT
NT, U
*Output, field, num=10, filter=filter_hfl
*ELEMENT OUTPUT,ELSET=WORK
HFL
*END STEP
**
**
*STEP
LET WORKPIECE COOL DOWN--II
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT
 ,9990.
*FIXED MASS SCALING,ELSET=WORK,FACTOR=4.E9
**
*DLOAD,OP=NEW
**
*END STEP