IPT 2018
Open Access
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

© O. Ivashtenko et al., published by EDP Sciences, 2019

Licence Creative CommonsThis 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 Miura-ori to vertically launch a standard Ping-Pong ball to a maximally possible height. We were also restricted to use only one uncut sheet of A4 paper (80 g/m2). 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, all-round compression, etc. Although both stretching and all-round compression can store a lot of energy, we will focus on the optimization of a one-layer Miura-ori, which operates solely on bending deformations.

1.1 Design of the folding

Miura-ori is one of the popular models known from the mechanical modeling of metamaterials [13]. 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 Miura-ori is a pattern consisting of equal parallelograms connected by folds (Fig. 2a and b). Besides a standard flat Miura-ori (Fig. 2a) there are multiple isomorphic variations [4,5], leading to a wide range of so called Miura-like structures. We call the standard Miura pattern flat, because Miura-like structures in a folded state may have different shapes, e.g., cylindrical. To construct a mathematical model of the Miura-ori folding pattern we parametrized it using 5 parameters: the length L1 and the height L2 of the paper sheet (which are fixed for our problem), the number N1 of horizontal folds, the number N2 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 Miura-like structures (Fig. 2c) that can be non-homogeneous or even non-periodic.

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= ϕ 2 IE 2l . \begin{equation*} W = \frac{{\phi}^2 IE}{2l}.\end{equation*}(1)

