Interactive Calculator

A shaft–hub connection transmits torque through two parallel mechanisms: friction from the interference fit pressure, and bonding shear from an adhesive film applied at the contact interface. The classical approach sizes the hub by experience or iteration — pick a wall thickness, check Tresca, adjust. Croccolo, De Agostinis & Vincenzi (2012) showed that for a given shaft diameter and material pair, an analytically optimal hub ratio \(Q_H = D_C / D_{He}\) exists that maximizes the torque-to-mass ratio of the joint. They derived a closed-form expression for this optimum that depends on only two dimensionless parameters: \(\phi\) (density–geometry ratio) and \(\chi\) (adhesive–strength ratio).

This tool implements that optimization as a six-step calculation pipeline. The pipeline takes shaft and hub materials, geometry, friction coefficient, and adhesive strength as input, and returns the optimal hub ratio \(Q_{H\_opt}\) , the corresponding Tresca-limited contact pressure \(p_C\) , the coupling geometry (\(D_C\) , \(D_{He}\) , \(L_C\) , \(Z\) ), and the design merit \(\text{DFLS}_T\) in kN·m/kg. A Python notebook ready to run on Google Colab is provided alongside the derivation, so every intermediate value can be verified.

Quick Example

A steel 39NiCrMo3 shaft must transmit a static torque \(T = 1.0\) kN·m through a hub in aluminum EN-AW6082. The allowable shaft shear stress is \(\tau_{allow} = 350\) MPa, the friction coefficient is \(\mu_T = 0.4\) . Three progressively more efficient designs:

Scenario (a)Scenario (b)Scenario (c)
ConfigurationSolid shaft, no adhesiveHollow shaft, no adhesiveHollow shaft + adhesive
\(Q_S\)00.70.7
\(\tau_{ad}\) (MPa)0010

Design merit: (a) 5.92 kN·m/kg → (b) 8.86 kN·m/kg → (c) 11.68 kN·m/kg

Going from a solid shaft to a hollow shaft with adhesive nearly doubles the torque-to-mass ratio. The pipeline quantifies exactly how much each design choice contributes.

Pipeline summary:

StepOperationKey output
1Shaft sizing from torsion\(D_C\)
2Normalizing parameters\(\phi\) , \(\chi\)
3Optimal hub ratio\(Q_{H\_opt}\)
4Tresca pressure limit\(p_C\)
5Joint geometry\(L_C\) , \(D_{He}\) , \(Z\)
6Design meritmass, \(\text{DFLS}_T\)

Pipeline Overview

The pipeline is fully closed-form — no iteration, no solver. Each step feeds the next with explicit formulas. The normalizing parameters \(\phi\) and \(\chi\) collapse all material and geometry combinations into a two-parameter family, so the same pipeline applies to any shaft–hub pair from magnesium to steel.


Step 1 — Shaft Sizing

Purpose. Determine the minimum coupling diameter \(D_C\) from the required torque and shaft material.

SymbolDescriptionUnit
\(T\)Required torqueN·mm
\(\tau_{allow}\)Allowable shaft shear stressMPa
\(Q_S = D_{Si}/D_C\)Shaft bore ratio (0 for solid shaft)

From the torsional section modulus of a hollow circular cross-section:

\[D_C = \left(\frac{16\,T}{\pi\,\tau_{allow}\,(1 - Q_S^3)}\right)^{1/3}\]

Note: the paper’s Eq. (30) uses \((1 - Q_S^3)\) rather than the exact \((1 - Q_S^4)\) . For \(Q_S = 0\) both forms are identical. For \(Q_S = 0.7\) the difference is about 5% on \(D_C\) . This tool follows the paper’s convention for verification consistency.

SymbolDescriptionUnit
\(D_C\)Coupling diametermm

Step 2 — Normalizing Parameters

Purpose. Reduce the material and geometry inputs to two dimensionless groups that fully characterize the optimization problem.

SymbolDescriptionUnit
\(Q_S\)Shaft bore ratio
\(\rho_S\) , \(\rho_H\)Shaft and hub densitieskg/mm³
\(\tau_{ad}\)Adhesive shear strengthMPa
\(S_{y\_H}\)Hub yield strengthMPa
\[\phi = (1 - Q_S^2)\,\frac{\rho_S}{\rho_H}\]
\[\chi = \frac{2\,\tau_{ad}}{S_{y\_H}}\]

