Start here

What This Field Is About

Two diagrams that place the work in context, for readers coming from outside impact mechanics. Both are taken from the PhD thesis behind the constitutive modeling stream.

Diagram showing the three parts of a numerical simulation: geometric representation, constitutive equation, and physical principles, with the flow stress equation highlighted inside the constitutive branch.
Figure 1. The three parts of a numerical simulation of a metal structure under impact, with the flow stress equation highlighted. Reproduced from the author's PhD thesis, Fig. 1-1.

Where the uncertainty actually sits

Any simulation of a metal structure under impact rests on three things: a geometric representation of the part, the physical conservation laws, and a constitutive equation describing how the material itself responds. The first two are mature and largely settled. The third is not.

Inside the constitutive branch, the shaded box is the flow stress equation, which sets the stress at which the metal yields and keeps deforming as a function of strain, strain rate, and temperature. In practice it is the largest single error source in an otherwise well-posed impact simulation. That is the reason this equation, rather than the solver or the mesh, is the subject of the constitutive modeling work on this site.

Diagram comparing three routes from characterized material data to a flow stress equation: the conventional fitted-equation route, the machine learning route, and the discrete-data route with data screening and SVD/CP decomposition.
Figure 2. Three routes from measured material data to a usable flow stress equation. The shaded branch on the right is the route developed in the thesis. Reproduced from the author's PhD thesis, Fig. 1-3.

Three ways to get from data to an equation

The conventional route on the left picks a fixed formula in advance, the Johnson-Cook model being the usual choice, and fits its constants to the measurements. It is fast and every solver accepts it, but the form of the equation is decided before the data is looked at, so whatever the material does that the formula cannot express is lost.

The machine learning route in the middle drops the fixed form and trains a network directly on the data. It fits well inside the tested range and gives little confidence outside it, which is exactly where impact problems live.

The third route keeps the discrete data as the primary object. Measurements are first screened for validity, then SVD or CP decomposition is used to separate the underlying structure from the noise, and an artificial neural network handles the non-ideal data whose strain rate was never really constant. This is the line of work that became TEDI-FS, described on the constitutive modeling page.

Source: Xianglin Huang, Determination of dynamic flow stress equation based on discrete material data, PhD thesis, Department of Mechanical and Aerospace Engineering, The University of Manchester, 2024. Supervisor: Prof. Q.M. Li.

Long-form

Technical Notes

Self-contained write-ups of the mathematics behind specific methods, written to be readable on their own rather than as supplements to a particular paper.

Linear Algebra Technical Note
Singular Value Decomposition of a Matrix
The SVD from the ground up: the factorization and its singular values and vectors, a constructive existence proof from the eigendecomposition of the Gram matrix, the four fundamental subspaces, the Eckart-Young-Mirsky optimal low-rank approximation theorem, the connection to PCA, and a fully worked 2x2 example, closing with the rank-1 legitimacy criterion for decoupled flow stress equations.
8 sections  ·  7 equations  ·  Aug 2026
Tensor Decomposition Technical Note
CP (CANDECOMP/PARAFAC) Decomposition of Order-3 Tensors
The canonical polyadic model as the tensor generalization of SVD: the rank-1 outer-product model, tensor rank, matricization and the Khatri-Rao product, the alternating least squares (CP-ALS) algorithm with its normal equations and MTTKRP kernel, rank selection, Kruskal's uniqueness condition, and a hand-verifiable worked example, with the relation to the non-negative variant used in the flow stress work.
8 sections  ·  9 equations  ·  Aug 2026
Tensor Decomposition Technical Note
Non-negative CP Decomposition of Order-3 Tensors
The mathematical principles of non-negative CP (CANDECOMP/PARAFAC) decomposition: the outer-product tensor model, matricization and the Khatri-Rao product, the non-negative least squares problem and its KKT conditions, the alternating multiplicative update algorithm, Kruskal's uniqueness condition, and worked Rank-1 to Rank-3 numerical examples that can be checked by hand.
8 sections  ·  15 equations  ·  Jul 2026
Dynamic Testing Technical Note
The Split Hopkinson Bar (SHPB)
How the split Hopkinson pressure bar measures material behaviour at 10^2 to 10^4 per second: the apparatus and one-dimensional wave principle, the two-wave and three-wave data-reduction equations, the stress-equilibrium and constant-strain-rate validity conditions, pulse shaping and geometric dispersion, and the tension and torsion variants of the technique.
8 sections  ·  4 equations  ·  diagram  ·  Aug 2026
Constitutive Models Reference Column
Common Constitutive Models in Impact Simulation
A reference column on the models used most often in crash, impact, and penetration analysis: the Johnson-Cook flow stress and fracture models and their Cowper-Symonds and Zerilli-Armstrong relatives for metals, the pressure-dependent Drucker-Prager model, the Deshpande-Fleck crushable foam model for cellular solids, and the Johnson-Holmquist JH-2 model for brittle ceramics, each with its yield surface, governing equations, and primary references.
8 sections  ·  8 equations  ·  Aug 2026
Machine Learning Reference Column
Neural Network Architectures: ANN, CNN, LSTM, Transformer
A detailed column on the four architectures used most often in engineering and signal analysis: the fully connected ANN/MLP with the backpropagation derivation, the CNN with the convolution and pooling operations, the LSTM recurrent network with its gate equations, and the Transformer with the scaled dot-product self-attention derivation, each with a schematic diagram and primary references.
8 sections  ·  9 equations  ·  4 diagrams  ·  Aug 2026
In Preparation
More notes to follow
Further notes are planned on the data-screening rules that precede any flow stress fitting, on wave dispersion correction in Hopkinson bar signals, and on acoustic emission feature extraction for battery health monitoring. They will be added here as they are written.
Planned
Terminology

