Pith. sign in

REVIEW 4 major objections 6 minor 1 cited by

Biaxial characterization of soft elastomers: experiments and data-adaptive configurational forces for fracture

T0 review · 4 major / 6 minor · reviewed 2026-08-07 · deepseek-v4-flash

Pith's one-line read The paper claims that the magnitude of the total crack-tip configurational force at crack onset gives a finite-strain fracture-toughness criterion for soft elastomers under equi-biaxial loading.

desk verdict Solid biaxial dataset and a coherent pipeline, but the headline fracture criterion is calibrated to the same onset data and is not yet a prediction. read the letter →

arxiv 2505.20244 v1 pith:3WUVJKT4 submitted 2025-05-26 cond-mat.soft cond-mat.mtrl-sci

classification cond-mat.softcond-mat.mtrl-sci MSC 74B2074R1074S05
keywords configurationalforcesfracturetoughnesssoftelastomersequi-biaxialloadinghyperelasticityJ-integraldata-drivenconstitutivemodelingdigitalimagecorrelation
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 argues that fracture onset in soft elastomers under equi-biaxial stretching can be predicted by a single computed number: the total configurational force at the crack tip. The authors run equi-biaxial experiments to rupture on five elastomers (Sylgard, three Elastosil blends, and VHB tape) and observe a wide behavioral spectrum, from brittle Sylgard to VHB with crack-tip strains above 150%. They combine these experiments with a data-adaptive B-spline hyperelastic model and a post-processing Configurational Force Method that estimates the J-integral without contour integration. At the experimentally observed crack onset, the critical configurational force values (0.0386 to 0.4228 N/m) rank the materials consistently with their measured total work of fracture. If the criterion holds, the critical configurational force offers a computationally efficient, finite-strain fracture-toughness parameter for soft solids.

What carries the argument

The framework rests on two components. First, a data-adaptive hyperelastic strain energy density is represented by additive cubic B-spline interpolation over the isochoric invariants \bar I_1 and \bar I_2, with interpolation points and values identified by nonlinear least-squares optimization against reaction forces and DIC displacement fields from the biaxial tests; the optimizer consistently set the \bar I_2 contribution to zero, leaving \bar I_1-only energy functions. Second, the Configurational Force Method is implemented as a post-processing algorithm: after solving the forward boundary-value problem with the identified energy, the Eshelby stress tensor is computed and integrated against the FE basis-function gradients to produce nodal configurational forces, and the nodal forces inside a 0.1 mm radius cylinder around the crack tip are summed to give the total crack-tip configurational force, which the paper equates with the J-integral.

What would settle it

Refine the finite-element mesh near the crack tip and vary the summation radius around 0.1 mm while recomputing the total configurational force at crack onset; if the reported critical values (0.0386 to 0.4228 N/m) change substantially, the proposed fracture criterion is not reliable. An independent check would evaluate the contour J-integral on a Pacman-shaped domain in the same virtual experiments and test whether it equals the summed nodal configurational force.

Watch

Extended reading notes

Core claim

The central claim is that the magnitude of the total crack-tip configurational force at crack onset serves as a measure of the material's fracture toughness, leading to a configurational-force-based fracture criterion: fracture onset occurs when the total configurational force reaches a critical value, F_CNF = F_CNF,C. This force is computed by summing nodal configurational forces derived from the Eshelby stress inside a small cylindrical region of radius 0.1 mm around the crack tip, and is presented as a computationally efficient estimate of the J-integral under finite strains and biaxial loading. The authors report F_CNF,C values of 0.0386 N/m for Elastosil 2:1, 0.0994 N/m for Elastosil 8:5, 0.2149 N/m for Elastosil 1:1, 0.1286 N/m for Sylgard, and 0.4228 N/m for VHB. They note that Sylgard's intermediate critical force combined with its low total work of fracture indicates brittle behavior, while VHB's high values mark a highly deformable fracture response. For the sideways-fracturing Elastosil 1:1, the method yields a forward-oriented configurational force, which the authors attribute to the isotropic hyperelastic model lacking anisotropic fracture resistance.

Load-bearing premise

The numerical estimate of the total configurational force equals the true J-integral for the chosen mesh and the 0.1 mm summation radius, an equality inherited from prior work rather than demonstrated by a convergence study in this paper.

Editorial extensions

If this is right

  • Fracture onset under biaxial loading could be predicted without contour integration or crack-growth simulation by stopping at the loading stage where the total configurational force reaches the material's critical value.
  • The reported critical values provide a finite-strain toughness ranking across a stiffness range of roughly 5 to 350 kPa, with VHB about eleven times tougher than Elastosil 2:1 by this measure.
  • The data-adaptive optimization identifies \bar I_1-only energy functions as sufficient for the biaxial constitutive response of these elastomers, supporting simpler calibrations for this class of materials.
  • The same virtual testbed could be reused for non-equi-biaxial loading paths and, with added anisotropic fracture resistance, for post-onset crack propagation predictions.

Reading between the lines

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

  • Because the paper states that nodal configurational forces depend on the FE mesh yet reports no mesh-convergence study, the specific F_CNF,C values should be treated as provisional until a denser-mesh verification stabilizes them.
  • A natural next experiment would be to apply the same B-spline calibration to uniaxial or shear data and test whether the resulting \bar I_1-only energy reproduces the biaxial force-displacement curves, which would probe the transferability of the data-adaptive model across loading modes.
  • The reported consistency between F_CNF,C and total work of fracture suggests that the method could rank fracture toughness in other soft solids, but the directional limitation seen in Elastosil 1:1 implies that the magnitude, not the crack-growth direction, is the trustworthy output for isotropic models.
  • Comparing F_CNF,C against an independent contour J-integral evaluated on the same finite-element solution would directly test whether the 0.1 mm summation radius is the load-bearing numerical choice.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

4 major / 6 minor

