## Abstract

Structural instability, once a catastrophic phenomenon to be avoided in engineering applications, is being harnessed to improve functionality of structures and materials, and has been a catalyst of substantial research in the field. One important application is to create functional metamaterials that deform their internal structure to adjust performance, resembling phase transformations in natural materials. In this paper, we propose a novel origami pattern, named the Shrimp pattern, with application to multi-phase architected metamaterials whose phase transition is achieved mechanically by snap-through. The Shrimp pattern consists of units that can be easily tessellated in two dimensions, either periodically with homogeneous local geometry or non-periodically with heterogeneous local geometries. We can use a few design parameters to program the unit cell to become either monostable or tune the energy barrier between the bistable states. By tessellating these unit cells into an architected metamaterial, we can create complex yet navigable energy landscapes, leading to multiple metastable phases of the material. As each phase has different geometries, the metamaterial can switch between different mechanical properties and shapes. The geometric origin of the multi-stable behavior implies that, conceptually, our designs are scale-independent, making them candidates for a variety of innovative applications, including reprogramable materials, reconfigurable acoustic waveguides, and microelectronic mechanical systems and energy storage systems.

## 1 Introduction

Elastic structures with two or more possible equilibria can transition from one state to another via a rapid process known as snap-through. Snap-through allows plants and animals to store elastic energy and release it suddenly to generate rapid and powerful motions. Through a suite of structural modifications of their raptorial appendages, mantis shrimps generate forceful predatory strikes in a fraction of a second, which is so fast that it causes cavitation in the water [1,2]. Venus flytraps snap their leaves together rapidly to capture insects, which is enabled by the doubly curved geometry of the leaves [3]. Ladybird beetles fold their wings in a non-rigid foldable origami pattern to store elastic energy when stowed, which allows quick release of their wings when they start flying [4]. Similarly, the earwig wings have incompatible folding patterns, which remain open by a bistable locking mechanism during flight and self-fold rapidly without muscular actuation [5].

Such multi-stable structures, if adopted as building blocks of architected materials, can lead to some interesting and unprecedented properties beyond conventional materials [6–9], such as energy trapping [10,11], elastic pattern formation [12,13], and guided transition waves [14]. In this work, we aim to design an architected material that can switch properties by multi-stable snapping [15–18], resembling the concept of phase transitions. Owing to the different structures at the atomic and molecular scale, the same material displays different properties in different phases. For natural materials, the transitions between phases are characterized by changes in thermodynamic states. The concept of a phase being a rearranged structural form of materials can be generalized at larger scales, which results in architected materials whose internal structures can be intentionally changed to display different properties on demand.

As suggested by examples found in nature [3–5], origami patterns can organize geometry of thin sheets and membranes to produce functionality by harnessing their inherent elastic instabilities [19–23]. There are a few known origami patterns that lead to structures with multi-stability. Examples include the square twist [24], the Kresling tube [25], and the hypar origami [26]. However, many of these patterns cannot be easily tessellated due to irregular geometry (e.g., Kresling tube), or if they can be tessellated, their deformation modes and multi-stable energy landscape differ significantly between the unit cell and tessellation, making their behavior difficult to predict (e.g., the square twist [24] and hypar origami pattern [26]).

In this paper, we present a multi-stable origami pattern, named the Shrimp pattern (see Fig. 2), whose geometry is tailored to accommodate tessellations easily. The pattern is inspired by the segmented structure of the shrimp, and hence named after it, as shown in Fig. 1. The Shrimp pattern combines features of both rigid and non-rigid origami, allowing the pattern to behave as a rigid origami [27] in a range of stable configurations, while the switching between different ranges of stable configurations has to go through an energy barrier by snap-through (i.e., non-rigid origami) [28]. The bistability of a unit cell of this pattern can be programed by a few geometric parameters. Hence, if used for architected materials, the Shrimp tessellated material would have multiple stable phases, and switch between phases is possible through snap-through, which is interpreted as mechanical phase transition. We further show that the Shrimp tessellated material exhibits different mechanical properties in its different phases. Moreover, because each Shrimp unit cell is bistable, by tessellating the Shrimp unit cells, the number of phases of the resultant architected material grows combinatorially.

