*HEADING
 FLAT SHEET ROLLING --- 2 ROLLERS, 20% REDUCTION IN EACH ROLL,
 3D ELEMENTS, ANALYTICAL RIGID SURFACE
*PREPRINT, MODEL=YES
*IMPORT, UPDATE=YES, STATE=YES, STEP=1, LIB=xx1_anls
SHEET,
**
*ELSET,ELSET=TOP,GEN
61,80,1
*ELSET,ELSET=END,GEN
20,80,20
*NODE
1000,.1196,.066,.0025
2000,.1796,.062,.0025
*BOUNDARY
1000,1,6
2000,1,6
*BOUNDARY
BOTTOM,YSYMM
ALL,ZSYMM
**
*NSET,NSET=CNODES
1000, 2000
**
**  Rigid body definition for rollers.
** 
*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, 1.E-3
*BOUNDARY,TYPE=VELOCITY,OP=MOD
1000,6,6,600.
2000,6,6,600.
*Contact
*Contact Inclusions, All Exterior
*CONTACT PROPERTY ASSIGNMENT
,,ROLL_TO_STEEL
*OUTPUT,FIELD,VAR=PRESELECT
*OUTPUT,FIELD,NUMBER INTERVAL=1
*ELEMENT OUTPUT,ELSET=QA_TEST
PEEQ,LE
*END STEP