1.1.11 The Hertz contact problem

Products: ABAQUS/Standard  ABAQUS/Explicit  

The Hertz contact problem (see Timoshenko and Goodier, 1951) provides a classic example for verifying the contact capabilities in ABAQUS. It also serves as an excellent illustration of the use of substructuring in ABAQUS/Standard for locally nonlinear cases (local surface contact). In addition, the problem is analyzed under dynamic conditions in ABAQUS/Standard to illustrate the use of contact surfaces in such cases.

The Hertz contact problem studied consists of two identical, infinitely long cylinders pressed into each other. The solution quantities of most interest are the pressure distribution on the contacting area, the size of the contact area, and the stresses near the contact area. The material behavior is assumed to be linear elastic, and geometric nonlinearities are ignored. Therefore, the only nonlinearity in the problem is the contact constraint.

Problem description

Contact modeling

Because of symmetry, the contact problem can be modeled as a deformable cylinder being pressed against a flat, rigid surface. Therefore, two contact surfaces are required: one (the slave surface in ABAQUS/Standard) on the deformable cylinder and the other (the master surface in ABAQUS/Standard) on the rigid body.

For illustrative purposes several different techniques are used to define the contacting surface pairs. The slave surface is defined by (1) grouping the free faces of elements in an element set that includes all elements in the region that potentially will come into contact (ABAQUS defines the faces automatically), (2) specifying the faces of the elements (or the element sets) in the contact region, or (3) identifying the nodes on the deformable body in the contact region that may come into contact. The master surface is defined by (1) specifying the faces of the rigid elements (or element sets) used to define the rigid body or (2) defining the rigid surface with the *SURFACE option in conjunction with the *RIGID BODY option. Any combination of these techniques can be used together.

By default, ABAQUS uses a finite-sliding contact formulation for modeling the interaction between contact pairs. The contacting surfaces undergo negligible sliding relative to each other, which makes this problem a candidate for the small-sliding contact option. The small-sliding option is invoked by including the SMALL SLIDING parameter on the *CONTACT PAIR option. For a discussion of small- versus finite-sliding contact, see Contact formulation for ABAQUS/Standard contact pairs, Section 29.2.2 of the ABAQUS Analysis User's Manual, or Contact formulation for ABAQUS/Explicit contact pairs, Section 29.4.4 of the ABAQUS Analysis User's Manual.

The surface contact formulation in ABAQUS/Standard gives an accurate solution for the contact area and pressure distribution between the surfaces because of the choice of integration scheme used. Irons and Ahmad (1980) suggest a Gaussian integration rule for calculating self-consistent areas for surface boundary condition problems, which for second-order elements can lead to oscillating results for the pressure distribution on the surface. Oden and Kikuchi explain why this behavior occurs (1980) and present the remedy of using Simpson's integration rule instead. This technique is used in ABAQUS/Standard, and no oscillations in the pressure distribution are found.

The default contact pair formulation in the normal direction in ABAQUS/Standard is hard contact, which gives strict enforcement of contact constraints. Some standard analyses of this problem are conducted with both hard and augmented Lagrangian contact to demonstrate that the default penalty stiffness chosen by the code does not affect stress results significantly. The augmented Lagrangian method is invoked by specifying the AUGMENTED LAGRANGE parameter on the *SURFACE BEHAVIOR option. The hard and augmented Lagrangian contact algorithms are described in Constraint enforcement methods for ABAQUS/Standard contact pairs, Section 29.2.3 of the ABAQUS Analysis User's Manual.

The default contact pair formulation in ABAQUS/Explicit is kinematic contact, which gives strict enforcement of contact constraints. (Note: the small-sliding contact option mentioned previously is available only with kinematic contact.) The explicit dynamic analyses of this problem are conducted with both kinematic and penalty contact to demonstrate that the penetration characteristic of the penalty method can affect stress results significantly in problems with displacement-controlled loading and purely elastic response. The penalty method is invoked by specifying the MECHANICAL CONSTRAINT=PENALTY parameter on the *CONTACT PAIR option. The kinematic and penalty contact algorithms are described in Contact formulation for ABAQUS/Explicit contact pairs, Section 29.4.4 of the ABAQUS Analysis User's Manual.

Substructure ABAQUS/Standard model

This type of contact problem is very suitable for analysis using the substructuring technique in ABAQUS/Standard, since the only nonlinearity in the problem is the contact condition, which is quite local. The cylinder can be defined as a substructure and, thus, reduced to a small number of retained degrees of freedom on the surface where contact may occur or where boundary conditions may be changed. During the iterative solution for contact only these external degrees of freedom on the substructure appear in the equations, thus substantially reducing the cost per iteration. Once the local nonlinearity has been resolved, the solution in the cylinder is recovered as a purely linear response to the known displacements at these retained degrees of freedom. This technique is particularly effective in this case because the rigid surface is flat and there is no friction on the surface; therefore, only the displacement component normal to the surface needs to be retained in the nonlinear iterations.

