Pith. sign in

REVIEW 2 major objections 6 minor 68 references

JAX-FEM-ANISO: Differentiable GPU-Accelerated Finite Element Framework for Inverse Identification of Finite-Strain Anisotropic Plasticity

T0 review · 2 major / 6 minor · reviewed 2026-07-14 · grok-4.5

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read A GPU-differentiable finite-element code recovers anisotropic metal plasticity parameters from one full-field displacement test, with multi-fold forward speedups.

desk verdict Solid systems paper: first end-to-end JAX GPU+AD finite-strain anisotropic plasticity (Hill + Yld2004-18p) with real Abaqus speedups and high-dimensional synthetic FEMU that works. read the letter →

arxiv 2606.17390 v2 pith:7NJ7BEAY submitted 2026-06-16 cs.CE

Deepak Sharma , Itzel Salgado , Lu Huang , Hui-Ping Wang , Jian Cao This is my paper · ORCID
classification cs.CE
keywords inversematerialcharacterizationfinite-strainanisotropicplasticitydifferentiablesimulationsJAX-FEMPDE-constrainedoptimizationGPUaccelerationautomaticdifferentiationfull-fieldidentification
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 claims that finite-strain anisotropic plasticity can be simulated and calibrated end-to-end on GPUs by making the entire finite-element workflow differentiable. Automatic differentiation supplies consistent material tangents and adjoint gradients of full-field displacement mismatch, so high-dimensional yield and hardening parameters—including spatially varying fields and Barlat-type coefficients—can be recovered by gradient optimization from a single information-rich specimen rather than many conventional tests. On a three-million-degree-of-freedom problem the same framework runs up to roughly nine times faster than a 24-core commercial baseline. A sympathetic reader cares because metal forming, welding, and additive manufacturing need both fast large-deformation solvers and reliable calibration of anisotropic models; the work argues that differentiability plus GPU throughput removes the two main bottlenecks that have kept one-shot, high-fidelity characterization impractical.

What carries the argument

End-to-end automatic differentiation through the constitutive return map, global Newton solver, and discrete adjoint of the time-stepping residual: custom Jacobian-vector products give consistent local tangents without manual derivation, while reverse-mode adjoints supply objective gradients whose cost is nearly independent of parameter dimension.

What would settle it

Apply the same adjoint pipeline to real digital-image-correlation displacements from a physical information-rich specimen and check whether recovered anisotropic parameters predict independent validation tests within engineering tolerance; failure of that transfer would overturn the single-test claim.

Watch

Extended reading notes

Core claim

A fully differentiable, GPU-accelerated finite-element implementation of finite-strain anisotropic elasto-plasticity (Hill–48 and Barlat Yld2004-18p) recovers anisotropic yield and hardening parameters—including spatially varying fields—from single information-rich full-field displacement tests via adjoint automatic-differentiation gradients, while delivering up to about 9.4× forward speed-up versus a 24-core commercial CPU baseline at roughly three million degrees of freedom.

Load-bearing premise

Every inverse-identification result uses synthetic displacement data generated by the same constitutive family and finite-element model (with optional added noise), not real full-field measurements with imperfect boundaries and model error.

Editorial extensions

If this is right

  • High-dimensional anisotropic yield surfaces become calibratable from one carefully designed full-field test instead of multi-test campaigns.
  • Spatially heterogeneous plasticity parameters induced by manufacturing can be identified by the same adjoint workflow.
  • Forward large-deformation anisotropic simulations at multi-million DOF become practical on a single high-end GPU.
  • Finite-difference sensitivity analysis is no longer required for gradient-based FEMU of complex plasticity models.

Reading between the lines

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

  • If synthetic-to-real transfer holds, the same pipeline could co-optimize specimen topology and loading paths for maximum parameter identifiability.
  • Adjoint gradients of full-field mismatch open a path to joint geometry-and-material design under manufacturing-induced heterogeneity.
  • The framework’s differentiability would also support uncertainty quantification over plasticity parameters once experimental noise models are included.
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

2 major / 6 minor

Summary. The manuscript presents JAX-FEM-ANISO, a fully differentiable GPU-accelerated finite-element framework built on JAX-FEM for forward simulation and inverse identification of finite-strain anisotropic plasticity (Hill–48 and Barlat Yld2004-18p). It implements Miehe’s modular logarithmic-strain formulation with F-bar kinematics, AD-derived consistent tangents and custom JVPs for local return mapping and spectral maps, GPU-resident residual/stiffness evaluation, sparse assembly, and linear solves (AMGX/CuDSS), plus discrete adjoint gradients for FEMU. Forward results are verified against Abaqus (flange drawing; single-element Barlat path). On ~3M DOF, a single H100 yields up to 9.4× wall-clock speedup versus a 24-core Abaqus CPU baseline. Adjoint gradients match central finite differences to <0.2% at h=10⁻⁶. Inverse studies on synthetic full-field displacements recover anisotropic yield/hardening parameters for information-rich tensile and cruciform specimens (with noise), two-region spatial heterogeneity, and a high-dimensional Barlat parameter set from a topology-optimized specimen.

