Research Article

## Folding of Tubular Waterbomb

## Abstract

Origami has recently emerged as a promising building block of mechanical metamaterials because it offers a purely geometric design approach independent of scale and constituent material. The folding mechanics of origami-inspired metamaterials, i.e., whether the deformation involves only rotation of crease lines (rigid origami) or both crease rotation and facet distortion (nonrigid origami), is critical for fine-tuning their mechanical properties yet very difficult to determine for origami patterns with complex behaviors. Here, we characterize the folding of tubular waterbomb using a combined kinematic and structural analysis. We for the first time uncover that a waterbomb tube can undergo a mixed mode involving both rigid origami motion and nonrigid structural deformation, and the transition between them can lead to a substantial change in the stiffness. Furthermore, we derive theoretically the range of geometric parameters for the transition to occur, which paves the road to program the mechanical properties of the waterbomb pattern. We expect that such analysis and design approach will be applicable to more general origami patterns to create innovative programmable metamaterials, serving for a wide range of applications including aerospace systems, soft robotics, morphing structures, and medical devices.

## 1. Introduction

Mechanical metamaterials are artificially designed structures that offer extreme and unusual, yet useful, mechanical properties determined by their structural and geometric configurations rather than only intrinsic material properties of their composing elements. Conventional mechanical metamaterials are often formed by quasi-1D rods or links [1–6]. Recently, origami has emerged as a promising building block of mechanical metamaterials with versatile functionalities and programmability [7–17], due to its capability of transforming a 2D crease pattern into a complex 3D sculpture, purely geometric traits independent of both scale and constituent materials, and ease of manufacturing [18, 19].

The mechanical properties of an origami metamaterial are primarily determined by its folding mechanics. When rigid patterns such as Miura-ori [7, 8] are utilized, the facets do not stretch or bend but only rotate about the creases, and the metamaterials behave like a kinematic mechanism during folding. On the other hand, nonrigid ones [15, 16] enable simultaneous crease rotation and facet distortion, resulting in structural deformation of the metamaterials. Difference in folding mechanics leads to distinct mechanical properties such as stiffness [15]. Here, we report a tubular waterbomb pattern with a transition between rigid origami motion and structural deformation, making it possible for programming the behavior of the derived metamaterial.

In origami, the waterbomb tube refers to an origami structure made from a crease pattern obtained by tessellation of the waterbomb bases. A typical waterbomb base is a six-crease pattern with two colinear mountain creases and four diagonal valley ones intersecting at a common vertex [20]. The base has been used to create many fascinating origami objects including the structure that is the focus of this article [21–23]. One of the most distinctive characteristics of the waterbomb tube is that it has a negative Poisson’s ratio: when compressed, both its length and radius get smaller. This has led to a number of notable practical applications such as an expandable medical stent graft [24], a transformable worm robot [25], and a deformable robot wheel [26]. Recently, the authors also obtained programmable stiffness and shape modulation in the waterbomb tube using a bar-and-hinge numerical model [27]. Despite that, its precise folding behaviors and mechanical properties have remained ambiguous. Therefore, in this paper, we aim to expose the exact folding mechanics of the waterbomb tube by means of kinematics and structural analysis.

## 2. Kinematic Modelling of Rigid Folding

Figure 1(a) illustrates the crease pattern of a waterbomb tube defined by four independent geometric parameters—width , sector angle , and the number of bases longitudinally, , and circumferentially, . When we join together the left and right edges of the pattern, we can obtain a waterbomb tube [21, 22]. We will illustrate the motion behavior of the waterbomb tube with a representative model. First, we create a waterbomb tube in the fully contracted configuration (Figure 1(b) ①), where the facets in the middle row collide. When we slightly expand the tube along its central axis, its radius increases as well as the length, and a uniform radius along the tube is obtained (we shall demonstrate later that such a configuration always exists) (Figure 1(b) ②). With further expansion, it develops a pineapple shape with closure at both ends (Figure 1(b) ③ and ④). Subsequently, it opens up its ends again (Figure 1(b) ⑤) and then regains a uniform radius (Figure 1(b) ⑥). After this, the tube can be only marginally deployed, and the change in shape is hardly noticeable.

