The mechanics of a multistable metamaterial that is comprised of a series of bistable unit cells can be derived from the model for a single unit cell. Consider the metamaterial with *n* unit cells in series (Fig. 2(a)) and the corresponding theoretical model shown in Fig. 2(b). In the theoretical model, the beams are connected by a rigid link. The bottom beam (beam 1) has clamped–clamped boundary conditions; the left and right sections of the other beams (beams 2, 3,…, *n*) are only allowed vertical (*y*-direction) displacement. The thicknesses and mode imperfection sizes of the curved part of each unit cell are *t*_{i} and *a*_{3,}_{i}, respectively. For all unit cells, the initial height of the curved parts is *h*. The displacements of the midpoint of the curved part of unit cell *i* is *u*_{i} (see Fig. 2(b)); the deformations of the unit cells can be expressed as *d*_{1} = *u*_{1}, $di=ui\u2212ui\u22121$ for *i* = 2,…, *n*, and the total deformation is $dtot=\u2211i=1ndi$. The total deformation, *d*_{tot}, the total strain energy, *U*_{tot}, and the loading force, *f*, are normalized with respect to the parameters of unit cell 1, i.e., the following definitions are used for the normalized force, $f\xaf$, the normalized deformation, $d\xaftot$, and the normalized strain energy, $U\xaftot$: