Issue 
Emergent Scientist
Volume 3, 2019
IPT 2018



Article Number  5  
Number of page(s)  8  
Section  Physics  
DOI  https://doi.org/10.1051/emsci/2019004  
Published online  04 June 2019 
Research Article
Origami launcher
^{1}
V. N. Karazin Kharkiv National University – Svobody Sq. 4, 61022 Kharkiv, Ukraine
^{2}
Department of Aerospace Engineering Sciences, University of Colorado at Boulder, 429 UCB, Boulder, CO 80309, USA
^{3}
A. Ya. Usikov Institute for Radiophysics and Electronics, National Academy of Sciences of Ukraine, 12 Proskura Str., Kharkiv 61085, Ukraine
^{*} email: polly.kofman@icloud.com
Received:
30
September
2018
Accepted:
5
April
2019
The article studies the elastic and locomotive properties of Miuraoritype paper origami. The mechanics of a single paper crease is studied experimentally, and its nonelastic properties turn out to be crucial. The entire origami construction is then described as a collection of individual creases, its capability to launch small objects is evaluated, and the equation of motion is found. Thus, the height of the launched ball is studied theoretically and experimentally as a function of governing parameters.
Key words: Origami / paper structures / elastic properties / Miuraori
© O. Ivashtenko et al., published by EDP Sciences, 2019
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
In our investigation, we will consider the properties of paper origami, designed to vertically launch small loads. This problem was proposed by the organizing committee of the X International Physicists’ Tournament. In its original version, we were to optimize a paper construction such as Miuraori to vertically launch a standard PingPong ball to a maximally possible height. We were also restricted to use only one uncut sheet of A4 paper (80 g/m^{2}). It had a fixed thickness, density and inherent anisotropy caused by the fibrous structure.
We define origami as a paper construction created by folding paper. The creation of creases is accompanied by damage to the material and plastic deformations. This process is irreversible, as can be seen from the photo (Fig. 1). However, once folded, origami may be subjected to further deformations that can be elastic in the vicinity of the stable position. Different types of these deformations are possible: bending, stretching, allround compression, etc. Although both stretching and allround compression can store a lot of energy, we will focus on the optimization of a onelayer Miuraori, which operates solely on bending deformations.
1.1 Design of the folding
Miuraori is one of the popular models known from the mechanical modeling of metamaterials [1–3]. The Miura fold is a rigid origami, i.e. each cell is flat and does not change either its shape or its linear dimensions during any transformations of the folding. That is why in our treatment we will consider only bending deformations in the edges and require that the energy stored in bending is much greater than the energy stored in other types of deformations. This principle also allows us to treat the construction as a sum of single folds. We note that multilayer structures will be disregarded because of stretching deformations appearing in them.
Standard Miuraori is a pattern consisting of equal parallelograms connected by folds (Fig. 2a and b). Besides a standard flat Miuraori (Fig. 2a) there are multiple isomorphic variations [4,5], leading to a wide range of so called Miuralike structures. We call the standard Miura pattern flat, because Miuralike structures in a folded state may have different shapes, e.g., cylindrical. To construct a mathematical model of the Miuraori folding pattern we parametrized it using 5 parameters: the length L_{1} and the height L_{2} of the paper sheet (which are fixed for our problem), the number N_{1} of horizontal folds, the number N_{2} of vertical folds, and the angle of construction α (which are variable). This set of parameters can describe the standard Miura pattern (Fig. 2b). Moreover, by varying angles and distances, one can construct many possible Miuralike structures (Fig. 2c) that can be nonhomogeneous or even nonperiodic.
2 Dead ends
In this section we try to evaluate the effectiveness of the origami launcher purely theoretically, by applying the linear theory of elasticity and avoiding the use of any phenomenological laws and experimentally determined dependencies. Due to the peculiarities of the paper structure and the presence of very large deformations, this theory could provide at most an estimate of the height reached by the ball.
From the linear theory of beams, we can write the expression for the elastic potential energy of one fold [6]: $$W=\frac{{\varphi}^{2}IE}{2l}.$$(1)
Here ϕ = π is the bending angle of the fold. E is the Young’s modulus of the paper. $l=\pi \frac{d}{2}$ denotes the length of the neutral line, with d being the paper thickness. $I=\frac{a{d}^{3}}{12}$ is the geometric moment of inertia of the paper’s crosssection, with $a=L{N}_{1}+\frac{{L}_{2}{N}_{2}}{\mathrm{sin}\alpha}$ being the total length of all the creases on the sheet of paper (Fig. 2). Substituting all these expressions into (1), we obtain: $$W=\left(L{N}_{1}+\frac{{L}_{2}{N}_{2}}{\mathrm{sin}\alpha}\right){d}^{2}\frac{E\pi}{12}.$$(2)
The height to which the ball rises can be estimated from the law of conservation of energy, by equating (2) to the kinetic energy of the system (14), thus obtaining the initial velocity of the ball v, and expressing the maximum height using the law of uniformly accelerated motion, ${h}_{\text{max}}=\frac{{v}^{2}}{2g}$. The answer is $${h}_{\text{max}}=\frac{({L}_{1}{N}_{1}\mathrm{sin}\alpha +{L}_{2}{N}_{2}){d}^{2}E\pi}{4\mathrm{sin}\alpha (3m+{m}_{0})g}\text{\hspace{0.17em}}.$$(3)
This expression correctly reproduces some features seen experimentally in Figure 8: linear dependence of h_{max} on N_{1} and N_{2}, monotonous decrease of h_{max} as a function of α. Still, assuming E ≈ 3 × 10^{9} Pa [7] and substituting the typical geometric properties from our experiments, we get the value h_{max} ≈ 1 km, which is 3 orders of magnitude too high. The reason for this discrepancy can be seen in Figure 1: the inner side of the fold does not compress but corrugates, which also causes a much smaller stretching of the outer side of the fold, with an E_{p} much smaller than (2). Moreover, the real curvature radius of the fold is much larger than the assumed value of $\frac{d}{2}$. To some extent, these deficiencies of the theory can be remedied by introducing into (1)–(2) an “effective thickness” much smaller than h and an “effective length” much larger than $l=\pi \frac{d}{2}$. The unavoidable need for such arbitrary fitting parameters and a scandalously large height reached by the ball, in our opinion, render the existence of a simple axiomatic theory highly implausible. We thus consider this unsuccessful attempt to be a motivation for the creation of the phenomenological theory, as presented in the following part of this article.
Fig. 1 Demonstration of modifications in the structure of the folded paper. (a) The first folding of an initially flat sheet of paper. One can see the corrugation of the internal layer. (b) The same sample bent backwards. Now the preliminary stretched layer forms many more creases. 
Fig. 2 The Miura pattern and parametrization. (a) Photo of the standard (flat) Miura origami. (b) The standard Miura scheme and its parameters: the length L_{1} and the height L_{2} of the sheet (fixed), the number N_{1} of horizontaland N_{2} of vertical creases, and the angle of construction α (variable). (c) The Miuralike structure scheme. All cells may be arbitrarily varied, but L_{1} and L_{2} are still fixed. 
3 Methods
3.1 Dynamics of a single paper crease
The properties of the paper are related to numerous elastic and frictional processes in its fibrous microstructure. The torque provided by the bent paper nonlinearly depends on the angle between the facets. Moreover there is a hysteresis (from the cycle of opening/closing the fold), so the properties of paper are processdepending. Moving layers are subjected to the internal friction that also depends on the motion rate. The physics behind all these values is very complicated and sensitive to the particular conditions of the experiment, such as the brand and the quality of the paper, the state of the environment, the folding method etc. That is why it is very difficult to study from microscopics the elastic properties of such a complex material as paper. However, it is possible to construct a simple phenomenological theory describing this system.
To do this, we decided firstly to study the dynamics of a single paper crease, to obtain its basic elastic characteristics from the experiment, and then to apply this result to describe the entire origami. Let us write the equation of motion for the system consisting of a single paper crease fixed by one face to the floor, with an additional load fixed to the edge of its moving face: $$I\ddot{\varphi}=M(\varphi ,\dot{\varphi}){M}_{\text{g}}(\varphi ),$$(4)
where I is the moment of inertia of the moving part of the paper and the load, $M(\varphi ,\dot{\varphi})$ is the effective torque due to the elastic and frictional forces in the crease, and M_{g} (ϕ) is the torque related to the gravity force acting on the load and the paper. In our approach, we assume that the effective torque in its turn consists of the static torque M_{0} and the dissipative term that depends on the angle ϕ between the two paper faces and the angular velocity ϕ,
Fig. 3 Scheme of the experimental technique with a single fold. The load (red rectangle) is put at the extremity of one face, while the other is clamped to the surface. Considering the faces to be rigid (indeed, the crease is much easier to bend that the undamaged paper, as we saw in the experiment), we track the temporal evolution ϕ(t) of the crease opening, starting from ϕ = 0. 
Here we assume that M_{0} is only a function of the angle ϕ, and that the dissipative term is proportional to the angular velocity $\dot{\varphi}$ (some analogy to a viscoelastic model) with the coefficient of proportionality μ depending in its turn on the angle ϕ. Substitutingthese expressions and treating the values of the moment of inertia and the gravitational term through the paper and the load parameters, we obtain the differential equation of motion: $$\begin{array}{l}\left(m{l}^{2}+{m}_{0}\frac{{L}^{2}}{3}\right)\ddot{\varphi}\\ \text{\hspace{1em}}={M}_{0}(\varphi )\mu (\varphi )\dot{\varphi}\left(ml+{m}_{0}\frac{L}{2}\right)g\mathrm{cos}\varphi ,\end{array}$$(6)
where m and l are respectively the mass and the lever arm of the load, while m_{0} and L are the mass and the length of the moving part of the paper. In this equation, the static torque M_{0} (ϕ) and the dynamical dissipative coefficient μ(ϕ) as functions of the angle are to be determined from the experimental study.
In our experiment (Fig. 3), we filmed the dynamics of a loaded paper crease, traced it using the program Tracker,^{1} thus obtaining the time dependence of the angle ϕ(t) for differentmasses of the load, numerically differentiated them and treated the angular velocity $\dot{\varphi}(t)$ and acceleration $\ddot{\varphi}(t)$ as functions of ϕ (Fig. 4). Before performing the numerical differentiation, the experimental data were subjected to smoothing, using spline interpolation. All the points in our experiment come from a single paper sample. These dependencies are reproducible, but they are very sensitive to any changes: paper should be taken from the same package, fold should be made in the same way by the same person (because different pressure during folding, etc. causes different damage to the material), environmental conditions such as humidity should also be the same, etc. Therefore, unfortunately, to repeat all the following computations in different circumstances one needs to obtain corresponding characteristics of a single bend once again.
We analyzed the obtained experimental data to find the empirical parameters of our model. For this purpose, we substituted $\dot{\varphi}(\varphi )$ and $\ddot{\varphi}(\varphi )$ into equation (6) and treated the effective torque in each case. When varying the mass of the load, the angular velocities are different for the same value of the angle. As such, we fixed the value of the angle, constructed the dependence of the effective torque on the angular velocity, and fitted these data using our linear assumption with M_{0} and μ as fitting parameters (Fig. 5a). The linearity of this dependence verifies our assumption regarding the viscoelastic form of the effective torque of $M(\varphi ,\dot{\varphi})$.
Repeating the same procedure for different fixed angles we received the static torque M_{0} (ϕ) and the dynamical dissipative coefficient μ(ϕ) as functions of the angle (Fig. 5b). After interpolation, these functions can be used to solve the equations of motion of complex origami structures consisting of a number of creases.
Fig. 4 Angular dependences of (a) the angular velocity and (b) the acceleration. We used the measured timedependences of the angle, the angular velocity and the angular acceleration to express the latter two as functions of the angle. Different curves correspond to different masses of the load. 
3.2 Operation of the construction as a whole
Relying on the obtained results for a single paper crease, we now can build a theoretical model for the operation of a complex origami structure. We treat the origami folding as a set of single creases, so we need only to sum up the work of all these creases, taking into consideration the construction properties. Let us consider the energy transformations of the system including the origami and the load
Fig. 5 Processing the experimental data. (a) Example of the dependence of the effective torque $M({\varphi}_{0},\dot{\varphi})$ on the angular velocity $\dot{\varphi}$ for the fixed opening angle ϕ_{0} = 0.4 radians. Variation of the effective torque at a fixed angle is achieved by changing the load mass. This dependence is fitted according to our assumption $M(\varphi ,\dot{\varphi})={M}_{0}(\varphi )+\mu (\varphi )\dot{\varphi}$ with M_{0} (ϕ) and μ(ϕ) as fitting parameters. So, from this linear fit we get M_{0}(ϕ_{0}) and μ(ϕ_{0}). The same procedure is then repeated for other values of ϕ_{0}. (b) The resulting empirical parameters of our model as functions of the angle ϕ of inclination of the paper. The upper plot is the static torque provided by the paper crease and the lower one is the dynamical dissipative parameter slowing down the motion of the paper. The latter is the coefficient of proportionality relating the effective torque and the angular velocity. It should be noted that in further calculations these values will be normalized according to the corresponding length of the crease. 
where ΔE_{p} and ΔE_{k} are respectively the increments of the potential and kinetic energy of both the paper construction and the load, and ΔE_{l} describes all other energy losses related to the load, e.g., air drag, rotation, etc. By a simple estimation (see Appendix A) it can be shown that these losses are negligible with respect to the other terms of the equation, so we will not take ΔE_{l} into consideration in our further computations.
Our task is to study the height of the load, so it would be reasonable to rewrite all of the equations in terms of the vertical coordinate x (Fig. 6a). To do this, we need to establish relations between the change of the coordinate x and that of each of the angles ϕ_{i}. The elementary work of the effective torque then reads $$\begin{array}{lll}\mathcal{M}(x,\dot{x})dx\hfill & =\hfill & {\sum}_{i,\text{\hspace{0.05em}}edges}{M}_{0}[{\varphi}_{i}(x)]\frac{d{\varphi}_{i}}{dx}dx\hfill \\ \hspace{0.17em}\hfill & \hfill & {\sum}_{i,\text{\hspace{0.05em}}edges}\mu [{\varphi}_{i}(x)]\dot{x}{\left(\frac{d{\varphi}_{i}}{dx}\right)}^{2}\text{}dx\text{\hspace{0.17em}},\hfill \end{array}$$(8)
where M_{0} and μ are numerical functions (Fig. 5b) normalized according to the length of the creases.
To find the coefficients of transformation $\frac{d{\varphi}_{i}}{dx}$, we need to express the angles ϕ_{i} in terms of the coordinate x. Firstly, we express the angles between the two particular facets in terms of the local displacement of one fold, and then we generalize this result to the entire construction with the total displacement x. In the case of a standard flat Miuraori, we have two different types of angles (Fig. 6b). We solve this geometric problem (see Appendix B) and obtain the following expressions: $${\varphi}_{1}=\mathrm{arcsin}\left[2\frac{x}{{L}_{2}}\text{\hspace{0.17em}}\sqrt[]{1\frac{{x}^{2}}{{L}_{2}^{2}}}\text{\hspace{0.17em}}\right],$$(9) $${\varphi}_{2}=\mathrm{arcsin}\left[\frac{2\mathrm{cos}\alpha \text{\hspace{0.17em}}\sqrt[]{\frac{{L}_{2}^{2}}{{x}^{2}}1}}{\frac{{L}_{2}^{2}}{{x}^{2}}{\mathrm{sin}}^{2}\alpha}\right].$$(10)
As can be seen, these expressions are related to the parameters of the particular pattern discussed in the previous section. Here, L_{2} is the height of the paper sheet and α is the construction angle (see Fig. 2b). For the standard flat Miura pattern we have only one type of cell, and hence only one construction angle. We also need to find the derivatives of the obtained relations. Differentiating (9) and (10) we get $$\begin{array}{lll}\frac{d{\varphi}_{1}}{dx}\hfill & =\hfill & \frac{2}{\sqrt[]{{L}_{2}^{2}{x}^{2}}},\hfill \\ \frac{d{\varphi}_{2}}{dx}\hfill & =\hfill & \frac{2\mathrm{cos}\alpha \text{\hspace{0.17em}}}{\sqrt[]{{L}_{2}^{2}{x}^{2}}\left(1\frac{{x}^{2}}{{L}_{2}^{2}}{\mathrm{sin}}^{2}\alpha \right)}.\hfill \end{array}$$
Using these expressions (9)–(12), $\mathcal{M}(x,\dot{x})$ can be treated in terms of x: $$\begin{array}{lll}\mathcal{M}(x,\dot{x})\hfill & =\hfill & {N}_{1}{L}_{1}(\text{}{M}_{0}[{\varphi}_{1}(x)]\frac{d{\varphi}_{1}}{dx}\text{}\mu [{\varphi}_{1}(x)]\dot{x}\text{}{\left(\text{}\frac{d{\varphi}_{1}}{dx}\text{}\right)}^{2}\text{})\hfill \\ \hspace{0.17em}\hfill & \hfill & +\frac{{N}_{2}{L}_{2}}{\mathrm{sin}\alpha}(\text{}{M}_{0}[{\varphi}_{2}(x)]\frac{d{\varphi}_{2}}{dx}\text{}\mu [{\varphi}_{2}(x)]\dot{x}\text{}{\left(\text{}\frac{d{\varphi}_{2}}{dx}\text{}\right)}^{2}\text{}),\hfill \end{array}$$(13)
where N_{1} and N_{2} are respectively the numbers of horizontal and vertical folds. The multipliers L_{1} and $\frac{{N}_{2}{L}_{2}}{\mathrm{sin}\alpha}$ come from calculating the length of the folds of a given type (see Fig. 2b).
In respect of the kinetic and potential energy terms E_{k} and E_{p}, they will consist of the energies of the load and of the origami itself. The expression for the total kinetic energy of the origami is complicated as it includes all of the inner motions of each particular crease. To simplify the calculations, we will neglect these inner motions and consider only the vertical velocity of the paper. Of course, this neglect may result in errors of our computations, but we accept it as the first approximation. Thus, the kinetic energy of the origami also can be easily treated in terms of the vertical motion.
We assume that the bottom of the origami remains at rest (we fix it to the table), its top moves at the speed of the ball v, and the vertical velocity of the paper changes linearly with height (in reality, it is not exactly so because of the increasing effective load on the lower folds due to the mass of the above paper). Then, averaging the squared velocity over the length of the origami, we get the following expression for the kinetic energy of the system: $${E}_{\text{k}}={E}_{\text{k},\text{ball}}+{E}_{\text{k},\text{origami}}=\frac{m{v}^{2}}{2}+\frac{1}{3}\frac{{m}_{0}{v}^{2}}{2},$$(14)
where m is the mass of the ball, m_{0} is the mass of the origami, andthe factor $\frac{1}{3}$ comes from the averaging. As for the potential energy, it can also be treated in terms of x as a sum of the two energies (using the same notations): $${E}_{\text{p}}={E}_{\text{p},\text{ball}}+{E}_{\text{p},\text{origami}}=mgx+{m}_{0}g\frac{x}{2}.$$(15)
Having the expressions for the kinetic energy (14), the potential energy (15), and the effective torque (13), as well as the relation for the energy transformations in the system (7), we can consider the energy change per unit time that provides a new equation of motion similar to (6): $$\left(m+\frac{{m}_{0}}{3}\right)\ddot{x}=\mathcal{M}(x,\dot{x})\left(m+\frac{{m}_{0}}{2}\right)g,$$(16)
where $\mathcal{M}(x,\dot{x})$ is given by (13).
Fig. 6 Miuraori launching a ball. (a) The construction disposed vertically with a ball on it. x is the vertical coordinate of the ball measured from the initial position of its center of mass, the origami being in the folded state. (b) A single crease within the entire construction. Each of the angles ϕ_{i} corresponds to one of the types of folds. For the case of a standard flat Miuraori presented in this figure, there are two different types of folds: near the side a and near the side b. 
Fig. 7 Timedependence of the vertical velocity ẋ of the load obtained from the numerical solutions of the equation (16) for different numbers N_{1} of horizontaland N_{2} of vertical folds. Green lines show N_{1} varying from2 to 20 with fixed N_{2} = 10, blue lines – vice versa. Naturally, when increasing the number of folds, the height of the ball also increases. The curves are cut off at the point of maximum, that is in the moment of separation from the construction. 
4 Results
The numerical solution of this equation gives the time dependence of the coordinate of the load and its vertical velocity (Fig. 7). It is influenced by all of the parameters of our system. As we recall, our task was to investigate the maximum height of the launched load, so we should find the velocity of the ball at the moment of its separation from the origami. Then we could easily convert it into the maximum potential energy and find the height of the ball. From our solutions, the maximal velocity of the ball can be obtained as a function of the governing parameters N_{1}, N_{2}, and the constructionangle α.
4.1 Experiment
To study the dependence of the maximum height of the ball on the proposed parameters, we carried out a set of experiments varying independently N_{1}, N_{2} and α. We drew and printed three series of patterns with two parameters fixed and the third one varying in some range. All of the obtained dependencesare shown in Figure 8 and compared to the corresponding theoretical curves. As can be seen from these plots, the height of the launched load increases with increasing numbers of horizontal and vertical folds. This result is quite logical as the lever arm of the folds that decreases with increasing number of folds is not important for the resulting speed of the load. The energy coming from the creases is determined by their number and their inner properties. As concerns the impact of the constructionangle (Fig. 8c), its increase results in a decrease of the height. Indeed, in this case the projection on the vertical axis of the force coming from the vertical creases becomes smaller and decreases the energy given to the load.
It should be noted that fixing some typical values of the parameters in the datasets does not allow us to reach the maximal height. To do this, we should perform the optimization while involving all of the parameters simultaneously. From our result, one can see that to increase the height of the ball it is necessary to increase the numbers of folds and to decrease the construction angle. Experimentally, we are naturally limited by our possibility to construct an origami with arbitrary parameters. When trying to reach the optimal parameters, the height approaches a value of the order of 50 cm. Our best result was achieved when we tried to make an origami corresponding to the extreme values of the parameters from our plots (Fig. 8). Here, the load becomes unstable on top of the Miuraori, so the vertical launch becomes complicated, and the full collection of statistics becomes cumbersome. Moreover, when the cells are so small, the impact of the human factor dramatically increases. The characteristic scale of construction errors becomes comparable to the size of the cells, so each inaccuracy results in significant changes to the launcher properties.
Fig. 8 Experimental dependences of the maximum height of the launched ball on the governing parameters: (a) on the number N_{1} of horizontalfolds (N_{2} = 7, α = 60°); (b) on the number N_{2} of vertical folds (N_{1} = 5, α = 60°); (c) on the construction angle α (N_{1} = 7, N_{2} = 7). Error bars in the plots are estimated approximately, to include the instrumental errors, the statistical distribution of repeated experiments, and the measurement inaccuracies related to the human factor. Each experiment was repeated several dozen times, and then from the initial data set only vertical launches were selected. Each data set is compared to the theoretical solution with the corresponding parameters. 
5 Discussion
It is impossible to disregard the fact that the standard flat Miuraori studied in our paper is not the best construction for vertically launching a ball. We investigated it in order to understand the general principles of its work and to study the influence of different factors.
Our method was based on experimentally measured paper characteristics. That is why all of our experiments were carried out with the same batch of paper. When we tried to take another A4 paper, we got a strikingly different result: the same folding launches a ball almost two times higher. So, our theoretical functions cannot be applied to describe another type of paper. However, our treatment allows us to introduce arbitrary initial curves and parameters, and to perform all of the calculations again. Moreover, we can use the relations between the angles and the vertical coordinate to describe other Miuralike foldings. In general, the proposed method is applicable to the investigation of any kind of structure that can be treated as a sum of individual creases.
For example, the simplest extension of a standard Miuraori is the cylindrical Miura folding (Fig. 9). The properties of such constructions have been studied in detail in [8]. This type of pattern is particularly interesting for our problem because it provides a stable vertical launch of the ball and allows one to effectively use the energy stored in the folds.
The geometry of the cylindric folding is more complex than that of the flat one. In order to be cylindric, it already includes two different construction angles related with each other and the number of facets of the polygonal cylinder base. Thereby, the number of different angles between the pairs of faces and their derivatives with respect to the vertical coordinate increases. These relations become more complicated and imply some geometric restrictions caused by the ability of the folding to be easily compressed and decompressed. For example, for some values of angles, it is impossible to compress the cylinder by a continuous movement. It requires a kind of “geometry breaking”: to be compressed, the folding needs to be selfcrossed, like the Möbius strip. As we are launching a ball, it is an important limitation in our problem, together with other one concerning breaking the rigidity of some cylindrical constructions. That is, having two stable states, the origami needs its faces to be bent when getting from one stable state to another. These elastic modes, as well as other properties of the cylindrical Miuraori are discussed in [8].
This type of folding is able to provide a sufficiently large height of the ball (up to about 60 cm, which is slightly higher than our results for a flat Miuraori) and can be equally well described by our theoretical model with the corresponding geometry.
Fig. 9 The cylindrical Miura folding. (a) Shape of the construction with the polygonal cylinder base. (b) Scheme of this folding. Onecan notice that there are two different construction angles. These angles are related to each other and the number of edges of the polygonal base. 
6 Conclusion
The elastic properties of the origami arise from the combination of its folds, so we first study the mechanics of an individual crease. Theoretically and experimentally, we determined the dynamics of a paper fold. In our experiments we noticed that paper has nonelastic and nonlinear properties. We saw the hysteresis loop on the static loading curve when opening and closing the fold. We experimentally measured the dependence of the opening angle, the angular velocity and the angular acceleration on time. Thereby, we found the torque provided by the paper fold and the dissipation in the fold depending on the angle and the angular velocity. Using this result, we solved the differential equation of motion and measured the dependence of the maximal angular velocity on the mass of the load.
With the aid of the result for a single fold, we found the velocity of the entire origami, the velocity of the ball and the height reached by the ball. To do this, we created an algorithm to calculate the angles between the facets of a Miuraori structure, and implemented it numerically.
In addition, we have made a number of experimental measurements, in which we studied the velocity and the coordinate reached by the ball as a function of variable parameters, such as the number of horizontal and vertical folds, and the angle of the constructionof the origami. Comparing the experiment data to the numerical solution, we saw that the discrepancy was small. Still, if we use paper with other characteristics, its rigid properties will be different, and we will need to go through the fitting and calibration process once again.
Our theoretical and experimental methods devised for the Miuraori can be used to study other similar constructions.
Acknowledgments
We are very grateful to our friend Alexander Kostenko for revising the style of the manuscript, as well as to the referees and editors of the Emergent Scientist, whose comments helped to substantially improve the quality of the article.
Appendix A Estimation of the energy loss of the load
The load loses energy by interacting with air. With typical orders of magnitude of all the values in our problem, the estimated value of the Reynolds number of a load is $$\text{Re}=\frac{\rho vD}{\eta}~\frac{1\frac{\text{kg}}{{\text{m}}^{3}}1\frac{\text{m}}{\text{s}}{10}^{2}\text{m}}{{10}^{5}\frac{\text{kg}}{\text{ms}}}~{10}^{3},$$(A.1)
so we will use the quadratic law to estimate the air drag force F_{dr}. The energy dissipation is equal to the work of this force over the height of the launched load (the load is considered to be a ball) $${E}_{\text{dis}}~\overline{{F}_{\text{dr}}}h~\frac{C\rho \pi {r}^{2}\overline{{v}^{2}}h}{2}~\frac{C\rho \pi {r}^{2}g{h}^{2}}{3}.$$(A.2)
The potential energy gain of the load is $${E}_{\text{p}}=mgh.$$(A.3)
The ratio of these two energies is (units of measurement are omitted, all the values are written out in the SI) $$\begin{array}{ll}\frac{{E}_{\text{dis}}}{{E}_{\text{p}}}\hfill & ~\frac{C\rho \pi {r}^{2}g{h}^{2}}{3mgh}~\frac{C\rho \pi {r}^{2}h}{3m}\hfill \\ \hfill & ~\frac{0.5\cdot \pi \cdot 1\cdot {10}^{4}\cdot 0.5}{3\cdot {10}^{3}}~0.03\hfill \end{array}$$(A.4)
Thus, the energy dissipated in the air is two orders of magnitude less than the typical energies in our problem. Finally, to estimate the energy of rotation of the load, we assume the load to be a ball making about one turn over the raising time, as was observed in our experiments. $$\begin{array}{ll}{E}_{\text{rot}}=\hfill & \frac{I{\omega}^{2}}{2}~\frac{m{r}^{2}{(2\pi /t)}^{2}}{3}\hfill \\ \hfill & ~\frac{m{r}^{2}{(2\pi )}^{2}/(2h/g)}{3}~\frac{2{\pi}^{2}m{r}^{2}g}{3h},\hfill \end{array}$$(A.5)
so the ratio of this energy and the kinetic energy of the translational motion (which is almost equal to its potential energy gain) is (values are also given in SI) $$\begin{array}{ll}\frac{{E}_{\text{rot}}}{{E}_{\text{p}}}\hfill & ~\frac{2{\pi}^{2}m{r}^{2}g}{3mg{h}^{2}}=\frac{2{\pi}^{2}{r}^{2}}{3{h}^{2}}\hfill \\ \hfill & ~\frac{2{\pi}^{2}\cdot {10}^{4}}{3\cdot {0.5}^{2}}~\mathrm{0.003.}\hfill \end{array}$$(A.6)
Thus, the energy lost to the load rotation is also small.
Appendix B Calculation of opening angles
To find the opening angles as functions of the vertical coordinate x let us consider a paper sheet with the length L_{1} and the height L_{2} folded using the standard Miura pattern with N_{1} horizontal creases (Fig. B.1). As the origami opens, the net displacement variable x increases (x = 0 for the completely closed origami). This total displacement is a sum of the local displacements Δx of each cell, and we consider them to be equal (Δx = x∕N_{1}). There are two different opening angles in this construction: ϕ_{1} near the “horizontal” crease and ϕ_{2} near the “vertical” crease.
Fig. B.1 A folded standard Miura sheet with the total displacement x distributed over N_{1} local displacements Δx. The opening angle ϕ_{1} is located forexample between faces 1 and 3, and ϕ_{2} is respectively between faces 2 and 3 (the numbers of the faces are indicated in yellow circles). 
Let us focus on one particular cell (Fig. B.2). It consists of parallelograms of sizes l_{2} and l_{1} (where l_{2} = L_{2}∕(N_{1} sinα)) and has the construction angle α. From the right triangle ΔA_{3}A_{4}B_{1} (Fig. B.2a) one can find A_{4}B_{1} = l_{2} sinα. The line segment A_{4}B_{2} shows the distance between twoparallel planes A_{1}A_{4}A_{5} and A_{2}A_{3}A_{6}, so it is equal to the local displacement Δx. The angle between the planes A_{1}A_{2}A_{3} and A_{2}A_{3}A_{6} is equal to a half of the opening angle ϕ_{1} and can be found from the right triangle ΔA_{4}B_{1}B_{2} (see also Fig. B.2b). So, $$\mathrm{sin}\frac{{\varphi}_{1}}{2}=\frac{\Delta x}{h\mathrm{sin}\alpha}=\frac{x}{{L}_{2}},$$(B.1)
which after a trigonometrical transformation results in equation (9).
From the right triangle ΔA_{1}A_{4}C_{2} with ∠A_{1}A_{4}C_{2} = α we have A_{1}C_{2} = lsinα, and A_{4}C_{2} = lcosα. At the horizontal plane (Fig. B.2c) let ∠A_{4}A_{3}B_{2} = γ. Then $\mathrm{sin}\gamma \text{}=\text{}\frac{\Delta x}{h}$. As A_{4} A_{5}  A_{3}A_{6}, we get ∠C_{1}A_{4}C_{2} = γ, so from the right triangle ΔC_{1}A_{4}C_{2} we obtain C_{1}C_{2} = lcosαtanγ. Finally, in the right triangle ΔA_{1}C_{2}C_{1} (Fig. B.2d) the angle ∠C_{1}A_{1}C_{2} is equal to one half of the opening angle ϕ_{2}, so one can find $$\begin{array}{ll}\mathrm{sin}\frac{{\varphi}_{2}}{2}\hfill & =\frac{l\mathrm{cos}\alpha \mathrm{tan}\gamma}{l\mathrm{sin}\alpha}=\frac{\mathrm{cot}\alpha \mathrm{sin}\gamma}{\sqrt{1{\mathrm{sin}}^{2}\gamma}}\hfill \\ \hfill & =\frac{\mathrm{cos}\alpha}{\sqrt{\frac{{L}_{2}^{2}}{{x}^{2}}{\mathrm{sin}}^{2}\alpha}},\hfill \end{array}$$(B.2)
and for the doubled argument one can find equation (10).
Fig. B.2 One Miura cell with the additional construction showing the two planes (a) and its separate elements (b, c, d). The planes A_{1}A_{4}A_{5} and A_{2} A_{3}A_{6} always remain parallel and in a real origami launcher are parallel to the ground. The local displacement Δx is the distance between these planes. 
References
 J.L. Silverberg, A.A. Evans, L. McLeod, R.C. Hayward, T. Hull, C.D. Santangelo, I. Cohen, Using origami design principles to fold reprogrammable mechanical metamaterials, Science 345, 647 (2014) [CrossRef] [Google Scholar]
 M. Schenk, S.D. Guest, Geometry of Miurafolded metamaterials, Proc. Natl. Acad. Sci. 110, 3276 (2013) [CrossRef] [Google Scholar]
 M. Schenk, S.D. Guest, Origami folding: a structural engineering approach, in Fifth International Meeting of Origami Science, Mathematics, and Education, 2011. [Google Scholar]
 W. Gilewski, J. Pelczynski, P. Stawarz, A comparative study of Origami inspired folded plates, Proc. Eng. 91, 220–225 (2014) [CrossRef] [Google Scholar]
 P. Sareh, S.D. Guest, Design of nonisomorphic symmetric descendants of the Miuraori, Smart Mater. Struct. 24, 085002 (2015) [CrossRef] [Google Scholar]
 S. Timoshenko, Strength of Materials. Part 1: Elementary Theory and Problems (1948), p. 296 [Google Scholar]
 A.R. Jones, An experimental investigation of the inplane elastic moduli of paper (1967), pp. 125–129 [Google Scholar]
 A. Reid, F. Lechenault, S. Rica, M. AddaBedia, Geometry and design of origami bellows with tunable response, Phys. Rev. E 95, 013002 (2017) [CrossRef] [Google Scholar]
