Research Article | Open Access

Volume 2021 |Article ID 6626727 | https://doi.org/10.1155/2021/6626727

Donghua Yin, Xueliang Zhang, Yonghui Chen, Guosheng Lan, Yanhui Wang, Shuhua Wen, "Statistical Model of Normal Contact Stiffness of Fixed Joint Surface during Unloading after First Loading", Advances in Materials Science and Engineering, vol. 2021, Article ID 6626727, 8 pages, 2021. https://doi.org/10.1155/2021/6626727

Revised07 Feb 2021
Accepted13 Feb 2021
Published29 Apr 2021

#### 1. Introduction

The machine tool is the mother of manufacturing industry, which is assembled by various components. There are some parts that contact each other, namely joint surface. Among them, the fixed joint surface is one of the widely existing joint surfaces. Research showed that  the ratio of contact stiffness of fixed joint surface to total stiffness of the whole machine was more than 60%. Therefore, it is of great significance to establish a more accurate contact stiffness model for the analysis of static and dynamic characteristics of the machine tool structure.

#### 2. Model of Normal Contact Stiffness of a Single Asperity

Figure 1 presents the deformation of a single asperity before and after contact with a rigid plane, where the dashed lines show the situation before deformation. The displacement of the rigid plane is the deformation of the asperity under the applied of the normal load p. is the curvature radius of the asperity summit.

##### 2.1. Model of Normal Contact Stiffness of a Single Asperity during Loading

The critical deformation of a single asperity , when it transforms from the elastic to the elastic-plastic deformation regime, is given by where is the hardness of the softer material, which is related to its yield strength by , is the hardness coefficient in the form , where is Poisson’s ratio of the softer material, And , the comprehensive elastic modulus, is given in , where and and and are the elastic modulus and Poisson’s ratios of the two materials, respectively. In this model, the rigid plane is smooth by , so the comprehensive elastic modulus can be simplified as .

When , a single asperity deforms elastically. According to Hertz theory , the normal contact load during loading can be expressed as follows:

Therefore, the critical contact load of a single asperity that marks the transition from the elastic to the elastic-plastic deformation regime can be expressed as follows:

According to equation (2), the normal contact stiffness of a single asperity during loading is given by

When , elastic-plastic deformation of a single asperity occurs. The normal contact load during loading can be expressed as follows:

Similarly, the normal contact stiffness of a single asperity during this loading is obtained as follows:

When , a single asperity has a completely plastic deformation, in which there is no stiffness.

##### 2.2. Model of Normal Contact Stiffness of a Single Asperity during Unloading

When , the normal contact load and the normal contact stiffness of the asperity during unloading can be expressed as follows:

According to , the relationship between normal contact load and unloading deformation of asperity iswhere the index of the plastic stage is given by .

The residual deformation and the radius of residual nonuniform curvature depend on the maximum load at the beginning of the unloading . According to relations (5) and (6) between contact load and deformation, the maximum load can be deduced as follows:

By substituting equations (12), (13) into equation , respectively, and making further differentiation, the unloading stiffness of an elastic-plastic deformable asperity can be deduced as follows:

In the same way, the plastically deformable asperity does not recover. There is no contact stiffness during unloading.

In , the ratio relation was given as follows:

#### 3. Statistical Model of Normal Contact Stiffness of Joint Surface

##### 3.1. Statistical Model of Normal Contact Stiffness of Joint Surface during Loading

Based on Greenwood and Williamson’s model (GW model), this paper assumes that there is no interaction between asperities, and all deformation is limited to the contacting asperities. So, the fixed joint surface is equivalent to the contact between a rough surface and a smooth rigid plane, as shown in Figure 3. is the height of asperities, is the distance between the mean of asperity heights and the rigid plane, is the distance between the mean of surface heights and the rigid plane, and is the distance between the mean of asperity heights and the mean of surface heights, satisfying the relation . The rough surface is isotropic, and its morphology is defined by three independent parameters: the area density of asperities , the ratio of the standard deviation of asperity heights to the standard deviation of surface heights , and the radius of curvature of asperity summit .

The relation of the ratio can be expressed aswhere is a dimensionless surface roughness parameter in the form:

Assuming that there are asperities on the nominal contact area , the expected number of contact asperities on the joint surface is given bywhere is the probability density function of the normal distribution of asperity heights.

The distribution function of dimensionless asperity heights is described by a dimensionless Gaussian standard probability density function in the form:

The dimensionless distance between the mean of asperity heights and the mean of surface heights is given by 

The random dimensionless interference of a single asperity can be expressed as follows:

In this paper, the plastic index form proposed by GW is adopted:

According to equation (1), the normal contact stiffness of the joint surface during loading is

The dimensionless form of equation (24) is

The normal contact stiffness of joint surface during unloading can be expressed as

The dimensionless form of the above equation iswhere .

#### 4. Simulation and Result Analysis of the Model

It can be seen from equation (25) that the dimensionless loading normal contact stiffness is a function of the standard deviation of surface heights , the radius of curvature at the initial summit of asperity , the dimensionless surface roughness parameter , the dimensionless surface mean separation , and so on. It can be seen from equation (27) that the dimensionless unloading normal contact stiffness is a function of the standard deviation of surface heights , the radius of curvature at the initial summit of asperity , the dimensionless surface roughness parameter , the dimensionless surface mean separation , the dimensionless residual deformation , and so on. And, it is not affected by residual nonuniform curvature radius . In the simulation analysis, the parameters are given such as the elastic modulus , Poisson’s ratio , the hardness , the radius of curvature at the initial summit of asperity , and the dimensionless surface roughness parameters and  (shown in Table 1). Equations (25) and (27) are simulated by using the data of each variable, and the corresponding results are shown in Figure 4.

 Number 1 0.0339 2 0.0476 3 0.0541 4 0.0601

