scholarly journals Multiscale Friction Simulation of Dry Polymer Contacts: Reaching Experimental Length Scales by Coupling Molecular Dynamics and Contact Mechanics

Author(s):  
Daniele Savio ◽  
Jannik Hamann ◽  
Pedro A. Romero ◽  
Christoph Klingshirn ◽  
Ravindrakumar Bactavatchalou ◽  
...  

Abstract This work elucidates friction in Poly-Ether-Ether-Ketone (PEEK) sliding contacts through multiscale simulations. At the nanoscale, non-reactive classical molecular dynamics (MD) simulations of dry and water-lubricated amorphous PEEK-PEEK interfaces are performed. During a short running-in phase, we observe structural transformations at the sliding interface that result in flattening of the initial nanotopographies accompanied by strong polymer chain alignment in the shearing direction. Our MD simulations reveal a linear pressure-dependence of the shear stress τMD (P,σH2o) [MPa]=0.18P + 50.5 - 1.25σH2o, where σH2o [nm-2] is the surface number density of adsorbed water molecules. This constitutive law is of central importance for our multiscale approach, since it forms a link between MD and elastoplastic contact mechanics calculations. An integration of τMD (P,σH2o) over the real area of contact yields a macroscopic friction coefficient μmacro (σH2o) that allows for a meaningful comparison with friction coefficients μexp≈0.5-0.7 which are in good agreement with the calculated dry friction coefficients μmacro(σH2o=0).For milder experimental loads, our multiscale model suggests that the lower friction states with μexp≈0.2 originate in the presence of physisorbed molecules (e.g. water), which significantly reduce interfacial adhesion.

Author(s):  
Peng-zhe Zhu ◽  
Hui Wang ◽  
Yuan-zhong Hu

Three-dimensional molecular dynamics (MD) simulations have been performed to investigate behaviors of nanoindentation and nano-scratch. The first case concerns the effects of material defect on the nanoindentation of nickel thin film. The defect is modeled by a spherical void embedded in the substrate and located under the surface of indentation. The simulation results reveal that compared to the case without defect, the presence of the void softens the material and allows for larger indentation depth at a given load. MD simulations are then performed for nano-scratch of single crystal copper, with emphasis on the effect of indenter shape (sharp and blunt) on the substrate deformation. The results show that the blunt indenter causes larger deformation region and much more dislocations at both the indentation and scratch stages. It is also found that during the scratching stage the blunt indenter results in larger chip volume in front of the indenter and gives rise to more friction than the sharp indenter. The scope of the simulations has been extended by introducing a multiscale model which couples MD simulations with Finite Element Method (FEM), and multiscale simulations are performed for two-dimensional nanoindentation of copper. The model has been validated by well-consistent load-depth curves obtained from both multiscale and full MD simulations, and by good continuity of deformation observed in the handshake region. The simulations also reveal that indenter radius and indentation velocity significantly affect the nanoindentation behavior. By use of multiscale method, the system size to be explored can be greatly expanded without increasing much computational cost.


2021 ◽  
Author(s):  
Maximilian Ries ◽  
Paul Steinmann ◽  
Sebastian Pfaller

Nano-filled polymers are becoming more and more important to meet the continuously growing requirements of modern engineering problems. The investigation of these composite materials at the molecular level, however, is either prohibitively expensive or just impossible. Multiscale approaches offer an elegant way to analyze such nanocomposites by significantly reducing computational costs compared to fully molecular simulations.When coupling different time and length scales, however, it is in particular important to ensure that the same material description is applied at each level of resolution.The Capriccio method, for instance, couples a particle domain modeled with molecular dynamics (MD) with a finite element based continuum description and has been used i.a. to investigate the effects of nano-sized silica additives embedded in atactic polystyrene (PS). However, a simple hyperelastic constitutive law is used so far for the continuum description which is not capable to fully match the behavior of the particle domain. To overcome this issue and to enable further optimization of the coupling scheme, the material model used for the continuum should be derived directly from pure MD simulations under thermodynamic conditions identical to those used by the Capriccio method.To this end, we analyze the material response of pure PS under uniaxial deformation using strain-controlled MD simulations. Analogously, we perform simulations under pure shear deformation to obtain a comprehensive understanding of the material behavior.As a result, the present PS shows viscoelastic characteristics for small strains, whereas viscoplasticity is observed for larger deformations. The insights gained and data generated are used to select a suitable material model whose parameters have to be identified in a subsequent parameter optimization.


Sign in / Sign up

Export Citation Format

Share Document