Quick answer For spherical contact, the Hertzian contact radius is a=(3F·R*/(4E*))^(1/3) (R*=equivalent radius of curvature, E*=equivalent elastic modulus), and the maximum contact pressure is pmax=3F/(2πa²). The maximum shear stress occurs below the surface at z≈0.48a and is ≈0.31·pmax, becoming the origin of fatigue failure.
Based on Hertz contact theory, compute contact pressure distribution, contact radius, and subsurface maximum shear stress for sphere and cylinder contacts in real time.
Parameter Settings
Contact Type
R₁ — Radius 1
mm
R₂ — Radius 2
mm
E₁ — Young's Modulus (Body 1)
GPa
E₂ — Young's Modulus (Body 2)
GPa
ν₁ — Poisson's Ratio 1
ν₂ — Poisson's Ratio 2
P — Applied Load
N
L — Contact Length (Cylinder)
mm
Results
—
μm
Contact Half-Width a
—
MPa
Max Contact Pressure p₀
—
MPa
Max Shear Stress τ_max
—
μm
τ_max Depth z*
Contact Pressure Distribution p(x)
Subsurface Stress Distribution σz, τmax vs Depth z
Hertz contact radius $a$ (m): $F$ is the applied load (N), $R^{*}$ is the equivalent radius of curvature (m), and $E^{*}$ is the equivalent elastic modulus (Pa).
$$p_0 = \frac{3F}{2\pi a^2}$$
Maximum contact pressure (Pa): the pressure at the center of the contact ellipse. For structural steel it can reach 1–3 GPa.
Indentation depth (m): the approach between the centers of the two spheres. A nonlinear elastic deformation proportional to the load raised to the 2/3 power.
What exactly is Hertz contact stress? I see it mentioned for bearings and gears, but what's happening at the microscopic level when two curved surfaces touch?
🎓
Basically, it's the localized stress that develops when two elastic bodies are pressed together. Unlike a uniform load, the pressure isn't spread evenly—it peaks at the center of the contact patch and drops to zero at the edges, forming a semi-ellipsoidal "pressure bulb." In practice, this is why a steel ball bearing can support a huge load on a tiny contact area without immediate failure. Try moving the "Applied Force" slider in the simulator above to see how the contact area grows and the peak pressure increases.
🙋
Wait, really? So the stress is actually under the surface? That seems counterintuitive. How do the material properties, like changing from steel to aluminum, affect this?
🎓
Exactly right! For many materials, the maximum shear stress occurs not at the surface but about half a contact radius below it. This is where subsurface cracks often start. The material's stiffness (Young's Modulus, $E$) and "squishiness" (Poisson's ratio, $\nu$) are critical. A softer material, like aluminum, will deform more, creating a larger contact area and thus lowering the peak pressure for the same force. You can test this by adjusting the $E_1$ and $E_2$ sliders for the two bodies in the simulator and watching the contact circle change size.
🙋
Okay, that makes sense. So what's the big deal with the "equivalent radius" ($R^*$) I see in the formulas? Why combine the curvatures of both bodies?
🎓
Great question! The equivalent radius tells you how "concentrated" the contact is. Think of it like the effective curvature of a single body pressing against a flat plane. If you press two identical spheres together, $R^*$ is just half of one sphere's radius—the contact is quite sharp. But if you press a small sphere against the inside of a huge, gently curved bowl, $R^*$ becomes large, meaning the contact patch is much wider and gentler. This is key for designing joints. Use the "Contact Type" selector to switch between sphere-sphere and sphere-flat contacts and see how $R^*$ changes instantly.
Physical Model & Key Equations
The core of Hertzian contact theory is calculating the size of the contact area and the maximum pressure at its center. For two spheres in contact, the contact patch is a circle. The theory assumes perfectly smooth, elastic materials and small, frictionless contact.
$$a = \sqrt[3]{\frac{3FR^{*}}{4E^{*}}}$$
Here, $a$ is the contact radius , $F$ is the applied normal force, $R^{*}$ is the equivalent radius of curvature, and $E^{*}$ is the equivalent Young's modulus. A larger $a$ means the load is spread over a bigger area.
The maximum contact pressure, $p_0$, occurs at the center of this circle. It is not simply force divided by area; due to the elliptical pressure distribution, it's 1.5 times the average pressure.
This peak Hertzian pressure is the critical value for predicting surface yielding or fatigue. The material properties are combined into the equivalent modulus $E^{*}$, and the geometries into the equivalent radius $R^{*}$, defined as:
For a sphere on a flat surface, $R_2 = \infty$, so $R^* = R_1$. A negative $R_2$ denotes a concave surface (like a ball in a socket).
Frequently Asked Questions
Select the "Sphere – Flat" (or "Cylinder – Flat") contact type from the dropdown. The flat side is then automatically treated as having an infinite radius of curvature and the R₂ input is hidden, so you only need to enter the sphere's (or cylinder's) radius as R₁.
The contact area may be too small relative to the load. Check whether the radius of curvature is extremely small or if a material with a very high Young's modulus has been set. Also, ensure that the unit system (e.g., N and mm, MPa) is consistent.
No. This simulator accepts only positive radii of curvature (1–500 mm), covering convex–convex and convex–flat contact. Conforming contact with a concave surface (e.g., a ball in a socket) would require a negative radius in the equivalent-radius formula R* and cannot be simulated directly with this tool.
The position (depth) of the maximum shear stress occurring just below the contact surface serves as an indicator for predicting the location of fatigue failure (such as spalling or pitting). For example, in the design of rolling bearings or gears, if inclusions exist at this depth, the risk of damage increases, making it a useful reference for material selection and heat treatment.
Real-World Applications
Ball & Roller Bearings: This is the classic application. Hertz theory is used to calculate the contact stress between rolling elements and raceways to predict fatigue life (known as the L10 life). Engineers select materials, hardness, and curvature to keep subsurface shear stresses below the material's endurance limit.
Gear Tooth Design: The contact between meshing gear teeth is modeled as contact between two cylinders (for spur gears). Hertz contact stress, often called "contact fatigue" or "pitting resistance," is a primary failure mode. Surface treatments like carburizing are used to create a hard case that resists this pressure.
Railway Wheels & Rails: The contact patch between a steel train wheel and the rail is a small ellipse. Hertz analysis helps understand wear patterns, plastic deformation (leading to "rail head checking"), and the forces that can cause derailment on curves.
Biomechanics & Prosthetics: In artificial hip or knee joints, the metal or ceramic ball contacts a polymer or ceramic socket. Hertz calculations, though simplified, help estimate contact pressures to minimize wear debris generation and prevent premature failure of the implant.
Common Misconceptions and Points to Note
When you start using this simulator, there are several pitfalls that beginners in CAE often fall into. A major misconception is the idea that "if you make the material stronger (using high-strength steel), it can handle any high load". While material strength is indeed important, Hertzian contact stress is a "contact pressure". For example, if you double the load, the contact radius $a \propto P^{1/3}$ increases only by about 1.26 times, but the maximum contact pressure $p_0$ jumps up by about 1.59 times, proportional to $P^{2/3}$. Often, simply improving material strength isn't enough, and a fundamental design review—increasing the shape (radius of curvature) to enlarge the contact area—is necessary.
Next, a point of caution regarding parameter settings. While the simulator treats a "flat surface" as having an infinite radius of curvature, in practice, perfectly flat surfaces hardly ever exist. For instance, in the contact between a ball bearing's ball and its race, the raceway groove's radius $R_2$ is set to a value slightly larger than the ball's radius $R_1$ (e.g., 1.05 times). If you consider this a "flat" surface, the calculated contact pressure will be significantly higher (i.e., overly conservative on the safe side) than in reality, potentially leading to over-design.
Finally, it's important to understand the assumptions of this tool. Hertz theory fundamentally assumes perfectly elastic bodies, smooth surfaces, and behavior within the elastic limit. This means that if the load is high enough to cause plastic deformation, if surfaces are rough, or if fatigue from repeated loading is the primary concern, these results are merely a "first approximation". For example, if a hardened steel surface has a tiny dent, stress concentration can occur there, potentially initiating failure at a location completely different from the calculated maximum shear stress point. You should treat simulation results as idealized "reference values", and it's crucial to maintain an attitude of verifying them with physical tests or more advanced CAE analysis.
Standard / reference: Hertz contact theory (H. Hertz, 1882; K.L. Johnson, Contact Mechanics, 1985). Sphere contact: radius \(a = (3FR^{*}/4E^{*})^{1/3}\), peak pressure \(p_0 = 3F/(2\pi a^2)\). Reduced modulus \(1/E^{*}=(1-\nu_1^2)/E_1+(1-\nu_2^2)/E_2\), reduced curvature \(1/R^{*}=1/R_1+1/R_2\). Subsurface max shear \(\tau_\text{max}\approx0.31\,p_0\) at \(z\approx0.48a\) (sphere), \(\approx0.30\,p_0\) at \(z\approx0.78a\) (cylinder).
Model assumptions: Linear-elastic, isotropic-homogeneous, frictionless, contact small relative to curvature (\(a \ll R\)), surface roughness neglected. Pressure is the semi-elliptical \(p(x)=p_0\sqrt{1-(x/a)^2}\). The subsurface stress profiles are first-order on-axis approximations.
Scope & limits: Defaults (\(R_1{=}R_2{=}20\,\text{mm}, E{=}200\,\text{GPa}, \nu{=}0.3, P{=}1000\,\text{N}\)) give \(a{=}408.7\,\mu\text{m},\ p_0{=}2.86\,\text{GPa}\). Not valid under plastic flow (\(p_0\) above yield), for \(a/R>0.1\), or when roughness is comparable to a. Real designs add dynamic, roughness and lubrication (EHL) corrections.
How to Use
Enter the radius of the first body (R1) in millimeters—for sphere-sphere contact, use radii; for sphere-flat, choose the "Sphere – Flat" contact type (R2 is then treated as infinite automatically).
Input elastic moduli (E1, E2) in GPa for each material: steel typically 200 GPa, aluminum 70 GPa, glass 70 GPa.
Specify the normal contact load in Newtons, then click Calculate to visualize pressure distribution and maximum subsurface stress depth (approximately 0.5 contact radius below surface).
Worked Example
Two steel spheres (E = 200 GPa) with radii R1 = 10 mm and R2 = 15 mm under a 5 kN normal load: the equivalent radius is R* = 6 mm and the contact radius is a ≈ 0.59 mm, giving a peak Hertzian pressure of approximately 6.9 GPa at the contact center. Maximum shear stress occurs at depth z ≈ 0.48a ≈ 0.28 mm, reaching about 2.1 GPa (≈ 0.31 p₀). This stress concentration is critical for bearing wear prediction and fatigue life assessment in rolling element bearings.
Practical Notes
For cylinder-cylinder contact (crossed axes), use equivalent radii: 1/Req = 1/R1 + 1/R2 to convert to sphere-equivalent problem.
Material pairing affects stress: steel-on-ceramic (200 GPa vs. 300 GPa) produces asymmetric subsurface stresses; use reduced modulus E* = 1/[(1-ν1²)/E1 + (1-ν2²)/E2].
Hertz theory assumes elastic deformation, frictionless contact, and smooth surfaces; real contacts with lubrication or roughness require traction stress analysis.