Deep Indentation of Hyperelastic Materials Reveals Tip Independent Parabolic Force Depth Response via Strain Energy Delocalization

Jimmy Hsia; Mattia Bacca; Mohammad Shojaeifard; Norman Fleck; Zhenwei Ma

arxiv: 2606.30971 · v2 · pith:QX2AJQBOnew · submitted 2026-06-29 · ❄️ cond-mat.soft

Deep Indentation of Hyperelastic Materials Reveals Tip Independent Parabolic Force Depth Response via Strain Energy Delocalization

Mohammad Shojaeifard , Zhenwei Ma , Jimmy Hsia , Norman Fleck , Mattia Bacca This is my paper

Pith reviewed 2026-07-02 20:05 UTC · model grok-4.3

classification ❄️ cond-mat.soft

keywords deep indentationhyperelastic materialsstrain energy delocalizationtip independenceparabolic force-depthOgden modelsoft materials

0 comments

The pith

Deep indentation of hyperelastic materials produces a parabolic force-depth law F = beta E D^2 that is independent of indenter tip shape.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper establishes that when indentation depth D greatly exceeds indenter radius R, flat and spherical indenters on hyperelastic materials converge to the identical force-depth relation F equals beta times Young's modulus E times D squared. This convergence occurs because strain energy delocalizes into a spheroidal domain whose linear size grows with D, so the activated volume scales as D cubed, stored elastic energy scales as E D cubed, and force as the derivative with respect to D scales as E D squared. A sympathetic reader would care because the prefactor beta depends primarily on the material's strain-stiffening response, is only mildly affected by friction, and is independent of tip radius or shape, enabling extraction of hyperelastic parameters from tests where contact geometry is uncertain or variable.

Core claim

At deep indentation depths D much larger than indenter radius R, flat-punch and spherical indenters on hyperelastic solids produce identical parabolic force-depth curves F = beta E D^2. The prefactor beta is independent of tip radius and shape, mildly sensitive to friction, and determined by the hyperelastic strain-stiffening response. This scaling originates from strain-energy delocalization: the region of significant strain-energy density expands as a spheroid whose linear size scales with D, so the activated volume grows as D cubed, stored energy U scales as E D cubed, and force F equals dU/dD scales as E D squared. Strain-energy density fields far from the contact collapse onto the Bouss

What carries the argument

strain-energy delocalization into a spheroidal domain whose size scales with indentation depth D

Load-bearing premise

Finite-element simulations with the chosen hyperelastic constitutive model accurately capture the dominant strain-energy distribution at large deformations.

What would settle it

Experiments on a well-characterized hyperelastic material that find a force-depth exponent different from 2 or a beta that depends strongly on tip shape at D much larger than R would falsify the tip-independent parabolic scaling.

read the original abstract

Indentation is a practical route for probing soft materials when standard tests are difficult, destructive, or cannot be performed in situ. Conventional indentation is usually interpreted in the shallow-depth regime, where the indentation depth D is small compared with the indenter radius R. In this limit, the response is controlled by local contact geometry and primarily identifies the small-strain Young's modulus E. Here, we show that at deep indentation, D >> R, flat and spherical indenters converge to the same parabolic force-depth law, F = beta E D^2. The coefficient beta is independent of indenter radius and tip shape, only mildly affected by interfacial friction, and controlled by the hyperelastic strain-stiffening response. Finite-element simulations show that this scaling arises from strain-energy delocalization: the region where SED/mu > 0.01 expands into a spheroidal domain whose size scales with D. The activated volume therefore scales as D^3, giving stored elastic energy U ~ E D^3 and force F = dU/dD ~ E D^2. Far from contact, the strain-energy-density fields collapse toward the Boussinesq far-field solution when distances are normalized by a = sqrt(F/E), which scales as D in the deep-indentation regime. These results provide a mechanistic basis for tip-shape independence and link beta to the Ogden strain-stiffening parameter alpha, enabling hyperelastic parameter extraction from deep-indentation data.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

