Pith. sign in

REVIEW 3 major objections 5 minor 26 references

Moving Collinear Cracks in a Prestressed Dry Sandy Medium Fracture Response under Traveling Punch Loads

T0 review · 3 major / 5 minor · reviewed 2026-08-01 · deepseek-v4-flash

Pith's one-line read The paper derives closed-form formulas for the Mode-I stress intensity factors and crack opening displacement of two collinear cracks moving in a pre-stressed dry sandy strip under a moving punch, with explicit finite-thickness corrections.

desk verdict The paper offers a new analytical configuration but solves a boundary-value problem different from the one it states, and its 'closed-form' SIFs contain unevaluated constants. read the letter →

arxiv 2607.17155 v1 pith:GNLSJLNZ submitted 2026-07-19 math-ph cond-mat.mtrl-scimath.MP

classification math-phcond-mat.mtrl-scimath.MP MSC 74R1045E05
keywords movingcollinearcracksdrysandymediuminitialstresspunchloadintensityfactorcrackopeningdisplacementsingularintegralequationsfiniteHilberttransform
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper sets out to prove that the fracture response of two collinear Griffith cracks (idealized line cracks) moving steadily in a pre-stressed dry sandy strip under a moving punch can be written in closed form. The target results are explicit formulas for the crack density functions, the Mode-I stress intensity factors at the inner and outer crack tips, and the crack opening displacement, each carried to the leading correction in the inverse square of the strip thickness. The model simultaneously includes crack speed, initial stress, sandiness, strip thickness, crack geometry, and punch position, so the formulas expose which combinations of parameters drive crack-tip intensification. A careful reader would care because this is the geometry of pavement, railway foundation, and underground excavation failures, and because an analytical result for this parameter combination has not previously been available. The route is steady-state transformation, Fourier transforms, coupled Cauchy-type singular integral equations, a thick-strip kernel reduction, and inversion by the finite Hilbert transform.

What carries the argument

The mechanism carrying the argument is the asymptotic kernel reduction of Section 4.1. The exact transform-domain kernel multipliers N₁₁, N₁₂, N₂₁ contain 1/sinh(r_j ζh) factors; for a thick strip the paper reduces them, using the series identities for 1/k² and 1/(2k+1)², to simple linear forms proportional to P s/h². That reduction changes two coupled Cauchy-type singular integral equations into four algebraic conditions, which are inverted in closed form by the finite Hilbert transform, an integral inversion technique for Cauchy principal-value equations. All physical parameters enter through the effective elastic coefficients D₁–D₅, the sandiness-dependent shear coefficient S = µ/χ, the i

What would settle it

Solve the coupled integral equations with the exact kernels numerically for the same parameters as Section 5 and compare the resulting stress intensity factors against the closed-form formulas; if the difference is not O(h⁻²) or does not vanish as h/e increases, the asymptotic reduction is refuted. A cheaper check is to expand the kernel multipliers around ζ = 0 and confirm that the 1/ζ poles cancel exactly, leaving the claimed coefficient P.

Watch

Extended reading notes

Core claim

The central claim is that the coupled problem can be solved, not just reduced, after a thick-strip simplification: for large h the regular kernels collapse to linear forms proportional to s/h² with one common material-loading constant P, and the density functions then split into an unbounded-strip part plus an h⁻² correction. Substitution into the near-tip limits gives the Mode-I stress intensity factor at the inner tip, K_I^(c) = (π/(2Q₁)) (1/(√c√(e²−c²))) (D₁ + D₃/h²), and an analogous closed-form expression at the outer tip, with punch loading and crack interaction entering through constants H₁ and H₂. The crack opening displacement follows as an elliptic-integral expression with the same

Load-bearing premise

The load-bearing premise is that the exact thick-strip kernels reduce to the simple O(h⁻²) linear forms with the single constant P using only two series identities; the paper does not exhibit the cancellation of the kernel poles at ζ = 0 or verify the reduction numerically, so if that reduction is incomplete the finite-thickness corrections in the stress-intensity-factor formulas are wrong.

Editorial extensions

