Tap a thin panel, such as an aircraft skin, and it rings at a set of well-defined frequencies. Each frequency belongs to a mode shape, a pattern of deflection the panel likes to vibrate in, and for small vibrations these modes are wonderfully well behaved. They are decoupled from one another, and any motion can be written as a sum of them.
Tap the panel harder and this stops being true. Once the deflection grows to something like the panel’s own thickness, bending it also requires stretching, and stretching a membrane is a far stiffer proposition than bending it. The panel’s stiffness is now dependent on the amount of deformation, which means that frequencies of vibration change with amplitude, modes begin to trade energy, and the tidy catalogue of independent mode shapes falls apart. One of the challenges is that engineers increasingly want to operate in this regime but predicting what happens is expensive. A finite element model of a real component may carry a million degrees of freedom, and charting its nonlinear dynamic behaviour can take weeks of computer time.
The standard escape route is a reduced-order model, or ROM. The insight behind this approach is that although the model may have a million degrees of freedom, the motion you actually care about does not weigh all degrees of freedom equally. Instead the motion is confined to a thin, curved surface threading through the million-dimensional space of possible deformed shapes. If you can find that surface, and describe it with two or three coordinates, then you can discard the rest and a simulation that took weeks now takes minutes.
There is a catch though. Commercial finite element packages do not hand you the nonlinear equations of motion but are a black box: you supply loads and the solver returns displacements. So-called indirect reduced-order methods have to find the surface using only what a black box can give you. One of the oldest of reduced-order methods is based on a simple premise: push the structure with static loads shaped like the mode shapes you intend to keep, record how the structure deforms and where it ends up, and then fit a polynomial through the force–displacement response. The set of shapes reachable this way is what we now call the stress manifold, and the ROM amounts to declaring that the only shapes the structure can ever adopt are the ones lying on that surface.
Applied force methods have a reputation for working well and for being held together with rules of thumb. Our two companion papers (1, 2) in Nonlinear Dynamics, written with PhD student Max de Bono and Tom Hill, try to replace the rules of thumb with something an engineer can actually check. They split the problem of trusting a ROM into two questions that the literature has long run together.
Have I drawn the surface accurately?
That is the subject of the first paper, which is all about verification. Many previously assumed that the reduced equations are cubic polynomials, on the reasonable-sounding grounds that the full equations are cubic. However, when you condense a cubic system down onto a lower-dimensional surface, what comes out is not necessarily cubic. The true functions are smooth, so a polynomial will approximate them locally, but “cubic” is a local approximation in exactly the same sense that “linear” is.
Once you appreciate this, verification becomes a well-posed question: raise the polynomial degree until the description of the surface stops changing. We set out an automated algorithm for doing this that adds load cases where they are needed and otherwise touches the finite element model as little as possible, since every call to it is the expensive part.
As a demonstration, a curved double arch with over a million degrees of freedom was reduced to a two-mode ROM and verified in sixty minutes on a desktop computer, reproducing an internal resonance present in the full-order system.
Is it the right surface?
This is validation, the subject of the second paper, and it is much harder. A ROM can be a beautifully accurate map of the wrong territory. The classic failure is internal resonance: as amplitude grows and frequencies drift, two modes can wander into a simple whole-number ratio, lock together and start trading energy back and forth. If the second mode is not among the ones you kept, your ROM has no ability to represent it and will report a periodic vibration that does not exist in the real structure.
The usual remedy is to check against full-order simulations, which defeats the point of building a ROM in the first place, or to build a larger ROM and see whether the answer moves, which is expensive and requires guessing which modes to add.
Our approach starts from a simple observation: to find out whether an orbit would survive in a larger ROM, you do not need the larger ROM. You only need its slope at the points you have already visited. Take a vibration predicted by the small ROM, and at every instant during the cycle attach the flat tangent plane of the bigger surface. Stitched together over one period, these planes form a time-varying validation manifold, and the question becomes whether the orbit, now permitted to drift, actually does so. That question turns out to be a set of linear differential equations with periodic coefficients, which is cheap to solve. If the sideways displacement stays at zero, the orbit exists in the larger ROM and the smaller set of modes was sufficient. If it does not, the solution names the modes you need to add. The stability of the same equations flags bifurcations the small ROM has missed.
The gradients this requires can be extracted while the original surface is being built, so the extra cost is modest. In the double arch again problem, a two-mode ROM of a 100,000 degree-of-freedom model, exhibiting internal resonance, was built and fully validated in under fifteen minutes. Refining the mesh to a million degrees of freedom pushed this to 170 minutes. The equivalent check by full-order simulation would take days. Crucially, nothing in the procedure assumes that an internal resonance exists, or which modes might be involved.
This “validation backbone” is a diagnostic tool and not a prediction. The further an orbit drifts off the original surface, the less the linearisation that was used to construct it can be trusted. Therefore, it tells you reliably that a result is wrong but it does not replace building the larger ROM to find out what is right.
Stepping back, the interesting question here is not really which reduction method is most elegant. Methods built on invariant manifolds are more rigorous and, dimension for dimension, more efficient. But dimensional efficiency is not the same as computational efficiency: a point on a stress manifold is cheap to find and comes with no awkward questions about existence or uniqueness, so a slightly larger reduced model can still be quicker to build and interrogate than a smaller, more sophisticated one. What has really restricted applied force methods in the past is not their mathematics but the absence of any principled way of knowing when you can trust what it produces.
Both papers are open access, and every result in them can be reproduced with AFR Tool written by Max de Bono.
