To determine the folding mechanics, we first build a kinematic model for the waterbomb tube. The pattern has three distinct vertex groups: central vertex and edge vertices and , the motion of which can be modeled as a spherical linkage with three degrees of freedom (DOFs) [28, 29]. As such, the tube is a network of these linkages, leading to a multi-DOF system (details in the supplementary material, S2). Thus, we reduce the DOF by making the following assumptions of symmetry based on our observation in Figure 1(b). First, the motion of the tube is symmetric about the equatorial plane that passes through the middle of the tube and divides it into two identical top and bottom halves. Second, all bases that are circumferentially placed in the same row have identical folding behavior. Finally, each base moves in a plane-symmetric way about a plane (presented as the red dot-dash line in Figure 2(a)) passing through the two mountain creases and the central axis of the tube. We first discuss the case when (i.e., the number of rows in the tube) is odd. Figure 2(a) presents a strip out of the origami pattern forming a tube with an odd number of rows. The equatorial plane passes through the center of the middle row defined as row 0. According to the above assumptions, linkage is symmetric about the equatorial plane and the plane passing vertex and tube axis, which reduces its DOF to one. Subsequently, the overall DOF of the tube becomes one, and its motion, described by the respective dihedral angles of all the linkages, can be found out through kinematic analysis (details in the supplementary material, S3).

Consider a particular example with , , and . As shown in Figure 2(b), is defined as the dihedral angle between two triangular facets that pass the top mountain crease in linkage . Taking as the input, we can obtain the other five dihedral angles of the linkage by kinematic analysis. Subsequently, we can determine all dihedral angles in the tube by using these known dihedral angles as inputs for adjacent linkages. To depict the extent of deployment of the tube, we define as the folding angle between the two largest triangular facets of a base on row 0 and the following equation holds: . The nondimensional radii of the vertices, , are plotted against in Figure 2(c), together with five representative configurations of the tube during deployment. We can draw three conclusions from the result. First, there exist two particular configurations, II with and IV with . On those configurations, the radii of all the vertices and become equal and so do those of vertices (), see the red dots on the curves. In other words, all the bases take the same geometric form at either configuration, resulting in a tube of a uniform radius. Thus, we prove theoretically that it is possible to construct a uniform waterbomb tube out of rigid sheet materials, provided that the amount of prefolding is correct. Second, the rigid-foldable range of the tube is bounded by two values of , and . The lower bound corresponds to the compactly folded configuration I: , i.e., two triangular facets on both sides of the common crease overlap entirely, whereas the upper bound is associated with the most expanded configuration V in which the upper sides of the bases on row 1 form a regular hexagon with a side length . The supplementary material (S3.B.1) provides detailed derivations on how both bounds are obtained. Third, the uniform radius and the bound configurations divide the tube deployment process with distinct shapes, i.e., a pineapple shape with the largest radius attained at crease in row 0, such as configuration III with , when , and a dogbone shape with the smallest radius reached at crease , when and .

## 3. Mechanism-Structure-Mechanism Transition

When we increase the number of rows, existing rows will retain their motion in the original tube with and drive concurrently the motion of newly added ones. However, the upper sides of bases on the end rows will form an -sided regular polygon earlier, resulting in the termination of the motion. For example, if we increase to 7, remains to be 60° while decreases to 144.24°. Figure 3(a) presents the motion sequence of the tube, where configurations I and V correspond to the minimum and maximum values of , respectively. Figure 3(b) plots the dihedral angles between adjacent bases of row () against . Figure 3(c) gives the radii of the vertices , , and , , , and , *vs.*. As in the case of , the curves intersect at two points (marked by red dots), which indicates that the tube also has a uniform radius at configurations II and IV when . Comparing it to the tube with , it is found that the values of at these two configurations are exactly the same. In fact, it can be proven kinematically that the uniform radius configurations are independent of the number of rows (supplementary material, S3.B.2). We can intuitively imagine that when all bases in the tube are in identical shape, more rows can be added to the tube in a geometrically compatible way.

