5 Johnson–Cook parameters. 3 strain rate regimes. 3 specimens instead of 12. Simultaneous identification from full DIC fields — with validated uncertainty quantification and independent test confirmation.
Inconel 625 is widely used in aerospace, nuclear, and defence applications — turbine casings, pressure vessels, high-velocity impact shields — where the material must be proven safe across a wide range of temperatures, strain rates, and stress states. Reliable simulation of these components depends entirely on the quality of the constitutive material card.
Classical calibration builds that card from a suite of uniaxial specimens: smooth tensile for A, B, n — then separate high-rate compression tests for C — then heated specimens for m. Each test is post-processed independently, the stress–strain curve is curve-fitted, and parameters are handed off sequentially. Each step propagates the errors of the previous one.
The core problem: every classical test delivers a single global scalar — one force–displacement curve. The stress and strain fields inside the specimen are assumed uniform and never measured. For a material like Inconel 625 where the hardening, rate sensitivity, and failure mode interact non-trivially, this assumption is a source of systematic bias that no amount of additional specimens resolves.
SIMFORGE applied the Virtual Fields Method (VFM) — a full-field inverse technique grounded in the Principle of Virtual Work — to identify all five Johnson–Cook parameters simultaneously from DIC strain fields acquired during a small number of purpose-designed heterogeneous tests.
Rather than assuming a uniform stress state, the method exploits the spatial gradient of the strain field as a rich mechanical signal. Virtual fields are constructed analytically to be maximally sensitive to each target parameter, directly conditioning the identification problem.
Specimen geometry was not arbitrary. A sensitivity-driven FE pre-study was used to select notch geometry such that the condition number of the identification system was minimised — guaranteeing parameter decoupling and numerical robustness before a single test was run.
| Parameter | Physical meaning | Regime |
|---|---|---|
| A | Initial yield stress | Quasi-static |
| B, n | Strain hardening coefficient & exponent | Quasi-static |
| C | Strain rate sensitivity | Intermediate + SHPB |
| m | Thermal softening exponent | High-temp campaign |
| ε̇₀ | Reference strain rate (fixed) | 10⁻³ s⁻¹ baseline |
VFM is built on the weak form of the equilibrium equation — the Principle of Virtual Work. In standard FEA, this principle is used to solve for displacements given a constitutive model. VFM inverts the problem: given measured displacements and strains (from DIC), it solves for the constitutive parameters.
For the Johnson–Cook model, the VFM cost function assembles a linear system in the parameter space by choosing virtual fields u* that are orthogonal to each parameter's sensitivity field. This yields a closed-form identification equation — no iterative FE model updating required.
The sensitivity-based virtual field approach requires that the strain field span the full parameter sensitivity space simultaneously. A smooth dog-bone specimen fails this requirement: its gauge section is nearly uniform, meaning H is near-singular and parameter cross-talk is severe.
For Inconel 625, a notched flat specimen with a 4 mm notch radius and 3 mm minimum width was selected after iterating on the FE sensitivity analysis. The objective function was the condition number κ(H) — minimised across candidate geometries:
The optimised geometry produces local equivalent strains up to ε̄ = 0.44 at the notch root before fracture — sufficient to activate the full hardening range — while maintaining quasi-uniform strain in the far-field gauge section for reference state comparison.
| Minimum width | 3 mm |
| Notch radius | 4 mm |
| Gauge length | 20 mm |
| Peak ε̄ (notch root) | 0.44 |
| κ(H) reduction vs smooth | 3.1× |
The full JC identification is structured across three strain rate regimes, each requiring a different experimental platform. The VFM cost function and parameter extraction engine are identical across all three — only the inertia treatment changes.
The 1–100 s⁻¹ window is systematically under-characterised in the literature: too fast for quasi-static machines, too slow to be treated as SHPB wave propagation. Classical approaches either ignore it or interpolate between two extremes. VFM treats it natively — the dynamic balance equation includes the inertia term explicitly from measured acceleration fields, making this the highest-fidelity identification window for the rate sensitivity parameter C.
The following parameter set was identified via VFM across the three strain rate regimes. Confidence intervals are derived from a Monte Carlo uncertainty propagation on the DIC displacement noise floor (σ_u ≈ 0.02 px RMS) — providing traceable bounds, not just point estimates.
| Parameter | Symbol | VFM Result | Literature | Regime |
|---|---|---|---|---|
| Yield stress | A | 485 ± 11 MPa | 450–520 | QS |
| Hardening coeff. | B | 612 ± 22 MPa | 580–650 | QS |
| Hardening exp. | n | 0.41 ± 0.02 | 0.38–0.45 | QS |
| Rate sensitivity | C | 0.013 ± 0.002 | 0.010–0.017 | Int. + SHPB |
| Thermal softening | m | 1.24 ± 0.08 | 1.10–1.40 | High-T |
| Reference rate | ε̇₀ | 10⁻³ s⁻¹ | — | Fixed |
Charpy-type notched impact specimen at 5 m/s impact velocity. FE simulation using VFM-identified JC parameters reproduced peak force within 4.2% and energy within 3.8%. Classical calibration on the same batch: 11% peak force error on the same validation geometry.
The performance advantage of VFM over classical calibration is not marginal — it is structural. Classical calibration extracts one scalar per test. VFM extracts 10,000–50,000 spatial data points per image. This fundamentally changes the information density of the identification problem.
For rare or expensive nickel alloys — single-crystal, ODS, AM-deposited — the 3× reduction in specimen count alone can justify the VFM approach on material cost. For high-rate aerospace qualification programmes, the combined time and accuracy improvement is the differentiator.
Classical identification of C relies on comparing two independently obtained stress–strain curves at different rates. Any systematic error in either test (friction, wave dispersion, alignment) transfers directly into C. VFM identifies C from the spatial strain rate gradient within a single specimen — the systematic errors cancel. For Inconel 625, the VFM-derived C was 18% lower than the classical estimate from the same batch.
The 1–100 s⁻¹ range is where inertia effects become measurable but wave mechanics do not dominate. VFM's explicit inertia term makes this regime natively tractable. The dynamic VFM balance provides richer constraint on both B/n and C simultaneously — in a single test — than either a quasi-static suite or a SHPB set provides individually.
The VFM system H·p = b is linear in the parameters (for the isothermal JC case). This means a Monte Carlo propagation of DIC noise through the identification system yields full confidence intervals on every parameter simultaneously — in minutes, not days. Classical methods typically report no UQ, or derive it from repeat specimen scatter, which conflates identification error with material variability.
Classical standards (ISO 6892, ASTM E8) prescribe smooth geometry specifically to achieve uniform stress — which is exactly what VFM does not want. The notched geometry used here is not a compromise; it is a precision instrument. The sensitivity-based FE pre-study that selects the notch radius is as important as the test itself — and it takes one morning, not a separate test campaign.
| Specimen design | FE sensitivity analysis to minimise κ(H) — notch geometry selected for parameter decoupling across all five JC parameters |
| Test programme | 3–4 specimens across quasi-static, intermediate (servo-hydraulic), and SHPB regimes — DIC acquisition at each |
| VFM identification | Full-field PVW cost function · sensitivity-based virtual fields · simultaneous 4-parameter isothermal identification + thermal softening m from separate high-T campaign |
| Uncertainty quantification | Monte Carlo propagation of DIC displacement noise (σ_u = 0.02 px) through H·p = b — 95% CI on all parameters |
| Validation | Independent Charpy-type impact test at 5 m/s · FE simulation vs experimental force-time history · 4.2% peak force error |
| Deliverable | JC material card for LS-DYNA / Abaqus / Nastran SOL 700 · parameter report with sensitivity maps · full identification traceability |
VFM is not a bolt-on service — it requires simultaneous expertise in inverse methods, DIC acquisition, high-rate mechanics, and uncertainty quantification. SIMFORGE combines all four in a single engagement. The deliverable is not just a number — it is a parameter set with traceable confidence bounds, validated against an independent test, ready to drop into your simulation solver.
Whether it's a nickel superalloy for bird-strike simulation, a titanium structure under ballistic loading, or a CFRP layup for crash analysis — if it needs a high-fidelity constitutive card with traceable uncertainty, we scope it in a 30-minute call and deliver a fixed-fee proposal within 48 hours.
Initial response · 24h · NDA available · Fixed-fee proposal within 48h