Summary. The manuscript combines original equi-biaxial fracture experiments on five soft elastomers (Elastosil P7670 at three mixing ratios, Sylgard 184, and VHB 4905) with a data-adaptive hyperelastic constitutive framework and a post-processing Configurational Force Method. The constitutive model represents the isochoric strain energy as B-spline functions of the invariants, calibrated by finite element model updating against DIC displacement fields and reaction forces. The authors then compute the total crack-tip configurational force as an estimator of the J-integral and report critical values at crack onset, proposing a fracture criterion based on a critical configurational force FCNF,C. The paper also provides a public dataset on Zenodo.

Significance. If the quantitative claims are robust, the paper would offer a useful finite-strain alternative to contour-integral J evaluation for soft elastomers under multiaxial loading, and the B-spline data-adaptive constitutive framework is a flexible contribution. The experimental dataset, including full-field DIC data and force-displacement curves for five materials, is a valuable community resource. The authors are transparent about limitations, including the ad hoc I1-only reduction, the directional shortcoming for Elastosil 1:1, and a discrepancy with an earlier value in [52]. However, the central fracture-criterion claim is not independently validated, and the numerical convergence of the reported FCNF,C values is not established.

major comments (4)
  1. [§4.2, Table 1] The proposed fracture criterion in §4.2 is not demonstrated as predictive. The critical values FCNF,C are computed by running each forward simulation up to the experimentally measured average displacement at crack onset (Table 1) and reading off the configurational force at that loading stage. Therefore, FCNF,C is a calibration output on the same experiments used to fit the constitutive model, not a prediction of onset. No test applies FCNF,C to a different geometry, crack length, or loading condition, and no independent toughness comparison is provided. I recommend either an out-of-sample predictive test or a revised framing that presents FCNF,C as a descriptive characterization rather than a fracture criterion.
  2. [§4.1, §4.2] The numerical robustness of FCNF,C is not established. The paper states in §4.2 that nodal configurational forces depend on the FE mesh, yet no mesh-convergence study is reported. Likewise, the choice of a 0.1 mm cylinder radius for summing spurious configurational forces (§4.1) is not justified by a radius-sensitivity study. Since the reported critical values and the material ranking depend on these two numerical choices, a convergence and sensitivity analysis (e.g., two additional mesh densities and radii) is needed before the quantitative values can be relied upon.
  3. [Table 1, §4.2] Uncertainty in the experimentally measured onset displacement is not propagated to FCNF,C. For example, VHB has onset displacement 48.8 ± 13.7 mm and Sylgard 5.6 ± 0.4 mm, but FCNF,C is reported as a single number computed at the mean onset displacement. Because the configurational force is a nonlinear function of displacement, the value at the mean need not equal the mean value, and the ordering of materials (e.g., Sylgard 0.1286 N/m vs Elastosil 1:1 0.2149 N/m) may not be robust. The authors should compute FCNF,C for each experimental repetition and report the distribution, or provide a sensitivity analysis.
  4. [§3.3, §3.4] The data-adaptive constitutive functions are validated only in-sample: the force and displacement errors in Figures 7–11 are on the same biaxial experiments used for calibration. Since the accuracy of the configurational forces depends directly on the constitutive model, an independent check (e.g., prediction of a uniaxial or different biaxial test, or comparison with an independent modulus or toughness measurement) would substantially strengthen the claim. The paper itself acknowledges in §3.4 that the I1-only reduction is ad hoc, which further motivates such a check.
minor comments (6)
  1. [§3 heading] The heading contains a typo: “postpocessing” should be “post-processing.”
  2. [§5.1] The phrase “fracture fracture” appears in the Discussion; one of the two words should be removed.
  3. [§3.3] The word “compreshensive” is misspelled; it should be “comprehensive.”
  4. [§4.2] The term “configurational-forced-based” should be “configurational-force-based.”
  5. [Figure 12] The caption lists panels (a), (b.1), and (b.2), but the text refers to panel (c) for the critical-force bar plot; the panel labels should be made consistent.
  6. [Eqs. (10)–(11)] The objective function is split across two numbered equations with dangling terms; it would be clearer to present it as a single equation or to explicitly reference the undisplayed remainder.

Circularity Check

1 steps flagged · score 6.0 of 10

The fracture criterion is calibrated at the experimentally measured onset displacement, so FCNF = FCNF,C 'predicts' onset by reproducing its own input; no independent validation is provided.

  1. fitted input called prediction [Section 4.2 (Results), Figure 12 caption, 'Total configurational force at crack onset yields a fracture criterion' paragraph]
    "The analysis is carried out up to the point at which fracture propagation begins, as determined experimentally from the average displacement at crack onset (see Table 1). ... The loading stage selected in the virtual experiment (computation) is exactly the one corresponding to the displacement at crack onset (average values in Table 1). ... This leads to what we term a configurational-forced-based fracture criterion, in which fracture onset (not complete rupture) is predicted when the magnitude of the total configurational force reaches a critical value, FCNF = FCNF,C."

    FCNF,C is not predicted or measured independently: it is read off from the simulation at the experimentally known onset displacement, and the paper states 'These values correspond to the final points of the force-displacement curves in Figure 12.b.' The criterion FCNF = FCNF,C therefore asserts that onset occurs when the simulation reaches the value it was forced to attain at the measured onset, so the 'prediction' reproduces its calibration input by construction for every material. The constitutive energy functions used in the same model were fitted to the same biaxial experiments, so ranking FCNF,C against the measured work of fracture is an internal-consistency check within one dataset, not an independent test.

full rationale