Glossary

Terms that appear throughout this site, defined in a sentence or two. Highlighted terms in the page text link directly to the corresponding entry here.

CONSTITUTIVE MODELING AND DYNAMIC TESTING

Flow stress
also written as the flow stress equation, or constitutive relation

The stress a metal sustains while it is plastically deforming, as a function of the accumulated plastic strain, the strain rate, and the temperature. It is the single material input that governs how a finite element solver such as ABAQUS or LS-DYNA predicts a crash, impact, or perforation event, so the accuracy of the flow stress description places a ceiling on the accuracy of the whole simulation.

Determining it is difficult because the three variables cannot be controlled independently in a single experiment: a high strain rate test heats the specimen through plastic work, so strain rate and temperature effects arrive coupled and have to be separated afterwards.

From our papers
Huang & Li, the four-paper series in International Journal of Impact Engineering: 173 (2023) 104453, 206 (2025) 105403, Part 1, 206 (2025) 105432, Part 2, and 219 (2027) 105863.
SHPB (Split Hopkinson Pressure Bar)
also called the Kolsky bar

The standard laboratory technique for measuring material response at strain rates of roughly 102 to 104 s-1. A short specimen is sandwiched between two long elastic bars; a striker generates a compressive stress wave in the incident bar, and the strains recorded on the two bars are used to reconstruct the stress, strain, and strain rate history of the specimen.

Its accuracy depends on assumptions that are only approximately met: stress uniformity across the specimen, one-dimensional wave propagation, and negligible friction and inertia. These are the reasons a raw SHPB curve cannot simply be fitted, and why data screening comes before any model identification.

From our papers
Huang & Li, Int. J. Impact Eng. 206 (2025) 105403, Part 1 and 206 (2025) 105432, Part 2, in which the flow stress equation is determined from screened Hopkinson bar data rather than from a fitted curve.
Taylor-Hopkinson impact test
a Taylor cylinder experiment instrumented with a Hopkinson bar

A cylindrical specimen is fired against a Hopkinson bar acting as the anvil, so that the deformed final shape of the specimen and the force history recorded by the bar are both available from the same experiment. It reaches strain rates above what a conventional SHPB can reach, but the deformation field is strongly non-uniform, so the flow stress is recovered by inverse identification: a finite element model of the test is run repeatedly until its predictions match the measurements.

From our papers
Huang & Li, Int. J. Impact Eng. 219 (2027) 105863, where Taylor-Hopkinson tests supply the very-high-rate measurements used to validate the extrapolation.
Strain rate regimes
quasi-static, intermediate, high, and very high strain rate

Strain rate is how fast deformation accumulates, in units of s-1. Different experimental techniques cover different bands: conventional testing machines reach about 10-4 to 100 s-1, servo-hydraulic machines the intermediate range, the Hopkinson bar roughly 102 to 104 s-1, and plate impact or laser-driven experiments beyond that.