If this is right

  • If the paper is correct, the Mode-I stress intensity factors at both crack tips of a thick sandy strip can be computed from closed-form expressions instead of solving singular integral equations.
  • The formulas quantify a specific asymmetry: the outer crack tip carries a larger stress intensity factor than the inner tip because the moving concentrated load sits closer to it.
  • The results show a common dynamic pattern: crack-tip intensity remains nearly flat until the crack speed approaches the effective shear-wave speed, then rises sharply across all material and geometric parameters.
  • The limiting cases V = 0, zero initial stress, and the zero-sandiness/punch-free configuration reproduce earlier simpler models, providing internal consistency checks.
  • The finite-thickness correction couples punch loading and inner-tip intensity through the constants H₁ and H₂, so the model captures the interaction between boundary loading and crack interaction at the first nontrivial order.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Beyond the paper: the thick-strip kernel reduction could be tested numerically by discretizing the exact integral equations; a mismatch that is not O(h⁻²) would map out where the asymptotic window ends.
  • Beyond the paper: because the sandiness and initial-stress effects enter only through effective coefficients, the closed-form structure should survive for any isotropic medium with a modified shear modulus; the geometrical form of the formulas is set by the crack configuration, not by the material model.
  • Beyond the paper: the same reduction-and-Hilbert-transform route could extend to periodic collinear crack arrays, where the two-crack ligament is replaced by a lattice spacing and the Hilbert-transform weights would need only modest changes.
  • Beyond the paper: the predicted non-monotonic inner-tip response as a function of crack length, with a minimum near c/e ≈ 0.10–0.15, is a specific signature that targeted simulations or experiments on crack interaction could confirm or refute.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper treats two steadily moving collinear Mode-I Griffith cracks in an initially stressed dry sandy strip of finite thickness, with concentrated crack-face loads and a moving punch pressure on the outer surfaces. The authors use a moving coordinate frame, Fourier transforms, and a large-thickness asymptotic reduction to convert the mixed boundary-value problem into coupled Cauchy singular integral equations (58)-(59), which they solve analytically with the finite Hilbert transform. The output consists of the crack-density functions (81)-(82), closed-form Mode-I stress intensity factors at the inner and outer crack tips (87)-(88), and a crack-opening-displacement expression (89), all to O(h^{-2}), together with parametric studies of the effects of crack speed, sandiness, initial stress, strip thickness, geometry, and loading position. The derivation is self-contained and the algebra connecting the integral equations to the asymptotic solution is internally consistent; the main weaknesses are in the statement of the boundary conditions and in the unproved large-thickness kernel reduction.

Significance. If the results are correct, this is a useful analytical benchmark for a fairly complicated configuration: two interacting moving cracks in a pre-stressed, sand-modified elastic strip under combined concentrated and distributed moving loads. The paper provides explicit, closed-form expressions rather than purely numerical results, and it contains several limiting-case checks, including a claimed reduction to a previously published monoclinic-strip problem. These are genuine strengths. However, the significance is conditional on resolving two load-bearing issues: the boundary-value problem actually solved differs from the one stated, and the asymptotic kernel formulas that generate the finite-thickness/punch corrections are asserted without derivation or numerical verification.

