*HEADING
 Adiabatic stress analysis in cylindrical billet problem using C3D10M
elements
*NODE
*ELEMENT,TYPE=C3D10M,ELSET=BILLET
*SOLID SECTION,ELSET=BILLET,MATERIAL=METAL
*MATERIAL,NAME=METAL
*ELASTIC
200.E9,.3
*PLASTIC
7.E8,0.00
3.7E9,10.0
*RATE DEPENDENT
*SPECIFIC HEAT
586.,
*DENSITY
7833.,
*CONDUCTIVITY
52.,
*EXPANSION
1.2E-5,
*INELASTIC HEAT FRACTION
******
*NSET,NSET=BNDZ,GEN
*NSET,NSET=BNDY,GEN
*NSET,NSET=BNDX,GEN
******
*NODE,NSET=RSNODE
*NSET,NSET=TOP,GEN
*NSET,NSET=TOP1,GEN
*ELSET,ELSET=CONTACT,GEN
******
*rigid body,analytical surface=RSURF,REF NODE=9999
*surface,TYPE=CYLINDER,NAME=RSURF
START,0.,0.
LINE,0.030,0.
*surface,NAME=TOPSURF
CONTACT,S1
*CONTACT PAIR, INTERACTION=SINTER,TYPE=SURFACE
TOPSURF,RSURF
*SURFACE INTERACTION,NAME=SINTER
*SURFACE BEHAVIOR,NO SEPARATION,PENALTY
,,0.1
*FRICTION
0.5
******
*BOUNDARY
BNDX,1
BNDY,2
BNDZ,3
******
*RESTART,WRITE,FREQ=10
*STEP,INC=200,AMPLITUDE=RAMP,NLGEOM,UNSYMM=YES
*STATIC,ADIABATIC
1.5,100.
*BOUNDARY
RSNODE,1,2
RSNODE,4,6
RSNODE,3,,-.009
*PRINT,CONTACT=YES
*MONITOR,NODE=9999,DOF=3
*EL PRINT, ELSET=BILLET,FREQUENCY=200
S,MISES
PEEQ,
PENER,
TEMP,
*NODE PRINT,FREQUENCY=100
U,
RF,
NT,
*CONTACT PRINT,SLAVE=TOPSURF,FREQUENCY=100
*NODE FILE,NSET=RSNODE
RF,
*OUTPUT,FIELD
*NODE OUTPUT,NSET=RSNODE
RF,
*CONTACT FILE,SLAVE=TOPSURF,FREQUENCY=100
*OUTPUT,FIELD,FREQ=100
*CONTACT OUTPUT,VARIABLE=PRESELECT,SLAVE=TOPSURF
*END STEP