The remainder of the article is organized as follows. Section 2 describes the geometry of the standard Shrimp pattern unit that tessellates on a plane. Section 3 explains how the energy landscape of a single Shrimp unit cell can be programed. Section 4 verifies the programable bistable behavior by means of a reduced-order bar-and-hinge numerical model. Section 5 elaborates on the creation of the multi-stable architected material by tessellating the Shrimp unit cells. More importantly, we demonstrate that the Shrimp tessellated material can switch property by mechanical phase transition. We summarize current work and propose future development in Sec. 6. Appendix A explains the strategy to tessellate the Shrimp pattern in curved configuration. Appendix B supplements the paper by offering tips for fabricating the non-developable Shrimp pattern.

## 2 Geometry of the Standard Shrimp Unit

*a*, front panel length

*b*, rear panel length

*c*, and tail length

*d*. The angular parameters include front panel angle

*α*

_{F}, rear panel angle

*α*

_{R}, and rib (triangle) panel angle

*α*

_{B}. For standard Shrimp pattern designs, edge length

*c*satisfies the following condition:

*a*lie in the

*xz*-plane, as shown in Fig. 2. Although within each unit cell the tail panels are not connected to the rib panels, in an assembly/tessellation of the Shrimp pattern, the rib panels are actually between the tail panels and front panels, and they should be attached to the last row of tail panels.

*ψ*≤ 2

*α*

_{B}, the Shrimp pattern allows rigid origami behavior with a single degree-of-freedom (DOF). The configuration is determined by any of the three edge inclination angles (

*ψ*/2) (=

*ψ*

_{2}),

*ϕ*

_{F}, and

*ϕ*

_{R}, which are related as follows:

*H*,

*H*

_{1},

*H*

_{2}), width (

*W*), and lengths (

*L*,

*L*

_{1},

*L*

_{2},

*L*

_{3}) of the Shrimp unit cell are given by

*ψ*≤ 2

*α*

_{B},

*O*

_{6}

*O*

_{7}is a mountain or valley fold. When

*O*

_{6}

*O*

_{7}is a mountain fold, the Shrimp pattern is said to be in its eggbox mode (

*E*), and the “+” sign shall be used in Eq. (21). When

*O*

_{6}

*O*

_{7}is a valley fold, the Shrimp pattern is said to be in its Miura mode (

*M*), and the “−” sign shall be used. This feature is inherited from the Morph pattern investigated in Ref. [29], as the Shrimp pattern shares the same geometry with the Morph pattern by deleting the tail and rib panels. Figure 4 shows how the geometric parameters (a)

*α*

_{F}, (b)

*α*

_{R}, (c)

*d*, and (d)

*b*/

*a*affect the Poisson’s ratio of a Shrimp unit. We note that the in limit of

*d*= 0, the value of

*ν*

_{yx}recovers that of the Morph pattern [29]. The above derivation only applies to the rigid origami range of the Shrimp pattern, when the folding process is periodic across all unit cells. However, in reality, due to the compliance of the panels, a true rigid origami behavior is rarely observed [28].

## 3 Creating Bistability by Breaking Connectedness of Configuration Space

*ϕ*and

*ψ*(=

*ψ*′ =

*ψ*″) are related by the following formula:

*ψ*′ and

*ϕ*.

*ψ*′ and

*ψ*″ cannot exceed 2

*α*

_{B}, and thus a certain range of folding is blocked under rigid origami assumption, as indicated by the dashed lines in Fig. 5. The solid parts of the curves refer to the remaining rigidly foldable ranges of the Shrimp pattern, which is separated into two disconnected parts. The two limit configurations of the blocked range of rigid folding (dashed lines in Fig. 5) in terms of

*ϕ*are given by

*ϕ*

*α*

_{F}.

In reality, the compliance of materials offers extra DOF that allow the Shrimp pattern to transition between the two disconnected rigid origami configurational ranges through non-rigid deformation. The tail panels are added to allow the blocked range to be overcome by gentle panel bending; otherwise, the blocked range of rigid folding can only be overcome by stretching the materials, which is likely to generate rupture or other types of irreversible damages. Together, the addition of rib panels and tail panels leads to bistable snapping between the two disconnected parts of the rigid origami configurations of the Shrimp pattern. The conditions for the emergence of bistability are discussed later.

