The pursuit of enhanced durability, fatigue resistance, and operational reliability in modern machinery places stringent demands on the surface integrity of critical components, among which spur and pinion gears are paramount. Mass finishing processes, particularly spindle-type barrel finishing, offer a versatile and cost-effective solution for improving surface topography, inducing beneficial residual stresses, and achieving isotropic surface textures on complex geometries like spur and pinion gears. However, the underlying mechanisms at the granular level, especially the particle action behavior on the intricate tooth flank surfaces, remain inadequately understood, hindering predictive process optimization. This work employs a Discrete Element Method (DEM) simulation framework to fundamentally investigate the dynamic interactions between abrasive media and the tooth surfaces of a spur and pinion gear during spindle finishing, elucidating the effects of key process parameters.
Fundamentals of Spindle Finishing and Numerical Modeling
The spindle finishing setup consists of a rotating cylindrical drum partially filled with abrasive media, into which the workpiece—a spur and pinion gear mounted on a shaft—is immersed. Both the drum and the workpiece shaft rotate concurrently, generating complex relative motion and interaction forces between the media and the gear tooth flanks. The process aims to uniformly treat surfaces through a combination of micro-cutting, peening, and plastic deformation induced by countless particle-workpiece contacts.

To model this stochastic process, the DEM approach is adopted, wherein the motion of each individual abrasive particle is tracked by solving Newton’s second law of motion. The contact forces between particles, and between particles and boundaries (drum, gear), are calculated using a Hertz-Mindlin contact model. For wear prediction on the gear surface, the Archard adhesive wear model is integrated, which posits that the wear volume is proportional to the normal contact force and the sliding distance. The incremental wear depth \(\Delta h\) on a surface element is given by:
$$ \Delta h = \frac{K P v}{H} \Delta t $$
where \(K\) is a dimensionless wear coefficient, \(P\) is the normal pressure at the contact, \(v\) is the relative sliding velocity, \(H\) is the surface hardness of the workpiece, and \(\Delta t\) is the time increment. This equation underpins the analysis, as the product \(P \cdot v\) is a key indicator of the local material removal rate on the flank of the spur and pinion gear.
The simulation model comprises a steel drum, spherical brown alumina abrasive particles (3 mm diameter), and a standard spur and pinion gear model. The material properties and contact parameters are summarized in the following tables.
| Material | Density, \(\rho\) (kg/m³) | Poisson’s Ratio, \(\varepsilon\) | Shear Modulus, G (MPa) |
|---|---|---|---|
| Drum (Steel) | 7850 | 0.300 | 7.94e3 |
| Particle (Alumina) | 2675 | 0.360 | 1.26e3 |
| Spur and pinion gear (40Cr) | 7870 | 0.277 | 8.08e3 |
| Interaction Pair | Coefficient of Restitution | Static Friction Coefficient | Rolling Friction Coefficient |
|---|---|---|---|
| Particle-Drum | 0.50 | 0.35 | 0.10 |
| Particle-Spur and pinion gear | 0.43 | 0.36 | 0.10 |
| Particle-Particle | 0.46 | 0.39 | 0.10 |
A parametric study is designed to isolate the influence of two critical process parameters on the particle action behavior on the spur and pinion gear tooth flank:
- Workpiece Immersion Depth (\(h_1\)): The vertical distance from the gear’s top face to the stationary media free surface. This primarily affects the hydrostatic pressure and the number of particles actively contacting the gear.
- Rotational Speed (\(n_1, n_2\)): The drum speed (\(n_1\)) and gear shaft speed (\(n_2\)), maintained at a fixed ratio (\(n_1:n_2 = 5:4\)). This primarily governs the kinetic energy and relative velocity of the media.
The parameter sets for simulation are: \(h_1 = [80, 110, 140]\) mm and \(n_1 = [12, 21, 30]\) rpm.
Analysis of Granular Flow Field and Particle Action Cycle
The introduction of the spur and pinion gear into the flowing media significantly disrupts the otherwise parabolic free surface. A pronounced pile-up of particles forms on the leading face of the gear (in the direction of drum rotation), while a low-pressure “wake” region devoid of particles trails behind it. This asymmetric pressure distribution is a fundamental characteristic influencing the processing.
Velocity vector analysis reveals the complex, three-dimensional flow path around the spur and pinion gear. Particles impacting the leading face are deflected upward, cascading over the gear’s top. A significant portion also flows laterally around the gear sides. Crucially, the flow within the tooth space itself exhibits a distinct, periodic cycle synchronized with the gear’s rotation. This cycle can be segmented into three sequential phases for a given tooth flank:
- Filling Phase: As the tooth rotates out of the wake, media rapidly flows into the tooth space from the top and the tip. Particle velocity relative to the flank is high and directed primarily into the root.
- Stable Filling (Dwell) Phase: The tooth space is fully packed with media and is positioned within the high-pressure pile-up region. Particles are subjected to intense confinement and shear, leading to sustained, high normal contact forces against the flank. The relative particle motion transitions to an upward sliding direction along the flank contour.
- Emptying Phase: As the gear continues to rotate, the tooth moves out of the pile-up region. Assisted by centrifugal force and gravity, the media rapidly evacuates the tooth space from the tip and bottom, leading to a sharp drop in contact forces.
This cyclical action is the core material removal mechanism for each flank of the spur and pinion gear. The dwell phase is identified as the dominant period for surface modification, where both contact forces and interaction durations are maximized.
Influence of Parameters on Contact Forces and Relative Velocity
The DEM simulation allows for the precise extraction of time-averaged normal contact force (\(F\)) and particle relative sliding velocity (\(v\)) on discrete patches of the spur and pinion gear tooth flank during the stable dwell phase. The effects of immersion depth and rotational speed are analyzed separately.
Effect of Immersion Depth (\(h_1\)): Increasing the immersion depth significantly amplifies the normal contact force on the tooth flank. This is attributed to the increased hydrostatic head and the greater mass of media contributing to the pile-up pressure on the gear. Quantitatively, a 75% increase in \(h_1\) (from 80 mm to 140 mm) results in an approximate 76% increase in the average normal contact force \(F\) on the flank. In contrast, the relative sliding velocity \(v\) shows only a minimal increase (≈4%) with depth, as it is more governed by kinematic factors than static pressure. The relationship can be summarized as:
$$ \Delta F \propto \Delta h_1, \quad \Delta v \approx f(\Delta h_1) $$
Effect of Rotational Speed (\(n_1, n_2\)): Increasing the rotational speed has a profound effect on the relative kinetic energy in the system. The relative sliding velocity \(v\) between the media and the spur and pinion gear flank increases nearly proportionally with speed. A 150% increase in \(n_1\) (from 12 rpm to 30 rpm) leads to a 148% increase in \(v\). The normal contact force \(F\), however, experiences a much smaller increase (≈18%) with speed, as the inertial impacts contribute less to the sustained pressure than the static/media weight loading does. This relationship is captured as:
$$ \Delta v \propto \Delta n_1, \quad \Delta F \approx f(\Delta n_1) $$
Spatial Non-Uniformity of Processing on the Spur and Pinion Gear
A critical finding from the simulation is the inherent non-uniformity of particle action across the tooth flank of the spur and pinion gear. This manifests in two primary dimensions:
- Profile Direction (Root-to-Tip): Due to the semi-enclosed geometry of the tooth space, particle motion and force transmission are hindered near the root. Consequently, both contact force \(F\) and relative velocity \(v\) are lower near the root and increase progressively towards the tip. Neither immersion depth nor speed significantly alters this inherent gradient along the profile.
- Axial Direction (Top-to-Bottom Face): A notable asymmetry exists between the upper and lower halves of the tooth flank. The upper flank experiences significantly greater action. On average, the normal contact force on the upper flank is 1.5 to 1.8 times that on the lower flank. Similarly, the relative sliding velocity is 1.35 to 1.45 times higher on the upper flank. This is due to the dominant media flow path being deflected over the top of the gear.
Importantly, the axial non-uniformity can be mitigated by process parameter selection. Increasing the immersion depth \(h_1\) effectively submerges a greater portion of the gear in the active media flow, reducing the disparity between the upper and lower regions. At a shallow depth (80 mm), the simulated force difference is maximal, whereas at a deep immersion (140 mm), the forces and velocities across the axial length become more uniform.
Experimental Validation and Surface Outcome Correlation
To validate the DEM findings, experimental strain measurements were conducted on a spur and pinion gear during actual spindle finishing under matching parameters. Strain gauges mounted at strategic locations (upper flank, top face, bottom face) recorded the dynamic stress signals induced by media contact.
The experimental stress data corroborated the simulation trends conclusively:
- Stress signals on the tooth flank showed clear periodicity, matching the predicted particle action cycle.
- Stresses increased modestly (17-28%) with increasing rotational speed, mirroring the simulated 18% increase in contact force.
- Stresses increased substantially (35-59%) with increasing immersion depth, aligning with the simulated 76% increase in force.
Furthermore, surface roughness (\(R_a\)) measurements were taken after processing. The results confirmed the parametric influence and the effect on uniformity:
- Increasing rotational speed \(n_1\) significantly improved roughness reduction across the gear, with the most pronounced effect on the bottom region.
- Increasing immersion depth \(h_1\) dramatically improved both the magnitude and the axial uniformity of roughness reduction. At \(h_1 = 80\) mm, the roughness reduction rate varied significantly from the top to the bottom of the flank (17% vs. 36%). At \(h_1 = 140\) mm, the rate became highly uniform (62% vs. 55%), confirming the simulation prediction that deep immersion homogenizes axial processing.
Conclusions
This integrated DEM simulation and experimental study provides a granular-level understanding of the particle action behavior during spindle finishing of spur and pinion gears. The key conclusions are:
- The processing action on an individual tooth flank of a spur and pinion gear is inherently periodic, comprising distinct filling, stable dwell, and emptying phases. The stable dwell phase is responsible for the majority of surface modification.
- Process parameters exert distinct and primary influences: Immersion depth (\(h_1\)) is the dominant control parameter for increasing the normal contact force on the flank. Rotational speed (\(n_1, n_2\)) is the dominant control parameter for increasing the relative sliding velocity between particles and the flank.
- Spatial non-uniformity is intrinsic to the process. The upper flank of a spur and pinion gear consistently experiences more intense action (1.5-1.8x higher force, 1.35-1.45x higher velocity) than the lower flank. A gradient also exists from the root to the tip of the tooth.
- Axial processing uniformity on the spur and pinion gear can be effectively improved by increasing the immersion depth. A greater depth submerges more of the gear’s axial length into the active media flow, equalizing the action between the upper and lower regions of the flank, as confirmed by surface roughness measurements. The profile-direction gradient, however, is less sensitive to parameter changes.
These insights establish a fundamental framework for modeling and predictively optimizing the spindle finishing process for spur and pinion gears, moving beyond empirical trial-and-error towards science-based process design for achieving superior and uniform surface integrity.
