Materials Science · Architected materials and metamaterials
Reconfigurable structures and instabilities
Hinges, folds, bistability or elastic buckling can create more than one accessible shape. A designed transition can change stiffness, permeability or wave response. Geometry allows functions beyond a fixed dense solid’s response.
Useful background: Composites, fibres and interfaces; Elasticity, plasticity and tensile strength. Each worked case explains the quantities it uses.
Back to the chapterAfter reading, you can
- Distinguish a useful structural transition from uncontrolled failure.
- Derive the force and tangent stiffness of a specified geometrically nonlinear elastic unit from its stored energy.
- Classify equilibria under a dead force, a prescribed displacement and an actuator with finite stiffness.
- Explain a trigger, retained state and return path without treating a static equilibrium curve as a time history.
- Account for work and released energy through a complete cycle, and identify evidence needed to distinguish useful reconfiguration from damage.
A structure can change its configuration
A shallow arch can be pressed through its flat position and remain inverted after the pressing force is removed. A different arch returns immediately. The useful question is not just whether a shape changes, but which state is stable under the actual supports and controller, how the transition is triggered, and what brings it back.
Hinges, folds, bistability or elastic buckling can create more than one accessible shape. A designed transition can change stiffness, permeability or wave response. Geometry allows functions beyond a fixed dense solid’s response.
A monostable unit has one stable configuration under a specified set of conditions. A bistable unit has two; multistability allows more. These are structural states, not automatically different chemical phases. Elastic members can support two shapes because their arrangement changes the energy landscape. Conversely, an apparently retained shape could involve friction, plastic strain or damage rather than a second elastic equilibrium.
Geometry, architecture and effective response explains why arrangements change stiffness, while Auxetic behaviour and mechanical lattices follows a particular rotation mechanism. This lesson instead follows one nonlinear unit through a transition and its control conditions. Reconfiguring a pore or a wave structure could change its function, but that function needs its own transport or wave evidence. Shape change alone supplies no permeability or band-gap value.
Build one explicit elastic energy landscape
Take an invented planar unit with fixed pin supports at (−a,0) and (+a,0) and a pin-jointed apex at (0,y). A frictionless vertical guide keeps the apex’s horizontal coordinate zero; it can provide a horizontal reaction but no vertical force. The signed height y is positive upward. The drawing uses the same length scale in both directions. Clearance permits either sign of y without contact.
Each member is an ideal axial spring with positive stiffness k in N/m and natural length L₀. Its force is k(L−L₀), positive in tension. Each stays straight and transmits no bending moment. The guide and this spring idealisation remove lateral and member-bending modes: their stability is not being proved. Gravity, friction, plasticity, thermal changes and fracture are omitted. No named constituent material is being qualified over this motion.
The member length is L(y)=√(a²+y²). Choose L₀=√(a²+h²), so the two mirror shapes y=+h and y=−h have the same natural member length. Each spring stores ½k(L−L₀)²; the two together store U(y)=k[L(y)−L₀]² joules. All spring constants and stored energies are nonnegative, even though the overall structure will have a negative force slope in part of its path.
For the numerical model use a=10.0 mm, h=2.00 mm and k=1000 N/m per member. L₀=10.198039 mm. Both ±h states have U=0. Flattening to y=0 shortens each spring to a and raises U to 39.2195 μJ. That is the barrier along this unloaded one-coordinate route, not a calibrated thermal activation energy or a work loss per cycle. A load changes the landscape.
Here U means elastic strain energy, not the molar Gibbs energy used to compare chemical phases. We will add the mechanical controller’s energy or load potential before judging equilibrium. Metastability and persistence provides the general distinction between a local minimum and a lower alternative; the calculation here states which mechanical motions are actually allowed.
A positive spring can produce a negative structural slope
Let Q be the externally applied apex force, positive upward and conjugate to y: a small displacement dy has external work Q dy. In quasistatic equilibrium Q=dU/dy. Since dL/dy=y/L, differentiation gives Q(y)=2k(L−L₀)y/L=2ky(1−L₀/L). The same result follows by resolving the two axial member forces vertically.
The tangent stiffness is the local slope K_t=dQ/dy=2k[1−L₀a²/L³], in N/m. It describes a small change of force per change of height about a specified configuration. It is not Q/y, a constituent Young’s modulus, or a material strength limit. Geometry, rather than a negative spring constant, creates its changing sign.
At y=±h, Q=0 and K_t=76.9231 N/m. At y=0, Q=0 but K_t=−39.6078 N/m. Thus zero force is not enough to identify a stable state. The first figure places the energy and force curves on the same signed-height axis so that a well, a crest and a changing force slope can be compared.
For a prescribed dead force Q̄, meaning a force whose magnitude and direction stay fixed during a perturbation, the total potential is Π_F(y)=U(y)−Q̄y. Its stationary points satisfy Q(y)=Q̄. For a small perturbation δy about an equilibrium, the leading energy change is ½K_t(δy)². Positive K_t gives a local minimum in this one-coordinate conservative model; negative K_t gives an unstable stationary point. This local test identifies a restoring equilibrium; without damping, a disturbed system can oscillate rather than settle.
If K_t is zero, the quadratic test alone is inconclusive: examine higher terms and the actual constraints. A local minimum is also not a comparison of every possible state or every mechanical mode. A real unit can bend out of plane or buckle a member before the predicted apex transition; the guide and axial model deliberately do not represent those possibilities.
Follow a trigger and a return path under force control
Start at y=+h with zero apex force. Reducing Q̄ below zero presses downward and moves the upper equilibrium toward the flat position. It loses its local minimum at a limit point, or fold, where K_t=0. Solving L³=L₀a² gives y_c=√[(L₀a²)^(2/3)−a²]. The upper fold is at y_c=1.147144 mm and Q̄=−30.1915 mN. Its mirror is at −y_c with +30.1915 mN.
A force just beyond the downward limit no longer has an upper equilibrium nearby. The apex can move rapidly toward the lower branch: snap-through. A static curve locates equilibria and the limiting loss of stability; it does not give a transition time, velocity, overshoot or transient peak stress. A numerical continuation method may trace the unstable middle branch without showing that a force-controlled specimen can occupy it.
For the force-cycle calculation below, stipulate slow changes along stable branches, a force held during each jump, and enough damping to settle after each jump. The limiting downward jump connects the upper fold to the other stable equilibrium at y≈−2.317007 mm at the same limiting force. This is a destination after settling in the ideal model, not a predicted dynamic trajectory.
Then return Q̄ to zero along the lower branch. The apex comes to y=−h and remains there in this model: removing the force does not force a return to the upper well. To restore the original shape by the mirror path, increase Q̄ upward to the other fold, pass its positive threshold, allow the reverse snap and settling, then reduce Q̄ to zero at +h. An alternative imposed-displacement path could return it differently.
The 39.2195 μJ unloaded barrier is not a constant barrier under those forces. The term −Q̄y tilts the potential; the departing minimum and adjacent crest merge at the fold. Nor does a demonstrated retained shape establish a useful lifetime, resistance to vibration or repeatable operation after many cycles. Those require observations of the actual device and environment.
Include the controller before interpreting a negative slope
Under ideal direct displacement control, the apex height y itself is prescribed. There is no free apex δy perturbation to which the dead-force test applies. The controller supplies the required reaction Q(y), including the negative-slope segment. This does not stabilise every bending or lateral mode in a real specimen, and it does not make the free-standing flat state stable after the constraint is removed.
A finite apparatus is more informative. Connect the apex to a driven stage through a positive linear vertical spring of stiffness K. Write the stage coordinate as Y after absorbing the spring’s constant natural-length offset, so its force is zero when y=Y. The commanded coordinate Y, not y, is prescribed. The apparatus force on the apex is K(Y−y), and its stored energy is ½K(Y−y)². The total potential for a fixed Y is Π_K(y;Y)=U(y)+½K(Y−y)².
Equilibrium requires Q(y)=K(Y−y). Stability against an allowed apex perturbation now uses ∂²Π_K/∂y²=K_t+K. A negative intrinsic K_t can therefore be part of a stable combined system if K_t+K>0. The positive controller contributes restoring curvature. Checking only the specimen’s force slope or only the commanded stage displacement would miss this distinction.
At Y=0 and y=0, K=20 N/m gives total curvature −19.6078 N/m, so that equilibrium is unstable. K=60 N/m instead gives +20.3922 N/m and a local minimum. The specimen’s own flat-configuration slope is unchanged. With K=20 and Y=0, the combined system has minima at y≈±1.396707 mm; these differ from ±h because the attached spring still exerts force.
In this ideal one-coordinate model K_t is most negative at y=0. Thus K>2k(L₀/a−1)=39.6078 N/m makes Π_K strictly convex in y and removes the fold in that model. Equality at this threshold needs higher-order analysis, and a real apparatus also has inertia, damping, bandwidth and additional degrees of freedom. A scalar inequality here is not a universal stability condition for a lattice or a solid with a negative modulus.
| Control at the flat configuration | Relevant test | What follows in this model |
|---|---|---|
| Dead force Q̄=0 | K_t=−39.6078 N/m | Unstable free apex equilibrium |
| Direct prescribed y=0 | Apex δy is not an allowed perturbation | Controller holds that coordinate; other modes untested |
| Prescribed Y=0; K=20 N/m | K_t+K=−19.6078 N/m | Flat combined equilibrium unstable |
| Prescribed Y=0; K=60 N/m | K_t+K=+20.3922 N/m | Flat combined equilibrium a local minimum |
A repeatable mechanism needs limits
Check trigger load, energy barriers, contact, friction and damage over cycles. A single snap-through does not establish reversible long-life operation. External actuation or stored elastic energy supplies the work; reconfiguration does not provide free energy.
For direct, quasistatic displacement control of the ideal conservative springs, Q=dU/dy and the work between endpoints is their energy difference. A complete path +h→−h→+h retraced without dissipation has ∮Q dy=0. Energy must be supplied while rising toward a barrier and can be returned while descending. A negative force or tangent does not provide energy from nowhere.
Force-controlled jumps with settling are different. Use the energy balance W_ext=ΔU+ΔT_kin+D, where T_kin is kinetic energy and D≥0 is mechanical energy removed from the resolved elastic storage and motion through damping or other dissipative channels in this account. For settled endpoints ΔT_kin=0. During the limiting downward jump, Q̄≈−0.0301915 N and Δy≈−0.00346415 m, so W_ext≈104.588 μJ. The spring energy changes by 4.47265−17.54487=−13.0722 μJ. The balance requires D≈117.660 μJ.
The force does positive work while the springs also release energy. The mirror jump has the same limiting loss. For the complete model force cycle returning to the same shape and rest state, net work is about 235.320 μJ. Smooth stable-branch portions are reversible here; the two assumed damped jumps account for the net loss. With no damping, the apex would not simply settle and this quasistatic hysteresis construction would not describe its motion.
The first figure’s force panel shows that idealised force cycle and the direct-displacement alternative. Its horizontal jumps mean a force held during a limiting snap, not zero duration or a resolved path through intermediate equilibria. A measured loop needs the force conjugate to the measured displacement and a defined system boundary. If an actuator spring is included, its stored energy and work at the moving stage must be accounted for; mixing stage force with a different displacement can give a false loss.
A loop that does not return the internal state or velocity to its starting value can include stored or kinetic energy changes. Even for a closed settled cycle, total dissipated work alone does not identify friction, viscoelasticity, damage or a particular geometric mode. Compare unloaded configuration, local strain, contact, trigger force, cycle work and the intended function over the actual observation interval.
Separate designed transitions from failure and rate effects
An elastomer sheet with rotating units and compliant ligaments can retain a changed arrangement after load release, while another pattern returns. Mechanical loading curves and observed shapes can be compared with finite-element predictions to test the proposed snap-through mechanism. The relation between unloading and a deliberate reset must be recorded; a label such as “recoverable” does not itself specify that protocol or a cyclic lifetime.
Real constituents can introduce their own time dependence. In a tested viscoelastic hole-patterned metamaterial, changing compression rate delayed or suppressed an internal pattern change under particular lateral confinement, changing the hysteretic response. That observation is evidence for a competition of loading and internal relaxation times in that system, not a universal rule that faster loading always gives more geometric snapping or the same loss.
Member buckling, joint bending, asymmetric imperfections and contact can introduce pathways absent from a guided one-coordinate model. Compression of an ideal spring is not proof that a real slender member can carry the same force without buckling or yielding. A negative structural force slope is also not evidence of fracture by itself: elastic instability can be designed, while accumulating damage can alter a superficially similar loop.
When several units are coupled, transitions may occur sequentially, localise, interact or follow different paths. Test the constraints and accessible modes of the assembly rather than multiplying one unit’s trigger force or stored energy without justification. If the useful function is a reconfigured acoustic or electromagnetic response, verify it in the retained state with its frequency, direction, losses and boundaries.
A useful reconfiguration claim therefore names the actuated coordinate, stable states, trigger and reset, complete energy account, observed function and repeatability interval. It remains conditional on the actual controller and environment. The model here explains how to make those questions quantitative without claiming a construction recipe, a qualified actuator or indefinite durability.
Worked examples
Work through the case
A bistable unit snaps to a second state under load. Must it return when the load is removed?
- The second state may be stable without that load.
- Returning requires another path or input.
- Inspect the full force–displacement cycle.
No. Bistability can retain the second state.
Zero force does not settle stability
For a=10.0 mm, h=2.00 mm and k=1000 N/m per member, compare y=+h, 0 and −h with no apex force. Find the stored energy and the local force slope.
- Convert lengths to metres. L₀=√(0.0100²+0.00200²)=0.010198039 m. At ±h both members have that natural length, so U=0 and Q=0.
- At y=0, L=a=0.0100 m. The two-spring energy is 1000(0.0100−0.010198039)²≈3.92195×10⁻⁵ J=39.2195 μJ. The signed force is still zero because the two compressed members are horizontal.
- K_t(±h)=2k h²/L₀²≈76.9231 N/m, whereas K_t(0)=2k(1−L₀/a)≈−39.6078 N/m. Under zero dead force, the outer equilibria are local minima and the flat equilibrium is unstable in the allowed apex coordinate.
- The model says nothing about a real member’s yield or bending stability. Its positive k and positive energy are compatible with the negative overall tangent caused by geometry.
Both ±2.00 mm states are stable zero-force wells in this one-coordinate model. The flat zero-force state is a 39.2195 μJ crest with a negative tangent, not a third stable unloaded shape.
The same specimen, two apparatus stiffnesses
At commanded Y=0, compare the flat equilibrium with K=20 and 60 N/m. For K=20, where are the other equilibria? Does moving a structure to y=0 prove it stays there after detachment?
- For y=0, Q=0 and the apparatus force K(Y−y)=0, so the force balance is satisfied for either K. Stability requires the total curvature, not force balance alone.
- Add K to −39.6078 N/m. K=20 gives −19.6078 N/m, unstable; K=60 gives +20.3922 N/m, a local minimum in the combined model.
- At Y=0, nonzero equilibria satisfy 2k(1−L₀/L)+K=0. Therefore L=L₀/[1+K/(2k)]. For K=20, L≈10.097068 mm and y=±√(L²−a²)≈±1.396707 mm. Their total curvature is +38.6519 N/m.
- Removing the apparatus removes its added curvature and any stored-energy contribution. The free flat equilibrium again has negative K_t. A held or stabilised configuration is not automatically a retained free-standing one.
The flat equilibrium is unstable with K=20 and locally stable with K=60 N/m. The K=20 combined system instead has stable mirror states near ±1.396707 mm at Y=0; detachment changes the problem.
Trace the state through a full force protocol
Start at the upper unloaded shape. Lower the controlled force quasistatically until snap-through, allow damping to settle, and then return the force to zero. Find the limiting trigger and retained state. What is required for the mirror return?
- Solve K_t=0: L_c=(L₀a²)^(1/3) and y_c=√(L_c²−a²)≈+1.147144 mm. Substitution in Q(y) gives −30.1915 mN. This is the limiting downward trigger, with motion starting after it is passed or suitably perturbed.
- At that limiting force the other stable equilibrium is near −2.317007 mm. Damping and the stipulated held force permit settling there; a static root does not specify the transient overshoot.
- Raising the force back to zero on the lower branch brings y to −2.00 mm. It remains in the second unloaded well rather than returning to +2.00 mm.
- Continue upward force to the mirror fold at −1.147144 mm and +30.1915 mN, pass the threshold, settle on the upper branch and reduce force to zero. This completes the idealised return protocol; an actual reset may instead use displacement or another actuator.
The limiting downward force is −30.1915 mN. After settling and unloading, the retained shape is y=−2.00 mm. Returning needs an appropriate reverse input, not merely removal of the original force.
Account for the energy of a snap before calling it loss
Use the limiting downward jump: y_i=+1.147144 mm, y_f=−2.317007 mm, Q̄=−30.191474 mN, U_i=17.544872 μJ and U_f=4.472651 μJ. The force is held and both endpoints are at rest. Find external work and D. Explain why D can exceed the initial spring energy.
- Δy≈−3.464151 mm=−0.003464151 m. Held-force work is Q̄Δy≈(−0.030191474)(−0.003464151)=0.000104588 J=104.588 μJ. Both force and displacement are downward, so the work is positive.
- The stored spring-energy change is U_f−U_i≈−13.072221 μJ. With no net endpoint kinetic energy, D=W_ext−ΔU≈104.587823−(−13.072221)=117.660 μJ.
- The force supplies work during the jump while the springs release some stored energy. Thus D is not limited to U_i, and it is not the unloaded 39.2195 μJ barrier.
- For the stated full symmetric force cycle, the mirror jump adds the same loss and the stable-branch work completes the return to the same U. Net external work is 2D≈235.320 μJ. The result depends on damped settling and the specified boundary/control; the conservative direct-displacement round trip would instead have zero loop work.
W_ext≈104.588 μJ and D≈117.660 μJ for this limiting settled jump. The complete symmetric force-cycle work is about 235.320 μJ; external work as well as stored energy feeds the dissipative channels.
Try the reasoning
Is all buckling necessarily destructive?
Show the answer and reasoning
No.
Controlled instability can be a designed mechanism.
Why measure cycles?
Show the answer and reasoning
Joints, friction and local strain can evolve.
One transition does not establish durability.
Keep a=10.0 mm and k=1000 N/m, but choose h=1.00 mm. Recalculate the two unloaded wells and the unloaded barrier. Does halving h merely halve that barrier?
Show the answer and reasoning
The wells are at y=±1.00 mm. L₀≈10.049876 mm and U(0)≈2.48758 μJ, about 6.3427% of the original barrier, not one half.
Use the changed natural length √(a²+h²), then k(a−L₀)². The nonlinear geometry makes the barrier highly sensitive to h. In the shallow limit L₀−a≈h²/(2a), so the barrier scales approximately as h⁴; that approximation explains the near one-sixteenth result without replacing the exact calculation.
Double each member’s k to 2000 N/m while retaining the original a,h and apparatus K=60 N/m. Is the combined flat equilibrium at Y=0 still a strict local minimum by the quadratic test?
Show the answer and reasoning
No. The intrinsic flat tangent doubles to −79.2156 N/m; adding 60 gives −19.2156 N/m. A strictly positive combined curvature there requires K>79.2156 N/m.
Increasing the positive member stiffness also increases the magnitude of the geometric negative slope. A controller that was sufficient for one specimen need not stabilise a stiffer one. This is the specified one-coordinate test, not a general lattice or active-control criterion.
Someone sets K exactly to 39.6078 N/m, the rounded magnitude of the original flat negative tangent, and declares the configuration stable because the sum is zero. What is missing?
Show the answer and reasoning
Zero quadratic curvature alone is inconclusive; use exact parameters and the next nonzero energy terms or a full perturbation analysis.
The strict positive-curvature sufficient test does not include equality. Rounding a threshold also does not produce an exact physical equality. In this conservative model the actual potential can be examined further, but neither “zero means stable” nor “zero always means unstable” follows from the second derivative alone.
A separate slow operation uses 0.40 mJ of mechanical work and ends at rest, but its internal stored energy is 0.15 mJ higher than at the start. There are no other work or heat inputs in the supplied energy account. What D follows, and is 0.40 mJ automatically dissipated work?
Show the answer and reasoning
D=0.40−0.15=0.25 mJ. The complete 0.40 mJ is not automatically dissipation.
Apply W_ext=ΔU+ΔT_kin+D with both endpoints at rest. A path ending in a different stored-energy state is not a closed cycle. Returning the external force to its starting value alone does not prove that every internal state has returned.
In a separate invented device record, complete settled-cycle work falls from 0.30 mJ in cycle 1 to 0.20 mJ in cycle 100; the trigger magnitude falls from 30 to 24 mN. The unloaded geometry still looks the same. Find the relative changes. Do they identify damage or establish a lifetime?
Show the answer and reasoning
Work falls by 33.333…% and trigger magnitude by 20%. Those two observations identify neither a unique cause nor a lifetime.
Use the first-cycle values as denominators: (0.30−0.20)/0.30 and (30−24)/30. Frictional conditioning, viscoelastic history, geometry or damage could contribute. A matching unloaded outline does not establish unchanged local strain, stiffness or useful performance, and two cycle observations do not justify indefinite extrapolation. These supplied observations are not predictions of the ideal arch.
A latched pore structure retains its expanded outline at zero external force. Must it have the same permeability on every cycle, or the acoustic gap predicted for an unrelated periodic lattice?
Show the answer and reasoning
No. Retained geometry is one structural observation, not a measurement of either function.
Actual pore connectivity, contacts, dimensions and flow conditions matter to permeability. Wave response additionally depends on frequency, mode, boundaries, damping and constituent properties. Measure or model the intended retained-state function under its own conditions.
What to carry forward
- Identify stable states and transition paths.
- Track energy, actuation and repeatability.
- Different control variables change the permitted perturbations and the total energy: force balance alone does not establish stability.
- A negative structural tangent can arise from positive elastic members; include apparatus stiffness and every relevant mode before transferring the result.
- A retained second well needs a reset path. A static fold identifies a limit, not the time or stresses of the dynamic snap.
- Count stored, kinetic and dissipated energy with the work at the correct boundary. Geometry recovery and a few closed loops do not establish cyclic durability.
References
Teaching resources and research supporting this lesson’s explanations and examples.
- Bistable Auxetic Mechanical Metamaterials Inspired by Ancient Geometric Motifs · 2016 preprintA. Rafsanjani and D. PasiniPrimary research preprint · selected passages
- Structural Mechanics · Lecture 10: Advanced Topic in Column BucklingTomasz Wierzbicki · MIT OpenCourseWareUniversity notes · selected snap-through and control discussion · PDF
- Aerospace Structures · Chapter 10: Structural Stability of Discrete Conservative Systems · 2022Eric R. Johnson · Virginia TechUniversity textbook chapter · selected energy and limit-point analysis · PDF
- Viscoelastic Snapping Metamaterials · Journal of Applied Mechanics 86, 111012 · 2019David M. J. Dykstra, Joris Busink, Bernard Ennis and Corentin CoulaisPrimary research paper · selected experimental and discussion passages · PDF