Numerical Study of the Optimum Fiber Content of Sealing Grease Using Discrete Element Method
Abstract
:1. Introduction
2. Discrete Element Modeling
2.1. Particle Flow Code
2.2. Discrete Element—Geometric Model
- The pressure column was modeled with a 40 mm × 40 mm × 50 mm rectangular. The sieve was modeled using a series of regularly arranged triangular unit walls with the screen opening size of 2.4 mm × 2.85 mm and 5 mm × 5 mm for sealing tests and rheology tests, respectively, as shown in Figure 3a.
- A total of 5533 powder spheres were generated using the “ball distribute” command, with a porosity of 0.6 and a radius of 1 mm (Figure 3b) [25,26,27]. The selected size of the powder spheres was small enough to extrude from the sieve openings, although it cannot be the same size as the actual particles at micron order.
- The model of single fiber was created by a series of spheres with a radius of 1.5 mm and an embedded length of 0.5 mm (Figure 3c) [9]. It looked similar to a pearl necklace where the pearls geometrically interpenetrate with each other. A parallel contact model was adopted between the spheres, which can guarantee that the fiber does not suffer a compressive fracture during compression. The number of spheres for every single fiber depended on the fiber’s aspect ratio, while the fiber concentration determined the number of fibers. A series of fibers with a certain number and random direction were generated in the pressure column using a nested loop structure, as shown in Figure 3d. The physical parameters of the models for the sealing grease simulations are provided in Table 1.
- The tackifier acts as a bond between fiber spheres and powder spheres, which was simulated by a contact model rather than spheres. The parameters of its contact model were determined by the type of tackifier, which is significant for the properties of sealing grease samples.
2.3. Computational Procedure of Pressure Resistant Sealing
2.4. The Criterion for Sealing State
2.5. Contact Element—Mechanical Modeling
2.5.1. Experiment of Rheology Test
2.5.2. Simulation of Rheology Test
2.6. Parametric Study-Viscosity and Fiber Aspect Ratio
3. Results and Discussion
3.1. The Effect of Viscosity Grade on Optimum Fiber Number
3.2. The Effect of Fiber’s Aspect Ratio on Optimum Fiber Number
4. Conclusions
- (1)
- Viscosity grade had a significant impact on the pumpability of sealing grease, which was verified by experiment and DEM simulation of rheology tests. However, the viscosity of sealing grease showed little influence on the optimum fiber number. Noting that the increase in viscosity can improve the sealing effect, apparently.
- (2)
- Compared with the influence of viscosity, the increase in fiber’s aspect ratio reduced the optimum fiber number, which was verified by the curves of unbalanced force and fiber area proportion of the DEM models. Furthermore, the variation curves of the unbalanced force and fiber area proportion were practically identical and can be divided into four stages. Although the optimum fiber number was dependent on fiber’s aspect ratio, the mass of fiber is related to the total number of fiber spheres. This study shows that the sealing grease sample made of the same fiber material and quality can reach the same seal state and seal effect.
- (3)
- The DEM models and simulation results provided a solid and convenient tool to predict the quality of sealing grease. Many other parameters affecting the optimum fiber content and sealing performance, such as fiber type and fiber gradation, should be further considered in the DEM simulation for a more complete parametric study for guiding and improving the sealing grease design with numerical simulation.
Author Contributions
Funding
Institutional Review Board Statement
Informed Consent Statement
Data Availability Statement
Acknowledgments
Conflicts of Interest
References
- Sharghi, M.; Chakeri, H.; Ozcelik, Y. Investigation into the effects of two component grout properties on surface settlements. Tunn. Undergr. Space Technol. 2017, 63, 205–216. [Google Scholar] [CrossRef]
- Lin, C.; Wu, S.; Zhang, Z. One Technique for Shield Tail Brushes’ Repair Using Liquid Nitrogen Freezing Method in the Artesian Aquifer. ICPTT 2011, 2011, 858–867. [Google Scholar] [CrossRef]
- Zhang, D.; Huang, Z.; Yin, Z.-Y.; Ran, L.; Huang, H.-W. Predicting the grouting effect on leakage-induced tunnels and ground response in saturated soils. Tunn. Undergr. Space Technol. 2017, 65, 76–90. [Google Scholar] [CrossRef]
- Wang, D.Q. Formulation Research and Property Charaterization of A Kind of Shield Tail Sealing Grease. Tunn. Constr. 2013, 33, 277–280. [Google Scholar] [CrossRef]
- Matsuda, Y.; Iwase, Y. Numerical simulation of rock fracture using three-dimensional extended discrete element method. Earth Planets Space 2002, 54, 367–378. [Google Scholar] [CrossRef] [Green Version]
- Potyondy, D.O.; Cundall, P.A. A bonded-particle model for rock. Int. J. Rock Mech. Min. Sci. 2004, 41, 1329–1364. [Google Scholar] [CrossRef]
- Tiwari, R.; Kumar, R.; Abouelregal, A.E. Analysis of a magneto-thermoelastic problem in a piezoelastic medium using the non-local memory-dependent heat conduction theory involving three phase lags. Mech. Time Depend. Mater. 2021, 1–17. [Google Scholar] [CrossRef]
- Tiwari, R.; Mukhopadhyay, S. Analysis of wave propagation in the presence of a continuous line heat source under heat transfer with memory dependent derivatives. Math. Mech. Solids 2017, 23, 820–834. [Google Scholar] [CrossRef]
- Zhang, D.; Tu, H.; Li, Y.; Weng, Y. Effect of fiber content and fiber length on the dynamic compressive properties of strain-hardening ultra-high performance concrete. Constr. Build. Mater. 2022, 328, 127024. [Google Scholar] [CrossRef]
- Kimura, H.; Imai, Y.; Yamamoto, Y. Study on Fiber Length Control for Ester-Based Lithium Soap Grease. Tribol. Trans. 2001, 44, 405–410. [Google Scholar] [CrossRef]
- Cundall, P.A.; Strack, O.D.L. A discrete numerical model for granular assemblies. Géotechnique 1979, 29, 47–65. [Google Scholar] [CrossRef]
- Leclerc, W.; Haddad, H.; Guessasma, M. On the suitability of a Discrete Element Method to simulate cracks initiation and propagation in heterogeneous media. Int. J. Solids Struct. 2017, 108, 98–114. [Google Scholar] [CrossRef]
- Le, B.; Dau, F.; Charles, J.; Iordanoff, I. Modeling damages and cracks growth in composite with a 3D discrete element method. Compos. Part B Eng. 2016, 91, 615–630. [Google Scholar] [CrossRef] [Green Version]
- Khattak, M.J.; Khattab, A. Modeling tensile response of fiber-reinforced polymer composites using discrete element method. Polym. Compos. 2013, 34, 877–886. [Google Scholar] [CrossRef]
- Yang, M.; Li, S.; Yao, Q. Mechanistic studies of initial deposition of fine adhesive particles on a fiber using discrete-element methods. Powder Technol. 2013, 248, 44–53. [Google Scholar] [CrossRef]
- Cai, R.-R.; Zhang, L.-Z. Modeling of dynamic deposition and filtration processes of airborne particles by a single fiber with a coupled lattice Boltzmann and discrete element method. Build. Environ. 2016, 106, 274–285. [Google Scholar] [CrossRef]
- Ganapathy, T.; Sathiskumar, R.; Sanjay, M.R.; Senthamaraikannan, P.; Saravanakumar, S.S.; Parameswaranpillai, J.; Siengchin, S. Effect of Graphene Powder on Banyan Aerial Root Fibers Reinforced Epoxy Composites. J. Nat. Fibers 2019, 18, 1029–1036. [Google Scholar] [CrossRef]
- Su, H.; Fu, Z.; Wen, Z.; Gao, A. Particle Flow Code Method-Based Meso-scale Identification for Seepage Failure of Soil Levee. Transp. Porous Media 2017, 119, 311–336. [Google Scholar] [CrossRef]
- He, J.; Li, X.; Li, S.; Yin, Y.; Qian, H. Study of seismic response of colluvium accumulation slope by particle flow code. Granul. Matter 2010, 12, 483–490. [Google Scholar] [CrossRef]
- Mehranpour, M.H.; Kulatilake, P.H. Improvements for the smooth joint contact model of the particle flow code and its applications. Comput. Geotech. 2017, 87, 163–177. [Google Scholar] [CrossRef]
- Ding, X.; Zhang, L. A new contact model to improve the simulated ratio of unconfined compressive strength to tensile strength in bonded particle models. Int. J. Rock Mech. Min. Sci. 2014, 69, 111–119. [Google Scholar] [CrossRef]
- Potyondy, D.O. Simulating stress corrosion with a bonded-particle model for rock. Int. J. Rock Mech. Min. Sci. 2007, 44, 677–691. [Google Scholar] [CrossRef]
- Gang, Y.B.; Wabf, J.F.; Yang, G.B. PFC Numerical Simulation on Anisotropic Properties of Cohesive Soil. J. Eng. Geol. 2009, 17, 638–642. [Google Scholar] [CrossRef]
- Khattab, A.; Khattak, M.J.; Fadhil, I.M. Micromechanical discrete element modeling of fiber reinforced polymer composites. Polym. Compos. 2011, 32, 1532–1540. [Google Scholar] [CrossRef]
- Wei, P.; Wei, Z.; Chen, Z.; He, Y.; Du, J. Thermal behavior in single track during selective laser melting of AlSi10Mg powder. Appl. Phys. A 2017, 123, 604. [Google Scholar] [CrossRef]
- Basinskas, G.; Sakai, M. Numerical study of the mixing efficiency of a ribbon mixer using the discrete element method. Powder Technol. 2016, 287, 380–394. [Google Scholar] [CrossRef]
- Tsunazawa, Y.; Shigeto, Y.; Tokoro, C.; Sakai, M. Numerical simulation of industrial die filling using the discrete element method. Chem. Eng. Sci. 2015, 138, 791–809. [Google Scholar] [CrossRef]
- Li, X.; Yang, Y.; Li, F. Comparative Study on Mechanical Properties of Sealing Grease Composed of Different Base Oils for Shield Tunnel. Materials 2020, 13, 692. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Wang, W.; Gu, W.; Liu, K.; Wang, F.; Tang, Z. DEM Simulation on the Startup Dynamic Process of a Plain Journal Bearing Lubricated by Granular Media. Tribol. Trans. 2014, 57, 198–205. [Google Scholar] [CrossRef]
Parameter | Dimension | Unit |
---|---|---|
Pressure column | 4 × 4 × 5 | cm |
Screen opening of seal test | 2.4 × 2.85 | mm |
Screen opening of rheology test | 5 × 5 | mm |
Radius of powder | 1 | mm |
Radius of fiber | 1.5 | mm |
Length of fiber (34) | 5 | cm |
Pressure (rheology test) | 1000 | Pa |
Pressure (sample loading) | 1000 | Pa |
Pressure (pressure resistant sealing test) | 3.5 | MPa |
Experiment | PB1300:PB2400 | 5:0 | 4:1 | 3:2 | 2:3 | 1:4 | 0:5 |
---|---|---|---|---|---|---|---|
Simulation | Normal stiffness (N/m) | 1 × 105 | 1 × 105 | 1 × 105 | 1 × 105 | 1 × 105 | 1 × 105 |
Friction force (N) | 1 × 102 | 1 × 102 | 1 × 102 | 1 × 102 | 1 × 102 | 1 × 102 | |
Normal critical damping ratio | 1 × 106 | 1.5 × 106 | 2 × 106 | 2.5 × 106 | 3 × 106 | 3.5 × 106 | |
Tangential critical damping ratio | 1 × 106 | 1.5 × 106 | 2 × 106 | 2.5 × 106 | 3 × 106 | 3.5 × 106 | |
Tensile strength (Pa) | 1 × 103 | 1 × 103 | 1 × 103 | 1 × 103 | 1 × 103 | 1 × 103 | |
Cohesion (Pa) | 1 × 102 | 1 × 102 | 1 × 102 | 1 × 102 | 1 × 102 | 1 × 102 | |
Viscosity Grade | 1 | 2 | 3 | 4 | 5 | 6 |
Fiber’s Aspect Ratio | Fiber Number | Number of Drops (Stability) | Pore Distribution at Cross Section |
---|---|---|---|
17 | 500 | 2078 | One Large pore at up left corner and one tiny pore at bottom right corner |
700 | 1673 | One tiny pore at up left corner | |
750 | 1300 | One tiny pore at up left corner | |
900 | 993 | No pores | |
20 | 500 | 1877 | One Large pore at up left corner |
600 | 1256 | Two tine pores at up left corner and up right corner | |
700 | 1190 | One tiny pore at up left corner | |
850 | 967 | No pores | |
24 | 400 | 1890 | One large pore at bottom left corner and two tiny pores at up right corner |
500 | 1760 | One large pore at up left corner and one tiny pore at bottom | |
600 | 1077 | One tiny pore at bottom | |
700 | 923 | No pores | |
28 | 400 | 1430 | Two tiny pores at bottom |
550 | 1296 | One tiny pore at up left corner | |
580 | 1162 | One tiny pore at up left coener | |
650 | 960 | No pores | |
34 | 350 | 1630 | Two tiny pores at bottom right corner |
400 | 1295 | One tiny pore at bottom | |
500 | 921 | One tiny pore at bottom | |
600 | 867 | No pores |
Fiber’s Aspect Ratio | Cross-Sectional Area | Fiber Area | Proportion of Fiber Area |
---|---|---|---|
17 | 52,269 | 81,859 | 80.08% |
20 | 80,052 | 62,264 | 77.91% |
24 | 62,300 | 49,035 | 78.71% |
28 | 85,069 | 66,796 | 78.52% |
34 | 70,754 | 56,851 | 80.35% |
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
Zhou, X.; Wei, Y.; Yang, Y.; Xu, P. Numerical Study of the Optimum Fiber Content of Sealing Grease Using Discrete Element Method. Materials 2022, 15, 3485. https://doi.org/10.3390/ma15103485
Zhou X, Wei Y, Yang Y, Xu P. Numerical Study of the Optimum Fiber Content of Sealing Grease Using Discrete Element Method. Materials. 2022; 15(10):3485. https://doi.org/10.3390/ma15103485
Chicago/Turabian StyleZhou, Xiong, Yingjie Wei, Yuyou Yang, and Pengfei Xu. 2022. "Numerical Study of the Optimum Fiber Content of Sealing Grease Using Discrete Element Method" Materials 15, no. 10: 3485. https://doi.org/10.3390/ma15103485
APA StyleZhou, X., Wei, Y., Yang, Y., & Xu, P. (2022). Numerical Study of the Optimum Fiber Content of Sealing Grease Using Discrete Element Method. Materials, 15(10), 3485. https://doi.org/10.3390/ma15103485