*HEADING
EXTRUSION OF A BAR WITH FRICTIONAL HEAT GENERATION 
ENHANCED HOURGLASS CONTROL
EXPLICIT [CAX4RT; DIE MODELED WITH RAX2 ELEMENTS]
*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
*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=ENHANCED
*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 
  5001, 0.250000000E+00,-0.180000000E+00
  5002, 0.866000000E-01,-0.180000000E+00
  5003, 0.666000000E-01,-0.180000000E+00
  5004, 0.666000000E-01,-0.170000000E+00
  5005, 0.666000000E-01,-0.150000000E+00
  5006, 0.666000000E-01,-0.122745000E+00
  5007, 0.792341000E-01,-0.810862000E-01
  5008, 0.873649000E-01,-0.689138000E-01
  5009, 0.999990000E-01,-0.272546000E-01
  5010, 0.999990000E-01, 0.000000000E+00
  5011, 0.999990000E-01, 0.300000000E+00
  5012, 0.999990000E-01, 0.310000000E+00
  5013, 0.200000000E+00, 0.410000000E+00
  5014, 0.239788000E+00,-0.180000000E+00
  5015, 0.229575000E+00,-0.180000000E+00
  5016, 0.219362000E+00,-0.180000000E+00
  5017, 0.209150000E+00,-0.180000000E+00
  5018, 0.198938000E+00,-0.180000000E+00
  5019, 0.188725000E+00,-0.180000000E+00
  5020, 0.178512000E+00,-0.180000000E+00
  5021, 0.168300000E+00,-0.180000000E+00
  5022, 0.158088000E+00,-0.180000000E+00
  5023, 0.147875000E+00,-0.180000000E+00
  5024, 0.137663000E+00,-0.180000000E+00
  5025, 0.127450000E+00,-0.180000000E+00
  5026, 0.117238000E+00,-0.180000000E+00
  5027, 0.107025000E+00,-0.180000000E+00
  5028, 0.968125000E-01,-0.180000000E+00
  5029, 0.766000000E-01,-0.180000000E+00
  5030, 0.666000000E-01,-0.160000000E+00
  5031, 0.666000000E-01,-0.140915000E+00
  5032, 0.666000000E-01,-0.131830000E+00
  5033, 0.674114000E-01,-0.111743000E+00
  5034, 0.698280000E-01,-0.100979000E+00
  5035, 0.737975000E-01,-0.906857000E-01
  5036, 0.928015000E-01,-0.593143000E-01
  5037, 0.967710000E-01,-0.490210000E-01
  5038, 0.991876000E-01,-0.382569000E-01
  5039, 0.999990000E-01,-0.181697000E-01
  5040, 0.999990000E-01,-0.908487000E-02
  5041, 0.999990000E-01, 0.100000000E-01
  5042, 0.999990000E-01, 0.200000000E-01
  5043, 0.999990000E-01, 0.300000000E-01
  5044, 0.999990000E-01, 0.400000000E-01
  5045, 0.999990000E-01, 0.500000000E-01
  5046, 0.999990000E-01, 0.600000000E-01
  5047, 0.999990000E-01, 0.700000000E-01
  5048, 0.999990000E-01, 0.800000000E-01
  5049, 0.999990000E-01, 0.900000000E-01
  5050, 0.999990000E-01, 0.100000000E+00
  5051, 0.999990000E-01, 0.110000000E+00
  5052, 0.999990000E-01, 0.120000000E+00
  5053, 0.999990000E-01, 0.130000000E+00
  5054, 0.999990000E-01, 0.140000000E+00
  5055, 0.999990000E-01, 0.150000000E+00
  5056, 0.999990000E-01, 0.160000000E+00
  5057, 0.999990000E-01, 0.170000000E+00
  5058, 0.999990000E-01, 0.180000000E+00
  5059, 0.999990000E-01, 0.190000000E+00
  5060, 0.999990000E-01, 0.200000000E+00
  5061, 0.999990000E-01, 0.210000000E+00
  5062, 0.999990000E-01, 0.220000000E+00
  5063, 0.999990000E-01, 0.230000000E+00
  5064, 0.999990000E-01, 0.240000000E+00
  5065, 0.999990000E-01, 0.250000000E+00
  5066, 0.999990000E-01, 0.260000000E+00
  5067, 0.999990000E-01, 0.270000000E+00
  5068, 0.999990000E-01, 0.280000000E+00
  5069, 0.999990000E-01, 0.290000000E+00
  5070, 0.107142000E+00, 0.317143000E+00
  5071, 0.114285000E+00, 0.324286000E+00
  5072, 0.121428000E+00, 0.331429000E+00
  5073, 0.128571000E+00, 0.338571000E+00
  5074, 0.135714000E+00, 0.345714000E+00
  5075, 0.142857000E+00, 0.352857000E+00
  5076, 0.149999000E+00, 0.360000000E+00
  5077, 0.157142000E+00, 0.367143000E+00
  5078, 0.164285000E+00, 0.374286000E+00
  5079, 0.171428000E+00, 0.381429000E+00
  5080, 0.178571000E+00, 0.388571000E+00
  5081, 0.185714000E+00, 0.395714000E+00
  5082, 0.192857000E+00, 0.402857000E+00
  5083, 0.200000000E+00, 0.000000000E+00