Cite this article as: Oryna Ivashtenko, Polina Kofman, Oleksiy Golubov, Zakhar Maizelis. Origami launcher, Emergent Scientist 3, 5 (2019)
All Figures
Fig. 1 Demonstration of modifications in the structure of the folded paper. (a) The first folding of an initially flat sheet of paper. One can see the corrugation of the internal layer. (b) The same sample bent backwards. Now the preliminary stretched layer forms many more creases. 

In the text 
Fig. 2 The Miura pattern and parametrization. (a) Photo of the standard (flat) Miura origami. (b) The standard Miura scheme and its parameters: the length L_{1} and the height L_{2} of the sheet (fixed), the number N_{1} of horizontaland N_{2} of vertical creases, and the angle of construction α (variable). (c) The Miuralike structure scheme. All cells may be arbitrarily varied, but L_{1} and L_{2} are still fixed. 

In the text 
Fig. 3 Scheme of the experimental technique with a single fold. The load (red rectangle) is put at the extremity of one face, while the other is clamped to the surface. Considering the faces to be rigid (indeed, the crease is much easier to bend that the undamaged paper, as we saw in the experiment), we track the temporal evolution ϕ(t) of the crease opening, starting from ϕ = 0. 

In the text 
Fig. 4 Angular dependences of (a) the angular velocity and (b) the acceleration. We used the measured timedependences of the angle, the angular velocity and the angular acceleration to express the latter two as functions of the angle. Different curves correspond to different masses of the load. 