Metals typically show a marked upturn in strength somewhere above about 103 s-1. Because no single technique spans the full spectrum, the behaviour in the very-high-rate band has to be extrapolated from lower-rate data, and doing that extrapolation defensibly is a central problem of the constitutive modeling stream.

From our papers
Huang & Li, Int. J. Impact Eng. 206 (2025) 105403, Part 1 for the rate dependence itself, and 219 (2027) 105863 for the extrapolation beyond the tested band.
DIF (Dynamic Increase Factor)
the strain rate sensitivity of strength, expressed as a ratio

The ratio of the dynamic strength at a given strain rate to the quasi-static strength of the same material. It is a compact way of reporting strain rate sensitivity, and it is the quantity most often plotted when comparing materials, but it is defined against a reference condition that is not always stated consistently in the literature, which makes cross-study comparison less straightforward than it appears.

Johnson-Cook model
JC; the multiplicative strain, strain rate, and temperature flow stress equation

The most widely used empirical flow stress equation in impact engineering, built into essentially every explicit finite element solver. It assumes the effects of strain hardening, strain rate, and thermal softening multiply together, each with its own small set of fitted constants.

Its convenience is also its limitation: the multiplicative separability is an assumption, not a measurement, and it is exactly this assumption that a data-driven treatment can test rather than presume.

From our papers
Tested rather than assumed in Huang & Li, Int. J. Impact Eng. 173 (2023) 104453.
TEDI-FS
Tensor-Decomposition-based Extrapolation and Dynamic Identification of Flow Stress

The framework developed in the constitutive modeling stream. Discrete experimental flow stress data over strain, strain rate, and temperature is arranged as an order-3 tensor and decomposed into non-negative one-dimensional latent modes, so that the strain rate dependence is carried by its own mode rather than being assumed to take a prescribed algebraic form. Those modes are then extrapolated beyond the calibrated domain under consistency constraints, and the extrapolated model is identified inversely against Taylor-Hopkinson impact experiments.

From our papers
Introduced in Huang & Li, Int. J. Impact Eng. 219 (2027) 105863, Paper 4 of the series.
MAPE (Mean Absolute Percentage Error)
the standard accuracy metric in this flow stress work

The average of the absolute differences between predicted and measured values, each expressed as a percentage of the measured value. It is used in preference to an absolute error metric because flow stress varies over a wide range across the test matrix, and a percentage measure weights the low-stress and high-stress regions comparably.

DATA-DRIVEN METHODS

CP decomposition
CANDECOMP/PARAFAC; also called canonical polyadic decomposition

A way of writing a multi-way array (a tensor) as a sum of a small number of simple components, each of which is the outer product of one vector per dimension. For flow stress data indexed by strain, strain rate, and temperature, this means expressing the whole three-dimensional dataset in terms of a few one-dimensional curves, one per variable per component.

Unlike matrix decompositions, CP decomposition is essentially unique under mild conditions without needing an artificial constraint such as orthogonality, which is what makes its components candidates for physical interpretation rather than mere basis functions.

From our papers
Huang & Li, Int. J. Impact Eng. 206 (2025) 105432, Part 2, where flow stress over strain, strain rate and temperature is treated as an order-3 tensor, and 219 (2027) 105863.
NN-CP (Non-negative CP decomposition)
CP decomposition with all factor matrices constrained to be non-negative

CP decomposition in which every factor is required to be non-negative, making it the three-way generalization of non-negative matrix factorization. Because the components can then only add and never cancel, each one is a physically admissible partial contribution rather than an arbitrary signed basis vector, which matters when the quantity being decomposed, such as stress, is itself non-negative.

It is the tensor-decomposition step underneath the TEDI-FS framework.

From our papers
The decomposition step in Huang & Li, Int. J. Impact Eng. 206 (2025) 105432, Part 2 and 219 (2027) 105863.
SVD (Singular Value Decomposition)
the matrix counterpart of tensor decomposition

The factorization of a matrix into orthogonal directions ordered by how much of the data each explains. Applied to two-way flow stress data, for instance stress over strain and strain rate at a fixed temperature, it separates the dominant trend from noise and gives a principled way of deciding how many components the data actually supports.

