pith. sign in

arxiv: 2509.01805 · v1 · pith:OZKFQEIWnew · submitted 2025-09-01 · 🧮 math-ph · cond-mat.soft· math.MP· nlin.PS

Wave-number lock-in in buckled elastic structures: an analogue to parametric instabilities

Pith reviewed 2026-05-22 13:39 UTC · model grok-4.3

classification 🧮 math-ph cond-mat.softmath.MPnlin.PS
keywords bucklingelastic beammodulated Winkler foundationwave-number lock-inparametric instabilitystatic analoguestructural instability
0
0 comments X

The pith

Buckling patterns in elastic beams on modulated foundations lock into periodic modes at specific wave numbers, just as driven dynamic systems do.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper establishes that a static elastic beam resting on a spatially modulated Winkler foundation exhibits wave-number lock-in, switching between quasi-periodic and periodic buckling modes as the modulation strength or wavelength changes. This mirrors the frequency lock-in seen in parametrically driven oscillators, where a periodic drive forces the response to synchronize at particular amplitudes and frequencies. A sympathetic reader would care because the finding shows that purely static structural systems can reproduce instability behaviors normally associated with time-dependent driving, opening routes to control buckling patterns through foundation design alone.

Core claim

The buckling patterns of an elastic beam resting on a modulated Winkler foundation display the same kind of frequency lock-in observed in dynamic systems; through simulations and experiments, compressed elastic strips with modulated height alternate between predictable quasi-periodic and periodic buckling modes.

What carries the argument

The spatially modulated Winkler foundation, whose periodic variation in stiffness selects and locks the buckling wave number.

Load-bearing premise

The observed alternation between quasi-periodic and periodic buckling modes is caused by the foundation modulation itself rather than by unaccounted experimental artifacts, simulation discretization effects, or boundary conditions.

What would settle it

Performing the same compression experiments and simulations on a uniform (unmodulated) foundation and observing the same alternation between periodic and quasi-periodic modes would falsify the claim that modulation drives the lock-in.

Figures

Figures reproduced from arXiv: 2509.01805 by Adel Djellouli, Arnaud Lazarus, Giada Risso, Helen E. Read, Katia Bertoldi.

Figure 1
Figure 1. Figure 1: Frequency lock-in in dynamics and statics (a) An inverted parametric pendulum. (b) Maximum normalized imaginary component of Floquet exponents sj for the inverted pendulum in the first Brillouin zone. (c) A modulated Winkler foundation. (d) Maximum normalized imaginary component of sj for the modulated Winkler foundation. where K¯ = K1/K0 and P¯ = 4π 2P/√ K0EI. Moreover, λ¯ = λ/λ0 and ¯x = x/λ0, with λ0 = … view at source ↗
Figure 2
Figure 2. Figure 2: Wavenumber response of the compressed elas￾tic strip (a) Schematic of our system. (b) When compressed along its bottom edge, the strip buckles out of plane. (c) Crit￾ical wavenumbers as a function of λmod at H1/H0 = 0.5. (d) Some example reconstructions; periodicity marked with grey dashed line. (e) Wavenumber lock-in tongues; the horizontal slice that generates 2c is shown in magenta. and 2λmod, respectiv… view at source ↗
Figure 3
Figure 3. Figure 3: b and Fig. 3c, respectively. All samples were com￾pressed to between 5–6.5% strain, which is much larger than the compressive strain needed for buckling onset (see Supplementary Information Sections VI for details). As expected, the sample with λmod = 8.2 mm exhibits a dis￾ordered postbuckling deformation, while the sample with λmod = 12.2 mm shows an ordered, apparently periodic postbuckling pattern. To q… view at source ↗
Figure 4
Figure 4. Figure 4: Schematic of an inverted pendulum with a shaking base. We study the stability and spectral response of [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Schematic of the modulated Winkler foundation. We study the buckling patters [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Floquet exponents for the inverted pendulum Maximum real (a) and imaginary (b) components of Floquet exponents within the first Brillouin zone for our inverted pendulum problem for a few values of non-dimensionalized amplitude, A¯, as a function of non-dimensionalized period, T /¯ (2π). (c) and (d) are max(Re(s) and max(Im(s)) in the modulation parameter space A¯ and T /¯ (2π); we can clearly see the emerg… view at source ↗
Figure 7
Figure 7. Figure 7: Floquet exponents as a function of loading. Re(s) (a, d), Im(s) (b, e) as a function of axial loading P¯ and output curves (c, f) for two values of period λ¯, both at K¯ = 0.4. We see the critical P¯cr occurs when a pair of si collapse to have real part equal to 0; this point is notated with a dashed line. The imaginary part of the Floquet exponents at the critical loading, P¯cr, tells us if our solutions … view at source ↗
Figure 8
Figure 8. Figure 8: Real, imaginary components of sn and the calculated critical applied load associated with each geometry. In (a, b, c) we show the (a) maximum real component, (b) maximum imaginary component, and (c) calculated critical load as a function of λ¯ for a few K¯ . In (d, e, f) we generalize these results in heat maps. Like for the inverted pendulum, we see the emergence of tongues where the maximum real componen… view at source ↗
Figure 9
Figure 9. Figure 9: Validation of unmodulated theory for finite length simulations (a) Deformation of the first linear buckling mode for H0 = 10 mm (b) DFT of that mode; we locate the wavenumber associated with the peak and calculate the wavelength of the mode. (c) Wavelength as a function of H0; we see that our FE simulations match theory. solutions are quasi-periodic. We plot the locations of our peaks (black) and our exten… view at source ↗
Figure 10
Figure 10. Figure 10: Validation of postbuckling simulations (a) Deformation of the top edge z(x) for a simulation λmod = 8.2 mm, H1/H0 = 0.5 with increasing strain. (b) We see that the locations of the peaks do not significantly move in the post￾buckling regime. (c) Our post-buckling simulations (L = 200 mm) have the same fundamental peak locations as linear buckling simulations (L = 1200 mm). In Fig. 10a, we plot several sna… view at source ↗
Figure 11
Figure 11. Figure 11: Validation of unmodulated theory for Bloch wave simulations (a) For a given geometry (here, H0 = 15 mm), we sweep over many wavenumbers ˜ν at different strains ϵ. (b) Each wavenumber ˜ν yields a 0-frequency response at a different strain, ϵcr; we fit a spline to determine the minimum of the underlying function. (c) Using the calculated ˜νcr, we calculate the critical wavelength and compare to the theoreti… view at source ↗
Figure 12
Figure 12. Figure 12: Modulated Bloch wave results (a) We search for both the critical strain ϵcr and the wavenumber ˜νcr for each geometry. (b) The deformation of the complex unit cell at that critical strain; we can use its deformation to reconstruct the infinite solution. (c) Comparison of finite linear buckling simulations with L = 1200 mm to Bloch wave simulations for H1/H0 = 0.25 [PITH_FULL_IMAGE:figures/full_fig_p023_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Snapshots of the main steps required to fabricate our structures. [PITH_FULL_IMAGE:figures/full_fig_p025_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Laser cut mold with engraved waving pattern on its edge. [PITH_FULL_IMAGE:figures/full_fig_p025_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Example of compression test of a structure with [PITH_FULL_IMAGE:figures/full_fig_p026_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: All experimental samples (green) compared with post-buckling simulations (blue) for [PITH_FULL_IMAGE:figures/full_fig_p027_16.png] view at source ↗
read the original abstract

Parametric instabilities are a known feature of periodically driven dynamic systems; at particular frequencies and amplitudes of the driving modulation, the system's quasi-periodic response undergoes a frequency lock-in, leading to a periodically unstable response. Here, we demonstrate an analogous phenomenon in a purely static context. We show that the buckling patterns of an elastic beam resting on a modulated Winkler foundation display the same kind of frequency lock-in observed in dynamic systems. Through simulations and experiments, we reveal that compressed elastic strips with modulated height alternate between predictable quasi-periodic and periodic buckling modes. Our findings uncover previously unexplored analogies between structural and dynamic instabilities, highlighting how even simple elastic structures can give rise to rich and intriguing behaviors.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript claims that buckling patterns of an elastic beam on a modulated Winkler foundation exhibit a static analogue of parametric frequency lock-in, with the system alternating between quasi-periodic and periodic buckling modes at specific modulation parameters, as demonstrated through finite-element simulations and physical experiments.

Significance. If the central claim is substantiated with controls that isolate the modulation effect, the work would establish a novel connection between dynamic parametric instabilities and static structural buckling, offering a framework for predicting and tuning wave-number selection in modulated elastic systems. This could inform design of adaptive structures or metamaterials, though the current presentation provides only qualitative observations without quantitative benchmarks.

major comments (2)
  1. [Simulation section] Simulation section: no mesh-convergence study, element-type variation, or boundary-condition sensitivity tests are reported. This directly bears on the central claim because the observed alternation between quasi-periodic and periodic modes could arise from discretization artifacts or clamping effects rather than modulation-induced lock-in, as noted in the stress-test concern.
  2. [Experimental results] Experimental results: the manuscript provides no quantitative data (measured wave-numbers, transition thresholds, error bars, or statistical measures) and no control experiments with uniform foundations. Without these, it is not possible to confirm that the mode alternation is caused by the foundation modulation itself rather than fabrication imperfections or generic selection mechanisms.
minor comments (2)
  1. [Abstract] The abstract states that the modes are 'predictable' but does not specify the prediction criterion or how it is derived from the modulation parameters.
  2. [Introduction] Notation for the modulation amplitude and wave-number is introduced without an explicit equation reference in the early sections, making it harder to follow the lock-in condition.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive feedback on our manuscript. We have carefully reviewed the comments regarding the simulation and experimental sections and provide detailed responses below. We believe these points can be addressed through targeted revisions that strengthen the evidence for the wave-number lock-in phenomenon without altering the core claims.

read point-by-point responses
  1. Referee: [Simulation section] Simulation section: no mesh-convergence study, element-type variation, or boundary-condition sensitivity tests are reported. This directly bears on the central claim because the observed alternation between quasi-periodic and periodic modes could arise from discretization artifacts or clamping effects rather than modulation-induced lock-in, as noted in the stress-test concern.

    Authors: We agree that explicit verification of numerical robustness is necessary to support the central claim. In the revised manuscript we will add a dedicated subsection reporting mesh-convergence studies (with at least three successively refined discretizations), results for both linear and quadratic elements, and boundary-condition sensitivity tests (including variations in clamping stiffness and domain length). These additional simulations confirm that the alternation between quasi-periodic and periodic buckling modes remains unchanged, indicating that the lock-in is not an artifact of discretization or boundary conditions. revision: yes

  2. Referee: [Experimental results] Experimental results: the manuscript provides no quantitative data (measured wave-numbers, transition thresholds, error bars, or statistical measures) and no control experiments with uniform foundations. Without these, it is not possible to confirm that the mode alternation is caused by the foundation modulation itself rather than fabrication imperfections or generic selection mechanisms.

    Authors: We acknowledge that the original presentation emphasized qualitative visualization of the mode alternation. To address this directly, the revised manuscript will incorporate quantitative measurements: wave-numbers extracted from digital image correlation of the buckled configurations, transition thresholds as functions of modulation amplitude, and error bars derived from at least five independent experimental realizations. We will also add control experiments performed on uniform (non-modulated) foundations under identical compression protocols, demonstrating that the periodic lock-in does not occur in the absence of modulation and thereby isolating the effect of the height variation. revision: yes

Circularity Check

0 steps flagged

No circularity: claims rest on independent simulations and experiments

full rationale

The paper demonstrates wave-number lock-in via new finite-element simulations and physical experiments on beams with modulated Winkler foundations. No equations, fitted parameters, or self-citations are shown that reduce the reported alternation between quasi-periodic and periodic modes to a prior definition or input by construction. The central analogy is presented as an observed phenomenon rather than a derived result that presupposes itself. The derivation chain is self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The claim rests on the Winkler foundation being a faithful model for the elastic support and on the modulation being the sole driver of the observed lock-in; no free parameters or new entities are mentioned in the abstract.

axioms (1)
  • domain assumption The Winkler foundation model accurately captures the restoring force from the modulated base.
    The abstract frames all results around an elastic beam on a modulated Winkler foundation.

pith-pipeline@v0.9.0 · 5663 in / 1169 out tokens · 35443 ms · 2026-05-22T13:39:07.858359+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

62 extracted references · 62 canonical work pages · 1 internal anchor

  1. [1]

    J. R. Sanmart´ ın Losada, American Journal of Physics52, 937 (1984)

  2. [2]

    A. A. Grandi, S. Proti` ere, and A. Lazarus, Extreme Me- chanics Letters , 101195 (2021)

  3. [3]

    Faraday, Philosophical transactions of the Royal Soci- ety of London , 299 (1831)

    M. Faraday, Philosophical transactions of the Royal Soci- ety of London , 299 (1831)

  4. [4]

    T. B. Benjamin and F. J. Ursell, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences 225, 505 (1954)

  5. [5]

    T. L. Heugel, M. Oscity, A. Eichler, O. Zilberberg, and R. Chitra, Physical review letters 123, 124301 (2019)

  6. [6]

    Zhang, P

    J. Zhang, P. W. Hess, A. Kyprianidis, P. Becker, A. Lee, J. Smith, G. Pagano, I.-D. Potirniche, A. C. Potter, A. Vishwanath, et al., Nature 543, 217 (2017)

  7. [7]

    Smith and J

    H. Smith and J. A. Blackburn, American Journal of Physics 60, 909 (1992)

  8. [8]

    Acheson, Nature 366, 215 (1993)

    D. Acheson, Nature 366, 215 (1993)

  9. [9]

    A. A. Grandi, S. Proti` ere, and A. Lazarus, Nonlinear Dy- namics 111, 12339 (2023)

  10. [10]

    Audoly and Y

    B. Audoly and Y. Pomeau, Elasticity and geometry: From hair curls to the nonlinear response of shells (Oxford press, 2010)

  11. [11]

    Z. P. Bazant and L. Cedolin, Stability of structures: elas- tic, inelastic, fracture and damage theories (World Scien- tific, 2010)

  12. [12]

    Krieger, Nature News 488, 146 (2012)

    K. Krieger, Nature News 488, 146 (2012)

  13. [13]

    P. M. Reis, Journal of Applied Mechanics 82, 111001 (2015)

  14. [14]

    J. Shim, C. Perdigou, E. R. Chen, K. Bertoldi, and P. M. Reis, Proceedings of the National Academy of Sciences 109, 5978 (2012)

  15. [15]

    Lazarus and P

    A. Lazarus and P. M. Reis, Advanced Engineering Mate- rials 17, 815 (2015)

  16. [16]

    Risso, M

    G. Risso, M. Sakovsky, and P. Ermanni, Composite Struc- tures 268, 113946 (2021)

  17. [17]

    Y. Yang, H. Read, M. Sbai, A. Zareei, A. E. Forte, D. Melancon, and K. Bertoldi, Advanced Materials 36, 2406611 (2024)

  18. [18]

    E. Lee, M. Zhang, Y. Cho, Y. Cui, J. Van der Spiegel, N. Engheta, and S. Yang, Advanced Materials 26, 4127 (2014)

  19. [19]

    Doar´ e and S

    O. Doar´ e and S. Michelin, Journal of Fluids and Struc- tures 27, 1357 (2011)

  20. [20]

    Michelin and O

    S. Michelin and O. Doar´ e, Journal of Fluid Mechanics 714, 489 (2013)

  21. [21]

    S. H. Kang, S. Shan, W. L. Noorduin, M. Khan, J. Aizen- berg, and K. Bertoldi, Advanced materials (2013)

  22. [22]

    B. Deng, S. Yu, A. E. Forte, V. Tournat, and K. Bertoldi, Proceedings of the National Academy of Sciences 117, 31002 (2020). 6

  23. [23]

    L. Jin, R. Khajehtourian, J. Mueller, A. Rafsanjani, V. Tournat, K. Bertoldi, and D. M. Kochmann, Pro- ceedings of the National Academy of Sciences 117, 2319 (2020)

  24. [24]

    Nadkarni, A

    N. Nadkarni, A. F. Arrieta, C. Chong, D. M. Kochmann, and C. Daraio, Physical review letters116, 244501 (2016)

  25. [25]

    F. Brau, H. Vandeparre, A. Sabbah, C. Poulard, A. Boudaoud, and P. Damman, Nature Physics 7, 56 (2011)

  26. [26]

    F. Brau, P. Damman, H. Diamant, and T. A. Witten, Soft Matter 9, 8177 (2013)

  27. [27]

    Lagrange and D

    R. Lagrange and D. Averbuch, International Journal of Mechanical Sciences 63, 48 (2012)

  28. [28]

    J. A. Richards, Analysis of periodically time-varying sys- tems (Springer Science & Business Media, 2012)

  29. [29]

    G. W. Hill, Acta Mathematica 8, 3 (1886)

  30. [30]

    Deconinck, F

    B. Deconinck, F. Kiyak, J. D. Carter, and J. N. Kutz, Mathematics and Computers in Simulation 74, 370 (2007)

  31. [31]

    Bentvelsen and A

    B. Bentvelsen and A. Lazarus, Nonlinear Dynamics 91, 1349 (2018)

  32. [32]

    Mora and A

    T. Mora and A. Boudaoud, The European Physical Jour- nal E 20, 119 (2006)

  33. [33]

    Si´ efert, N

    E. Si´ efert, N. Cattaud, E. Reyssat, B. Roman, and J. Bico, Physical review letters 127, 168002 (2021)

  34. [34]

    Brillouin, Wave propagation in periodic structures, 1st ed

    L. Brillouin, Wave propagation in periodic structures, 1st ed. (McGraw-Hill New York, 1946)

  35. [35]

    Chen and J

    X. Chen and J. W. Hutchinson, Journal of Applied Me- chanics 71, 597 (2004)

  36. [36]

    H. Lee, J. G. Bae, W. B. Lee, and H. Yoon, Soft matter 13, 8357 (2017)

  37. [37]

    D. Wu, Y. Huang, Q. Zhang, P. Wang, Y. Pei, Z. Zhao, and D. Fang, Journal of the Mechanics and Physics of Solids 162, 104838 (2022)

  38. [38]

    L¨ owe, X

    C. L¨ owe, X. Zhang, and G. Kovacs, Advanced engineering materials 7, 361 (2005)

  39. [39]

    Florijn, C

    B. Florijn, C. Coulais, and M. van Hecke, Programmable mechanical metamaterials, Vol. 113 (APS, 2014)

  40. [40]

    Ghatak, M

    A. Ghatak, M. Brandenbourger, J. Van Wezel, and C. Coulais, Proceedings of the National Academy of Sci- ences 117, 29561 (2020)

  41. [41]

    Optimal dynamical stabilization

    A. Lazarus and E. Tr´ elat, arXiv preprint arXiv:2508.00006 (2025)

  42. [42]

    Lazarus and E

    A. Lazarus and E. Tr´ elat, arXiv preprint arXiv:2507.13125 (2025)

  43. [43]

    Lazarus, C

    A. Lazarus, C. Maurini, and S. Neukirch, in Extremely Deformable Structures (Springer, 2015) Chap. 1, pp. 1– 53

  44. [44]

    R. A. Calico and W. E. Wiesel, Journal of Guidance, Con- trol, and Dynamics 7, 671 (1984)

  45. [45]

    J. Zhou, T. Hagiwara, and M. Araki, Systems & control letters 53, 141 (2004)

  46. [46]

    Liu and K

    J. Liu and K. Bertoldi, Proc. R. Soc. A 471, 20150493 (2015)

  47. [47]

    ˚Aberg and P

    M. ˚Aberg and P. Gudmundson, The Journal of the Acous- tical Society of America 102, 2007 (1997). 7 Appendix A: Parametric Pendulum We begin by considering the inverted pendulum shown in Fig. 1 of the main text, or Fig. 4 here. Fig. 4. Schematic of an inverted pendulum with a shaking base. We study the stability and spectral response of θ(t) as a function...

  48. [48]

    frequency lock-in

    Floquet theory for the Inverted Pendulum Here, we apply Floquet theory to the inverted pendulum considered in Section A. It follows from Eq. (A13) that the J matrix for this system is given by J(τ) =   0 1 −1 − ¯A cos ¯Ωτ − ¯C   . (C17) Following the procedure described in the above section, we rewrite J(τ) as a Fourier expansion to obtain   0 1 −1 ...

  49. [49]

    Floquet theory for the Periodic Winkler Foundation Here, we apply Floquet theory to the periodic Winkler foundation considered in Section. A. It follows from Eq. (B14) that the J matrix for this system is given by 15 J(¯x) =   0 1 0 0 0 0 1 0 0 0 0 1 −16π4 − 16π4 ¯K cos ¯Ω¯x 0 − ¯P 0   . (C20) Following the procedure described in the above s...

  50. [50]

    Since the period of modulation, λmod will not always divide evenly into the total length, some simulations have a non-integer number of periods; this did not affect our results

    Finite size simulations We simulate finite size samples over a variety of lengths, ranging from L = 200 mm (matching the experimental samples) to L = 1200 mm. Since the period of modulation, λmod will not always divide evenly into the total length, some simulations have a non-integer number of periods; this did not affect our results. To impose uniform co...

  51. [51]

    To study the buckling of our structure, we first compress our unit cell to some strain ϵ, then investigate the frequency response of small-amplitude waves of linear wavenumber ˜ν

    Bloch wave simulations We also investigate the behavior of infinite periodic strips using simulations of a single unit cell [46]; in the modulated case, this unit cell has length λmod, while in the unmodulated case, its length is arbitrary; we denote the length of our unit cell as L∗. To study the buckling of our structure, we first compress our unit cell...

  52. [52]

    Mold Preparation for the modulated strip on a thin block: We laser-cut four acrylic parts with 6.35 mm thickness to form the mold for the modulated strip on a thin block, as illustrated in Fig. 13a. The black part (part #2) has the desired wavy pattern engraved (highlighted in Fig. 14). The depth of the mold can be tailored by adjusting the laser power du...

  53. [53]

    Casting the modulated strip on a thin block: Parts #1 and #2 are stacked and aligned using four pins (Fig. 13b). The assembly is then filled with uncured PVS elastomer (Fig. 13c)

  54. [54]

    13d) are placed on top of the assembly, and pressure clamps are used to ensure uniform thickness (Fig

    Closure and Curing: Parts #3 and #4 (Fig. 13d) are placed on top of the assembly, and pressure clamps are used to ensure uniform thickness (Fig. 13e). The elastomer is cured at room temperature (25 ◦C) for 30 minutes before demolding. At this point, the modulated strip on a thin block is ready

  55. [55]

    Additionally, two top and bottom acrylic plates are laser-cut to enclose the box

    Mold Preparation for thickening the block: A box with two open faces is 3D printed using a BambuLab X1 Carbon printer. Additionally, two top and bottom acrylic plates are laser-cut to enclose the box

  56. [56]

    Final casting and curing: The cured modulated strip on a thin block is placed at the center of the bottom plate (Fig. 13f). The curing box is assembled, and the uncured PVS elastomer mixture is poured inside (Fig. 13g). The assembly is cured at room temperature (25 ◦C) for 30 minutes

  57. [57]

    Finishing T ouches: The edges of the wavy pattern are manually colored black using a Sharpie. 25 Fig. 13. Snapshots of the main steps required to fabricate our structures. Fig. 14. Laser cut mold with engraved waving pattern on its edge. 26 Appendix F: Testing Our samples are compressed using a universal testing machine (Instron 5969, outfitted with a 253...

  58. [58]

    We define these locations to be ( x, z0(x))

    Post-process the undeformed sample to calculate the location of the top edge in the reference configuration. We define these locations to be ( x, z0(x))

  59. [59]

    We define this to be ( xfinal, z0(xfinal))

    Post-process the deformed edge to calculate the location of the top edge in the current configuration. We define this to be ( xfinal, z0(xfinal))

  60. [60]

    Using the known global deformation, ux, project the deformed coordinates ( xfinal, zdeformed(xfinal)) back into the reference frame to calculate ( x, zdeformed(x))

  61. [61]

    Calculate z(x) = zdeformed(x)) − z0(xfinal)) 27

  62. [62]

    Calculate the DFT of z(x) with nfft = 214 All our experimental results can be seen in Fig. 16. We find strong agreement for λmod = 8.2 mm and 12.2 mm, and moderate agreement for all other samples. Fig. 16. All experimental samples (green) compared with post-buckling simulations (blue) for λmod = 8.2, 9.4, 9.8, 10.2, and 12.2 mm. DATA AVAILABILITY STATEMEN...