*ELEMENT, TYPE=RAX2,ELSET=EL_DIE                                 
  5001,  5001,  5014
  5002,  5014,  5015
  5003,  5015,  5016
  5004,  5016,  5017
  5005,  5017,  5018
  5006,  5018,  5019
  5007,  5019,  5020
  5008,  5020,  5021
  5009,  5021,  5022
  5010,  5022,  5023
  5011,  5023,  5024
  5012,  5024,  5025
  5013,  5025,  5026
  5014,  5026,  5027
  5015,  5027,  5028
  5016,  5028,  5002
  5017,  5002,  5029
  5018,  5029,  5003
  5019,  5003,  5004
  5020,  5004,  5030
  5021,  5030,  5005
  5022,  5005,  5031
  5023,  5031,  5032
  5024,  5032,  5006
  5025,  5006,  5033
  5026,  5033,  5034
  5027,  5034,  5035
  5028,  5035,  5007
  5029,  5007,  5008
  5030,  5008,  5036
  5031,  5036,  5037
  5032,  5037,  5038
  5033,  5038,  5009
  5034,  5009,  5039
  5035,  5039,  5040
  5036,  5040,  5010
  5037,  5010,  5041
  5038,  5041,  5042
  5039,  5042,  5043
  5040,  5043,  5044
  5041,  5044,  5045
  5042,  5045,  5046
  5043,  5046,  5047
  5044,  5047,  5048
  5045,  5048,  5049
  5046,  5049,  5050
  5047,  5050,  5051
  5048,  5051,  5052
  5049,  5052,  5053
  5050,  5053,  5054
  5051,  5054,  5055
  5052,  5055,  5056
  5053,  5056,  5057
  5054,  5057,  5058
  5055,  5058,  5059
  5056,  5059,  5060
  5057,  5060,  5061
  5058,  5061,  5062
  5059,  5062,  5063
  5060,  5063,  5064
  5061,  5064,  5065
  5062,  5065,  5066
  5063,  5066,  5067
  5064,  5067,  5068
  5065,  5068,  5069
  5066,  5069,  5011
  5067,  5011,  5012
  5068,  5012,  5070
  5069,  5070,  5071
  5070,  5071,  5072
  5071,  5072,  5073
  5072,  5073,  5074
  5073,  5074,  5075
  5074,  5075,  5076
  5075,  5076,  5077
  5076,  5077,  5078
  5077,  5078,  5079
  5078,  5079,  5080
  5079,  5080,  5081
  5080,  5081,  5082
  5081,  5082,  5013
*NSET, NSET=REFNODE
 5083,
*ELEMENT,TYPE=HEATCAP,ELSET=CAP
99001,5083
*HEATCAP,ELSET=CAP
1.0,
**
*NSET,NSET=TOP,GENERATE
61,2061,100
*NSET,NSET=ALL
1,2061,1
*INITIAL CONDITIONS,TYPE=TEMPERATURE
ALL,20. 
REFNODE,20.
*NSET,NSET=NFILEOUT
AXIS,OUTSIDE,TOP
*ELSET,ELSET=EFILEOUT
BOT,SIDE,TOP
*SURFACE,TYPE=ELEMENT,NAME=DIE
EL_DIE,SPOS
*SURFACE,TYPE=ELEMENT, NAME=BAR
BOT, S1
SIDE,S2
*RIGID BODY, REF NODE=5083, ELSET=EL_DIE, ISOTHERMAL=YES
**
**
*ELSET, ELSET=QA_TEST
EFILEOUT,
*NSET, NSET=QA_TEST
NFILEOUT,
*STEP
STABILIZE WORKPIECE INSIDE DIE
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT
 ,1.
*FIXED MASS SCALING,ELSET=WORK,FACTOR=1.E5
*CONTACT PAIR, INTERACTION=CONTACT
BAR,DIE
*SURFACE INTERACTION,NAME=CONTACT
*FRICTION
0.1 , 
*GAP HEAT GENERATION
1.0,0.5
*GAP CONDUCTANCE,PRESSURE
0.0,0.0
0.0,10.0
*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,TIME INTERVAL=0.5,VAR=PRESELECT
**
*OUTPUT, FIELD, TIME MARKS=YES, NUMBER INTER=1
*ELEMENT OUTPUT, ELSET=QA_TEST

PEEQ,
*NODE OUTPUT, NSET=QA_TEST
U,
*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
*FILE OUTPUT, NUM=1
*EL FILE,ELSET=EFILEOUT
PE, 
*NODE FILE,NSET=NFILEOUT
NT, 
*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,NUM=1, VAR=PRESELECT
*OUTPUT,HISTORY,TIME INTERVAL=.5,VAR=PRESELECT
** OUTPUT REQUEST FOR QA
*OUTPUT, FIELD, TIME MARKS=YES, NUMBER INTER=1
*ELEMENT OUTPUT, ELSET=QA_TEST

PEEQ,
*NODE OUTPUT, NSET=QA_TEST
U, 
*OUTPUT,HISTORY,TIME INTERVAL=10,VAR=PRESELECT
**
*END STEP
**
**
*STEP
LET WORKPIECE COOL DOWN--II
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT
 ,9990.
*FIXED MASS SCALING,ELSET=WORK,FACTOR=4.E9
**
*DLOAD,OP=NEW
** Output request for QA
*OUTPUT, FIELD, TIME MARKS=YES, NUMBER INTER=1
*ELEMENT OUTPUT, ELSET=QA_TEST

PEEQ,
*NODE OUTPUT, NSET=QA_TEST
U, 
*OUTPUT,HISTORY,TIME INTERVAL=9990,VAR=PRESELECT
**
*FILE OUTPUT, NUM=1
*EL FILE,ELSET=EFILEOUT
PE, 
*NODE FILE,NSET=NFILEOUT
NT, 
*OUTPUT,FIELD,NUM=1, VAR=PRESELECT
*OUTPUT,HISTORY,TIME INTERVAL=500,VAR=PRESELECT
*END STEP