Spiral bevel gears transmit motion between intersecting shafts and are widely used in aerospace, automotive, and industrial drives. Their tooth flanks are spatial curved surfaces generated on a conical pitch surface, and the contact path is diagonal and sensitive to small changes in topography. Surface roughness, grinding marks, burrs, and edge condition influence the lubricant film, contact-pressure distribution, vibration, noise, and fatigue resistance. Grinding can produce accurate geometry, but it often leaves directional marks and micro-burrs; shot peening improves residual stress but does not remove those features. Mass finishing is therefore attractive as a terminal process because it can reach tooth roots and flanks, remove burrs, and modify the surface texture without changing the macro-geometry.
Parallel-spindle barrel finishing is a mass-finishing variant in which the workpiece axis is parallel to the drum axis. The workpiece and drum rotate in opposite directions while the gear is immersed in a flowing abrasive medium. Material removal is driven by the relative velocity between the abrasive and the tooth surface, together with the local confinement of the media in the tooth space. For spiral bevel gears, the curved flanks and tapered disc geometry make the flow field strongly non-uniform, so a purely kinematic estimate is not sufficient to predict where material is removed. The process therefore needs to be analyzed at two levels: a global motion model for cutting speed, and a discrete-element model for particle-scale wear distribution.

Kinematics of parallel-spindle barrel finishing
Consider a point A on the workpiece at radius r from the workpiece axis. The workpiece axis is offset from the drum axis by distance R. In a frame in which the drum is fixed, the workpiece axis orbits around the drum axis at angular velocity Ω, while the workpiece spins about its own axis at angular velocity ω in the opposite direction. The position of point A can be written as
$$
\begin{aligned}
x &= r\sin(\theta_1-\theta_2)-R\sin\theta_1,\\
y &= -r\cos(\theta_1-\theta_2)+R\cos\theta_1,
\end{aligned}
$$
where \(\theta_1=\Omega t\) and \(\theta_2=\omega t\). Differentiating and taking the magnitude of the velocity vector gives the relative cutting speed:
$$
v=\left[r^{2}(\Omega-\omega)^{2}+\Omega^{2}R^{2}+2Rr\,\Omega(\Omega-\omega)\cos\theta\right]^{1/2},
$$
where θ is the instantaneous phase angle of the point in the relative motion. The expression shows that the cutting speed increases with workpiece radius r, drum speed Ω, and the speed difference between the drum and the workpiece. For a representative set of conditions—drum speed \(N=45\ \mathrm{r/min}\), workpiece speed \(n=-45\) to \(-135\ \mathrm{r/min}\), center distance \(R=0.41\ \mathrm{m}\), and workpiece radius \(r=0.05\) to \(0.1\ \mathrm{m}\)—the calculated speed rises as either the radius or the spindle speed increases. The outer region of the gear and higher spindle speeds therefore experience greater cutting energy.
This kinematic result provides a first-order control variable, but it treats the abrasive medium as a continuum moving with the drum. It does not resolve particle collisions, tooth-space shadowing, or the clustering of media on one flank. Those effects determine the local wear field and must be captured with a discrete-element simulation.
Discrete-element wear simulation
The finishing process was simulated in EDEM using a Hertz–Mindlin contact model coupled with Archard wear. The drum and disc-type gear geometry were imported and scaled to keep the particle count and time step manageable. The drum rotated at \(-50\ \mathrm{r/min}\) and the workpiece at \(+115\ \mathrm{r/min}\), with the sign convention that opposite rotation corresponds to opposite signs. The simulation used 200 000 spherical abrasive particles with a mean diameter of 2 mm and a normal size distribution. The abrasive was assigned a Poisson’s ratio of 0.25 and a shear modulus of 129 GPa; the gear was assigned a Poisson’s ratio of 0.30 and a shear modulus of 70 GPa. Rotational motion started at 2 s, the total simulated time was 4 s, the time step was 0.02 s, and the grid size was set to \(3R_{\min}\), where \(R_{\min}\) is the smallest particle radius.
| Contact pair | Restitution coefficient | Static friction coefficient | Rolling friction coefficient |
|---|---|---|---|
| Abrasive–abrasive | 0.50 | 0.50 | 0.15 |
| Abrasive–gear | 0.30 | 0.55 | 0.01 |
The simulated wear field was not uniform. Because the disc-type spiral bevel gear has an inverted-cone shape, the small end showed greater wear than the large end. Wear also increased with processing time for every spindle speed tested. When the workpiece speed was set to 115, 75, and 50 r/min, the highest speed produced the largest wear volume. This trend agrees with the kinematic model: a larger speed difference increases the relative velocity between the abrasive and the tooth surface, raising the collision energy available for abrasive wear. The simulation was short and geometrically scaled, so the absolute wear depth should not be treated as a production-life prediction. Its value lies in ranking the process variables and revealing the spatial bias of the wear field.
Experimental conditions and measurement
The finishing trials used a disc-type aviation spiral bevel gear with module \(m=4.4\ \mathrm{mm}\), \(z_2=41\) teeth, and a mean spiral angle \(\beta=35^\circ\). The gear material was 20Cr2Ni4A, carburized and quenched, form-ground, and shot-peened, with a mean surface hardness of 59.2 HRC. The drum diameter was 1260 mm. The abrasive medium consisted of 2 mm corundum-based spherical media, and the fill ratio was 70% of the drum volume. The compound was a mixture of water and polishing liquid at a fixed ratio and fixed flow rate. The gear was mounted on a fixture with its axis parallel to the drum axis; the center distance between the workpiece and drum axes was 410 mm, and the immersion depth was 150 mm. The drum and workpiece rotated in opposite directions, and the forward and reverse rotation times were kept equal so that both flanks were exposed under comparable conditions.
Surface roughness was measured with a contact profilometer capable of measuring curved surfaces. Readings were taken every 6 min at the same location, and three measurements were averaged. A dedicated fixture with a small indexing plate was used to reposition the disc-type gear accurately and reduce measurement scatter.
Surface-quality response
The arithmetic mean roughness \(R_a\) decreased from 0.8 μm to approximately 0.4 μm. Before finishing, the convex flank was rougher than the concave flank. During processing, the concave flank roughness decreased more rapidly. In the narrow tooth space, particles tend to cluster and form vortex-like circulation, which increases the number of abrasive collisions on the concave side and raises the local removal rate. The roughness curves flattened after the major peaks had been removed, indicating that the process should be stopped before excessive edge radiusing occurs. The small end also showed a higher roughness-reduction rate than the large end, consistent with the simulated wear bias.
| Characteristic | Before finishing | After finishing |
|---|---|---|
| Arithmetic mean roughness \(R_a\) | 0.8 μm | ~0.4 μm |
| Edges and burrs | Burrs at tooth tip and outer diameter | Burrs removed; edges rounded |
| Texture | Directional grinding marks | Marks reduced; more isotropic |
After finishing, burrs at the tooth tip and outer diameter were removed, edges were rounded, and grinding marks on the flanks and in the tooth roots were reduced. The surface appeared brighter. At 140× magnification under a digital three-dimensional video microscope, the finished surface showed a more isotropic texture than the ground surface. This change is relevant to running-in and lubrication, because an isotropic texture can distribute contact stress and retain lubricant more uniformly than a strongly directional grinding texture.
Mechanistic interpretation and engineering implications
The results support a two-stage view of barrel finishing. In the first stage, abrasive impacts remove burrs and flatten the highest asperities, so roughness drops quickly. In the second stage, the surface is smoother and the removal rate falls; continued processing mainly rounds edges and may begin to alter the contact pattern. For spiral bevel gears, the contact pattern is sensitive to edge geometry, so finishing time should be optimized against roughness, edge radius, and contact-pattern requirements rather than roughness alone.
Cutting speed is a useful but incomplete control variable. The kinematic model predicts that increasing the spindle speed or the workpiece radius raises the relative speed. The discrete-element results show that the same variables also change the spatial distribution of wear. In production, fixture orientation, immersion depth, media size, and cycle time can be adjusted to compensate for the small-end and concave-flank bias. Media shape and size distribution are also important: spherical media were used in these trials, whereas angular or mixed media would change cutting efficiency and surface texture.
The evidence here is limited to roughness and topography. Residual stress, hardness, and contact fatigue were not measured, and the trials covered one gear geometry and one media specification. A complete assessment would require contact-pattern checks, residual-stress measurements, and bench or FZG-type fatigue testing. The simulation also used a short physical time and scaled geometry; longer runs with calibrated wear coefficients would be needed for quantitative prediction. Within those limits, the combination of kinematics, DEM wear simulation, and finishing trials provides a coherent basis for applying parallel-spindle barrel finishing to spiral bevel gears.
Conclusions
Parallel-spindle barrel finishing reduced the roughness of a disc-type spiral bevel gear from \(R_a=0.8\ \mu\mathrm{m}\) to about \(0.4\ \mu\mathrm{m}\), removed burrs at the tooth tip and outer diameter, rounded edges, and produced a more isotropic surface texture. The kinematic model identifies workpiece radius, drum speed, and spindle speed as first-order variables controlling cutting speed. DEM simulations using the Hertz–Mindlin contact model with Archard wear reproduce the main trends: wear increases with spindle speed and is higher at the small end of the inverted-cone gear. The experiments confirm a faster initial roughness reduction on the concave flank, followed by a plateau as the major asperities are removed. Uniformity and edge control remain the key engineering issues for industrial application.
