Next Article in Journal
Synthesis, X-ray Crystal Structure, and Photochromism of a Sandwich-Type Mono-Aluminum Complex Composed of Two Tri-Lacunary α-Dawson-Type Polyoxotungstates
Next Article in Special Issue
Surface Roughness Effects on the Broadband Reflection for Refractory Metals and Polar Dielectrics
Previous Article in Journal
Green Synthesis and Characterization of Pullulan Mediated Silver Nanoparticles through Ultraviolet Irradiation
Previous Article in Special Issue
Surface-Enhanced Raman Spectroscopy (SERS) Study Using Oblique Angle Deposition of Ag Using Different Substrates
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Experimental Verification and Comparative Analysis of Equivalent Methods on Metal’s Fixed Joint Interface

1
College of Intelligence Science and Technology, National University of Defense Technology, Changsha 410073, China
2
School of Mechanical Engineering, Hunan International Economics University, Changsha 410205, China
*
Author to whom correspondence should be addressed.
Materials 2019, 12(15), 2381; https://doi.org/10.3390/ma12152381
Submission received: 27 June 2019 / Revised: 22 July 2019 / Accepted: 24 July 2019 / Published: 26 July 2019
(This article belongs to the Special Issue Surface and Interface Engineering)

Abstract

:
In order to effectively improve the dynamic characteristics of the fixed metal joint interface, it is important to establish a correct equivalent model of the metal joint interface. In this paper, three equivalent methods for simulating the metal joint interface are analyzed, including the virtual material method, spring damping method, finite element method, and verification by modal experiment. First, according to the contact mechanics model of the constructed metal joint interface, the physical properties of the three-dimensional models of the fixed joint interface are assigned in the ANSYS software. Then, three methods are used for the modal analysis and compared with a modal experiment. The results show that the modal shapes of the three theoretical methods are consistent with those of the experimental modes. The first five natural frequencies obtained by the virtual material method are closest to the experimental natural frequencies, and the errors are within 10%. The errors of the other two methods are between 9% and 39%. Therefore, the virtual material method is a better equivalent method of the metal joint interface.

Graphical Abstract

1. Introduction

The metal surface machined by the machine tool is a rough surface composed of many microscopic irregular peaks and valleys. Studies have shown that [1,2,3,4], although the surface morphology produced by a machine tool shows complex and random features in local minute details, there are still certain rules to follow. The machined metal surface has self-affine fractal features that can be described by fractal geometry. Therefore, the metal’s fixed joint interface consists essentially of two rough surfaces with fractal properties. The development of contact mechanics for rough surfaces originated from the ideal smooth surface frictionless elastic contact theory proposed by Hertz in 1881, but contact mechanics developed rapidly within decades. Based on the Hertz contact theory, Majumdar and Bhushan [5] proposed the elastoplastic contact fractal geometry theory of rough surface with scale independence, the MB model, in 1991. The MB model uses the fractal parameters to quantitatively express the contact area and contact load of the metal’s fixed joint interface.
Through the study of the joint surface performance, scholars found that 50% of the system stiffness and more than 90% of the system damping come from the fixed joint interface [6]. Therefore, the selection of the exact equivalent method of the metal joint interface is of great significance for engineering applications. The precise equivalent method of the joint interface can solve the engineering practicality problem of the dynamic model. The current equivalent methods of fixed metal joint interface include the virtual material method (VMM), spring damping method (SDM) and finite element method (FEM). The virtual material method is the most popular equivalent method of the joint interface in recent years.
The VMM, which is also called the equivalent gap method [7], was proposed by Tian [8]. Because the joint surface has the characteristics of gap, friction, wear and hysteresis, the method is proposed from the microscopic point of view, considering the elastoplastic deformation of the real contact surface, the material variation characteristics, thermal characteristics and timeliness of the contact metal. The virtual material is composed of a rough metal surface and a medium of contact surface, which is the main deformation area. The method uses the fractal theory [9] to establish the virtual material model of the joint interface, and to calculate the parameters of the model from the microscopic point of view, such as the Poisson ratio, density and elastic modulus of the virtual material. The length and width of the virtual material are the length and width of the joint, respectively. The thickness of the virtual material is determined according to the actual gap between the two surfaces. The virtual material and the fixed joint are fixedly connected on both sides, and finally applied to the finite element analysis. Zhao [10] improved the VMM, and adopted the improved method to simulate the mechanical joint interface. As a result, the accuracy of the results was improved. Huang [11] applied the virtual material method to simulate the bolt joint interface in the modal analysis. The SDM is used to characterize the stiffness and damping between mechanical joints by providing spring damping elements on the nodes of joint interfaces [12] to achieve the equivalent mechanical joint interface. The selection of the nodes depends on the shape of the fixed joint interface, the connection mode of the fixed joint interface, and the area of the fixed joint interface. The modal analysis of the method in the free state proves the practicability and accuracy of the method.
The accurate analysis of the dynamic characteristics of the joint interface is of great significance for the optimization of the dynamic performance of the whole machine. Therefore, it is necessary to select a high-precision equivalent method of the joint interface. In this paper, the modal analysis of the joint interface is carried out by a modal experiment. The qualitative comparison of the mode shape and the quantitative comparison of the natural frequency are used to verify the feasibility of the joint equivalent method and to select a better equivalent method of the metal joint interface.

