A spiral bevel gear pair transmits rotation between intersecting shafts with a high load capacity, smooth conjugate motion and low noise, which is why it dominates final drives in automobiles, rotorcraft and aircraft engine accessories, rolling mills, power plant equipment and rotary drilling rigs. Its weak point is dynamic: mesh stiffness varies through the cycle, loaded teeth enter and leave contact progressively, and transmission error excites the housing and the shafting. Manufacturing and assembly error add a second, often larger, excitation. When the centre of mass of a gear does not coincide with its axis of rotation, the resulting eccentricity generates centrifugal and inertial forces in the plane of rotation, the contact point drifts away from its theoretical position, and the vibration of the whole train rises.
Two routes to lower vibration
Countermeasures divide into active and passive families. Active measures attack the source: better tooth flank finishing, revised gear parameters, and profile or lead modification that compensates for loaded deflection. They are effective but computationally demanding, expensive and slow to iterate, because each variant requires a new cutter or a new grinding program. Passive measures accept the excitation and dissipate the resulting kinetic energy in an added device — a squeeze-film damper, a damping ring shrunk onto the gear web, or wire-rope isolators.
Particle damping belongs to the passive family and is unusually cheap. A set of loose particles is sealed inside a cavity that is rigidly attached to the vibrating structure, forming a closed, strongly nonlinear coupled oscillator. Energy leaves the primary structure through two channels: inelastic impacts among particles and between particles and the cavity wall, and friction at the same contacts. Because dissipation does not depend on a single tuned frequency, particle dampers are broadband, and their geometry can be adapted freely to the installation space. Reported implementations include flexible fabric-wrapped “bean bag” dampers placed in the cavity of a main generator end cover on a diesel locomotive, rectangular box dampers bolted to the free end of an L-shaped cantilever, cube dampers on free-free beams verified by finite element analysis and test, and flat box dampers magnetically attached to the frame of a high-speed train power pack. The damper can also be the structure itself: symmetric grooves machined at the hub of a wheel-like part, sub-cavities inside a hollow shaft, or a cavity bored into the front of a boring bar.
Inside a gearbox the established practice is to mount the damper on the housing, i.e. in a stationary field. Work on rotating gears is thinner. One line of research machines through-holes in the gear web as cavities and closes them with baffles; the discrete element method is then used to trace particle motion under centrifugal load and to explain the dissipation mechanism for cylindrical gears and for spiral bevel gears. A related design is a particle dry-friction primary/secondary damping ring seated in a groove cut specifically into the bevel gear web. Locating the damper on the gear puts it close to the source, but every such solution changes the gear body, requires sealing components and mounting holes, and therefore raises doubts about tooth-root strength while adding manufacturing cost. A damper carried on the gear shaft avoids all of this, and it is this configuration that is examined below.
Dynamic model of a spiral bevel gear pair with eccentric mass
The pair is represented by a lumped-parameter model with the origin at the intersection of the two shaft axes. The mesh is collapsed into a spring-damper element with time-varying stiffness k(t) and mean mesh damping cm; bearings are elastic supports with stiffness kij and damping cij (i = x, y, z; j = 1, 2), assumed isotropic; tooth friction is neglected, and the static transmission error is omitted. Each gear then retains three translational coordinates and one rotation about its own axis, giving eight degrees of freedom:
$$\begin{aligned}
m_1\ddot{x}_1 + c_{x1}\dot{x}_1 + k_{x1}x_1 &= -F_x\\
m_1\ddot{y}_1 + c_{y1}\dot{y}_1 + k_{y1}y_1 &= -F_y + F_{ly1} – F_{gy1}\\
m_1\ddot{z}_1 + c_{z1}\dot{z}_1 + k_{z1}z_1 &= -F_z + F_{lz1} – F_{gz1}\\
I_1\ddot{\theta}_1 + F_m r_1 &= T_1\\
m_2\ddot{x}_2 + c_{x2}\dot{x}_2 + k_{x2}x_2 &= F_x + F_{lx2} – F_{gx2}\\
m_2\ddot{y}_2 + c_{y2}\dot{y}_2 + k_{y2}y_2 &= F_y\\
m_2\ddot{z}_2 + c_{z2}\dot{z}_2 + k_{z2}z_2 &= F_z + F_{lz2} – F_{gz2}\\
I_2\ddot{\theta}_2 – F_m r_2 &= -T_2
\end{aligned}$$
where m1, m2 are equivalent masses including the shaft sections, I1, I2 are moments of inertia, T1, T2 are driving and load torque, Fx, Fy, Fz are the components of the normal dynamic mesh force, Fm is the normal mesh force itself, and Flij and Fgij are the radial components of centrifugal and inertial force produced by the eccentricity.
Because a spiral bevel pair meshes progressively and its contact ratio is large, the stiffness does not jump as it does in a spur pair; it oscillates about a mean value, so a mean stiffness plus a first-harmonic term is adequate:
$$k(t) = k_m + k_{kl}\cos\left(l\omega_n t + \psi_{kl}\right)$$
with l = 1, km the mean mesh stiffness, ωn the mesh frequency and the initial phase ψkl taken as zero. With an eccentricity present, the relative displacement along the line of action becomes
$$\lambda_m = \left[a_1\left(x_1 + r_{p1}\theta_1 – e_1\cos\theta_1\right) – a_2\left(x_2 + r_{p2}\theta_2 – e_2\cos\theta_2\right)\right] + \left[b_1\left(y_1 + r_{p1}\theta_1 – e_1\sin\theta_1\right) – b_2\left(y_2 + r_{p2}\theta_2 – e_2\sin\theta_2\right)\right] – \left(c_1 z_1 – c_2 z_2\right)$$
where θ1 = ω1t − θp and θ2 = ω2t − θg are the eccentricity angles, and the geometric coefficients are
$$a_q = \cos\delta_q\sin\alpha_n,\quad b_q = \cos\delta_q\cos\alpha_n\sin\beta_q + \sin\delta_q\sin\alpha_n$$
$$c_q = \sin\delta_q\cos\alpha_n\sin\beta_q + \cos\delta_q\sin\alpha_n,\quad r_{pq} = \frac{r_q}{\cos\alpha_q}$$
with δq the pitch cone angle, βq the mean spiral angle, αn the normal pressure angle, rq the base circle radius and αq the mesh angle of wheel q. The normal dynamic mesh force and its coordinate components follow as
$$F_m = k(t)\lambda_m + c_m\dot{\lambda}_m$$
$$\begin{aligned}
F_x &= -F_m\left(\sin\alpha_n\cos\delta_2 + \cos\alpha_n\sin\beta_2\sin\delta_2\right)\\
F_y &= F_m\left(\sin\alpha_n\sin\delta_2 – \cos\alpha_n\sin\beta_2\cos\delta_2\right)\\
F_z &= -F_m\cos\alpha_n\cos\beta_2
\end{aligned}$$
The eccentricity adds, in each plane of rotation, centrifugal components that scale with the square of the angular velocity,
$$\begin{aligned}
F_{ly1} &= m_1 e_1\left(\omega_1 + \dot{\theta}_p\right)^2\cos\theta_1, & F_{lz1} &= m_1 e_1\left(\omega_1 + \dot{\theta}_p\right)^2\sin\theta_1\\
F_{lx2} &= m_2 e_2\left(\omega_2 + \dot{\theta}_g\right)^2\cos\theta_2, & F_{lz2} &= m_2 e_2\left(\omega_2 + \dot{\theta}_g\right)^2\sin\theta_2
\end{aligned}$$
and inertial components that scale with the angular acceleration of the torsional coordinate,
$$\begin{aligned}
F_{gy1} &= m_1 e_1\ddot{\theta}_1\cos\theta_1, & F_{gz1} &= m_1 e_1\ddot{\theta}_1\sin\theta_1\\
F_{gx2} &= m_2 e_2\ddot{\theta}_2\cos\theta_2, & F_{gz2} &= m_2 e_2\ddot{\theta}_2\sin\theta_2
\end{aligned}$$
Simulated response to eccentricity
The model was integrated numerically with a Runge–Kutta scheme for a 27/79 tooth pair whose geometry is listed in Table 1, driven at 60 kW and 3000 r/min against a 200 N·m load. Identical eccentricities of 50, 100 and 200 μm were assigned to both members so that the effect of a single parameter could be isolated.
| Parameter | Pinion | Gear |
|---|---|---|
| Number of teeth z | 27 | 79 |
| Hand | left | right |
| Face width b (mm) | 22.8 | 22.8 |
| Pitch cone angle δ (°) | 18.87 | 71.13 |
| Shaft angle (°) | 90 | 90 |
| Transverse module m (mm) | 3.15 | 3.15 |
| Mean spiral angle β (°) | 35 | 35 |
| Pitch diameter d (mm) | 85.05 | 248.85 |
| Tip circle diameter df (mm) | 92.17 | 249.88 |
| Normal pressure angle α (°) | 20 | 20 |
Three quantities were tracked: the relative displacement along the line of action, the normal dynamic mesh force, and the radial vibration of the pinion. All three grow with the eccentricity. The mechanism is geometric rather than frictional — the eccentricity acts as a cam, displacing the contact from its theoretical position and thereby enlarging the relative displacement along the normal; that displacement, multiplied by the time-varying stiffness, produces a larger dynamic mesh force, which in turn feeds back into the vibration. The response remains periodic, but its waveform is shaped jointly by the mesh force, the centrifugal force and the torsional inertia term, so the peak amplitude is not simply proportional to e. Even the smallest offset examined, 50 μm, is enough to shift the vibration level measurably, which is the practical justification for a damper that can be retrofitted without touching the tooth geometry.
A shaft-mounted sector-cavity damper
For a rotating gear the particle damper has a property that stationary applications do not share: the centrifugal acceleration, which at a few thousand revolutions per minute is orders of magnitude larger than gravity, presses the particles against the outer wall and converts the filling from a loose granular gas into a compacted bed sliding on the wall. The equivalent stiffness and damping of such a bed vary with time and with speed, the number of particles taking part in energy exchange fluctuates, and individual collisions are unpredictable. There is no reliable closed-form description, so the parameter set has to be identified experimentally.
The hardware adopted is a disc split into two symmetric halves with six equal sector cavities, sized so that no cavity becomes a dead zone where particles settle, and shaped to maximize the frequency of particle–particle and particle–wall contact. Each cavity shell carries a threaded hole for filling and sealing; the flange carries two through-holes so that the halves can be bolted together around the shaft. The body is additively manufactured in PLA thermoplastic and weighs 243.5 g. Its symmetry and its small mass relative to the shaft keep the static balance of the rotor effectively unchanged, and each cavity is loaded with the same particle mass weighed on a precision balance so that no residual unbalance is introduced. As speed rises, the charge redistributes along the inner circumferential surface of the cavity under centrifugal load.
Test bench and performance metrics
The test bench comprised a drive motor, the spiral bevel gear pair, a load motor, a variable-frequency drive, the damper, and the instrumentation chain. Two SJ4-1 high-precision eddy-current sensors, held by magnetic-base stands on a cast-iron platform, measured the radial and axial displacement of the pinion; the signals were acquired by a DH8302 dynamic signal analyser and processed in DHDAS software. The damper was fitted to the pinion shaft rather than the gear shaft, since the pinion runs faster and a higher surface speed of the particle bed was expected to dissipate more energy.
Two dimensionless indices summarize the outcome. The energy dissipation coefficient compares the radial kinetic energy of the shaft with and without the damper,
$$H = \frac{E_0 – E_z}{E_0},\qquad E_0 = \frac{1}{2}m_1 v_0^2,\qquad E_z = \frac{1}{2}m_1 v_z^2$$
where v0 is the mean radial vibration velocity without the damper and vz the value with it. The relative displacement reduction rate compares peak-to-peak amplitudes,
$$S = \frac{A_0 – A_z}{A_0}$$
with A0 and Az the peak-to-peak radial displacement of the pinion without and with the damper. A larger H or S means stronger suppression.
Filling ratio, particle size and material
The three parameters that govern performance were swept one at a time over an input speed range of 600–1500 r/min. Filling ratio was varied from 30% to 90% with 4 mm steel spheres; sphere diameter from 2 mm to 6 mm at 80% filling; and material — carbon steel, aluminium and resin — at 4 mm and 80% filling. Both H and S follow the same non-monotonic pattern in every sweep, and both rise monotonically with rotational speed.
| Parameter | Range examined | Fixed conditions | Best setting | Observed behaviour |
|---|---|---|---|---|
| Filling ratio | 30–90% | steel, d = 4 mm, 600–1500 r/min | 80% | H and S peak at 80%, lower on both sides |
| Particle diameter | 2–6 mm | steel, 80% filling, 600–1500 r/min | 4 mm | H and S peak at 4 mm, lower on both sides |
| Particle material | carbon steel, aluminium, resin | d = 4 mm, 80% filling, 600–1500 r/min | carbon steel | steel > aluminium > resin at every speed |
| Rotational speed | 600–1500 r/min | any of the above | highest tested | dissipation increases with speed |
The physical reading is straightforward. A sparse filling gives too few contacts to absorb energy, while an over-filled cavity forces the particles into a tightly packed bed in which relative motion, and hence the number of effective collisions, falls again; there is an optimum near 80%. Small spheres raise the collision count but carry little momentum each, so the energy transferred per impact is small; large spheres reduce the collision frequency and leave the bed insufficiently agitated. The 4 mm optimum balances the two effects. Material acts mainly through density and restitution: carbon steel delivers the largest momentum and the highest contact forces for a given impact velocity, resin the smallest, with aluminium in between.
Representative time histories quantify the effect. At 1000 r/min with an optimally configured steel charge, the peak-to-peak radial displacement of the pinion fell from 0.278 mm to 0.227 mm, a reduction of 17.5%. At 1200 r/min the corresponding pair was 0.304 mm and 0.257 mm (16.4%), and at 1500 r/min 0.314 mm and 0.255 mm (18.7%). The speed dependence is consistent with the collision-frequency argument: a faster rotor agitates the bed more violently and pushes more of the charge into the sliding, dissipative regime.
| Test condition | Input speed (r/min) | Without damper (mm) | With damper (mm) | Reduction S (%) |
|---|---|---|---|---|
| Balanced rotor | 1000 | 0.278 | 0.227 | 17.5 |
| Balanced rotor | 1200 | 0.304 | 0.257 | 16.4 |
| Balanced rotor | 1500 | 0.314 | 0.255 | 18.7 |
| 37.5 g magnet at 50 mm radius on the gear | 1000 | 0.308 | 0.245 | 20.4 |
Suppression with a deliberately imposed eccentric mass
To confirm that the damper still works when the excitation it was designed for is actually present, an eccentric mass was added to the transmission by the mass-addition method: a magnet of 37.5 g was attached 50 mm from the axis of the driven gear, which corresponds to an unbalance of 1875 g·mm, or 1.875 × 10−3 kg·m. The driven gear was chosen because the pinion end face is too small to hold the magnet securely, and because the eccentricity of either member excites the pair once the two wheels are in mesh. With the pinion shaft running at 1000 r/min and the damper fitted with its optimal charge, the peak-to-peak radial displacement of the pinion fell from 0.308 mm to 0.245 mm, a reduction of 20.4% — a larger relative gain than the 17.5% obtained at the same speed on the balanced rotor, which is what a dissipation mechanism that scales with vibration amplitude would predict.
Design implications and open questions
Three conclusions follow from the model and the measurements. First, eccentricity is a first-order excitation for a spiral bevel pair: the relative displacement along the line of action, the dynamic mesh force and the radial vibration all increase with the offset, and the response stays periodic but not proportionally so, because mesh stiffness, centrifugal force and torsional inertia act together. Second, a disc damper with equal sector cavities, bolted to the gear shaft, is a workable way of introducing particle damping into a centrifugal field without modifying the gear body — the feature that distinguishes it from web-hole and damping-ring solutions, which require cavities, seals and grooves cut into the gear itself. Third, the dissipation capacity grows with shaft speed, so the device is most useful exactly where gear whine and dynamic tooth load are most severe.
The numbers also mark the limits. A 17–20% reduction in peak-to-peak displacement is useful but modest, and it was obtained on a bench with a single damper geometry, a single charge mass and one pair of gears; the optimum at 80% filling and 4 mm steel spheres is specific to these six cavities and this 243.5 g body, since fill ratio and particle size are normalized against cavity volume and wall speed respectively, and a differently proportioned damper will need its own sweep. The added 243.5 g on the pinion shaft must be accounted for in the rotordynamics of the actual machine, even though the symmetry of the body keeps static balance intact. Nothing in the tests addresses long-term behaviour: PLA is sensitive to temperature and to lubricating oil, and a production version would need a metal or engineering-polymer body, together with evidence that the charge does not compact irreversibly or wear through the cavity wall over service life. Extending the speed range well beyond 1500 r/min, where centrifugal acceleration completely dominates and the particle bed may simply lock against the wall, is the next question worth answering, because that is the regime in which a spiral bevel stage in a real driveline spends most of its life.