ISO/TS 16281 — Bearing Internal Load Distribution, Contact Angle & Stiffness Calculator
ISO/TS 16281 — Rolling bearings — Methods for calculating the modified reference rating life for universally loaded bearings
Page last updated
ISO/TS 16281 is the technical specification that looks inside a rolling bearing. Where ISO 281 works from a single equivalent load P and treats the bearing as one spring with a catalogue number on it, ISO/TS 16281 solves the bearing's internal equilibrium: the rolling elements are displaced against the raceways until the loads they carry balance the applied radial force, axial force and tilting moment. From that solution come the quantities an equivalent-load method has no counterpart to — the load on each individual rolling element, the contact angle each one actually operates at, the extent of the load zone, and the bearing's stiffness.
That distinction matters because rolling-contact fatigue life goes as the inverse cube of the contact load, so the single worst-loaded element governs. A few arc-minutes of misalignment — an ordinary mounting error — concentrates the load zone onto fewer elements and raises that worst contact sharply. An equivalent-load calculation is structurally incapable of showing it.
MechanixCalc runs the internal load-distribution method for ball bearings on the bearing's own internal geometry, with the solve performed on the server rather than in the browser. It reports the per-element load, the loaded contact-angle range, the load-zone extent, and the radial, axial and tilting stiffness — the last of which is the number a shaft-deflection or rotordynamic model needs and cannot get from a catalogue.
What ISO/TS 16281 covers
- Internal load distribution — the load carried by each rolling element under combined radial, axial and moment loading
- The contact angle under load, which opens under thrust and varies around the race under combined load
- Operating clearance or interference after mounting, and its effect on the extent of the load zone
- Misalignment between the rings, imposed or equilibrated against an applied tilting moment
- Bearing stiffness — radial, axial and tilting — as the tangent of the load–displacement relation at the operating point
- The modified reference rating life L_nmref, which integrates the contact stress over the raceway — part of the standard that MechanixCalc does NOT implement
Governing formulas
Qⱼ = Kn · δⱼ^1.5where Qⱼ = load on rolling element j (N); δⱼ = total elastic approach at that element, inner and outer contacts summed (mm); Kn = the load–deflection constant of the ball/raceway pair (N/mm^1.5), computed from the general elliptical Hertz solution for the bearing's own curvatures and combined in series across the two raceways. The exponent is 3/2 for a ball (point contact); a roller's line contact uses 10/9 and is outside MechanixCalc's implementation.
δⱼ = √( (A·sin α₀ + δa + θ·(Dpw/2)·cos ψⱼ)² + (A·cos α₀ + δr·cos ψⱼ − Pd/2)² ) − Awhere A = (fᵢ + f₀ − 1)·Dw = the free distance between the raceway groove curvature centres (mm); fᵢ, f₀ = osculation, groove radius divided by ball diameter; Dw = ball diameter (mm); Dpw = pitch diameter (mm); α₀ = free contact angle; ψⱼ = azimuth of element j from the load line; δr, δa = radial and axial displacement of one ring relative to the other (mm); θ = misalignment (rad); Pd = diametral clearance (mm, negative for interference). A negative δⱼ means the element is out of contact and carries nothing — a rolling element cannot pull.
tan αⱼ = (A·sin α₀ + δa + θ·(Dpw/2)·cos ψⱼ) / (A·cos α₀ + δr·cos ψⱼ − Pd/2)where αⱼ = the contact angle element j actually operates at, which is NOT the free angle α₀. Under thrust the angle opens: on a 20-ball 15° angular-contact geometry carrying 4 kN of pure thrust it reaches about 21.4°. A model that holds α at α₀ understates both the axial capacity and the axial stiffness.
Fr = Σ Qⱼ·cos αⱼ·cos ψⱼ · Fa = Σ Qⱼ·sin αⱼ · M = Σ Qⱼ·sin αⱼ·(Dpw/2)·cos ψⱼwhere Three equations in the three unknowns (δr, δa, θ), solved simultaneously. The summations run over all Z rolling elements. Because the element loads derive from an elastic potential, the tangent of this system — the bearing's stiffness matrix — is symmetric.
Fr = Z·Q_max·J_r(ε) with ε = ½·(1 − (Pd/2)/δr)where ε = load-zone parameter (0.5 at zero clearance, below it with clearance, above it under preload); J_r(ε) = the radial load-distribution integral. At ε = 0.5 the classical result is Q_max ≈ 4.37·Fr/Z. ε describes an extent only while the load zone is smaller than the full race.
Frequently asked questions
What does ISO/TS 16281 add over ISO 281?
ISO 281 reduces the applied loads to one equivalent load P and computes life from the ratio C/P. That is the right method for selecting a bearing from a catalogue, and it cannot see inside the bearing. ISO/TS 16281 solves the internal equilibrium instead, so it produces the load on each rolling element, the contact angle each one runs at, the size of the load zone, and the bearing's stiffness. Because life goes as the inverse cube of contact load, the worst-loaded element governs — and only the internal method can tell you what it carries.
Does MechanixCalc compute the ISO/TS 16281 modified reference rating life?
No. MechanixCalc implements the standard's internal load-distribution method and reports the per-element load, the contact angles, the load zone and the bearing stiffness. The modified reference rating life L_nmref — which integrates the contact stress along the raceways and applies the standard's life-modification treatment — is not implemented. The life figures on the bearing calculator are ISO 281 L10 and Lnmh. Treat the internal load distribution as design insight, not as an ISO/TS 16281 compliance calculation.
Why does misalignment matter so much?
Misalignment tilts one ring relative to the other, which shrinks the load zone and concentrates the load on fewer rolling elements. On a 20-ball angular-contact geometry under combined radial and axial load, one milliradian — about 3.4 arc-minutes, an ordinary mounting error — takes the load zone from eleven loaded elements to nine and raises the worst contact by roughly 14 %. Because rolling-contact fatigue life goes as Q⁻³, that is about a third of the life. That figure is for a shaft stiff enough to IMPOSE the misalignment, which is the usual case and the calculator's default; a bearing free to tilt instead finds its own alignment, and there the same milliradian can lower the worst contact rather than raise it. Which boundary condition applies is a real modelling choice, and the calculator names the one it used on every result. An equivalent-load calculation shows no change under either.
What internal geometry do I need to run this?
The number of rolling elements, the ball diameter, the pitch diameter, the free contact angle, the raceway osculation (groove radius divided by ball diameter) and the mounted diametral clearance. None of those appear on the catalogue page a bearing's load ratings come from, so MechanixCalc can estimate them from the bore and outside diameter to get you started. That estimate is shown, not hidden, and it is editable: a real ball count depends on the series, the cage and the filling method, and two bearings with identical boundary dimensions routinely differ. Replace the estimate with the bearing's own geometry before relying on the result.
Why does the calculator report a bearing stiffness?
Because nothing else can supply it. A shaft-deflection model, a rotordynamic model and a thermal-growth model all need the bearing's radial, axial and tilting stiffness as an input, and a bearing catalogue publishes load ratings rather than stiffness. The stiffness is the tangent of the load–displacement relation at the operating point, so it falls straight out of the internal equilibrium — and it changes with load, clearance and preload, which is exactly why a single tabulated figure would not do.
Why is the contact angle not the one in the bearing designation?
The angle in the designation is the FREE contact angle, measured with no load on the bearing. Under thrust the rings shift axially and the contact point climbs the groove, so the operating angle opens — on a 15° row under 4 kN of pure thrust it reaches about 21.4°. Under combined load it also varies from element to element around the race. Holding the angle at its free value understates the axial capacity and the axial stiffness.
Does it handle roller bearings?
No. MechanixCalc's implementation covers ball bearings, where each contact is a Hertzian point contact with a 3/2 load-deflection exponent. A roller's line contact needs the roller to be divided into laminae along its length, with the load distribution solved along that length as well as around the race — that is a different calculation and it is not implemented. The ISO 281 life calculation on the bearing page covers roller bearings as usual.
What is left out, and which way does it err?
Centrifugal loading of the rolling elements is not modelled. It throws the balls outward, so it ADDS to the outer-raceway contact — roughly 10 % of the worst contact at 20 000 rpm on a 90 mm pitch diameter, and more above that. That omission is not conservative for high-speed duty, which is why every result says so. The rings are also treated as rigid: for a stiff housing that is conservative, because a compliant ring would spread the load over more elements, but for a thin-section ring in a soft or out-of-round housing, ovalisation concentrates the load instead. Gyroscopic and cage effects and any thermal change in clearance are outside the model.
Related standards
Run a ISO/TS 16281 calculation on your own numbers
Free 30-minute preview — no sign-up. A free 14-day account trial unlocks every tool and the branded PDF report.
Open the Bearing AnalysisUsing MechanixCalc at work? See plans & pricing — one subscription unlocks all 50 calculators, PDF reports and saved projects.