\(\phi\) measures the relative weight contribution of the shaft: high \(\phi\) means the shaft dominates the joint mass (heavy material, solid section), pushing the optimum toward thinner hubs. \(\chi\) measures how much “free” shear capacity the adhesive adds relative to the hub’s yield limit. Two problems with identical \((\phi, \chi, \mu)\) have identical \(Q_{H\_opt}\) regardless of the actual materials and diameters.

SymbolDescriptionUnit
\(\phi\)Density–geometry parameter
\(\chi\)Adhesive–strength parameter

Step 3 — Optimal Hub Ratio

Purpose. Find the hub aspect ratio \(Q_H\) that maximizes the torque-to-mass ratio \(\text{DFLS}_T\) .

SymbolDescriptionUnit
\(\phi\) , \(\chi\)Normalizing parameters
\(\mu\)Friction coefficient

Without adhesive (\(\chi = 0\) ), Eq. (25) of the paper:

\[Q_{H\_opt} = \sqrt{\frac{1 - \sqrt{\phi}}{1 - \phi}}\]

With adhesive (\(\chi > 0\) ), Eq. (29) of the paper:

\[Q_{H\_opt} = \sqrt{\frac{1 - \sqrt{\dfrac{1}{1-\chi^2}\left[\phi - \dfrac{\chi}{\mu}\sqrt{1-\chi^2}\,(1-\phi)\right]}}{1 - \phi}}\]

For \(\chi = 0\) the second form reduces identically to the first. For \(\phi \to 1\) the denominator vanishes but the limit is finite: \(Q_{H\_opt} \to 1/\sqrt{2} \approx 0.707\) .

The approximation \(\mu^2 \ll 1/(1-Q_H^2)^2\) is used in deriving both forms. The paper’s Table 2 confirms the error on \(Q_{H\_opt}\) stays below 5% for all practical \((\phi, \mu)\) combinations.

SymbolDescriptionUnit
\(Q_{H\_opt}\)Optimal hub ratio

Step 4 — Tresca Pressure Limit

Purpose. Compute the maximum contact pressure the hub can sustain at the optimal \(Q_H\) before yielding at the inner surface.

SymbolDescriptionUnit
\(Q_{H\_opt}\)Optimal hub ratio
\(\mu\)Friction coefficient
\(\tau_{ad}\)Adhesive shear strengthMPa
\(S_{y\_H}\)Hub yield strengthMPa

The Tresca criterion at the hub bore, combined with friction and adhesive shear on the interface, gives a quadratic in \(p_C\) . Defining:

\[a = 4\left(\frac{1}{(1-Q_H^2)^2} + \mu^2\right), \quad b = 8\,\mu\,\tau_{ad}, \quad c = 4\,\tau_{ad}^2 - S_{y\_H}^2\]
\[p_C = \frac{-b + \sqrt{b^2 - 4ac}}{2a}\]

For \(\tau_{ad} = 0\) this simplifies to:

\[p_C = \frac{S_{y\_H}}{2\sqrt{\dfrac{1}{(1-Q_H^2)^2} + \mu^2}}\]
SymbolDescriptionUnit
\(p_C\)Tresca-limited contact pressureMPa

Step 5 — Joint Geometry

Purpose. Size the coupling length, hub outer diameter, and required interference.

SymbolDescriptionUnit
\(T\)Required torqueN·mm
\(D_C\)Coupling diametermm
\(p_C\)Contact pressureMPa
\(\mu\) , \(\tau_{ad}\)Friction coefficient, adhesive strength—, MPa
\(Q_{H\_opt}\) , \(Q_S\)Hub and shaft bore ratios
\(E_S\) , \(\nu_S\) , \(E_H\) , \(\nu_H\)Elastic constantsMPa, —

Coupling length from torque equilibrium, Eq. (6):

\[L_C = \frac{2\,T}{(\mu\,p_C + \tau_{ad})\,\pi\,D_C^2}\]

Hub outer diameter:

\[D_{He} = \frac{D_C}{Q_{H\_opt}}\]

Required diametral interference from Lamé compatibility, Eq. (3):