2. Materials and Methods

Two identical square plates are selected as experimental specimens. The specimen, made of 45#steel, has 200 mm length, 6 mm thickness, 7850 kg/m3 density, 3.2 µm surface roughness, 209 Gpa elastic modulus, and 0.3 Poisson ratio. The two plates are pressed by a pressure of 0.625 MPa, as shown in Figure 1. Three methods which are equivalent to the joint interface of the above conditions are as follows:

2.1. Virtual Material Method

The VMD divides the three-dimensional model of the two boards into upper, middle, and lower boards (numbers 1, 2, and 3, respectively) as shown in Figure 2. The thickness of the board 1 and board 3 is h 1 = h 2 = 5.6   mm . Their physical properties are unchanged. The board 2 is the virtual material. By theoretical calculation, the thickness of the board 2 is h = 0.8   mm [8], and the formula for calculating the elastic modulus of the virtual material can be written as [10]:
E = 2 D 4.5 E * D φ 1 D 2 G 1 D a L D 2 3 π 1 + D 2 ( ln   γ ) 1 2 A 0 ( a L 1 2 a C 1 2 )
where D and G are the fractal parameters of the surface profile of the joint interface, φ describes the domain extension factor for the micro-contact size distribution associated with D , γ is the scaling parameter, a L is the truncated area of the largest elastic micro-contact, a c is the critical truncated area demarcating the elastic and plastic deformation regimes, A 0 is the nominal contact area and E * is the equivalent elastic modulus. According to [8], D and G are determined by the roughness of the surface of the joint interface. According to [10], a L is only affected by the surface pressure exerted on the board when the physical properties of the material are unchanged. The formula of the equivalent elastic modulus can be written as [13]:
1 E * = ( 1 ν 1 2 ) E 1 + ( 1 ν 2 2 ) E 2
where ν 1 and E 1 are the Poisson’s ratio and elastic modulus of the board 1 and board 2, and ν 2 and E 2 are the Poisson’s ratio and elastic modulus of the board 3.
The normal stiffness value between the joint interfaces and the shear modulus of the virtual material can be described as [14]
K n = E A 0 / h
G τ = K τ h / A 0
where K τ is the tangential contact stiffness. According to [15], the ratio of the tangential contact stiffness to the normal contact stiffness is 0.25 to 0.35. Without a loss of generality, K τ / K n is taken as 0.35.
The Poisson’s ratio and density of the virtual material can be described, respectively, as [10,16]:
υ = E 2 G τ 1
ρ = ρ 1 h 1 + ρ 2 h 2 h 1 + h 2
where ρ 1 and ρ 2 are the density of board 1 and board 2, respectively. Through the above calculation, E = 0.3554   Gpa , ν = 0.429 and ρ = 7850   kg / m 3 .

