We develop an approach that combines the power of nonlinear dynamics with the evolution equations for the mobile and immobile dislocation densities and force to explain force fluctuations in nanoindentation experiments. The model includes nucleation, multiplication, and propagation thresholds for mobile dislocations, and other well known dislocation transformation mechanisms. The model predicts all the generic features of nanoindentation such as the Hertzian elastic branch followed by several force drops of decreasing magnitudes, and residual plasticity after unloading. The stress corresponding to the elastic force maximum is close to the yield stress of an ideal solid. The predicted values for all the quantities are close to those reported by experiments. Our model allows us to address the indentation-size effect including the ambiguity in defining the hardness in the force drop dominated regime. At large indentation depths, the hardness remains nearly constant with a marginal decreasing trend.

• Received 6 December 2013

DOI:https://doi.org/10.1103/PhysRevB.90.094104