7.1 Hertzian contact theory

Hertzian contact theory explains how curved surfaces interact under load, predicting the stress distributions and deformations that result. It's foundational to friction and wear analysis because nearly every mechanical system with moving parts involves curved surfaces pressing together, from bearings and gears to wheel-rail interfaces.

Heinrich Hertz developed this theory in 1882 while studying Newton's rings (optical interference patterns between curved glass surfaces). The theory assumes smooth, frictionless surfaces undergoing only elastic deformation. Those assumptions limit its direct applicability, but Hertzian theory remains the starting point for virtually all contact mechanics analysis.

Fundamentals of Hertzian contact

Hertzian contact theory predicts what happens when two curved bodies are pressed together: how large the contact area becomes, how pressure distributes across it, and how the materials deform. These predictions feed directly into wear rate estimates, fatigue life calculations, and friction analysis for mechanical components.

Elastic deformation principles

Elastic deformation is the temporary shape change a material undergoes when loaded. Remove the load, and the material returns to its original shape. Hertzian theory deals exclusively with this regime.

The governing relationship is Hooke's law: stress is proportional to strain within the elastic limit. Two material properties control elastic behavior:

• Elastic modulus (Young's modulus, $E$): How stiff the material is. Higher $E$ means less deformation for a given stress.
• Poisson's ratio ($\nu$): How much a material contracts laterally when stretched axially. Most metals fall between 0.25 and 0.35.

All the energy stored during elastic deformation is recoverable when the load is removed. This is a key distinction from plastic deformation, where energy is dissipated permanently.

Assumptions and limitations

Hertzian theory relies on several simplifying assumptions. You need to know these because they tell you when the theory applies and when it doesn't:

• Surfaces are perfectly smooth (no roughness)
• No friction or adhesion between contacting bodies
• Only elastic deformation occurs (stresses stay below yield strength)
• Deformations are small compared to the overall dimensions of the bodies
• Materials are isotropic (same properties in all directions) and homogeneous
• No time-dependent behavior (no creep or viscoelasticity)

When any of these assumptions is significantly violated, you need to move beyond classical Hertzian theory. Real surfaces are rough, friction is always present, and plastic deformation often occurs at high loads. The "Advanced topics" section at the end covers some of these extensions.

Historical development

• Hertz (1882): Developed the original theory for contact between two spheres or a sphere and a plane.
• Boussinesq (1885): Contributed solutions for point loads on elastic half-spaces, providing mathematical tools that underpin contact mechanics.
• Sneddon (1965): Extended the theory for axisymmetric indentation problems (e.g., flat punch, cone, and paraboloid indenters).
• Later researchers generalized the theory to cylinders, ellipsoids, and other geometries.
• Modern approaches use finite element analysis (FEA) and computational methods to handle complex geometries and non-linear behavior that analytical solutions can't address.

Contact geometry

The shape of the contact between two bodies determines how stress distributes and how wear progresses. Hertzian theory distinguishes three fundamental contact types, each with different pressure profiles and engineering implications.

Point contact

Point contact occurs when two curved surfaces touch at what is initially a single point. Under load, this point expands into a small circular or elliptical area.

• Typical examples: ball bearings, spherical joints, ball-on-flat test configurations
• Contact pressures are high and concentrated at the center
• Deformation is highly localized
• Common in precision instruments and rolling element bearings where loads are transmitted through small regions

Line contact

Line contact results when two cylindrical surfaces meet along a line. Under load, this line expands into a narrow rectangular strip.

• Typical examples: roller bearings, gear teeth, cam-follower systems
• Stress distributes more uniformly along the contact length compared to point contact
• For the same total load, contact pressures are lower than in point contact because the load spreads over a longer region
• Preferred in applications where better load distribution and reduced wear are needed

Area contact

Area contact occurs between conforming surfaces that already share a region of contact even before load is applied.

• Typical examples: flat bearings, seals, gaskets, journal bearings
• Pressure distribution is more uniform than point or line contact
• Contact pressures are the lowest of the three types for a given load
• Provides the best stability and load-bearing capacity, but Hertzian theory is least applicable here because the "small contact area" assumption breaks down

Stress distribution

Understanding how stresses distribute through and beneath the contact area is critical for predicting fatigue, wear, and failure. Three stress measures are particularly important.

Normal stress

Normal stress acts perpendicular to the contact surface. In Hertzian contact, it's highest at the center and drops to zero at the edges, following a semi-ellipsoidal profile.

For a circular contact area, the normal stress at the surface is:

$\sigma_z = -p_0 \sqrt{1 - (r/a)^2}$

where $p_0$ is the maximum contact pressure, $r$ is the radial distance from the center, and $a$ is the contact radius. The negative sign indicates compression.

This stress profile determines overall deformation and is the first thing to check when evaluating whether a contact will remain elastic or begin to yield.

Shear stress

Shear stress acts parallel to the contact surface and is responsible for much of the damage in contacting components.

• The maximum shear stress occurs below the surface, not at it. For a frictionless circular contact, this maximum is roughly at a depth of $0.48a$.
• Subsurface shear stress drives crack initiation and propagation, which is why bearing failures often start beneath the surface.
• Repeated shear stress cycling causes fatigue damage over time, even when individual load cycles are well within elastic limits.

von Mises stress

The von Mises stress combines all normal and shear stress components into a single scalar value that predicts yielding in ductile materials:

$\sigma_{vm} = \sqrt{\frac{1}{2}[(\sigma_1 - \sigma_2)^2 + (\sigma_2 - \sigma_3)^2 + (\sigma_3 - \sigma_1)^2]}$

where $\sigma_1$, $\sigma_2$, and $\sigma_3$ are the principal stresses. If $\sigma_{vm}$ exceeds the material's yield strength, plastic deformation begins. This criterion is used to determine safety factors and to guide material selection for contacting components.

Contact pressure

Contact pressure is the force per unit area at the interface between two bodies. Its magnitude and distribution directly affect wear rates, friction behavior, and fatigue life.

Maximum contact pressure

For a circular (sphere-on-plane) contact, the maximum pressure occurs at the center and is given by:

$p_0 = \frac{3F}{2\pi a^2}$

where $F$ is the applied normal load and $a$ is the contact radius. This is 1.5 times the mean pressure over the contact area.

Maximum contact pressure is the critical parameter for checking whether the contact remains elastic. It depends on load, material stiffness, and the radii of curvature of the contacting bodies. Smaller radii and stiffer materials produce higher peak pressures.

Pressure distribution

The pressure across a circular contact follows a semi-ellipsoidal profile:

$p(r) = p_0 \sqrt{1 - (r/a)^2}$

This means pressure is highest at the center and smoothly drops to zero at the edge of the contact circle. For conforming area contacts, the distribution is more uniform. Surface treatments or coatings can modify the effective stiffness and alter the pressure profile.

Effect of load

The relationship between applied load and contact geometry is nonlinear. For circular point contact, the contact radius is:

$a = \sqrt[3]{\frac{3FR}{4E^*}}$

where $R$ is the equivalent radius of curvature and $E^*$ is the equivalent elastic modulus.

Because of the cube root, the contact area grows as $F^{2/3}$. Doubling the load does not double the contact area; it increases it by a factor of about 1.59. This nonlinearity is a defining feature of Hertzian contact and means that contact pressure increases more slowly than you might expect as load rises.

Contact area

The contact area is the region where two bodies physically touch. Its size and shape govern how load, pressure, and stress are distributed.

Circular vs. elliptical contact

• Circular contact occurs between identical spheres, or a sphere and a flat plane. The contact area is $A = \pi a^2$.
• Elliptical contact results when the contacting bodies have different curvatures in two perpendicular directions, such as a ball in a grooved raceway. The contact area is $A = \pi ab$, where $a$ and $b$ are the semi-major and semi-minor axes.

Elliptical contacts are extremely common in practice. Ball bearings, for instance, produce elliptical contact patches because the ball's curvature differs from the raceway's curvature in the axial and circumferential directions.

Calculation methods

For simple geometries, closed-form Hertzian equations give exact results. The contact radius for sphere-on-plane is:

$a = \sqrt[3]{\frac{3FR}{4E^*}}$

For elliptical contacts, the calculation requires solving transcendental equations involving elliptic integrals, or using published approximation formulas and tabulated coefficients. For complex or non-standard geometries, finite element analysis (FEA) is the standard approach and can also account for surface roughness and non-Hertzian effects.

Influence of material properties

The equivalent elastic modulus $E^*$ captures the combined stiffness of both contacting materials:

$\frac{1}{E^*} = \frac{1-\nu_1^2}{E_1} + \frac{1-\nu_2^2}{E_2}$

Softer materials (lower $E$) produce larger contact areas and lower peak pressures for the same load. When two materials with very different stiffnesses are in contact, the softer one dominates the deformation. Surface coatings can change the effective properties at the contact, which is one reason hard coatings (like TiN on steel) are used to reduce wear.

Deformation analysis

Deformation analysis examines how the shapes of contacting bodies change under load. This matters for predicting conformity changes, wear evolution, and failure modes.

Surface displacement

Surface displacement describes how much each point on the surface moves vertically under the contact pressure. Maximum displacement occurs at the center of the contact. Outside the contact area, the surface still deflects, but the displacement decreases with distance.

The displacement field for a circular contact involves complex expressions from elasticity theory. In practice, these displacements are small (micrometers for typical metal contacts), but they determine the real contact area and pressure distribution. Experimental measurement techniques include optical interferometry and surface profilometry.

Subsurface deformation

The stresses transmitted through the contact don't just affect the surface. They penetrate into the bulk material, and the maximum shear stress typically occurs at a depth below the surface (around $0.48a$ for frictionless contact).

This is why contact fatigue cracks often initiate subsurface and propagate toward the surface, eventually causing pitting or spalling. The depth of maximum shear stress shifts depending on the friction coefficient at the surface: friction pulls the maximum shear stress location closer to the surface.

Elastic vs. plastic deformation

The transition from elastic to plastic deformation occurs when the von Mises stress first exceeds the material's yield strength. This typically happens subsurface, at the location of maximum shear stress.

• At low loads, deformation is fully elastic and reversible.
• As load increases, a plastic zone forms beneath the surface while the surrounding material remains elastic (contained plasticity).
• At very high loads, the plastic zone reaches the surface and the contact becomes fully plastic.
• Repeated elastic-plastic cycling causes fatigue damage, residual stress buildup, and progressive changes in surface topography.

Analyzing plastic deformation requires nonlinear FEA because the closed-form Hertzian equations only cover the elastic regime.

Hertzian equations

The Hertzian equations are the mathematical core of the theory. They provide closed-form solutions for contact area, pressure, and deformation for idealized geometries.

Derivation of key equations

The derivation starts from classical elasticity theory and proceeds through several steps:

1. Model each contacting body as an elastic half-space (valid when the contact area is small relative to body dimensions).
2. Express the surface displacement due to a distributed pressure using Boussinesq's solution for a point load on a half-space.
3. Apply the geometric constraint that the gap between the undeformed surfaces must equal the combined surface displacement within the contact area.
4. Apply boundary conditions: pressure is zero outside the contact area, and there's no interpenetration of surfaces.
5. Solve the resulting integral equation to obtain the pressure distribution, contact radius, and displacement field.

The mathematics involves potential theory, integral transforms, and superposition. The result is the semi-ellipsoidal pressure distribution and the cube-root load-displacement relationships that define Hertzian contact.

Application to different geometries

Sphere-on-plane (circular contact):

• Pressure distribution: $p(r) = p_0 \sqrt{1 - (r/a)^2}$
• Contact radius: $a = \sqrt[3]{\frac{3FR}{4E^*}}$
• Maximum pressure: $p_0 = \frac{3F}{2\pi a^2}$

Cylinder-on-plane (line contact):

• Pressure distribution: $p(x) = p_0 \sqrt{1 - (x/a)^2}$
• Here $a$ is the half-width of the contact strip, and the pressure profile has the same semi-elliptical shape but extends uniformly along the cylinder length.

Elliptical contact (e.g., crossed cylinders, ball in grooved raceway) requires numerical solutions or tabulated approximation coefficients because the equations involve elliptic integrals.

Limitations and extensions

• Valid only for small elastic deformations of homogeneous, isotropic materials
• Does not account for friction, adhesion, or surface roughness
• Extensions for viscoelastic materials allow analysis of polymers and biological tissues
• Layered structure models handle coated or surface-treated components
• Non-Hertzian theories (e.g., Greenwood-Williamson for rough surfaces, JKR/DMT for adhesion) address scenarios where classical assumptions fail
• FEA handles arbitrary geometries, nonlinear materials, and coupled thermal-mechanical effects

Material considerations

Material properties directly control contact area size, pressure magnitude, and deformation behavior. Selecting the right materials for contacting surfaces is one of the most impactful design decisions in tribology.

Elastic modulus

The elastic modulus ($E$) measures resistance to elastic deformation. Steel has $E \approx 200$ GPa, aluminum about 70 GPa, and polymers can be below 5 GPa.

• Higher $E$ means smaller contact areas and higher contact pressures for a given load.
• The modulus controls the overall stiffness and load-bearing capacity of the contact.
• For anisotropic materials like composites, the modulus depends on direction, complicating the analysis.
• $E$ generally decreases with increasing temperature, which matters for high-temperature applications.

Poisson's ratio

Poisson's ratio ($\nu$) describes how much a material contracts laterally when stretched. It typically ranges from about 0.1 (cork) to nearly 0.5 (rubber, which is essentially incompressible).

Most metals fall between 0.25 and 0.35. Poisson's ratio affects the lateral expansion under compressive contact loads and influences the subsurface stress distribution. It appears in the equivalent modulus formula, so it does affect contact area and pressure calculations, though its influence is smaller than that of $E$.

Influence on contact behavior

The equivalent elastic modulus $E^*$ combines the properties of both contacting materials into a single parameter:

$\frac{1}{E^*} = \frac{1-\nu_1^2}{E_1} + \frac{1-\nu_2^2}{E_2}$

This means the softer material dominates the contact behavior. A steel ball on a polymer surface behaves very differently from a steel ball on a steel surface, even under the same load. The polymer contact will have a much larger area and lower peak pressure.

Beyond elastic properties, the yield strength determines when plastic deformation begins, hardness correlates with wear resistance, and toughness affects fatigue crack resistance. Surface treatments like carburizing, nitriding, or hard coatings can improve contact performance without changing the bulk material.

Applications in engineering

Bearings and gears

Rolling element bearings are the most direct application of Hertzian theory. Each ball or roller creates a Hertzian contact patch on the raceway, and the contact stress at these patches determines bearing life.

• Bearing life prediction (e.g., the Lundberg-Palmgren model) relies on Hertzian contact stress calculations.
• Gear tooth contact is modeled as line contact (spur gears) or elliptical contact (helical gears) to optimize tooth profiles and predict pitting fatigue.
• Lubrication film thickness in elastohydrodynamic lubrication (EHL) depends on the Hertzian contact pressure and deformation.
• Common failure modes like pitting, spalling, and micropitting are all driven by Hertzian contact stresses.

Wheel-rail contact

Railway engineering relies heavily on contact mechanics. The wheel-rail contact patch is typically elliptical, roughly 1-2 cm across, and transmits enormous loads.

• Hertzian theory predicts the contact patch size, shape, and pressure distribution.
• Wear patterns on wheels and rails are analyzed using contact stress predictions.
• Rolling contact fatigue causes rail defects like head checks and squats, driven by repeated Hertzian stress cycling.
• For more accurate predictions in cases with significant creep or conformal contact, non-Hertzian methods are used.

MEMS devices

At the microscale, Hertzian contact mechanics still applies but with additional considerations.

• Microswitches and microactuators involve tiny contact areas where Hertzian pressures can be extremely high.
• Adhesion forces (van der Waals, capillary) become significant relative to applied loads at small scales, often requiring JKR or DMT models rather than pure Hertzian theory.
• Surface roughness effects are proportionally larger at the microscale.
• Stiction (static friction preventing release) is a major reliability concern in MEMS, and contact mechanics analysis helps mitigate it.

Experimental validation

Theoretical predictions need experimental verification, especially when pushing beyond the idealized Hertzian assumptions.

Photoelasticity techniques

Photoelasticity exploits the fact that certain transparent materials become birefringent (they split light into two polarized beams) when stressed. By loading a transparent model and viewing it through polarized light, you can directly see the stress field as colored fringe patterns.

• Provides a full-field visualization of stress distributions, not just point measurements
• Allows direct observation of where maximum shear stresses occur
• Useful for validating stress predictions in complex geometries
• Limited to transparent materials with suitable photoelastic properties (e.g., epoxy resins, polycarbonate)

Ultrasonic methods

Ultrasonic techniques use high-frequency sound waves to probe the contact interface in opaque, real-world components.

• When an ultrasonic wave hits a contact interface, the proportion reflected depends on the contact pressure at that point.
• By scanning across the interface, you can map the pressure distribution and real contact area.
• This is a nondestructive technique that works on actual engineering components during operation.
• Requires careful calibration to convert reflection coefficients into pressure values.

Numerical simulations

Finite Element Analysis (FEA) is the standard computational tool for contact problems that go beyond simple Hertzian cases.

• Can model complex geometries, nonlinear materials, friction, adhesion, and surface roughness
• Provides complete stress, strain, and displacement fields throughout both contacting bodies
• Used to optimize designs and predict failure modes before building prototypes
• Results should always be validated against analytical solutions (for simple cases) or experimental data to ensure accuracy

Advanced topics

Classical Hertzian theory covers idealized cases. Real engineering contacts involve complications that require extended theories.

Non-Hertzian contact

Non-Hertzian contact analysis is needed when one or more Hertzian assumptions are significantly violated:

• Conformal contacts where the contact area is large relative to body dimensions (e.g., journal bearings)
• Large deformations that violate the small-strain assumption
• Edge effects in finite-sized bodies where the half-space assumption breaks down
• Layered or coated materials where properties vary with depth

These problems generally require numerical methods, though some analytical solutions exist for specific cases.

Rough surface effects

Real surfaces are never perfectly smooth. Surface roughness has major implications for contact mechanics:

• The real contact area (where asperities actually touch) is typically a small fraction of the apparent contact area. For metal surfaces under moderate loads, real contact area might be only 1-10% of apparent area.
• Contact occurs at discrete asperity tips, creating local pressures much higher than the nominal Hertzian prediction.
• The Greenwood-Williamson model treats rough surface contact statistically, modeling asperities as a distribution of spherical contacts.
• Roughness strongly influences friction, wear, lubrication effectiveness, and electrical/thermal contact resistance.

Adhesion and friction influences

At small scales or with soft materials, adhesion forces become comparable to applied loads and can no longer be ignored.

• The JKR (Johnson-Kendall-Roberts) model applies when adhesion is strong and materials are compliant. It predicts a contact area larger than Hertzian theory and a finite pull-off force.
• The DMT (Derjaguin-Muller-Toporov) model applies when adhesion is weaker and materials are stiffer. It predicts Hertzian contact area but adds an adhesive tensile force outside the contact.
• Friction introduces tangential tractions that alter the stress field. Under partial slip conditions, the contact edge experiences microslip while the center remains stuck, creating an annular slip zone that drives fretting fatigue.
• Coupled normal-tangential analysis is needed for accurate predictions when friction is significant.