Significance. If the reported accuracy and performance hold, this is a substantial systems and methods contribution for computational solid mechanics and Material Testing 2.0: it removes manual consistent-tangent derivation for advanced anisotropic models, makes high-dimensional FEMU gradients practical via reverse-mode AD, and demonstrates GPU scaling on million-DOF nonlinear problems. Strengths include careful Abaqus verification, explicit FD–AD gradient checks with the expected step-size U-curve, transparent reporting of Voce (Q,b) non-uniqueness while matching hardening curves, and progressive inverse demos up to Barlat Yld2004-18p and multi-material fields. These establish differentiable GPU FEM as a credible engine for simulation and synthetic-data inverse design; real DIC transfer remains future work as stated in §4.

major comments (2)
  1. Abstract and §4 claim a “practical high-throughput engine for … characterization” and that information-rich single tests can replace many conventional experiments. All inverse results (§3.3–3.6) use same-model synthetic displacements (optional i.i.d. Gaussian noise), not real DIC with imperfect BCs, incomplete observations, or model discrepancy—explicitly deferred to future work in §4. The algorithmic claims (AD gradients, GPU speedups, synthetic recovery) are supported; the characterization-engine phrasing overreaches the evidence. Please qualify abstract/conclusions/keywords so “single-test characterization” is framed as demonstrated on synthetic full-field data and as a pathway toward DIC-based MT2.0, not as an established experimental replacement.
  2. §3.5 titles and abstract language refer to “spatially varying material properties,” but the demonstration is two discrete material regions with piecewise-constant parameters (12 scalars), not a continuous field (e.g., nodal or element-wise parameter map). That is still a useful multi-material test and shows AD scaling with parameter count, but it does not yet establish continuous spatial-field identification. Clarify the claim (two-region heterogeneity vs. field reconstruction) and, if continuous fields are intended as a contribution, either add a small field example or soften the wording.
minor comments (6)
  1. Typos/grammar: “largly unexplored” (Introduction); “memeory demands,” “stiffness materix,” “uisng JAX” (§2.1); “V oce” spacing throughout (should be Voce); “non-finite solutions” is fine but check consistency of “finite-difference” hyphenation.
  2. Table 4: Abaqus 24-core reference is on Windows 11 while other Abaqus runs are on RHEL; note this more prominently when interpreting the 9.4× figure so readers do not over-interpret cross-platform wall-clock ratios.
  3. §3.2 gradient check: direction vector Di=0.1 is somewhat arbitrary; a brief note that relative component-wise errors (Table 6) are the primary verification is enough.
  4. §3.6 Barlat inverse: recovered coefficients differ from ground truth while yield surfaces match—good—but state explicitly in the table caption that coefficient non-uniqueness is expected and that surface/hardening agreement is the success metric.
  5. Figures 5, 14, 18, 26: ensure colorbars and deformed-mesh scales are readable in print; some multi-panel layouts are dense.
  6. Related work: briefly position against other differentiable FEM/plasticity efforts beyond JAX-FEM/JAX-CPFEM (e.g., other AD-enabled continuum codes) so novelty of the anisotropic finite-strain + inverse stack is sharper.

Circularity Check

1 steps flagged · score 1.0 of 10

No load-bearing circularity; inverse recoveries are standard same-model synthetic verification, gradients cross-checked by independent FD, and JAX-FEM base is an engineering substrate not a uniqueness claim.

  1. self citation load bearing [Abstract / §1 / §2.2.5 (JAX-FEM base)]
    "Built on JAX-FEM, the framework exploits modern accelerator architectures by parallelizing the three major computational bottlenecks in nonlinear FEM... JAX-FEM provides a general-purpose, GPU-accelerated differentiable finite element framework..."

    The implementation re-uses the authors’ own prior JAX-FEM library for elemental residual/stiffness kernels, sparse assembly and linear solves. This is ordinary engineering reuse of a software substrate rather than a load-bearing uniqueness or ansatz claim; the anisotropic constitutive update, custom JVPs, adjoint formulation and inverse results are independently derived and externally verified against Abaqus and finite differences. Flagged only as a minor self-citation that does not force any scientific conclusion.

full rationale

The paper implements Miehe’s modular logarithmic-strain finite-strain plasticity (Hill–48 and Yld2004-18p) inside an existing differentiable GPU FEM library, supplies consistent tangents and adjoints via custom JVPs and reverse-mode AD, and verifies forward solutions against Abaqus and objective gradients against central finite differences (absolute differences <0.2 % at h=10^{-6}). Inverse examples recover the very parameters used to generate the synthetic full-field displacements (with optional additive Gaussian noise); this is ordinary method verification, not a derivation that forces the result by construction, and the paper itself reports the well-known non-uniqueness of the (Q,b) hardening pair. The sole minor self-reference is the reuse of the authors’ prior JAX-FEM infrastructure for residual/stiffness evaluation and assembly; that infrastructure is an engineering substrate, not a uniqueness theorem or ansatz that closes the logical loop. Real DIC transfer is explicitly left to future work. Consequently the derivation chain is self-contained against external benchmarks and exhibits no circular reduction of the kind enumerated in the analyzer rules.

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

