pith. sign in

arxiv: 2604.25919 · v1 · submitted 2026-02-01 · 🧮 math.GM

Computational modeling of crack-tip fields in transversely isotropic strain-limiting solids subjected to piecewise linear slope loads

S. M. Mallikarjunaiah , Saugata Ghosh , Dambaru Bhatta This is my paper

Pith reviewed 2026-05-16 08:34 UTC · model grok-4.3

classification 🧮 math.GM
keywords strain-limiting solidstransversely isotropic materialscrack-tip fieldspiecewise linear loadsfinite element methodnonlinear constitutive relationsfracture mechanics
0
0 comments X

The pith

Piecewise linear slope loads allow computation of finite crack-tip fields in transversely isotropic strain-limiting solids via a nonlinear stress-strain relation.

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

The paper studies crack-tip fields inside transversely isotropic elastic solids whose strain is limited by a nonlinear relation to stress. Boundary loads are imposed so that vertical displacement changes linearly with opposite slopes on each half of the top and bottom edges. This setup is solved as a quasi-linear elliptic problem in the displacement field by a continuous Galerkin finite element method together with Picard iteration. Results are obtained for two fiber orientations. The work shows that these loads form a practical way to examine how strain-limiting behavior interacts with crack mechanics without producing infinite strains.

Core claim

A nonlinear constitutive relation that ties Cauchy stress directly to linearized strain removes non-physical strain singularities at the crack tip in a transversely isotropic solid; when this relation is paired with piecewise linear slope boundary conditions on displacement, a continuous Galerkin finite element discretization with Picard linearization yields bounded fields for two distinct fiber orientations.

What carries the argument

The nonlinear constitutive framework relating Cauchy stress to linearized strain, which removes strain singularities at the crack tip and converts the problem into a quasi-linear elliptic boundary-value problem in displacement.

If this is right

  • Vertical displacement prescribed with opposite linear slopes on each half of the top and bottom boundaries produces realistic loading for crack studies.
  • Computed displacement fields remain finite at the crack tip for both fiber orientations examined.
  • The continuous Galerkin method with Picard iteration converges for the quasi-linear system arising from the nonlinear constitutive law.
  • The configuration supplies a flexible platform for exploring strain-limiting effects on crack-tip mechanics in transversely isotropic media.

Where Pith is reading between the lines

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

  • The same nonlinear law and loading strategy could be applied to other classes of anisotropic materials to test whether singularities are suppressed.
  • Time-dependent or cyclic versions of the piecewise slope loads might reveal fatigue behavior in strain-limiting solids.
  • Direct comparison of the computed fields against laboratory measurements on fiber-reinforced polymers would provide an independent check on the model.

Load-bearing premise

The nonlinear relation between Cauchy stress and linearized strain is sufficient to eliminate non-physical strain singularities at the crack tip.

What would settle it

A finite element run or physical test that produces unbounded strain values at the crack tip under the stated piecewise linear slope loads would refute the central claim.

Figures

Figures reproduced from arXiv: 2604.25919 by Dambaru Bhatta, Saugata Ghosh, S. M. Mallikarjunaiah.