The single clear circular step is the fracture criterion: Section 4.2 defines FCNF,C as the total crack-tip configurational force evaluated at the experimentally measured average onset displacement (Table 1), then describes onset 'is predicted' when FCNF = FCNF,C. Because the critical value is extracted exactly where onset was already known, the rule's output duplicates its input for all five materials; this is a fitted input called a prediction. The constitutive modeling itself is not circular: the B-spline energy functions are calibrated by FEMU against pre-crack reaction forces and DIC displacement fields, with fit quality checked (Figures 7-11), and the configurational force genuinely depends on those fitted functions (footnote 13 acknowledges the VHB value changed with the improved model, indicating real dependence). The J-integral equivalence is imported from the authors' own prior work [52], but the same equivalence is established in the wider configurational-mechanics literature (refs. [53, 55]), so that self-citation is corroborated rather than load-bearing in a way that would raise the score. The missing mesh-convergence study and the fixed 0.1 mm summation radius ('Nodal configurational forces are ... a consequence of the FE spatial discretization') are correctness risks rather than circular steps, and the acknowledged failure to predict the sideways fracture direction for Elastosil 1:1 is an honest falsification of the direction claim, not a tautology. The underlying configurational-force framework retains independent mathematical content, so the circularity is partial: the central 'predictive' criterion reduces by construction to the measured onset it claims to forecast, with no independent validation on other geometries, crack lengths, or loading conditions, warranting score 6.

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

The central new numerical content is a set of fitted B-spline control values and domain boundaries; these are legitimate constitutive parameters but they are fitted to the same data used to evaluate the fracture criterion. No new physical entities are introduced.

free parameters (5)
  • B-spline interpolation values w(1)_i (I1-only strain energy function) = Not tabulated; shown as curves in Figs. 7-11a
    These are the principal unknowns in the FEMU calibration (Eqs. 10-13), fitted to reaction forces and DIC displacement fields for each of the five materials.
  • Last interpolation point I1^(n1) per material (domain boundary) = Varies by material: 5.76 (Elastosil 2:1), 5.06 (8:5), 4.35 (1:1), 3.07 (Sylgard), 9.77 (VHB)
    The interpolation domain is also optimized, with a KDE-based regularization pushing it toward the 95th percentile of sampled invariants (Eq. 11, Section B).
  • Cylinder radius for summing spurious configurational forces = 0.1 mm
    Chosen in Section 4.1 for all materials; no mesh or radius convergence study is provided, and nodal configurational forces are stated to be mesh-dependent.
  • Objective function weights omega_f and omega_u = 1/max(f_exp) and 1/max(u_exp)
    Normalization scheme in Eq. (12) to combine force and displacement residuals.
  • ROI dimensions for sampling DIC displacements = 20 mm x 20 mm, except 17 mm x 17 mm for Elastosil 8:5
    Manual choice of calibration region around the crack; its influence on the identified energy functions is not studied.
assumptions (3)
  • domain assumption The elastomers are incompressible and isotropic hyperelastic; volumetric behavior is enforced by a Lagrange multiplier p (Eqs. 4-6).
    Used to set up the forward boundary value problem (Section 3.1). Sideways fracture of Elastosil 1:1 suggests deformation-induced anisotropy that the isotropic model cannot capture.
  • standard math The configurational force at the crack tip equals the J-integral under finite strains.
    Assumed from prior work [52] and used to interpret F_CNF as J (Sections 1.1 and 4.1). This equivalence is not re-derived here.
  • domain assumption The DIC displacement fields and machine reaction forces are accurate and the FE model reproduces the experimental boundary conditions with a 45-degree notch.
    All calibration and fracture metrics rest on this data-fidelity assumption (Section 3.3).

how reviews work

0 comments
Cite this review

Pith. "Pith review of Biaxial characterization of soft elastomers: experiments and data-adaptive configurational forces for fracture." pith.science (2026). https://pith.science/paper/3WUVJKT4

@misc{pith2026250520244,
  author       = {Pith},
  title        = {Pith review of: Biaxial characterization of soft elastomers: experiments and data-adaptive configurational forces for fracture},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/3WUVJKT4}},
  note         = {Machine review of arXiv:2505.20244}
}
read the original abstract

Understanding the fracture mechanics of soft solids remains a fundamental challenge due to their complex, nonlinear responses under large deformations. While multiaxial loading is key to probing their mechanical behavior, the role of such loading in fracture processes is still poorly understood. Here, we present a combined experimental-computational framework to investigate fracture in soft elastomers under equi-biaxial loading. We report original equi-biaxial quasi-static experiments on five elastomeric materials, revealing a spectrum of material and fracture behavior, from brittle-like to highly deformable response with crack tip strains exceeding 150 %. Motivated by these observations, we develop a hybrid computational testbed that mirrors the experimental setup and enables virtual biaxial tests. Central to this framework are two components: a data-adaptive formulation of hyperelastic energy functions that flexibly captures material behavior, and a post-processing implementation of the Configurational Force Method, providing a computationally efficient estimate of the J-integral at the crack tip. Our data-adaptive framework for hyperelastic energy functions proves versatility to capture with high accuracy the hyperelastic behavior observed in the biaxial experiments. This is important because accurately capturing the constitutive behaviour of soft solids is key for a reliable application of the Configurational Force Method to soft solids. In the limit of crack onset, a critical value of the crack tip configurational force allows for a criterion of fracture toughness. Together, our experimental, theoretical, and computational contributions offer a new paradigm for characterizing and designing soft materials with tailored fracture properties.

Figures

Figures reproduced from arXiv: 2505.20244 by the authors.