Deep indentation yields tip-independent parabolic F ~ E D^2 from strain-energy delocalization into a D-scaled spheroid in Ogden FE simulations, but lacks reported mesh or domain checks.

read the letter

The main result is that at D much larger than indenter radius, flat and spherical tips on hyperelastic solids both give F = beta E D^2, with beta set by the strain-stiffening parameter rather than tip geometry. Finite-element runs show the region of significant strain-energy density expanding as a spheroid whose size tracks D, so activated volume scales as D^3, stored energy as E D^3, and force as its derivative with D.

What is new is the explicit connection between that volume scaling and the observed tip independence, plus the far-field collapse onto the Boussinesq solution when distances are normalized by sqrt(F/E). The simulations make the mechanism concrete and tie beta directly to the Ogden alpha.

The work is internally consistent on its own terms and gives a usable scaling for large-deformation indentation. The delocalization picture is straightforward once the volume growth is accepted.

The soft spot is that all evidence comes from FE with one constitutive model and no reported mesh-convergence or outer-domain-size studies in the deep regime. If the rate of energy spread depends on boundary conditions or discretization, the claimed independence of tip shape could change. Linking beta to alpha is also partly by construction, since alpha is a fitted constant, so the result functions more as a mapping than an independent derivation.

This is for people who interpret or model indentation on soft gels, elastomers, or tissues at large strains. It supplies both a practical law and a physical account that can guide experiments. The thinking is clear and the subfield engagement is standard.

Send it for peer review. The scaling and mechanism are worth referee scrutiny even if more validation details will be required.

Referee Report

2 major / 0 minor

Summary. The manuscript claims that at deep indentation (D ≫ R) of hyperelastic materials, flat and spherical indenters converge to the same parabolic force-depth law F = β E D². The prefactor β is independent of indenter radius and tip shape, only mildly affected by friction, and controlled by the hyperelastic strain-stiffening response. Finite-element simulations with an Ogden model attribute this to strain-energy delocalization into a spheroidal domain whose size scales with D, yielding activated volume ~ D³, stored energy U ~ E D³, and F = dU/dD ~ E D². Far-field strain-energy-density fields collapse to the Boussinesq solution when distances are normalized by a = sqrt(F/E) ~ D. The work links β to the Ogden α parameter for hyperelastic parameter extraction.

Significance. If the reported scaling and delocalization mechanism are robustly validated, the result would provide a mechanistic basis for tip-shape-independent characterization of hyperelastic materials at large depths, simplifying in-situ testing of soft matter where conventional shallow-indentation analysis fails. The physical picture of D³ volume scaling is a clear strength, though its dependence on the specific constitutive model limits immediate generality.

major comments (2)

[Abstract/Methods] Abstract and Methods: The central claim that the parabolic scaling and tip independence arise from strain-energy delocalization into a D-scaling spheroid rests entirely on finite-element simulations, yet no mesh-convergence studies, domain-size sensitivity checks, or quantitative validation (error bars, comparison to analytic limits) are reported for the deep-indentation regime. This is load-bearing because the D³ volume scaling and resulting β could be artifacts of the chosen Ogden model, far-field boundary conditions, or discretization, directly undermining the reported independence from tip shape.
[Abstract] Abstract: The statement that β is 'controlled by' and 'linked to' the Ogden strain-stiffening parameter α risks circularity, since α is itself a fitted constitutive parameter; the manuscript should demonstrate whether this link is an independent derivation from the delocalization mechanism or merely a re-expression of the model fit.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments. We address each major point below and will revise the manuscript to strengthen the validation of the simulation results and clarify the interpretation of the β-α relationship.

read point-by-point responses

