The continuous advancement of industrial intelligence has led to a growing demand for high-precision, small-modulus gears, placing increasingly stringent requirements on the production efficiency and machining accuracy of high-speed gear shaping machines, which are specialized equipment for manufacturing such components.
A critical limitation of traditional gear shaping machine spindles operating under high stroking conditions is the occurrence of dry friction and seizure phenomena, severely restricting their operational efficiency and broader application. The adoption of hydrostatic spindles, renowned for their high stiffness, exceptional stability, low friction, and minimal wear, presents a viable solution to significantly enhance the productivity of these machines. However, even with hydrostatic technology, persistent issues of localized friction and wear can emerge during sustained high-speed stroking operations. The hydrostatic spindle is the core component determining machining quality. Consequently, optimizing its structure to improve performance has become a key research focus. This study specifically investigates the influence of incorporating an axial groove structure on the working performance of a hydrostatic spindle for high-speed gear shaping applications, based on theoretical analysis and numerical simulation under various operating conditions.

The primary motion of a hydrostatic spindle in a gear shaping machine is a high-speed reciprocating stroking motion combined with a slow rotary indexing motion. Given the substantial difference in linear speeds between these two motions, the slow rotation can be neglected for analysis of the dynamic pressure effects, allowing the system to be approximated as a high-speed cylindrical hydrostatic guideway. During the cutting process, the spindle shaft is subjected to intermittent radial loads, the magnitude of which varies with different radial feed movements. The hydrostatic oil film provides the necessary supporting force to the shaft and simultaneously acts as a fluid medium transmitting this radial load to the stationary sleeve, potentially causing deformation.
Theoretical Foundation and Boundary Conditions
The load-carrying capacity of the hydrostatic oil film is a paramount parameter for spindle performance. When a radial load is applied during machining, the shaft centerline displaces relative to the sleeve, creating an eccentricity and resulting in a non-uniform oil film thickness. The pressurized section of the film generates a supporting reaction force, preventing direct metallic contact between the shaft and the sleeve.
The pressure within the hydrostatic pocket can be described by the equation for a capillary-compensated bearing:
$$p=\frac{p_s}{1+\lambda_0(h/h_0)^3}$$
where \(p_s\) is the supply pressure, \(h_0\) is the initial (design) film thickness, \(h\) is the minimum film thickness under load, and \(\lambda_0\) is the initial resistance ratio.
The eccentricity ratio of the hydrostatic spindle is defined as:
$$\varepsilon = \frac{e}{h_0} = 1 – \frac{h}{h_0}$$
where \(e\) is the relative radial displacement between the shaft and sleeve under load.
Thus, the load-carrying capacity \(F\) of the oil film is:
$$F = p A_e = \frac{p_s A_e}{1+\lambda_0(h/h_0)^3} = \frac{p_s A_e}{1+\lambda_0(1-\varepsilon)^3}$$
where \(A_e\) is the effective bearing area of the hydrostatic pocket.
The stiffness \(j\) of the oil film, which is the rate of change of force with respect to film thickness, is derived as:
$$j = \frac{\partial F}{-\partial h} = \frac{p_s A_e [3\lambda_0(1-\varepsilon)^2]}{h_0 [1+\lambda_0(1-\varepsilon)^3]^2}$$
This equation shows that the stiffness varies with the eccentricity ratio and hence with the minimum film thickness.
Since no direct contact occurs in a hydrostatic bearing, the dominant motion resistance arises from viscous shear within the oil film. The viscous resistance force \(F_n\) can be estimated by:
$$F_n = \mu A_F \frac{v}{h} + \mu A_Y \frac{v}{H}$$
where \(\mu\) is the dynamic viscosity of the oil, \(A_F\) is the area of the sealing lands, \(A_Y\) is the area of the hydrostatic pocket, \(H\) is the pocket depth, and \(v\) is the relative velocity. For a gear shaping machine spindle, the average linear stroking speed \(V\) is related to the strokes per minute \(n\) and the stroke length \(L\) (0.05 m in this case) by:
$$V = \frac{2Ln}{60}$$
This study focuses on a condition where the average unidirectional stroking speed is 2 m/s.
Numerical Analysis: Influence of Axial Groove Presence
Impact on Load Capacity and Stiffness
Finite element models of the oil film domain, both with and without an axial groove on the circumferential sealing land, were established. A constant pressure supply with capillary restrictors at the inlets was assumed. Factors such as assembly precision, surface roughness, and thermal effects were initially neglected. The simulation parameters are summarized in Table 1.
| Parameter | Symbol | Value | Unit |
|---|---|---|---|
| Shaft Diameter | d | 90 | mm |
| Radial Clearance | h0 | 0.02 | mm |
| Oil Pocket Axial Length | L | 100 | mm |
| Axial Land Width | l | 10 | mm |
| Oil Density | ρ | 880 | kg/m³ |
| Oil Specific Heat Capacity | C | 1884 | J/(kg·K) |
| Oil Pocket Circumferential Width | B | 51 | mm |
| Circumferential Land Width | b | 19 | mm |
| Oil Pocket Angle | θ | 60 | ° |
| Pocket Depth | H | 2 | mm |
| Oil Dynamic Viscosity | μ | 0.045 | Pa·s |
| Thermal Conductivity | K | 0.132 | W/(m·K) |
Simulations were conducted for various eccentricity ratios (ε = 0 to 0.5) to represent different radial load levels. The pressure distribution for ε = 0.3 revealed that pressure dropped sharply from the circumferential sealing land towards the outlet, while remaining nearly uniform within the pocket. Without an axial groove, the pressure zone extended across the loaded pocket and part of the sealing land. With the groove, pressure became more concentrated within the pocket itself, effectively reducing pressure on the adjacent sealing land.
The calculated load capacity as a function of eccentricity is summarized in Table 2. The load capacity increases with eccentricity, but the rate of increase diminishes. The presence of the axial groove enhances load capacity by approximately a factor of 1.5 at low eccentricity and about 1.3 at higher eccentricity.
| Eccentricity Ratio (ε) | Load Capacity (No Groove) [N] | Load Capacity (With Groove) [N] | Enhancement Factor |
|---|---|---|---|
| 0.1 | ~4500 | ~6750 | 1.50 |
| 0.3 | ~7000 | ~9800 | 1.40 |
| 0.5 | ~8500 | ~11000 | 1.29 |
The oil film stiffness, calculated using the derivative of the load-eccentricity data, decreases with increasing eccentricity. At the same eccentricity, the configuration without a groove exhibits marginally higher stiffness. This difference is more pronounced at lower eccentricities and diminishes as eccentricity increases.
Impact on Viscous Resistance and Flow
The viscous shear force distribution for V=2 m/s shows that shear stress is highest in the thin film regions of the sealing lands and significantly lower within the deep pocket. The shear force distribution is asymmetric relative to the stroking direction. Monitoring planes on the sealing lands revealed a highly non-uniform flow distribution under eccentric conditions, with flow at the maximum film thickness location being up to 5 times greater than at the minimum film thickness location on the leading land. The axial groove structure increases the total flow across the sealing lands and reduces the effective area subject to high shear, thereby lowering the overall viscous resistance force for the spindle, as shown in Table 3 for different eccentricities.
| Eccentricity Ratio (ε) | Viscous Resistance (No Groove) [N] | Viscous Resistance (With Groove) [N] | Reduction |
|---|---|---|---|
| 0.1 | ~125 | ~118 | ~5.6% |
| 0.3 | ~132 | ~122 | ~7.6% |
| 0.5 | ~140 | ~130 | ~7.1% |
Thermal Behavior During Continuous Stroking
In high-speed gear shaping, viscous shear heating within the oil film can lead to temperature rise. Transient simulations over multiple strokes were performed. The temperature distribution is primarily concentrated on the sealing lands where the film is thinnest and shear is highest. The temperature distribution is uneven, with the sealing land opposite the direction of motion experiencing higher temperatures due to lower flushing flow. The axial groove reduces the area of concentrated shear heating and promotes better flow, aiding in heat evacuation. Over ten consecutive strokes, the average oil film temperature with the groove stabilized around 28°C with minor fluctuations, whereas without the groove, it showed a gradual increasing trend.
Impact on Sleeve Deformation (Fluid-Structure Interaction)
A one-way fluid-structure interaction analysis was conducted to evaluate the deformation of the bronze alloy sleeve under fluid pressure. The location and magnitude of maximum deformation differ between the two designs, as summarized in Table 4. Without a groove, the maximum deformation shifts from the circumferential land at low eccentricity towards the pocket center at high eccentricity. With a groove, the maximum deformation consistently occurs at the pocket center. At low eccentricity (ε < 0.3), the sleeve without a groove deforms more (up to 40% greater). At high eccentricity (ε > 0.3), the sleeve with a groove deforms slightly more (up to 15% greater).
| Condition | Low Eccentricity (ε ~0.1) | High Eccentricity (ε ~0.5) | Deformation Trend |
|---|---|---|---|
| No Axial Groove | Higher deformation (~40% more than with groove). Location: Circumferential land. | Lower deformation. Location: Pocket center & land. | Increases by factor of ~1.2 from low to high ε. |
| With Axial Groove | Lower deformation. Location: Pocket center. | Higher deformation (~15% more than without groove). Location: Pocket center. | Increases by factor of ~1.9 from low to high ε. |
Influence of Axial Groove Cross-Sectional Shape
Further simulations investigated the effect of different groove shapes (rectangular, circular, triangular, trapezoidal) with equal projected area on the sealing land. The supply pressure was set at 2 MPa and the instantaneous speed at 2 m/s. The results for load capacity and stiffness are summarized in Tables 5 and 6.
| Eccentricity Ratio (ε) | Rectangular | Triangular | Circular | Trapezoidal |
|---|---|---|---|---|
| 0.1 | High | Highest | Medium | Medium |
| 0.3 | Highest | High | Low | Medium |
| 0.5 | Highest | Medium | Low | High |
While all shapes show increasing load capacity with eccentricity, the influence of shape is relatively minor. The triangular groove offers the best performance at very low eccentricity, while the rectangular groove generally provides the highest capacity at moderate to high eccentricity.
| Eccentricity Ratio (ε) | Stiffness Ranking (Highest to Lowest) | Key Observation |
|---|---|---|
| 0.1 | Triangular > Rectangular > Trapezoidal > Circular | Shape has significant effect (~10% diff. between highest & lowest). |
| 0.3 | Trapezoidal > Rectangular > Triangular > Circular | Shape has notable effect. |
| 0.5 | Rectangular > Trapezoidal > Triangular > Circular | Shape effect is present but relatively smaller. |
The groove shape has a more pronounced effect on stiffness, particularly at lower eccentricities (ε = 0.1 and 0.3). The triangular groove provides the highest initial stiffness, but its stiffness degrades fastest with increasing eccentricity. The rectangular and trapezoidal grooves maintain more consistent performance across the eccentricity range.
Experimental Validation
To validate the modeling approach, the established finite element method was applied to a benchmark case from literature: a cylindrical hydrostatic guideway. The experimental setup and parameters from a known study were used to create a corresponding simulation model for a single guideway under a 150 kg load. The key comparison is shown in Table 7.
| Stiffness Value Source | Oil Film Stiffness [N/μm] | Deviation from Experiment |
|---|---|---|
| Experimental Measurement | 160 | 0% (Baseline) |
| Theoretical Calculation (Literature) | 166.4 | +3.8% |
| Present Simulation Model | 163.4 | +2.0% |
The close agreement between the simulation result and the experimental measurement (within 2%) confirms the validity and accuracy of the finite element model and methodology used in this study for analyzing hydrostatic systems relevant to gear shaping machine components.
Conclusions
This investigation into the effect of an axial groove structure on the performance of a hydrostatic spindle for high-speed gear shaping machines yields the following conclusions:
- The incorporation of an axial groove on the circumferential sealing land significantly improves pressure distribution, concentrating it within the hydrostatic pocket. This enhances the oil film’s load-carrying capacity by a factor of 1.3 to 1.5 and provides superior stiffness at lower eccentricities compared to the no-groove design.
- The groove structure reduces the effective area of high viscous shear, increases lubricant flow across the bearing lands, and consequently lowers the total viscous resistance during the high-speed stroking motion characteristic of gear shaping operations.
- Regarding thermal management, the axial groove mitigates viscous heating by reducing the zone of concentrated shear on the sealing land and enhancing convective heat removal via increased flow. This leads to a stabilized average oil film temperature over multiple strokes, which is critical for maintaining consistent performance and oil viscosity in a demanding gear shaping environment.
- Fluid-structure interaction analysis reveals that the groove alters the deformation pattern of the support sleeve. While the no-groove design exhibits larger deformation at low loads, the grooved design shows slightly higher deformation at very high loads, with the maximum deformation consistently located at the pocket center.
- The cross-sectional shape of the axial groove has a moderate influence on performance. While its effect on load capacity is relatively small, the choice of shape (rectangular, triangular, trapezoidal, circular) has a more notable impact on oil film stiffness, especially at lower eccentricities relevant to precision machining in gear shaping. A rectangular groove offers a good balance of high load capacity and stable stiffness across the operating range.
This study provides valuable insights and a validated numerical framework for optimizing the structural parameters of hydrostatic spindles, directly contributing to the advancement of high-performance, reliable gear shaping machines. Future work will focus on incorporating transient thermal effects, surface roughness, and the combined impact of stroking and slow rotation to further refine the performance predictions for these critical machine tool components.