2.2. Spring Damping Method and Finite Element Method

In the SDM model, a set of twelve spring damping units connects to the four corners of the joint interface, so that the joint interface is evenly stressed. The schematic diagram of the model is shown in Figure 3a. At each node, there are three sets of spring damping units. The normal and tangential total stiffness and damping of the spring are the normal and tangential total stiffness and damping of the joint, respectively. The value of each spring and damping can be expressed, respectively, as:
K i x = K i y = K τ / 4
K i z = K n / 4
C i x = C i y = C τ / 4
C i z = C n / 4
where K τ is the tangential contact stiffness, K n is the normal stiffness, C τ is the tangential damping and C n is the normal damping. According to [17], the stiffness value and damping value of the spring damping unit can be obtained: K x = K y = 8.53 e 7   N / m , K z = 2.44 e 8   N / m , C x = C y = 5.22 e 5   N · s / m and C z = 4.62 e 8   N · s / m .
The FEM connects the three-dimensional model of the joint interface through the bond connection method in the ANSYS software, as shown in Figure 3b, so that the two boards become a whole. The original physical properties of the board were used, and the finite element analysis was performed. Therefore, the elastic modulus E 0 , the Poisson ratio ν 0 and the density ρ 0 of the board of the FEM are as follows: E 0 = 209   Gpa , ν 0 = 0.3 and ρ 0 = 7850   kg / m 3 .

3. Simulation, Modal Experiment and Results

3.1. Modeling and Simulation

The above three models are modeled in three dimensions through the CATIA software. The model was imported into the ANSYS or ABAQUS software, the boundary condition was set to the free state, and the physical properties of the model were assigned. The mesh diagram of the model of the FEM is shown in Figure 4. The grid has a total of 129,955 nodes and 63,599 elements. The six holes in the Figure 4 are bolt holes. A force was applied around the bolt holes so that the pressure on the joint interface is equal to 0.625 MPa. Finally, the modal analysis is performed. The first five natural frequencies are shown in Table 1, and the mode shapes are shown in Table 2.

3.2. Modal Experiment

The modal experiment (ME) uses an exciter to excite and laser scans to acquire data. The natural frequencies and mode shapes of the test piece are calculated and analyzed based on the measured response data. By comparing the results obtained by finite element analysis and modal experiments, the optimal equivalent method of the joint interface can be obtained.
First, a torque wrench was used to fasten two identical boards together by bolts, as shown in Figure 5a. Surface pressure is applied to the plates through the bolts. In order to measure the natural frequencies of the joint interface, a boundary-less condition is selected. However, it cannot be fully realized in the experiment. Therefore, the experiment chooses two elastic ropes to hang the two boards in front of the exciter. The specimens are vibrated by the exciter. The laser is used to detect the modal shapes and natural frequencies of the specimen. The final data is collected on the computer for signal processing and identification. The natural frequencies and mode shapes of the first five orders are obtained. The experimental setup is shown in Figure 5b.

3.3. Results

The first five natural frequencies of the modal experiment are 660.9 Hz, 1132.8 Hz, 1325.0 Hz, 1564.1 Hz, 1564.1 Hz, 1635.1 Hz, respectively, as shown in Figure 6. The first five modal shapes obtained by the experiment are shown in Table 2.

4. Discussion