Furthermore, the tube undergoes nonrigid folding within a region of between ° at configuration III_{L} and at III_{R} (supplementary material, S3.B.3). This is clearly demonstrated by the intermediate configuration III in Figure 3(a), where the central vertices of the bases on the two end rows collide with each other. The reason is that the dihedral angle in this range, indicates interference among the facets on row 3, which is not permitted in rigid origami. This conclusion is further supported by the fact that in the same range as shown in Figure 3(c). Therefore, the tube has rigid origami motion only within two distinct regions of : and . At configurations III_{L} and III_{R}, the ends of the cylinder are closed, and the tube becomes a concealed volume. It has been shown that structural deformation is required for any change in a concealed volume [30]. Thus, in order to move a physical tube specimen from one rigid-foldable range to the other, the component sheet material has to deform, and the tube works as a structure instead of a mechanism. This kind of folding behavior is referred to as *mechanism-structure-mechanism transition*.

The mechanism-structure-mechanism transition leads to a dramatic variation in stiffness for waterbomb-based structures and metamaterials. Within the rigid origami regime, the tube has a low stiffness determined by the torsional stiffness of the creases. When entering the structural range, the stiffness will be significantly increased due to facet deformation taken place in the tube for its shape change. This feature is demonstrated through a structural analysis of a tube using the finite element method. The model had identical , , and with that in Figure 3(a), and the facets were set to be 2 orders stiffer than the creases to distinguish the deformation of these two components [31–33] (supplementary material, S4). The simulation started from near configurations III_{R} and terminated at just beyond configuration III_{L}. The folding process of the tube is shown in Figure 3(e), from which physical contact and deformation of the facets at the ends are clearly seen at configuration III. The elastic strain energy of the tube is plotted in Figure 3(f). To manifest the effect of facet deformation, a fictional mechanism motion was also simulated for the tube by allowing the facets to freely penetrate into each other. In the mechanism mode, most energy is stored in the creases and distributed almost linearly. In the structural mode, both facets and creases deform. The energy difference between the two acts as an indicator of the level of deviation from rigid origami motion. The larger it is, the more extra deformation is required to enable the tube to move. It is also worth pointing out that the energy difference is not strictly zero at configuration III_{L} where the tube resumes rigid folding. This is due to some localized residual deformation near the vertices in contact in the structural mode, without affecting the global folding behaviour of the tube. Compared with the fictional mechanism mode, a 110% higher maximum strain energy is required to fold the tube through the structural range, indicating a larger stiffness. The precise stiffness variation depends on the sheet material of the tube, but when that is known, it is possible to predict the external loading that will cause the tube to move from one rigid origami range to the other. Notice that the stiffness property of this tubular waterbomb can also be revealed with the bar-and-hinge model such as the Merlin code [31] and the Origami Contact Simulator [34], or the smooth fold model [35].

## 4. Programmability of Folding Mechanics

Using the established kinematic model, we can program the existence and range of the mechanism-structure-mechanism transition by varying geometric parameters of the pattern. Take and as an example. The radius of the vertex and *vs*. for various is presented in Figure 3(d), where we can find the transition occurs when . When , the tube conducts a pure rigid origami motion from to during deployment. Since and are also related to , the motion range shrinks with . If , the tube experiences only one curtailed rigid origami range due to physical interference. For instance, when , the rigid motion range is .

The folding behavior of the waterbomb tube can also be programmed by varying the number of bases in a row. To demonstrate this, consider and , but changing from 4 to 20. We find from the theoretical model that rigid folding occurs only when . The upper and lower limits, and , are plotted against in Figure 4(a), together with the transitional structural range (highlighted by red lines) if it exists. The result indicates that the mechanism-structure-mechanism transition consistently appears when . However, when is relatively large, in the current geometric setup, the mechanism motion range in the neighborhood of becomes very narrow, and practically, the tube can be considered reaching a stable configuration rather than another mechanism motion range. Another interesting phenomenon when is large is that the tube can form a spherical shape that was used as an origami wheel [36] or artificial muscles [37]. The deployment process of such a tube with , , and is presented in Figure 4(b). Our analysis shows that when reaches 120.75°, the two triangles in yellow and purple, respectively, meet and overlap at the shaded area, which ceases the rigid folding process before the tube ends are closed.

