*HEADING
HERTZ CONTACT PROBLEM 3D:
ELEM: C3D20R, C3D20, R3D4    INTERFACE: SMALL SLIDING
*RESTART,WRITE,OVERLAY,FREQUENCY=10
*NODE
*NODE,NSET=RIGNODE
*NODE,SYSTEM=C
*NGEN,LINE=C
*NGEN
*NSET,NSET=NX,GENERATE
*ELEMENT,TYPE=C3D20R, ELSET=CYL
*ELEMENT,TYPE=C3D20,ELSET=CYL
*ELSET,ELSET=TOP
*ELSET,ELSET=BOT
*EQUATION
*MPC
QUADRATIC,82,81,83,85
QUADRATIC,84,81,83,85
C BIQUAD,1083,81,85,2085,2081,83,1085,2083,1081
QUADRATIC,2082,2081,2083,2085
QUADRATIC,2084,2081,2083,2085
QUADRATIC,122,121,123,125
QUADRATIC,124,121,123,125
C BIQUAD,1123,121,125,2125,2121,123,1125,2123,1121
QUADRATIC,2122,2121,2123,2125
QUADRATIC,2124,2121,2123,2125
QUADRATIC,126,125,127,129
QUADRATIC,128,125,127,129
C BIQUAD,1127,125,129,2129,2125,127,1129,2127,1125
QUADRATIC,2126,2125,2127,2129
QUADRATIC,2128,2125,2127,2129
QUADRATIC,130,129,131,133
QUADRATIC,132,129,131,133
C BIQUAD,1131,129,133,2133,2129,131,1133,2131,1129
QUADRATIC,2130,2129,2131,2133
QUADRATIC,2132,2129,2131,2133
QUADRATIC,189,209,169,129
QUADRATIC,149,209,169,129
QUADRATIC,2189,2209,2169,2129
QUADRATIC,2149,2209,2169,2129
C BIQUAD,1169,129,209,2209,2129,169,1209,2169,1129
*ELSET,ELSET=REF1
*ELSET,ELSET=REF2
**
** Contact definition
**
*ELEMENT,TYPE=R3D4,ELSET=ERIGID
*RIGID BODY,ELSET=ERIGID,REF NODE=10000
*surface,NAME=RIGID
ERIGID,SPOS
*surface,NAME=CYLBOT
BOT,S4
*CONTACT PAIR, INTERACTION=HERTZ,SMALL SLIDING, TYPE=SURFACE TO SURFACE
CYLBOT,RIGID
*SURFACE INTERACTION,NAME=HERTZ
*SURFACE BEHAVIOR, AUGMENTED LAGRANGE
*NSET,NSET=ZSYM,GENERATE
*BOUNDARY
ZSYM,3,3
NX,1,1
10000, 1,6
*SOLID SECTION,MATERIAL=A1, ELSET=CYL
*MATERIAL,NAME=A1
*ELASTIC
30.E6,0.3
*STEP,INC=10
*STATIC,DIRECT
.1,1.
*BOUNDARY
1,2,,-0.4
*PRINT,CONTACT=YES
*EL PRINT,FREQUENCY=2
S11,E11
*NODE PRINT,FREQUENCY=10
U,
RF,
*OUTPUT,FIELD,FREQ=10
*NODE OUTPUT
RF,
*OUTPUT,FIELD,FREQ=5
*CONTACT OUTPUT,VARIABLE=PRESELECT,MASTER=RIGID,SLAVE=CYLBOT
*END STEP