*HEADING
Thick plate rolling: 3-Dimensional, ABAQUS/Explicit
SECTION CONTROLS USED (KINEMA=CENTROID, HOURGLASS=STIFFNESS)
Add *integrated output to output section force on contact pair surface SURF1.
SOF and SOM should match CFT and CMT on surface SURF1.
*RESTART,WRITE,NUM=4
***********************************************
*NODE
** Bar
1,        0.,     0.
801,      0.,     0.020
47,     -00.092,  0.
847,    -00.092,  0.020
24001,    0.,     0.   , -0.020
24801,    0.,     0.020, -0.020
24047,  -00.092,  0.   , -0.020
24847,  -00.092,  0.020, -0.020
*NGEN,NSET=BOT1   
1,47,1
*NGEN,NSET=TOP1
801,847,1
*NGEN,NSET=BOT2   
24001,24047,1
*NGEN,NSET=TOP2
24801,24847,1
*NFILL,NSET=ZSYMM
BOT1,TOP1,8,100
*NFILL,NSET=SIDE
BOT2,TOP2,8,100
*NFILL,NSET=BAR
ZSYMM,SIDE,8,3000
*NSET,NSET=BOTTOM,GEN
1,47,1
3001,3047,1
6001,6047,1
9001,9047,1
12001,12047,1
15001,15047,1
18001,18047,1
21001,21047,1
24001,24047,1
**
** Roller, Radius = .170 m
**
*NODE
** Reference node
30000,    0.0409 ,  0.185 , 0.010
31000,    0.0409 ,  0.185 ,-0.030
** Nodes forming roll
30001 ,  -0.1291 ,  0.185 , 0.010
30090 ,   0.0409 ,  0.015 , 0.010
32001 ,  -0.1291 ,  0.185 ,-0.030
32090 ,   0.0409 ,  0.015 ,-0.030
**
*NGEN,NSET=ROLL1,LINE=C
30001,30090,1,30000
*NGEN,NSET=ROLL2,LINE=C
32001,32090,1,31000
**
***********************************************
***** Bar
**
*ELEMENT,TYPE=C3D8R,ELSET=METAL
1,2,1,3001,3002,102,101,3101,3102
*ELGEN,ELSET=METAL
1,8,100,100,46,1,1,8,3000,1000
**
*ELSET,ELSET=TOP,GEN
701,746,1
1701,1746,1
2701,2746,1
3701,3746,1
4701,4746,1
5701,5746,1
6701,6746,1
7701,7746,1
*ELSET,ELSET=BACK,GEN
46,746,100
1046,1746,100
2046,2746,100
3046,3746,100
4046,4746,100
5046,5746,100
6046,6746,100
7046,7746,100
*ELSET,ELSET=SIDE,GEN
7001,7046,1
7101,7146,1
7201,7246,1
7301,7346,1
7401,7446,1
7501,7546,1
7601,7646,1
7701,7746,1
**
*SOLID SECTION,ELSET=METAL,MAT=C15, CONTROL=C
1., 
*SECTION CONTROLS,KINEMA=CENTROID,HOURGLASS=STIFFNESS,NAME=C
**
*MATERIAL,NAME=C15
*ELASTIC
1.5E11,.3
*PLASTIC
168.72E06,0
219.33E06,0.1
272.02E06,0.2
308.53E06,0.3
337.37E06,0.4
361.58E06,0.5
382.65E06,0.6
401.42E06,0.7
418.42E06,0.8
434.01E06,0.9
448.45E06,1.0
*DENSITY
7850, 
**
***********************************************
****   Roll
**
*ELEMENT,TYPE=R3D4,ELSET=ROLLER
10001,30001,30002,32002,32001
*ELGEN,ELSET=ROLLER
10001,89,1,1
**
*INITIAL CONDITIONS,TYPE=VELOCITY
BAR,1,1.0367
**
***************** Step 1
*SURFACE,TYPE=ELEMENT, NAME=SURF1
TOP,S2
SIDE,S5
*SURFACE,TYPE=ELEMENT, NAME=SURF2
ROLLER,SNEG
*RIGID BODY,ELSET=ROLLER,REF NODE=30000
*STEP
*DYNAMIC,EXPLICIT
,0.089286
*FIXED MASS SCALING,FACTOR=2758.5,ELSET=METAL
**
*BOUNDARY
BOTTOM,2,2
ZSYMM,ZSYMM
30000,1,5
**
*BOUNDARY,TYPE=VELOCITY
**30000,6,6,330.0
30000,6,6,6.2832
**
*SURFACE INTERACTION,NAME=FRICT
*FRICTION
 0.3,
*CONTACT PAIR,INTERACTION=FRICT,CPSET=TEST
SURF1,SURF2
**
**
*FILE OUTPUT,TIMEMARKS=YES,NUM=1
*EL FILE, ELSET=QA_TEST
PEEQ,
*NODE FILE, NSET=QA_TEST
U, 
*ENERGY FILE
*****
*ELSET, ELSET=QA_TEST
SIDE,
*NSET, NSET=QA_TEST
SIDE,
*OUTPUT, FIELD, TIME MARKS=YES, NUMBER INTERVAL=1
*ELEMENT OUTPUT, ELSET=QA_TEST
PEEQ, 
*NODE OUTPUT, NSET=QA_TEST
U, 
*OUTPUT, FIELD, TIME MARKS=YES, NUMBER INTERVAL=1
*ELEMENT OUTPUT
PEEQ,
*NODE OUTPUT
U,
*OUTPUT, HIST, TIME INT=0.02
*ENERGY OUTPUT, VAR=PRESELECT 
*INTEGRATED OUTPUT, SURFACE=SURF1
SOF,SOM,
*CONTACT OUTPUT, CPSET=TEST
CFT,CMT
****
******
*END STEP
**********************************************