*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, TYPE=SURFACE TO SURFACE
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,MASTER=RIGID,SLAVE=CYLBOT
*END STEP