It can be seen from Figure 4 that the normal contact stiffness of joint surface during loading and unloading is a nonlinear function of the mean surface separation and decreases with the increase of the mean surface separation. When the plastic index is smaller, the contact between asperities is more elastic, so the normal contact stiffness curves of joint surface during loading and unloading are close. When the plastic index is larger, the plastic deformation cannot recover due to the large proportion of the plastically deformed asperities, so the normal contact stiffness decreases rapidly during unloading.

#### 5. Conclusions

Unfortunately, our paper did have some limitations and shortcomings, which will be verified by supplementary experiments in the future.

#### Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

#### Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

#### Acknowledgments

This work was supported by Shanxi Provincial Natural Science Foundation of China (Grant no. 201901D111248), Shanxi Provincial “1331” Engineering Key Discipline Construction Project of China, and Shanxi Provincial Graduate Education Innovation Project of China (Grant No. 2020BY112).

1. Y. Ren and C. F. Beards, “Identification of “effective” linear joints using coupling and joint identification techniques,” Journal of Vibration and Acoustics, vol. 120, no. 2, pp. 331–338, 1998. View at: Publisher Site | Google Scholar
2. Y.-Y. Lin and C. Y. Hui, “Mechanics of contact and adhesion between viscoelastic spheres: an analysis of hysteresis during loading and unloading,” Journal of Polymer Science Part B: Polymer Physics, vol. 40, no. 9, pp. 772–793, 2002. View at: Publisher Site | Google Scholar
3. S. L. Yan and L. Y. Li, “Finite element analysis of cyclic indentation of an elastic-perfectly plastic half-space by a rigid sphere,” Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, vol. 217, no. 5, pp. 505–514, 2003. View at: Publisher Site | Google Scholar
4. R. E. Jones, “Models for contact loading and unloading of a rough surface,” International Journal of Engineering Science, vol. 42, no. 17-18, pp. 1931–1947, 2004. View at: Publisher Site | Google Scholar
5. I. Etsion, Y. Kligerman, and Y. Kadin, “Unloading of an elastic-plastic loaded spherical contact,” International Journal of Solids and Structures, vol. 42, no. 13, pp. 3716–3729, 2005. View at: Publisher Site | Google Scholar
6. Y. Kadin, Y. Kligerman, and I. Etsion, “Unloading an elastic-plastic contact of rough surfaces,” Journal of the Mechanics and Physics of Solids, vol. 54, no. 12, pp. 2652–2674, 2006. View at: Publisher Site | Google Scholar
7. Y. Zait, Y. Kligerman, and I. Etsion, “Unloading of an elastic-plastic spherical contact under stick contact condition,” International Journal of Solids and Structures, vol. 47, no. 7-8, pp. 990–997, 2010. View at: Publisher Site | Google Scholar
8. J. H. Zhao, S. Nagao, and Z. L. Zhang, “Loading and unloading of a spherical contact: from elastic to elastic-perfectly plastic materials,” International Journal of Mechanical Sciences, vol. 56, no. 1, pp. 70–76, 2012. View at: Publisher Site | Google Scholar
9. W. Fu, L. T. Lou, and Z. Q. Gao, “Theoretical model for the contact stiffness and damping of mechanical joint surface,” Journal of Mechanical Engineering, vol. 53, no. 9, pp. 073–082, 2017. View at: Publisher Site | Google Scholar
10. R. An, X. L. Zhang, and Y. S. Wang, “Statistical model of unloading normal contact stiffness of joint surface,” Journal of Taiyuan University of Science and Technology, vol. 38, no. 5, pp. 373–377, 2017. View at: Google Scholar
11. H. Nassiraei, “Local joint flexibility of CHS X-joints reinforced with collar plates in jacket structures subjected to axial load,” Applied Ocean Research, vol. 93, pp. 1–10, 2019. View at: Publisher Site | Google Scholar
12. H. Nassiraei, “Local joint flexibility of CHS T/Y-connections strengthened with collar plate under in-plane bending load: parametric study of geometrical effects and design formulation,” Ocean Engineering, vol. 202, pp. 1–11, 2020. View at: Publisher Site | Google Scholar
13. H. Nassiraei, “Geometrical effects on the LJF of tubular T/Y-joints with doubler plate in offshore wind turbines,” Ships and Offshore Structures, vol. 9, pp. 1–11, 2020. View at: Publisher Site | Google Scholar
14. L. Kogut and I. Etsion, “Elastic-plastic contact analysis of a sphere and a rigid flat,” Journal of Applied Mechanics, vol. 69, no. 5, pp. 657–662, 2002. View at: Publisher Site | Google Scholar
15. D. Tabor, The Hardness of Metals, Clarendon Press, Oxford, UK, 1951.
16. K. L. Johnson, Contact Mechanics, Cambridge University Press, Cambridge, UK, 1985.
17. L. Li, A. J. Cai, and L. G. Cai, “Micro-contact model of bolted-joint interface,” Journal of Mechanical Engineering, vol. 52, no. 7, pp. 205–212, 2016. View at: Publisher Site | Google Scholar