\[Z = p_C\,D_C\left[\frac{1}{E_H}\left(\frac{1+Q_H^2}{1-Q_H^2}+\nu_H\right) + \frac{1}{E_S}\left(\frac{1+Q_S^2}{1-Q_S^2}-\nu_S\right)\right]\]
SymbolDescriptionUnit
\(L_C\)Coupling lengthmm
\(D_{He}\)Hub outer diametermm
\(Z\)Required diametral interferencemm

Step 6 — Design Merit

Purpose. Compute the joint mass and the design function for lightweight structures.

SymbolDescriptionUnit
\(D_C\) , \(D_{Si}\) , \(D_{He}\)Coupling, bore, and hub diametersmm
\(L_C\)Coupling lengthmm
\(\rho_S\) , \(\rho_H\)Shaft and hub densitieskg/mm³
\(T\)Required torqueN·m
\[m = \frac{\pi}{4}\left[(D_C^2 - D_{Si}^2)\,\rho_S + (D_{He}^2 - D_C^2)\,\rho_H\right] L_C\]
\[\text{DFLS}_T = \frac{T}{m} \quad [\text{Nm/kg}]\]

\(\text{DFLS}_T\) does not depend on \(D_C\) : at fixed materials and ratios, a joint transmits torque with the same mass efficiency regardless of absolute scale. This is a structural property of the torque case (where \(T \propto D_C^2\) and \(m \propto D_C^2\) ).

SymbolDescriptionUnit
\(m\)Joint masskg
\(\text{DFLS}_T\)Torque-to-mass ratioN·m/kg

Numerical Case — Croccolo et al. (2012), §4.1

Input

PropertyShaft (39NiCrMo3 Steel)Hub (EN-AW6082 Al)
\(\rho\) (kg/mm³)\(7.87 \times 10^{-6}\)\(2.75 \times 10^{-6}\)
\(E\) (MPa)207 00069 000
\(\nu\)0.290.33
\(S_y\) (MPa)304

\(T = 1000\) N·m, \(\tau_{allow} = 350\) MPa, \(\mu_T = 0.4\) .

Scenario (a) — Solid shaft, no adhesive

\(Q_S = 0\) , \(\tau_{ad} = 0\) MPa.

Step 1. \(D_C = (16 \times 10^6 / (\pi \cdot 350 \cdot 1))^{1/3} = 24.4\) mm. Paper rounds to 25 mm.

Step 2. \(\phi = (1 - 0) \cdot 7.87/2.75 = 2.862\) . \(\chi = 0\) .

Step 3. \(Q_{H\_opt} = \sqrt{(1 - \sqrt{2.862})/(1 - 2.862)} = 0.610\) . Paper: 0.61. ✓

Step 4. \(p_C = 304 / (2\sqrt{1/(1-0.61^2)^2 + 0.16}) = 92.6\) MPa. Paper: 93 MPa. ✓

Step 5. \(L_C = 2 \times 10^6 / (0.4 \cdot 93 \cdot \pi \cdot 25^2) = 27.4\) mm. Paper: 28 mm. ✓. \(D_{He} = 25/0.61 = 41.0\) mm. \(Z = 0.093\) mm.

Step 6. \(m = 168.9\) g. \(\text{DFLS}_T = 1000/0.169 = 5920\) N·m/kg ≈ 5.92 kN·m/kg. Paper: 5.9. ✓

Scenario (b) — Hollow shaft, no adhesive

\(Q_S = 0.7\) , \(\tau_{ad} = 0\) MPa.

Step 1. \(D_C = (16 \times 10^6 / (\pi \cdot 350 \cdot 0.657))^{1/3} = 28.1\) mm. Paper: 28 mm. ✓

Step 2. \(\phi = (1 - 0.49) \cdot 2.862 = 1.460\) . \(\chi = 0\) .

Step 3. \(Q_{H\_opt} = \sqrt{(1 - \sqrt{1.46})/(1 - 1.46)} = 0.673\) . Paper: 0.67. ✓

Step 4. \(p_C = 81.2\) MPa. Paper: 82 MPa. ✓

Step 5. \(L_C = 24.8\) mm. Paper: 25 mm. ✓. \(D_{He} = 41.6\) mm. \(Z = 0.128\) mm.