major comments (3)
  1. [§2.4 and §3.3, Eqs. (19), (23), (49)-(54)] The stated boundary-value problem and the problem actually solved are not the same. Eq. (19) requires u1(x1,h)=0 for |x1|<c and |x1|>e, and Eq. (23) similarly requires u1(x1,0)=0 outside the cracks. The transform solution never imposes these conditions: the amplitude relations (49)-(54) enforce σ12=0 at x2=0 and x2=h and use the representations (50)-(51) to enforce u2=0 on the relevant intervals, but U(ζ,h) and U(ζ,0) remain unconstrained. Consequently the SIE system (58)-(59) and all subsequent SIF/COD formulas correspond to the smooth-punch/Mode-I problem with u2=0 and σ12=0 prescribed, not to the clamped conditions in (19) and (23). In fact, for a Mode-I crack the correct mid-plane condition is u2=0 outside the cracks, not u1=0; so Eq. (23) is not a Mode-I condition either. The manuscript must either delete the u1=0 conditions from the problem statement or extend the transform solutio
  2. [§4.1, Eqs. (70)-(72)] The large-thickness asymptotic reduction is the only source of the finite-thickness and punch-coupling corrections in the final SIFs, yet it is asserted without derivation. The exact kernel multipliers in (63)-(67) contain F1, F2, F3, F4 built from 1/sinh(r_j ζh), which near ζ=0 behave as 1/(r_j ζh), giving 1/ζ poles. For the leading terms to be O(h^{-2}) and linear in s as claimed in (70)-(72), these poles must cancel after combining Π_j, Λ_j, r_j, and the finite parts must collapse to the single constant P in (73). None of this cancellation is exhibited, and no numerical check of (70)-(72) is reported. Because H1 and H2 in (83)-(84) and the h^{-2} contributions to K_I^(c) and K_I^(e) in (87)-(88) depend directly on these forms, this is a load-bearing gap. Please provide the derivation or an appendix, or verify the asymptotic kernels numerically against the exact integral representation
  3. [§4.5.3, Eq. (106)] The 'sandiness-free' limiting case is mis-specified. The constitutive relation (4) contains the factor 2µ/χ with χ>1, and §2.2 explicitly states that χ=1 recovers the classical isotropic response. Setting χ=0 in Eq. (106), and subsequently in (107)-(111), makes the effective shear modulus singular and is not a well-defined limit. If the intended consistency check is a sandiness-free medium, the parameter should be χ=1 throughout this subsection; otherwise the claimed reduction to the corresponding case in [9] is not established. This also propagates into the formulas for H_1^(b) and H_2^(b) in (114)-(115).
minor comments (5)
  1. [§1.6] The second objective states that the solution uses the 'Schmidt method', but the paper actually uses the finite Hilbert transform technique. Please correct the wording.
  2. [§4.2, after Eq. (82)] The sentence defining the coefficients H1 and H2 ('The coefficients H1 and H2 ... are given by') appears twice in succession. Please remove the duplication.
  3. [§2.3 and §4.2] The symbol D_j is used for two different sets of quantities: the effective coefficients in Eq. (31) and the integration constants in Eqs. (81)-(82). This is confusing and should be changed, e.g., to A_j for the integration constants.
  4. [§5, Table 1] The quantity Q_s/Q_c is repeatedly called the 'tangential loading ratio', but Q_s is the normal punch pressure in Eq. (80) and Q_c is the concentrated normal crack-face load. A neutral term such as 'punch-to-concentrated load ratio' would be more accurate.
  5. [§4.1, Eq. (73)] The constant P is defined in (73) as a combination of the Π_j, Λ_j, r_j, and C, but the text does not explain why the same P appears in all three kernel asymptotics (70)-(72). A sentence indicating the algebraic origin of this common factor would help.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the SIF/COD results are derived from the stated boundary-value problem, not fitted or defined in terms of the targets.

full rationale

The derivation is self-contained. The constitutive law (4) is taken from external reference [21]; the moving-frame field equations (15)-(16) follow from it. The transform-domain amplitudes are reduced using the stated boundary conditions, and the auxiliary density functions g1,g2 are determined from the coupled Cauchy singular integral equations (58)-(59) together with the consistency conditions (68), not from the target SIFs. The large-thickness kernel reduction (70)-(72) is an asymptotic approximation asserted after the series identities (69); although the pole cancellations are not exhibited and the reduction is not numerically verified, this is a completeness/correctness risk, not a circular step, because the claimed kernel forms are not defined in terms of the final stress-intensity-factor outputs. The closed-form SIFs (87)-(88) and COD (89) are obtained by solving those integral equations via the finite Hilbert transform; the constants D_j are integration constants fixed by the supplementary conditions (68), not fitted parameters. The only self-citation is the consistency check in Sec. 4.5.3 claiming reduction to the authors' prior model [9]; that comparison is not load-bearing in the derivation, which does not assume the SIF/COD values from [9]. A separate concern that the stated upper-surface condition u1=0 (Eq. 19) is not enforced in Sec. 3.3 is a boundary-condition correctness issue, not a circularity.

