Next Article in Journal
Tensile Strength of Class F Fly Ash and Fly Ash with Bentonite Addition as a Material for Earth Structures
Next Article in Special Issue
High-Temperature Tribological Behavior of HDPE Composites Reinforced by Short Carbon Fiber under Water-Lubricated Conditions
Previous Article in Journal
Special Issue: “Surface Engineering and Coating Technologies for Corrosion and Tribocorrosion Resistance”
Previous Article in Special Issue
Impact Abrasive Wear Property of CrAlN/TiSiN Multilayer Coating at Elevated Temperatures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Experimental and Numerical Analysis on the Impact Wear Behavior of TP316H Steel

1
Key Lab of Advanced Technologies of Materials, Tribology Research Institute, Southwest Jiao Tong University, Chengdu 610031, China
2
DEC Academy of Science and Technology Co., Ltd., Chengdu 611731, China
*
Author to whom correspondence should be addressed.
Materials 2022, 15(8), 2881; https://doi.org/10.3390/ma15082881
Submission received: 25 March 2022 / Revised: 7 April 2022 / Accepted: 12 April 2022 / Published: 14 April 2022

Abstract

:
In this work, the contact force model and experiment methods were used to study the dynamic response and impact wear behavior of TP316H steel. The Flore model and the classic Hertz model were selected for comparison with the experimental results, and the model was revised according to the section parameters of the TP316H tube. The results show that there is a large difference between the models without considering the effect of structural stiffness on the impact system and the test results, whereas the revised model has a good agreement. With the rise in impact mass, the coefficient of restitution increases from 0.65 to 0.78, whereas the energy dissipation and wear volume decrease. Spalling, delamination, plastic deformation, and oxidative wear are the main impact wear mechanism of TP316H steel.

1. Introduction

