*HEADING
 HERTZ CONTACT PROBLEM 3D:             
 ELEM: C3D20R, C3D20, R3D4    INTERFACE: SMALL SLIDING
*RESTART,WRITE,OVERLAY,FREQUENCY=10
*NODE
   1, 0.,0.
   5, 2.,0.
   9, 5.5,0.
  13, 10.,0.
  41, 0.,-3.
  45, 2.,-3.
  49, 5.2,-3.
  81, 0.,-5.5
  85, 2.,-5.5
  89, 4.,-5.5
 121, 0.,-7.5
 129, 2.,-7.5
 133, 3.5,-7.5
 101, 0.,-6.5
 103, 1.,-6.5
 105, 2.,-6.5
 107, 3.75,-6.5
 201, 0.,-10.
2001,.0,.0,1.
1001,.0,.0,.5
2005,2.,.0,1.
1005,2.,.0,.5
2009,5.5,0.,1.
1009,5.5,0.,.5
2013,10.,0.,1.
1013,10.,0.,.5
2041,0.,-3.,1.
1041,0.,-3.,.5
2045,2.,-3.,1.
1045,2.,-3.,.5
2049,5.2,-3.,1.
1049,5.2,-3.,.5
2081,.0,-5.5,1.
1081,.0,-5.5,.5
2085,2.,-5.5,1.
1085,2.,-5.5,.5
2089,4.,-5.5,1.
1089,4.,-5.5,.5
2121,0.,-7.5,1.
1121,0.,-7.5,.5
2129,2.,-7.5,1.
1129,2.,-7.5,.5
2133,3.5,-7.5,1.
1133,3.5,-7.5,.5
2101,0.,-6.5,1.
1101,0.,-6.5,.5
2103,1.,-6.5,1.
1103,1.,-6.5,.5
2105,2.,-6.5,1.
1105,2.,-6.5,.5
2107,3.75,-6.5,1.
1107,3.75,-6.5,.5
2201,.0,-10.,1.
1201,.0,-10.,.5
*NODE,NSET=RIGNODE
5000,-.5,-10,-.5
5001,-.5,-10,1.5
5002,2.5,-10,-.5
5003,2.5,-10,1.5
10000,1.,-11,1
*NODE,SYSTEM=C
 33,10.,-15.
 53,10.,-30.
 73,10.,-37.5
 93,10.,-45.
113,10.,-55.
213,10.,-65.
209,10.,-78.
1033,10.,-15.,.5
2033,10.,-15.,1.
1053,10.,-30.,.5
2053,10.,-30.,1.
1073,10.,-37.5,.5
2073,10.,-37.5,1.
1093,10.,-45.,.5
2093,10.,-45.,1.
1113,10.,-55.,.5
2113,10.,-55.,1.
1213,10.,-65.,.5
2213,10.,-65.,1.
1209,10.,-78.,.5
2209,10.,-78.,1.
*NGEN,LINE=C
 209,213
 201,209
1209,1213,1,1001
2209,2213,1,2001
1201,1209,1,1001
2201,2209,1,2001
*NGEN
 1,5,2
 5,9,2
 9,13,2
 41,45,2
 45,49,2
 49,53,2
 81,85
 85,89,2
 89,93,2
 121,129
 129,133
 133,213,40
 131,211,40
 1,41,20
 5,45,20
 9,49,20
 41,81,20
 45,85,20
 49,89,20
 81,121,20
 121,201,20
 123,203,20
 125,205,20
 127,207,20
 129,209,20
 161,163
 163,165
 165,167
 167,169