The comparison between the first five-order mode shapes of the modal experiment and those of the simulation is shown in Table 2. Compared with the experimental mode shape diagrams, the mode shapes of the three equivalent methods of the joint interface are roughly similar to the experimental mode shapes. However, the second-order mode shape obtained by the SDM is significantly different from the mode shape obtained by the experiment. Furthermore, the mode shapes obtained by simulation analysis are slightly different from those obtained by the experiment in the edge area of the mode shape, especially the fifth mode shape. The reason for this is that the test area is the red border area, as shown in Figure 5a, and the mode shapes outside the red frame cannot be measured. Therefore, it can be proved that the VMD and the FEM are more suitable for simulating the metal’s fixed joint interface.
The contrast and error of the first five natural frequencies obtained by simulation and modal experiment are shown in Table 3 and Figure 7. By comparing the first five natural frequencies obtained by the modal experiment, it is found that the error of the natural frequencies obtained by the VMM is the smallest, being within 10%. The errors of the natural frequencies obtained by the other two equivalent methods of the joint interface are between 9% and 39%. The results prove that the VMM is more realistic than the other two methods. It can be clearly seen from Figure 7 that, as the order increases, the difference between the natural frequency obtained by the simulation and the modal experiment is larger. Therefore, it can be explained that for a high-order modal simulation, the equivalent methods for simulating the metal joint interface will lose their authenticity.
The three equivalent methods of the joint interface all have disadvantages. The VMD and the SDM are the joint interface equivalent models established based on the fractal contact theory. By modelling the contact mechanics of the joint interface, the physical properties of the virtual material, and the stiffness and damping values of the spring damping unit, are obtained respectively. However, the fractal contact theory has certain limitations. The contact mechanics model ignores the interaction between the asperities. When the squeezing force is small, this assumption is approximately true. However, when the squeezing force is relatively large, the interaction between asperities cannot be ignored. Therefore, it will affect the establishment of the equivalent method of the fixed joint interface. The VMM is also affected by the thickness of the virtual material. If the thickness of the virtual material is too large or too small, the VMM will lose some authenticity. The SDM cannot accurately determine the number and position of the spring damping unit. Generally, the number and position of the spring damping unit are determined according to the number of nodes of the joint interface, but it also causes a certain error. For example, the second-order mode shape of the SDM does not match the experimental mode shape. The FEM ignores the tiny structure between the joint interface so that the modal simulation brings a big error. Additionally, the simulation environment cannot be exactly the same as the experiment. Therefore, there are certain errors in the equivalent methods of the joint interface.
By analyzing the error of the results obtained by the modal experiment and the modal simulation, the following aspects can be improved. First, in the basic theory part, it is necessary to consider the interaction between the asperities of the fixed joint interface to establish a more precise contact mechanics model. Second, the multi-variable modal simulation is carried out on the thickness of the virtual material and the number and position of the spring damping unit, and then compared with the modal experiment to obtain a more accurate equivalent method on the metal’s fixed joint interface.

5. Conclusions

In this paper, the modal analysis of the VMM, the FEM and the SDM was carried out by the finite element simulation software. The first five natural frequencies and mode shapes of the three joint interface equivalent models were obtained, and then the modal experiment was performed. Finally, the experimental results were compared with the results of the simulation. The conclusions can be summarized as follows:
  • The modal shapes of the modal experiment are roughly similar to the modal shapes obtained by the three equivalent methods of the joint interface. It is proved that all three methods can simulate the microscopic metal joint interface.
  • The natural frequencies obtained by three equivalent methods of the metal joint interface have obvious deviations from the results of the modal experiment. The error of the natural frequency value obtained by the VMM is the smallest, being within 10%. Obviously, it is lower than the errors of the other two methods, which proves that the VMM is the closest to the real situation among the three equivalent methods of the metal joint interface.

Author Contributions

Modal experiment was completed by R.H. and K.Z. under the supervision of G.L. and J.G. The simulation was completed by R.H. under the guidance of M.Z., R.H. and K.Z. wrote the paper.

Funding

This research was funded by the National Natural Science Foundation of China, Grant number 51775552 and Key Research Project of The Education Department of Hunan Province, Grant number 15A107.

Acknowledgments