In the text 
Fig. 5 Processing the experimental data. (a) Example of the dependence of the effective torque $M({\varphi}_{0},\dot{\varphi})$ on the angular velocity $\dot{\varphi}$ for the fixed opening angle ϕ_{0} = 0.4 radians. Variation of the effective torque at a fixed angle is achieved by changing the load mass. This dependence is fitted according to our assumption $M(\varphi ,\dot{\varphi})={M}_{0}(\varphi )+\mu (\varphi )\dot{\varphi}$ with M_{0} (ϕ) and μ(ϕ) as fitting parameters. So, from this linear fit we get M_{0}(ϕ_{0}) and μ(ϕ_{0}). The same procedure is then repeated for other values of ϕ_{0}. (b) The resulting empirical parameters of our model as functions of the angle ϕ of inclination of the paper. The upper plot is the static torque provided by the paper crease and the lower one is the dynamical dissipative parameter slowing down the motion of the paper. The latter is the coefficient of proportionality relating the effective torque and the angular velocity. It should be noted that in further calculations these values will be normalized according to the corresponding length of the crease. 

In the text 
Fig. 6 Miuraori launching a ball. (a) The construction disposed vertically with a ball on it. x is the vertical coordinate of the ball measured from the initial position of its center of mass, the origami being in the folded state. (b) A single crease within the entire construction. Each of the angles ϕ_{i} corresponds to one of the types of folds. For the case of a standard flat Miuraori presented in this figure, there are two different types of folds: near the side a and near the side b. 

In the text 
Fig. 7 Timedependence of the vertical velocity ẋ of the load obtained from the numerical solutions of the equation (16) for different numbers N_{1} of horizontaland N_{2} of vertical folds. Green lines show N_{1} varying from2 to 20 with fixed N_{2} = 10, blue lines – vice versa. Naturally, when increasing the number of folds, the height of the ball also increases. The curves are cut off at the point of maximum, that is in the moment of separation from the construction. 

In the text 
Fig. 8 Experimental dependences of the maximum height of the launched ball on the governing parameters: (a) on the number N_{1} of horizontalfolds (N_{2} = 7, α = 60°); (b) on the number N_{2} of vertical folds (N_{1} = 5, α = 60°); (c) on the construction angle α (N_{1} = 7, N_{2} = 7). Error bars in the plots are estimated approximately, to include the instrumental errors, the statistical distribution of repeated experiments, and the measurement inaccuracies related to the human factor. Each experiment was repeated several dozen times, and then from the initial data set only vertical launches were selected. Each data set is compared to the theoretical solution with the corresponding parameters. 

In the text 
Fig. 9 The cylindrical Miura folding. (a) Shape of the construction with the polygonal cylinder base. (b) Scheme of this folding. Onecan notice that there are two different construction angles. These angles are related to each other and the number of edges of the polygonal base. 

In the text 
Fig. B.1 A folded standard Miura sheet with the total displacement x distributed over N_{1} local displacements Δx. The opening angle ϕ_{1} is located forexample between faces 1 and 3, and ϕ_{2} is respectively between faces 2 and 3 (the numbers of the faces are indicated in yellow circles). 

In the text 
Fig. B.2 One Miura cell with the additional construction showing the two planes (a) and its separate elements (b, c, d). The planes A_{1}A_{4}A_{5} and A_{2} A_{3}A_{6} always remain parallel and in a real origami launcher are parallel to the ground. The local displacement Δx is the distance between these planes. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.