Figure 1
Figure 1. Experimental setup for biaxial stretch tests. (a) Close-up view of the grips used in the experimental testbed. Each grip consists of three independent clamps designed to hold rectangular specimens. A sliding grip system allows the gripping fingers to translate perpendicular to the loading direction, i.e., in the direction of the edge of the sample. To that end, the gripping fingers can freely move in this direction.… view at source ↗
Figure 2
Figure 2. Experimental results for the biaxial experiments. The two independent axes move identically at a displacement rate of 0.85 mm s−1 , which renders a strain rate of 0.01 s−1 . (a.1-5) Force displacement results for biaxial experiments and for each of the axes of the machine (blue and orange colored-lines). The results for each experimental repetition are marked with ⃝1 , ⃝2 , ⃝3 , and ⃝4 so that they can be easily ide… view at source ↗
Figure 3
Figure 3. Experimental total work of fracture and strain at failure for the biaxial experiments. (a) Work of fracture for each material calculated as the total work performed by both axes during the deformation of the sample until full propagation of the crack divided by the initial volume of the sample. The volume of the samples is 85 mm × 85 mm × 2 mm (thickness) for all samples except for VHB. For VHB, thickness is 0.5 mm.… view at source ↗
Figures from the paper (14 more)
Figure 4
Figure 4. Figure 4: Strain fields at crack onset. The fields are engineering strain fields computed from the Lagrange strain tensor. The fields correspond to the loading stages described in [PITH_FULL_IMAGE:figures/full_fig_p013_4.png]
Figure 5
Figure 5. Figure 5: Schematics of the data-adaptive hyperelastic constitutive modeling framework. The invariants I¯1, I¯2 are discretized with interpolation points (black circles) and the corresponding interpolation values (blue circles) represent unknown material parameters. Interpolatio…
Figure 6
Figure 6. Figure 6: Region of interest used to sample the experimental displacement fields to the FE mesh and calibrate data-adaptive strain energy functions. A square region of interest (ROI) with dimensions 20 mm x 20 mm in the vicinity of the crack contour is defined. The ROI is chosen…
Figure 7
Figure 7. Figure 7: Data-adaptive strain energy function for Elastosil 2:1 obtained from biaxial experimental data: reaction forces and full-field DIC measurements. (a) Data-adaptive hyperelastic strain energy density as a function of the first invariant, I¯1. The energy function is recon…
Figure 8
Figure 8. Figure 8: Data-adaptive strain energy function for Elastosil 8:5 obtained from biaxial experimental data: reaction forces and full-field DIC measurements. (a) Data-adaptive hyperelastic strain energy density as a function of the first invariant, I¯1. The energy function is recon…
Figure 9
Figure 9. Figure 9: Data-adaptive strain energy function for Elastosil 1:1 obtained from biaxial experimental data: reaction forces and full-field DIC measurements. (a) Data-adaptive hyperelastic strain energy density as a function of the first invariant, I¯1. The energy function is recon…
Figure 10
Figure 10. Figure 10: Data-adaptive strain energy function for Sylgard obtained from biaxial experimental data: reaction forces and full-field DIC measurements. (a) Data-adaptive hyperelastic strain energy density as a function of the first invariant, I¯1. The energy function is reconstruc…
Figure 11
Figure 11. Figure 11: Data-adaptive strain energy function for VHB obtained from biaxial experimental data: reaction forces and full-field DIC measurements. (a) Data-adaptive hyperelastic strain energy density as a function of the first invariant, I¯1. The energy function is reconstructed …
Figure 12
Figure 12. Figure 12: Results for data-adaptive configurational forces ante crack onset and fracture criterion based on their critical values. (a) Nodal configurational forces at the crack tip and spurious configurational forces at the crack tip vicinity are depicted for Elastosil 2:1, Ela…
Figure 13
Figure 13. Figure 13: Strain fields at onset of crack growth. The fields are engineering strain fields computed from the Lagrange strain tensor. The fields correspond to the loading stages described in [PITH_FULL_IMAGE:figures/full_fig_p035_13.png]
Figure 14
Figure 14. Figure 14: Strain fields at onset of crack growth. The fields are engineering strain fields computed from the Lagrange strain tensor. The fields correspond to the loading stages described in [PITH_FULL_IMAGE:figures/full_fig_p036_14.png]
Figure 15
Figure 15. Figure 15: Strain fields at onset of crack growth. The fields are engineering strain fields computed from the Lagrange strain tensor. The fields correspond to the loading stages described in [PITH_FULL_IMAGE:figures/full_fig_p037_15.png]
Figure 16
Figure 16. Figure 16: Kernel Density Estimation (KDE) employed to obtain a differentiable representation of the probability density function (PDF) for the sampled invariants I¯1, I¯2. At each optimization iteration, the sampled invariants (¯I1, ¯I2) at the quadrature points of the FE-mesh …
Figure 17
Figure 17. Figure 17: Sampled invariants predicted by the FE forward simulation using the identified data-adaptive hyperelastic strain energy functions. The vertical line marks the 95th-percentile of the invariants. 40 [PITH_FULL_IMAGE:figures/full_fig_p040_17.png]

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Configurational forces explain echelon cracks in soft materials

    cond-mat.soft 2025-07 conditional novelty 6.0 of 10

    Configurational force computations on idealized crack geometries show that facet tilting, local binormal shear, and facet spacing govern echelon crack growth and coalescence in soft materials.

Reference graph

Works this paper leans on

117 extracted references · 70 canonical work pages · cited by 1 Pith paper

  1. [52]

    Moreno-Mateos, M. A. & Steinmann, P. Configurational force method enables fracture assessment in soft materials. Journal of the Mechanics and Physics of Solids 186, 105602 (2024)

  2. [1]

    Designing toughness and strength for soft materials

    Zhao, X. Designing toughness and strength for soft materials. Proceedings of the National Academy of Sciences 114, 8138–8140 (2017)

  3. [2]

    A., Hossain, M., Steinmann, P

    Moreno-Mateos, M. A., Hossain, M., Steinmann, P. & Garcia-Gonzalez, D. Hybrid magnetorheological elastomers enable versatile soft actuators. npj Computational Materials 8, 162 (2022)

  4. [3]

    U., Khalid, M

    Arif, Z. U., Khalid, M. Y ., Tariq, A., Hossain, M. & Umer, R. 3d printing of stimuli-responsive hydrogel materials: Literature review and emerging applications. Giant 17, 100209 (2024)

  5. [4]

    Treloar, L. R. G. Stress-strain data for vulcanised rubber under various types of deformation. Trans. Faraday Soc. 40, 59–70 (1944)

  6. [5]

    G., Green, A

    James, A. G., Green, A. & Simpson, G. M. Strain energy functions of rubber. i. characterization of gum vulcanizates. Journal of Applied Polymer Science 19, 2033–2058 (1975)

  7. [6]

    & Kawai, H

    Kawabata, S., Matsuda, M., Tei, K. & Kawai, H. Experimental survey of the strain energy density function of isoprene rubber vulcanizate. Macromolecules 14, 154–162 (1981)

  8. [7]

    Pancheri, F. Q. & Dorfmann, L. Strain-controlled biaxial tension of natural rubber: New experimental data. Rubber Chemistry and Technology 87, 120–138 (2014)