As illustrated in Fig. 6(a), the folding energy of the system increases as the pattern is folded away from the initial configuration, but the bending energy only appears within the blocked range of rigid folding, i.e., when *ψ* > 2*α*_{B}. While the bending energy gradually vanishes beyond the range of non-rigid folding, the folding energy keeps increasing. Therefore, if the stored energy becomes smaller than the peak energy after the non-rigid deformation, the Shrimp pattern will experience bistability. We denote the initial configuration as configuration 0, the peak energy state as configuration 1, and the state after non-rigid folding as configuration 2, i.e., the metastable state (see Fig. 6(b)).

### 3.1 Estimating Peak Energy During the Non-Rigid Folding Range.

Although the actual non-rigid deformation of the Shrimp pattern during the transitioning is a complex combination of bending, folding, and stretching, we may consider a simplified deformation process to estimate the peak energy by assuming isometric deformations. We assume that only the front and tail panels will bend along one of their diagonals, as demonstrated in Fig. 6(a). Therefore, *ψ* can be different from *ψ*′ and *ψ*″ and can exceed 2*α*_{B}, as depicted in Fig. 6(b)〈0〉.

*δ*

_{F}and the bending angle of the tail panels as

*δ*

_{T}. We define

*ρ*as the rotation angle of edge

*O*

_{1}

*O*

_{2}(see Fig. 6(b)〈1〉), which can be approximated by

*O*

_{5}, using spherical trigonometry, we obtain

*O*

_{8}, we obtain

*ρ*is largest when

*ψ*= 2

*α*

_{R}, and consequently,

*δ*

_{F}and

*δ*

_{T}become maximal. Therefore, the total bending energy reaches peak when

*ψ*= 2

*α*

_{R}, when the two rear panels become coplanar. According to Eqs. (26) and (28), the maximal bending angles are given by

*ψ*= 2

*α*

_{R}is calculated as

*ψ*

_{0}and

*ϕ*

_{0}. The symbol $k0F$ denotes the linear folding stiffness. The deformed angles at configuration 1 can be estimated by taking

*ψ*

_{1}= 2

*α*

_{R}and

*ϕ*

_{R,1}= 0 using the same set of equations. The total energy at peak is a combination of folding and bending energy:

### 3.2 Estimating Stored Energy After Snapping.

*γ*

_{FR},

*γ*

_{RR}, and

*γ*

_{RT}experience deformation, and the other folding hinges stay at the same angles as in the initial configuration. The three dihedral angles can be derived based on Eqs. (35) and (36), using Eqs. (14)–(20). Then we can compute the stored energy at configuration 2 as

If it exists, the bistable state with locally minimal energy happens before the assumed configuration 2, thus Eq. (38) overestimates the minimal energy after the non-rigid folding. However, because Eq. (34) also overestimates the peak energy, when we use Δ*U* = *U*_{1} − *U*_{2} to determine the existence of a metastable state, the error shall be reduced.

### 3.3 Determining Existence of Metastable State.

As illustrated in Fig. 6, if Δ*U* = *U*_{1} − *U*_{2} > 0, a local minimum of stored energy exists besides the initial configuration (i.e., the metastable state), and thus the Shrimp pattern exhibits bistability. This energy difference is affected by the geometry of the pattern and material properties. Since we assume isometric deformations, according to Eqs. (34) and (38), the material-associated parameter that influences the value of Δ*U* is $k0B/k0F$, the ratio between bending stiffness and folding stiffness.

Figures 7(a) and 7(b) show the variation of Δ*U* with respect to the pairs of geometric parameters (*α*_{F}, *α*_{R}), and (*a*, *d*). We observe that *α*_{R} and *a* have positive influence on Δ*U*, while *d* has negative influence. We also find that *α*_{F} does not have significant influence on Δ*U*. The effect of $k0B/k0F$ is plotted in Fig. 7(c) in pair with *α*_{R}. The emergence of the bistability is defined by the contour of Δ*U* = 0, which is plotted in Fig. 7(d) as a function of *α*_{R}, *d*, and $k0B/k0F$. When $k0B/k0F$ is larger, the Shrimp pattern exhibits bistability over a wider range of geometries.

## 4 Verifying Bistable Behavior by Numerical Simulations