All information that is relevant to the substructure generation must be given within the *SUBSTRUCTURE GENERATE step, including the degrees of freedom that will be retained in the *RETAINED NODAL DOFS option. The substructure creation and usage cannot be included in the same input file. Only one substructure can be generated per input file. Any number of unit load cases can be defined for the substructure by using the *SUBSTRUCTURE LOAD CASE option. Although this feature is not necessary in this example, it is used in one of the input files for verification purposes.

Substructures are introduced into an analysis model by the *ELEMENT option, where the element number and nodes are defined for each usage of each substructure. Node and element numbers within a substructure and at the usage level are independent—the same node and element numbers can be reused in different substructures and on the usage level. It is also possible to refer to a substructure several times if the structure has identical sections. Thus, once a substructure has been created, it is used just as a standard element type.

Results and discussion

Dynamic analysis in ABAQUS/Standard

A simple dynamic example is created in ABAQUS/Standard by giving the cylinder a uniform initial velocity with the contact conditions all open. This represents the experiment of dropping the cylinder onto a rigid, flat floor under a gravity field.

The impact algorithm used in ABAQUS/Standard for dynamic contact is based on the assumption that, when any contact occurs, the total momentum of the bodies remains unchanged while the points that are contacting will acquire the same velocity instantaneously. In this example the cylinder contacts a rigid surface, which implies that each contacting point will suddenly have zero vertical velocity. This means that a compressive stress wave will emanate from the contacting point and will travel back into the cylinder. After some time this will cause the cylinder to rebound.

It is important to understand that the ABAQUS/Standard dynamic contact algorithm is a “locally perfectly plastic impact” algorithm, as described above, which gives excellent results when it is used correctly. However, it is readily seen that, if the cylinder were modeled as a concentrated mass, with one vertical degree of freedom, the algorithm would imply that the cylinder stops instantaneously when it hits the rigid surface. In reality neither the cylinder nor the surface it hits are rigid: stress waves are started in each. Enough of this detail must be modeled for the results to be meaningful. In this example the cylinder itself is modeled in reasonable detail to capture at least the overall dynamic behavior. If the physical problem from which the example has been developed is that of two cylinders with equal and opposite velocities, this solution is probably useful. If the physical problem is that of a single cylinder hitting a flat surface, it may be necessary to include some elements to model the material below the surface (and the propagation of energy into that domain), unless that material is very dense so that this propagation can be neglected.

Input files

References

Figures

Figure 1.1.11–1 Mesh for the Hertz contact example, ABAQUS/Standard.

Figure 1.1.11–2 Contact pressure (normalized pressure) versus position (normalized distance) for the Hertz contact (no friction) example, ABAQUS/Standard.

Figure 1.1.11–3 Mises stress distribution for the Hertz contact problem, ABAQUS/Standard.

Figure 1.1.11–4 Displaced configuration for the Hertz contact problem, ABAQUS/Standard.

Figure 1.1.11–5 Mesh for the Hertz contact example using CPE4R elements, ABAQUS/Explicit.

Figure 1.1.11–6 Mesh for the Hertz contact example using CPE6M elements, ABAQUS/Explicit.

Figure 1.1.11–7 Contact pressure contour for the Hertz contact problem using CPE4R elements and kinematic contact, ABAQUS/Explicit.

Figure 1.1.11–8 Contact pressure contour for the Hertz contact problem using CPE4R elements and penalty contact, ABAQUS/Explicit.

Figure 1.1.11–9 Mises stress distribution for the Hertz contact problem using CPE4R elements and kinematic contact, ABAQUS/Explicit.

Figure 1.1.11–10 Mises stress distribution for the Hertz contact problem using CPE4R elements and penalty contact, ABAQUS/Explicit.

Figure 1.1.11–11 Displaced configuration for the Hertz contact problem using CPE4R elements and kinematic contact, ABAQUS/Explicit.

Figure 1.1.11–12 Displaced configuration for the Hertz contact problem using CPE4R elements and penalty contact, ABAQUS/Explicit.

Figure 1.1.11–13 Contact pressure contour for the Hertz contact problem using CPE6M elements and kinematic contact, ABAQUS/Explicit.

Figure 1.1.11–14 Mises stress distribution for the Hertz contact problem using CPE6M elements and kinematic contact, ABAQUS/Explicit.

Figure 1.1.11–15 Displaced configuration for the Hertz contact problem using CPE6M elements and kinematic contact, ABAQUS/Explicit.