The central claims rest on standard continuum mechanics and FEM practice plus engineering modeling choices (additive logarithmic-strain modular plasticity, F-bar, Voce hardening, associative flow, synthetic DIC-like data). No new physical entities are postulated. Free parameters are the constitutive coefficients being identified and a few algorithmic tolerances/step sizes; they are either optimization unknowns or standard numerical controls, not hidden fits that manufacture the forward verification.

free parameters (4)
  • Constitutive parameter vector θ (Hill r_ij, Voce σ0/Q/b, Barlat c′/c′′ coefficients, m)
    These are the inverse unknowns; ground-truth values are chosen for synthetic experiments and recovered by optimization. Claims of accurate recovery depend on these synthetic targets and bounds (±50% for Barlat demo).
  • Finite-difference step size h and AD regularization δ for spectral map
    h is used only for gradient verification; spectral perturbation δ≪1 is an implementation choice to make eigh differentiable. Both affect numerical gradient quality if poorly chosen.
  • L-BFGS-B convergence tolerances, parameter normalization bounds, noise scale δ_noise
    Optimizer and synthetic-noise settings (Tables 9–12) influence reported iteration counts and recovered errors; not physics free parameters but they shape the inverse results shown.
  • Adaptive time-step bounds (Δt_initial, Δt_min, Δt_max) and line-search β, η
    Numerical controls for Newton robustness; affect forward convergence and thus inverse feasibility, especially for Barlat.
assumptions (5)
  • domain assumption Finite-strain anisotropic plasticity is adequately represented by Miehe’s additive logarithmic-strain modular framework with associative flow and isotropic Voce hardening.
    §2.2 adopts this modular preprocessor/constitutive/postprocessor split; inverse claims assume model form is correct for the material of interest.
  • domain assumption Quasi-static equilibrium without body forces/inertia; F-bar kinematics control volumetric locking.
    Governing equations (1)–(2) and F-bar §2.2.4; standard solid-mechanics modeling choices.
  • standard math Discrete adjoint / reverse-mode AD through the time-discrete residual yields correct dJ/dθ for path-dependent plasticity when custom VJPs use the implicit-function theorem at local and global Newton levels.
    §2.3 and Appendix A; standard adjoint FEMU theory plus JAX custom_vjp practice.
  • ad hoc to paper Synthetic full-field displacements (optionally with i.i.d. Gaussian noise) are a sufficient proxy to demonstrate inverse identification performance relevant to DIC-based Material Testing 2.0.
    All inverse sections §3.3–3.6 use synthetic data; conclusions defer real DIC. This is the load-bearing experimental-transfer assumption.
  • domain assumption Automatic differentiation through spectral logarithmic strain can be made well-defined via small diagonal regularization or analytic Seth–Hill JVPs.
    §2.2.5; required for end-to-end differentiability of the geometric pre/postprocessors.

how reviews work

0 comments
Cite this review

Pith. "Pith review of JAX-FEM-ANISO: Differentiable GPU-Accelerated Finite Element Framework for Inverse Identification of Finite-Strain Anisotropic Plasticity." pith.science (2026). https://pith.science/paper/7NJ7BEAY

@misc{pith2026260617390,
  author       = {Pith},
  title        = {Pith review of: JAX-FEM-ANISO: Differentiable GPU-Accelerated Finite Element Framework for Inverse Identification of Finite-Strain Anisotropic Plasticity},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/7NJ7BEAY}},
  note         = {Machine review of arXiv:2606.17390}
}
abstract

We present a fully differentiable, GPU-accelerated finite element framework JAX-FEM-ANISO, for forward simulation and inverse parameter identification of finite-strain anisotropic plasticity. Built on JAX-FEM, the framework exploits modern accelerator architectures by parallelizing the three major computational bottlenecks in nonlinear FEM: elemental weak-form and tangent-stiffness evaluation, global sparse matrix assembly, and sparse linear solution. For a large-scale forward problem with 3 million degrees of freedom, JAX-FEM-ANISO on a single NVIDIA H100 GPU achieves up to 9.4$\times$ speed-up over a 24-core CPU Abaqus baseline. Automatic differentiation is applied through the constitutive update and solver workflow, providing consistent Jacobians for complex constitutive models without manual derivation and accurate gradients for PDE-constrained inverse analysis. Compared with finite differences, the JAX-AD gradients avoid step-size sensitivity and provide the required sensitivities at substantially lower computational cost. For inverse characterization, we combine information-rich, topology-optimized heterogeneous specimens with full-field displacement data to identify advanced constitutive model parameters from a single test, replacing what would otherwise require many conventional experiments. We demonstrate accurate recovery of anisotropic yield and hardening parameters in progressively challenging settings, including uniform and spatially varying material properties. The resulting AD-based formulation enables efficient optimization in high-dimensional parameter spaces where finite-difference approaches are computationally infeasible. These results establish differentiable, GPU-accelerated FEM as a practical high-throughput engine for simulation, characterization, and optimization workflows in advanced manufacturing.

Figures

Figures reproduced from arXiv: 2606.17390 by the authors.

Figure 1
Figure 1. Overview of the proposed differentiable GPU-accelerated framework for gradient-based parameter identification of finite-strain anisotropic plasticity: (a) parameterization of the constitutive model (M) by a set of material parameters (θ), (b) forward simulation using differentiable FEM, (c) objective function evaluation by comparing the mismatch of the simulated and ground truth displacement fields, (d) gradient com… view at source ↗
Figure 2
Figure 2. Schematic of an elastoplastic solid body with displacement and traction boundary conditions. 2.2.2 Additive Lagrangian Modular Approach Finite-strain plasticity with material anisotropy is essential for modeling sheet metals subjected to complex loading paths. Several algorithmic challenges arise at large strains, including the choice of strain measure, objective stress mapping, robust return mapping for path-depend… view at source ↗
Figure 3
Figure 3. Overview of the differentiable GPU-accelerated forward-solve implementation. while CuDSS is employed as a fallback when the iterative solve does not converge within a prescribed number of iterations. A comparative assessment of these solver choices in terms of computational performance is provided in Appendix D. Automatic time-stepping scheme. To improve robustness in strongly nonlinear loading regimes, the forward … view at source ↗
Figures from the paper (25 more)
Figure 4
Figure 4. Figure 4: Schematic diagram of the thin circular flange considered for the drawing simulations. All dimensions are in mm [PITH_FULL_IMAGE:figures/full_fig_p013_4.png]
Figure 5
Figure 5. Figure 5: Comparison of equivalent plastic strain in flange drawing simulations between JAX-FEM and Abaqus. The deformed meshes are shown at radial displacements of u = 25, 50, and 75 mm. Results correspond to anisotropic responses characterized by (a) δ = 3.4641 and (b) δ = 0.8…
Figure 6
Figure 6. Figure 6: Verification setup for the Barlat Yld2004-18p model. (a) Unit-cube single-element benchmark with the applied boundary conditions. (b) Prescribed displacement history imposed on the highlighted nodes [PITH_FULL_IMAGE:figures/full_fig_p015_6.png]
Figure 7
Figure 7. Figure 7: Comparison of the stress plots obtained using JAX-FEM and Abaqus for the Barlat Yld2004-18p verification test. 3.1.3 Computational Performance Benchmarks To assess the computational performance and scalability of JAX-FEM, we consider the unit-cube geometry shown in [P…
Figure 8
Figure 8. Figure 8: Schematic of the unit-cube example with loading and boundary conditions used for the forward-simulation benchmarks [PITH_FULL_IMAGE:figures/full_fig_p017_8.png]
Figure 9
Figure 9. Figure 9: Speed-up comparison between JAX-FEM and Abaqus for the unit-cube benchmark with 100 × 100 × 100 elements, using the Abaqus CPU (24-core) case as the reference configuration. The Abaqus simulations are performed on different CPU configurations ranging from 1 core to 128…
Figure 10
Figure 10. Figure 10: Speed-up comparison between JAX-FEM and Abaqus for the larger unit-cube benchmark with approximately 11.2 million degrees of freedom, using the Abaqus CPU (24-core) case as the reference configuration [PITH_FULL_IMAGE:figures/full_fig_p018_10.png]
Figure 11
Figure 11. Figure 11: Loading and boundary conditions for the unit-cube domain considered for gradient verification [PITH_FULL_IMAGE:figures/full_fig_p020_11.png]
Figure 12
Figure 12. Figure 12: Finite-difference error check EFD,check as a function of step size h [PITH_FULL_IMAGE:figures/full_fig_p020_12.png]
Figure 13
Figure 13. Figure 13: Information-rich tensile specimen geometry adapted from Zhang et al. [19]. (a) (b) [PITH_FULL_IMAGE:figures/full_fig_p022_13.png]
Figure 14
Figure 14. Figure 14: (a) Distribution of equivalent plastic strain and (b) stress-state distribution for integration points subjected to plastic deformation. Both axes are normalized by the equivalent Hill stress; σ11 and σ22 are aligned with the rolling direction (RD) and transverse dire…
Figure 15
Figure 15. Figure 15: Objective-function histories until convergence of the optimization [PITH_FULL_IMAGE:figures/full_fig_p024_15.png]
Figure 16
Figure 16. Figure 16: Comparison of hardening curves for different material-parameter initializations [PITH_FULL_IMAGE:figures/full_fig_p024_16.png]
Figure 17
Figure 17. Figure 17: Cruciform specimen geometry with the considered loading and boundary conditions [PITH_FULL_IMAGE:figures/full_fig_p026_17.png]
Figure 18
Figure 18. Figure 18: Plastic-strain distribution at the final load step [PITH_FULL_IMAGE:figures/full_fig_p027_18.png]
Figure 19
Figure 19. Figure 19: Visualization of noisy synthetic data used for inverse identification at different noise levels [PITH_FULL_IMAGE:figures/full_fig_p027_19.png]
Figure 20
Figure 20. Figure 20: Convergence history for inverse identification using noise-free synthetic data. (a) (b) [PITH_FULL_IMAGE:figures/full_fig_p028_20.png]
Figure 21
Figure 21. Figure 21: Convergence histories for inverse identification using noisy synthetic data: (a) objective function values and (b) reduction in objective function values measured as the difference between the initial and current objective function values during optimization. 3.5 Inve…
Figure 22
Figure 22. Figure 22: (a) Cruciform geometry with spatially varying material properties; (b) distribution of plastic strain at the final loading step [PITH_FULL_IMAGE:figures/full_fig_p029_22.png]
Figure 23
Figure 23. Figure 23: Convergence history for identifying spatially varying material properties [PITH_FULL_IMAGE:figures/full_fig_p030_23.png]
Figure 24
Figure 24. Figure 24: Comparison of inverse-identified hardening curves with ground-truth results [PITH_FULL_IMAGE:figures/full_fig_p030_24.png]
Figure 25
Figure 25. Figure 25: Topology-optimized specimen geometry used for inverse identification of Barlat’s Yld2004-18p model, including the applied boundary conditions and loading direction. The specimen thickness is set to 1.0 mm. All dimensions are in mm [PITH_FULL_IMAGE:figures/full_fig_p0…
Figure 26
Figure 26. Figure 26: (a) Equivalent plastic-strain distribution with the deformed geometry at the final time step, (b) principal strain-state diagram at the final time step, and (c) the corresponding principal stress-state diagram with yield surfaces overlaid. Only material points undergo…
Figure 27
Figure 27. Figure 27: Inverse identification of Barlat’s model: (a) comparison of reference and identified yield surfaces in the principal-stress plane for two hardening states (α = 0 and α = 0.2), and (b) comparison of reference and identified hardening curves [PITH_FULL_IMAGE:figures/fu…
Figure 28
Figure 28. Figure 28: Inverse optimization for Barlat material parameter identification: (a) convergence history of the optimization, and (b) comparison of forward-pass simulation times with and without warm-start initialization. 4 Conclusions This work presented an end-to-end differentiab…

Review history (2 revisions) →

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

68 extracted references · 1 linked inside Pith

  1. [1]

    J. C. Simo and T. J. R. Hughes.Computational Inelasticity. Springer, New York, NY , 1998

    1998

  2. [2]

    Belytschko, W

    T. Belytschko, W. K. Liu, and B. Moran.Nonlinear Finite Elements for Continua and Structures. John Wiley & Sons, Chichester, UK, 2000

    2000

  3. [3]

    R. Hill. A theory of the yielding and plastic flow of anisotropic metals.Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences, 193(1033):281–297, 1948

    1948

  4. [4]

    Barlat, J

    F. Barlat, J. C. Brem, J. W. Yoon, K. Chung, R. E. Dick, D. J. Lege, F. Pourboghrat, S.-H. Choi, and E. Chu. Plane stress yield function for aluminum alloy sheets—part 1: theory.International Journal of Plasticity, 19(9):1297–1319, 2003

    2003

  5. [5]

    Barlat, H

    F. Barlat, H. Aretz, J. W. Yoon, M. E. Karabin, J. C. Brem, and R. E. Dick. Linear transformation-based anisotropic yield functions.International Journal of Plasticity, 21(5):1009–1039, 2005

    2005

  6. [6]

    E. A. de Souza Neto, D. Peri´c, and D. R. J. Owen.Computational Methods for Plasticity: Theory and Applications. John Wiley & Sons, Chichester, UK, 2008

    2008

  7. [7]

    Michaleris, D

    P. Michaleris, D. A. Tortorelli, and C. A. Vidal. Tangent operators and design sensitivity formulations for transient non- linear coupled problems with applications to elastoplasticity.International Journal for Numerical Methods in Engineering, 37(14):2471–2499, 1994

    1994

  8. [8]

    Miehe, N

    C. Miehe, N. Apel, and M. Lambrecht. Anisotropic additive plasticity in the logarithmic strain space: modular kinematic formulation and implementation based on incremental minimization principles for standard materials.Computer Methods in Applied Mechanics and Engineering, 191:5383–5425, 2002

    2002

Show all 68 references
  1. [9]

    N. Aravas. On the numerical integration of a class of pressure-dependent plasticity models.International Journal for Numerical Methods in Engineering, 24(7):1395–1416, 1987

    1987

  2. [10]

    Güner, C

    A. Güner, C. Soyarslan, A. Brosius, and A. E. Tekkaya. Characterization of anisotropy of sheet metals employing inhomogeneous strain fields for Yld2000-2D yield function.International Journal of Solids and Structures, 49(25):3517– 3527, 2012

    2012

  3. [11]

    J. H. Kim, F. Barlat, F. Pierron, and M. G. Lee. Determination of anisotropic plastic constitutive parameters using the virtual fields method.Experimental Mechanics, 54(7):1189–1204, 2014

    2014

  4. [12]

    M. A. Sutton, J.-J. Orteu, and H. Schreier.Image Correlation for Shape, Motion and Deformation Measurements: Basic Concepts, Theory and Applications. Springer, New York, NY , 2009

    2009

  5. [13]

    Bertin, F

    M. Bertin, F. Hild, and S. Roux. On the identifiability of Hill-1948 plasticity model with a single biaxial test on very thin sheet.Strain, 53(5):e12233, 2017

    1948

  6. [14]

    Coppieters, T

    S. Coppieters, T. Hakoyama, D. Debruyne, S. Takahashi, and T. Kuwabara. Inverse yield locus identification of sheet metal using a complex cruciform in biaxial tension and digital image correlation.Proceedings, 60(2):382, 2018

    2018

  7. [15]

    Zhang, M

    Y . Zhang, M. Rossi, R. Pettermann, and C. Sommitsch. Inverse identification of plastic anisotropy through multiple non-conventional mechanical experiments.International Journal of Solids and Structures, 285:112534, 2023

    2023

  8. [16]

    Pierron and M

    F. Pierron and M. Grédiac. Towards material testing 2.0. a review of test design for identification of constitutive parameters from full-field measurements.Strain, 57(1):e12370, 2021

    2021

  9. [17]

    J. M. P. Martins, A. Andrade-Campos, and S. Thuillier. Comparison of inverse identification strategies for constitutive mechanical models using full-field measurements.International Journal of Mechanical Sciences, 145:330–345, 2018

    2018

  10. [18]

    Marek, F

    A. Marek, F. M. Davis, M. Rossi, and F. Pierron. Extension of the sensitivity-based virtual fields to large deformation anisotropic plasticity.International Journal of Material Forming, 12(3):457–476, 2019

    2019

  11. [19]

    Zhang, M

    Y . Zhang, M. Rossi, J. Reisch, R. Pettermann, and C. Sommitsch. Enhancing the information-richness of sheet metal specimens for inverse identification of plastic anisotropy through strain fields.International Journal of Mechanical Sciences, 214:106891, 2022

    2022

  12. [20]

    Conde, Y

    M. Conde, Y . Zhang, J. Henriques, S. Coppieters, and A. Andrade-Campos. Design and validation of a heterogeneous interior notched specimen for inverse material parameter identification.Finite Elements in Analysis and Design, 214:103866, 2023

    2023

  13. [21]

    Gonçalves, S

    M. Gonçalves, S. Thuillier, and A. Andrade-Campos. Inverse identification of anisotropic plasticity model parameters using FEMU and a heterogeneous test. InESAFORM 2025, volume 54, pages 1548–1557, 2025

    2025

  14. [22]

    D. T. Seidl and B. N. Granzow. Calibration of elastoplastic constitutive model parameters from full-field data with automatic differentiation-based sensitivities.International Journal for Numerical Methods in Engineering, 123(1):69–100, 2022

    2022

  15. [23]

    Kumar, D

    S. Kumar, D. T. Seidl, B. N. Granzow, J. Yang, and J. N. Fuhg. A comparative study of calibration techniques for finite strain elastoplasticity: Numerically-exact sensitivities for FEMU and VFM.Computer Methods in Applied Mechanics and Engineering, 444:118159, 2025

    2025

  16. [24]

    C. C. Margossian. A review of automatic differentiation and its efficient implementation.Wiley Interdisciplinary Reviews: Data Mining and Knowledge Discovery, 9(4):e1305, 2019

    2019

  17. [25]

    Rothe and S

    S. Rothe and S. Hartmann. Automatic differentiation for stress and consistent tangent computation.Archive of Applied Mechanics, 85(8):1103–1125, 2015

    2015

  18. [26]

    Dummer, M

    A. Dummer, M. Neuner, P. Gamnitzer, and G. Hofstetter. Robust and efficient implementation of finite strain generalized continuum models for material failure: Analytical, numerical, and automatic differentiation with hyper-dual numbers. Computer Methods in Applied Mechanics an...

    2024

  19. [27]

    Y . Jia, W. Li, and X. S. Zhang. Multimaterial topology optimization of elastoplastic composite structures.Journal of the Mechanics and Physics of Solids, 196:106018, 2025

    2025

  20. [28]

    Jia and X

    Y . Jia and X. S. Zhang. Multimaterial topology optimization for finite strain elastoplasticity: theory, methods, and applications. Computer Methods in Applied Mechanics and Engineering, 449:118445, 2026

    2026

  21. [29]

    Bradbury, R

    J. Bradbury, R. Frostig, P. Hawkins, M. J. Johnson, C. Leary, D. Maclaurin, G. Necula, A. Paszke, J. Vander- Plas, S. Wanderman-Milne, and Q. Zhang. JAX: composable transformations of Python+NumPy programs. http: //github.com/google/jax, 2018

    2018

  22. [30]

    Paszke, S

    A. Paszke, S. Gross, F. Massa, A. Lerer, J. Bradbury, G. Chanan, et al. PyTorch: An imperative style, high-performance deep learning library.Advances in Neural Information Processing Systems, 32, 2019

    2019

  23. [31]

    M. Macklin. Warp: A high-performance python framework for GPU simulation and graphics. InNVIDIA GPU Technology Conference (GTC), volume 3, 2022

    2022

  24. [32]

    NVIDIA H100 Tensor Core GPU datasheet

    NVIDIA Corporation. NVIDIA H100 Tensor Core GPU datasheet. https://resources.nvidia.com/ en-us-hopper-architecture/nvidia-tensor-core-gpu-datasheet, 2023

    2023

  25. [33]

    Intel Xeon 6952P Processor specifications

    Intel Corporation. Intel Xeon 6952P Processor specifications. https://www.intel.com/content/www/us/en/ products/sku/241643/intel-xeon-6952p-processor-480m-cache-2-10-ghz/specifications.html, 2026

    2026

  26. [34]

    Williams, A

    S. Williams, A. Waterman, and D. Patterson. Roofline: an insightful visual performance model for multicore architectures. Communications of the ACM, 52(4):65–76, 2009

    2009

  27. [35]

    Martínez-Frutos, P

    J. Martínez-Frutos, P. J. Martínez-Castejón, and D. Herrero-Pérez. Fine-grained GPU implementation of assembly-free iterative solver for finite element problems.Computers&Structures, 157:9–18, 2015

    2015

  28. [36]

    Aissa, T

    M. Aissa, T. Verstraete, and C. Vuik. Toward a GPU-aware comparison of explicit and implicit CFD simulations on structured meshes.Computers&Mathematics with Applications, 74(1):201–217, 2017

    2017

  29. [37]

    Stavroulakis, D

    G. Stavroulakis, D. G. Giovanis, V . Papadopoulos, and M. Papadrakakis. A GPU domain decomposition solution for spectral stochastic finite element method.Computer Methods in Applied Mechanics and Engineering, 327:392–410, 2017

    2017

  30. [38]

    Ghysels and R

    P. Ghysels and R. Synk. High performance sparse multifrontal solvers on modern GPUs.Parallel Computing, 110:102897, 2022

    2022

  31. [39]

    F. Li, F. Zou, and J. Rao. A multi-GPU and CUDA-aware MPI-based spectral element formulation for ultrasonic wave propagation in solid media.Ultrasonics, 134:107049, 2023

    2023

  32. [40]

    Kiran, D

    U. Kiran, D. Sharma, and S. S. Gautam. Development of GPU-based matrix-free strategies for large-scale elastoplasticity analysis using conjugate gradient solver.International Journal for Numerical Methods in Engineering, 125(7):e7421, 2024

    2024

  33. [41]

    S. Hong, G. Jang, and W. K. Jeong. MG-FIM: A multi-GPU fast iterative method using adaptive domain decomposition. SIAM Journal on Scientific Computing, 44(1):C54–C76, 2022

    2022

  34. [42]

    S. Liao, A. Golgoon, M. Mozaffar, and J. Cao. Efficient GPU-accelerated thermomechanical solver for residual stress prediction in additive manufacturing.Computational Mechanics, 71(5):879–893, 2023

    2023

  35. [43]

    E. A. Träff, A. Rydahl, S. Karlsson, O. Sigmund, and N. Aage. Simple and efficient GPU accelerated topology optimisation: Codes and applications.Computer Methods in Applied Mechanics and Engineering, 410:116043, 2023

    2023

  36. [44]

    Kiran, D

    U. Kiran, D. Sharma, and S. S. Gautam. An efficient framework for matrix-free SpMV computation on GPU for elastoplastic problems.Mathematics and Computers in Simulation, 216:318–346, 2024

    2024

  37. [45]

    Karatarakis, P

    A. Karatarakis, P. Karakitsios, and M. Papadrakakis. GPU accelerated computation of the isogeometric analysis stiffness matrix.Computer Methods in Applied Mechanics and Engineering, 269:334–355, 2014

    2014

  38. [46]

    T. Xue, S. Liao, Z. Gan, C. Park, X. Xie, W. K. Liu, and J. Cao. JAX-FEM: A differentiable GPU-accelerated 3d finite element solver for automatic inverse design and mechanistic data science.Computer Physics Communications, 291:108802, 2023

    2023

  39. [47]

    F. Hu, S. Niezgoda, T. Xue, and J. Cao. Efficient GPU-computing simulation platform JAX-CPFEM for differentiable crystal plasticity finite element method.npj Computational Materials, 11(1):46, 2025

    2025

  40. [48]

    F. Hu, R. Zhou, K. Ryou, R. Zha, S. Niezgoda, T. Xue, and J. Cao. An efficient graphical processing unit-accelerated calibration of crystal plasticity model parameters by multi-objective optimization with automatic differentiation-based sensitivities.Journal of Applied Mechani...

    2026

  41. [49]

    B. P. Ferreira and M. A. Bessa. Automatically differentiable model updating (ADiMU): conventional, hybrid, and neural network material model discovery including history-dependency.arXiv preprint arXiv:2505.07801, 2025

    2025 arXiv

  42. [50]

    Knapik, S

    S. Knapik, S. Deng, L. Wang, A. Chandrasekhar, W. K. Liu, J. Cao, and W. Chen. Co-design of geometry and thermal-elastic gradient alloy distribution with temperature-dependent material properties.Structural and Multidisciplinary Optimization, 68(7):127, 2025

    2025

  43. [51]

    J. P. Leonor and G. J. Wagner. GO-MELT: GPU-optimized multilevel execution of LPBF thermal simulations.Computer Methods in Applied Mechanics and Engineering, 426:116977, 2024

    2024

  44. [52]

    AmgX: GPU accelerated linear solvers.https://developer.nvidia.com/amgx, 2026

    NVIDIA Corporation. AmgX: GPU accelerated linear solvers.https://developer.nvidia.com/amgx, 2026

    2026

  45. [53]

    E. A. de Souza Neto, D. Peric, and D. R. Owen.Computational Methods for Plasticity: Theory and Applications. John Wiley & Sons, Chichester, UK, 2011

    2011

  46. [54]

    W. M. Scherzinger. A return mapping algorithm for isotropic and anisotropic plasticity models using a line search method. Computer Methods in Applied Mechanics and Engineering, 317:526–553, 2017

    2017

  47. [55]

    Scalet and F

    G. Scalet and F. Auricchio. Computational methods for elastoplasticity: An overview of conventional and less-conventional approaches.Archives of Computational Methods in Engineering, 25(3):545–589, 2018

    2018

  48. [56]

    Alfano, L

    G. Alfano, L. Rosati, and N. Valoroso. A tangent–secant approach to rate-independent elastoplasticity: formulations and computational issues.Computer Methods in Applied Mechanics and Engineering, 179(3–4):379–405, 1999

    1999

  49. [57]

    Miehe and M

    C. Miehe and M. Lambrecht. Algorithms for computation of stresses and elasticity moduli in terms of Seth–Hill’s family of generalized strain tensors.Communications in Numerical Methods in Engineering, 17(5):337–353, 2001

    2001

  50. [58]

    Z. Fu, T. J. Lewis, R. M. Kirby, and R. T. Whitaker. Architecting the finite element method pipeline for the GPU.Journal of Computational and Applied Mathematics, 257:195–211, 2014

    2014

  51. [59]

    Meng, B.-L

    H.-T. Meng, B.-L. Nie, S. Wong, C. Macon, and J.-M. Jin. GPU accelerated finite-element computation for electromagnetic analysis.IEEE Antennas and Propagation Magazine, 56(2), 2014

    2014

  52. [60]

    Kiran, D

    U. Kiran, D. Sharma, and S. S. Gautam. GPU-warp based finite element matrices generation and assembly using coloring method.Journal of Computational Design and Engineering, 6(5):705–718, 2019

    2019

  53. [61]

    M. B. Giles and N. A. Pierce. An introduction to the adjoint approach to design.Flow, Turbulence and Combustion, 65:393–415, 2000

    2000

  54. [62]

    Xu and E

    K. Xu and E. Darve. Physics constrained learning for data-driven inverse modeling from sparse observations.Journal of Computational Physics, 453:110938, 2022

    2022

  55. [63]

    Papadopoulos and J

    P. Papadopoulos and J. Lu. A general framework for the numerical solution of problems in finite elasto-plasticity.Computer Methods in Applied Mechanics and Engineering, 159(1–2):1–18, 1998

    1998

  56. [64]

    unified material model driver for plasticity (UMMDp)

    H. Takizawa, K. Oide, K. Suzuki, T. Yamanashi, T. Inoue, T. Ida, T. Nagai, and T. Kuwabara. Development of the user subroutine library “unified material model driver for plasticity (UMMDp)” for various anisotropic yield functions.Journal of Physics: Conference Series, 1063:012...

    2018

  57. [65]

    Tammas-Williams and I

    S. Tammas-Williams and I. Todd. Design for additive manufacturing with site-specific properties in metals and alloys. Scripta Materialia, 135:105–110, 2017

    2017

  58. [66]

    Gonçalves, A

    M. Gonçalves, A. Andrade-Campos, and B. Barroqueiro. On the design of mechanical heterogeneous specimens using multilevel topology optimization.Advances in Engineering Software, 175:103314, 2023

    2023

  59. [67]

    Gonçalves, M

    M. Gonçalves, M. G. Oliveira, S. Thuillier, and A. Andrade-Campos. Key performance indicators for heterogeneous mechanical tests.International Journal of Mechanical Sciences, 264:108821, 2024

    2024

  60. [68]

    Efficient and modular implicit differentiation.Advances in neural information processing systems, 35:5230–5242, 2022

    Mathieu Blondel, Quentin Berthet, Marco Cuturi, Roy Frostig, Stephan Hoyer, Felipe Llinares-López, Fabian Pedregosa, and Jean-Philippe Vert. Efficient and modular implicit differentiation.Advances in neural information processing systems, 35:5230–5242, 2022

    2022

Pith tools

Reviewed July 14, 2026 · model on record in the stance chip above.