1001,1005,2
2001,2005,2
1005,1009,2
2005,2009,2
1009,1013,2
2009,2013,2
1041,1045,2
2041,2045,2
1045,1049,2
2045,2049,2
1049,1053,2
2049,2053,2
1081,1085
2081,2085
1085,1089,2
2085,2089,2
1089,1093,2
2089,2093,2
1121,1129
2121,2129
1129,1133
2129,2133
1133,1213,40
2133,2213,40
1131,1211,40
2131,2211,40
1001,1041,20
2001,2041,20
1005,1045,20
2005,2045,20
1009,1049,20
2009,2049,20
1041,1081,20
2041,2081,20
1045,1085,20
2045,2085,20
1049,1089,20
2049,2089,20
1081,1121,20
2081,2121,20
1121,1201,20
2121,2201,20
1123,1203,20
2123,2203,20
1125,1205,20
2125,2205,20
1127,1207,20
2127,2207,20
1129,1209,20
2129,2209,20
1161,1163
2161,2163
1163,1165
2163,2165
1165,1167
2165,2167
1167,1169
2167,2169
*NSET,NSET=NX,GENERATE
1,201,20
1001,1201,40
2001,2201,20
*ELEMENT,TYPE=C3D20R, ELSET=CYL
1,1,41,45,5,2001,2041,2045,2005,21,43,25,3,2021,2043,2025,
2003,1001,1041,1045,1005
2,5,45,49,9,2005,2045,2049,2009,25,47,29,7,2025,2047,2029,
2007,1005,1045,1049,1009
3,9,49,53,13,2009,2049,2053,2013,29,51,33,11,2029,2051,2033,
2011,1009,1049,1053,1013
4,41,81,85,45,2041,2081,2085,2045,61,83,65,43,2061,2083,2065,
2043,1041,1081,1085,1045
5,45,85,89,49,2045,2085,2089,2049,65,87,69,47,2065,2087,2069,
2047,1045,1085,1089,1049
6,49,89,93,53,2049,2089,2093,2053,69,91,73,51,2069,2091,2073,
2051,1049,1089,1093,1053
7,81,121,125,83,2081,2121,2125,2083,101,123,103,82,2101,2123,2103,
2082,1081,1121,1125,1083
8,83,125,129,85,2083,2125,2129,2085,103,127,105,84,2103,2127,2105,
2084,1083,1125,1129,1085
9,85,129,133,89,2085,2129,2133,2089,105,131,107,87,2105,2131,2107,
2087,1085,1129,1133,1089
10,89,133,213,93,2089,2133,2213,2093,107,173,113,91,2107,2173,2113,
2091,1089,1133,1213,1093
11,121,161,163,123,2121,2161,2163,2123,141,162,143,122,2141,2162,2143,
2122,1121,1161,1163,1123
12,123,163,165,125,2123,2163,2165,2125,143,164,145,124,2143,2164,2145,
2124,1123,1163,1165,1125
13,125,165,167,127,2125,2165,2167,2127,145,166,147,126,2145,2166,2147,
2126,1125,1165,1167,1127
14,127,167,169,129,2127,2167,2169,2129,147,168,149,128,2147,2168,2149,
2128,1127,1167,1169,1129
19,129,209,211,131,2129,2209,2211,2131,169,210,171,130,2169,2210,2171,
2130,1129,1209,1211,1131
20,131,211,213,133,2131,2211,2213,2133,171,212,173,132,2171,2212,2173,
2132,1131,1211,1213,1133
*ELEMENT,TYPE=C3D20,ELSET=CYL
15,161,201,203,163,2161,2201,2203,2163,181,202,183,162,2181,2202,2183,
2162,1161,1201,1203,1163
16,163,203,205,165,2163,2203,2205,2165,183,204,185,164,2183,2204,2185,
2164,1163,1203,1205,1165
17,165,205,207,167,2165,2205,2207,2167,185,206,187,166,2185,2206,2187,
2166,1165,1205,1207,1167
18,167,207,209,169,2167,2207,2209,2169,187,208,189,168,2187,2208,2189,
2168,1167,1207,1209,1169
*ELSET,ELSET=TOP
 1,2,3,4,5,6,7,8,9,10,11,12,13,14,19,20
*ELSET,ELSET=BOT
 15,16,17,18
*EQUATION
 2, 
 3,2,-1.,1,2,1.
 2, 
 5,2,-1.,1,2,1.
 2, 
 7,2,-1.,1,2,1.
 2, 
 9,2,-1.,1,2,1.
 2, 
11,2,-1.,1,2,1.
 2, 
13,2,-1.,1,2,1.
2, 
2001,2,-1.,1,2,1.
2, 
2003,2,-1.,1,2,1.
2, 
2005,2,-1.,1,2,1.
2, 
2007,2,-1.,1,2,1.
2, 
2009,2,-1.,1,2,1.
2, 
2011,2,-1.,1,2,1.
2, 
2013,2,-1.,1,2,1.
2, 
1001,2,-1.,1,2,1.
2, 
1005,2,-1.,1,2,1.
2, 
1009,2,-1.,1,2,1.
2, 
1013,2,-1.,1,2,1.
*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
11,12,13,14,19,20
*ELSET,ELSET=REF2
7,8,9
**
** Contact definition
**
*ELEMENT,TYPE=R3D4,ELSET=ERIGID
100,5000,5001,5003,5002
*RIGID BODY,ELSET=ERIGID,REF NODE=10000
*surface,NAME=RIGID
ERIGID,SPOS
*surface,NAME=CYLBOT
BOT,S4
*CONTACT PAIR, INTERACTION=HERTZ,SMALL SLIDING
CYLBOT,RIGID
*SURFACE INTERACTION,NAME=HERTZ
*SURFACE BEHAVIOR, AUGMENTED LAGRANGE
*NSET,NSET=ZSYM,GENERATE
1,300,1
2001,2300,1
*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,MAIN=RIGID,SECONDARY=CYLBOT
*END STEP