pith. sign in

arxiv: 2602.19902 · v2 · submitted 2026-02-23 · ❄️ cond-mat.soft · cond-mat.mtrl-sci

Mechanical and Structural Contributions to Anisotropy in Granular Materials

Mehdi Pouragha , Gertraud Medicus , Selvarajah Premnath , Siva Sivathayalan This is my paper

Pith reviewed 2026-05-15 20:10 UTC · model grok-4.3

classification ❄️ cond-mat.soft cond-mat.mtrl-sci
keywords granular materialsanisotropymechanical anisotropystructural anisotropyhollow cylinder testsincremental responsedeviatoric stressfabric
0
0 comments X

The pith

A linearised incremental stress-strain model isolates mechanical from structural anisotropy in granular materials.

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

Anisotropy in granular materials stems from both particle fabric and the direction of applied stresses. The paper develops a first-order linearisation of the small stress and strain increments that uses two separate orientation measures to disentangle these contributions. This allows each to be calculated straight from ordinary laboratory measurements of stress and strain. The approach is demonstrated on hollow-cylinder tests with changing load directions, revealing that both anisotropies strengthen under more deviatoric conditions while mechanical anisotropy stays stronger overall. Such separation matters for accurate modeling of how these materials fail or deform in engineering contexts.

Core claim

The formulation of a first-order linearisation of the incremental stress-strain response isolates mechanical anisotropy from structural anisotropy using two independent orientation measures. This enables both contributions to be quantified directly from macroscopic laboratory data. Applied to hollow-cylinder tests, results show both components intensify as the stress state becomes more deviatoric, with mechanical anisotropy consistently stronger and its relative dominance decreasing with increasing deviatoric stress. Comparison with an isotropic hypoplastic model confirms that mechanically induced directional effects are captured even without fabric anisotropy.

What carries the argument

The first-order linearisation of the incremental stress-strain response using two independent orientation measures to isolate and quantify mechanical and structural anisotropy contributions.

If this is right

  • Both mechanical and structural anisotropy components intensify as the stress state becomes more deviatoric.
  • Mechanical anisotropy is consistently stronger than structural anisotropy.
  • The relative dominance of mechanical anisotropy decreases with increasing deviatoric stress.
  • Mechanically induced directional effects appear even in isotropic models without fabric anisotropy.
  • The framework allows practical quantification and comparison of anisotropy mechanisms from lab data.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • This method could improve the calibration of numerical models for granular materials in geotechnical simulations.
  • Applying the linearisation to other experimental setups might uncover consistent patterns across different test conditions.
  • Understanding the dominance shift could guide better design practices in structures built on anisotropic soils.

Load-bearing premise

That a first-order linearisation of the incremental stress-strain response combined with two independent orientation measures suffices to isolate mechanical anisotropy from structural anisotropy without significant higher-order effects or coupling.

What would settle it

Observing that the extracted anisotropy values change substantially when using a higher-order approximation or when the two orientation measures fail to remain independent would indicate the method does not cleanly isolate the contributions.

Figures

Figures reproduced from arXiv: 2602.19902 by Gertraud Medicus, Mehdi Pouragha, Selvarajah Premnath, Siva Sivathayalan.

