*HEADING
- AXISYMMETRIC UPSETTING PROBLEM
- WITH RIGID SURFACE
- SECTION CONTROLS USED (HOURGLASS=STIFFNESS)
- FAST UPSETTING
- PENALTY CONTACT
*NODE
1,
13,.01
1201,0.,.015
1213,.01,.015
*NGEN,NSET=MIDDLE
1,13
*NGEN,NSET=TOP
1201,1213
*NSET, NSET=SIDE, GEN
13, 1213, 100
*NFILL
MIDDLE,TOP,12,100
*NSET,NSET=AXIS,GEN
1,1201,100
*ELEMENT,TYPE=cax4rt,ELSET=BILLET
1,1,2,102,101
*ELGEN,ELSET=BILLET
1,12,1,1,12,100,100
*NODE, NSET=NRIGID
2003,0.01,.02
*SOLID SECTION,ELSET=BILLET,MATERIAL=METAL,CONTROL=B
*SECTION CONTROLS, HOURGLASS=STIFFNESS, NAME=B
*MATERIAL,NAME=METAL
*ELASTIC
200.E9,.3
*PLASTIC
7.E8,0.00
3.7E9,10.0
*RATE DEPENDENT
40.,5.
*SPECIFIC HEAT
586., 
*DENSITY
7833., 
*CONDUCTIVITY
52., 
*EXPANSION
1.2E-5, 
*INELASTIC HEAT FRACTION
0.9, 
*BOUNDARY
MIDDLE,2
AXIS,1
2003,1
2003,3,6
*SURFACE, NAME = BILLET
TOP,S3
SIDE,S2
*SURFACE, NAME=RIGID, TYPE=SEGMENTS
START, 0.02,.015
LINE, 0.00,.015
*RIGID BODY, REF NODE=2003, ANALYTICAL SURFACE =RIGID
*ELSET,ELSET=SIDE,GEN
12,1112,100
*STEP
*DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT
 ,0.1
*FIXED MASS SCALING,ELSET=BILLET,FACTOR=1.E+4
*BOUNDARY,TYPE=VELOCITY
2003,2,,-9.e-2
*ELSET,ELSET=TOP,GEN
1101,1112,1
*SURFACE INTERACTION,NAME=RIG_BILL
*FRICTION
 1.0,
*CONTACT PAIR,INTERACTION=RIG_BILL,
MECHANICAL CONSTRAINT=PENALTY
RIGID,BILLET
*MONITOR,NODE=2003,DOF=2
*FILE OUTPUT,NUM=1
*EL FILE, ELSET=QA_TEST
PEEQ,
*NODE FILE, NSET=NRIGID
U,RF
*NODE FILE, NSET=QA_TEST
NT,
*OUTPUT,FIELD,VAR=PRESELECT,NUM=5
*OUTPUT,HISTORY
*NODE OUTPUT,NSET=NRIGID
U2,RF2
*ENERGY OUTPUT
ALLIE,ALLKE,ALLAE,ALLPD
*****
*NSET, NSET=QA_TEST_REFN
NRIGID,
*ELSET, ELSET=QA_TEST
SIDE,
*NSET, NSET=QA_TEST
SIDE,
*OUTPUT, FIELD, TIME MARKS=YES, NUMBER INTERVAL=1
*ELEMENT OUTPUT, ELSET=QA_TEST

PEEQ,
*NODE OUTPUT, NSET=QA_TEST
U,
*NODE OUTPUT, NSET=QA_TEST_REFN
RF,
*OUTPUT, HIST, FREQ=9999
*ENERGY OUTPUT, VAR=PRESELECT
******
*END STEP