REVIEW 4 major objections 5 minor 19 references
Localised stress and strain distribution in sliding
T0 review · 4 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read Simulated and measured strain fields match around a sliding nanoscratch, exposing a residual shear band that explains a dislocation line in copper.
desk verdict New residual strain maps of a nanoscratch with a plausible mechanistic story, but the extracted-FIB-lamella relaxation problem and calibrated-model circularity weaken the central quantitative claim. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The machinery is a coupled experimental-simulation loop: nanoscratch testing with a Berkovich tip on (001) copper, HR-EBSD cross-correlation to map residual elastic strain at the surface and in two cross-sectional planes, and a physically based crystal-plasticity user material for Abaqus that solves the same nanoscratch in three dimensions with constitutive parameters calibrated on the same experiments. The interpretive key is the simulated shear-stress distribution, in particular the $\sigma_{13}$ and $\sigma_{23}$ components in the cross-sectional planes: comparing the loaded indentation, loaded scratch, and unloaded scratch isolates the residual positive shear band beneath the scratch wake.
What would settle it
Perform the same CPFE simulation with constitutive parameters obtained from independent tests on the same crystal, for example tensile or microcompression tests on [100] copper, rather than from the scratch itself, and compare the predicted residual elastic strain fields to the HR-EBSD maps. If the field-by-field agreement disappears, the match rests on calibration rather than on the model's physical content. A second check is to scratch at a different normal load, predict the depth of the zero-shear line between the positive and negative bands, and search by cross-sectioning and EBSD for the dislocation trace at that depth.
Extended reading notes
Core claim
The paper argues that the residual elastic strain fields measured by HR-EBSD around a 3 mN Berkovich nanoscratch in (001) single-crystal copper match, in sign and magnitude, the elastic strain fields from a three-dimensional physically based crystal-plasticity finite-element simulation. On the basis of that match, it reads the simulated stress history: during sliding, large shear stresses form ahead of the contact; after unloading, a band of positive shear stress remains beneath the wake, between two opposite-signed shear bands. The paper identifies this residual positive band as the feature that pushes dislocations down ahead of the contact and pulls them back in the wake, so they arrest along the line of zero shear stress that is the dislocation trace observed in earlier sliding-friction work.
Load-bearing premise
The paper's central comparison assumes that residual elastic strains measured by HR-EBSD after unloading can be matched field-for-field to the elastic strains in the CPFE solution, with the CPFE constitutive parameters already calibrated on the same nanoscratch experiments.
Editorial extensions
If this is right
- The residual stress state around a single sliding asperity in copper is specified in three dimensions: an inner compressive zone, edge pile-up, and opposing shear bands below the track.
- The shear-stress history explains why a dislocation trace line forms below a sliding surface at a fixed depth rather than uniformly: dislocations are driven down ahead of the contact and pulled back in the wake, arresting at the zero-shear line.
- The edge pile-up and lattice rotation asymmetry distinguish sliding from indentation, tying the higher scratch hardness of copper to the smaller effective contact area and the stress concentration ahead of the leading edge.
- Because the CPFE model captures the measured residual fields, it can be used to predict where strain localisation will occur at different loads and tip geometries, giving a route toward predicting the transition from plastic redistribution to material detachment.
- Subsurface stress gradients below a sliding contact are steeper than the elastic analytical solution predicts, indicating that plastic deformation concentrates the stress and drives tribolayer formation.
Reading between the lines
- A natural extension is to change the sliding direction on the same crystal: if the residual positive shear band follows the crystallographic slip geometry rather than the contact geometry, the predicted dislocation-line depth and position should shift accordingly.
- Because the validation is conducted at one load, tip orientation, and material, a testable extension is to vary the normal load and confirm that the residual shear band and the observed dislocation line move deeper together.
- The same stress-gradient argument could be carried to brittle or quasi-brittle materials, where the analogous residual shear band might instead nucleate subsurface cracks; searching for such a band at the onset of micro-cracking would test the transferability.
- A direct experimental signature of the proposed mechanism would be the disappearance or migration of the dislocation line after annealing or after repeated scratching, matching the simulated evolution of the residual positive shear band.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports a combined experimental and CPFE study of the residual elastic strain and stress fields around a 3 mN nanoscratch in (001) single-crystal copper produced by a Berkovich indenter sliding along [100]. HR-EBSD is used to map elastic strains on the deformed surface and on two cross-sectional planes, and a 3D crystal plasticity finite element model is used to simulate the loaded and unloaded states. The authors claim that the sign and magnitude of the simulated strain fields correlate well with the measured fields, and they use the simulated stress fields to identify a residual subsurface band of positive shear stress after unloading that they relate to the self-organised dislocation line observed in sliding friction.
Significance. If the central claim is correct, this is a useful three-dimensional picture of the residual stress state produced by a micrometre-scale sliding asperity, and it links a measurable stress gradient to a proposed dislocation self-organisation mechanism. The paper's strengths are its physically based CPFE treatment, multi-plane HR-EBSD measurements, and explicit separation of loaded and unloaded stress states. The main limitations are that the model is calibrated on the same experiments being compared, the HR-EBSD cross-sections may have been measured on relaxed FIB lamellae, the quantitative agreement is not assessed, and the highest-strain regions are excluded. These issues mean that, as presented, the result is a plausible consistency check rather than a validated prediction.
major comments (4)
- [Section 2 (experimental description; Figs. 2-3)] The text states that '3 µm thick lamella were prepared and extracted using a focused ion beam (FIB)' before HR-EBSD mapping. If these lamellae were fully lifted out, the new free surfaces allow elastic relaxation of the residual stress field, so the measured elastic strains are not directly comparable, field by field, with the semi-infinite bulk CPFE solution. The manuscript contains no relaxation correction and does not state whether the cross-sections remained attached to the bulk during EBSD acquisition. This is load-bearing for the Section 2 claim that simulated and measured strain fields correlate well and for the validation of the residual stress band in Figure 4j; please clarify the specimen state or add a relaxation correction.
- [Section 2 (CPFE calibration)] The CPFE constitutive laws were calibrated against the same nanoscratch experiments, with all details in reference [11]. The simulated fields in Figures 1-4 are therefore not independent predictions but outputs of a model fitted to the very data used for comparison. Use of the same dataset for calibration and validation makes the agreement a consistency check; to support the stronger statement that the model 'accurately capture[s]' the measured fields, the authors should provide a hold-out comparison (e.g., a different load, scratch orientation, or scratch length) or explicitly re-frame the claim as a qualitative consistency check.
- [Section 2 and Figures 1-3] The agreement between simulation and HR-EBSD is assessed only visually: no quantitative metric (e.g., correlation coefficient, mean absolute error, or strain-uncertainty-weighted residual) is reported, no error bars are given for the HR-EBSD strains, and the most heavily deformed subsurface region is excluded because of pattern quality. The statement 'the sign and magnitude of the simulated strain fields correlate well' needs a quantitative basis, and the authors should state explicitly how the exclusion of the near-apex region bounds the comparison.
- [Figure 4 and discussion (dislocation self-organisation)] The mechanistic conclusion that the residual positive shear band 'validated by the residual HR-EBSD measurement' causes the dislocation self-organisation line involves a causal step that is not demonstrated: the HR-EBSD comparison is made in a lower-strain region and the model is calibrated to the same experiment. The authors should either compare with an independent dislocation-level simulation (such as the DDP analysis in reference [9]) or present an explicit stress-based criterion and show that the predicted line position matches observations across more than one load or geometry.
minor comments (5)
- [Abstract] The abstract's claim that the model 'can accurately capture the measured elastic and plastic strain fields' is stronger than the evidence presented; consider 'is consistent with' or add a quantitative measure.
- [Throughout] There are typos: 'unkown' on page 2, 'inhmogenous' on page 8, and 'The stress distribution ... is given for the unloaded scratch is given in h., j., l.' in the Figure 4 caption.
- [Figure 4] The colour scale is labelled in GPa; state whether the HR-EBSD stress maps use the same elastic stiffness constants as the CPFE model and how the HR-EBSD zero-strain reference pattern is defined.
- [Section 2] Because reference [11] contains all experimental and constitutive details, the manuscript should include a supplementary table with scratch parameters, tip orientation, copper elastic constants, and calibrated slip-system parameters so that the results can be reproduced without consulting [11].
- [Figures 1-4] The schematics would benefit from marking the scratch direction explicitly on every plot and from indicating the location of the excluded high-distortion region, so that the reader can see which parts of the field are being compared.
Circularity Check
The central 'validation' reduces to a calibrated-model consistency check: the CPFE model was calibrated on the same nanoscratch experiments whose strain fields are then presented as correlating with simulation.
-
fitted input called prediction
[Section 2, 'Experimental and simulation' paragraph, sentence preceding Figures 1–3 and the 'correlate well' statement]
"The constitutive laws in the model were calibrated against the experiments. Full details of the nanoscratch experiments, cross-section preparation, EBSD/HR-EBSD analysis, and the constitutive laws and parameters used in the CPFE model are provided in [11] ... In all cases the sign and magnitude of the simulated strain fields correlate well with those measured experimentally."
The simulated elastic strain fields in Figures 1–3 and the residual shear-stress band in Figure 4 are outputs of a CPFE model whose constitutive laws were calibrated against the same nanoscratch experiments that produced the HR-EBSD measurements. The statement 'correlate well' therefore describes how well the calibrated model reproduces its own calibration data or companion data from the same test, not an independent prediction. No withheld data set or parameter-free test is offered, so the central validation is a consistency check rather than a first-principles result.
-
self citation load bearing
[Section 2, reference to [11] for all model parameters and calibration details]
"Full details of the nanoscratch experiments, cross-section preparation, EBSD/HR-EBSD analysis, and the constitutive laws and parameters used in the CPFE model are provided in [11] where we reported the lattice rotation fields from these experiments."
The load-bearing premise that the CPFE model faithfully represents the deformation is delegated entirely to the authors' own prior paper [11], which contains both the experiments and the calibration. The present paper supplies no independent verification, code reproduction, or parameter-free derivation, so the chain of support closes on the same data and authors: the calibrated model from [11] is applied to the experiments from [11] and then reported as agreement. This makes the apparent corroboration of the simulated stress interpretation partially self-referential.
full rationale
The paper contains two related circularity elements. First, the comparison between simulated and measured elastic strain fields is presented as a successful correlation, but the simulation model was explicitly calibrated against the same experiments. This is a fitted-input-called-prediction pattern: the model output is not a prediction of a withheld quantity but a reproduction of data used in calibration, so the agreement in Figures 1–4 is expected to a significant degree. Second, all constitutive parameters and calibration details are imported from the authors' prior paper [11], which is also based on the same experiments; this self-citation is load-bearing because the present paper's central interpretation of the residual shear-stress band depends on the fidelity of that externally supplied, unverified calibration. A separate experimental concern—that the HR-EBSD cross-sections were obtained from extracted 3 um FIB lamellae, whose free surfaces relax the residual stress field—is a validity threat to the field-by-field comparison, but it is a correctness issue rather than a circularity issue and does not affect this score. No further circularity patterns, such as renaming known results or importing uniqueness theorems, are present. The derivation does not literally define the output in terms of the input, but the central predictive claim is weakened to a consistency check, which warrants a score of 6.
Assumptions & free parameters
free parameters (1)
- CPFE constitutive parameters =
not stated in this paper (calibrated in [11])
assumptions (4)
- domain assumption Crystal plasticity constitutive equations (Dunne et al. UMAT) describe single-crystal copper deformation
- domain assumption HR-EBSD cross-correlation measures elastic strain relative to a reference pattern with sufficient sensitivity
- domain assumption Berkovich tip geometry and edge-forward orientation in simulation match experiment
- domain assumption Residual elastic strain fields from HR-EBSD are comparable to unloaded CPFE elastic strains
Cite this review
Pith. "Pith review of Localised stress and strain distribution in sliding." pith.science (2026). https://pith.science/paper/HXE5P47Z
@misc{pith2026241115815,
author = {Pith},
title = {Pith review of: Localised stress and strain distribution in sliding},
year = {2026},
howpublished = {\url{https://pith.science/paper/HXE5P47Z}},
note = {Machine review of arXiv:2411.15815}
}
read the original abstract
In this paper, we present a comprehensive analysis of the contact mechanics associated with a micron-sized sliding asperity, which plays a crucial role in the abrasive wear processes. Utilising nanoscratch testing, we experimentally investigate the deformation and employ High-Resolution Electron Backscatter Diffraction (HR-EBSD) to characterise the resulting strain fields at various locations in the residual nanoscratch. To simulate these experiments, we utilise a physically-based Crystal Plasticity Finite Element (CPFE) model, enabling a three-dimensional simulation that can accurately capture the measured elastic and plastic strain fields around the sliding contact. This knowledge serves as a foundation from which we may be able to discern the multi-physical processes governing micro-scale wear phenomena.
Figures
Reference graph
Works this paper leans on
- [11]
-
[9]
Y. Xu, F. Ruebeling, D. S. Balint, C. Greiner, D. Dini, On the origin of microstructural discontinuities in sliding contacts: A discrete dislocation plasticity analysis, International Journal of Plasticity 138 (2021) 102942
work page 2021
-
[1]
Jost, Economic impact of tribology, Proc mechanical failures preven- tion group (1976) 117–139
P. Jost, Economic impact of tribology, Proc mechanical failures preven- tion group (1976) 117–139
work page 1976
-
[2]
K. Holmberg, A. Erdemir, Influence of tribology on global energy con- sumption, costs and emissions, Friction 5 (2017) 263–284
work page 2017
-
[3]
S. Lim, M. Ashby, J. Brunton, Wear-rate transitions and their relation- ship to wear mechanisms, Acta metallurgica 35 (6) (1987) 1343–1348
work page 1987
-
[4]
S. C. Lim, M. Ashby, Overview no. 55 wear-mechanism maps, Acta metallurgica 35 (1) (1987) 1–24
work page 1987
-
[5]
F. P. Bowden, D. Tabor, The friction and lubrication of solids, Vol. 1, Oxford university press, 2001
work page 2001
- [6]
Show all 19 references
-
[7]
Greiner, Z
C. Greiner, Z. Liu, R. Schneider, L. Pastewka, P. Gumbsch, The origin of surface microstructure evolution in sliding friction, Scripta Materialia (2018)
2018
-
[8]
Greiner, Z
C. Greiner, Z. Liu, L. Strassberger, P. Gumbsch, Sequence of stages in the microstructure evolution in copper under mild reciprocating tribo- logical loading, ACS applied materials & interfaces 8 (24) (2016) 15809– 15819
2016
-
[10]
Brazil, G
O. Brazil, G. M. Pharr, Direct observation of partial interface slip in micrometre-scale single asperity contacts, Tribology International 155 (2021) 106776. 9
2021
-
[12]
Kareer, X
A. Kareer, X. Hou, N. M. Jennett, S. V. Hainsworth, The existence of a lateral size effect and the relationship between indentation and scratch hardness in copper, Philosophical Magazine 96 (32-34) (2016) 3396–3413
2016
-
[13]
Kareer, X
A. Kareer, X. Hou, N. M. Jennett, S. Hainsworth, The interaction be- tween lateral size effect and grain size when scratching polycrystalline copper using a berkovich indenter, Philosophical Magazine 96 (32-34) (2016) 3414–3429
2016
-
[14]
G. Hamilton, Explicit equations for the stresses beneath a sliding spheri- cal contact, Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science 197 (1) (1983) 53–59
1983
-
[15]
A. J. Wilkinson, G. Meaden, D. J. Dingley, High-resolution elastic strain measurement from electron backscatter diffraction patterns: New levels of sensitivity, Ultramicroscopy 106 (4-5) (2006) 307–313
2006
-
[16]
A. J. Wilkinson, D. Randman, Determination of elastic strain fields and geometrically necessary dislocation distributions near nanoindents using electron back scatter diffraction, Philosophical magazine 90 (9) (2010) 1159–1177
2010
-
[17]
Britton, A
T. Britton, A. J. Wilkinson, High resolution electron backscatter diffrac- tion measurements of elastic strain variations in the presence of larger lattice rotations, Ultramicroscopy 114 (2012) 82–95
2012
-
[18]
Dunne, D
F. Dunne, D. Rugg, A. Walker, Lengthscale-dependent, elastically anisotropic, physically-based hcp crystal plasticity: Application to cold- dwell fatigue in ti alloys, International Journal of Plasticity 23 (6) (2007) 1061–1083
2007
-
[19]
Dunne, R
F. Dunne, R. Kiwanuka, A. Wilkinson, Crystal plasticity analysis of micro-deformation, lattice rotation and geometrically necessary disloca- tion density, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences 468 (2145) (2012) 2509–2531. 10
2012
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.