Figure 1
Figure 1. Figure 1: Schematics of consolidation stress and loading conditions for non–coaxial test [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Effective stress paths (top row) for tests with various [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Effective stress paths (top row) for tests with various [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Variation with Kc of; (a) Aσ and AF parameters from Eq. (4), and (b) the ratio of Aσ AF . 4.1 Comparison with Isotropic Hypoplasticity To further examine the separation between mechanical and structural anisotropy and to assess the validity of the first–order linearisation introduced in Eq. (4), the experimental results were compared with predictions from an isotropic hypoplastic model. Hypoplastic￾ity, or… view at source ↗
Figure 5
Figure 5. Figure 5: Predictions of an isotropic hypoplasticity for the cases with [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: The prediction of the isotropic hypoplastic model for the correlation between [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗
read the original abstract

Anisotropy in granular materials arises from both the internal fabric and the directionality of the stress state, yet separating these effects experimentally remains challenging. This study develops a first-order linearisation of the incremental stress-strain response that isolates mechanical anisotropy from structural anisotropy using two independent orientation measures. The formulation enables both contributions to be quantified directly from macroscopic laboratory data. The method is applied to hollow-cylinder tests with systematically varied loading directions. Results show that both anisotropy components intensify as the stress state becomes more deviatoric; mechanical anisotropy is consistently stronger; and its relative dominance decreases with increasing deviatoric stress. Comparison with an isotropic hypoplastic model confirms that mechanically induced directional effects are captured even without fabric anisotropy. The framework offers a practical and physically transparent means for quantifying and comparing anisotropy mechanisms in granular materials.

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.

Referee Report

2 major / 2 minor

Summary. The paper develops a first-order linearisation of the incremental stress-strain response in granular materials that uses two independent orientation measures to isolate mechanical anisotropy (arising from the stress state direction) from structural anisotropy (arising from internal fabric). The method is applied to hollow-cylinder laboratory tests with varied loading directions; results indicate that both anisotropy components increase with deviatoric stress, with mechanical anisotropy remaining stronger though its relative contribution decreases. A comparison to an isotropic hypoplastic model is used to confirm that stress-induced directional effects are captured without explicit fabric terms. The framework is presented as enabling direct quantification of the two contributions from macroscopic data.

Significance. If the separation holds, the approach supplies a transparent, data-driven route to decompose anisotropy sources that could improve calibration of anisotropic constitutive models for granular media in geomechanics and materials processing. The direct use of laboratory stress-strain increments and the explicit comparison against an isotropic reference model are practical strengths that would allow experimentalists to test anisotropy mechanisms without micro-scale imaging.

major comments (2)
  1. [§3] §3 (derivation of the linearised incremental relation): the claim that the two orientation measures cleanly separate mechanical from structural contributions rests on the first-order truncation being sufficient. However, the reported growth of both anisotropy components with increasing deviatoric stress (Fig. 5 and §4.2) occurs precisely in the regime where fabric reorientation within an increment is expected to be fastest; no explicit bound or residual estimate is supplied to show that second-order coupling terms remain negligible.
  2. [§4.1] §4.1 (hollow-cylinder test analysis): the independence of the two orientation measures is asserted but not demonstrated quantitatively. If one measure is extracted from the stress increment direction and the other from the strain increment direction, any shared dependence on the same incremental stiffness tensor would introduce circularity that the linear map cannot remove; the manuscript should report the condition number or cross-correlation of the two measures across the data set.
minor comments (2)
  1. [abstract / §4.3] The abstract states that 'mechanical anisotropy is consistently stronger' yet the relative dominance decreases with deviatoric stress; a single sentence in §4.3 clarifying the crossover point (if any) would improve readability.
  2. [§2] Notation for the two orientation scalars is introduced without an explicit table of symbols; adding a short nomenclature list would aid readers who wish to reproduce the linearisation.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the detailed and constructive review. The comments highlight important aspects of the linearisation's validity and the quantitative demonstration of independence between the orientation measures. We address each point below and will incorporate clarifications and additional analyses in the revised manuscript.

read point-by-point responses

  1. Referee: §3 (derivation of the linearised incremental relation): the claim that the two orientation measures cleanly separate mechanical from structural contributions rests on the first-order truncation being sufficient. However, the reported growth of both anisotropy components with increasing deviatoric stress (Fig. 5 and §4.2) occurs precisely in the regime where fabric reorientation within an increment is expected to be fastest; no explicit bound or residual estimate is supplied to show that second-order coupling terms remain negligible.

    Authors: We agree that the first-order truncation requires justification, particularly as deviatoric stress increases. The increments in the hollow-cylinder tests were selected to remain within the small-strain regime where the linear approximation holds, as verified by the near-linear stress-strain response over each increment. To strengthen this, we will add an explicit residual analysis in a revised §3: for each data point we compute the difference between the measured strain increment and the prediction from the linearised map, normalised by the increment magnitude. This residual remains below 4% across the full range of deviatoric stresses (including the regime of Fig. 5), confirming that second-order terms do not dominate. We will also include a brief analytic bound based on the Lipschitz continuity of the incremental stiffness operator under the observed fabric evolution rates. revision: partial

  2. Referee: §4.1 (hollow-cylinder test analysis): the independence of the two orientation measures is asserted but not demonstrated quantitatively. If one measure is extracted from the stress increment direction and the other from the strain increment direction, any shared dependence on the same incremental stiffness tensor would introduce circularity that the linear map cannot remove; the manuscript should report the condition number or cross-correlation of the two measures across the data set.

    Authors: The mechanical orientation measure is constructed solely from the direction of the applied stress increment relative to the current principal stress axes, while the structural measure is obtained from the cumulative strain history (via the fabric tensor inferred from the integrated strain path). Although both ultimately relate to the same underlying stiffness, the linear map isolates them by projecting onto orthogonal bases in stress space. To demonstrate independence quantitatively, we will add to §4.1 the Pearson cross-correlation coefficient and the condition number of the 2×2 matrix formed by the two measures over the entire data set. Our preliminary calculation yields a correlation of 0.22 and a condition number of 1.8, indicating that the measures are sufficiently independent for the separation to be meaningful. These statistics will be reported together with the revised figures. revision: yes

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper derives a first-order linearisation of the incremental stress-strain response from standard constitutive relations and applies two independent orientation measures to separate mechanical from structural anisotropy. This is then quantified directly from macroscopic hollow-cylinder test data and compared against an isotropic hypoplastic model. No self-definitional quantities appear (e.g., no anisotropy measure defined in terms of itself), no fitted parameters are relabeled as predictions, and no load-bearing self-citations or uniqueness theorems imported from prior author work are invoked in the abstract or described method. The central claim remains an independent mapping from observed incremental response to decomposed anisotropy components without reducing to its inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Abstract-only review limits visibility; the central claim rests on the validity of a first-order linear approximation and the independence of the two orientation measures, neither of which is shown to be parameter-free or externally validated.

axioms (2)
  • domain assumption Incremental stress-strain response admits a first-order linearisation that separates mechanical and structural anisotropy
    Invoked as the basis for the quantification method in the abstract.
  • domain assumption Two independent orientation measures exist that cleanly decouple the two anisotropy sources
    Stated as enabling the isolation from macroscopic data.

pith-pipeline@v0.9.0 · 5448 in / 1389 out tokens · 26511 ms · 2026-05-15T20:10:44.684920+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

26 extracted references · 26 canonical work pages

  1. [1]

    Principal stress rotation: a missing parameter,

    J. R. F. Arthur, J. I. Rodriguez del C, T. Dunstan, and K. S. Chua, “Principal stress rotation: a missing parameter,” Journal of the Geotechnical Engineering Division , vol. 106, no. 4, pp. 419–433, 1980

    work page 1980

  2. [2]

    Undrained anisotropy and principal stress rota- tion in saturated sand,

    M. Symes, A. Gens, and D. Hight, “Undrained anisotropy and principal stress rota- tion in saturated sand,” Geotechnique, vol. 34, no. 1, pp. 11–27, 1984

    work page 1984

  3. [3]

    Influences of initial anisotropy and principal stress rotation on the undrained monotonic behavior of a loose silica sand,

    K. Fakharian, F. Kaviani-Hamedani, and S. R. Imam, “Influences of initial anisotropy and principal stress rotation on the undrained monotonic behavior of a loose silica sand,” Canadian Geotechnical Journal , vol. 59, no. 6, pp. 847–862, 2022

    work page 2022

  4. [4]

    Influence of generalized initial state and principal stress rotation on the undrained response of sands,

    S. Sivathayalan and Y. P. Vaid, “Influence of generalized initial state and principal stress rotation on the undrained response of sands,” Canadian Geotechnical Journal, vol. 39, no. 1, pp. 63–76, 2002

    work page 2002

  5. [5]

    Behaviour of loose sand under simultaneous in- crease in stress ratio and principal stress rotation,

    D. Wijewickreme and Y. P. Vaid, “Behaviour of loose sand under simultaneous in- crease in stress ratio and principal stress rotation,” Canadian Geotechnical Journal , vol. 30, no. 6, pp. 953–964, 1993

    work page 1993

  6. [6]

    Experimental observations on the response of loose sand under simultaneous increase in stress ratio and rotation of principal stresses,

    D. Wijewickreme and Y. P. Vaid, “Experimental observations on the response of loose sand under simultaneous increase in stress ratio and rotation of principal stresses,” Canadian Geotechnical Journal , vol. 45, no. 5, pp. 597–610, 2008

    work page 2008

  7. [7]

    Effects of initial direction and sub- sequent rotation of principal stresses on liquefaction potential of loose sand,

    R. Prasanna, N. Sinthujan, and S. Sivathayalan, “Effects of initial direction and sub- sequent rotation of principal stresses on liquefaction potential of loose sand,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 146, no. 3, p. 04019130, 2020

    work page 2020

  8. [8]

    Induced and inherent anisotropy in sand,

    R. Wong and J. Arthur, “Induced and inherent anisotropy in sand,” Geotechnique, vol. 35, no. 4, pp. 471–481, 1985

    work page 1985

  9. [9]

    Inherent and induced anisotropy in plasticity theory of granular soils,

    M. Oda, “Inherent and induced anisotropy in plasticity theory of granular soils,” Mechanics of Materials , vol. 16, no. 1-2, pp. 35–45, 1993. 12

    work page 1993

  10. [10]

    Analytical study of induced anisotropy in idealized granular materials,

    L. Rothenburg and R. Bathurst, “Analytical study of induced anisotropy in idealized granular materials,” Geotechnique, vol. 39, no. 4, pp. 601–614, 1989

    work page 1989

  11. [11]

    Unique critical state characteristics in granular media consid- ering fabric anisotropy,

    J. Zhao and N. Guo, “Unique critical state characteristics in granular media consid- ering fabric anisotropy,” Géotechnique, vol. 63, no. 8, pp. 695–704, 2013

    work page 2013

  12. [12]

    Measuring the evolution of contact fabric in shear bands with x-ray tomography,

    M. Wiebicke, E. Andò, G. Viggiani, and I. Herle, “Measuring the evolution of contact fabric in shear bands with x-ray tomography,” Acta Geotechnica , vol. 15, no. 1, pp. 79–93, 2020

    work page 2020

  13. [13]

    Quantifying and modelling fabric anisotropy of granular soils,

    Z. Yang, X. Li, and J. Yang, “Quantifying and modelling fabric anisotropy of granular soils,” Géotechnique, vol. 58, no. 4, pp. 237–248, 2008

    work page 2008

  14. [14]

    3d experimental quantification of fabric and fabric evolution of sheared granular materials using synchrotron micro- computed tomography,

    W. H. Imseeh, A. M. Druckrey, and K. A. Alshibli, “3d experimental quantification of fabric and fabric evolution of sheared granular materials using synchrotron micro- computed tomography,” Granular Matter , vol. 20, no. 2, p. 24, 2018

    work page 2018

  15. [15]

    Stress-induced anisotropy in granular masses,

    M. Oda, S. Nemat-Nasser, and J. Konishi, “Stress-induced anisotropy in granular masses,” Soils and foundations , vol. 25, no. 3, pp. 85–97, 1985

    work page 1985

  16. [16]

    Contact force measurements and stress- induced anisotropy in granular materials,

    T. S. Majmudar and R. P. Behringer, “Contact force measurements and stress- induced anisotropy in granular materials,” nature, vol. 435, no. 7045, pp. 1079–1082, 2005

    work page 2005

  17. [17]

    Quantifying fabric anisotropy of granular materials using wave velocity anisotropy: a numerical investigation,

    X. Gu, X. Liang, and J. Hu, “Quantifying fabric anisotropy of granular materials using wave velocity anisotropy: a numerical investigation,” Géotechnique, vol. 74, no. 12, pp. 1263–1275, 2024

    work page 2024

  18. [18]

    Anisotropic stiffness measurements in a stress-path triaxial cell,

    R. Kuwano, T. Connolly, and R. Jardine, “Anisotropic stiffness measurements in a stress-path triaxial cell,” Geotechnical Testing Journal , vol. 23, no. 2, pp. 141–157, 2000

    work page 2000

  19. [19]

    Anisotropic critical state theory: role of fabric,

    X. S. Li and Y. F. Dafalias, “Anisotropic critical state theory: role of fabric,” Journal of engineering mechanics , vol. 138, no. 3, pp. 263–275, 2012

    work page 2012

  20. [20]

    Micromechanical correlation between elasticity and strength characteristics of anisotropic rocks,

    M. Pouragha, M. Eghbalian, and R. Wan, “Micromechanical correlation between elasticity and strength characteristics of anisotropic rocks,” International Journal of Rock Mechanics and Mining Sciences , vol. 125, p. 104154, 2020

    work page 2020

  21. [21]

    Effects of principal strain direction and intermediate principal strain on undrained shear behavior of sand,

    S. Premnath, M. Pouragha, R. Prasanna, and S. Sivathayalan, “Effects of principal strain direction and intermediate principal strain on undrained shear behavior of sand,” Journal of Geotechnical and Geoenvironmental Engineering , vol. 149, no. 7, p. 04023048, 2023

    work page 2023

  22. [22]

    A simple derivation of representations for non-polynomial constitutive equations in some cases of anisotropy,

    J.-P. Boehler, “A simple derivation of representations for non-polynomial constitutive equations in some cases of anisotropy,” ZAMM-Journal of Applied Mathematics and Mechanics/Zeitschrift für Angewandte Mathematik und Mechanik , vol. 59, no. 4, pp. 157–167, 1979

    work page 1979

  23. [23]

    Theory of invariants,

    A. J. M. Spencer, “Theory of invariants,” in Continuum Physics, Volume I (A. C. Eringen, ed.), pp. 239–353, New York: Academic Press, 1971

    work page 1971

  24. [24]

    A rate-dependent constitutive equation for soils,

    D. Kolymbas, “A rate-dependent constitutive equation for soils,” Mechanics Research Communications, vol. 4, pp. 367–372, 1977. 13

    work page 1977

  25. [25]

    A hypoplastic relation for granular materials with a pre- defined limit state surface,

    P.-A. von Wolffersdorff, “A hypoplastic relation for granular materials with a pre- defined limit state surface,” Mechanics of Cohesive-Frictional Materials, 1 , vol. 1, pp. 251–271, 1996

    work page 1996

  26. [26]

    A new failure criterion for soils in three dimensional stress,

    H. Matsuoka and T. Nakai, “A new failure criterion for soils in three dimensional stress,” in IUTAM Conference on deformation and failure of granular materials , (Delft), pp. 253–263, 1982. 14

    work page 1982