*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,CONTROLS=HGC
*SECTION CONTROLS,NAME=HGC,HOURGLASS=COMBINED
*MATERIAL,NAME=METAL
*ELASTIC
210.E9,.3
*PLASTIC
250.E6,0.0
300.E6,1.0
*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, direct
1.E-8, 2.E-3
*BOUNDARY,TYPE=VELOCITY
1000,6,6,600.
*CONTACT 
*CONTACT INCLUSIONS
,
*CONTACT PROPERTY ASSIGNMENT
,,ROLL_TO_STEEL
*RESTART,WRITE,NUM=1
*OUTPUT,FIELD,NUM=8
*ELEMENT OUTPUT
S,E
*NODE OUTPUT
U,RF
*CONTACT OUTPUT
FSLIP,FSLIPR,CSTRESS
*FILE OUTPUT,NUMBER INTERVAL=1
*EL FILE,ELSET=END
PEEQ,
*NODE FILE,NSET=CNODES
RF,
***********************
*ELSET, ELSET=QA_TEST
END,
*OUTPUT,FIELD,VAR=PRESELECT
*OUTPUT,FIELD,NUMBER INTERVAL=1
*ELEMENT OUTPUT,ELSET=QA_TEST
PEEQ,LE
************************
**
**  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, direct
1.E-8, 1.E-3
*BOUNDARY,TYPE=VELOCITY,OP=MOD
2000,6,6,600.
*RESTART,TIMEMARKS=NO,WRITE,NUM=2
*END STEP