Particle Damping for Eccentric-Mass Vibration in Spiral Bevel Gear Drives

Spiral bevel gears carry high torque at right angles with smooth, quiet motion, which is why they dominate final drives in automotive axles, aerospace gearboxes, power transmission equipment and drilling machinery. Their load capacity comes at a price: mesh stiffness varies through the contact cycle, the contact ratio is large but not constant, and any offset between the geometric centre and the rotational axis turns the gear into a rotating unbalance source. Eccentricity arising from mounting error, bore clearance or manufacturing scatter generates centrifugal and inertial forces in the plane of rotation, and those forces add directly to the dynamic mesh force. The result is a periodic radial vibration whose amplitude scales with the eccentricity rather than with the transmitted load.

Remedies divide into active and passive families. Active approaches attack the source—grinding accuracy, macro-geometry optimisation, profile modification—and are effective but computationally heavy, slow to iterate and expensive to implement on a production gear set. Passive approaches absorb or dissipate vibration energy with an add-on device: squeeze-film dampers, damping rings fitted to the gear web, or wire-rope elements. Particle damping belongs to this second family and has a distinct advantage: a cavity filled with loose granules behaves as a strongly nonlinear, broadband absorber, because energy leaves the structure through inelastic collisions and sliding friction between particles and between particles and the cavity wall.

Most gear applications of particle damping place the device on the gearbox housing, where the vibration field is steady and gravity sets the particle bed. Mounting inside the rotating gear itself—drilling pockets in the web, sealing them with cover plates, or fitting split damping rings into machined grooves—brings the particles close to the source but weakens the web, adds sealing hardware and raises manufacturing cost. A damper that clamps onto the gear shaft avoids all three penalties, at the price of operating in a centrifugal field where the particle bed is pressed against the outer wall and the effective stiffness and damping of the device become time-varying and difficult to describe analytically. This article examines the mechanism of eccentricity-induced vibration in a spiral bevel pair, describes a shaft-mounted sector-cavity damper, and reviews the experimental evidence on how filling ratio, granule size and granule material govern energy dissipation in the rotating frame.

Dynamic model of an eccentric spiral bevel pair

The pair is represented as a lumped-parameter system in a Cartesian frame whose origin is the intersection of the two shaft axes. The mesh is a spring–damper element with time-varying stiffness k(t) and mean mesh damping cm; bearings are isotropic elastic supports with stiffness kij and damping cij (i = x, y, z; j = 1, 2). Tooth friction and static transmission error are neglected, leaving each gear with three translational coordinates plus rotation about its own axis, giving eight degrees of freedom for the pair:

$$
\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}
$$

Because a spiral bevel pair engages progressively over a large contact ratio, the mesh stiffness shows no step discontinuity of the spur-gear kind; it fluctuates about a mean value and is approximated by a first-order harmonic:

$$k(t) = k_m + k_{kl}\cos\!\left(l\omega_n t + \psi_{kl}\right), \qquad l = 1,\ \psi_{kl} = 0$$

Eccentricity enters through the relative displacement along the line of action. With the pinion centre offset by e1 and the wheel by e2, the projection of the mesh-point motion acquires additional cosine and sine terms:

$$
\begin{aligned}
\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)
\end{aligned}
$$

$$a_q = \cos\delta_q\sin\alpha_n, \qquad 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, \qquad r_{pq} = r_q / \cos\alpha_n$$

where q = 1, 2 labels pinion and wheel, δq is the pitch cone angle, βq the mid-face spiral angle and αn the normal pressure angle. The normal dynamic mesh force and its 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 offset also produces a centrifugal force Fl in the plane of rotation and an inertial force Fg driven by torsional acceleration, both of which project onto the radial directions of the respective gear:

$$
\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}
$$

$$
\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}
$$