Figure 1
Figure 1. Figure 1: Schematic representation of the problem geometry and boundary conditions. A rectan [PITH_FULL_IMAGE:figures/full_fig_p009_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Evolution of the adaptive mesh refinement. Left: The initial coarse mesh (Cycle 0). [PITH_FULL_IMAGE:figures/full_fig_p013_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Advanced stages of adaptive mesh refinement. Left: The mesh at Cycle 8. Right: [PITH_FULL_IMAGE:figures/full_fig_p013_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Convergence of the vertical displacement component [PITH_FULL_IMAGE:figures/full_fig_p014_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Evolution of crack-tip field profiles along the crack line ( [PITH_FULL_IMAGE:figures/full_fig_p015_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Comparison of the growth rates for crack-tip stress and strain across adaptive refinement [PITH_FULL_IMAGE:figures/full_fig_p016_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Parametric study of the strain-limiting parameter [PITH_FULL_IMAGE:figures/full_fig_p017_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Influence of the parameter α on the crack-tip fields for fibers aligned with the x-axis (Case I), with fixed β = 0.02. Higher values of α lead to a marked amplification of stress, strain, and strain energy density, indicating sharper gradients and reduced strain-limiting effectiveness near the tip. Figure (8) illustrates the influence of the parameter α. In contrast to the behavior observed for β , increas… view at source ↗
Figure 9
Figure 9. Figure 9: Contour plots of the axial strain (left) and axial stress (right) for Case I ( [PITH_FULL_IMAGE:figures/full_fig_p018_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Displacement contours for Case I (α = 0.5, β = 0.05). Left: Horizontal displacement u1, showing mild variation. Right: Vertical displacement u2, displaying a sharp antisymmetric transition across the horizontal mid-line induced by the piecewise boundary loading. Figure (10) presents the contour plots of the horizontal displacement (u1), and vertical dis￾placement (u2). The horizontal displacement field va… view at source ↗
Figure 11
Figure 11. Figure 11: Elevated surface view of the vertical displacement [PITH_FULL_IMAGE:figures/full_fig_p020_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Influence of the parameter β on the crack-tip fields for fibers aligned with the y-axis (Case II). Similar to Case I, increasing β significantly suppresses the peak magnitudes of stress, strain, and energy density, reflecting the material’s enhanced resistance to crack opening. Figure (12) shows the effect of varying β when the fibers are aligned with the y-axis. In this configuration, the preferred direc… view at source ↗
Figure 13
Figure 13. Figure 13: Influence of the parameter α on the crack-tip fields for fibers aligned with the y-axis (Case II), with β = 0.02. Consistent with Case I, higher α values lead to amplified crack-tip quantities and a reduction in the effective strain-limiting behavior. Figure (13) presents the influence of α for the case of fibers aligned with the y-axis. For fixed β, higher values of α lead to a marked amplification of th… view at source ↗
Figure 14
Figure 14. Figure 14: Contour plots of axial strain (left) and axial stress (right) for Case II ( [PITH_FULL_IMAGE:figures/full_fig_p022_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Displacement contours for Case II (α = 0.5, β = 0.05). Left: Horizontal displacement u1. Right: Vertical displacement u2. The vertical displacement field is dominant and exhibits a clear sign reversal across the mid-height, driven by the boundary loading. Figure (15) displays the contour plots of the horizontal displacement (u1) and the vertical displacement (u2) for the case in which the fibers are align… view at source ↗
Figure 16
Figure 16. Figure 16: Elevated surface view of the vertical displacement [PITH_FULL_IMAGE:figures/full_fig_p024_16.png] view at source ↗
read the original abstract

Crack-tip fields within a transversely isotropic strain-limiting elastic body are investigated under the influence of piecewise linear slope boundary loads. The mechanical response is characterized via a nonlinear constitutive framework relating the Cauchy stress to the linearized strain, by which non-physical strain singularities at the crack tip are eliminated. The governing system is formulated as a quasi-linear elliptic boundary value problem in terms of the displacement field and is solved utilizing a continuous Galerkin finite element method coupled with a Picard linearization scheme. Boundary conditions are prescribed such that the vertical displacement varies piecewise linearly along the top and bottom edges, exhibiting opposite slopes on each half of the boundary. Numerical results are derived for two distinct fiber orientations. It is demonstrated that piecewise slope loads provide a flexible and realistic configuration for elucidating the interplay between strain-limiting behavior and crack-tip mechanics in transversely isotropic media.

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 investigates crack-tip fields in transversely isotropic strain-limiting elastic solids subjected to piecewise linear slope boundary loads. It employs a nonlinear constitutive relation between Cauchy stress and linearized strain to eliminate non-physical singularities at the crack tip, formulates the problem as a quasi-linear elliptic BVP in displacements, and solves it via continuous Galerkin FEM with Picard iteration. Numerical results are presented for two fiber orientations, with the claim that such loads provide a flexible configuration for studying the interplay between strain-limiting behavior and crack-tip mechanics.

Significance. If the numerical evidence holds, the work would offer a computational demonstration of bounded crack-tip strains in anisotropic strain-limiting materials under non-standard boundary conditions, potentially aiding modeling in fracture mechanics where classical linear elasticity predicts singularities.

major comments (2)
  1. [Numerical Results] The central claim that the nonlinear constitutive framework eliminates non-physical strain singularities rests on the numerical results, yet no mesh-convergence study or error estimates are reported for the crack-tip strains (or displacements) under the piecewise linear slope loads for either fiber orientation. Without such verification as h→0, it remains possible that the computed fields exhibit mesh-dependent behavior, undermining the demonstration.
  2. [Numerical Results] No validation against known analytical limits or benchmark problems (e.g., uniform tension or standard mode-I loading in the isotropic limit) is provided to confirm that the FEM-Picard scheme recovers expected behavior away from the crack tip or for the chosen constitutive parameters.
minor comments (2)
  1. [Introduction] The abstract and introduction should explicitly state the specific form of the nonlinear constitutive law (e.g., the function relating stress to strain) and the values of material parameters used in the computations.
  2. [Figures] Figure captions and axis labels should include units and clarify the fiber orientations (e.g., angle with respect to the crack plane) for reproducibility.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments, which highlight important aspects of numerical verification. We address each major comment below and will incorporate the suggested additions into the revised manuscript.

read point-by-point responses

  1. Referee: [Numerical Results] The central claim that the nonlinear constitutive framework eliminates non-physical strain singularities rests on the numerical results, yet no mesh-convergence study or error estimates are reported for the crack-tip strains (or displacements) under the piecewise linear slope loads for either fiber orientation. Without such verification as h→0, it remains possible that the computed fields exhibit mesh-dependent behavior, undermining the demonstration.

    Authors: We agree that mesh-convergence verification is necessary to support the numerical demonstration. In the revised version we will add a dedicated subsection presenting a systematic h-refinement study for both fiber orientations. This will include tabulated values and plots of the maximum crack-tip strain and selected displacement components versus mesh size (or degrees of freedom), together with observed convergence rates, confirming that the reported fields remain bounded and stable as h→0. revision: yes

  2. Referee: [Numerical Results] No validation against known analytical limits or benchmark problems (e.g., uniform tension or standard mode-I loading in the isotropic limit) is provided to confirm that the FEM-Picard scheme recovers expected behavior away from the crack tip or for the chosen constitutive parameters.

    Authors: We accept that explicit validation strengthens the credibility of the computational framework. The revised manuscript will include a new validation subsection that (i) recovers the isotropic linear-elastic solution under uniform far-field tension by setting the anisotropy parameters to zero and (ii) compares the computed fields under standard mode-I loading against the known asymptotic behavior away from the tip. These tests will be performed with the same Picard iteration and constitutive parameters (in the isotropic limit) to demonstrate that the scheme reproduces expected linear-elastic behavior outside the strain-limiting zone. revision: yes

Circularity Check

0 steps flagged

No circularity: numerical solution of constitutive model is independent of reported fields

full rationale

The paper formulates a quasi-linear elliptic BVP from a nonlinear stress-strain relation chosen to bound strains at the crack tip, then discretizes it with continuous Galerkin FEM plus Picard iteration. No equation reduces a reported field to a fitted parameter or self-referential definition; the piecewise-linear slope BCs and fiber orientations are external inputs, and the boundedness result follows from the model properties rather than from the numerical output itself. No self-citation chains or ansatz smuggling appear in the derivation steps described.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Only the abstract is available; the ledger is therefore limited to the constitutive assumption stated there.

axioms (1)
  • domain assumption A nonlinear constitutive relation between Cauchy stress and linearized strain eliminates non-physical strain singularities at the crack tip.
    Explicitly invoked in the abstract as the mechanical characterization that removes singularities.

pith-pipeline@v0.9.0 · 5453 in / 1204 out tokens · 64697 ms · 2026-05-16T08:34:13.487148+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

23 extracted references · 23 canonical work pages

  1. [1]

    T. L. Anderson.Fracture Mechanics: Fundamentals and Applications. CRC Press, 2005

    work page 2005

  2. [2]

    The deal.II library, version 9.5.Journal of Numerical Mathematics, 31(3):231–246, aug 2023

    Daniel Arndt, Wolfgang Bangerth, Maximilian Bergbauer, Marco Feder, Marc Fehling, Jo- hannes Heinz, Timo Heister, Luca Heltai, Martin Kronbichler, Matthias Maier, Peter Munch, Jean-Paul Pelteret, Bruno Turcksin, David Wells, and Stefano Zampini. The deal.II library, version 9.5.Journal of Numerical Mathematics, 31(3):231–246, aug 2023. Preprint available ...

    work page 2023

  3. [3]

    The deal.II finite element library: Design, features, and insights.Computers & Mathematics with Appli- cations, 81:407–422, 2021

    Daniel Arndt, Wolfgang Bangerth, Denis Davydov, Timo Heister, Luca Heltai, Martin Kro- nbichler, Matthias Maier, Jean-Paul Pelteret, Bruno Turcksin, and David Wells. The deal.II finite element library: Design, features, and insights.Computers & Mathematics with Appli- cations, 81:407–422, 2021

    work page 2021

  4. [4]

    L. Beck, M. Bulíček, J. Málek, and E. Süli. On the existence of integrable solutions to nonlinear elliptic systems and variational problems with linear growth.Archive for Rational Mechanics and Analysis, 225(2):717–769, 2017. 25

    work page 2017

  5. [5]

    Bulíček, J

    M. Bulíček, J. Málek, K. R. Rajagopal, and E. Süli. On elastic solids with limiting small strain: modelling and analysis.EMS Surveys in Mathematical Sciences, 1(2):283–332, 2014

    work page 2014

  6. [6]

    Bustamante

    R. Bustamante. Some topics on a new class of elastic bodies.Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 465(2105):1377–1392, 2009

    work page 2009

  7. [7]

    Gou and S

    K. Gou and S. M. Mallikarjunaiah. Finite element study of v-shaped crack-tip fields in a three-dimensional nonlinear strain-limiting elastic body.Mathematics and Mechanics of Solids, 28(10):2155–2176, 2023

    work page 2023

  8. [8]

    K. Gou, S. M. Mallikarjunaiah, K. R. Rajagopal, and J. R. Walton. Modeling fracture in the context of a strain-limiting theory of elasticity: A single plane-strain crack.International Journal of Engineering Science, 88:73–82, 2015

    work page 2015

  9. [9]

    M. E. Gurtin.An Introduction to Continuum Mechanics. Academic Press, 1981

    work page 1981

  10. [10]

    C. E. Inglis. Stresses in a plate due to the presence of cracks and sharp corners.Transactions of the Royal Institute of Naval Architectes, 60:219–241, 1913

    work page 1913

  11. [11]

    H. Itou, V. A. Kovtunenko, and K. R. Rajagopal. On the states of stress and strain adjacent to a crack in a strain-limiting viscoelastic body.Mathematics and Mechanics of Solids, 23(3):433– 444, 2018

    work page 2018

  12. [12]

    Kulvait, J

    V. Kulvait, J. Málek, and K. R. Rajagopal. Anti-plane stress state of a plate with a v-notch for a new class of elastic solids.International Journal of Fracture, 179(1-2):59–73, 2013

    work page 2013

  13. [13]

    S. Lee, H. C. Yoon, and S. M. Mallikarjunaiah. Finite element simulation of quasi-static tensile fracture in nonlinear strain-limiting solids with the phase-field approach.Journal of Computational and Applied Mathematics, 399:113715, 2022

    work page 2022

  14. [14]

    A. E. H. Love.A treatise on the mathematical theory of elasticity. Cambridge university press, 2013

    work page 2013

  15. [15]

    S. M. Mallikarjunaiah.On Two Theories for Brittle Fracture: Modeling and Direct Numerical Simulations. PhD thesis, Texas A&M University, 2015

    work page 2015

  16. [16]

    S. M. Mallikarjunaiah and J. R. Walton. On the direct numerical simulation of plane-strain fractureinaclassofstrain-limitinganisotropicelasticbodies.International Journal of Fracture, 192(2):217–232, Apr 2015

    work page 2015

  17. [17]

    K. R. Rajagopal. On implicit constitutive theories.Applications of Mathematics, 48(4):279– 319, 2003

    work page 2003

  18. [18]

    K. R. Rajagopal. The elasticity of elasticity.Zeitschrift fÃŒr Angewandte Mathematik und Physik (ZAMP), 58(2):309–317, 2007

    work page 2007

  19. [19]

    K. R. Rajagopal. Non-linear elastic bodies exhibiting limiting small strain.Mathematics and Mechanics of Solids, 16(1):122–139, 2011

    work page 2011

  20. [20]

    K. R. Rajagopal and J. R. Walton. Modeling fracture in the context of a strain-limiting theory of elasticity: a single anti-plane shear crack.International journal of fracture, 169(1):39–48, 2011. 26

    work page 2011

  21. [21]

    H. C. Yoon, S. Lee, and S. M. Mallikarjunaiah. Quasi-static anti-plane shear crack propagation in nonlinear strain-limiting elastic solids using phase-field approach.International Journal of Fracture, 227(2):153–172, 2021

    work page 2021

  22. [22]

    H. C. Yoon and S. M. Mallikarjunaiah. A finite-element discretization of some boundary value problems for nonlinear strain-limiting elastic bodies.Mathematics and Mechanics of Solids, 27(2):281–307, 2022

    work page 2022

  23. [23]

    H. C. Yoon, K. K. Vasudeva, and S. M. Mallikarjunaiah. Finite element model for a coupled thermo-mechanical system in nonlinear strain-limiting thermoelastic body.Communications in Nonlinear Science and Numerical Simulation, 108:106262, 2022. 27

    work page 2022