In this paper, molecular dynamic (MD) simulation was employed to simulate the AFM tip-based nanoscratching process on multilayer graphene to investigate graphene deformation and coefficient of friction (COF). MD Simulation results showed cross-linking structures were created at the interface of two-layer graphene during nanoindentation and nanoscratching. The loading force had a characteristic of a periodic wave crest-trough transition. The COF of the diamond tip fluctuated at around 0.15 during nanoscratching process. © 2021 Proceedings of the 21st International Conference of the European Society for Precision Engineering and Nanotechnology, EUSPEN 2021. All rights reserved.