The contact-impact occurs in the mechanical system with clearance fits under the influence of external forces. Macroscopically, the impact process is a complex physical phenomenon, accompanied by the characteristics of very short duration, high-contact force, rapid energy dissipation, and large speed changes of the bodies [1]. At the microscopic level (contact interface), repeated impacts cause tribo-chemical reactions between the two bodies, which eventually lead to oxidation, wear, fatigue, and material transfer [2]. The contact-impact is a complex physical and chemical process under actual working conditions (such as high temperature or fluid). The contact-impact and the wear, fatigue, and cracks caused by it have been extensively studied [3,4,5,6,7].
The early contact force model was proposed by Hertz in 1882 [8], which is a completely elastic contact force model and does not consider the energy dissipation during the impact process. Based on Hertz’s contact theory, researchers improved the contact force model by adding either a linear spring connected to a linear damper or a nonlinear viscoelastic element [9,10]. The new models take into account the energy dissipation and have been well researched and developed [11,12,13]. At present, the consensus view is that the contact-impact process can be regarded as two stages (Figure 1b,c [11]: namely the compression, approaching or loading period, and the restitution, separating or unloading period [14,15]. In the compression stage, the contact force and deformation gradually increase from zero to the maximum value. The acceleration and velocity decrease to zero (Figure 1c). In return, the end of the compression stage is the beginning of the restitution stage, the contact force and deformation gradually decrease to zero when the two bodies are separated. Meanwhile, the velocity of the reverse will gradually increase (Figure 1c) [16]. In this process, energy dissipation due to internal damping (through vibration, heat, sound, etc.) usually occurs [17], as shown in Figure 1a,b. It can be defined by Newton’s Law of restitution, that is, the ratio of the velocity after and before the impact, which is also expressed as the coefficient of restitution [18]. The elastic strain energy stored by the material in the compression stage will finally become the kinetic energy of the rebound [19]. Although the theoretical contact force model has been extensively studied, there are few reports on its application to impact wear. This study combines the contact force model and impact wear test to study the dynamic response and wear behavior of the impact interface.
In the heat exchangers of nuclear power, the flow of the heat transfer medium will cause the vibration of the structure, such as the heat transfer tube, against the anti-vibration bar [20]. In the past few decades, scholars have carried out a lot of research on the impact wear of heat transfer tubes, including different contact forms, environments, and materials (Inconel 600, Inconel 690, and 316L) [21,22,23,24]. With the development of nuclear power technology, more and more attention is paid to candidate materials for the fourth generation reactor. The wear properties of 2.25Cr-1Mo steel, a candidate material for heat transfer tubes in a sodium-cooled fast reactor (SFR), are reported [25,26]. Another candidate material for the SFR is TP316H [27], which is used for heat exchangers [28]. There are few studies about the wear of TP316H [29,30]. Thus, this study takes the TP316H steel as the sample and combines existing models and experimental data to develop a revised contact force model applied to the thin-walled tube and discuss its correctness. The wear mechanism and energy dissipation during the impact wear process are discussed. This study also compared the differences between different models and test results, which is not only significant for the development of the contact force model but also provides data support for the engineering practice of TP316H steel.

2. Materials and Methods

2.1. Materials and Experimental Procedure

In this study, the impact wear of TP316H steel against 316H (the tube support plate material) was carried out under dry conditions. The main component of TP316H steel is Fe, and the mass percentages of other components are C(0.04–0.10), Mn(2.0), P(≤0.04), S(≤0.3), Si(≤0.075), Cr(16.0–18.0), Ni(11.0–14.0), and Mo(2.0–3.0). Table 1 show the mechanical properties of the tribo-pairs [29,30].
The self-developed impact wear test rig with controllable kinetic energy was used for this study, as shown in Figure 1a [31]. It can conveniently realize the impact wear test of different velocities, impact mass, and temperature. The impact velocity and contact force can be recorded in real-time. The energy dissipation and relative deformation between the contact bodies can be analyzed by the kinetic energy theorem and the integral of velocity against time, respectively. A detailed equipment introduction can be referred to in our previous work [31]. Table 2 shows the detailed experimental parameters. The initial impact kinetic energy was 3.5 mJ, the variable was the impact mass, and the impact frequency and the number of cycles were 3 Hz, 104. After the test, scanning electron microscopy (SEM, JSM 7800F, JEOL, Tokyo, Japan), energy-dispersive X-ray spectrometry (EDS, Aztec X-Max 80, Oxford Instrument, Abingdon, UK), and a 3D profiler were used to characterize the morphology, element distribution, and wear volume of the wear scars.

2.2. Contact Force Model for Impact Wear

P. A. Engel and R. G. Bayer used the Hertz contact force model to analyze the impact wear in 1974 [32]. A ball with a mass of m1 is assumed to impact a fixed plane, as shown in Figure 2a. The equation of motion for the ball is written in terms of contact force F N and the relative deformation δ [32]:
m δ ¨   +   F N   =   0 ,
where the δ ¨ is the acceleration of the ball, F N is based on an idealized analysis of the impact of two solid spheres, as shown in Figure 2b, and is expressed as [8]:
F N   =   K δ n ,
The n is the nonlinear power exponent (usually equal to 3/2). And K is the generalized parameter, which depends on the properties of the materials and the radii of the spheres [8]:
K   =   4 3 h 1   +   h 2     R 1 R 2 R 1   +   R 2 1 / 2 ,
where hl is defined for the l-th body as [8]:
h l   =   1     v l 2 E l ,   l   =   1 ,   2 .
and the variables vl, El is the Poisson’s ratio and Young’s modulus of the two spheres, respectively.
Even though the Hertz contact force model has defects, it was widely used until the 1960s. Then the Kelvin−Voigt model was proposed, taking into account the energy transfer (dissipation) during the impact, and is modeled as a linear spring and a linear damper element, as shown in Figure 2c. The contact force is expressed as [9]:
F N   =   K δ   +   D δ ˙ ,
The right side of the formula represents the elastic force and energy dissipation, where D is the hysteresis coefficient and δ ˙ is the relative velocity. Although Formula (5) considers energy dissipation, its constitutive law has some weaknesses. At the beginning and end of the impact process, the deformation and the contact force should be null, but the damping component in Formula (5) will make the contact force discontinuous at these two moments. Subsequently, Hunt and Crossley gave a new contact model with a nonlinear viscous-elastic element [10]:
F N   =   K δ n   +   χ δ n δ ˙ ,
where the χ is the hysteresis damping factor given by [10]:
χ   =   3 K 1     c r 2 δ ˙ ,
where cr is the coefficient of restitution, and δ ˙ is the initial relative velocity. As a result, the expression for this contact force model has the following form [10]:
F N   =   K δ n ( 1   +   3 1     c r 2 δ ˙ δ ˙ ) ,
Scholars expanded and developed more models on this basis. Among them, the method of comparing the kinetic energy lost during the impact with the energy obtained by the integration of the contact force on the deformation is widely accepted. The Lankarani−Nikravesh model was proposed in 1989, and its hysteresis damping factor is expressed as [33]:
χ   =   3 K 1     c r 2 4 δ ˙ ,
The meaning of the variables in Formula (9) is consistent with the previous formula. The contact force can be expressed as [33]:
F N   =   K δ n ( 1   +   3 1     c r 2 4 δ ˙ δ ˙ ) ,
The Flore model was proposed in 2010, and its hysteresis damping factor is expressed as [17]:
χ   =   8 K 1     c r 5 c r δ ˙ ,
the contact force can be expressed as [17]:
F N   =   K δ n ( 1   +   8 1     c r 5 c r δ ˙ δ ˙ ) ,
In recent years, the comparison between various contact force models has been widely reported [1,14,15,17]. For “medium” values of the coefficient of restitution, the Flore model is better at predicting the properties of impact energy dissipation [14]. The Flore model and classic Hertz model were selected to compare with the experimental results in this study.

3. Results

3.1. Impact Dynamic Response of Thin-Walled Tube

Table 3 lists the experimental results and calculation parameters under varying impact mass. The pre-impact and post-impact velocity are recorded through the experiment, and the generalized parameter (K), coefficient of restitution (cr), and hysteresis damping factor ( χ ) are calculated. These pre-parameters and Formulas (1) and (12) are written as subroutines of MATLAB and solved.
Figure 3 shows the curve of relative velocity versus contact time under varying impact mass. According to the kinetic energy theorem, to keep the initial impact kinetic energy consistent, the pre-impact velocity must decrease as the impact mass increases, as shown in Figure 3a. The results show that post-impact velocity also gradually decreases with the increase in the impact mass, but the magnitude of the decrease is less than the initial impact velocity (Table 3). The results of the Flore model agree with the experimental results because the coefficient of restitution is the same, as shown in Figure 3. The coefficient of restitution of the Hertz model is 1, and its pre-impact velocity is equal to the post-impact velocity, as shown in Figure 3b.
Figure 4 and Figure 5 show the curve of contact force and relative deformation versus contact time under varying impact mass. The experimental results show a slight difference between the first and last impact, as shown in Figure 3a, Figure 4a and Figure 5a. Because the repeated impact process changes the contact area, the accumulation of wear debris simultaneously makes the contact interface change [31]. The numerical results show a significant difference between the Hertz model and Flore model, as shown in Figure 3b, Figure 4b and Figure 5b. This is because the coefficient of restitution causes the Flore model to have smaller results than the Hertz model [1,14,15]. The Flore model has a larger contact force and smaller relative deformation than the experimental results, as shown in Figure 4 and Figure 5. Generally, a structure with higher rigidity is less likely to deform under the same test parameters, and its contact force will also be larger [34].
The contact force and relative deformation of the Fore model and the experiment all form a closed ring, and its area decreases as the impact mass increases, shown in Figure 6. It represents the work carried out by the damping force, that is, the energy dissipation during the impact [1]. The Hertz model does not take into consideration that the damping force, the contact force-relative deformation curve, is a line. In summary, the Flore model, which considers energy dissipation, fits the experimental results better than the Hertz model. However, there is still a significant difference between the Flore model and the impact wear test.

3.2. A Revised Contact Force Model for Thin-Walled Tubes

The current contact force models are developed based on the ideal impact of two solid spheres (Figure 2b). The sample used in this study is a thin-walled TP316H steel tube (wall thickness is 1.2 mm). The rigidity of a solid sphere is greater than that of a thin-walled tube structure of the same material. The dynamic response of the impact interface is not only related to the material but also the structure. A material with a greater hardness [34] or a structure with a greater rigidity (such as a tube with a larger wall thickness or a tube with a smaller span) has a greater contact force and a smaller relative deformation [35,36].
The bending stiffness of a solid cylinder or a tube can be expressed as [37]:
M   =   E I ,
where E is Young’s modulus, and I is the section second moment of inertia. For a solid cylinder, I can be expressed as [37]:
I 1   =   π D 4 4 ,
where D is the diameter of the cylinder, and for a tube, I can be expressed as [38]:
I 2   =   π ( D 4     d 4 ) 4 ,
where D, d is the outer and inner diameter of the tube.
The bending stiffness of the tube with the same outer diameter is smaller than that of a solid cylinder. A bending stiffness coefficient is introduced to evaluate this difference and is expressed as:
ε   =   M 2 M 1   =   D 4     d 4 D 4 ,
ε is a constant and is related to the cross-sectional size of the tube.
The right side of the Flore model can be applied to the tube by multiplying the bending stiffness coefficient, the revised model is expressed as:
F N   =   ε K δ n ( 1   +   8 1     c r 5 c r δ ˙ δ ˙ ) ,
Figure 7 shows the comparison between the revised model and the experimental results under varying impact mass. The contact force and relative deformation of the original Flore model are quite different from the experimental results (Figure 4 and Figure 5), whereas the revised model is almost the same as the experimental results, as shown in Figure 7a,c. Meanwhile, the velocity change curve of the revised model is also consistent with the experiment and Flore model. It shows that the revised model makes up the difference and maintains the consistency between the original Flore model and the experiment.
To further analyze the old and revised models and experiment results, Figure 8 and Figure 9 show the two critical parameters in impact wear: peak contact force and maximum relative deformation. The Hertz model that does not consider energy dissipation has the largest peak contact force, and the Flore model that considers energy dissipation is second, as shown in Figure 8. The maximum relative deformation is smallest in the Flore model, followed by the Hertz model, as shown in Figure 9. Since the model does not consider the influence of the structure on the overall stiffness, it has a greater stiffness, that is, a greater ability to resist deformation. Therefore, the contact force of the model is larger than that of the experimental results, and the relative deformation is smaller, as shown in Figure 8 and Figure 9. The revised model is very close to the experiment in terms of peak contact force and maximum relative deformation, as shown in Figure 8 and Figure 9. It shows that the revised model, which considers the structure of the impact body, is more suitable for analyzing the impact wear of thin-walled tubes than the Flore model.
With the same initial kinetic energy, the peak contact force and the maximum relative deformation of the model and test results gradually decrease with the increase in the impact mass, as shown in Figure 8 and Figure 9. The contact time also decreases gradually, as shown in Figure 4, Figure 5 and Figure 7a–c. Inertia is the ability of an object to maintain a state of motion, and the magnitude of inertia is reflected by mass [39]. The increase in mass increases the difficulty of changing its state of motion, and the time for the compression and restitution period also increases. This results in a reduction in the slopes of the force-time and relative deformation-time curves during the compression and restitution periods, as shown in Figure 4, Figure 5 and Figure 7a–c. At the end of the compression period (i.e., exhaustion of the initial kinetic energy), the peak point of the curve decreases as the slope decreases. This is consistent with other reported results [34,39].

3.3. Impact Wear Behavior of TP316H

Figure 10 shows the cross-sectional profile and 3D morphology of the wear scar under varying impact mass. Under different impact masses, the wear area is visible and significantly lower than the substrate. The surface of the wear scar is relatively smooth, and the section has some tiny jagged contours (Figure 10), which may be the accumulation of wear debris or the loss of material [35,36]. As the impact mass increases, the wear depth, area, and volume gradually decrease, as shown in Figure 10 and Figure 11.
Figure 12 shows the SEM images of the wear scar with an impact mass of 900 g. There is an accumulation of wear debris, delamination, and generation of cracks near the boundary of the wear scar, as shown in Figure 12b. Material plastic flow and obvious material spalling pits are also found inside of the wear scar, as shown in Figure 12c,d. The SEM cross-section micrograph of the wear scar, made near the white line in Figure 12a, is shown in Figure 13. The enlarged view of the cross-section clearly shows the accumulated wear debris and some material spalling pits, which are consistent with the cross-sectional profile of Figure 10 and the surface topography of Figure 12. Elemental analysis found that the area where the wear debris accumulated was rich in oxygen, indicating that the impact wear process was accompanied by oxidation of the material [34], as shown in Figure 13c.

4. Discussion

4.1. The Validity of the Revised Model

The existing contact force model is based on the impact between two solid spheres, and although it considers the effect of the properties of the material, it ignores the effect of structural stiffness on the impact system. The overall stiffness of the two-body impact system is closely related to its structure. The stiffness determines its ability to resist deformation, and the structure with higher stiffness shows a larger peak contact force in the impact test under the same parameters [35,36]. The contact force model is an expression of the contact force, so the most direct way to reduce the contact force is to multiply it by a factor (less than 1). This coefficient must be related to the structure to have broad applicability. Bending stiffness represents the ability of a part to resist bending deformation, and it is a good way to explain the difference between a thin-walled tube and a solid cylinder. As described in Section 3.2, the bending stiffness coefficient is the ratio of the bending stiffness of the thin-walled tube to the cylinder. The Flore model multiplied by the bending stiffness coefficient is in good agreement with the experimental results, as shown in Figure 7.
In addition, to further verify the applicability of the revised model, we conducted another set of experiments. The impact mass remains at 700 g and the initial impact kinetic energy changes (2.24, 350, and 5.4 mJ). The experimental results are also in good agreement with the revised model results, as shown in Figure 14. It shows that the revised model has a wide range of validity for the impact wear of thin-walled tubes. The maximum peak contact force and the maximum relative deformation increase with the initial impact kinetic energy, which is consistent with other studies [34,35,36,39].

4.2. Energy Dissipation during Impact

The impact wear process is inevitably accompanied by energy dissipation. Evaluating impact wear by calculating the energy loss during impact has been widely used, and the wear volume is proportional to the energy loss [6,35,40]. There are two methods to calculate the energy dissipation, one is to use the kinetic energy theorem, and the formula can be expressed as [22]:
Δ E   =   1 2 m v 1 2     v 2 2 ,
where m is the impact mass, v1, v2 are the pre-impact and post-impact velocity, respectively.
Another method is to calculate the mechanical work carried out by the damping force, which can be written as [33]:
Δ E   =   χ δ 1.5 δ ˙ d δ ,
That is the area enclosed by the contact force-relative deformation curve, as shown in Figure 1b.
Figure 15 shows the energy dissipation of the experiment and different contact force models under varying impact mass. Since the restitution coefficient in the model comes from the experimental results, the energy dissipation of the revised model, the Flore model, and the experimental results are relatively consistent, and the work carried out by the damping force is almost equal to the kinetic energy dissipation, as shown in Figure 15. This further verifies the correctness of the Flore model and the revised model. As the impact mass increases, the lost energy decreases gradually. Most of the energy lost is used for wear and plastic deformation of material [31,34], and a small part is lost in other forms such as sound and heat [17]. On the one hand, the coefficient of restitution increases with the increase in the impact mass (Table 3), then the kinetic energy loss decreases gradually [14,15]. The increase in mass (inertia) increases the difficulty for the impact block to change its state of motion [39]. When the initial kinetic energy is the same, with the rise in impact mass, the time of the contact force in the compression period becomes longer, and its peak force decreases. Thus, the area enclosed by the contact force and relative deformation will decrease, and the work carried out by the damping force decreases. On the other hand, as the impact mass increases, the momentum loss also decreases, as shown in Figure 15. The loss of momentum represents the product of the contact force and time during the impact process, that is, the impulse [41,42]. As the impact mass increases, the decrease in impulse is also responsible for the decrease in energy loss. In summary, with the rise in impact mass, the energy loss in the impact process decreases, resulting in less energy being used for wear, so the wear volume decreases. This further explains the results in Section 3.3.

5. Conclusions

In this study, the dynamic response and wear behavior of the impact interface of thin-walled TP316H steel were studied by combining the contact force model and experiment, and the main conclusions are as follows:
(a)
The revised contact force model is in good agreement with the experimental results and is more suitable for studying the dynamic response of the thin-walled tube impact interface than the existing model. This provides a reference for applying the contact force model to impact wear tests.
(b)
With the rise in impact mass, the coefficient of restitution increases from 0.65 to 0.78, whereas the peak contact force and maximum relative deformation decrease, and the momentum loss also decreases.
(c)
When the kinetic energy is equal, the energy dissipation and wear volume decrease as the impact mass increases. The impact wear mechanism of TP316H steel is mainly spalling, delamination, plastic deformation, and oxidative wear.

Author Contributions

Conceptualization, Q.-h.Y.; Data curation, L.-W.W.; Formal analysis, S.-L.M.; Investigation, F.Z.; Methodology, Q.-h.Y.; Project administration, S.-L.M.; Software, F.Z.; Validation, Z.-B.C.; Writing—original draft, X.-d.C.; Writing—review and editing, K.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Sichuan Science and Technology Planning Project: 2022JDJQ0019.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Skrinjar, L.; Slavič, J.; Boltežar, M. A review of continuous contact-force models in multibody dynamics. Int. J. Mech. Sci. 2018, 145, 171–187. [Google Scholar] [CrossRef]
  2. Lewis, R. A modelling technique for predicting compound impact wear. Wear 2007, 262, 1516–1521. [Google Scholar] [CrossRef]
  3. Zhong, L.; Feng, G.-H.; Zhu, H.-M.; Wang, G.-R.; Chen, W.-B. Structural Optimization Analysis and Erosion Resistance Study of Frac Ball Seat. Surf. Technol. 2021, 50, 213–228. [Google Scholar]
  4. Li, L.; Ge, P.-Q. Analysis of SPH and FEM Coupling Simulation of Abrasive Grain Scratching Single Crystal SiC. Surf. Technol. 2021, 50, 44–53. [Google Scholar]
  5. Cai, Z.B.; Li, Z.Y.; Yin, M.G.; Zhu, M.H.; Zhou, Z.R. A review of fretting study on nuclear power equipment. Tribol. Int. 2020, 144, 106095. [Google Scholar] [CrossRef]
  6. Souilliart, T.; Rigaud, E.; Le Bot, A.; Phalippou, C. Energy-based wear law for oblique impacts in dry environment. Tribol. Int. 2017, 105, 241–249. [Google Scholar] [CrossRef] [Green Version]
  7. Zalzalah, M.; Lewis, R.; Slatter, T. A new predictive model for normal and compound impact wear. Wear 2021, 480, 203954. [Google Scholar] [CrossRef]
  8. Johnson, K.L. One hundred years of Hertz contact. Proc. Inst. Mech. Eng. 1982, 196, 363–378. [Google Scholar] [CrossRef]
  9. Goldsmith, W. Impact: The Theory and Physical Behavior of Colliding Solids; Edward Arnold: Manchester, CT, USA, 1960. [Google Scholar]
  10. Hunt, K.H.; Crossley, F.R.E. Coefficient of restitution interpreted as damping in vibroimpact. J. Appl. Mech. Trans. ASME 1975, 42, 440–445. [Google Scholar] [CrossRef]
  11. Marhefka, D.W.; Orin, D.E. A compliant contact model with nonlinear damping for simulation of robotic systems. IEEE Trans. Syst. Man Cybern. Paart A-Syst. Hum. 1999, 29, 566–572. [Google Scholar] [CrossRef]
  12. Dubowsky, S.; Gardner, T. Design and analysis of multilink flexible mechanisms with multiple clearance connections. J. Eng. Ind.-Trans. ASME 1977, 99, 88–96. [Google Scholar] [CrossRef]
  13. Gonthier, Y.; McPhee, J.; Lange, C.; Piedboeuf, J.-C. A regularized contact model with asymmetric damping and dwell-time dependent friction. Multibody Syst. Dyn. 2004, 11, 209–233. [Google Scholar] [CrossRef]
  14. Alves, J.; Peixinho, N.; da Silva, M.T.; Flores, P.; Lankarani, H.M. A comparative study of the viscoelastic constitutive models for frictionless contact interfaces in solids. Mech. Mach. Theory 2015, 85, 172–188. [Google Scholar] [CrossRef]
  15. Machado, M.; Moreira, P.; Flores, P.; Lankarani, H.M. Compliant contact force models in multibody dynamics: Evolution of the Hertz contact theory. Mech. Mach. Theory 2012, 53, 99–121. [Google Scholar] [CrossRef]
  16. Pereira, C.M.; Ramalho, A.L.; Ambrósio, J.A. A critical overview of internal and external cylinder contact force models. Nonlinear Dyn. 2010, 63, 681–697. [Google Scholar] [CrossRef] [Green Version]
  17. Flores, P.; Machado, M.; Silva, M.T.; Martins, J.M. On the continuous contact force models for soft materials in multibody dynamics. Multibody Syst. Dyn. 2010, 25, 357–375. [Google Scholar] [CrossRef]
  18. Lee, T.W.; Wang, A. On the dynamics of intermittent-motion mechanisms. Part 1: Dynamic model and response. J. Mech. Trans. Automation. 1983, 105, 534–540. [Google Scholar] [CrossRef]
  19. Gilardi, G.; Sharf, I. Literature survey of contact dynamics modelling. Mech. Mach. Theory 2002, 37, 1213–1239. [Google Scholar] [CrossRef]
  20. Souilliart, T.; Rigaud, E.; Le Bot, A.; Phalippou, C. An energy approach for impact wear in water environment. Wear 2017, 376–377, 738–746. [Google Scholar] [CrossRef] [Green Version]
  21. Hong, S.M.; Kim, I.S. Impact fretting wear of alloy 690 tubes at 25 °C and 290 °C. Wear 2005, 259, 356–360. [Google Scholar] [CrossRef]
  22. Cai, Z.; Chen, Z.; Sun, Y.; Jin, J.; Peng, J.; Zhu, M. Development of a novel cycling impact–sliding wear rig to investigate the complex friction motion. Friction 2017, 7, 32–43. [Google Scholar] [CrossRef] [Green Version]
  23. Kayali, Y. The corrosion and wear behavior of TiN and TiAlN coated AISI 316 L stainless steel. Prot. Met. Phys. Chem. Surf. 2014, 50, 412–419. [Google Scholar] [CrossRef]
  24. Li, H.; Ramezani, M.; Li, M.; Ma, C.; Wang, J. Tribological performance of selective laser melted 316L stainless steel. Tribol. Int. 2018, 128, 121–129. [Google Scholar] [CrossRef]
  25. Chen, X.D.; Wang, L.W.; Yang, L.Y.; Tang, R.; Cai, Z.B. Effects of Temperature on the Fretting Wear Behavior of 2.25Cr-1Mo Tubes against Gr5C12 Rods. Materials 2020, 13, 3388. [Google Scholar] [CrossRef]
  26. Yu, Y.Q.; Chen, X.D.; Yang, L.Y.; Guan, X.M.; Cai, Z.B.; Wang, L.W. Investigation on fretting wear behavior of 2.25 Cr-1Mo tube in water at various temperatures. Wear 2021, 476, 203727. [Google Scholar] [CrossRef]
  27. Čiripová, L.; Falat, L.; Ševc, P.; Vojtko, M.; Džupon, M. Ageing effects on room-temperature tensile properties and fracture behavior of quenched and tempered T92/TP316H dissimilar welded joints with Ni-based weld metal. Metals 2018, 8, 791. [Google Scholar] [CrossRef] [Green Version]
  28. Blach, J.; Falat, L. The influence of thermal exposure and hydrogen charging on the notch tensile properties and fracture behaviour of dissimilar T91/TP316H weldments. High Temp. Mater. Process. 2014, 33, 329–337. [Google Scholar] [CrossRef]
  29. Chen, X.-D.; Feng, S.; Wang, L.-W.; Tang, R.; Zhang, F.; Ming, S.-L.; Cai, Z.-B. Effect of QPQ on the fretting wear behavior of TP316H steel at varying temperatures in liquid sodium. J. Nucl. Mater. 2022, 562, 153583. [Google Scholar] [CrossRef]
  30. Cai, Z.-B.; Chen, X.-D.; Rui, T.; Wang, L.-W.; Liu, X.-H.; Zhang, F. Investigation of the fretting corrosion mechanism of QPQ-treated TP316H steel in liquid sodium at 450 °C. Corros. Sci. 2022, 201, 110282. [Google Scholar] [CrossRef]
  31. Chen, X.-D.; Wang, L.-W.; Yang, L.-Y.; Tang, R.; Yu, Y.-Q.; Cai, Z.-B. Investigation on the impact wear behavior of 2.25 Cr–1Mo steel at elevated temperature. Wear 2021, 476, 203740. [Google Scholar] [CrossRef]
  32. Engel, P.; Bayer, R. The wear process between normally impacting elastic bodies. J. Lubr. Tech. 1974, 96, 595–604. [Google Scholar] [CrossRef]
  33. Lankarani, H.M.; Nikravesh, P.E. A contact force model with hysteresis damping for impact analysis of multibody systems. J. Mech. Des. 1990, 112, 369–376. [Google Scholar] [CrossRef]
  34. Wang, Z.; Cai, Z.-B.; Chen, Z.-Q.; Sun, Y.; Zhu, M.-H. Low-Velocity Impact Wear Behavior of Ball-to-Flat Contact Under Constant Kinetic Energy. J. Mater. Eng. Perform. 2017, 26, 5669–5679. [Google Scholar] [CrossRef]
  35. Sun, Y.; Cai, Z.-B.; Chen, Z.-q.; Qian, H.; Tang, L.-C.; Xie, Y.; Zhou, Z.-R.; Zhu, M.-h. Impact fretting wear of Inconel 690 tube with different supporting structure under cycling low kinetic energy. Wear 2017, 376–377, 625–633. [Google Scholar] [CrossRef]
  36. Lin, Y.-W.; Cai, Z.-B.; Chen, Z.-Q.; Qian, H.; Tang, L.-C.; Xie, Y.-C.; Zhu, M.-H. Influence of diameter–thickness ratio on alloy Zr-4 tube under low-energy impact fretting wear. Mater. Today Commun. 2016, 8, 79–90. [Google Scholar] [CrossRef] [Green Version]
  37. Berinskii, I.; Krivtsov, A.; Kudarova, A. Bending stiffness of a graphene sheet. Phys. Mesomech. 2014, 17, 356–364. [Google Scholar] [CrossRef]
  38. Shadmehri, F.; Derisi, B.; Hoa, S.V. On bending stiffness of composite tubes. Compos. Struct. 2011, 93, 2173–2179. [Google Scholar] [CrossRef]
  39. Wang, Z.; Cai, Z.-B.; Sun, Y.; Peng, J.-F.; Zhu, M.-H. Low velocity impact wear behavior of MoS2/Pb nanocomposite coating under controlled kinetic energy. Surf. Coat. Technol. 2017, 326, 53–62. [Google Scholar] [CrossRef]
  40. Abdollah, M.F.B.; Yamaguchi, Y.; Akao, T.; Inayoshi, N.; Tokoroyama, T.; Umehara, N. The Effect of Maximum Normal Impact Load, Absorbed Energy, and Contact Impulse, on the Impact Crater Volume/Depth of DLC Coating. Tribol. Online 2011, 6, 257–264. [Google Scholar] [CrossRef] [Green Version]
  41. Brach, R.M. Impact dynamics with applications to solid particle erosion. Int. J. Impact Eng. 1988, 7, 37–53. [Google Scholar] [CrossRef]
  42. Brach, R.M. Classical planar impact theory and the tip impact of a slender rod. Int. J. Impact Eng. 1993, 13, 21–33. [Google Scholar] [CrossRef]
Figure 1. Test equipment (a): (1) voice coil motor, (2) damping punch, (3) guide rail, (4) impact block, (5) force sensor, (6) fixture, (7) sample, (8) tri-bro pair, (9) velocity sensor; and the dynamic response of an impact event: (b) contact force versus deformation, (c) contact force and velocity versus time.
Figure 1. Test equipment (a): (1) voice coil motor, (2) damping punch, (3) guide rail, (4) impact block, (5) force sensor, (6) fixture, (7) sample, (8) tri-bro pair, (9) velocity sensor; and the dynamic response of an impact event: (b) contact force versus deformation, (c) contact force and velocity versus time.
Materials 15 02881 g001
Figure 2. Schematic diagram of impact wear (a), idealized contact-impact between two spheres (b), and equivalent system (c).
Figure 2. Schematic diagram of impact wear (a), idealized contact-impact between two spheres (b), and equivalent system (c).
Materials 15 02881 g002
Figure 3. Relative velocity vs time under varying impact mass.
Figure 3. Relative velocity vs time under varying impact mass.
Materials 15 02881 g003
Figure 4. Contact force vs time under varying impact mass.
Figure 4. Contact force vs time under varying impact mass.
Materials 15 02881 g004
Figure 5. Relative deformation vs time under varying impact mass.
Figure 5. Relative deformation vs time under varying impact mass.
Materials 15 02881 g005
Figure 6. Contact force vs relative deformation under varying impact mass.
Figure 6. Contact force vs relative deformation under varying impact mass.
Materials 15 02881 g006
Figure 7. Dynamic response of impact wear with varying impact mass, experimental and revised model results.
Figure 7. Dynamic response of impact wear with varying impact mass, experimental and revised model results.
Materials 15 02881 g007
Figure 8. Peak contact force during a contact−impact event with varying impact mass.
Figure 8. Peak contact force during a contact−impact event with varying impact mass.
Materials 15 02881 g008
Figure 9. Maximum deformation during a contact-impact event with varying impact mass.
Figure 9. Maximum deformation during a contact-impact event with varying impact mass.
Materials 15 02881 g009
Figure 10. Cross-section profiles and 3D topography of worn scars of varying impact mass.
Figure 10. Cross-section profiles and 3D topography of worn scars of varying impact mass.
Materials 15 02881 g010
Figure 11. Wear volume, area, and depth of TP316H.
Figure 11. Wear volume, area, and depth of TP316H.
Materials 15 02881 g011
Figure 12. SEM micrographs of wear scar (a) and enlarged view (bd) under the impact mass of 900 g.
Figure 12. SEM micrographs of wear scar (a) and enlarged view (bd) under the impact mass of 900 g.
Materials 15 02881 g012
Figure 13. SEM cross-section micrographs of wear scar (ae) and EDS (c) under the impact mass of 900 g.
Figure 13. SEM cross-section micrographs of wear scar (ae) and EDS (c) under the impact mass of 900 g.
Materials 15 02881 g013
Figure 14. Dynamic response of contact-impact: experimental and numerical results of varied impact energy.
Figure 14. Dynamic response of contact-impact: experimental and numerical results of varied impact energy.
Materials 15 02881 g014
Figure 15. Energy dissipation and momentum loss during a contact-impact event with varying impact mass.
Figure 15. Energy dissipation and momentum loss during a contact-impact event with varying impact mass.
Materials 15 02881 g015
Table 1. Mechanical properties of tribo-pairs.
Table 1. Mechanical properties of tribo-pairs.
MaterialsHV(0.1)Poisson’s RatioYoung’s ModulusSurface Roughness
TP316H174.630.29112.98 GpaRa = 0.12 μm
316H172.550.2889.11 GpaRa = 0.13 μm
Table 2. Experimental parameters for the impact wear test.
Table 2. Experimental parameters for the impact wear test.
Impact Wear TestSample Geometry
Tube materialTP316HTube
Rod material316HOuter diameter16 mm
TemperatureRoom temperatureInner diameter13.6 mm
Impact mass500 g, 600 g, 700 g, 800 g, 900 gLength30 mm
Impact energy3.50 mJRod
Impact frequency3 HzDiameter10 mm
Number of cycles104Length20 mm
Table 3. Experimental results and calculation parameters under varying impact mass.
Table 3. Experimental results and calculation parameters under varying impact mass.
Test VariablesThe Experimental DataThe Calculated Data
Pre-Impact Velocity (v1)
mm/s
Post-Impact Velocity (v2)
mm/s
Kcr χ
500 g117.976.44.0 × 1090.652.9 × 1010
600 g107.973.80.682.8 × 1010
700 g99.673.00.732.4 × 1010
800 g93.070.00.752.3 × 1010
900 g89.769.80.782.0 × 1010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, X.-d.; Wang, L.-W.; Yu, Q.-h.; Zhang, F.; Mo, K.; Ming, S.-L.; Cai, Z.-B. Experimental and Numerical Analysis on the Impact Wear Behavior of TP316H Steel. Materials 2022, 15, 2881. https://doi.org/10.3390/ma15082881

AMA Style

Chen X-d, Wang L-W, Yu Q-h, Zhang F, Mo K, Ming S-L, Cai Z-B. Experimental and Numerical Analysis on the Impact Wear Behavior of TP316H Steel. Materials. 2022; 15(8):2881. https://doi.org/10.3390/ma15082881

Chicago/Turabian Style

Chen, Xu-dong, Li-Wen Wang, Qi-hang Yu, Fan Zhang, Kun Mo, Shi-Lin Ming, and Zhen-Bing Cai. 2022. "Experimental and Numerical Analysis on the Impact Wear Behavior of TP316H Steel" Materials 15, no. 8: 2881. https://doi.org/10.3390/ma15082881

APA Style

Chen, X. -d., Wang, L. -W., Yu, Q. -h., Zhang, F., Mo, K., Ming, S. -L., & Cai, Z. -B. (2022). Experimental and Numerical Analysis on the Impact Wear Behavior of TP316H Steel. Materials, 15(8), 2881. https://doi.org/10.3390/ma15082881

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop