Development of an inverse identification method for identifying constitutive parameters by metaheuristic optimization algorithm: Application to hyperelastic materials
Pith reviewed 2026-05-25 02:06 UTC · model grok-4.3
The pith
An adapted particle swarm optimizer identifies hyperelastic constitutive parameters by matching full kinematic fields from one cruciform test.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The constitutive parameters of hyperelastic models can be identified by minimizing the difference between experimental and numerical kinematic fields using an adaptation of the Particle Swarm Optimization algorithm based on the PageRank algorithm. The experimental field comes from digital image correlation on the surface of a cruciform specimen loaded in equibiaxial tension; the numerical field is generated by a finite-element reconstruction of the same specimen and loading. The optimization loop updates the parameters inside the hyperelastic constitutive law until the two fields agree.
What carries the argument
An adaptation of the Particle Swarm Optimization algorithm based on the PageRank algorithm that drives iterative adjustment of constitutive parameters to reduce mismatch between measured and simulated kinematic fields.
If this is right
- Constitutive parameters are obtained from a single heterogeneous test instead of three separate homogeneous tests.
- Full-field data from digital image correlation supplies the large number of mechanical states needed for identification.
- The cost function directly penalizes differences in the entire measured displacement field rather than averaged stress-strain points.
- The PageRank-based swarm update rule guides the search through the parameter space without requiring gradient information.
Where Pith is reading between the lines
- The method could be tested on other nonlinear constitutive models by swapping only the material subroutine inside the same finite-element reconstruction.
- If measurement noise is spatially correlated, the cost function might need an additional weighting step before the swarm can converge reliably.
- Running the optimization on data collected during the test itself could support closed-loop control of the loading to improve parameter identifiability.
Load-bearing premise
The finite-element model is assumed to faithfully reproduce the experimental boundary conditions, geometry, and measurement noise so that any remaining mismatch is due only to incorrect constitutive parameters.
What would settle it
Apply the parameters identified from the cruciform test to predict the response of a different specimen geometry or loading condition and compare the predicted surface displacements against new experimental measurements.
read the original abstract
In the present study, a numerical method based on a metaheuristic parametric algorithm has been developed to identify the constitutive parameters of hyperelastic models, by using FE simulations and full kinematic field measurements. The full kinematic field is measured at the surface of a cruciform specimen submitted to equibiaxial tension. The sample is reconstructed by FE to obtain the numerical kinematic field to be compared with the experimental one. The constitutive parameters used in the numerical model are then modified through the optimization process, for the numerical kinematic field to fit with the experimental one. The cost function is then formulated as the minimization of the difference between these two kinematic fields. The optimization algorithm is an adaptation of the Particle Swarm Optimization algorithm, based on the PageRank algorithm used by the famous search engine Google. INTRODUCTION The constitutive parameters of hyperelastic models are generally identified from three homogeneous tests, basically the uniaxial tension, the pure shear and the equibiaxial tension. From about 10 years, an alternative methodology has been developed [1, 2, 3, 4], and consists in performing only one heterogeneous test as long as the field is sufficiently heterogeneous. This is tipically the case when a multiaxial loading is applied to a 3 branch or a 4-branch cruciform specimen, which induces a large number of mechanical states at the specimen surface. The induced heterogeneity is generally analysed through the distribution of the biaxiality ratio and the maximal eigen value of the strain. The Digital Image Correlation (DIC) technique is generally used to retrieve the different mechanical states induces, and provides the full kinematic field at the specimen surface, i.e. a large number of experimental data to be analysed to identify the constitutive parameters of the behaviour model considered.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims to develop a numerical inverse identification method for hyperelastic constitutive parameters. It uses full kinematic fields from DIC on a cruciform specimen under equibiaxial tension, compares them to FE-simulated fields, and minimizes the difference via an adapted Particle Swarm Optimization algorithm based on the PageRank algorithm. The cost function is the kinematic-field discrepancy, with the method positioned as an alternative to multiple homogeneous tests.
Significance. If the central claim holds, the work would contribute an algorithmic variant of PSO for inverse problems in hyperelasticity and demonstrate parameter recovery from a single heterogeneous test. The approach leverages existing DIC and FE tools, but its significance depends on demonstrated accuracy rather than the workflow description alone.
major comments (2)
- Abstract: the workflow is described but supplies no quantitative validation, convergence data, comparison against known parameter sets, or error metrics; without these it is impossible to verify whether the optimization recovers the intended parameters.
- Introduction (paragraph on FE reconstruction and cost-function formulation): the claim that parameter identification follows from minimizing the kinematic discrepancy assumes the FE model reproduces experimental boundary conditions, geometry, and measurement noise to a precision finer than the cost-function sensitivity to the constitutive parameters. No error-propagation study, synthetic-data recovery test with controlled perturbations, or residual analysis is provided to support this assumption for the chosen cruciform geometry and DIC data.
minor comments (1)
- Introduction: 'tipically' is a typo and should read 'typically'.
Simulated Author's Rebuttal
We thank the referee for the detailed and constructive report. The two major comments correctly identify gaps in quantitative validation and supporting analyses. We will revise the manuscript to address both points directly.
read point-by-point responses
-
Referee: Abstract: the workflow is described but supplies no quantitative validation, convergence data, comparison against known parameter sets, or error metrics; without these it is impossible to verify whether the optimization recovers the intended parameters.
Authors: We agree that the abstract lacks quantitative results. The revised abstract will report the recovered constitutive parameters, the final cost-function value, convergence iterations, and comparison against reference values obtained from homogeneous tests. revision: yes
-
Referee: Introduction (paragraph on FE reconstruction and cost-function formulation): the claim that parameter identification follows from minimizing the kinematic discrepancy assumes the FE model reproduces experimental boundary conditions, geometry, and measurement noise to a precision finer than the cost-function sensitivity to the constitutive parameters. No error-propagation study, synthetic-data recovery test with controlled perturbations, or residual analysis is provided to support this assumption for the chosen cruciform geometry and DIC data.
Authors: We accept the criticism. The revised manuscript will add a dedicated validation subsection containing (i) synthetic-data recovery tests with added Gaussian noise matching the DIC uncertainty and (ii) a brief error-propagation analysis showing that the chosen cost function remains sensitive to the constitutive parameters under realistic perturbations of boundary conditions and geometry. revision: yes
Circularity Check
No significant circularity; standard external optimization loop around independent data.
full rationale
The paper presents an inverse identification procedure that minimizes a cost function defined as the difference between experimental DIC kinematic fields and FE-simulated fields, with parameters updated via an adapted PSO algorithm. This constitutes an external numerical optimization around independently measured data and standard FE modeling; the cost function and update rule are not defined in terms of the target parameters themselves. No self-definitional steps, fitted inputs renamed as predictions, or load-bearing self-citations appear in the derivation. The central claim remains an independent algorithmic method rather than a tautology.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The finite-element reconstruction accurately captures the experimental boundary conditions and measurement noise.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.