Referee: [Abstract/Methods] Abstract and Methods: The central claim that the parabolic scaling and tip independence arise from strain-energy delocalization into a D-scaling spheroid rests entirely on finite-element simulations, yet no mesh-convergence studies, domain-size sensitivity checks, or quantitative validation (error bars, comparison to analytic limits) are reported for the deep-indentation regime. This is load-bearing because the D³ volume scaling and resulting β could be artifacts of the chosen Ogden model, far-field boundary conditions, or discretization, directly undermining the reported independence from tip shape.

Authors: We agree that the absence of explicit mesh-convergence and domain-size studies is a limitation in the current manuscript. In the revision we will add a dedicated subsection (or supplementary material) reporting mesh refinement results, domain-size variations, and quantitative error metrics for the deep-indentation regime, together with direct comparisons to the far-field Boussinesq solution. These additions will confirm that the reported D³ scaling and tip-independent β are robust. revision: yes
Referee: [Abstract] Abstract: The statement that β is 'controlled by' and 'linked to' the Ogden strain-stiffening parameter α risks circularity, since α is itself a fitted constitutive parameter; the manuscript should demonstrate whether this link is an independent derivation from the delocalization mechanism or merely a re-expression of the model fit.

Authors: The β(α) relation is obtained by systematically varying α in otherwise identical simulations and observing the resulting change in β; we interpret this as a consequence of how strain-stiffening modulates the delocalized energy volume. To remove any appearance of circularity we will revise the abstract and discussion to state explicitly that the link is empirical from the simulations and to frame its use for parameter extraction as a practical outcome rather than an independent analytic derivation. revision: partial

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper's central derivation proceeds from finite-element simulations that demonstrate strain-energy delocalization into a spheroidal domain whose linear size scales with D, producing activated volume ~ D^3, stored energy U ~ E D^3, and thus F = dU/dD ~ E D^2. This scaling is obtained directly from the simulated fields rather than by algebraic rearrangement of an input fit or by self-citation. The reported link between beta and the Ogden alpha parameter is an observed outcome of the chosen constitutive model within those simulations, not a definitional equivalence or a prediction forced by prior fitting of the same quantity. No load-bearing step reduces to its own inputs by construction, and the far-field collapse check is an independent consistency test against the Boussinesq solution.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim rests on the Ogden hyperelastic constitutive model and on the numerical observation that activated volume scales with D cubed; both are introduced without independent experimental confirmation in the abstract.

free parameters (2)

beta
Coefficient in the parabolic law F = beta E D^2; stated to be controlled by alpha but no explicit fitting procedure is given in the abstract.
alpha
Ogden strain-stiffening exponent that controls beta; a standard material parameter fitted to data.

axioms (2)

domain assumption Hyperelastic response is adequately described by the Ogden model with a single strain-stiffening parameter alpha.
Invoked to link the observed beta to material behavior.
domain assumption Finite-element discretization and boundary conditions faithfully reproduce the far-field Boussinesq solution when distances are normalized by sqrt(F/E).
Used to confirm delocalization mechanism.

pith-pipeline@v0.9.1-grok · 5809 in / 1433 out tokens · 25521 ms · 2026-07-02T20:05:05.076309+00:00 · methodology

Review history (2 revisions) →

discussion (0)

Reference graph

Works this paper leans on

5 extracted references · 1 canonical work pages · 1 internal anchor

[1]

2 also illustrates how material parameters are extracted from the normalized force–depth plots

Fig. 2 also illustrates how material parameters are extracted from the normalized force–depth plots. Extrapolating the shallow linear-elastic fits to 𝐷=𝑅 gives 𝐹/𝐷%=16𝐸/9 for spherical indentation, from Eq. (2), and 𝐹/𝐷%=8𝐸/3 for flat-punch indentation, from Eq. (3). These intercepts provide the small-strain elastic modulus 𝐸. Once 𝐸 is determined, the ho...

2025
[2]

Discussion and Conclusions A key limitation of the proposed indentation theory is the assumption that the indented sample behaves like a hyperelastic half-space. The parabolic deep-indentation regime requires both the lateral and vertical sample dimensions to remain large compared with the indenter radius 𝑅, the indentation depth 𝐷, and thus the strain-en...

