Mixed Lubrication Modeling of Multi-Lip Reciprocating Seals Based on Elastohydrodynamic Lubrication Theory
Abstract
:1. Introduction
2. Theoretical Models
2.1. Geometrical Model
2.2. Solid Mechanics
2.3. Fluid Mechanics
2.4. Contact Mechanics
2.5. Deformation Mechanics
2.6. Check Convergence and Auxiliary Calculation
2.6.1. Check Convergence
2.6.2. Auxiliary Calculation
2.7. Computational Scheme
- Input the fluid property parameters (including fluid viscosity μ, pressure–viscosity coefficient α, etc.), the operating conditions parameters (including fluid sealed pressure Ps, piston rod speed U) and the seal surface parameters (including RMS roughness σ, asperity density η, etc.) in the MATLAB software.
- Based on the basic structural parameters of the DAS sealing system, the large-scale finite element general analysis software ANSYS is used for solid mechanics calculation. The static contact pressure Psc, contact length Lx, and deformation coefficient matrix K in the sealing zone are obtained, respectively.
- The initial static film thickness Hs obtained by Psc was used to solve the Reynolds equation in fluid mechanics, until Pf converges, otherwise Φ is updated.
- The G–W contact model in contact mechanics was solved by Hs, and the asperity contact pressure Pc was obtained.
- Compare the sum of the fluid film pressure Pf and the asperity contact pressure Pc with the static contact pressure Psc to determine whether the three forces are balanced. If they are not balanced, the film thickness H is updated using the deformation coefficient matrix in the deformation mechanics and iterated repeatedly until the three forces are balanced.
- Auxiliary calculation of the flow through the three seal lips, whether , , and are equal. If not equal, return to step 3 to update the fluid pressures Pi1 and Pi2 between the inter-lip. Iterate the loop again until .
- The program calculation is completed and the results are output.
3. Simulation Parameters Setup
4. Results and Discussion
4.1. Macro Seal Characteristics
4.2. Micro Seal Characteristics
4.2.1. Lubrication Characteristics
4.2.2. Effects of Sealed Pressure and Piston Speed
4.2.3. Effects of Sealed Pressure and Surface Roughness
4.2.4. Effects of Piston Speed and Surface Roughness
5. Experimental Study
5.1. Experimental Apparatus
5.2. Measurement Principle
5.3. Experimental Results
6. Conclusions
- (1)
- In mixed lubrication state, the asperity contact pressure increases obviously with the increase of sealed pressure, while the film thickness decreases, and the asperity contact pressure of the main lip is the largest. It shows that during the high-pressure service of the hydraulic cylinder, the main lip will produce more wear and accelerate the failure of the seal. With the increase of the piston rod speed, the asperity contact pressure decreases relatively, and the increase of the film thickness is not obvious, while the asperity contact pressure and the film thickness both increase with the increase of the seal surface roughness.
- (2)
- With the increase of piston rod speed and seal surface roughness, the leakage increases. When the seal surface roughness is within a certain range, the leakage decreases significantly with the sealed pressure increase. Therefore, to minimize internal leakage, the hydraulic cylinder piston seal surface roughness should be selected reasonably under the condition that the hydraulic cylinder meets the requirements of working conditions.
- (3)
- When the pre-compression amount of the seal is constant, the friction force is mainly determined by the sealed pressure. With the increase of sealed pressure, the friction force of the sealing zone increases significantly. Although increasing the piston rod speed and seal surface roughness can increase the film thickness, the viscous friction force is relatively weak, so the piston rod speed and seal surface roughness have little influence on the total friction force.
Author Contributions
Funding
Data Availability Statement
Conflicts of Interest
Nomenclature
C10, C01, d | Moony–Rivlin coefficients |
Ds | Seal outside diameter |
ds | Seal inside diameter |
E | Elastic modulus |
F | Cavitation index |
Ff | Friction force |
FL | Load force of the hydraulic cylinder |
f | Friction coefficient for asperity contact |
h | Fluid film thickness |
H | Dimensionless average film thickness, h/σ |
Hs | Dimensionless static film thickness, hs/σ |
HT | Dimensionless average truncated film thickness, hT/σ |
I1, I2 | Green–Lagrange strain invariant |
J | Elastic volume ratio |
K | Influence coefficient matrix |
L | Contact length of sealing zone |
Lstroke | Piston/rod stroke length |
pa | Ambient pressure |
Pc | Dimensionless asperity contact pressure, pc/E |
Pf | Dimensionless fluid film pressure, pf/pa |
Dimensionless fluid pressure in current iterated step | |
Dimensionless fluid pressure in last iterated step | |
Pi1, Pi2 | Pressure in the first and second inter-lip region |
Pl | Dimensionless pressure in low-pressure side, pl/pa |
Ps | Dimensionless sealed pressure, ps/pa |
Psc | Dimensionless static contact pressure, psc/E |
Q | Flow rate per stroke |
Dimensionless flow rate, | |
R | Radius of asperities |
U | Dimensionless piston rod speed, |
WM−R | Strain-energy function |
Ws | Width of the seal |
Dimensionless coordinate parallel to the fluid film thickness, x/L | |
z | Dimensionless coordinate normal to the fluid film thickness |
Dimensionless pressure–viscosity coefficient, αpa | |
ξ | |
Φ | Fluid pressure/density function |
ϕf, ϕfs, ϕfp | Shear stress factors |
ϕs.c.x, ϕxx | flow factor |
η | Asperity density |
μ0 | Viscosity at atmospheric pressure |
v | Poisson’s ratio |
ρf | Fluid density |
Dimensionless density, ρ/ρf | |
Dimensionless RMS roughness of the seal, σR1/3η2/3 | |
Dimensionless viscous shear stress, | |
Dimensionless asperity shear stress, |
Appendix A
- Fluid Mechanics
- Contact mechanics
References
- Nikas, G. Eighty years of research on hydraulic reciprocating seals: Review of tribological studies and related topics since the 1930s. Proc. Inst. Mech. Eng. Part J J. Eng. Tribol. 2010, 224, 1–23. [Google Scholar] [CrossRef]
- White, C.; Denny, D.; Britain, G. The Sealing Mechanism of Flexible Packings; HM Stationery Office: London, UK, 1948; pp. 35–62. [Google Scholar]
- Karaszkiewicz, A. Hydrodynamic lubrication of rubber seals for reciprocating motion; leakage of seals with an O-ring. Tribol. Int. 1988, 21, 361–367. [Google Scholar] [CrossRef]
- Kaneta, M.; Todoroki, H.; Nishikawa, H.; Kanzaki, Y.; Kawahara, Y. Tribology of flexible seals for reciprocating motion. J. Tribol. 2000, 122, 787–795. [Google Scholar] [CrossRef]
- Öngün, Y.; André, M.; Bartel, D.; Deters, L. An axisymmetric hydrodynamic interface element for finite-element computations of mixed lubrication in rubber seals. Proc. Inst. Mech. Eng. Part J J. Eng. Tribol. 2008, 222, 471–481. [Google Scholar] [CrossRef]
- Peng, C.; Ouyang, X.; Zhu, Y.; Guo, S.; Zhou, Q.H.; Yang, H.Y. Investigation into the influence of stretching on reciprocating rod seals based on a novel 3-D model vs axisymmetric model. Tribol. Int. 2018, 117, 1–14. [Google Scholar] [CrossRef]
- Yuan, X.; Wang, J.; Lian, Z.; Wang, G. Partial lubrication modeling of reciprocating rod seals based on a developed EHL method. Tribol. Int. 2021, 153, 106585. [Google Scholar] [CrossRef]
- Salant, R.; Yang, B.; Thatte, A. Simulation of hydraulic seals. Proc. Inst. Mech. Eng. Part J J. Eng. Tribol. 2010, 224, 865–876. [Google Scholar] [CrossRef]
- Crudu, M.; Fatu, A.; Cananau, S.; Hajjam, M.; Pascu, A.; Cristescu, C. A numerical and experimental friction analysis of reciprocating hydraulic ‘U’ rod seals. Proc. Inst. Mech. Eng. Part J J. Eng. Tribol. 2012, 226, 785–794. [Google Scholar] [CrossRef]
- Huang, Y.; Salant, R.F. Numerical analysis of a hydraulic rod seal: Flooded vs. starved conditions. Tribol. Int. 2015, 92, 577–584. [Google Scholar] [CrossRef]
- Gadari, M.E.; Hajjam, M. Effect of the Grooved Rod on the Friction Force of U-Cup Hydraulic Rod Seal with Rough Lip. Tribol. Trans. 2018, 61, 661–670. [Google Scholar] [CrossRef]
- Wang, J.; Li, Y.; Lian, Z. Numerical Investigations on the Sealing Performance of a Reciprocating Seal Based on the Inverse Lubrication Method. J. Tribol. 2019, 141, 112201. [Google Scholar] [CrossRef]
- Nikas, G. Elastohydrodynamics and Mechanics of Rectangular Elastomeric Seals for Reciprocating Piston Rods. J. Tribol. 2003, 125, 60–69. [Google Scholar] [CrossRef]
- Nikas, G.; Almond, R.; Burridge, G. Experimental Study of Leakage and Friction of Rectangular, Elastomeric Hydraulic Seals for Reciprocating Motion from −54 to +135 °C and Pressures from 3.4 to 34.5 MPa. Tribol. Trans. 2014, 57, 846–865. [Google Scholar] [CrossRef]
- Bhaumik, S.; Kumaraswamy, A.; Guruprasad, S.; Bhandari, P. Investigation of friction in rectangular Nitrile-Butadiene Rubber (NBR) hydraulic rod seals for defence applications. J. Mech. Sci. Technol. 2015, 29, 4793–4799. [Google Scholar] [CrossRef]
- Mao, J.; Wang, W.; Liu, Y. Experimental and theoretical investigation on the sealing performance of the combined seals for reciprocating rod. J. Mech. Sci. Technol. 2012, 26, 1765–1772. [Google Scholar] [CrossRef]
- Peng, C.; Guo, S.; Ouyang, X.; Zhou, Q.; Yang, H. Mixed Lubrication Modeling of Reciprocating Seals Based on a Developed Multiple-Grid Method. Tribol. Trans. 2018, 61, 1151–1161. [Google Scholar] [CrossRef]
- Xiang, C.; Guo, F.; Jia, X.; Wang, Y.; Huang, X. Thermo-elastohydrodynamic mixed-lubrication model for reciprocating rod seals. Tribol. Int. 2019, 140, 105894. [Google Scholar] [CrossRef]
- Wang, B.; Peng, X.; Meng, X. Simulation of the effects of non-Newtonian fluid on the behavior of a step hydraulic rod seal based on a power law fluid model. J. Zhejiang Univ. Sci. A 2018, 19, 824–842. [Google Scholar] [CrossRef]
- Ran, H.; Wang, S.; Liu, D. A multiscale wear model for reciprocating rod stepseal under mixed lubricating conditions based on linear elasticity. Proc. Inst. Mech. Eng. Part J J. Eng. Tribol. 2020, 235, 161–180. [Google Scholar] [CrossRef]
- Yang, B.; Salant, R. A Numerical Model of a Reciprocating Rod Seal With a Secondary Lip. Tribol. Trans. 2008, 51, 119–127. [Google Scholar] [CrossRef]
- Prati, E.; Strozzi, A. A Study of the Elastohydrodynamic Problem in Rectangular Elastomeric Seals. J. Tribol. 1984, 106, 505–512. [Google Scholar] [CrossRef]
- Dragoni, E.; Strozzi, A. Analysis of an Unpressurized, Laterally Restrained, Elastomeric O-Ring Seal. J. Tribol. 1988, 110, 193–200. [Google Scholar] [CrossRef]
- Ogden, R. Large Deformation Isotropic Elasticity: On the Correlation of Theory and Experiment for Compressible Rubberlike Solids. Proc. R. Soc. Lond. A 1972, 328, 567–583. [Google Scholar] [CrossRef]
- Rivlin, R. Large elastic deformations of isotropic materials IV. further developments of the general theory. Philos. Trans. R. Soc. Lond. A 1948, 241, 379–397. [Google Scholar]
- Zhang, J.; Han, C.; Liang, Z. Physics of failure analysis of power section assembly for positive displacement motor. J. Loss Prev. Process Ind. 2016, 44, 414–423. [Google Scholar] [CrossRef]
- Payvar, P.; Salant, R. A Computational Method for Cavitation in a Wavy Mechanical Seal. J. Tribol. 1992, 114, 199–204. [Google Scholar] [CrossRef]
- Salant, R.; Maser, N.; Yang, B. Numerical Model of a Reciprocating Hydraulic Rod Seal. J. Tribol. 2007, 129, 91–97. [Google Scholar] [CrossRef]
- Patir, N.; Cheng, H. An Average Flow Model for Determining Effects of Three-Dimensional Roughness on Partial Hydrodynamic Lubrication. J. Lubr. Technol. 1978, 100, 12–17. [Google Scholar] [CrossRef]
- Patir, N.; Cheng, H. Application of Average Flow Model to Lubrication Between Rough Sliding Surfaces. J. Lubr. Technol. 1979, 101, 220–229. [Google Scholar] [CrossRef]
- Greenwood, J.; Williamson, J. Contact of nominally flat surfaces. Proc. R. Soc. Lond. A 1966, 295, 300–319. [Google Scholar]
- EI Gadari, M.; Fatu, A.; Hajjam, M. Shaft roughness effect on elasto-hydrodynamic lubrication of rotary lip seals: Experimentation and numerical simulation. Tribol. Int. 2015, 88, 218–227. [Google Scholar] [CrossRef]
- Streator, J. A model of mixed lubrication with capillary effects. Tribol. Ser. 2002, 40, 121–128. [Google Scholar]
- Schmidt, T.; André, M.; Poll, G. A transient 2D-finite-element approach for the simulation of mixed lubrication effects of reciprocating hydraulic rod seals. Tribol. Int. 2010, 43, 1775–1785. [Google Scholar] [CrossRef]
Type | References | Method | Conclusion |
---|---|---|---|
O-ring | Öngün et al. [5] | EHL | The hydrodynamic interface element formula of seal mixed lubrication is realized in ABAQUS framework, but the calculation is very time consuming. |
Peng et al. [6] | EHL | Based on Salant’s model, a 3-D finite element model is established to analyze the effect of stretching on sealing performance. | |
Yuan et al. [7] | EHL | The elastic deformation equation is used to solve the elastic deformation of seals caused by fluid pressure in mixed lubrication model. | |
U-cup | Salant et al. [8] | EHL | The effects of cavitation and roughness on mixed lubrication are considered in 1-D Reynolds equation for the first time, and the theory is widely used. |
Crudu et al. [9] | IHL | Since the effect of surface roughness is ignored in the numerical model, the mixed lubrication effect cannot be reflected. | |
Huang et al. [10] | EHL | A mixed model of the boundary conditions of flooded fluid and starved fluid is established, and the lubricating region of the liquid film under the condition of flooded fluid is determined. | |
Gadari et al. [11] | EHL | Considering the effect of rod and seal surface roughness on the “transient” elastohydrodynamics in a 1-D numerical model, it is found that rod roughness increases the thickness of the fluid film. | |
Wang et al. [12] | IHL | Combined with the Greenwood–Williamson (G–W) contact model, a mixed lubrication model considering seal roughness is established, which overcomes the friction caused by the traditional fluid dynamics only calculating the fluid viscosity. | |
Rectangular ring | Nikas et al. [13,14] | EHL | It is found that the roughness has little effect on the film thickness. However, the model is limited to the sealing performance of rectangular seal and cannot adapt to the seal with complex geometric structure. |
Bhaumik et al. [15] | IHL | The friction force of the combined G–W contact model and the full film model, it is found that the friction force calculated by the combined G–W contact model is closer to the real result. | |
C-shape & T-shape | Mao et al. [16] | IHL | The full film method is used to compare the leakage of the two combined seals, and it is found that the rod speed had a significant effect on the seal leakage of the combined seal, but the model does not consider the effect of roughness. |
VL | Peng et al. [17] | EHL | The average pressure and interpolation procedure are used in the numerical program to solve the problems of mesh refinement and poor convergence. |
Xiang et al. [18] | EHL | A thermo-elastohydrodynamic mixed-lubrication model is established, and the simulation results show that the temperature of the sealing area changes with time, but the model ignores the transient effect. | |
Step | Wang et al. [19] | EHL | Based on the mixed lubrication model, the effects of non-Newtonian and Newtonian fluids on sealing performance are studied. |
Ran et al. [20] | EHL | The multiscale wear numerical program is established, but the program ignores that the wear will change the seal geometry and lead to the change of fluid film. |
Parameter Type | Value |
---|---|
Material of seal | NBR |
Elastic modulus of seal, Es (MPa) | 14 |
Poisson’s ratio of seal, vs | 0.499 |
Asperity radius, R (μm) | 1.0 |
Asperity density, η (m−2) | 1013 |
Inside diameter of seal, ds (mm) | 47 |
Outside diameter of seal, Ds (mm) | 63 |
Width of seal, Ws (mm) | 12 |
Material of backup ring | TPE |
Elastic modulus of backup ring, Eb (MPa) | 300 |
Poisson’s ratio of backup ring, vb | 0.470 |
Material of wearing ring | POM |
Elastic modulus of wearing ring, Ew (MPa) | 2600 |
Poisson’s ratio of wearing ring, νw | 0.380 |
Material of hydraulic and cylinder | C45E4 |
Diameter of piston groove, Dg (mm) | 47 |
Insider diameter of cylinder, Dc (mm) | 63 |
Diameter of piston rod, Dr (mm) | 36 |
Elastic modulus of piston groove and cylinder, Eg, Ec (MPa) | 2.1 × 105 |
Poisson’s ratio of piston groove and cylinder, νg, νc | 0.3 |
Stroke length of hydraulic cylinder, Lstroke (mm) | 300 |
Empirical friction coefficient, f | 0.2 |
Oil type | L-HM 46 |
Density of oil, ρl (Kg·m−3) | 842 |
Pressure–viscosity coefficient, α (Pa−1) | 2 × 10−8 |
Reference viscosity, μ0 (Pa·s) | 0.0387 |
Parameter Type | Value |
---|---|
Sealed pressure, ps (MPa) | 2–10 |
Piston rod extension speed, u (m/s) | 0.05–0.5 |
Seal RMS roughness, σ (μm) | 0.6–1.4 |
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Cheng, D.; Gu, L.; Sun, Y. Mixed Lubrication Modeling of Multi-Lip Reciprocating Seals Based on Elastohydrodynamic Lubrication Theory. Machines 2022, 10, 483. https://doi.org/10.3390/machines10060483
Cheng D, Gu L, Sun Y. Mixed Lubrication Modeling of Multi-Lip Reciprocating Seals Based on Elastohydrodynamic Lubrication Theory. Machines. 2022; 10(6):483. https://doi.org/10.3390/machines10060483
Chicago/Turabian StyleCheng, Donghong, Lichen Gu, and Yu Sun. 2022. "Mixed Lubrication Modeling of Multi-Lip Reciprocating Seals Based on Elastohydrodynamic Lubrication Theory" Machines 10, no. 6: 483. https://doi.org/10.3390/machines10060483
APA StyleCheng, D., Gu, L., & Sun, Y. (2022). Mixed Lubrication Modeling of Multi-Lip Reciprocating Seals Based on Elastohydrodynamic Lubrication Theory. Machines, 10(6), 483. https://doi.org/10.3390/machines10060483