This work was carried out during a project, funded by National Natural Science Foundation of China (Grant: 51775552) and Key Research Project of The Education Department of Hunan Province (Grant: 15A107)

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Majumdar, A.; Tien, C.L. Fractal Characterization and Simulation of Rough Surfaces. Wear 1990, 136, 313–327. [Google Scholar] [CrossRef]
  2. Gagnepain, J.J.; Roques-Carmes, C. Fractal Approach to Two-Dimensional and Three-Dimensional Surface Roughness. Wear 1990, 109, 114–119. [Google Scholar] [CrossRef]
  3. Lopez, J.; Hansali, G.; Zahouani, H. 3D Fractal-Based Characterization for Engineered Surface Topography. Int. J. Mach. Tool Manuf. 1995, 35, 211–217. [Google Scholar] [CrossRef]
  4. Marotta, E.E.; Fletcher, L.S.; Dietz, T.A. Thermal Contact Resistance Modeling of Non-Flat, Roughened Surfaces with Non-Metallic Coatings. J. Heat Transf. 2000, 123, 11–23. [Google Scholar] [CrossRef]
  5. Majumdar, A.; Bhushan, B. Role of Fractal Geometry in Roughness Characterization and Contact Mechanics of Surfaces. J. Tribol. 1990, 112, 205–216. [Google Scholar] [CrossRef]
  6. Ibrahim, R.A.; Pettit, C.L. Uncertainties and dynamic problems of bolted joints and other fasteners. J. Sound Vib. 2005, 279, 857–936. [Google Scholar] [CrossRef]
  7. Li, X.; Liang, Y.; Guo, H.; Ju, X.; Wen, B. Study on equivalent model of generalized clearance of joint surface. J. Vib. Eng. 2014, 27, 25–32. [Google Scholar]
  8. Tian, H.; Li, B.; Liu, H.; Mao, K.; Peng, F.; Huang, X. A new method of virtual material hypothesis-based dynamic modeling on fixed joint interface in machine tools. Int. J. Mach. Tool Manuf. 2011, 51, 239–249. [Google Scholar] [CrossRef]
  9. Yuan, Y.; Gan, L.; Liu, K.; Yang, X. Elastoplastic Contact Mechanics Model of Rough Surface Based on Fractal Theory. Chin. J. Mech. Eng. 2017, 30, 207–215. [Google Scholar] [CrossRef]
  10. Zhao, Y.; Yang, C.; Cai, L.; Shi, W.; Liu, Z. Surface contact stress-based nonlinear virtual material method for dynamic analysis of bolted joint of machine tool. Precis. Eng. 2016, 43, 230–240. [Google Scholar] [CrossRef]
  11. Huang, K.; Jin, J. Research on Bolt Preload Simulation Based on Virtual Material Method. Mach. Des. Manuf. 2012, 8, 148–150. [Google Scholar]
  12. Jana, T.; Anirban, M.; Prasanta, S. Dynamic contact interactions of fractal surfaces. Appl. Surf. Sci. 2017, 392, 872–882. [Google Scholar] [CrossRef]
  13. Zhang, K.; Li, G.; Gong, J.; Zhang, M. Normal contact stiffness of rough surfaces considering oblique asperity contact. Adv. Mech. Eng. 2019, 11, 1–14. [Google Scholar] [CrossRef]
  14. Pan, W.; Li, X.; Wang, L.; Guo, N.; Mu, J. A normal contact stiffness fractal prediction model of dry-friction rough surface and experimental verification. Eur. J. Mech. A-Solid. 2017, 66, 94–102. [Google Scholar] [CrossRef]
  15. Mayer, M.H.; Gaul, L. Segment-to-segment contact elements for modelling joint interfaces in finite element analysis. Mech. Syst. Signal Process. 2007, 21, 724–734. [Google Scholar] [CrossRef]
  16. Ye, H.; Huang, Y.; Li, P.; Li, Y.; Bai, L. Virtual material parameter acquisition based on the basic characteristics of the bolt joint interfaces. Tribol. Int. 2016, 95, 109–117. [Google Scholar] [CrossRef]
  17. Yoshimura, M. Computer-aided design improvement of machine tool structure incorporating joint dynamics data. In Annals of the CIRP; CIRP: Paris, French, 1979; Volume 28, pp. 241–246. [Google Scholar]