Show all 117 references
  1. [8]

    Ahmad, D., Sahu, S. K. & Patra, K. Fracture toughness, hysteresis and stretchability of dielectric elastomers under equibiaxial and biaxial loading. Polymer Testing79, 106038 (2019)

  2. [9]

    N., Mai, T.-T., Urayama, K., Gong, J

    Khiêm, V . N., Mai, T.-T., Urayama, K., Gong, J. P. & Itskov, M. A multiaxial theory of double network hydrogels. Macromolecules 52, 5937–5947 (2019). Doi: 10.1021/acs.macromol.9b01044

  3. [10]

    & Hossain, M

    Esmaeili, A., George, D., Masters, I. & Hossain, M. Biaxial experimental characterizations of soft polymers: A review. Polymer Testing128, 108246 (2023)

  4. [11]

    Griffith, A. A. The phenomena of rupture and flow in solids. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences 221, 163–198 (1921)

  5. [12]

    Analysis of stresses and strains near the end of a crack traversing a plate

    Irwin, G. Analysis of stresses and strains near the end of a crack traversing a plate. Journal of Applied Mechanics 24, 361–364 (1957)

  6. [13]

    Cherepanov, G. P. Crack propagation in continuous media. Journal of Applied Mathematics and Mechanics 31, 503–512 (1967)

  7. [14]

    Rice, J. R. A path independent integral and the approximate analysis of strain concentration by notches and cracks. Journal of Applied Mechanics 35, 379 (1968)

  8. [15]

    F., Moran, B

    Shih, C. F., Moran, B. & Nakamura, T. Energy release rate along a three-dimensional crack front in a thermally stressed body. International Journal of Fracture30, 79–102 (1986). 42

  9. [16]

    Unstable crack propagation in metals: Cleavage and fast fracture.Proceedings of the Crack Propagation Symposium 1, paper 84 (1961)

    Wells, A. Unstable crack propagation in metals: Cleavage and fast fracture.Proceedings of the Crack Propagation Symposium 1, paper 84 (1961). Crack Tip Opening Displacement (CTOD)

  10. [17]

    & Sih, G

    Erdogan, F. & Sih, G. C. On the crack extension in plates under plane loading and transverse shear. Journal of Basic Engineering 85, 519–525 (1963)

  11. [18]

    Sih, G. C. Strain-energy-density factor applied to mixed mode crack problems. International Journal of Fracture 10, 305–321 (1974)

  12. [19]

    & Underwood, J

    Hussain, M., Pu, S. & Underwood, J. Strain energy release rate for a crack under combined mode i and mode ii. National Symposium on Fracture Mechanics 2–2–27 (1974)

  13. [20]

    Strifors, H. C. A generalized force measure of conditions at crack tips. International Journal of Solids and Structures 10, 1389–1404 (1974)

  14. [21]

    Y ., Jagota, A., Bennison, S

    Hui, C. Y ., Jagota, A., Bennison, S. J. & Londono, J. D. Crack blunting and the strength of soft elastic solids. Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences 459, 1489–1516 (2003)

  15. [22]

    & Long, R

    Qi, Y ., Zou, Z., Xiao, J. & Long, R. Mapping the nonlinear crack tip deformation field in soft elastomer with a particle tracking method. Journal of the Mechanics and Physics of Solids 125, 326–346 (2019)

  16. [23]

    & Long, R

    Lu, Y ., Qi, Y ., Tenardi, M. & Long, R. Mixed-mode fracture in a soft elastomer.Extreme Mechanics Letters 48, 101380 (2021)

  17. [24]

    Long, R., Krishnan, V . R. & Hui, C. Y . Finite strain analysis of crack tip fields in incompressible hyperelastic solids loaded in plane stress. Journal of the Mechanics and Physics of Solids 59, 672–695 (2011)

  18. [25]

    & Hui, C.-Y

    Long, R. & Hui, C.-Y . Crack tip fields in soft elastic solids subjected to large quasi-static deformation — a review. Extreme Mechanics Letters 4, 131–155 (2015)

  19. [26]

    Eshelby, J. D. The force on an elastic singularity. Philosophical Transactions of the Royal Society of London 244 (1951)

  20. [27]

    Application of material forces to hyperelastostatic fracture mechanics

    Steinmann, P. Application of material forces to hyperelastostatic fracture mechanics. i. continuum mechanical setting. International Journal of Solids and Structures 37, 7371–7391 (2000)

  21. [28]

    Configurational forces: are they needed? Mechanics Research Communications 29, 513–519 (2002)

    Podio-Guidugli, P. Configurational forces: are they needed? Mechanics Research Communications 29, 513–519 (2002)

  22. [29]

    & Schmidt, I

    Gross, D., Kolling, S., Mueller, R. & Schmidt, I. Configurational forces and their application in solid mechanics. European Journal of Mechanics - A/Solids 22, 669–692 (2003)

  23. [30]

    On boundary potential energies in deformational and configurational mechanics

    Steinmann, P. On boundary potential energies in deformational and configurational mechanics. Journal of the Mechanics and Physics of Solids 56, 772–800 (2008). 43

  24. [31]

    & Denzer, R

    Steinmann, P., Scherer, M. & Denzer, R. Secret and joy of configurational mechanics: From foundations in continuum mechanics to applications in computational mechanics. ZAMM - Journal of Applied Mathematics and Mechanics / Zeitschrift für Angewandte Mathematik und Mechanik 89,...

  25. [32]

    Spatial and Material Forces in Nonlinear Continuum Mechanics, vol

    Steinmann, P. Spatial and Material Forces in Nonlinear Continuum Mechanics, vol. 272 (Springer International Publishing, 2022)

  26. [33]

    Francfort, G. A. & Marigo, J. J. Revisiting brittle fracture as an energy minimization problem. Journal of the Mechanics and Physics of Solids 46, 1319–1342 (1998)

  27. [34]

    Schreyer, H. L. Analytical solutions for nonlinear strain-gradient softening and localization. Journal of Applied Mechanics 57, 522–528 (1990)

  28. [35]

    & Geers, M

    Peerlings, R., de Borst, R., Brekelmans, W. & Geers, M. Gradient-enhanced damage modelling of concrete fracture. Mechanics of Cohesive-frictional Materials 3, 323–342 (1998)

  29. [36]

    Nonlocal models for damage and fracture: Comparison of approaches

    Jirásek, M. Nonlocal models for damage and fracture: Comparison of approaches. International Journal of Solids and Structures 35, 4133–4145 (1998)

  30. [37]

    H., Wang, H

    Song, J. H., Wang, H. & Belytschko, T. A comparative study on finite element methods for dynamic fracture. Computational Mechanics 42, 239–250 (2008)

  31. [38]

    Yielding of steel sheets containing slits

    Dugdale, D. Yielding of steel sheets containing slits. Journal of the Mechanics and Physics of Solids 8, 100–104 (1960)

  32. [39]

    & Belytschko, T

    Moës, N., Dolbow, J. & Belytschko, T. A finite element method for crack growth without remeshing.International Journal for Numerical Methods in Engineering (1999)

  33. [40]

    Y ., Gu, L

    Belytschko, T., Lu, Y . Y ., Gu, L. & Tabbara, M. Element-free galerkin methods for static and dynamic fracture. International Journal of Solids and Structures 32, 2547–2570 (1995)

  34. [41]

    & Hofacker, M

    Miehe, C., Welschinger, F. & Hofacker, M. Thermodynamically consistent phase-field models of fracture: Variational principles and multi-field fe implementations. International Journal for Numerical Methods in Engineering 83, 1273–1311 (2010)

  35. [42]

    Kumar, A., Bourdin, B., Francfort, G. A. & Lopez-Pamies, O. Revisiting nucleation in the phase-field approach to brittle fracture. Journal of the Mechanics and Physics of Solids 142, 104027 (2020)

  36. [43]

    Lo, Y .-S., Hughes, T. J. & Landis, C. M. Phase-field fracture modeling for large structures. Journal of the Mechanics and Physics of Solids 105118 (2022)

  37. [44]

    A., Hossain, M., Steinmann, P

    Moreno-Mateos, M. A., Hossain, M., Steinmann, P. & Garcia-Gonzalez, D. Hard magnetics in ultra-soft magnetorheological elastomers enhance fracture toughness and delay crack propagation.Journal of the Mechanics and Physics of Solids 173, 105232 (2023)

  38. [45]

    Moreno-Mateos, M. A. & Steinmann, P. Crosslinking degree variations enable programming and controlling soft fracture via sideways cracking. npj Computational Materials 10, 282 (2024). 44

  39. [46]

    & Kästner, M

    Dammaß, F., Schab, D., Rohm, H. & Kästner, M. Rate- and temperature-dependent ductile-to-brittle fracture transition: Experimental investigation and phase-field analysis for toffee. Engineering Fracture Mechanics 297, 109878 (2024)

  40. [47]

    Hu, X. & Li, S. A stress-intensity-factor-driven phase field modeling of mixed mode fracture.Computer Methods in Applied Mechanics and Engineering 443, 118058 (2025)

  41. [48]

    & Ortiz, M

    Schmidt, B., Fraternali, F. & Ortiz, M. Eigenfracture: An eigendeformation approach to variational fracture. https://doi.org/10.1137/080712568 7, 1237–1266 (2009)

  42. [49]

    & Kaliske, M

    Storm, J. & Kaliske, M. The origin of the energy split in phase-field fracture and eigenfracture. PAMM 23, e202300295 (2023)

  43. [50]

    Silling, S. A. Reformulation of elasticity theory for discontinuities and long-range forces. Journal of the Mechanics and Physics of Solids 48, 175–209 (2000)

  44. [51]

    Javili, A., Ekiz, E., McBride, A. T. & Steinmann, P. Continuum-kinematics-inspired peridynamics: Thermo- mechanical problems. Continuum Mechanics and Thermodynamics 33, 2039–2063 (2021)

  45. [53]

    & Barth, F

    Steinmann, P., Ackermann, D. & Barth, F. Application of material forces to hyperelastostatic fracture mechanics. ii. computational setting. International Journal of Solids and Structures 38, 5509–5526 (2001)

  46. [54]

    & Haslach, H

    Kienzler, R., Herrmann, G. & Haslach, H. Mechanics in material space: With applications to defect and fracture mechanics. Applied Mechanics Reviews 55, B23–B24 (2002)

  47. [55]

    Denzer, R., Barth, F. J. & Steinmann, P. Studies in elastic fracture mechanics based on the material force method. International Journal for Numerical Methods in Engineering 58, 1817–1835 (2003)

  48. [56]

    Gurtin, M. E. & Podio-Guidugli, P. Configurational forces and the basic laws for crack propagation. Journal of the Mechanics and Physics of Solids 44, 905–927 (1996)

  49. [57]

    M., Xu, Y ., Zhang, P

    Zhou, R., Coombs, W. M., Xu, Y ., Zhang, P. & Wang, L. G. A configurational force-based material point method for crack propagation modelling in 2d. Theoretical and Applied Fracture Mechanics117, 103186 (2022)

  50. [58]

    Curved crack propagation based on configurational forces

    Schütte, H. Curved crack propagation based on configurational forces. Computational Materials Science 46, 642–646 (2009)

  51. [59]

    M., Pletz, M

    Frankl, S. M., Pletz, M. & Schuecker, C. Improved concept for iterative crack propagation using configurational forces for targeted angle correction. Engineering Fracture Mechanics 266, 108403 (2022)

  52. [60]

    & Ricoeur, A

    Schmitz, K. & Ricoeur, A. Theoretical and computational aspects of configurational forces in three-dimensional crack problems. International Journal of Solids and Structures 282, 112456 (2023)

  53. [61]

    Serrao, P. H. & Kozinov, S. Evaluation of configurational/material forces in strain gradient elasticity theory. Mechanics of Materials 203, 105240 (2025). 45

  54. [62]

    Guo, Y . & Li, Q. Material configurational forces applied to mixed mode crack propagation. Theoretical and Applied Fracture Mechanics 89, 147–157 (2017)

  55. [63]

    Liu, R., Hou, J. & Li, Q. Material configurational forces applied to mixed-mode fatigue crack propagation and life prediction in elastic-plastic material. International Journal of Fatigue 134, 105467 (2020)

  56. [64]

    & Müller, R

    Yan, S., Schlüter, A. & Müller, R. Configurational forces in a phase field model for fatigue. PAMM 22, e202200034 (2023)

  57. [65]

    D., Govindjee, S., Klein, P

    Nguyen, T. D., Govindjee, S., Klein, P. A. & Gao, H. A material force method for inelastic fracture mechanics. Journal of the Mechanics and Physics of Solids 53, 91–121 (2005)

  58. [66]

    & Müller, R

    Näser, B., Kaliske, M. & Müller, R. Material forces for inelastic models at large strains: Application to fracture mechanics. Computational Mechanics 40, 1005–1013 (2007)

  59. [67]

    & Näser, B

    Kaliske, M., Netzker, C. & Näser, B. Evaluation of crack-driving forces at finite viscoelasticity: Theory and experiment. IUTAM Bookseries 17, 193–202 (2009)

  60. [68]

    & Bhashyam, G

    Özenç, K., Kaliske, M., Lin, G. & Bhashyam, G. Evaluation of energy contributions in elasto-plastic fracture: A review of the configurational force approach. Engineering Fracture Mechanics 115, 137–153 (2014)

  61. [69]

    & Kegl, M

    Kolednik, O., Tiwari, A., Posch, C. & Kegl, M. Configurational force based analysis of creep crack growth. International Journal of Fracture236, 175–199 (2022)

  62. [70]

    & Menzel, A

    Denzer, R. & Menzel, A. Configurational forces for quasi-incompressible large strain electro-viscoelasticity – application to fracture mechanics. European Journal of Mechanics - A/Solids 48, 3–15 (2014)

  63. [71]

    A., Mehnert, M

    Moreno-Mateos, M. A., Mehnert, M. & Steinmann, P. Electro-mechanical actuation modulates fracture perfor- mance of soft dielectric elastomers. International Journal of Engineering Science 195, 104008 (2024)

  64. [72]

    Serrao, P. H. & Kozinov, S. Configurational forces for defect evolution in flexoelectricity. Engineering Fracture Mechanics 110903 (2025)

  65. [73]

    & Javili, A

    Steinmann, P., de Villiers, A., McBride, A. & Javili, A. Configurational peridynamics. Mechanics of Materials 185, 104751 (2023)

  66. [74]

    & Possart, G

    Steinmann, P., Hossain, M. & Possart, G. Hyperelastic models for rubber-like materials: consistent tangent operators and suitability for Treloar’s data. Archive of Applied Mechanics 82, 1183–1217 (2012)

  67. [75]

    & Wriggers, P

    Ricker, A. & Wriggers, P. Systematic Fitting and Comparison of Hyperelastic Continuum Models for Elastomers. Archives of Computational Methods in Engineering 30, 2257–2288 (2023)

  68. [76]

    A., Açan, A

    Dal, H., Denli, F. A., Açan, A. K. & Kaliske, M. Data-driven hyperelasticity, part i: A canonical isotropic formulation for rubberlike materials. Journal of the Mechanics and Physics of Solids 179, 105381 (2023)

  69. [77]

    Z., Açan, A

    Tikeno˘gulları, O. Z., Açan, A. K., Kuhl, E. & Dal, H. Data-driven hyperelasticity, part ii: A canonical framework for anisotropic soft biological tissues. Journal of the Mechanics and Physics of Solids 181, 105453 (2023). 46

  70. [78]

    & Ortiz, M

    Kirchdoerfer, T. & Ortiz, M. Data-driven computational mechanics. Computer Methods in Applied Mechanics and Engineering 304, 81–101 (2016)

  71. [79]

    & Itskov, M

    Abdusalamov, R., Hillgärtner, M. & Itskov, M. Automatic generation of interpretable hyperelastic material models by symbolic regression. International Journal for Numerical Methods in Engineering 124, 2093–2104 (2023)

  72. [80]

    & Sun, W

    Bahmani, B. & Sun, W. Physics-constrained symbolic model discovery for polyconvex incompressible hypere- lastic materials. International Journal for Numerical Methods in Engineering 125, e7473 (2024)

  73. [81]

    L., Jones, R

    Frankel, A. L., Jones, R. E. & Swiler, L. P. Tensor Basis Gaussian Process Models of Hyperelastic Materials. ArXiv abs/1912.10872 (2019)

  74. [82]

    & Gil, A

    Ellmer, N., Ortigosa, R., Martínez-Frutos, J. & Gil, A. J. Gradient enhanced gaussian process regression for constitutive modelling in finite strain hyperelasticity. Computer Methods in Applied Mechanics and Engineering 418, 116547 (2024)

  75. [83]

    & Kuhl, E

    Holthusen, H., Lamm, L., Brepols, T., Reese, S. & Kuhl, E. Theory and implementation of inelastic constitutive artificial neural networks. Computer Methods in Applied Mechanics and Engineering 428, 117063 (2024)

  76. [84]

    M., Long, K

    Hamel, C. M., Long, K. N. & Kramer, S. L. B. Calibrating constitutive models with full-field data via physics informed neural networks. Strain 59, e12431 (2023)

  77. [85]

    Dammaß, F., Kalina, K. A. & Kaestner, M. Neural Networks Meet Phase-Field: A Hybrid Fracture Model (2024)

  78. [86]

    Klein, D. K. et al. Neural networks meet hyperelasticity: A monotonic approach (2025). ArXiv:2501.02670 [cs]

  79. [87]

    Ortigosa, R. et al. A generalized theory for physics-augmented neural networks in finite strain thermo-electro- mechanics. Computer Methods in Applied Mechanics and Engineering 437, 117741 (2025)

  80. [88]

    Linka, K. et al. Constitutive artificial neural networks: A fast and general approach to predictive data-driven constitutive modeling by deep learning. Journal of Computational Physics 429, 110010 (2021)

  81. [89]

    & Kuhl, E

    Linka, K. & Kuhl, E. A new family of Constitutive Artificial Neural Networks towards automated model discovery. Computer Methods in Applied Mechanics and Engineering 403, 115731 (2023)

  82. [90]

    Linden, L. et al. Neural networks meet hyperelasticity: A guide to enforcing physics. Journal of the Mechanics and Physics of Solids 179, 105363 (2023)

  83. [91]

    Fuhg, J. N. et al. A Review on Data-Driven Constitutive Laws for Solids. Archives of Computational Methods in Engineering 32, 1841–1883 (2025)

  84. [92]

    Römer, U. et al. Reduced and All-at-Once Approaches for Model Calibration and Discovery in Computational Solid Mechanics (2024). ArXiv:2404.16980 [cs]

  85. [93]

    A., Orteu, J

    Sutton, M. A., Orteu, J. J. & Schreier, H. Image correlation for shape, motion and deformation measurements: basic concepts, theory and applications (Springer Science & Business Media, 2009). 47

  86. [94]

    & Grédiac, M

    Pierron, F. & Grédiac, M. The Virtual Fields Method: Extracting Constitutive Mechanical Parameters from Full-field Deformation Measurements (Springer, New York, NY , 2012)

  87. [95]

    Flaschel, M. et al. Automated discovery of interpretable hyperelastic material models for human brain tissue with EUCLID. Journal of the Mechanics and Physics of Solids 180, 105404 (2023)

  88. [96]

    Joshi, A. et al. Bayesian-EUCLID: Discovering hyperelastic material laws with uncertainties.Computer Methods in Applied Mechanics and Engineering 398, 115225 (2022)

  89. [97]

    Thakolkaran, P. et al. NN-EUCLID: Deep-learning hyperelasticity without stress data. Journal of the Mechanics and Physics of Solids 169, 105076 (2022)

  90. [98]

    Wiesheier, S., Moreno-Mateos, M. A. & Steinmann, P. Versatile data-adaptive hyperelastic energy functions for soft materials. Computer Methods in Applied Mechanics and Engineering 430, 117208 (2024)

  91. [99]

    & Steinmann, P

    Wiesheier, S., Mergheim, J. & Steinmann, P. Discrete data-adaptive approximation of hyperelastic energy functions. Computer Methods in Applied Mechanics and Engineering 416, 116366 (2023)

  92. [100]

    A theory of large elastic deformation

    Mooney, M. A theory of large elastic deformation. Journal of applied physics 11, 582–592 (1940). Publisher: American Institute of Physics

  93. [101]

    Characterization of elastic properties of carbon-black-filled rubber vulcanizates

    Yeoh, O. Characterization of elastic properties of carbon-black-filled rubber vulcanizates. Rubber chemistry and technology 63, 792–805 (1990)

  94. [102]

    d. Boor, C. A Practical Guide to Splines (Springer Verlag, 1978)

  95. [103]

    & Tiller, W

    Piegl, L. & Tiller, W. The NURBS Book (Springer-Verlag, New York, NY , USA, 1996), second edn

  96. [104]

    Dammaß, F., Kalina, K. A. & Kästner, M. When invariants matter: The role of i1 and i2 in neural network models of incompressible hyperelasticity (2025)

  97. [105]

    Arndt, D. et al. The deal.II Library, Version 9.4, 2022

  98. [106]

    Arndt, D. et al. The deal.II finite element library: Design, features, and insights. Computers & Mathematics with Applications 81, 407–422 (2021)

  99. [107]

    & Herrmann, G

    Kienzler, R. & Herrmann, G. On the properties of the eshelby tensor. Acta Mechanica 125, 73–91 (1997)

  100. [108]

    Eshelby, J. D. The elastic energy-momentum tensor. Journal of Elasticity 5, 321–335 (1975)

  101. [109]

    Eshelby, J. D. Energy relations and the energy-momentum tensor in continuum mechanics. Fundamental Contributions to the Continuum Theory of Evolving Phase Interfaces in Solids 82–119 (1999)

  102. [110]

    & Hackl, K

    Stumpf, H., Makowski, J. & Hackl, K. Configurational forces and couples in fracture mechanics accounting for microstructures and dissipation. International Journal of Solids and Structures 47, 2380–2389 (2010)

  103. [111]

    & Kaliske, M

    Behnke, R., Berger, T. & Kaliske, M. Numerical modeling of time- and temperature-dependent strain-induced crystallization in rubber. International Journal of Solids and Structures 141-142, 15–34 (2018). 48

  104. [112]

    & Bouklas, N

    Mulderrig, J., Li, B. & Bouklas, N. Affine and non-affine microsphere models for chain scission in polydisperse elastomer networks. Mechanics of Materials 160, 103857 (2021)

  105. [113]

    Arunachala, P. K. et al. A multiscale anisotropic polymer network model coupled with phase field fracture. International Journal for Numerical Methods in Engineering 125, e7488 (2024)

  106. [114]

    & Grasinger, M

    Mulderrig, J., Buche, M. & Grasinger, M. Polydisperse polymer networks with irregular topologies. Arxiv (2025)

  107. [115]

    & Tang, J

    Li, C., Yang, H., Suo, Z. & Tang, J. Fatigue-resistant elastomers. Journal of the Mechanics and Physics of Solids 134, 103751 (2020)

  108. [116]

    & Govindjee, S

    Reese, S. & Govindjee, S. A theory of finite viscoelasticity and numerical aspects. International Journal of Solids and Structures 35, 3455–3482 (1998)

  109. [117]

    C., Copas, J

    Atkinson, A. C., Copas, J. B., Pierce, D. A., Schervish, M. J. & Titterington, D. M. OXFORD STATISTICAL SCIENCE SERIES. OXFORD STATISTICAL SCIENCE SERIES . 49

Pith tools

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