To verify the analytical analysis about the bistable behavior of the Shrimp pattern unit cell, we conduct numerical simulations using the bar-and-hinge model. The implementation in this work is based on the N5B8 discretization scheme [31]. The basic idea of the N5B8 discretization is explained in Refs. [31,32]. We assume that Young’s modulus *E* = 1 GPa, Poisson’s ratio *ν* = 0.3, and sheet thickness *t* = 0.5 mm. The stiffness of folding hinges are assigned through scalable formulas as in Ref. [31]. We take the length scale factor *L** to be 30*L*_{F}, relative to the length of each folding hinge, which yields an average ratio between bending and folding stiffness per length $k0B/k0F\u224880$.

The results from numerical simulations of compressive folding are presented in Fig. 8. The boundary conditions for the numerical analysis are shown in Fig. 8(a). In this study, we fix some geometric parameters as $\alpha F=70deg$, *d* = 10 mm, *a* = 20 mm, and *b* = 25 mm. When $\alpha R>49deg$, the bistable behavior starts to appear. The displacement *u* measures the movement of node *O*_{8} in the direction of loading (i.e., −*y*). We normalize the displacement *u* by *L*_{2}, such that (1) when *u*/*L*_{2} = 1, the rear panels are approximately coplanar, close to configuration 1 in the analytical derivation and (2) when *u*/*L*_{2} = 2, the rear panels snap inside the front panels, close to configuration 2 in the analytical derivation.

