*HEADING
 FLAT SHEET ROLLING --- 2 ROLLERS, 20% REDUCTION IN EACH ROLL,
 3D ELEMENTS.
**
**  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 roller number 1.
** 
*NODE,NSET=ROLLER1
1000,.1196,.066,.0025
3000,.1196,.066,-.005
3001,.1196,.116,-.005
3091,.0696,.066,-.005
3181,.1196,.016,-.005
4000,.1196,.066,.010
4001,.1196,.116,.010
4091,.0696,.066,.010
4181,.1196,.016,.010
*NGEN,NSET=ROLLER1,LINE=C
3001,3091,1,3000
3091,3181,1,3000
4001,4091,1,4000
4091,4181,1,4000
*ELEMENT,TYPE=R3D4,ELSET=ROLLER1
3001, 3001,3004,4004,4001
*ELGEN,ELSET=ROLLER1
3001, 60,3,1
*BOUNDARY
1000,1,6
**
**  Rigid roller number 2.
** 
*NODE,NSET=ROLLER2
2000,.1796,.062,.0025
5000,.1796,.062,-.005
5001,.1796,.112,-.005
5091,.1296,.062,-.005
5181,.1796,.012,-.005
6000,.1796,.062,.010
6001,.1796,.112,.010
6091,.1296,.062,.010
6181,.1796,.012,.010
*NGEN,NSET=ROLLER2,LINE=C
5001,5091,1,5000
5091,5181,1,5000
6001,6091,1,6000
6091,6181,1,6000
*ELEMENT,TYPE=R3D4,ELSET=ROLLER2
5001, 5001,5004,6004,6001
*ELGEN,ELSET=ROLLER2
5001, 60,3,1
*BOUNDARY
2000,1,6
**
*NSET,NSET=CNODES
 1000, 2000
*ELSET,ELSET=ELHIST
 76,
**
*SURFACE,TYPE=ELEMENT,NAME=ROLLER1
ROLLER1,SPOS
*SURFACE,TYPE=ELEMENT,NAME=ROLLER2
ROLLER2,SPOS
*SURFACE,TYPE=ELEMENT,NAME=SHEET
TOP,S5
END,S4
*RIGID BODY,ELSET=ROLLER1,REF NODE=1000
*RIGID BODY,ELSET=ROLLER2,REF NODE=2000
**
*STEP
*DYNAMIC,EXPLICIT,ADIABATIC
,2.E-3
*BOUNDARY,TYPE=VELOCITY
1000,6,6,600.
*SURFACE INTERACTION,NAME=ROLL_TO_STEEL
*FRICTION,TAUMAX=144.3E6
 .3,
*CONTACT PAIR,INTERACTION=ROLL_TO_STEEL
ROLLER1,SHEET
**
*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,
*NSET,NSET=QA_TEST,ELSET=QA_TEST
*OUTPUT,FIELD,NUMBER INTERVAL=1,TIMEMARKS=YES
*NODE OUTPUT,NSET=QA_TEST
U,
************************
**
**  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
,2.E-3
*RESTART,TIMEMARKS=NO,WRITE,NUM=2
*BOUNDARY,TYPE=VELOCITY,OP=MOD
2000,6,6,600.
*CONTACT PAIR,INTERACTION=ROLL_TO_STEEL
ROLLER2,SHEET
*END STEP