Table 1. Geometry of the spiral bevel pair used in the model and tests
Parameter Pinion Wheel
Number of teeth z 27 79
Hand of spiral left right
Face width b (mm) 22.8 22.8
Pitch cone angle δ (°) 18.87 71.13
Shaft angle (°) 90 90
Transverse module (mm) 3.15 3.15
Mid-face spiral angle β (°) 35 35
Pitch diameter (mm) 85.05 248.85
Tip diameter (mm) 92.17 249.88
Normal pressure angle αn (°) 20 20

What eccentricity does to the response

Numerical integration of the eight equations by a Runge–Kutta scheme gives the response of the pair at 60 kW input, 3 000 r/min pinion speed and 200 N·m load torque, with identical eccentricities assigned to pinion and wheel at 50, 100 and 200 µm. Three quantities respond monotonically to the offset. The relative displacement along the mesh normal grows with eccentricity, because the offset makes the actual contact point depart from its theoretical position—a cam-like effect that lengthens the effective line of action. The normal dynamic mesh force grows in step, since the enlarged relative displacement is multiplied by the time-varying stiffness. The pinion radial vibration displacement grows as well, and it fluctuates periodically rather than randomly, because the mesh force, the centrifugal force and the inertial force all vary at once-per-revolution and mesh-order frequencies.

The practical reading of these trends is that eccentricity does not merely add a once-per-revolution component to an otherwise acceptable spectrum; it raises the dynamic mesh force itself, so the tooth contact and the bearing loads are affected simultaneously. Suppression therefore has to remove energy from the shaft rather than filter a single harmonic.

A shaft-mounted sector-cavity damper

The device is a disc split into two symmetric halves, each carrying six equal sector cavities, clamped to the gear shaft through two through-holes in the flange. Nothing on the gear body is modified, so web strength and tooth geometry are untouched, and the damper can be retrofitted. Each cavity has a threaded port for filling and sealing. The six-fold symmetry distributes mass evenly, limits particle piling in any one pocket under gravity, and increases the number of particle–wall and particle–particle impacts per revolution.

The housing was produced by 3D printing in PLA thermoplastic and weighs 243.5 g. Because the geometry is symmetric and the mass is small relative to the shaft assembly, static balance is essentially unaffected; during filling, each cavity is charged with an equal mass weighed on a precision balance to avoid introducing a new unbalance. As shaft speed rises, the granules are thrown outward and spread over the inner circumferential surface of the housing, forming a bed that is continuously re-excited by the shaft vibration. The damper was installed on the pinion shaft, which runs faster than the wheel shaft and therefore experiences the higher centrifugal acceleration and the larger collision frequency.

Spiral bevel gear with curved teeth on a tapered hub

Test rig and dissipation metrics

The experimental system comprised a drive motor, the spiral bevel gearbox, a loading motor, a frequency converter, the damper, and the instrumentation chain. Two SJ4-1 eddy-current displacement sensors, held on magnetic stands fixed to a cast-iron platform, measured the pinion radial and axial vibration; signals were acquired by a DH8302 dynamic signal analyser and processed in DHDAS software. Two dimensionless indices summarise the outcome. The energy dissipation coefficient compares the mean kinetic energy of the gear shaft with and without the damper, and the relative displacement reduction compares peak-to-peak vibration amplitudes:

$$H = \frac{E_0 – E_z}{E_0}, \qquad E = \frac{1}{2}m_1 v^2$$

$$S = \frac{A_0 – A_z}{A_0}$$

Here E0 and A0 are the kinetic energy and peak-to-peak amplitude of the undamped shaft, Ez and Az the values with the particle damper installed, and v the mean radial vibration velocity. A larger H or S means more energy diverted into the granular bed.

Selecting filling ratio, granule size and material