2025
[3]

Fifty shades of brain: a review on the mechanical testing and modeling of brain tissue

[Budday2020] Budday S, Ovaert TC, Holzapfel GA, Steinmann P, Kuhl E. Fifty shades of brain: a review on the mechanical testing and modeling of brain tissue. Archives of Computational Methods in Engineering. 2020 Sep; 27: 1187–1230. [Delaine-Smith2016] Delaine-Smith RM, Burney S, Balkwill FR, Knight MM. Experimental validation of a flat punch indentation m...

2020
[4]

Hollow Needle Puncture Mechanics for Biopsy Sampling

[Joodaki2018] Joodaki H, Panzer MB. Skin mechanical properties and modeling: A review. Proceedings of the Institution of Mechanical Engineers, Part H: Journal of Engineering in Medicine. 2018 Apr;232(4):323-43. [Lohr2022] Lohr MJ, Sugerman GP, Kakaletsis S, Lejeune E, Rausch MK. An introduction to the Ogden model in biomechanics: benefits, implementation ...

work page internal anchor Pith review Pith/arXiv arXiv 2018
[5]

Spherical indentation method for determining the constitutive parameters of hyperelastic soft materials

17 [Zhang2014] Zhang MG, Cao YP, Li GY, Feng XQ. Spherical indentation method for determining the constitutive parameters of hyperelastic soft materials. Biomechanics and Modeling in Mechanobiology. 2014;13(1):1–11

2014

[1] [1]

2 also illustrates how material parameters are extracted from the normalized force–depth plots

Fig. 2 also illustrates how material parameters are extracted from the normalized force–depth plots. Extrapolating the shallow linear-elastic fits to 𝐷=𝑅 gives 𝐹/𝐷%=16𝐸/9 for spherical indentation, from Eq. (2), and 𝐹/𝐷%=8𝐸/3 for flat-punch indentation, from Eq. (3). These intercepts provide the small-strain elastic modulus 𝐸. Once 𝐸 is determined, the ho...

2025

[2] [2]

Discussion and Conclusions A key limitation of the proposed indentation theory is the assumption that the indented sample behaves like a hyperelastic half-space. The parabolic deep-indentation regime requires both the lateral and vertical sample dimensions to remain large compared with the indenter radius 𝑅, the indentation depth 𝐷, and thus the strain-en...

2025

[3] [3]

Fifty shades of brain: a review on the mechanical testing and modeling of brain tissue

[Budday2020] Budday S, Ovaert TC, Holzapfel GA, Steinmann P, Kuhl E. Fifty shades of brain: a review on the mechanical testing and modeling of brain tissue. Archives of Computational Methods in Engineering. 2020 Sep; 27: 1187–1230. [Delaine-Smith2016] Delaine-Smith RM, Burney S, Balkwill FR, Knight MM. Experimental validation of a flat punch indentation m...

2020

[4] [4]

Hollow Needle Puncture Mechanics for Biopsy Sampling

[Joodaki2018] Joodaki H, Panzer MB. Skin mechanical properties and modeling: A review. Proceedings of the Institution of Mechanical Engineers, Part H: Journal of Engineering in Medicine. 2018 Apr;232(4):323-43. [Lohr2022] Lohr MJ, Sugerman GP, Kakaletsis S, Lejeune E, Rausch MK. An introduction to the Ogden model in biomechanics: benefits, implementation ...

work page internal anchor Pith review Pith/arXiv arXiv 2018

[5] [5]

Spherical indentation method for determining the constitutive parameters of hyperelastic soft materials

17 [Zhang2014] Zhang MG, Cao YP, Li GY, Feng XQ. Spherical indentation method for determining the constitutive parameters of hyperelastic soft materials. Biomechanics and Modeling in Mechanobiology. 2014;13(1):1–11

2014