REVIEW 3 major objections 5 minor 32 references
Mechanical behaviour of brain-skull interface (meninges) under shear loading through experiment and finite element modelling: Preliminary results
T0 review · 3 major / 5 minor · reviewed 2026-08-03 · deepseek-v4-flash
Pith's one-line read Cohesive layer captures brain-skull interface shear failure
desk verdict Useful new protocol and shear data for the brain-skull interface, but the reported normal traction is not identifiable from shear-only tests. 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 central mechanism is a cohesive zone model (CZM) implementing a traction-separation law with a maximum nominal stress damage-initiation criterion and energy-based damage evolution. This represents the brain-skull interface as a layer of cohesive elements. The brain tissue itself is modeled as a second-order Ogden hyperelastic material whose subject-specific parameters are calibrated first on brain-only samples, then reused in the brain-skull complex simulations. MRI-based segmentation provides the true three-dimensional sample geometry for the finite element meshes, avoiding idealization of sample shape.
What would settle it
Run the same shear protocol on a larger set of sheep brain-skull samples while independently measuring the interface elastic moduli (e.g., small-strain normal and shear tests on intact meninges). If the measured moduli deviate substantially from 61 kPa and 11 kPa, or if calibrated maximum tractions scatter well outside 2.8-3.4 kPa (normal) and 1.8-2.1 kPa (tangential), the claim that the cohesive layer captures interface behavior would be contradicted.
Extended reading notes
Core claim
The paper claims that a cohesive layer with a traction-separation law captures the force-displacement behavior and damage initiation of the brain-skull interface under shear. Failure was observed to initiate within the subarachnoid space, with the pia mater remaining attached to the brain and the dura mater to the skull. The calibrated cohesive properties — maximum normal traction 2.8-3.4 kPa, maximum tangential traction 1.8-2.1 kPa, and fracture energy 0.48-0.7 N/m — were consistent across the three brain-skull complex samples tested. The authors also note that after interface failure, a gradual force decrease occurred in experiments but not in the model, because adhesion between pia and du
Load-bearing premise
The calibrated cohesive tractions and fracture energy assume the elastic moduli of the brain-skull interface (61 kPa normal, 11 kPa shear) from bovine pia-arachnoid literature apply to sheep meninges at the tested strain rate; if those moduli are wrong, the calibrated failure parameters are systematically biased.
Editorial extensions
If this is right
- Computational head models could replace idealized brain-skull contact conditions with a cohesive layer whose parameters are experimentally derived, improving injury-prediction biofidelity.
- The consistent calibrated tractions across samples suggest a reproducible failure envelope for meninges under shear, at least in sheep.
- Failure localization in the subarachnoid space indicates where injury models should place the separation site.
- The framework can be extended to tension and compression to obtain a complete, experimentally grounded description of the brain-skull interface.
- The uncaptured gradual post-failure force decrease points to the need for a pia-dura adhesion or contact model after interface failure.
Reading between the lines
- If these values hold under larger sampling, they imply the meninges are dramatically weaker in shear than in in-plane tension, a distinction head-injury models may need to respect.
- A natural testable extension is to add a contact/friction law between pia and dura after cohesive failure; the paper's own data show this would reproduce the gradual force decay.
- The sheep-to-human transfer is qualitative at best; human meninges are thicker and differ in CSF drainage, so scaling or direct human measurements would be needed before adopting these numbers clinically.
- The uncalibrated elastic moduli taken from bovine pia-arachnoid literature are the main systematic risk; measuring E_nn, E_ss, E_tt on the same samples would strengthen the calibrated failure parameters.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This manuscript presents an experimental-computational framework to characterize the mechanical behavior of the ovine brain-skull interface (meninges) under shear loading. Brain tissue samples and brain-skull complex samples were extracted from three sheep cadaver heads, tested in shear at nominally 0.02/s, and imaged by MRI to build subject-specific finite element meshes. A second-order Ogden hyperelastic model was calibrated to the brain tissue force-displacement curves, and a cohesive zone model (linear elastic traction-separation, maximum stress damage initiation, energy-based evolution) was used for the brain-skull interface, with elastic moduli taken from the literature (Jin et al., bovine pia-arachnoid) and damage parameters calibrated to the brain-skull complex shear tests. The authors report maximum normal tractions of 2.8–3.4 kPa and maximum tangential tractions of 1.8–2.1 kPa, and conclude that a cohesive layer captures the force-displacement and damage initiation of the brain-skull interface.
Significance. If the quantitative results were sound, this would be a useful step toward replacing arbitrary boundary conditions in head FE models with data-derived descriptions of the brain-skull interface. The study has notable strengths: it tests the intact brain-skull complex rather than isolated meningeal layers, uses MRI-derived sample geometries, verifies the subarachnoid-space failure location with camera recordings, and explicitly acknowledges the post-failure mismatch. However, the headline normal-traction values are not identifiable from the shear-only experiments, and the brain tissue validation is in-sample. The framework is promising, but the quantitative claims need substantial revision before the results can be relied upon.
major comments (3)
- [§2.4.1, Eqs. (4)–(5), Table 3] Under the shear loading described in §2.3, the brain-skull interface is approximately parallel to the loading direction. In the uncoupled traction–separation law (Eq. 4), t_n = E_nn ε_n; with no normal separation across the cohesive layer, ε_n ≈ 0, so t_n does not affect the simulated force–displacement response. The damage-initiation criterion (Eq. 5) is then controlled solely by t_s/t_s0. The SLSQP calibration therefore has no sensitivity to t_n, and the values 2.8–3.4 kPa in Table 3 are arbitrary optimizer outcomes, not measured brain-skull interface properties. This directly undermines the abstract's quantitative summary and the claim to have determined the maximum normal traction. Please remove t_n from the headline results or support it with an independent normal-traction experiment/simulation.
- [§2.4.1, Table 3] The elastic moduli of the brain-skull interface, E_nn=61 kPa and E_ss=E_tt=11 kPa, are taken from Jin et al. [15] for the bovine pia-arachnoid complex and used without calibration for the ovine intact meninges at the tested strain rate. Since E_ss governs the pre-failure shear stiffness of the cohesive layer and directly participates in the force balance, an incorrect value will bias the calibrated tangential tractions and fracture energy. The paper should include a sensitivity/uncertainty analysis, or at least a quantitative justification for transferability from bovine pia-arachnoid to ovine whole meninges.
- [§3, Fig. 5, Table 2] The brain tissue Ogden parameters are calibrated by minimizing the difference between model and experimental force–displacement curves for the same samples, and then the agreement is reported as 'confirm the validity of our selection' of the hyperelastic model. This is in-sample validation; it demonstrates consistency but not predictive validity. Please rephrase the claim and, if feasible, provide a hold-out or cross-validation check to strengthen the conclusion.
minor comments (5)
- [Abstract / §3] The abstract states that a cohesive layer 'captures the force-displacement and damage initiation,' but Fig. 6 shows a rapid force drop to zero in the model versus a gradual decrease in experiments after failure. The authors acknowledge this in §3, so the abstract should be qualified to make clear the agreement holds only up to interface failure.
- [Eq. (1)] The second-order Ogden strain energy includes compressibility terms with D_i, but Table 2 reports only μ1, α1, μ2, α2, and μ0; D1 is not given. Since the paper uses a nearly incompressible formulation (Poisson's ratio 0.49), the missing D value prevents full reproduction of the constitutive model.
- [§2.4] No mesh convergence study is reported for the 0.5 mm element size. Given the importance of mesh density on cohesive element response and peak tractions, a brief convergence check or at least a justification of the chosen element size would strengthen the results.
- [Table 3] The column 'Failure force difference percentage (%)' is not defined in the text. It is unclear whether this is the maximum force difference, the difference at the failure point, or the difference in peak force. Please define it explicitly.
- [§2.3 / Fig. 2] The sample coordinate system and the orientation of the brain-skull interface relative to the loading direction are described verbally but could be clearer in Fig. 2. Adding labels for the interface plane and the shear loading direction would help the reader verify the mode-mix assumptions.
Circularity Check
Cohesive-layer agreement is an in-sample fit, and the headline normal traction is unidentifiable from shear-only tests.
-
fitted input called prediction
[Abstract; §2.4.1; §3 (Fig. 6, Table 3)]
"Our results indicate that a cohesive layer captures the force-displacement and damage initiation of the brain-skull interface. ... the nominal tractions (t_n and t_s = t_t, as defined in Equation 4) and fracture energy (G) were the cohesive properties of the brain-skull interface calibrated to match experimental force-displacement results of the brain-skull complex samples under shear loading."
The cohesive parameters (t_s/t_t, G) are determined by SLSQP minimization of the difference between the FE and experimental force-displacement curves; the agreement in Fig. 6 is therefore a least-squares fit residual, not an out-of-sample prediction. The independent content is limited to camera-observed failure location/kinematics that were not part of the objective function. Thus the abstract's 'captures' claim is partially supported only by the fitting data, with a modest independent check.
-
other
[§2.4.1, Eqs. 4–5; Table 3]
"The nominal tractions (t_n and t_s = t_t, as defined in Equation 4) and fracture energy (G) were the cohesive properties of the brain-skull interface calibrated to match experimental force-displacement results of the brain-skull complex samples under shear loading."
The sample is oriented so the brain-skull interface is parallel to the shear loading direction, giving pure tangential loading with normal nominal strain ε_n ≈ 0. Equation 4 then gives t_n = E_nn ε_n ≈ 0, and in the damage-initiation criterion, Equation 5, the term <t_n>/t_n0 is zero. Hence t_n never contributes to the simulated force-displacement response and cannot be identified by the SLSQP calibration. The reported maximum normal tractions of 2.8–3.4 kPa are therefore artifacts of the optimizer rather than properties constrained by the experiments; the headline quantitative result is disconnected from the data by construction.
full rationale
The paper is explicit that the cohesive parameters and the brain-tissue Ogden parameters are obtained by inverse FE calibration against the measured force-displacement curves; this is standard parameter identification rather than a hidden tautology, and the use of MRI-derived geometry adds genuine content. However, the abstract's statement that a cohesive layer 'captures' the force-displacement and damage initiation is an in-sample validation: the same curves used to fit the parameters are then used to demonstrate agreement. This reduces the force-displacement portion of the claim to the minimization itself. A second, independent issue is that the maximum normal traction t_n cannot be identified from shear tests with the interface parallel to the loading direction, because Eq. 4 and the damage criterion Eq. 5 eliminate t_n from the response; the reported 2.8–3.4 kPa values are therefore unsupported by the experiments. This is partly a correctness/identifiability problem rather than a textbook circularity, which is why no step is classified as self-citation or uniqueness smuggling; the self-citations to the authors' prior work are consistency checks, not load-bearing. The acknowledged limitations (six samples, post-failure force-drop mismatch, pia-dura adhesion not modelled) are stated honestly. Overall, the central claim has some independent support from the camera-observed subarachnoid failure location, but two construction-level issues make a partial circularity score appropriate.
Assumptions & free parameters
free parameters (5)
- Ogden shear moduli μ1, μ2 (brain tissue) =
B1: 800, 386.7 Pa; B2: 1210.8, 466.4 Pa; B3: 821.6, 599.9 Pa
- Cohesive maximum normal traction t_n =
3.0, 3.4, 2.8 kPa (S1-S3)
- Cohesive maximum tangential traction t_s=t_t =
2.1, 1.9, 1.8 kPa (S1-S3)
- Cohesive fracture energy G =
0.48, 0.54, 0.70 N/m (S1-S3)
- Interface elastic moduli E_nn, E_ss, E_tt =
E_nn=61 kPa, E_ss=E_tt=11 kPa
assumptions (7)
- domain assumption Sheep brain tissue is nearly incompressible with Poisson's ratio ν=0.49.
- domain assumption The second-order Ogden model (N=2) is an adequate constitutive form for brain tissue under shear up to γ=0.3.
- ad hoc to paper A single cohesive layer with linear elastic traction-separation, maximum stress damage initiation (Eq. 5), and energy-based evolution represents the brain-skull interface.
- ad hoc to paper Literature values for the pia-arachnoid complex (Jin et al.: E_nn=61 kPa, E_ss=E_tt=11 kPa) are valid for the whole ovine meninges at the tested strain rate.
- domain assumption Brain tissue parameters identified from the isolated tissue sample (Part 2) are transferable to the brain tissue within the brain-skull complex sample (Part 1) of the same head.
- domain assumption Cyanoacrylate glue creates a rigid, no-slip bond on loading and base platens.
- domain assumption The optimisation finds a unique, globally optimal parameter set.
Cite this review
Pith. "Pith review of Mechanical behaviour of brain-skull interface (meninges) under shear loading through experiment and finite element modelling: Preliminary results." pith.science (2026). https://pith.science/paper/U44QRHKS
@misc{pith2026251208425,
author = {Pith},
title = {Pith review of: Mechanical behaviour of brain-skull interface (meninges) under shear loading through experiment and finite element modelling: Preliminary results},
year = {2026},
howpublished = {\url{https://pith.science/paper/U44QRHKS}},
note = {Machine review of arXiv:2512.08425}
}
read the original abstract
The brain-skull interface (meninges) plays a critical role in governing brain motion during head impacts, yet computational models often simplify this interface using idealized contact conditions due to limited experimental data. This study presents an improved protocol combining experimental testing and computational modelling to determine the mechanical properties of the brain-skull interface under shear loading. Brain tissue and brain-skull complex samples were extracted from sheep cadaver heads and subjected to shear loading. Magnetic resonance imaging (MRI) was used to obtain accurate 3D geometries of the samples, which were then used to create computational grids (meshes) for simulation of the experiments using finite element (FE) models to determine subject-specific properties of the brain tissue and brain-skull interface. A second-order Ogden hyperelastic model was used for the brain tissue, and a cohesive layer was employed to model the brain-skull interface. Our results indicate that a cohesive layer captures the force-displacement and damage initiation of the brain-skull interface. The calibrated cohesive properties showed consistent patterns across samples, with maximum normal tractions ranging from 2.8-3.4 kPa and maximum tangential tractions from 1.8-2.1 kPa. This framework provides a foundation for improving the biofidelity of computational head models used in injury prediction and neurosurgical planning by replacing arbitrary boundary conditions with formulations derived from experimental data on brain-skull interface (meninges) biomechanical behaviour.
Reference graph
Works this paper leans on
-
[15]
Annals of biomedical engineering, 2014
Jin, X., et al., Constitutive modeling of pia–arachnoid complex. Annals of biomedical engineering, 2014. 42: p. 812-821
2014
-
[1]
International Journal for Numerical Methods in Biomedical Engineering, 2019
Miller, K., et al., Biomechanical modeling and computer simulation of the brain during neurosurgery. International Journal for Numerical Methods in Biomedical Engineering, 2019. 35(10): p. e3250
2019
-
[2]
Annals of biomedical engineering, 2013
Mostayed, A., et al., Biomechanical model as a registration tool for image -guided neurosurgery: evaluation against BSpline registration. Annals of biomedical engineering, 2013. 41: p. 2409-2425
2013
-
[3]
arXiv preprint arXiv:2402.14303, 2024
Huynh, A.T., et al., Open Meshed Anatomy: Towards a comprehensive finite element hexahedral mesh derived from open atlases. arXiv preprint arXiv:2402.14303, 2024
arXiv 2024
-
[4]
2011: Springer
Miller, K., Biomechanics of the Brain. 2011: Springer
2011
-
[5]
Scientific Reports, 2022
Yang, S., et al., Assessment of brain injury characterization and influence of modeling approaches. Scientific Reports, 2022. 12(1): p. 13597
2022
-
[6]
Biomechanics and Modeling in Mechanobiology,
Wang, F., et al., Prediction of brain deformations and risk of traumatic brain injury due to closed-head impact: quantitative analysis of the effects of boundary conditions 15 and brain tissue constitutive model. Biomechanics and Modeling in Mechanobiology,
-
[7]
subdural
Haines, D.E., H.L. Harkey, and O. Al-Mefty, The “subdural” space: a new look at an outdated concept. Neurosurgery, 1993. 32(1): p. 111-120
1993
Show all 32 references
-
[8]
Journal of Neurotrauma, 2021
Walsh, D.R., et al., Mechanical properties of the cranial meninges: a systematic review. Journal of Neurotrauma, 2021. 38(13): p. 1748-1761
2021
-
[9]
American Journal of Physiology - Regulatory, Integrative and Comparative Physiology, 1998
Boulton, M., et al., Raised intracranial pressure increases CSF drainage through arachnoid villi and extracranial lymphatics. American Journal of Physiology - Regulatory, Integrative and Comparative Physiology, 1998. 275(3): p. R889-R896
1998
-
[10]
Comparative medicine, 2020
Kinaci, A., et al., Histologic comparison of the dura mater among species. Comparative medicine, 2020. 70(2): p. 170-175
2020
-
[11]
and J.N.J
Banstola, A. and J.N.J. Reynolds, Mapping sheep to human brain: The need for a sheep brain atlas. Frontiers in veterinary science, 2022. 9: p. 961413
2022
-
[12]
PLoS One, 2014
Hoffmann, A., et al., The ovine cerebral venous system: comparative anatomy, visualization, and implications for translational research. PLoS One, 2014. 9(4): p. e92990
2014
-
[13]
Journal of biomechanics, 2005
Miller, K., Method of testing very soft biological tissues in compression. Journal of biomechanics, 2005. 38(1): p. 153-158
2005
-
[14]
Miller, K. and K. Chinzei, Mechanical properties of brain tissue in tension. Journal of biomechanics, 2002. 35(4): p. 483-490
2002
-
[16]
Fluids and Barriers of the CNS,
Conley Natividad, G., et al., Ex-vivo quantification of ovine pia arachnoid complex biomechanical properties under uniaxial tension. Fluids and Barriers of the CNS,
-
[17]
Acta biomaterialia, 2021
Walsh, D.R., et al., Mechanical characterisation of the human dura mater, falx cerebri and superior sagittal sinus. Acta biomaterialia, 2021. 134: p. 388-400
2021
-
[18]
Acta biomaterialia, 2017
MacManus, D.B., et al., Protection of cortex by overlying meninges tissue during dynamic indentation of the adolescent brain. Acta biomaterialia, 2017. 57: p. 384-394
2017
-
[19]
Towards Determining Mechanical Properties of Brain -Skull Interface Under Tension and Compression
Arzemanzadeh, S., et al. Towards Determining Mechanical Properties of Brain -Skull Interface Under Tension and Compression . in Computational Biomechanics for Medicine: Progress in Research and Applications. 2025. Springer
2025
-
[20]
Destrade, and M.D
Rashid, B., M. Destrade, and M.D. Gilchrist, Influence of preservation temperature on the measured mechanical properties of brain tissue. Journal of biomechanics, 2013. 46(7): p. 1276-1281
2013
-
[21]
34(5): p
Miller, K., How to test very soft biological tissues in extension? Journal of Biomechanics, 2001. 34(5): p. 651-657
2001
-
[22]
Miller, K. and K. Chinzei, Constitutive modelling of brain tissue: experiment and theory. Journal of biomechanics, 1997. 30(11-12): p. 1115-1121
1997
-
[23]
Journal of biomechanics, 2013
Jin, X., et al., A comprehensive experimental study on material properties of human brain tissue. Journal of biomechanics, 2013. 46(16): p. 2795-2801
2013
-
[24]
Shepherd, and D.M
Li, W., D.E.T. Shepherd, and D.M. Espino, Investigation of the compressive viscoelastic properties of brain tissue under time and frequency dependent loading conditions. Annals of Biomedical Engineering, 2021. 49(12): p. 3737-3747. 16
2021
-
[25]
Magnetic resonance imaging, 2012
Fedorov, A., et al., 3D Slicer as an image computing platform for the Quantitative Imaging Network. Magnetic resonance imaging, 2012. 30(9): p. 1323-1341
2012
-
[26]
Yang, K.H. and H. Mao, Modelling of the brain for injury simulation and prevention. Biomechanics of the Brain, 2019: p. 97-133
2019
-
[27]
Forschungsbericht- Deutsche Forschungs - und Versuchsanstalt fur Luft - und Raumfahrt, 1988
Kraft, D., A software package for sequential quadratic programming. Forschungsbericht- Deutsche Forschungs - und Versuchsanstalt fur Luft - und Raumfahrt, 1988
1988
-
[28]
Park, K. and G.H. Paulino, Cohesive zone models: a critical review of traction - separation relationships across fracture surfaces. Applied Mechanics Reviews, 2011. 64(6): p. 060802
2011
-
[29]
The Journal of Adhesion, 2020
Jaillon, A., et al., Mode I cohesive zone model parameters identification and comparison of measurement techniques for robustness to the law shape evaluation. The Journal of Adhesion, 2020. 96(1-4): p. 272-299
2020
-
[30]
Archives of Computational Methods in Engineering, 2019
Budday, S., et al., Fifty shades of brain: a review on the mechanical testing and modeling of brain tissue. Archives of Computational Methods in Engineering, 2019
2019
-
[31]
Mechanical properties of brain –skull interface in compression
Agrawal, S., et al. Mechanical properties of brain –skull interface in compression . in Computational Biomechanics for Medicine: New Approaches and New Applications
-
[32]
Acta biomaterialia, 2017
Budday, S., et al., Mechanical characterization of human brain tissue. Acta biomaterialia, 2017. 48: p. 319-340
2017
Reviewed August 3, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.