So far, we have only discussed the situation of being odd. Similar behavior exists when is even. The equatorial row of the tube no longer exists in this case, and the equatorial plane passes through the midpoints of the creases linking vertices and . However, the tube remains a single DOF system as the two directly adjacent rows above and below the equatorial plane must behave the same under symmetrical assumptions, and they subsequently drive the motion of the remaining rows. Moreover, the tube also has a pair of uniform radius configurations identical to its odd row counterpart with the same and , albeit differs (supplementary material, S3.C). It means that the uniform radius configurations of the tube are solely decided by parameters and . It is not related to .

## 5. Conclusions

We have uncovered the true folding mechanics of the tubular waterbomb and its dependence on pattern geometric parameters. Through a rigorous kinematic analysis, we have demonstrated that some waterbomb tubes are capable of rigid origami motion, whereas others will experience what we refer to as a mechanism-structure-mechanism transition. And a structural analysis has revealed a significant increase in stiffness when the tube transforms into the structural range. Furthermore, we have derived theoretically the correlation between the occurrence and range of the mechanism-structure-mechanism transition and geometric parameters of the pattern, making accurate programming of the mechanical properties readily achievable. Thus, this work can not only facilitate the development of mechanical metamaterials making use of the intriguing properties of the tubular waterbomb but also provide an analysis framework for novel programmable metamaterials with wide engineering applications such as soft robotics, morphing structures, and medical devices. To adopt our analytical model where symmetry assumption was made in these applications, the key is to maintain the symmetry of the pattern during motion. Considering that the thickness of the sheets in the waterbomb pattern cannot be ignored in physical applications, the analytical model should be adjusted to the thick-panel origami model [23]. The motion symmetry would be satisfied automatically by structural constraints introduced by the panel thickness.

## Conflicts of Interest

The authors declare no competing financial interests.

## Authors’ Contributions

Jiayao Ma and Huijuan Feng contributed equally to this work. J. Ma, Y. Chen, and Z. You initiated and directed the research. H. Feng and Y. Chen carried out the kinematic analysis. J. Ma and D. Hou carried out the numerical simulation. All the authors prepared all the figures and supplementary information. J. Ma, H. Feng, Y. Chen, and Z. You wrote the manuscript.

## Acknowledgments

This work was supported by the National Natural Science Foundation of China (Projects 51825503, 51721003, and 51575377) and the Air Force Office of Scientific Research (FA9550-16-1-0339). During the course of this research, Z. You was appointed as a visiting professor at Tianjin University.

## Supplementary Materials

*Supplementary 1. *Supplementary Information: derivations, figures, and equations.

*Supplementary 2. *Movie S1: animation of Figure 2(c).

*Supplementary 3. *Movie S2: animation of Figure 3(a).

*Supplementary 4. *Movie S3: animation of Figure 4(b).

## References