Assumptions & free parameters 3 free parameters · 7 assumptions · 0 invented entities

The model has no data-fitted parameters: every input is either a standard elastic constant (λ, μ, ρ from [21]) or a hand-chosen study value. The ledger is dominated by domain assumptions inherited from the literature, plus one assumption specific to this paper — the asserted O(h⁻²) kernel reduction (Eqs. 70-72) — which is load-bearing and unverified. No new entities (particles, forces, conserved quantities) are introduced; the 'sandiness parameter' χ is an inherited scalar rescaling of the shear modulus.

free parameters (3)
  • Sandiness parameter χ = 1.30 (reference value, from [21])
    Rescales the shear modulus in Eq. (4). Not fitted in this paper but inherited from the cited constitutive model [21]; varied parametrically in §5.1.1.
  • Initial stresses σ⁰₁₁, σ⁰₂₂ = 1.0×10⁹ Pa each (reference values)
    Hand-chosen reference values (Table 1), not from field data; varied parametrically in §5.1.2.
  • Loading ratio Qs/Qc and load position x₀/e = Qs/Qc = 2.1, x₀/e = 0.96 (Table 1); x₀/e = 0.5 (Table 2)
    Hand-chosen reference loading parameters for the numerical study; not fitted to any data.
assumptions (7)
  • domain assumption Steady-state assumption: cracks, crack-face loads, and surface punch translate together at constant speed V, so ∂/∂t → −V∂/∂x₁ (Eq. 14)
    Converts the transient problem into a stationary BVP in the moving frame; standard Yoffe-type reduction but restricts the model to uniform, steady long-time propagation.
  • domain assumption Small-on-large incremental equations with initial stress: σij,j + (u_i,k Σ⁰_kj),j = ρ ü_i (Eq. 11)
    Linearized incremental motion about a uniform biaxial pre-stress; assumes the initial stress enters only through the incremental equilibrium. Standard but strong.
  • domain assumption Sandiness constitutive law (Eq. 4): σij = λ εkk δij + 2μ/χ εij, χ > 1
    Adopted from [21]; algebraically identical to isotropic elasticity with effective shear modulus μ/χ, so the granular-physics content is a single scalar rescaling.
  • ad hoc to paper Large-thickness kernel asymptotics (Eqs. 70-72): Kij(x,s) = O(h⁻²) linear-in-s kernels with a common constant P
    The load-bearing reduction to the coupled O(h⁻²) system (77)-(78); asserted with only the series identities (69), and the cancellation of the 1/ζ poles of F₁-F₄ is not demonstrated.
  • domain assumption Distinct real characteristic roots r₁ ≠ r₂ (Eq. 40)
    The modal basis (41)-(42) requires distinct roots. At V = 0 the characteristic equation is a perfect square with r₁ = r₂ = 1, so the claimed stationary limit (§4.5.1) is not a genuine limit of the modal solution.
  • standard math Finite Hilbert transform inversion technique (as in [23])
    Used to invert the Cauchy-type equations (75)-(78); cited, not re-derived.
  • domain assumption Plane strain, infinitesimal deformation, no body forces (Section 2.1)
    Standard restrictions for this class of fracture models.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Moving Collinear Cracks in a Prestressed Dry Sandy Medium Fracture Response under Traveling Punch Loads." pith.science (2026). https://pith.science/paper/GNLSJLNZ

@misc{pith2026260717155,
  author       = {Pith},
  title        = {Pith review of: Moving Collinear Cracks in a Prestressed Dry Sandy Medium Fracture Response under Traveling Punch Loads},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/GNLSJLNZ}},
  note         = {Machine review of arXiv:2607.17155}
}
read the original abstract