Step 6. \(m = 112.9\) g. \(\text{DFLS}_T = \) 8.86 kN·m/kg. Paper: 8.9. ✓

Scenario (c) — Hollow shaft + adhesive

\(Q_S = 0.7\) , \(\tau_{ad} = 10\) MPa.

Step 2. \(\phi = 1.460\) (same as b). \(\chi = 2 \cdot 10 / 304 = 0.0658\) .

Step 3. \(Q_{H\_opt} = 0.725\) . Paper: 0.73. ✓

Step 4. \(p_C = 69.8\) MPa. Paper: 70 MPa. ✓

Step 5. \(L_C = 21.3\) mm. Paper: 21 mm. ✓. \(D_{He} = 38.7\) mm. \(Z = 0.126\) mm.

Step 6. \(m = 85.6\) g. \(\text{DFLS}_T = \) 11.68 kN·m/kg. Paper: 11.7. ✓

Summary

(a) Solid, dry(b) Hollow, dry(c) Hollow + adhesive
\(\phi\)2.8621.4601.460
\(\chi\)000.066
\(Q_{H\_opt}\)0.6100.6730.725
\(p_C\) (MPa)92.681.269.8
\(D_{He}\) (mm)41.041.638.6
\(L_C\) (mm)27.424.821.3
mass (g)168.9112.985.6
\(\text{DFLS}_T\) (kN·m/kg)5.928.8611.68

All values within 3% of the paper. The 4.2% discrepancy on \(p_C\) for scenario (c) is structural: the paper applies the simplified Tresca form without \(\mu^2\) , while the calculator uses the full quadratic. \(\text{DFLS}_T\) is unaffected (0.2% deviation).


DFLS vs Hub Ratio

The chart in the calculator shows \(\text{DFLS}_T\) as a function of \(Q_H\) for the three canonical scenarios. Each curve rises from zero (infinitely thick hub), peaks at the analytic optimum \(Q_{H\_opt}\) , and falls back to zero (infinitely thin hub, Tresca forces \(p_C \to 0\) ). The adhesive shifts the peak rightward (thinner hub is optimal) and upward (higher merit). The three faded reference curves remain visible as context while the cyan curve tracks the current input configuration.


How to Use the Notebook

The Python notebook implements the full six-step pipeline and runs directly on Google Colab with no installation required (numpy and matplotlib are preinstalled).

Open in Colab

Modify only the cell marked ✏️ Input. The parameters are:

ParameterSymbolUnitNote
Shaft densityrho_Skg/mm³e.g. \(7.87 \times 10^{-6}\) for steel
Shaft Young’s modulusE_SMPa
Shaft Poisson’s rationu_S
Hub densityrho_Hkg/mm³
Hub Young’s modulusE_HMPa
Hub Poisson’s rationu_H
Hub yield strengthSy_HMPaTresca-governing
Shaft bore ratioQ_S0 for solid shaft; ≤ 0.8
Required torqueTN·mStatic
Allowable shaft sheartau_allowMPa
Friction coefficientmuTypically 0.1–0.4
Adhesive shear strengthtau_adMPa0 for dry joint

Known limitations. The model assumes plane stress, axisymmetric geometry, and linear-elastic isotropic materials (\(\sigma_{VM} < \sigma_y\) must hold). Hub yielding at the inner surface governs (Tresca criterion). The approximation \(\mu^2 \ll 1/(1-Q_H^2)^2\) is embedded in the \(Q_{H\_opt}\) formulas (error < 5% per Table 2 of the paper). Static torque only — no fatigue, no dynamic loads. Shaft buckling for hollow shafts (\(Q_S > 0.8\) ) is not checked. The paper’s Eq. (30) uses \((1 - Q_S^3)\) for the torsional section modulus instead of the exact \((1 - Q_S^4)\) .


References

[1] Croccolo, D., De Agostinis, M. & Vincenzi, N. (2012). Design and optimization of shaft–hub hybrid joints for lightweight structures. International Journal of Mechanical Sciences, 56(1), 77–85.

[2] Croccolo, D. & Vincenzi, N. (2009). A generalized theory for shaft–hub couplings. Proc. IMechE Part C: J. Mechanical Engineering Science, 223, 2231–2239.