- Y. Chen, T. Li, F. Scarpa, and L. Wang, “Lattice metamaterials with mechanically tunable Poisson’s ratio for vibration control,”
*Physical Review Applied*, vol. 7, no. 2, article 024012, 2017. View at Publisher · View at Scopus · View at Google Scholar - X. Zheng, H. Lee, T. H. Weisgraber et al., “Ultralight, ultrastiff mechanical metamaterials,”
*Science*, vol. 344, no. 6190, pp. 1373–1377, 2014. View at Publisher · View at Scopus · View at Google Scholar - C. Coulais, D. Sounas, and A. Alù, “Static non-reciprocity in mechanical metamaterials,”
*Nature*, vol. 542, no. 7642, pp. 461–464, 2017. View at Publisher · View at Scopus · View at Google Scholar - C. Coulais, A. Sabbadini, F. Vink, and M. van Hecke, “Multi-step self-guided pathways for shape-changing metamaterials,”
*Nature*, vol. 561, no. 7724, pp. 512–515, 2018. View at Publisher · View at Scopus · View at Google Scholar - T. Frenzel, M. Kadic, and M. Wegener, “Three-dimensional mechanical metamaterials with a twist,”
*Science*, vol. 358, no. 6366, pp. 1072–1074, 2017. View at Publisher · View at Scopus · View at Google Scholar - L. Wu, Z. Dong, H. Du, C. Li, N. X. Fang, and Y. Song, “Bioinspired ultra-low adhesive energy interface for continuous 3D printing: reducing curing induced adhesion,”
*Research*, vol. 2018, article 4795604, 10 pages, 2018. View at Publisher · View at Google Scholar - M. Schenk and S. D. Guest, “Geometry of Miura-folded metamaterials,”
*Proceedings of the National Academy of Sciences of the United States of America*, vol. 110, no. 9, pp. 3276–3281, 2013. View at Publisher · View at Scopus · View at Google Scholar - Z. Y. Wei, Z. V. Guo, L. Dudte, H. Y. Liang, and L. Mahadevan, “Geometric mechanics of periodic pleated origami,”
*Physical Review Letters*, vol. 110, no. 21, article 215501, 2013. View at Publisher · View at Scopus · View at Google Scholar - S. Kamrava, D. Mousanezhad, H. Ebrahimi, R. Ghosh, and A. Vaziri, “Origami-based cellular metamaterial with auxetic, bistable, and self-locking properties,”
*Scientific Reports*, vol. 7, no. 1, article 46046, 2017. View at Publisher · View at Scopus · View at Google Scholar - S. Waitukaitis, R. Menaut, B. G. G. Chen, and M. van Hecke, “Origami multistability: from single vertices to metasheets,”
*Physical Review Letters*, vol. 114, no. 5, article 055503, 2015. View at Publisher · View at Scopus · View at Google Scholar - E. T. Filipov, T. Tachi, and G. H. Paulino, “Origami tubes assembled into stiff, yet reconfigurable structures and metamaterials,”
*Proceedings of the National Academy of Sciences of the United States of America*, vol. 112, no. 40, pp. 12321–12326, 2015. View at Publisher · View at Scopus · View at Google Scholar - V. Brunck, F. Lechenault, A. Reid, and M. Adda-Bedia, “Elastic theory of origami-based metamaterials,”
*Physical Review E*, vol. 93, no. 3, article 033005, 2016. View at Publisher · View at Scopus · View at Google Scholar - H. Fang, S. C. A. Chu, Y. Xia, and K. W. Wang, “Programmable self-locking origami mechanical metamaterials,”
*Advanced Materials*, vol. 30, no. 15, article 1706311, 2018. View at Publisher · View at Scopus · View at Google Scholar - H. Yasuda, C. Chong, E. G. Charalampidis, P. G. Kevrekidis, and J. Yang, “Formation of rarefaction waves in origami-based metamaterials,”
*Physical Review E*, vol. 93, no. 4, article 043004, 2016. View at Publisher · View at Scopus · View at Google Scholar - J. L. Silverberg, A. A. Evans, L. McLeod et al., “Using origami design principles to fold reprogrammable mechanical metamaterials,”
*Science*, vol. 345, no. 6197, pp. 647–650, 2014. View at Publisher · View at Scopus · View at Google Scholar - J. Ma, J. Song, and Y. Chen, “An origami-inspired structure with graded stiffness,”
*International Journal of Mechanical Sciences*, vol. 136, pp. 134–142, 2018. View at Publisher · View at Scopus · View at Google Scholar - S. Ren, J. Wang, C. Song et al., “Single-step organization of plasmonic gold metamaterials with self-assembled DNA nanostructures,”
*Research*, vol. 2019, article 7403580, 10 pages, 2019. View at Publisher · View at Google Scholar - N. Bassik, G. M. Stern, and D. H. Gracias, “Microassembly based on hands free origami with bidirectional curvature,”
*Applied Physics Letters*, vol. 95, no. 9, article 091901, 2009. View at Publisher · View at Scopus · View at Google Scholar - C. D. Onal, M. T. Tolley, R. J. Wood, and D. Rus, “Origami-inspired printed robots,”
*IEEE/ASME Transactions on Mechatronics*, vol. 20, no. 5, pp. 2214–2221, 2014. View at Publisher · View at Scopus · View at Google Scholar - S. Randlett,
*The Art of Origami: Paper Folding, Traditional and Modern*, E.P. Dutton, 1961. - S. Fujimoto and M. Nishiwaki,
*Sojo Suru Origami Asobi no Shotai (Invitation to Creative Origami Playing)*, 1982, Asahi Culture Center. - B. Kresling, “Plant "design": mechanical simulations of growth patterens and bionics,”
*Biomimetics*, vol. 3, pp. 105–120, 1996. View at Google Scholar - Y. Chen, H. Feng, J. Ma, R. Peng, and Z. You, “Symmetric waterbomb origami,”
*Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences*, vol. 472, no. 2190, article 20150846, 2016. View at Publisher · View at Scopus · View at Google Scholar - K. Kuribayashi, K. Tsuchiya, Z. You et al., “Self-deployable origami stent grafts as a biomedical application of ni-rich tini shape memory alloy foil,”
*Materials Science and Engineering: A*, vol. 419, no. 1-2, pp. 131–137, 2006. View at Publisher · View at Scopus · View at Google Scholar - C. D. Onal, R. J. Wood, and D. Rus, “An origami-inspired approach to worm robots,”
*IEEE/ASME Transactions on Mechatronics*, vol. 18, no. 2, pp. 430–438, 2012. View at Publisher · View at Scopus · View at Google Scholar - D. Y. Lee, J. S. Kim, S. R. Kim, J. S. Koh, and K. J. Cho, “The deformable wheel robot using magic-ball origami structure,” in
*Proceedings of the ASME 2013 International Design Engineering Technical Conferences and Computers and Information in Engineering Conference. Volume 6B: 37th Mechanisms and Robotics Conference*, Portland, Oregon, USA, 2013. View at Publisher · View at Scopus · View at Google Scholar - T. Mukhopadhyay, J. Ma, H. Feng et al., “Programmable stiffness and shape modulation in origami materials: emergence of a distant actuation feature,”
*Applied Materials Today*, vol. 19, article 100537, 2020. View at Publisher · View at Google Scholar - J. S. Dai and J. Rees Jones, “Mobility in metamorphic mechanisms of foldable/erectable kinds,”
*Journal of Mechanical Design*, vol. 121, no. 3, pp. 375–382, 1999. View at Publisher · View at Scopus · View at Google Scholar - C. H. Chiang,
*Kinematics of Spherical Mechanisms*, Krieger Publishing, 2000. - R. Connelly, I. Sabitov, and A. Walz, “The bellows conjecture,”
*Contributions to Algebra and Geometry*, vol. 38, no. 1, pp. 1–10, 1997. View at Google Scholar - K. Liu and G. H. Paulino, “Nonlinear mechanics of non-rigid origami: an efficient computational approach,”
*Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences*, vol. 473, no. 2206, article 20170348, 2017. View at Publisher · View at Scopus · View at Google Scholar - E. T. Filipov, K. Liu, T. Tachi, M. Schenk, and G. H. Paulino, “Bar and hinge models for scalable analysis of origami,”
*International Journal of Solids and Structures*, vol. 124, pp. 26–45, 2017. View at Publisher · View at Scopus · View at Google Scholar - J. L. Silverberg, J. H. Na, A. A. Evans et al., “Origami structures with a critical transition to bistability arising from hidden degrees of freedom,”
*Nature Materials*, vol. 14, no. 4, pp. 389–393, 2015. View at Publisher · View at Scopus · View at Google Scholar - Y. Zhu and E. T. Filipov, “An efficient numerical approach for simulating contact in origami assemblages,”
*Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences*, vol. 475, no. 2230, article 20190366, 2019. View at Publisher · View at Google Scholar - E. A. P. Hernandez, D. J. Hartl, E. Akleman, and D. C. Lagoudas, “Modeling and analysis of origami structures with smooth folds,”
*Computer-Aided Design*, vol. 78, pp. 93–106, 2016. View at Publisher · View at Scopus · View at Google Scholar - D. Y. Lee, S. R. Kim, J. S. Kim, J. J. Park, and K. J. Cho, “Origami wheel transformer: a variable-diameter wheel drive robot using an origami structure,”
*Soft Robotics*, vol. 4, no. 2, pp. 163–180, 2017. View at Publisher · View at Scopus · View at Google Scholar - S. Li, D. M. Vogt, D. Rus, and R. J. Wood, “Fluid-driven origami-inspired artificial muscles,”
*Proceedings of the National Academy of Sciences of the United States of America*, vol. 114, no. 50, pp. 13132–13137, 2017. View at Publisher · View at Scopus · View at Google Scholar