Figure 1. Dimensions of the specimen.
Figure 1. Dimensions of the specimen.
Materials 12 02381 g001
Figure 2. Schematic diagram of the three-dimensional model of VMM.
Figure 2. Schematic diagram of the three-dimensional model of VMM.
Materials 12 02381 g002
Figure 3. Schematic diagram of three-dimensional model of (a) theSDM and (b) the FEM.
Figure 3. Schematic diagram of three-dimensional model of (a) theSDM and (b) the FEM.
Materials 12 02381 g003
Figure 4. Grid diagram of the three-dimensional model of the FEM.
Figure 4. Grid diagram of the three-dimensional model of the FEM.
Materials 12 02381 g004
Figure 5. (a) Experimental fasten boards, and (b) experimental setup.
Figure 5. (a) Experimental fasten boards, and (b) experimental setup.
Materials 12 02381 g005
Figure 6. Curve of frequency response of specimen by experiment.
Figure 6. Curve of frequency response of specimen by experiment.
Materials 12 02381 g006
Figure 7. Comparison of the first five frequencies obtained by simulation and experiment.
Figure 7. Comparison of the first five frequencies obtained by simulation and experiment.
Materials 12 02381 g007
Table 1. The first five natural frequencies of the three equivalent methods of the joint interface.
Table 1. The first five natural frequencies of the three equivalent methods of the joint interface.
Methodsf1 (Hz)f2 (Hz)f3 (Hz)f4 (Hz)f5 (Hz)
VMM 71401049.81228.81707.11790.1
SDM855.81243.31540.62166.42237.8
FEM8671270.31504.22144.12211.7
Table 2. Comparison of the first five mode shapes obtained by the simulation and modal experiment.
Table 2. Comparison of the first five mode shapes obtained by the simulation and modal experiment.
OrderMEVMMFEMSDM
1 Materials 12 02381 i001 Materials 12 02381 i002 Materials 12 02381 i003 Materials 12 02381 i004
2 Materials 12 02381 i005 Materials 12 02381 i006 Materials 12 02381 i007 Materials 12 02381 i008
3 Materials 12 02381 i009 Materials 12 02381 i010 Materials 12 02381 i011 Materials 12 02381 i012
4 Materials 12 02381 i013 Materials 12 02381 i014 Materials 12 02381 i015 Materials 12 02381 i016
5 Materials 12 02381 i017 Materials 12 02381 i018 Materials 12 02381 i019 Materials 12 02381 i020
Table 3. Errors of the first five natural frequencies obtained by simulation and experiment.
Table 3. Errors of the first five natural frequencies obtained by simulation and experiment.
Order12345
Error of VMM8.12%7.32%7.26%9.14%9.93%
Error of FEM31.18%12.14%13.52%37.08%35.26%
Error of SDM29.49%9.75%16.27%38.51%36.86%

Share and Cite

MDPI and ACS Style

Han, R.; Li, G.; Gong, J.; Zhang, M.; Zhang, K. Experimental Verification and Comparative Analysis of Equivalent Methods on Metal’s Fixed Joint Interface. Materials 2019, 12, 2381. https://doi.org/10.3390/ma12152381

AMA Style

Han R, Li G, Gong J, Zhang M, Zhang K. Experimental Verification and Comparative Analysis of Equivalent Methods on Metal’s Fixed Joint Interface. Materials. 2019; 12(15):2381. https://doi.org/10.3390/ma12152381

Chicago/Turabian Style

Han, Renxiu, Guoxi Li, Jingzhong Gong, Meng Zhang, and Kai Zhang. 2019. "Experimental Verification and Comparative Analysis of Equivalent Methods on Metal’s Fixed Joint Interface" Materials 12, no. 15: 2381. https://doi.org/10.3390/ma12152381

APA Style

Han, R., Li, G., Gong, J., Zhang, M., & Zhang, K. (2019). Experimental Verification and Comparative Analysis of Equivalent Methods on Metal’s Fixed Joint Interface. Materials, 12(15), 2381. https://doi.org/10.3390/ma12152381

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