Its factors are orthogonal by construction and can be negative, so they are useful for compression and denoising but harder to read physically than NN-CP factors.

From our papers
Huang & Li, Int. J. Impact Eng. 173 (2023) 104453, where the singular values of two-way flow stress data are used to test whether a decoupled equation is legitimate for a given material.
ANN (Artificial Neural Network)
used here as a surrogate model, not as a black-box replacement for mechanics

A function approximator trained on data. In this work it is used in a deliberately narrow way: to reconstruct and interpolate a response surface that has already been given structure by decomposition, and to build fast surrogates for quantities such as perforation velocity that would otherwise each require a separate finite element run.

A network trained directly on raw stress-strain points would fit the data and explain nothing; combining it with decomposition keeps the strain rate and temperature dependence visible as separate one-dimensional factors.

From our papers
Reconstruction step in Huang & Li, Int. J. Impact Eng. 206 (2025) 105432, Part 2. Used as a surrogate model in Deng, Yang & Huang, Thin-Walled Structures 203 (2024) 112161 and Deng et al., Thin-Walled Structures 216 Part B (2025) 113685.
CNN-BiLSTM
convolutional layers followed by a bidirectional long short-term memory network

A deep learning architecture for signals that have both local shape and temporal order. The convolutional layers pick out local waveform features, and the bidirectional recurrent layers relate them along the time axis in both directions. It is used in the battery stream to classify damage states from acoustic emission recordings, where the diagnostic information lies both in the shape of individual bursts and in how the bursts follow one another.

From our papers

BATTERY SAFETY AND HEALTH MANAGEMENT

PHM (Prognostics and Health Management)
condition monitoring, diagnosis, and remaining-life prediction

The discipline concerned with inferring the internal condition of a system from measurable signals, diagnosing what has gone wrong, and estimating how much useful life remains. For lithium-ion batteries the practical question is whether a cell that has survived a mechanical impact without any visible external damage is still safe to use, which cannot be answered from voltage and capacity alone.

From our papers
Huang & Kim, Journal of Mechanical Science and Technology 40(4) (2026) 2405-2415, an invited review of the field.
Acoustic emission (AE)
passive listening to the elastic waves released by internal damage

When a material cracks, delaminates, or slips internally, it releases a transient elastic wave that a surface-mounted piezoelectric sensor can pick up. The measurement is passive, so it detects damage as it happens rather than after the fact, and it reaches internal events that external inspection cannot see.

In a battery cell the useful signals, such as separator tearing or electrode fracture, arrive mixed with electrochemical and environmental noise, so the sensing has to be paired with a classifier trained to tell them apart.

From our papers
SOC (State of Charge)
the remaining charge as a fraction of full capacity

How full the cell is, from 0 to 100 percent. It matters mechanically as well as electrically: the amount of lithium held in the electrodes changes their stiffness and the stress already stored in the cell, so a cell at high SOC responds differently to the same impact than the same cell at low SOC, and it carries more energy to release if it does fail.

From our papers
Internal short circuit (ISC)
a direct electrical path formed inside the cell between the two electrodes

Mechanical loading can tear or crush the thin polymer separator that keeps the electrodes apart, letting them touch. The resulting short circuit discharges the cell through a very small contact area, generating heat locally and quickly. It is the usual mechanical trigger for thermal runaway, and it can develop with a delay of minutes to hours after an impact that left no visible external damage.

From our papers
Thermal runaway
self-sustaining exothermic failure of a cell

Above a threshold temperature, the chemical reactions inside a lithium-ion cell begin to generate heat faster than the cell can shed it. The temperature then rises on its own, accelerating the reactions further, and ends in venting, fire, or rupture. Because the process is self-sustaining once started, safety engineering is concerned mainly with the conditions that initiate it, which is why the mechanical threshold for impact damage is the quantity worth measuring.

From our papers
Multiphysics FEM
finite element modeling with coupled mechanical, electrical, and thermal fields

A finite element model in which more than one physical field is solved together and each affects the others. For a battery under impact this means the mechanical deformation determines where the separator fails, that failure sets up an electrical short, the short generates heat, and the heat in turn changes the material properties feeding back into the mechanics. Treating any one of these fields alone gives the deformation without the consequence, or the consequence without its cause.

From our papers