The three parameters were swept one at a time over a pinion speed range of 600–1 500 r/min. A low filling ratio leaves too few contacts per unit time; an excessive one compacts the bed so that relative motion between granules is suppressed. Both limits appear in the measurements as a rise and fall of H and S with filling ratio, with a maximum at 80% for steel spheres of 4 mm diameter. Granule size behaves the same way: small spheres pack densely and collide often but carry little mass per impact, while large spheres collide rarely; the optimum was 4 mm over the 2–6 mm range examined. Material separates the three candidates unambiguously, with carbon steel giving the highest dissipation coefficient, aluminium intermediate and resin the lowest at every speed tested—a ranking consistent with the loss mechanisms being dominated by friction and inelastic impact, both of which scale with granule mass and contact stiffness.

Table 2. Parameter sweep for the shaft-mounted damper (pinion speed 600–1 500 r/min)
Design factor Range examined Optimum Observed behaviour
Filling ratio 30–90% 80% H and S rise then fall with filling ratio; both increase with speed at fixed filling ratio
Granule diameter 2–6 mm 4 mm Maximum H and S at 4 mm; both increase with speed at fixed diameter
Granule material carbon steel, aluminium, resin carbon steel Steel > aluminium > resin at every speed; dissipation increases with speed for a given material

The speed dependence is a consistent feature of the data: at fixed filling ratio, granule size or material, both indices grow as the pinion accelerates from 600 to 1 500 r/min, indicating that the collision frequency between granules and between granules and the cavity wall rises faster than any stiffening of the centrifugally compacted bed. For a gearbox operating at high speed, this is the favourable direction.

Table 3. Peak-to-peak pinion radial displacement with and without the damper (steel granules, 4 mm, 80% filling)
Test case Pinion speed (r/min) A0 (mm) Az (mm) S (%)
Filling-ratio series, 80% 1 000 0.278 0.227 17.5
Granule-size series, 4 mm 1 200 0.304 0.257 16.4
Material series, carbon steel 1 500 0.314 0.255 18.7
Deliberate eccentric mass, 37.5 g at 50 mm 1 000 0.308 0.245 20.4

Validation on a pair with a deliberate eccentric mass

The optimised configuration—carbon steel spheres, 4 mm, 80% filling—was then tested on a pair carrying a known unbalance. Because the pinion face is too small to hold a magnet securely, a 37.5 g magnet was attached to the face of the 79-tooth wheel at a radius of 50 mm from its axis, producing an unbalance of 1 875 g·mm that excites the whole mesh as the pair rotates. At an input speed of 1 000 r/min the wheel turns at about 342 r/min, so the resulting centrifugal force is only of the order of 2.4 N in rigid-body rotation; the measured effect is nevertheless clear, with the pinion radial peak-to-peak amplitude falling from 0.308 mm to 0.245 mm, a reduction of 20.4%. The reduction exceeds the 17.5% obtained at the same speed without the added unbalance, which suggests the granular bed draws proportionally more energy from a response containing a strong once-per-revolution component.

Interpretation and scope

Two mechanisms explain the observed behaviour. First, the granules act as a large number of small, dry-friction contacts that convert shaft kinetic energy into heat through impact and sliding; because the contacts are numerous and their orientations random, dissipation is broadband rather than tuned to one frequency. Second, centrifugation organises the bed against the outer wall, so the device self-aligns with the vibration direction instead of relying on gravity, which is why performance improves with speed rather than degrading.

The evidence also delimits the claim. Equivalent stiffness and damping of a granular bed are time-varying and cannot be written in closed form, so the damper was characterised experimentally and inserted into the system model only through the measured reduction indices; no coupled simulation of particles plus gears is offered. Validation under the deliberate eccentric mass was performed at a single speed, 1 000 r/min, and with the damper in one location—the pinion shaft. The remaining questions are quantitative: how the reduction scales between 1 000 and 1 500 r/min under sustained unbalance, whether a second damper on the wheel shaft adds measurable benefit, and how the PLA housing tolerates long-term service in an oil-lubricated gearbox. Within those limits, the shaft-mounted sector-cavity damper is a low-cost retrofit that reduces eccentric-mass-induced vibration of spiral bevel drives by roughly one-fifth without altering the gear body or its load path.

Scroll to Top