The negative forces in Fig. 8(b) during the bistable snapping display a long (negative) plateau without a strong peak, which behaves quite differently from most bistable origami structures, including the hypar origami [26]. The lack of strong negative peak force indicates that the snapping of the Shrimp pattern is a gentle process. We observe a stiffening effect occurring immediately after the snapping, as indicated by the sharp increase of forces when *u*/*L*_{2} > 2 in Fig. 8(b). As expected in the analytical analysis, if bistability appears, the stored energy of the Shrimp pattern experiences a peak near *u*/*L*_{2} = 1 and a local minimum around *u*/*L*_{2} = 2 (see Fig. 8(c)). As shown in Fig. 8(d), the analytical estimations of the energy difference Δ*U* agree well the numerically computed energy barriers. We record the angles of *ψ*′ during the simulations, as shown in Fig. 8(e). Compared with Fig. 5, we can clearly see that *ψ*′ is kinematically constrained to be below 2*α*_{B} ($=90deg$).

The numerical simulations verify that the energy landscape of the Shrimp pattern is programable, through geometric variations of the rear panel angle *α*_{R}. As *α*_{R} increases, both the peak force and magnitude of energy barrier (i.e., Δ*U*) increase accordingly. Such monotonic relationship allows for easy interpretation and manipulation of the mechanical properties of the Shrimp pattern.

## 5 Metastable Phases and Phase Transition of the Shrimp Tessellated Material

The programable bistability of the Shrimp pattern unit cell allows us to create architected materials with metastable phases by tessellating the shrimp unit cells. The compatibility condition for tessellation is that: in the *x*-direction, angle *ψ*′ of one unit must match *ψ*″ of the next unit; and in the *y*-direction, angle *ϕ*, and tail length *d* of adjacent units must be equal. Figure 9(a) shows an example of such designs when we tessellate three shrimp units in the *x*-direction to form a chain.

The state of least energy (when *U* = 0) of the Shrimp tessellation occurs when all three units are in their eggbox mode. Denote *ψ*_{0} as angle *ψ* at the state of least energy. We consider three variants of this tessellation: (1) $\psi 0=66deg$, $\alpha R=53deg$; (2) $\psi 0=86deg$, $\alpha R=53deg$; and (3) $\psi 0=86deg$, $\alpha R,AB=55deg$, $\alpha R,BC=53deg$, and $\alpha R,CD=51deg$; while we fix $\alpha F=70deg$, $\alpha B=45deg$, *a* = 20 mm, *b* = 25 mm, and *d* = 10 mm for all three variants. To analyze the behavior of this chain of Shrimp units, we consider each unit as a nonlinear one-dimensional spring, as shown in the left half of Fig. 9(a). This is referred to as the chain model in Fig. 9. Owing to the additional kinematic constraints implicitly induced by neighbor units in the tessellation, we update the boundary conditions used in Sec. 4 by removing the support on vertex *O*_{6} as it seems to be over restrictive. The compatibility between adjacent units is enforced by requiring that the angular difference between *ψ*′ of one unit and *ψ*″ of the next unit to be less than 0.01*π*. This chain tessellation has one zero energy phase and seven metastable phases. We consider a scenario when the whole chain is compressed by an amount of *u*_{A} = 2 *L*_{R,AB}. Qualitatively, we expect that there are three possible metastable phases, encoded by (E,E,M), (E,M,E), and (M,E,E), as illustrated in the right part of Fig. 9(a). The three major DOFs of the system are given by *u*_{A}, *u*_{B}, and *u*_{C}, as shown in Fig. 9(a). Based on the displacement–energy relationship revealed by numerical analysis of each individual unit, we can draw the complete energy landscape of the tessellation to reveal its three metastable states, as shown in Figs. 9(b)–9(d) for the three different variants of the chain tessellation. We can see that when the three units are the same, as of variants (1) and (2), the energy landscape is symmetric. However, when the three units are different in geometry, the energy landscape becomes asymmetric, and the energy state associated with unit CD in its Miura mode has a higher energy compared to the other two metastable states. This is because unit CD has a smaller *α*_{R}, as well as *L*_{R}, and hence at *u*_{A} = 2 *L*_{R,AB}, the whole chain needs to be further compressed beyond unit CD’s metastable state, which causes the system to store extra elastic energy.

To verify the predictions made by the simplified chain model, we perform simulation on a full model considering all three units. We make sure that in the full model, each unit is subject to equivalent boundary conditions as in the single unit simulation. We first find the local equilibrium of the chain tessellation when unit AB is at its Miura mode. We compute an initial rigid origami configuration as unit AB in its Miura mode (configuration 2) and units BC and CD in their eggbox modes (configuration 2), which is approximately the geometry of the metastable state (E,E,M). However, the zero energy states of the folding angles are set to be the same as all three units in their eggbox mode. This leads to some unbalanced forces in the chain, which is then intermediately resolved at the first increment of the nonlinear structural analysis as the algorithm seeks for equilibrium at the end of each increment [28]. Next, a displacement load *u*_{B} is applied while *u*_{A} is fixed to transition the tessellation from metastable state (E,E,M) to (E,M,E). This simulation results in a path on the energy landscape from one local minimum to another, which is then compared with the chain model predictions, as shown on the right sides of Figs. 9(b)–9(d). We find that the chain model is quite accurate, especially near the metastable states. Deviating from the stable states, the chain model appears to be stiffer than the full mode, because its DOFs are much less than the full model.

In each of the different phases, the Shrimp tessellated material displays different mechanical properties. Figure 10 shows how Poisson’s effect behaves differently for the variant (1) at different phases. Moreover, as the three units are the same, some phases behave the same way. Therefore, we can only find four distinct branches on Poisson’s ratio diagram. However, variant (3) has three different units, and hence its eight phases all behave differently, as shown in Fig. 10(b). Such a feature allows us to reprogram the material after it is manufactured, through mechanical phase transitions (i.e., snap-through), as illustrated in Fig. 11.

## 6 Concluding Remarks

The Shrimp origami pattern enables architected metamaterials with programable metastable phases through mechanical phase transitions. The Shrimp pattern unit is created by modifying the geometry of a rigid foldable pattern to block a range of rigid folding kinematics, which can only be overcome by non-rigid origami deformation. The non-rigid deformation separates the rigid folding kinematics of a Shrimp pattern into two disconnected ranges. Thus, the programable Shrimp pattern can display both rigid origami and non-rigid origami behavior. We study the connection between the geometry of a shrimp unit cell and its bistable behavior. An analytical estimation function is derived to approximate the energy barrier between the two stable states of a shrimp unit cell. The bistable Shrimp unit can be tessellated to form architected metamaterials. We show an example of the Shrimp tessellated materials with multiple stable configurations, which are interpreted as metastable phases. To transition from one phase to another, the tessellation must overcome an energy barrier through mechanical snapping. In each of the different phases, the metamaterial displays different mechanical properties, characterized here by Poisson’s ratio (yet other properties can also be explored). In Appendix B, we discuss approaches to fabricate such non-developable patterns.

Future investigation includes conducting experiments on the Shrimp pattern to investigate the performance of samples provided by different manufacturing techniques. In addition, the 3D stacking of the Shrimp pattern should be explored towards applications as smart multi-phase cellular materials.

## Acknowledgment

We thank the support from the US National Science Foundation (NSF) through Grant No. 1538830, the Japan Science and Technology Agency (JST) PRESTO through Grant No. JPMJPR1927, and the Raymond Allen Jones Chair at Georgia Tech. We are grateful to Emily D. Sanders for her helpful comments and suggestions to this paper.

## Conflict of Interest

There are no conflicts of interest.

## Data Availability Statement

The authors attest that all data for this study are included in the paper. Data provided by a third party are listed in Acknowledgements.

### Appendix A: Generalization of the Shrimp Pattern

Varying the geometry of the Shrimp pattern unit cell allows it to tessellate curved chains, such as the shrimp shape shown in Fig. 1. This is done by creating angular deficit on each of the shrimp units in an assemblage. Here, we illustrate the strategy by looking at the flat folded states of the Shrimp patten, as shown in Fig. 12.

In the flat folded state of a standard Shrimp pattern unit cell (eggbox mode), the points *O*_{5} and *O*_{8} lie on a line parallel to the line passing through *O*_{1} and *O*_{4}. These two parallel lines can be regarded as concentric circles of infinite radius (Fig. 12(a)). When we generalize the radii of the two concentric circles from infinite to finite values, we can construct a Shrimp unit cell that leads to a curved tessellation. As shown in Fig. 12(b), in the generalized Shrimp pattern unit cell, the lines of *O*_{1}*O*_{5}, *O*_{4}*O*_{8} are perpendicular to the circles passing through *O*_{1}*O*_{4} and *O*_{5}*O*_{8}, instead of the straight lines in the standard version. Now both lines of *O*_{1}*O*_{5}, *O*_{4}*O*_{8} are pointing to the shared center of the two circles. We let the extended line of *O*_{2}*O*_{6} to pass through the center, and let this line to bisect the angle formed by lines of *O*_{1}*O*_{5} and *O*_{4}*O*_{8}. As in the standard version, *O*_{1}*O*_{2} is parallel to *O*_{5}*O*_{6}, while the panel angle *α*_{F} is still a free design parameter. Following this design protocol, the front panel *O*_{1}*O*_{2}*O*_{6}*O*_{5} becomes an obtuse trapezoid instead of a parallelogram. We also keep the rectangle shape of the tail panel *O*_{3}*O*_{4}*O*_{8}*O*_{7}. As a result, the rear panel *O*_{2}*O*_{3}*O*_{7}*O*_{6} becomes a general quadrilateral. The generalized Shrimp pattern unit cell forms an angular deficit *θ*, such that a 1D tessellation leads to a curvature of *θ*/*L*, where *L* = ‖*O*_{1}*O*_{4}‖. When all unit cells in a tessellation are folded flat into the Miura mode, the curvature disappears and the tessellation becomes straight, as shown in Fig. 12(b). The Miura mode flat folded state can be obtained from the eggbox mode by two reflections about *O*_{2}*O*_{6} and *O*_{3}*O*_{7}.

### Appendix B: Fabrication of the Shrimp Pattern

The Shrimp pattern contains non-developable vertices, thus it cannot be folded from a single piece of flat sheet, like the Miura-ori. We fabricate the Shrimp pattern by gluing thin sheets or 3D printing.

We adopt a double layer technique to glue stripes of thin sheets, which is demonstrated in Fig. 13. The advantage of this strategy is that we obtain uniform thickness over the panels (double layers) and creases (single layer). The creases are treated by perforation or etching to reduce its folding stiffness. As we can observe from Fig. 7, to enable bistability, the Shrimp pattern needs to have a large $k0B/k0F$ ratio. By making the panels thicker, the double layer technique further increases the ratio of $k0B/k0F$.

Additive manufacturing technologies provide another avenue to realize the Shrimp pattern design. An automatic algorithm is written in Rhino/Grasshopper to generate a printable solid given an origami design. The printable solid is thicker within panels and thinner along folding hinges in order to differentiate their stiffness. The thin hinges and thick panels can also be printed with different materials. The procedure is outlined in Fig. 14.