Here ϕ = π is the bending angle of the fold. E is the Young’s modulus of the paper. l=π d 2 $l=\pi\frac{d}{2}$ denotes the length of the neutral line, with d being the paper thickness. I= a d 3 12 $I=\frac{a d^3}{12}$ is the geometric moment of inertia of the paper’s cross-section, with a=L N 1 + L 2 N 2 sinα $a=L N_1+ \frac{L_2 N_2}{\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=( L N 1 + L 2 N 2 sinα ) d 2 Eπ 12 . \begin{equation*} W = \left(L N_1+ \frac{L_2 N_2}{\sin \alpha}\right) d^2 \frac{E \pi}{12}.\end{equation*}(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 max = v 2 2g $h_{\textrm{max}}=\frac{v^2}{2g}$. The answer is h max = ( L 1 N 1 sinα+ L 2 N 2 ) d 2 Eπ 4sinα(3m+ m 0 )g . \begin{equation*} h_{\textrm{max}} = \frac{(L_1N_1\sin{\alpha}+L_2N_2)d^2E \pi}{4\sin{\alpha}(3m+m_0)g}\,. \end{equation*}(3)

This expression correctly reproduces some features seen experimentally in Figure 8: linear dependence of hmax on N1 and N2, monotonous decrease of hmax as a function of α. Still, assuming E ≈ 3 × 109 Pa [7] and substituting the typical geometric properties from our experiments, we get the value hmax ≈ 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 Ep much smaller than (2). Moreover, the real curvature radius of the fold is much larger than the assumed value of d 2 $\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=π d 2 $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.

thumbnail 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.

thumbnail 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 L1 and the height L2 of the sheet (fixed), the number N1 of horizontaland N2 of vertical creases, and the angle of construction α (variable). (c) The Miura-like structure scheme. All cells may be arbitrarily varied, but L1 and L2 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 micro-structure. 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 process-depending. 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 ϕ ¨ =M(ϕ, ϕ ˙ ) M g (ϕ), \begin{equation*} I\ddot{\phi} = M(\phi, \dot{\phi})-M_{\textrm{g}}(\phi),\end{equation*}(4)

where I is the moment of inertia of the moving part of the paper and the load, M(ϕ, ϕ ˙ ) $M(\phi, \dot{\phi})$ is the effective torque due to the elastic and frictional forces in the crease, and Mg (ϕ) 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 M0 and the dissipative term that depends on the angle ϕ between the two paper faces and the angular velocity ϕ,

thumbnail 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.

M(ϕ, ϕ ˙ )= M 0 (ϕ)+μ(ϕ) ϕ ˙ . \begin{equation*} M(\phi, \dot{\phi}) = M_0(\phi)+\mu(\phi)\dot{\phi}.\end{equation*}(5)

Here we assume that M0 is only a function of the angle ϕ, and that the dissipative term is proportional to the angular velocity ϕ ˙ $\dot{\phi}$ (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: ( m l 2 + m 0 L 2 3 ) ϕ ¨ = M 0 (ϕ)μ(ϕ) ϕ ˙ ( ml+ m 0 L 2 )gcosϕ, \begin{align*}&\left(ml^2+m_0\frac{L^2}{3}\right)\ddot{\phi}\nonumber\\ &\quad=M_0(\phi)-\mu(\phi)\dot{\phi}-\left(ml+m_0\frac{L}{2}\right)g\cos \phi, \end{align*}(6)

where m and l are respectively the mass and the lever arm of the load, while m0 and L are the mass and the length of the moving part of the paper. In this equation, the static torque M0 (ϕ) 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 ϕ ˙ (t) $\dot{\phi}(t)$ and acceleration ϕ ¨ (t) $\ddot{\phi}(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{\phi}(\phi)$ and ϕ ¨ (ϕ) $\ddot{\phi}(\phi)$ 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 M0 and μ as fitting parameters (Fig. 5a). The linearity of this dependence verifies our assumption regarding the viscoelastic form of the effective torque of M(ϕ, ϕ ˙ ) $M(\phi, \dot{\phi})$.

Repeating the same procedure for different fixed angles we received the static torque M0 (ϕ) 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.

thumbnail Fig. 4

Angular dependences of (a) the angular velocity and (b) the acceleration. We used the measured time-dependences 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

thumbnail Fig. 5

Processing the experimental data. (a) Example of the dependence of the effective torque M( ϕ 0 , ϕ ˙ ) $M(\phi_0, \dot{\phi})$ on the angular velocity ϕ ˙ $\dot{\phi}$ 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(ϕ, ϕ ˙ )= M 0 (ϕ)+μ(ϕ) ϕ ˙ $M(\phi,\dot{\phi})=M_0(\phi)+\mu(\phi)\dot{\phi}$ with M0 (ϕ) and μ(ϕ) as fitting parameters. So, from this linear fit we get M0(ϕ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.

i,edges M( ϕ i , ϕ i . )Δ ϕ i =Δ E p +Δ E k +Δ E l , \begin{align*} \sum_{i, \, edges} M(\phi_i, \dot{\phi_i})\Delta\phi_i=\Delta E_{\textrm{p}}+\Delta E_{\textrm{k}}+\Delta E_{\textrm{l}},\end{align*}(7)

where ΔEp and ΔEk are respectively the increments of the potential and kinetic energy of both the paper construction and the load, and ΔEl 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 ΔEl 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 M(x, x ˙ )dx = i,edges M 0 [ ϕ i (x)] d ϕ i dx dx i,edges μ[ ϕ i (x)] x ˙ ( d ϕ i dx ) 2 dx, \begin{eqnarray*} \mathcal{M}(x,\dot{x}) dx&=&\sum_{i, \,edges} M_0[\phi_i(x)]\frac{d\phi_i}{dx} dx\nonumber\\* &&-\sum_{i, \,edges} \mu[\phi_i(x)]\dot{x}\left(\frac{d\phi_i}{dx}\right)^2\!dx \,,\end{eqnarray*}(8)

where M0 and μ are numerical functions (Fig. 5b) normalized according to the length of the creases.

To find the coefficients of transformation d ϕ i dx $\frac{d\phi_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 Miura-ori, we have two different types of angles (Fig. 6b). We solve this geometric problem (see Appendix B) and obtain the following expressions: ϕ 1 =arcsin[ 2 x L 2 1 x 2 L 2 2 ], \begin{equation*} \phi_1 = \arcsin\left[2 \frac{x}{L_2}\; \sqrt[]{1-\frac{x^2}{L_2^2}}\,\right],\end{equation*}(9) ϕ 2 =arcsin[ 2cosα L 2 2 x 2 1 L 2 2 x 2 sin 2 α ]. \begin{equation*} \phi_2 = \arcsin\left[\frac{2 \cos\alpha\ \sqrt[]{\frac{L_2^2}{x^2}-1}}{\frac{L_2^2}{x^2}-\sin^2\alpha}\right].\end{equation*}(10)

As can be seen, these expressions are related to the parameters of the particular pattern discussed in the previous section. Here, L2 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 d ϕ 1 dx = 2 L 2 2 x 2 , d ϕ 2 dx = 2cosα L 2 2 x 2 ( 1 x 2 L 2 2 sin 2 α ) . \begin{eqnarray} \frac{d\phi_1}{dx} &=& \frac{2}{\sqrt[]{L_2^2-x^2}},\\\frac{d\phi_2}{dx} &=& \frac{2 \cos\alpha\ }{\sqrt[]{L_2^2-x^2}\left(1-\frac{x^2}{L_2^2}\sin^2\alpha\right)}.\end{eqnarray}

Using these expressions (9)–(12), M(x, x ˙ ) $\mathcal{M}(x,\dot{x})$ can be treated in terms of x: M(x, x ˙ ) = N 1 L 1 ( M 0 [ ϕ 1 (x)] d ϕ 1 dx μ[ ϕ 1 (x)] x ˙ ( d ϕ 1 dx ) 2 ) + N 2 L 2 sinα ( M 0 [ ϕ 2 (x)] d ϕ 2 dx μ[ ϕ 2 (x)] x ˙ ( d ϕ 2 dx ) 2 ), \begin{eqnarray*}\mathcal{M}(x,\dot{x})&=&N_1L_1\Bigg(\!M_0[\phi_1(x)]\frac{d\phi_1}{dx}\!-\mu[\phi_1(x)]\dot{x}\!\left(\!\frac{d\phi_1}{dx}\!\right)^2\!\Bigg)\nonumber\\* &&+\frac{N_2L_2}{\sin \alpha}\Bigg(\!M_0[\phi_2(x)]\frac{d\phi_2}{dx}\!-\mu[\phi_2(x)]\dot{x}\!\left(\!\frac{d\phi_2}{dx}\!\right)^2\!\Bigg),\nonumber\\ \end{eqnarray*}(13)

where N1 and N2 are respectively the numbers of horizontal and vertical folds. The multipliers L1 and N 2 L 2 sinα $\frac{N_2L_2}{\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 Ek and Ep, 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 k = E k,ball + E k,origami = m v 2 2 + 1 3 m 0 v 2 2 , \begin{equation*} E_{\textrm{k}}=E_{\textrm{k,ball}}+E_{\textrm{k,origami}}=\frac{m v^2}{2}+\frac{1}{3} \frac{m_0 v^2}{2},\end{equation*}(14)

where m is the mass of the ball, m0 is the mass of the origami, andthe factor 1 3 $\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 p = E p,ball + E p,origami =mgx+ m 0 g x 2 . \begin{align*} E_{\textrm{p}}=E_{\textrm{p,ball}}+E_{\textrm{p,origami}}=m g x+m_0 g\frac{x}{2}.\end{align*}(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): ( m+ m 0 3 ) x ¨ =M(x, x ˙ )( m+ m 0 2 )g, \begin{align*} \left(m+\frac{m_0}{3}\right)\ddot{x}=\mathcal{M}(x,\dot{x})-\left(m+\frac{m_0}{2}\right)g,\end{align*}(16)

where M(x, x ˙ ) $\mathcal{M}(x,\dot{x})$ is given by (13).

thumbnail Fig. 6

Miura-ori 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 Miura-ori presented in this figure, there are two different types of folds: near the side a and near the side b.

thumbnail Fig. 7

Time-dependence of the vertical velocity of the load obtained from the numerical solutions of the equation (16) for different numbers N1 of horizontaland N2 of vertical folds. Green lines show N1 varying from2 to 20 with fixed N2 = 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 N1, N2, 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 N1, N2 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 Miura-ori, 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.

thumbnail Fig. 8

Experimental dependences of the maximum height of the launched ball on the governing parameters: (a) on the number N1 of horizontalfolds (N2 = 7, α = 60°); (b) on the number N2 of vertical folds (N1 = 5, α = 60°); (c) on the construction angle α (N1 = 7, N2 = 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 Miura-ori 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 Miura-like 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 Miura-ori 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 self-crossed, 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 Miura-ori 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 Miura-ori) and can be equally well described by our theoretical model with the corresponding geometry.

thumbnail 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 Miura-ori 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 Miura-ori 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 Re= ρvD η ~ 1 kg m 3 1 m s 10 2 m 10 5 kg ms ~ 10 3 , \begin{equation*} \textrm{Re}=\frac{\rho v D}{\eta}\sim \frac{1\frac{\textrm{kg}}{{\textrm{m}^3}} 1 \frac{\textrm{m}}{\textrm{s}}10^{-2} \textrm{m}}{10^{-5}\frac{\textrm{kg}}{\textrm{m s}}}\sim10^3, \end{equation*}(A.1)

so we will use the quadratic law to estimate the air drag force Fdr. 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 dis ~ F dr ¯ h~ Cρπ r 2 v 2 ¯ h 2 ~ Cρπ r 2 g h 2 3 . \begin{equation*} E_{\textrm{dis}}\sim \overline{F_{\textrm{dr}}}h \sim\frac{C \rho \pi r^2 \overline{v^2} h}{2} \sim \frac{C\rho \pi r^2 g h^2}{3}. \end{equation*}(A.2)

The potential energy gain of the load is E p =mgh. \begin{equation*} E_{\textrm{p}}=m g h. \end{equation*}(A.3)

The ratio of these two energies is (units of measurement are omitted, all the values are written out in the SI) E dis E p ~ Cρπ r 2 g h 2 3mgh ~ Cρπ r 2 h 3m ~ 0.5π1 10 4 0.5 3 10 3 ~0.03 \begin{align*} \nonumber \frac{E_{\textrm{dis}}}{E_{\textrm{p}}}& \sim \frac{C \rho \pi r^2 g h^2}{3 m g h} \sim \frac{C \rho \pi r^2 h}{3 m}\\[6pt] &\sim \frac{0.5\cdot \pi \cdot 1\cdot 10^{-4}\cdot 0.5}{3\cdot 10^{-3}}\sim0.03 \end{align*}(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. E rot = I ω 2 2 ~ m r 2 (2π/t) 2 3 ~ m r 2 (2π) 2 /(2h/g) 3 ~ 2 π 2 m r 2 g 3h , \begin{align*} \nonumber E_{\textrm{rot}}=&\frac{I\omega^2}{2} \sim\frac{m r^2 (2 \pi/t)^2}{3}\\[6pt] &\sim\frac{m r^2 (2 \pi)^2/(2h/g)}{3} \sim\frac{2\pi^2 m r^2 g}{3h}, \end{align*}(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) E rot E p ~ 2 π 2 m r 2 g 3mg h 2 = 2 π 2 r 2 3 h 2 ~ 2 π 2 10 4 3 0.5 2 ~0.003. \begin{align*} \nonumber \frac{E_{\textrm{rot}}}{E_{\textrm{p}}}&\sim \frac{2\pi^2 m r^2 g}{3 m g h^2}=\frac{2\pi^2 r^2}{3h^2}\\[6pt] &\sim\frac{2\pi^2\cdot 10^{-4}}{3\cdot0.5^2}\sim0.003. \end{align*}(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 L1 and the height L2 folded using the standard Miura pattern with N1 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 = xN1). There are two different opening angles in this construction: ϕ1 near the “horizontal” crease and ϕ2 near the “vertical” crease.

thumbnail Fig. B.1

A folded standard Miura sheet with the total displacement x distributed over N1 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 l2 and l1 (where l2 = L2∕(N1 sinα)) and has the construction angle α. From the right triangle ΔA3A4B1 (Fig. B.2a) one can find A4B1 = l2 sinα. The line segment A4B2 shows the distance between twoparallel planes A1A4A5 and A2A3A6, so it is equal to the local displacement Δx. The angle between the planes A1A2A3 and A2A3A6 is equal to a half of the opening angle ϕ1 and can be found from the right triangle ΔA4B1B2 (see also Fig. B.2b). So, sin ϕ 1 2 = Δx hsinα = x L 2 , \begin{equation*} \sin{\frac{\phi_1}{2}}=\frac{\Delta x}{h \sin{\alpha}}=\frac{x}{L_2}, \end{equation*}(B.1)

which after a trigonometrical transformation results in equation (9).

From the right triangle ΔA1A4C2 with ∠A1A4C2 = α we have A1C2 = lsinα, and A4C2 = lcosα. At the horizontal plane (Fig. B.2c) let ∠A4A3B2 = γ. Then sinγ= Δx h $\sin{\gamma}\!=\!\frac{\Delta x}{h}$. As A4 A5 || A3A6, we get ∠C1A4C2 = γ, so from the right triangle ΔC1A4C2 we obtain C1C2 = lcosαtanγ. Finally, in the right triangle ΔA1C2C1 (Fig. B.2d) the angle ∠C1A1C2 is equal to one half of the opening angle ϕ2, so one can find sin ϕ 2 2 = lcosαtanγ lsinα = cotαsinγ 1 sin 2 γ = cosα L 2 2 x 2 sin 2 α , \begin{align*} \nonumber\sin{\frac{\phi_2}{2}}&=\frac{l \cos{\alpha}\tan{\gamma}}{l\sin{\alpha}}=\frac{\cot{\alpha}\sin{\gamma}}{\sqrt{1-\sin^2\gamma}}\\ &=\frac{\cos{\alpha}}{\sqrt{\frac{L_2^2}{x^2}-\sin^2\alpha}}, \end{align*}(B.2)

and for the doubled argument one can find equation (10).

thumbnail Fig. B.2

One Miura cell with the additional construction showing the two planes (a) and its separate elements (b, c, d). The planes A1A4A5 and A2 A3A6 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

  1. 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]
  2. M. Schenk, S.D. Guest, Geometry of Miura-folded metamaterials, Proc. Natl. Acad. Sci. 110, 3276 (2013) [CrossRef] [Google Scholar]
  3. M. Schenk, S.D. Guest, Origami folding: a structural engineering approach, in Fifth International Meeting of Origami Science, Mathematics, and Education, 2011. [Google Scholar]
  4. W. Gilewski, J. Pelczynski, P. Stawarz, A comparative study of Origami inspired folded plates, Proc. Eng. 91, 220–225 (2014) [CrossRef] [Google Scholar]
  5. P. Sareh, S.D. Guest, Design of non-isomorphic symmetric descendants of the Miura-ori, Smart Mater. Struct. 24, 085002 (2015) [CrossRef] [Google Scholar]
  6. S. Timoshenko, Strength of Materials. Part 1: Elementary Theory and Problems (1948), p. 296 [Google Scholar]
  7. A.R. Jones, An experimental investigation of the in-plane elastic moduli of paper (1967), pp. 125–129 [Google Scholar]
  8. A. Reid, F. Lechenault, S. Rica, M. Adda-Bedia, 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

thumbnail 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
thumbnail 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 L1 and the height L2 of the sheet (fixed), the number N1 of horizontaland N2 of vertical creases, and the angle of construction α (variable). (c) The Miura-like structure scheme. All cells may be arbitrarily varied, but L1 and L2 are still fixed.

In the text
thumbnail 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
thumbnail Fig. 4

Angular dependences of (a) the angular velocity and (b) the acceleration. We used the measured time-dependences 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
thumbnail Fig. 5

Processing the experimental data. (a) Example of the dependence of the effective torque M( ϕ 0 , ϕ ˙ ) $M(\phi_0, \dot{\phi})$ on the angular velocity ϕ ˙ $\dot{\phi}$ 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(ϕ, ϕ ˙ )= M 0 (ϕ)+μ(ϕ) ϕ ˙ $M(\phi,\dot{\phi})=M_0(\phi)+\mu(\phi)\dot{\phi}$ with M0 (ϕ) and μ(ϕ) as fitting parameters. So, from this linear fit we get M0(ϕ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
thumbnail Fig. 6

Miura-ori 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 Miura-ori presented in this figure, there are two different types of folds: near the side a and near the side b.

In the text
thumbnail Fig. 7

Time-dependence of the vertical velocity of the load obtained from the numerical solutions of the equation (16) for different numbers N1 of horizontaland N2 of vertical folds. Green lines show N1 varying from2 to 20 with fixed N2 = 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
thumbnail Fig. 8

Experimental dependences of the maximum height of the launched ball on the governing parameters: (a) on the number N1 of horizontalfolds (N2 = 7, α = 60°); (b) on the number N2 of vertical folds (N1 = 5, α = 60°); (c) on the construction angle α (N1 = 7, N2 = 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
thumbnail 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
thumbnail Fig. B.1

A folded standard Miura sheet with the total displacement x distributed over N1 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
thumbnail Fig. B.2

One Miura cell with the additional construction showing the two planes (a) and its separate elements (b, c, d). The planes A1A4A5 and A2 A3A6 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 (full-text 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 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.