Next Article in Journal
Influence of Heat Treatment on the Corrosion Resistance of Aluminum-Copper Coating
Next Article in Special Issue
Research on Shape Control Characteristics of Non-oriented Silicon Steel for UCMW Cold Rolling Mill
Previous Article in Journal
Bidirectional Response of Weak-Axis End-plate Moment Connections: Numerical Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Concurrent Multiscale Simulations of Rough Lubricated Contact of Aluminum Single Crystal

1
State Key Laboratory of Tribology, Tsinghua University, Beijing 100084, China
2
School of Mechanical, Materials, Mechatronic and Biomedical Engineering, University of Wollongong, Wollongong, NSW 2522, Australia
3
School of Mechanical, Electronic and Control Engineering, Beijing Jiaotong University, Beijing 100044, China
*
Author to whom correspondence should be addressed.
Metals 2020, 10(7), 965; https://doi.org/10.3390/met10070965
Submission received: 4 June 2020 / Revised: 1 July 2020 / Accepted: 6 July 2020 / Published: 17 July 2020

Abstract

:
In this paper, a concurrent multiscale simulation strategy coupling atomistic and continuum models was proposed to investigate the three-dimensional contact responses of aluminum single crystal under both dry and lubricated conditions. The Hertz contact is performed by using both the multiscale and full molecular dynamics (MD) simulations for validation. From the contact area, kinetic energy and stress continuity aspects, the multiscale model shows good accuracy. It can also save at least five times the computational time compared with the full MD simulations for the same domain size. Furthermore, the results of lubricated contact show that the lubricant molecules could effectively cover the contact surfaces; thereby separating the aluminum surfaces and bearing the support loads. Moreover, the surface topography could be protected by the thin film formed by the lubricant molecules. It has been found that the contact area decreases obviously with increasing the magnitude of load under both dry and lubricated contacts. Besides, a decrease in contact area is also seen when the number of lubricant molecules increases. The present study has confirmed that the dimension of lubricated contacts could be greatly expanded during the simulation using the proposed multiscale method without sacrificing too much computational time and accuracy.

1. Introduction

Tribology is an interdisciplinary field of research that studies the adhesion, friction, lubrication and wear of contacts from the nano-scale to macro-scale, which is quite complicated, as the processes are usually governed by physical mechanisms at different length scales and can be affected by many factors, such as load, temperature, velocity, environment and so on [1,2,3,4,5,6,7]. For contact mechanics, in particular, not only the long-range elastic and plastic deformation play roles, but also the atomic scale roughness at the interface, reported by Luan and Robbins, has critical effects on the contact area and stress [4]. They revealed the breakdown of continuum models, and intrigued new investigations for a number of contacting problems at the nanoscale, which are still open for in-depth analyses. Over the last 40 years, the molecular dynamics (MD) simulations were widely used to explore the mechanisms of contact, friction, wear and lubrication at the nanoscale for various materials [8,9,10,11,12]. As a complement to conventional experiments, MD simulation can produce some detailed observations which cannot be obtained by experimental instruments. Despite the advancing development of computer science, MD still suffers from the temporal and spatial limitations. It is still very challenging to slow the timescale of MD down to realistic vales, although many efforts have been devoted to approaching such a target, such as accelerated molecular dynamics [13] and multiscale simulations [14,15,16]. The finite element method (FEM), based on the framework of continuum mechanics, has also been widely used to study the elastic and elastic–plastic contact between rough surfaces and deformations [17,18,19,20]. However, capturing the accurate information at the atomic scale is difficult by using FEM due to incorrect constitutive law. The development of the coupling methods that bridge the atomistic and continuum models is propelled to overwhelm the limitations of MD and FEM [14,15,16].
In general, the multiscale model could be divided into three regions: the critical region, the coarse region and the transition region [21,22,23]. The critical region, composed of atoms, is in charge of the local details of deformation, such as the contact area. The deformation of the coarse region, which is usually far from the critical region, could be described by the continuum models. The region, where the critical and coarse regions overlap, is called the transition region or the handshaking region. It is very critical to complete the information exchange between the two regions, and also to properly remove the unphysically spurious wave at this region [24,25]. Based on such philosophy, several multiscale methods have been established, such as the quasicontinuum method (QCM), bridging scale method (BSM), bridging domain method (BSD) and the hybrid simulation method (HSM). The QCM proposed by Tadmor et al. [21] is one of the pioneering works in multiscale coupling. It has a well-defined energy formulation for both critical and coarse regions. However, the element size in the critical region needs to be refined down to the atomic level. The BSM proposed by Wagner and Liu [22] reduced the atomistic calculations to a small domain with a description of the effects of the eliminated atoms via a time history kernel technique. The BDM proposed by Xiao and Belytschko [23] applied the Lagrange multipliers to impose displacement or velocity compatibility between the atomistic and continuum in a transition domain. For both the BSM and BDM, there is no need to refine the mesh to the atomic size. Particularly, as there is no mathematical framework, it is easier to implement the HSM proposed by Luan et al. [26] than other multiscale methods. Furthermore, the HSM yields the convincing results for the investigations of the contact and friction at the nanoscale benchmarked by pure MD simulations. Therefore, in this study, the HSM was used to study the 3D rough lubricated contact of the single crystal aluminum.
Recently, the HSM has been used for nanoindentation [27] and dry contact between rough surfaces [28,29]. Revealing the behaviors of confined thin lubricant layers between metal surfaces could help understand mixed lubrication. Using the HSM could provide more statistics data and thereby more accurate information compared to the MD simulations [30,31]. Furthermore, the MD region in the HSM can handle the complicated behaviors at the lubricating interfaces, while the coarse region can properly accommodate the elastic deformation from the contact zone [26]. In the previous MD simulations, however, the very limited substrate height with the rigid boundary increased the system rigidity, which could result in the inappropriate accumulation of kinetic energy at the contact zone [32]. To the best of our knowledge, multiscale methods have not been used to perform the simulations of mixed lubrication. In this paper, we used the HSM to study the behaviors of the hexadecane molecules confined within the aluminum tribo-surfaces, which have attracted attention for many years due to their light-weight nature, good formability and low cost [33]. Moreover, Al single crystal as a model material was chosen for comparisons. The changes in the root mean square (RMS) roughness, contact area and surface pressure under dry and lubricated contacts are presented.

