*HEADING
 FLAT SHEET ROLLING --- 2 ROLLERS, 20% REDUCTION IN EACH ROLL,
 3D ELEMENTS, ANALYTICAL RIGID SURFACE
**
**  Flat sheet.
** 
*NODE
   1,0.,0.,0.
  21,.1,0.,0.
 401,0.,.02,0.
 421,.1,.02,0.
1001,0.,0.,.005
1021,.1,0.,.005
1401,0.,.02,.005
1421,.1,.02,.005
*NGEN,NSET=BOTTOM
   1,  21,1
1001,1021,1
*NGEN,NSET=TOP
 401, 421,1
1401,1421,1
*NFILL,NSET=ALL
BOTTOM,TOP,4,100
*ELEMENT,TYPE=C3D8R
1, 1,2,102,101, 1001,1002,1102,1101
*ELGEN,ELSET=SHEET
1, 20,1,1, 4,100,20
*ELSET,ELSET=TOP,GEN
61,80,1
*ELSET,ELSET=END,GEN
20,80,20
**
*SOLID SECTION,ELSET=SHEET,MATERIAL=METAL
*MATERIAL,NAME=METAL
*ELASTIC
210.E9,.3
*PLASTIC
250.E6, 
*DENSITY
7850., 
*INELASTIC HEAT FRACTION
.9, 
*SPECIFIC HEAT
460.46, 
*EXPANSION
1.E-5, 
*BOUNDARY
BOTTOM,YSYMM
ALL,ZSYMM
*INITIAL CONDITIONS, TYPE=VELOCITY
ALL,1,27.6
*INITIAL CONDITIONS, TYPE=TEMPERATURE
ALL,292.
**
**  Rigid body definition for rollers.
** 
*NODE
1000,.1196,.066,.0025
2000,.1796,.062,.0025
*BOUNDARY
1000,1,6
2000,1,6
**
*NSET,NSET=CNODES
1000, 2000
*ELSET,ELSET=ELHIST
76,
**
** Rigid Surface definition
**
*SURFACE, NAME=ROLLER1,TYPE=CYLINDER
.1196,.066,.010 ,.1296,.066,.010 
.1196,.066,-.005 
START, 0.05,0.0
CIRCL, 0.0,-0.05
CIRCL, -0.05,0.0
CIRCL, 0.0, 0.05
CIRCL, 0.05,0.0
*SURFACE, NAME=ROLLER2,TYPE=CYLINDER
.1796,.062,.010,.1896,.062,.010
.1796,.062,-.005
START, 0.05,0.0
CIRCL, 0.0,-0.05
CIRCL, -0.05,0.0
CIRCL, 0.0, 0.05
CIRCL, 0.05,0.0
*SURFACE,TYPE=ELEMENT,NAME=SHEET
TOP,S5
END,S4
*RIGID BODY,REF NODE=1000, ANALYTICAL SURFACE =ROLLER1
*RIGID BODY,REF NODE=2000, ANALYTICAL SURFACE =ROLLER2
**
*SURFACE INTERACTION,NAME=ROLL_TO_STEEL
*FRICTION
 .3,
*STEP
*DYNAMIC,EXPLICIT,ADIABATIC,ELEMENTBYELEMENT
,2.E-3
*BOUNDARY,TYPE=VELOCITY
1000,6,6,600.
**
*CONTACT 
*CONTACT INCLUSIONS
,
*CONTACT PROPERTY ASSIGNMENT
,,ROLL_TO_STEEL
*RESTART,TIMEMARKS=NO,WRITE,NUM=1
*OUTPUT,FIELD,NUM=5
*ELEMENT OUTPUT
S,E,
*NODE OUTPUT
U,RF
*CONTACT OUTPUT
FSLIP,FSLIPR,CSTRESS
*FILE OUTPUT,TIMEMARKS=YES,NUMBER INTERVAL=1
*EL FILE,ELSET=END
PEEQ,
*NODE FILE,NSET=CNODES
RF,
***********************
*ELSET, ELSET=QA_TEST
END,
*OUTPUT,FIELD,NUMBER INTERVAL=1,TIMEMARKS=YES
*ELEMENT OUTPUT,ELSET=QA_TEST
PEEQ
************************
**
**  History output.
**
***OUTPUT,HISTORY,TIME INTERVAL=1.E-4
***NODE OUTPUT,NSET=CNODES
**UR3,RM3
**RF,
***ELEMENT OUTPUT,ELSET=ELHIST
**MISES,PEEQ,PE,SDV
***ENERGY OUTPUT
***ALLKE,ALLSE,ALLWK,ALLPD,ALLIE,ALLVD,ETOTAL
*END STEP
**
*STEP
*****************
*DYNAMIC,EXPLICIT,ADIABATIC,ELEMENTBYELEMENT
,2.E-3
*BOUNDARY,TYPE=VELOCITY,OP=MOD
2000,6,6,600.
**
*RESTART,TIMEMARKS=NO,WRITE,NUM=2
*END STEP