The dynamic fracture behaviour of two moving collinear Griffith cracks in an initially stressed dry sandy medium subjected to concentrated crack-face loading and moving punch pressure is investigated. A moving coordinate system transforms the transient problem into a steady state formulation, while the effects of initial stress and sandiness are incorporated into the governing equations. Fourier integral transforms are employed to obtain the characteristic equation and the transformed traction displacement relations. The crack face, symmetry, and outer surface conditions reduce the problem to coupled Cauchy type singular integral equations. For a sufficiently thick strip, an asymptotic kernel reduction is developed, and the resulting equations are solved analytically using the finite Hilbert transform. Closed form expressions are derived for the crack density functions, Mode I stress intensity factors at the inner and outer crack tips, and the crack opening displacement. Several limiting cases are recovered from the general formulation, providing analytical consistency checks. The results demonstrate the combined influence of crack speed, crack geometry, initial stress, sandiness, and moving punch loading on crack tip intensification and crack opening. The proposed analytical framework provides useful insights for fracture assessment and the design of transportation infrastructure, geotechnical systems, underground excavations, and other engineering structures involving dry sandy media subjected to moving loads.

Figures

Figures reproduced from arXiv: 2607.17155 by the authors.

Figure 1
Figure 1. Schematic configuration of the problem showing two moving collinear [PITH_FULL_IMAGE:figures/full_fig_p008_1.png] view at source ↗
Figure 2
Figure 2. Variation of the normalized Mode-I stress intensity factor with crack [PITH_FULL_IMAGE:figures/full_fig_p024_2.png] view at source ↗
Figure 3
Figure 3. Variation of the normalized Mode-I stress intensity factors with the [PITH_FULL_IMAGE:figures/full_fig_p025_3.png] view at source ↗
Figures from the paper (18 more)
Figure 4
Figure 4. Figure 4: Variation of the normalized Mode-I stress intensity factors at the outer [PITH_FULL_IMAGE:figures/full_fig_p026_4.png]
Figure 5
Figure 5. Figure 5: Variation of the normalized Mode-I stress intensity factors at the outer [PITH_FULL_IMAGE:figures/full_fig_p026_5.png]
Figure 6
Figure 6. Figure 6: Variation of the normalized Mode-I stress intensity factors at the outer [PITH_FULL_IMAGE:figures/full_fig_p027_6.png]
Figure 7
Figure 7. Figure 7: Contour plots of the normalized Mode-I stress intensity factors as func [PITH_FULL_IMAGE:figures/full_fig_p029_7.png]
Figure 8
Figure 8. Figure 8: Contour plots of the normalized Mode-I stress intensity factors as func [PITH_FULL_IMAGE:figures/full_fig_p030_8.png]
Figure 9
Figure 9. Figure 9: Contour plots of the normalized Mode-I stress intensity factors as func [PITH_FULL_IMAGE:figures/full_fig_p031_9.png]
Figure 10
Figure 10. Figure 10: Contour plots of the normalized Mode-I stress intensity factors as [PITH_FULL_IMAGE:figures/full_fig_p031_10.png]
Figure 11
Figure 11. Figure 11: Contour plots of the normalized Mode-I stress intensity factors as [PITH_FULL_IMAGE:figures/full_fig_p032_11.png]
Figure 12
Figure 12. Figure 12: Variation of the normalized Mode-I stress intensity factors with the [PITH_FULL_IMAGE:figures/full_fig_p035_12.png]
Figure 13
Figure 13. Figure 13: Variation of the normalized Mode-I stress intensity factors with the [PITH_FULL_IMAGE:figures/full_fig_p035_13.png]
Figure 14
Figure 14. Figure 14: Variation of the normalized Mode-I stress intensity factors with the [PITH_FULL_IMAGE:figures/full_fig_p036_14.png]
Figure 15
Figure 15. Figure 15: Variation of the normalized Mode-I stress intensity factors with the [PITH_FULL_IMAGE:figures/full_fig_p036_15.png]
Figure 16
Figure 16. Figure 16: Variation of the normalized Mode-I stress intensity factors with the [PITH_FULL_IMAGE:figures/full_fig_p037_16.png]
Figure 17
Figure 17. Figure 17: Contour representation of the normalized Mode-I stress intensity fac [PITH_FULL_IMAGE:figures/full_fig_p038_17.png]
Figure 18
Figure 18. Figure 18: Contour representation of the normalized Mode-I stress intensity fac [PITH_FULL_IMAGE:figures/full_fig_p039_18.png]
Figure 19
Figure 19. Figure 19: Contour representation of the normalized Mode-I stress intensity fac [PITH_FULL_IMAGE:figures/full_fig_p040_19.png]
Figure 20
Figure 20. Figure 20: Contour representation of the normalized Mode-I stress intensity fac [PITH_FULL_IMAGE:figures/full_fig_p040_20.png]
Figure 21
Figure 21. Figure 21: Contour representation of the normalized Mode-I stress intensity fac [PITH_FULL_IMAGE:figures/full_fig_p041_21.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

26 extracted references

  1. [9]

    Fatigue & Fracture of Engineering Materials & Structures 49(1), 284–304 (2026)

    Diksha, ., Chaudhary, S., Othmani, C., Kumar Sharma, P.: Dynamic fracture analysis of moving mode i collinear cracks in monoclinic crystalline strip: An analytical approach using hilbert transform. Fatigue & Fracture of Engineering Materials & Structures 49(1), 284–304 (2026)

  2. [1]

    Physical Review E60(2), 2366 (1999)

    Adda-Bedia, M., Arias, R., Amar, M.B., Lund, F.: Generalized griffith criterion for dynamic fracture and the stability of crack motion at high velocities. Physical Review E60(2), 2366 (1999)

  3. [2]

    Applied Mathematics and Mechanics38(9), 1313–1328 (2017)

    Alam, P., Kundu, S., Gupta, S.: Dispersion and attenuation of torsional wave in a viscoelastic layer bonded between a layer and a half-space of dry sandy media. Applied Mathematics and Mechanics38(9), 1313–1328 (2017)

  4. [3]

    In: ARMA US Rock Mechanics/Geomechanics Symposium, p

    Amir Hosseini, M., Tahmasebi, P.: A new approach to model irregular particles breakage with deploying image analysis techniques on discrete element method. In: ARMA US Rock Mechanics/Geomechanics Symposium, p. D022S023R009. ARMA (2024)

  5. [4]

    Powder Technology162(3), 208– 229 (2006)

    Campbell, C.S.: Granular material flows–an overview. Powder Technology162(3), 208– 229 (2006)

  6. [5]

    i: Theory

    Chang, C.S., Chang, Y., Kabir, M.G.: Micromechanics modeling for stress-strain behav- ior of granular soils. i: Theory. Journal of geotechnical engineering118(12), 1959–1974 (1992)

  7. [6]

    Applied Mathematical Modelling140, 115916 (2025) 42 Diksha et al

    Chaudhary, S., Sharma, P.K., et al.: Crack dynamics in rotating, initially stressed ma- terial strip: A mathematical approach. Applied Mathematical Modelling140, 115916 (2025) 42 Diksha et al

  8. [7]

    International Journal of Solids and Structures43(25-26), 7880–7890 (2006)

    Das, S.: Interaction of moving interface collinear griffith cracks under antiplane shear. International Journal of Solids and Structures43(25-26), 7880–7890 (2006)

Show all 26 references
  1. [8]

    Mathematical Methods in the Applied Sciences (2026)

    Das, S., Murshid AK, M., Mondal, S.: Magnetoelastic plane wave propagation in an initially stressed anisotropic strip through collinear mode i griffith cracks with punch loading. Mathematical Methods in the Applied Sciences (2026)

  2. [10]

    International Journal for Numerical and Analytical Methods in Geomechanics40(1), 105–116 (2016)

    Druckrey, A.M., Alshibli, K.A.: 3d finite element modeling of sand particle fracture based on in situ x-ray synchrotron imaging. International Journal for Numerical and Analytical Methods in Geomechanics40(1), 105–116 (2016)

  3. [11]

    Proceedings of the Royal So- ciety A: Mathematical, Physical and Engineering Sciences463(2087), 3021–3035 (2007)

    Einav, I.: Fracture propagation in brittle granular matter. Proceedings of the Royal So- ciety A: Mathematical, Physical and Engineering Sciences463(2087), 3021–3035 (2007)

  4. [12]

    Cambridge university press (1998)

    Freund, L.B., Freud, L.: Dynamic fracture mechanics. Cambridge university press (1998)

  5. [13]

    European Journal of Soil Science46(4), 591–599 (1995)

    Hallett, P., Dexter, A., Seville, J.: The application of fracture mechanics to crack prop- agation in dry soil. European Journal of Soil Science46(4), 591–599 (1995)

  6. [14]

    Tumu Gongcheng Xuebao/China Civil Engineering Journal49(7), 9–35 (2016)

    Huang, M., Yao, Y., Yin, Z., Liu, E., Lei, H.: An overview on elementary mechani- cal behaviors, constitutive modeling and failure criterion of soils. Tumu Gongcheng Xuebao/China Civil Engineering Journal49(7), 9–35 (2016)

  7. [15]

    CRC press (2020)

    Iwashita, K., Oda, M.: Mechanics of granular materials: an introduction. CRC press (2020)

  8. [16]

    Acta Mechanica227(12), 3671–3686 (2016)

    Li, Y.: Multiple collinear griffith cracks in a one-dimensional hexagonal quasicrystalline layer. Acta Mechanica227(12), 3671–3686 (2016)

  9. [17]

    Theoretical and applied fracture mechanics 56(1), 49–57 (2011)

    Monfared, M., Ayatollahi, M.: Dynamic stress intensity factors of multiple cracks in a functionally graded orthotropic half-plane. Theoretical and applied fracture mechanics 56(1), 49–57 (2011)

  10. [18]

    Journal of Geophysical Research: Solid Earth125(1), e2019JB018604 (2020)

    Muqtadir, A., Al-Dughaimi, S., Dvorkin, J.: Deformation of granular aggregates: Static and dynamic bulk moduli. Journal of Geophysical Research: Solid Earth125(1), e2019JB018604 (2020)

  11. [19]

    In: Crack Dynamics in Metallic Materials, pp

    Nilsson, F.: Dynamic fracture theory. In: Crack Dynamics in Metallic Materials, pp. 1–68. Springer (1990)

  12. [20]

    Journal of Vibration Engineering & Technologies14(4), 199 (2026)

    Rahul, A.K., et al.: Love wave propagation across slip interfaces in a viscoelastic or- thotropic layer bounded by dry sandy and transversely isotropic media. Journal of Vibration Engineering & Technologies14(4), 199 (2026)

  13. [21]

    Waves in Random and Complex Media 35(6), 10803–10820 (2025)

    Singh, A.K., Singh, A.K.: Mathematical study on the propagation of griffith crack in a dry sandy strip subjected to punch pressure. Waves in Random and Complex Media 35(6), 10803–10820 (2025)

  14. [22]

    Journal of Thermal Stresses44(4), 433–455 (2020)

    Singh, R., Das, S.: Investigation of interactions among collinear griffith cracks situated in a functionally graded medium under thermo-mechanical loading. Journal of Thermal Stresses44(4), 433–455 (2020)

  15. [23]

    Proceedings of the Royal Society of Edinburgh Section A: Mathematics68(4), 309–321 (1970)

    Srivastava, K., Lowengrub, M.: Xx.—finite hilbert transform technique for triple integral equations with trigonometric kernels. Proceedings of the Royal Society of Edinburgh Section A: Mathematics68(4), 309–321 (1970)

  16. [24]

    Progress in materials science138, 101157 (2023)

    Tahmasebi, P.: A state-of-the-art review of experimental and computational studies of granular materials: Properties, advances, challenges, and future directions. Progress in materials science138, 101157 (2023)

  17. [25]

    Advances in Civil Engineering2018(1), 5459307 (2018)

    Wang, C., Zhao, Y., Zhao, Y., Wan, W.: Study on the interaction of collinear cracks and wing cracks and cracking behavior of rock under uniaxial compression. Advances in Civil Engineering2018(1), 5459307 (2018)

  18. [26]

    the moving griffith crack

    Yoffe, E.H.: Lxxv. the moving griffith crack. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science42(330), 739–750 (1951)

Pith tools

Reviewed August 1, 2026 · model on record in the stance chip above.