2. Methodology

2.1. Model Set-Up

The multiscale (MU) model for the contact between a rigid flat plane and a deformable substrate with a self-affine surface is shown in Figure 1. The substrate was built from Face-center-cubic (FCC) crystal aluminum with a (001) top surface. The lattice constant a of the crystal at 300 K is 0.4045 nm. The substrate dimension is 80a × 80a × 80a, which approximately corresponds to 32.36 nm × 32.36 nm × 32.36 nm. This size is at least 8 times larger than the previous MD simulations [30,34]. Except for the lubricant molecules, the total number of atoms ranges from 670,000 to 700,000 for different surface roughness (0.2 nm, 0.5 nm and 0.8 nm). The specific RMS values were chosen and the rough surfaces were created by using the random midpoint displacement (RMD) algorithm, following the previous MD simulations [34,35]. The corresponding RMS roughness of the generated atomic surface was 0.204, 0.492 and 0.771 nm, respectively. The atomic region of the substrate consists of three layers: (2) the deformable layer, (4) the thermostat layer and (3) the pad layer. The contact occurs among the rigid plane, lubricant and deformable layer. The thermostatting technique was applied to the atoms in the thermostat layer to maintain the constant temperature of 300 K for the system. Along the x and y directions, a periodic boundary condition was applied to both the atomic and FEM regions. The bottom boundary of the FEM region was fixed while the top surface of the deformable layer was free. The FEM region consists of 33,600 elements and 6615 nodes, corresponding to 1,433,600 atoms. A constitutive anisotropic law describes the relation between stress and strain for each integration point. The elastic constants C11 = 114.0, C12 = 61.6, C44 = 31.6 (unit: GPa) were taken from Mishin’s work [36].
The key ideas of the coupling strategy for the HSM is described here. The pad layer as a part of the FEM region provides a physical boundary for the atomic region. The displacements of each atom (e.g., uβ) in the pad layer are interpolated from the nodal displacements [26] as shown in Equation (1):
u β = j = 1 m nodes N j d j ,
where N is the shape function of atom β corresponding to the nodes in each element. mnodes is the number of the nodes in each element. As can be seen in Figure 1, the top boundary of the FEM region lies in the thermostat layer. The nodal displacements on this boundary (e.g., di) were obtained by averaging the displacements of atoms within the average spheres, as shown in Equation (2). The radius of each sphere is 4a, corresponding to the smallest size of the element edge. The four-node tetrahedron elements were used to divide the continuum region:
d i = α = 1 n atoms ω α u α
where natoms is the number of atoms in the average circle for node i. The weight ω of each atom within the averaging circle with a radius of rav can be calculated by a weighting function [37]:
ω ( d ) = { 1 3 d 2 + 2 d 3 , d < 1 0 , d 1
where d = r/rav and r is the distance between the atom within the average circle and the node.
The timestep is 2 fs. Before loading, the substrate is allowed to relax for equilibrium. The system has 0.4 ns to equilibrate after applying a load. To study the effects of different pressure on the lubricated contact, 0.05, 0.15 and 0.25 GPa were applied on the rigid plane, respectively. During the loading, at an interval of 200 timesteps, the forces and stresses were recorded. The calculation time of the multiscale simulation ranges from 74 to 98 h due to the different number of atoms. For the different cases here, the number of FEM nodes are not changed.

2.2. Force Field

As the plastic deformation on the rough surface with the atoms reshuffling or dislocations usually involves during compression and shearing between the rough surfaces, the EAM potential [38] was used to model the interactions between the aluminum atoms in the substrate. The atomic interactions between the rigid plane and the substrate were calculated by the Lennard–Jones (LJ) potential, which is characterized by the distance and energy parameters σ and ε, respectively. The LJ potential was used to simulate the aluminum oxide surface by reducing ε [34]. A united-atom (UA) model [39] was used to model the hexadecanes molecules. The UA model simplified all the CHX groups with pseudo carbon atoms. Linear chain molecules of hexadecane with the chemical formula C16H34 were chosen as the lubricant liquid. A chain of hexadecane is made of connected CH3 and CH2 groups as shown in Figure 2a. The optimization of the conformation of the hexadecane molecule in Figure 2b is attributed to the intra-molecular interactions, which include bond stretching, angle bending and the dihedral angle torsion given in Equations (1)–(3). To study the effects of the amount of lubricant on the lubricated contact, the different number of C16H34 molecules was added (0–4000) between the aluminum tribo-surfaces:
E bond = k b ( r     r 0 ) 2 ,
E bend = k θ ( θ     θ 0 ) 2 ,
E torsion = c 0   +   c 1 [ 1   +   cos ϕ ]   +   c 2 [ 1     cos ( 2 ϕ ) ]   +   c 3 [ 1   +   cos ( 3 ϕ ) ] ,
The corresponding parameters in Equations (4)–(6) are taken from Reference [39]. The inter-molecular and liquid–solid interactions are described by the LJ potential with a cutoff distance of 1.4 nm. The LJ parameters for the liquid–solid interactions were estimated via the Lorentz–Berthelot mixing rules (Table 1).

2.3. Contact Area

Calculations for the real contact area are of importance at the atomic scale, especially when the dimension of the devices is reduced down to the nanometer. The real contact area has an important effect on heat conduction, friction and adhesion. In continuum mechanics, the contact area can be obtained based on the edge of the contact zone. It is difficult to accurately define the atomic area and even the atomic contact. The atomic contact is determined by either repulsive force [28,40] or distance [35,41]. Here, the contact is established when the distance of two atoms is less than or equals to 0.5 nm [30], and the real contact area is estimated from a projected rectangular cuboid method [34]. In this method, the n×n grids as a bottom face of the rectangular cuboids are plotted on the x–y plane, on which all the atoms of the deformable layer are projected. Those atoms are divided into the corresponding grids according to their x and y coordinates. The length of each cuboid is the maximum z coordinate values of the atoms within it. Similar cuboids can be plotted for the rigid plane. The top faces of the deformable layer cuboids and the bottom faces of the rigid plane cuboids are used to determine the atomic contact.

2.4. Multiscale Strategy Valication

Here, we validate our MU model by performing a 3D Hertz contact, benchmarked with full MD simulations, following the the Anciaux and Molinali’s work [32]. Details of the model setup for the Hertz contact is shown in Figure S1 and demonstrated in the supplementary materials. The contact area, kinetic energy and stress continuality were considered. Figure S2a shows the relationship between the contact area and the normalized indent displacement. The contact region consists of the atoms on the stepped tip and is close to the area of a circle. The contact area was determined by averaging the distances between the center of the circle and the outermost atoms. In terms of the contact area, the MU predicts the MD well. It is found from Figure S2a that the results between the two models are very close. From the energy point of view, the MU works well based on the results in Figure S2b. Although there are a few deviating points, the kinetic energies of MU and MD are close to each other. In the MU model, the atomic, or the virial stress was calculated for the MD region while the Cauchy stress for the FEM region. The stress component σzz of both the MU model and the MD model when normalized indent displacement is 1.31a was shown in Figure S2c. First and foremost, the HSM can ensure stress continuity after comparing σzz of both the FEM region and the MD region (dotted square). Another important result shown in Figure S2c is that the MU model can predict the full MD model well.
From the contact area, kinetic energy and stress continuity aspects, our MU simulations show good accuracy compared with the full MD simulations. More importantly, it can save at least 5 times the computational time compared with the full MD simulations for the same domain size. With relatively good computational cost and accuracy, the MU model was here extended to study lubricated contacts.

3. Results and Discussion

3.1. RMS Roughness

Before compression, the rigid plane was above the substrate with a distance larger than the cut-off radius of the interaction potential. As mentioned above, the original RMS roughness of the generated atomic surface was 0.204, 0.492 and 0.771 nm, respectively. The atoms at the surface lack the full neighbors, and are attracted by the atoms at the subsurface. Therefore, the relaxing process at a temperature of 300 K modified the surface roughness with RMS = 0.198, 0.484 and 0.760 nm, respectively. Although the RMS roughness of the surfaces changes, it is still close to the required values. After relaxation, the top plane was forced to move downwards under constant pressure. Without any resisting force, the top plane will impact on the substrate with a high speed. This is not acceptable, since it will cause significant plastic deformation on the substrate. Therefore, the top plane was forced to move at a very low speed until it touched the deformable layer. The system heights were calculated during the whole compression. In Figure 3, the relationship between the system height and simulation time for dry and lubricated contact is shown. It indicates that all the simulations reached equilibrium after 0.4 ns. The heights in the lubricated cases were higher than the dry counterparts, as shown in Figure 3a,b. The thickness of the thin film formed by the hexadecane molecules contributed to the system height. In Figure 3c, two thousand hexadecane molecules had no effect on the system height when the surface roughness was 0.8 nm (indicated by the dashed circle in the Figure), since there was not enough quantity to fill the cavity between the two contacting surfaces. Meanwhile, four thousand hexadecane molecules lifted the top plane by a distance of 0.4 nm as indicated by the dashed arrow in Figure 3c. The system height is also related to the chain length of the lubricant molecules [42]. The effective separation for the contacting surfaces is dependent on the added amount of lubricant and surface roughness. A few lubricant molecules can work well for relatively flat surfaces, while much more molecules are required to fill a cavity in rougher surfaces and effectively separate them. The volume of the hexadecane molecule is also related to the rotation of C–C bonds—thereby affecting the system heights in Figure 3.
When the contact occurs, the top plane will flatten the substrate asperities as shown in Figure 4b,c. The RMS roughness of flattened surface changes as the load increases. The RMS of the flattened surface is denoted as RMSf, while that of the initial surface RMSu. The change of RMS, therefore, is defined as ∆RMS = RMSf − RMSu. Firstly, it indicates that RMSf decreases with the increasing loads for three surfaces, and this relationship is independent of RMSu. However, the percentage of change in RMS roughness is related to RMSu. There was a much smaller change for the rougher surfaces (solid triangle and square in Figure 4a), because the contact for those rougher surfaces occurred on a few high asperities which support the loads. It is the rest part on the surface that dominated the roughness. However, the surface RMS changed more for RMS = 0.2 nm, because much more asperities on the substrate contacted with the top plane were simultaneously flattened. The peak-to-valley (PV) distances [34] indicate roughly the number of the contacting asperities. The PV distance here is determined by the difference between the maximum peak height to the maximum valley depth of the atomistic surface. The PV of the atomic surfaces for RMS = 0.2, 0.5 and 0.8 nm are 0.172, 0.3285 and 0.514 nm, respectively. The height of the highest asperity for RMS = 0.2 nm is approximately 0.86 nm. Therefore, many asperities support the top plane simultaneously. Figure 5 shows the effects of the molecule number on the RMS change. Zheng et al. [42] demonstrated that the appropriate amount of lubricant is dependent on the RMS roughness. The increase in the molecule number reduces the change in RMS roughness or asperity flattening as shown in Figure 5. After compression, the lubricant molecules filled the valleys and covered the surface as shown in the insets. With 4000 molecules, the surface could be entirely covered. The asperities were covered by the lubricant, and were not contacted by the top plane. Consequently, the RMS changed only by about 1.8%. In addition, Wu et al. [43] indicated that the lubricant viscosity has an important effect on the changes in the surface roughness. In their work, the smaller reduction in the surface roughness was obtained for the lubricant with a higher viscosity. The viscosity in Wu’s work plays a similar role as the lubricant amount here. A well-formed lubricant film reduced the direct metal–metal contact, and protected the surfaces. For the lubricated contact, the new surface was formed by the lubricant molecules and substrate atoms. The RMS roughness of the new surface was 0.33 and 0.06 nm for 3000 and 4000 molecules, respectively. They are both much lower than the original roughness of 0.8 nm.

3.2. Contact Area

The effects of roughness and lubricant on the contact area calculated by the project method will be studied in this section. Figure 4a indicates that the contact area goes up with the increasing loads for the surfaces with three different RMS. Under the load of 0.05 GPa, the contact area of RMS = 0.2 nm is close to 50%, which is much higher than 4% and 3.5% of the RMS = 0.5 and 0.8 nm. As mentioned in the previous section, the heights of the asperities on the surface of RMS = 0.2 nm are relatively low, so a number of them were in contact simultaneously. For the rougher surface, however, only few asperities supported the load applied from the top plane. With the increasing loads, the contacting asperity deformed and their heights decreased while the lower non-contacting asperities were joined to support the top plane. Therefore, the contact area of the three surfaces was increased in Figure 4a. Considering the different topography due to the different random seed when generating the self-affine surfaces [35] and adhesive strength [34], the results could deviate a little. The relationship between the contact area and the loads obtained in this work qualitatively agrees with the published literature [44].
For the lubricated contacts, the contact area still increased with the loads at the range of 0.05 GPa (red diamond) to 0.25 GPa (purple square) as shown in Figure 6a. However, the contact area reduced with the increasing number of molecules. The lubricants could separate the substrate from the top plane. The separating gap depends on the lubricant amount. It indicates that under 0.25 GPa the contact area of 0 and 1000 molecules is almost identical while the contact area of 4000 molecules was reduced by 75%. However, it was not enough for the 2000 molecules to fill the valleys on the rough surface. From the point of view of normal force, it could also provide some clues. The normal force on the top plane of lubricated contact was divided into two parts: normal force FS from the substrate and normal force FL from the lubricant. The sum of FS and FL is equal to the normal force F0 on the top plane of dry contact (no lubricant molecules). Figure 6b shows the normal force on the top plane supported by dry and lubricated contacts as a function of molecule number. In the figure, the support force FS/F0 from the asperities roughly decreased with the larger number of molecules, which is consistent with the trend in the contact area. At the same time, the normal force FL/F0 from the lubricant increased with the higher number of molecules. The decrease in the normal force FS/F0 under each load indicates that the number of the contacting asperities dropped. At 0.25 GPa, more lubricant endured the larger force. When 4000 hexadecane molecules were used, FL reached 87.17 nN under 0.25 GPa, which corresponds to 85.45% of FS in dry contact. The project method can be used to calculate the volume occupied by the hexadecane molecules—thereby estimating its density. For RMS = 0.8 nm with 4000 molecules under 0.15 and 0.25 GPa, the hexadecane has the density of 6.9 and 7.2 g/cm3, respectively. The estimations are very close to 7.7 g/cm3 obtained from the previous experiments [45].

3.3. Surface Pressure

Surface pressure was calculated based on the projection method in Section 2.3. There are other approaches for calculating surface pressure, but for the comparisons the projection method was used. The virial stress of each atom was calculated at every 20 ps. After the cuboids were obtained, the surface pressure was estimated by averaging the virial stress along the z direction of all atoms within each cuboid. Before that, we show in Figure 7 the deformation of the FEM region due to the applied loads from the top plane on the substrate. Before compression, there is almost no deformation at the coarse region as shown in Figure 7a. With the compression of the top plane, the deformation is small in this region as shown in Figure 7b, as the FEM region is far from the surface. However, more deformation occurred in the FEM region, since the asperities were flattened in the corresponding area on the surface, indicating that the coarse region can properly accommodate the elastic deformation from the contact zone.
Figure 8 shows the pressure distribution for the unlubricated and lubricated contacts. To obtain the distribution histogram, the surface pressure in each cuboid was binned with an interval of 0.05 GPa. The shape of the probability distribution in this work was in agreement with the previous works [30,34]. It seems that the system size and the roughness of the top surface have little effect on the shape of the probability distribution. In Figure 8a, the sum of the negative pressure distribution is 88.1%, 80.0% and 66.79% for 0.05, 0.15 and 0.25 GPa, respectively. Meanwhile, in Figure 8b, the obtained values are 89.2%, 72.8% and 63.4%, respectively. More than half of the surface area mainly suffered the attractive forces. In addition, for both dry and lubricated surfaces, the probability sum of the positive pressure goes up as the loads increase as shown in Figure 8. Further compression caused the increase in the repulsive forces between the top plane, lubricants and substrate. As demonstrated by Cheng et al. [46], the lubricant molecules spread the pressure over a larger contact area. The shifting of probability with the increasing load for the lubricated contacts is much more than that for the dry contacts, indicating the significant role of lubricant molecules in bearing the loads.

4. Conclusions

In the present study, a concurrent multiscale model that combines the MD and FEM methods was proposed and successfully applied to study the three-dimensional rough lubricated contact of aluminum single crystal. We performed the Hertz contact to compare the MU model with the full MD model. Our MU simulations show good accuracy from the contact area, kinetic energy and stress continuity aspects, and it can save at least five times the computational time compared with the full MD simulations for the same domain size. With relatively good computational cost and accuracy, the MU model here has been extended to study lubricated contacts. A systematic investigation considering the effects of the lubricant amount, surface roughness and loads during contact in mixed lubrication at the atomic scale was conducted. It was found that the sufficient lubricant molecules could fill the cavity and cover the large portion of contacting surfaces and effectively separate the tribo-pair. Meanwhile, the insufficient lubricant molecules can work well for relatively flat surfaces. In addition, the results show that the increasing number of lubricant molecules leads to smaller changes in surface RMS under the studied loads, varying from 0.05 to 0.25 GPa, indicating that the thin film formed by the lubricant molecules is beneficial to maintain the initial surface topography. For both the lubricated and unlubricated contacts, the contact area decreases with the increasing loads. As the tribo-pair gap increases with the increasing lubricant amount, the contact area decreases, and the larger portion of loads is supported by the lubricant molecules. The changes in pressure distribution for the different loads in the lubricated contacts confirm the contributions of the lubricant molecules to the load-bearing.

Supplementary Materials

The following are available online at https://www.mdpi.com/2075-4701/10/7/965/s1, Figure S1: Geometry of the multiscale (a) and full MD (b) simulations, Figure S2: Comparisons between MU and full MD simulations: (a) Normalized contact area as a function of normalized indent displacement. (b) Kinetic energy near the substrate surface as a function normalized indent displacement. (c) Changes of stress component σzz in a line along the z direction when normalized indent displacement is 1.31a. Indent displacement and Z are normalized by lattice constant of a.

Author Contributions

Conceptualization, software, visualization, writing—original draft, J.Z.; methodology, J.Z. and L.S.; validation, writing—review and editing, L.S. and Z.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Fundamental Research Funds for the Central Universities (grant number 2020JBM059), the National Natural Science Foundation of China (grant number 51905296) and Australian Research Council (ARC) Discovery Early Career Researcher Award (DECRA) project (DE180100124).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Persson, B.N.J.; Albohr, O.; Tartaglino, U.; Volokitin, A.I.; Tosatti, E. On the nature of surface roughness with application to contact mechanics, sealing, rubber friction and adhesion. J. Phys. Condens. Matter 2004, 17, R1–R62. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Bhushan, B. Contact mechanics of rough surfaces in tribology: Multiple asperity contact. Tribol. Lett. 1998, 4, 1–35. [Google Scholar] [CrossRef]
  3. Manini, N.; Mistura, G.; Paolicelli, G.; Tosatti, E.; Vanossi, A. Current trends in the physics of nanoscale friction. Adv. Phys. X 2017, 2, 569–590. [Google Scholar] [CrossRef] [Green Version]
  4. Luan, B.; Robbins, M.O. The breakdown of continuum models for mechanical contacts. Nature 2005, 435, 929–932. [Google Scholar] [CrossRef]
  5. Deng, G.; Tieu, A.K.; Lan, X.; Su, L.; Wang, L.; Zhu, Q.; Zhu, H. Effects of normal load and velocity on the dry sliding tribological behaviour of CoCrFeNiMo0.2 high entropy alloy. Tribol. Int. 2020, 144, 106116. [Google Scholar] [CrossRef]
  6. Ta, T.D.; Tieu, A.K.; Zhu, H.; Zhu, Q.; Kosasih, P.B.; Zhang, J.; Deng, G. Tribological Behavior of Aqueous Copolymer Lubricant in Mixed Lubrication Regime. ACS Appl. Mater. Interfaces 2016, 8, 5641–5652. [Google Scholar] [CrossRef]
  7. Deng, G.Y.; Tieu, A.K.; Su, L.H.; Zhu, H.T.; Reid, M.; Zhu, Q.; Kong, C. Microstructural study and residual stress measurement of a hot rolling work roll material during isothermal oxidation. Int. J. Adv. Manuf. Technol. 2019, 102, 2107–2118. [Google Scholar] [CrossRef]
  8. Zhang, J.; Xu, Q.; Gao, L.; Ma, T.; Qiu, M.; Hu, Y.; Wang, H.; Luo, J. A molecular dynamics study of lubricating mechanism of graphene nanoflakes embedded in Cu-based nanocomposite. Appl. Surf. Sci. 2020, 511, 145620. [Google Scholar] [CrossRef]
  9. Ma, T.-B.; Hu, Y.-Z.; Wang, H. Molecular dynamics simulation of shear-induced graphitization of amorphous carbon films. Carbon 2009, 47, 1953–1957. [Google Scholar] [CrossRef]
  10. Kong, N.; Wei, B.; Zhuang, Y.; Zhang, J.; Li, H.; Wang, B. Effect of Compressive Prestrain on the Anti-Pressure and Anti-Wear Performance of Monolayer MoS2: A Molecular Dynamics Study. Nanomaterials 2020, 10, 275. [Google Scholar] [CrossRef] [Green Version]
  11. Zhao, X.; Lu, C.; Tieu, A.K.; Zhan, L.; Huang, M.; Su, L.; Pei, L.; Zhang, L. Deformation twinning and dislocation processes in nanotwinned copper by molecular dynamics simulations. Comp. Mater. Sci. 2018, 142, 59–71. [Google Scholar] [CrossRef] [Green Version]
  12. Zhang, L.; Shibuta, Y.; Lu, C.; Huang, X. Interaction between nano-voids and migrating grain boundary by molecular dynamics simulation. Acta Mater. 2019, 173, 206–224. [Google Scholar] [CrossRef]
  13. Li, Q.; Dong, Y.; Perez, D.; Martini, A.; Carpick, R.W. Speed Dependence of Atomic Stick-Slip Friction in Optimally Matched Experiments and Molecular Dynamics Simulations. Phys. Rev. Lett. 2011, 106, 126101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Zhang, J.; Michal, G.; Tieu, A.K.; Zhu, H.T.; Deng, G.Y. Hertz Contact at the Nanoscale with a 3D Multiscale Model. Appl. Mech. Mater. 2016, 846, 306–311. [Google Scholar] [CrossRef]
  15. Liu, W.K.; Karpov, E.G.; Zhang, S.; Park, H.S. An introduction to computational nanomechanics and materials. Comput. Method Appl. Mech. Eng. 2004, 193, 1529–1578. [Google Scholar] [CrossRef]
  16. Michal, G.; Lu, C.; Kiet, T.A. Multiscale model of elastic nanocontacts. Comp. Mater. Sci. 2014, 81, 98–103. [Google Scholar] [CrossRef]
  17. Hyun, S.; Pei, L.; Molinari, J.F.; Robbins, M.O. Finite-Element analysis of contact between elastic self-affine surfaces. Phys. Rev. E 2004, 70, 026117. [Google Scholar] [CrossRef] [Green Version]
  18. Pei, L.; Hyun, S.; Molinari, J.F.; Robbins, M.O. Finite element modeling of elasto-Plastic contact between rough surfaces. J. Mech. Phys. Solids 2005, 53, 2385–2409. [Google Scholar] [CrossRef]
  19. Deng, G.Y.; Tieu, A.K.; Su, L.H.; Zhu, H.T.; Zhu, Q.; Zamri, W.F.H.; Kong, C. Characterizing deformation behaviour of an oxidized high speed steel: Effects of nanoindentation depth, friction and oxide scale porosity. Int. J. Mech. Sci. 2019, 155, 267–285. [Google Scholar] [CrossRef]
  20. Deng, G.Y.; Zhu, H.T.; Tieu, A.K.; Su, L.H.; Reid, M.; Zhang, L.; Wei, P.T.; Zhao, X.; Wang, H.; Zhang, J.; et al. Theoretical and experimental investigation of thermal and oxidation behaviours of a high speed steel work roll during hot rolling. Int. J. Mech. Sci. 2017, 131–132, 811–826. [Google Scholar] [CrossRef] [Green Version]
  21. Tadmor, E.B.; Ortiz, M.; Phillips, R. Quasicontinuum analysis of defects in solids. Philos. Mag. A 1996, 73, 1529–1563. [Google Scholar] [CrossRef]
  22. Wagner, G.J.; Liu, W.K. Coupling of atomistic and continuum simulations using a bridging scale decomposition. J. Comput. Phys. 2003, 190, 249–274. [Google Scholar] [CrossRef]
  23. Xiao, S.P.; Belytschko, T. A bridging domain method for coupling continua with molecular dynamics. Comput. Method Appl. Mech. Eng. 2004, 193, 1645–1669. [Google Scholar] [CrossRef]
  24. Park, H.S.; Liu, W.K. An introduction and tutorial on multiple-scale analysis in solids. Comput. Method Appl. Mech. Eng. 2004, 193, 1733–1772. [Google Scholar] [CrossRef]
  25. Zhang, J.; Tieu, K.; Michal, G.; Zhu, H.; Zhang, L.; Su, L.; Deng, G.; Wang, H. A damping boundary condition for atomistic-continuum coupling. Chin. Phys. B 2017, 26, 68702. [Google Scholar] [CrossRef]
  26. Luan, B.Q.; Hyun, S.; Molinari, J.F.; Bernstein, N.; Robbins, M.O. Multiscale modeling of two-Dimensional contacts. Phys. Rev. E 2006, 74, 046710. [Google Scholar] [CrossRef]
  27. Zhu, P.; Hu, Y.; Fang, F.; Wang, H. Multiscale simulations of nanoindentation and nanoscratch of single crystal copper. Appl. Surf. Sci. 2012, 258, 4624–4631. [Google Scholar] [CrossRef]
  28. Luan, B.; Robbins, M.O. Hybrid Atomistic/Continuum Study of Contact and Friction between Rough Solids. Tribol. Lett. 2009, 36, 1–16. [Google Scholar] [CrossRef]
  29. Zhang, J.; Zhang, L.; Tieu, A.K.; Michal, G.; Zhu, H.T.; Deng, G.Y. Finite-Temperature Multiscale Simulations for 3D Nanoscale Contacts. Appl. Mech. Mater. 2016, 846, 288–293. [Google Scholar] [CrossRef]
  30. Zheng, X.; Zhu, H.; Kiet, T.A.; Kosasih, B. A molecular dynamics simulation of 3D rough lubricated contact. Tribol. Int. 2013, 67, 217–221. [Google Scholar] [CrossRef]
  31. Zheng, X.; Zhu, H.; Tieu, A.K.; Kosasih, B. Roughness and Lubricant Effect on 3D Atomic Asperity Contact. Tribol. Lett. 2013, 53, 215–223. [Google Scholar] [CrossRef]
  32. Anciaux, G.; Molinari, J.-F. Contact mechanics at the nanoscale, a 3D multiscale approach. Int. J. Numer. Methods Eng. 2009, 79, 1041–1067. [Google Scholar] [CrossRef]
  33. Su, L.; Deng, G.; Luzin, V.; Wang, H.; Wang, Z.; Yu, H.; Li, H.; Tieu, A.K. Effect of cryogenic temperature equal channel angular pressing on microstructure, bulk texture and tensile properties of AA1050. Mater. Sci. Eng. A 2020, 780, 139190. [Google Scholar] [CrossRef]
  34. Spijker, P.; Anciaux, G.; Molinari, J.-F. The effect of loading on surface roughness at the atomistic level. Comput. Mech. 2012, 50, 273–283. [Google Scholar] [CrossRef]
  35. Spijker, P.; Anciaux, G.; Molinari, J.-F. Relations between roughness, temperature and dry sliding friction at the atomic scale. Tribol. Int. 2013, 59, 222–229. [Google Scholar] [CrossRef]
  36. Mishin, Y.; Farkas, D.; Mehl, M.J.; Papaconstantopoulos, D.A. Interatomic potentials for monoatomic metals from experimental data and ab initio calculations. Phys. Rev. B 1999, 59, 3393–3407. [Google Scholar] [CrossRef] [Green Version]
  37. Miller, R.E.; Tadmor, E.B. A unified framework and performance benchmark of fourteen multiscale atomistic/continuum coupling methods. Model. Simul. Mater. Sci. Eng. 2009, 17, 053001. [Google Scholar] [CrossRef]
  38. Foiles, S.M.; Baskes, M.I.; Daw, M.S. Embedded-atom-method functions for the fcc metals Cu, Ag, Au, Ni, Pd, Pt, and their alloys. Phys. Rev. B 1986, 33, 7983–7991. [Google Scholar] [CrossRef]
  39. Martin, M.G.; Siepmann, J.I. Transferable Potentials for Phase Equilibria. 1. United-Atom Description of n-Alkanes. J. Phys. Chem. B 1998, 102, 2569–2577. [Google Scholar] [CrossRef]
  40. Mo, Y.; Turner, K.T.; Szlufarska, I. Friction laws at the nanoscale. Nature 2009, 457, 1116–1119. [Google Scholar] [CrossRef]
  41. Anciaux, G.; Ramisetti, S.B.; Molinari, J.F. A finite temperature bridging domain method for MD-FE coupling and application to a contact problem. Comput. Method Appl. Mech. Eng. 2012, 205–208, 204–212. [Google Scholar] [CrossRef] [Green Version]
  42. Zheng, X.; Zhu, H.; Kosasih, B.; Tieu, A.K. A molecular dynamics simulation of boundary lubrication: The effect of n-alkanes chain length and normal load. Wear 2013, 301, 62–69. [Google Scholar] [CrossRef]
  43. Wu, C.; Zhang, L.; Li, S.; Jiang, Z.; Qu, P. A unified method for characterizing multiple lubrication regimes involving plastic deformation of surface asperities. Tribol. Int. 2016, 100, 70–83. [Google Scholar] [CrossRef]
  44. Campañá, C.; Müser, M.H. Contact mechanics of real vs. randomly rough surfaces: A Green’s function molecular dynamics study. Europhys. Lett. 2007, 77, 38005. [Google Scholar] [CrossRef]
  45. Tanaka, H.; Yamaki, Y.; Kato, M. Solubility of carbon dioxide in pentadecane, hexadecane, and pentadecane + hexadecane. J. Chem. Eng. Data 1993, 38, 386–388. [Google Scholar] [CrossRef]
  46. Cheng, S.; Luan, B.; Robbins, M.O. Contact and friction of nanoasperities: Effects of adsorbed monolayers. Phys. Rev. E 2010, 81, 016102. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Schematic of the model set-up. a = 0.4045 nm is the lattice constant of the single crystal aluminum.
Figure 1. Schematic of the model set-up. a = 0.4045 nm is the lattice constant of the single crystal aluminum.
Metals 10 00965 g001
Figure 2. Schematic of one hexadecane molecule’s (a) initial structure and (b) after the optimization of conformation. The pink and yellow atoms indicate the united atoms of CH2 and CH3.
Figure 2. Schematic of one hexadecane molecule’s (a) initial structure and (b) after the optimization of conformation. The pink and yellow atoms indicate the united atoms of CH2 and CH3.
Metals 10 00965 g002
Figure 3. Variation of the system heights with the time for the root mean square (RMS) = (a) 0.2 nm, (b) 0.5 nm, (c) 0.8 nm under 0.25 MPa. The red lines indicate dry contact while the green and blue lines indicate lubricated contact with 2000 and 4000 hexadecane molecules, respectively.
Figure 3. Variation of the system heights with the time for the root mean square (RMS) = (a) 0.2 nm, (b) 0.5 nm, (c) 0.8 nm under 0.25 MPa. The red lines indicate dry contact while the green and blue lines indicate lubricated contact with 2000 and 4000 hexadecane molecules, respectively.
Metals 10 00965 g003
Figure 4. (a) RMS change ∆RMS and contact area A/A0 of dry contact as a function of pressure. The circle, square and triangle indicate RMS = 0.2, 0.5 and 0.8 nm, respectively. Changes in the topography before (b) and after (c) loading for the surfaces with RMS = 0.2, 0.5 and 0.8 nm, respectively.
Figure 4. (a) RMS change ∆RMS and contact area A/A0 of dry contact as a function of pressure. The circle, square and triangle indicate RMS = 0.2, 0.5 and 0.8 nm, respectively. Changes in the topography before (b) and after (c) loading for the surfaces with RMS = 0.2, 0.5 and 0.8 nm, respectively.
Metals 10 00965 g004
Figure 5. Variations of the RMS change ∆RMS with the molecule number under 0.25 GPa. The RMS roughness of the surface is 0.8 nm. Insets show the surface topography and the hexadecane molecule distribution after loading.
Figure 5. Variations of the RMS change ∆RMS with the molecule number under 0.25 GPa. The RMS roughness of the surface is 0.8 nm. Insets show the surface topography and the hexadecane molecule distribution after loading.
Metals 10 00965 g005
Figure 6. Contact area (a) and the normalized normal force (b) as a function of the lubricant molecule number for surface with a roughness RMS = 0.8 nm. The molecule number of zero indicates dry contact. The load is 0.25 GPa in (b).
Figure 6. Contact area (a) and the normalized normal force (b) as a function of the lubricant molecule number for surface with a roughness RMS = 0.8 nm. The molecule number of zero indicates dry contact. The load is 0.25 GPa in (b).
Metals 10 00965 g006
Figure 7. Deformation of the FEM region (a) before compression and (b) after compression for dry contact.
Figure 7. Deformation of the FEM region (a) before compression and (b) after compression for dry contact.
Metals 10 00965 g007
Figure 8. Pressure distribution in the dry contact (a) and 4000 molecules lubricated contact (b) under 0.05 GPa (red line), 0.15 GPa (green line) and 0.25 GPa (blue line).
Figure 8. Pressure distribution in the dry contact (a) and 4000 molecules lubricated contact (b) under 0.05 GPa (red line), 0.15 GPa (green line) and 0.25 GPa (blue line).
Metals 10 00965 g008
Table 1. Parameters for the Lennard–Jones (LJ) potentials for the non-bonded interactions.
Table 1. Parameters for the Lennard–Jones (LJ) potentials for the non-bonded interactions.
Neighborsσ (nm)ε (eV)
Al Al a0.25960.010677
Al CH 2 c0.32730.006506
AlCH 3 c0.31730.009496
CH 2 CH 2 b0.39500.003964
CH 3 CH 3 b0.37500.008445
CH 3 CH 2 c0.38500.005786
a taken from Ref. [34]. b taken from Ref. [39]. c calculated from the Lorentz–Berthelot rules.

Share and Cite

MDPI and ACS Style

Zhang, J.; Su, L.; Wang, Z. Concurrent Multiscale Simulations of Rough Lubricated Contact of Aluminum Single Crystal. Metals 2020, 10, 965. https://doi.org/10.3390/met10070965

AMA Style

Zhang J, Su L, Wang Z. Concurrent Multiscale Simulations of Rough Lubricated Contact of Aluminum Single Crystal. Metals. 2020; 10(7):965. https://doi.org/10.3390/met10070965

Chicago/Turabian Style

Zhang, Jie, Lihong Su, and Zhongnan Wang. 2020. "Concurrent Multiscale Simulations of Rough Lubricated Contact of Aluminum Single Crystal" Metals 10, no. 7: 965. https://doi.org/10.3390/met10070965

APA Style

Zhang, J., Su, L., & Wang, Z. (2020). Concurrent Multiscale Simulations of Rough Lubricated Contact of Aluminum Single Crystal. Metals, 10(7), 965. https://doi.org/10.3390/met10070965

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