Original Research ARTICLE
Curvature Tuning in Folded Strips Through Hyperstatic Applied Rotations
- Japan Agency for Marine-Earth Science and Technology, Department of Mathematical Science and Advanced Technology, Yokohama Institute for Earth Sciences, Yokohama, Japan
Folding a strip of paper generates extremely localized plastic strains. The relaxation of the residual stresses results in a ridge that joins two flat faces at an angle known as the dihedral angle. When constrained isostatically, the strip will be at its undeformed roof-like state. Instead, if confined hyperstatically, the flat faces will undergo bending. We demonstrate that the generated curvatures can change their sign with appropriate rotations applied at the ends. We use Euler's theory of the Elastica and a shooting method to match the applied rotations at the boundaries. We also consider a constitutive model for the discontinuous rotation that takes into account the initial dihedral angle and the rotational stiffness of the fold. We show that the curvatures on the left and the right of the fold change according to a law also confirmed by the Euler-Bernoulli beam theory for small displacements and rotations. For opposite applied rotations, the fold disappears in the limit of zero rotational stiffness; instead, for applied rotations of the same sign, there exists a theoretical non-zero critical rotational stiffness that neutralizes the fold. Below such critical value, the fold can mutate, for example, from a mountain to a valley fold.
Folds in thin strips are extremely localized curvatures that occur over a small length-scale. In the limit of such length to zero, the shape of the strip transforms into one with a sharp corner (Lechenault et al., 2014). Such mechanical instabilities in thin strips are often the consequence of an over-constraining environment and material mismatch (Tanaka et al., 1987; Bowden et al., 1998), such as films on foundations.
Examples of patterns generated in confined conditions (Wang and Zhao, 2015) include wrinkles (Huang et al., 2004, 2005; Guvendiren et al., 2010), creases (Hong et al., 2009; Lechenault et al., 2014), blisters, folds (Conti and Maggi, 2007; Pocivavsek et al., 2008; Bayart et al., 2014), crumples (Vliegenthart and Gompper, 2006; Deboeuf et al., 2013), and ridges (Lobkovsky et al., 1995; Lobkovsky, 1996; Lobkovsky and Witten, 1997; Audoly and Pomeau, 2010). The transitions from one to the other manifest as a consequence of the change in the boundary conditions (Holmes and Crosby, 2010; Jin et al., 2015).
Several researchers explored the mechanical snap-back and snap-through instabilities of thin structures constrained by controlled boundary conditions, both on displacements and rotations.
In Beharic et al. (2014), the authors examined the bi-stable snap-through behavior of a beam with both positive and negative curvatures. To achieve such configuration, they confined the beam with applied symmetrical rotations at both ends and end-to-end horizontal span smaller than the arc-length. By deriving a strain energy landscape, they were able to determine a critical angle at which the deformation becomes mono-stable. Such critical angle depends monotonically only on the (compression) ratio between the arc-length and the horizontal span and not on the material or thickness, (Gomez et al., 2017) by studying a very similar system, demonstrated that the speed of the snapping is lower than the prediction based on the speed of sound; in fact, even without dissipation, they discovered the existence of a slow-down near criticality.
Cazzolli and Dal Corso (2018) studied the snapping mechanisms for elastic strips without folds, generalizing the work of Beharic et al. (2014). They derived an energy release map of all the possible rotations (symmetrical and non-symmetrical) for which the strip exhibits a snap mechanism. Other recent works on snapping in over-confined thin strips include (Yu and Hanna, 2019), where they used the Kirchhoff equations for anisotropic rods to explain the stability under compression, shear, and symmetric clamping: they revealed a series of boundary conditions that can generate stable configurations. Sano and Wada (2018) instead explored a range of controlled asymmetrical boundary conditions (hinged-clamped) that generate snapping.
From the literature survey, there appears to be a considerable amount of research in confined thin sheets, especially in determining their stability. Very little yet exists on confined folded sheet. By fold, in this paper, we mean a sharp corner in the deformation (Lobkovsky, 1996; Lechenault et al., 2014).
Despite the abundance of solutions of the Euler's Elastica (Bigoni, 2012; Manning, 2014), there seems to be a lack of closed-form (or even numerical) solutions of the Elastica containing folds as we intend in this paper. Related works are (Dado et al., 2004) and (Phungpaingam and Chucheepsakul, 2018), where they studied the post-buckling response of a cantilevered column (fixed-free ends) with a rotational spring in the middle, subjected to a concentrated horizontal force at the free end. However, there is no specific connection to folds, nor the spring contains a residual dihedral angle.
A sharp corner translates into a discontinuity of the rotation of the strip. Using a numerical meshfree method initially developed for cracks, our recently published work (Barbieri et al., 2019) presented a discretization method to model thick plates with multiple folds with infinite rotational stiffness.
In this paper, we are particularly interested in how constraining a thin strip with a fold with finite rotational stiffness leads to non-zero curvatures. Also, we study how such curvatures can be tuned by modifying the applied rotations at the boundaries (Figure 1). The question of whether such shapes are stable or unstable remains a topic that requires further investigation.
Figure 1. Curvature induced by hyperstatic confinement, and change in curvature (tuning) through applied rotations. (A) Folds and curvatures induced by applied rotations θ0 and θ1. (B) Signs of the curvatures.
2. Numerical Model Using Closed-Form Solutions of the Elastica
To examine which boundary conditions create curvatures in folded strips, we carried out semi-analytical simulations. Assuming invariance in the width, we modeled the strip as a plate under cylindrical bending. Therefore, instead of the Föppl-von Kármán equations, we modeled the strip as a planar, linearly elastic, unshearable, and inextensible rod, according to the Euler's theory of the Elastica.
with boundary conditions
where (·)′ = d/dS, with 0 ≤ S ≤ L being the curvilinear coordinate, μ is the curvature, θ is the rotation, x and y the Cartesian coordinates of the deformation, L the length of the rod, B = EI/(1 − ν2) the bending stiffness, E the Young modulus, I the second moment of area of the cross-section, ν the Poisson ratio, H0 the applied horizontal force at S = 0, V0 the applied vertical force at S = 0, M0 the bending moment at S = 0, and θ0 the rotation at S = 0, with positions x0, y0 at S = 0 (Figure 2A).
Figure 2. The Elastica with a fold at S = Sf. (A) The Elastica with a fold and its boundary conditions. (B) Dihedral angle β0, rest opening angle φ0 and discontinuity 〚θ〛 in the rotation.
We pass to the dimensionless form using the following normalizations:
and the following system of non-linear ordinary differential equations [the is removed for ease of notation]
with boundary conditions
where , with . Following standard techniques (Bigoni, 2012) of solution for Equation (5), we get
with θF0 such that
where am is the Jacobi amplitude function and
with F is the elliptic integral of the first kind and
The Cartesian coordinates of the deformation are
where dn is the Jacobi dn function
where E is the elliptic integral of the second kind, Imax = I(ϕmax) and is the number of times from 0 to S when ϕ′ changes sign. Furthermore,
The strain energy of the rod is
Let us now assume that a fold exists at 0 < Sf < 1. The fold is a discontinuity 〚·〛 on the rotation, with the following constitutive model (Lechenault et al., 2014)
with the curvilinear coordinate immediately to the right of Sf and the one immediately on the left (Figure 2B). The quantity κ > 0 is the rotational stiffness of the fold, and M(Sf) is the bending moment, assumed continuous in the absence of concentrated moments applied at Sf. The angle φ0 is the rest opening angle. When φ0 = 0, there is no fold; when φ0 < 0 the fold is called a mountain and when φ0 > 0 the fold is called a valley (Figure 1). The linear relationship between bending moment and dihedral angle was showed experimentally in Lechenault et al. (2014) and Jules et al. (2019) to hold for a wide range of dihedral angles. In their experiments, Lechenault and co-workers used the setup in Figure 3. The sample is clamped at one end, with a dynamometer attached to the other end. Both ends lie in the same plane. The right end can translate as the force pulls the specimen. The sample is illuminated by a laser, and the successive deformations are captured by a digital camera. In this way, it is possible to measure the dihedral angle β0. The moment at the fold is calculated as the product of the force measured by the dynamometer and the maximum height of the profile y(Sf).
Figure 3. Experimental setup used in Lechenault et al. (2014) for the characterization of the rotational stiffness of the fold (adapted).
The geometrical meaning of φ0 is shown in Figure 2B. The rest opening angle is related to the dihedral angle β0 by the simple relation
The rest opening angle is a mechanical property related to the yield stress of the material of the strip (Lechenault et al., 2014): to create a sharp fold, one needs to apply an extreme deformation that generates localized and irreversible strains. Such strains, in turn, lead to high stress concentrations beyond the yield stress of the material. An estimate of the dihedral angle is (Lechenault et al., 2014)
where ϵp is the plastic strain, σY the yield stress and E the Young modulus of the material composing the strip. In addition, φ0 appears to be independent from the thickness of the strip (Lechenault et al., 2014).
To tune the curvatures, we will further assume that at least one of the following quantities is assigned:
which renders one or all of H0, V0, M0 unknowns to be determined. The procedure to compute such unknowns is a shooting method (Press et al., 2007). Firstly, the solutions (7), (15), (16), and (22) are applied from 0 ≤ S ≤ Sf: let us call this solution μ−, θ−, x−, and y−. Then, the values θ(1), x(1), and y(1) are computed by applying the solution (7), (15), and (16) with the following boundary conditions:
the procedure is iterated with a non-linear solver until θ(1), x(1), and y(1) computed in this way match the conditions (27), within a certain tolerance.
In this section we report the results of the Elastica model for strips containing folds. We examine both isostatic and hyperstatic boundary conditions. In all the numerical calculations, we set the origin of the Cartesian system to the origin of the rod, therefore x0 = y0 = 0.
3.1. Isostatic Boundary Conditions: no Curvature
Let us examine the case where only y(1) is assigned
Figure 4 shows that no curvature is generated under conditions (30), regardless of the position of the fold and rotational stiffness. This configuration corresponds to the undeformed state of the folded strip when is at rest. This undeformed state depends only on the position of the fold Sf and the rest opening angle φ0.
Figure 4. Deformations for Sf = 1/5 and κ = 2, and for Sf = 1/2 and κ = ∞, both with rest opening angle φ = −60° with y1 = 0 (isostatic confinement). Blue continuous line: Elastica, red dashed line: Euler- Bernoulli (Equations 34).
3.2. Hyperstatic Boundary Conditions
In addition to condition (29), we consider cases where the boundary conditions are provided by a couple of rotations applied respectively at the start and the end of the strip:
We also impose that the end lies in the same plane as the start of the strip, therefore y1 = 0. Under boundary conditions (31),
Firstly, we consider the case of a fold with infinite rotational stiffness. By changing the applied rotations, we investigate the changes in sign of the curvatures at both sides of the fold (Figure 1B). Secondly, we repeat the same calculations for a finite rotational stiffness.
Thirdly, we isolate the effect of the rotational stiffness by varying κ for different strips under the strips under the same couple of applied rotations and the same position of the fold. We consider anti-symmetrical applied rotations.
3.2.1. Variable Applied Rotations and Rotational Stiffness
Figure 5 shows the deformations for a sheet under the same applied rotation, infinite rotational stiffness, same fold position and same rest angle.
Figure 5. Non-uniqueness of the deformation for a sheet with an infinitely stiff fold positioned at Sf = 0.5, with rest angle and applied rotations , .
There exist multiple solutions, depending on the initial moment M0 and shear force V0. This non-uniqueness is a consequence of the non-linearity of the Elastica: in fact, the solutions (15) and (16) are periodic in space, with the period depending on M0, H0 and V0 (Equation 13).
In the proceeding of the paper, we will always refer to the solution with the lowest strain energy. In fact, Equation (23) states that for the same θ0 and θ1, the strain energy grows with M0 and F0
Having resolved this disambiguation, we now examine the changes in sign of the curvatures (Figure 1B). Figure 6A shows the curvatures at both sides of the fold for varying applied rotations at both ends. We considered θ0 and θ1 ranging from −|φ0| to |φ0|, with an increment of 1°. The fold has position S0 = 1/2 and infinite rotational stiffness. With respect to the signs of the curvatures, there exist three regions: one where both curvatures are positive, one where both are negative and one where the curvatures have different signs.
Figure 6. Signs of the curvatures at the left M0 = μ(0) and at the right μ(1) of the fold. The rest angle is and the position is Sf = 1/2. (A) Infinite rotational stiffness (κ = ∞). (B) Finite rotational stiffness (κ = 1/2).
Figure 6B shows a map of the signs of the curvatures for the same strip with reduced rotational stiffness κ = 1/2. We can still identify three domains, as in Figure 6A, except that the domains with the same curvatures have decreased, while the sub-regions with opposite curvature have increased. Regardless of the rotational stiffness, the three regions meet at the same point, which corresponds to the isostatic boundary conditions.
Figures 7A,B show instead the dimensionless shear force at S = 0. Both graphs are identical, with the zero shear force line being straight with equation θ1 = −θ0.
Figure 7. Sign of the shear force V0 for a fold with rest angle φ = −60° and position Sf = 1/2. (A) Infinite rotational stiffness (κ = ∞). (B) Finite rotational stiffness (κ = 1/2).
3.2.2. Variable Rotational Stiffness for Symmetrical and Anti-symmetrical Applied Rotations
Figure 8 shows the deformations for a strip with a fold at Sf = 1/2, rest angle and variable rotational stiffness κ. The boundary conditions are given by a fixed y-coordinate at the end y1 = 0 and symmetrical applied rotations . The curvature μ(S) is uniform throughout the rod, and decreases with the decrease in κ. As expected, the dihedral angle β0 increases as κ decreases. When κ reaches a critical value κc, the discontinuity on the rotation disappears; below κc, the fold transitions from a mountain into a valley. In the limit of zero rotational stiffness, the fold becomes a perfect hinge: the deformation is equivalent to one of a flat rod, under isostatic boundary conditions, containing a fold with a rest opening angle equal to twice the applied rotation at the end.
Figure 8. Deformations for variable rotational stiffness for symmetrical applied rotations , with rest angle .
Instead, when the applied rotations are equal (anti-symmetrical rotations), the limit configuration for zero rotational stiffness is the solution of the Elastica with no folds (Figure 9).
Figure 9. Deformations for variable rotational stiffness for anti-symmetrical applied rotations , with rest angle .
In this section, we derive a simplified analytical model to explain the numerical results obtained in section 3.
4.1. Influence of the Isostatic Confinement
In this case, the solution of the system (5) is
with δ being the Dirac delta function, the Heaviside function and
is the ramp function. The solutions in Equation (34) correspond to a piece-wise roof-like deformation, where each side of the fold is flat.
The angle θ0 is given by
meaning that the deformation depends only upon the position of the fold Sf and the rest angle φ0, as anticipated in Figure 4.
4.2. Influence of the Applied Rotations
For simplicity, we will also assume x1 unassigned, making H0 = 0. In addition, we will assume, in first instance, that the fold is infinitely stiff, meaning κ → ∞. In this case, the Elastica coincides with the classic Euler-Bernoulli beam equation:
and the solutions for the beam are
Assigning the following boundary conditions
the two unknowns V0 and M0 are given by
For Sf = 1/2, the curvature will be uniform (V0 = 0) if
regardless of the rest angle of the fold. Such straight line in Equation (41) is observable in both Figures 7A,B.
Therefore, the curvature will be identically zero if both V0 = M0 = 0, which happens for
The curvature in Equation (38a) is a linear function in S. Therefore, by inspecting the signs of the values at S = 0 and S = 1, it is possible to discern whether the beam has the same curvature (either positive or negative) or opposite curvatures:
The values of μ(0) and μ(1) are
Figure 10. Curvature signs due to hyperstatic confinement through applied rotations. The fold has infinite rotational stiffness. Comparison between the Elastica model (dots) and Euler-Bernoulli beam (continuous straight lines). Pink dots: positive curvatures; green dots: negative curvatures; blue dots: opposite curvatures.
4.3. Influence of the Rotational Stiffness
Interestingly, the values in Equation (42) hold also for folds with finite rotational stiffness. In fact, when κ is finite, the solutions are
with given by Equation (40a) and by Equation (40b). By substituting Equations (42) into (46), the values of and also become zero.
In addition, for Sf = 1/2, θ0 = −θ1, also , which means that for opposite applied rotations, the curvature will be uniform also for folds with finite rotational stiffness.
The Euler-Bernoulli theory for soft folds in Equation (46) confirms the shrinkage of the domains with same-sign curvatures, as showed in Figure 11.
Figure 11. Curvature signs due to hyperstatic confinement through applied rotations. The fold has finite low rotational stiffness κ = 1/2. Comparison between the Elastica model (dots) and Euler-Bernoulli beam (continuous straight lines). Pink dots: positive curvatures; green dots: negative curvatures; blue dots: opposite curvatures.
The denominator in Equation (46) is always positive
preventing the appearance of any singularities.
The discontinuity in rotation at the fold is given by
and is showed in Figure 12 for and .
Figure 12. Discontinuity on the rotation at the fold for variable rotational stiffness. Comparison between values from Figure 8 and Equation (48).
The discontinuity becomes zero for
For Sf = 1/2:
When θ0 is of the same sign of φ0, then kc > 0 as in Figures 8, 13 and there exist a physical critical value of the rotational stiffness. Otherwise, for θ0 = −θ1 and φ0 of the opposite sign, then there is no critical value, and the solution is an undeformed roof-like deformation, as determined by Equation (42).
Figure 13. Discontinuity on the rotation at the fold for variable rotational stiffness. Comparison between values from Figure 9 and Equation (48).
It must be remarked, however, that even though a critical rotational stiffness for fold disappearance exists, the linear relationship in (24) might not hold for small rotational stiffness.
A thin sheet containing a fold assumes a roof-like undeformed state when confined isostatically. Instead, if confined hyperstatically, it will possess one or more curvatures. Such curvatures change according to the rotations applied at the ends of the sheet. Using the Euler's Elastica theory, we calculated the values of the applied rotations corresponding to a switch in the curvatures' signs.
We considered a fold as a discontinuity in the rotations, equal to the supplementary angle to its dihedral angle. The mechanical behavior of the fold depends upon two material properties: the rest opening angle and the rotational stiffness. The rest opening angle is a material property closely related to the yield stress; the rotational stiffness contributes to increasing the dihedral angle of the fold when subjected to a bending moment.
For infinite rotational stiffness, the dihedral angle depends only on the rest opening angle. For the signs of the curvatures before and after the fold, we can identify three regions: one where both curvatures are positive, one where they are both negative and one when they are opposite in sign. We presented a map of such occurrences according to the applied rotations; for small rotations and displacements, the numerical results agree with a simplified model based on the Euler-Bernoulli beam. Also, this map confirmed the presence of values of applied rotations for which the curvature is uniform throughout the sheet; in particular, there is one couple of applied rotations when such uniform curvature is zero. These rotations correspond to the isostatic confinement.
For folds with finite rotational stiffness, the two regions with same-sign curvatures shrink. Interestingly, there exists a critical value of the rotational stiffness when the fold disappears, meaning that the discontinuity on the rotation becomes zero. In particular, when the applied rotations are opposite in signs and the fold is a mountain, the critical value of the rotational stiffness is different from zero, and represents a transition from mountain to valley. In the limit of zero rotational stiffness, the sheet transforms into an undeformed valley with rest angle equal to twice the applied rotation at the end.
All datasets generated for this study are included in the manuscript and/or the supplementary files.
EB designed the research, conducted the analysis, and wrote the manuscript.
EB's work was supported by the Japan Society for the Promotion of Science JSPS KAKENHI Grant Number JP18K18065.
Conflict of Interest Statement
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Barbieri, E., Ventura, L., Grignoli, D., and Bilotti, E. (2019). A meshless method for the nonlinear von Kármán plate with multiple folds of complex shape - a bridge between cracks and folds. Comput. Mech. 1–19. doi: 10.1007/s00466-019-01671-w
Bowden, N., Brittain, S., Evans, A. G., Hutchinson, J. W., and Whitesides, G. M. (1998). Spontaneous formation of ordered structures in thin films of metals supported on an elastomeric polymer. Nature 393, 146. doi: 10.1038/30193
Dado, M., Al-Sadder, S., and Abuzeid, O. (2004). Post-buckling behavior of two elastica columns linked with a rotational spring. Int. J. Nonlinear Mech. 39, 1579–1587. doi: 10.1016/j.ijnonlinmec.2004.01.003
Guvendiren, M., Burdick, J. A., and Yang, S. (2010). Solvent induced transition from wrinkles to creases in thin film gels with depth-wise crosslinking gradients. Soft Matter 6:5795. doi: 10.1039/c0sm00317d
Manning, R. S. (2014). A catalogue of stable equilibria of planar extensible or inextensible elastic rods for all possible dirichlet boundary conditions. J. Elasticity 115, 105–130. doi: 10.1007/s10659-013-9449-y
Phungpaingam, B., and Chucheepsakul, S. (2018). Postbuckling behavior of variable-arc-length elastica connected with a rotational spring joint including the effect of configurational force. Meccanica 53, 2619–2636. doi: 10.1007/s11012-018-0847-x
Keywords: fold, crease, ridge, Elastica, discontinuous rotations
Citation: Barbieri E (2019) Curvature Tuning in Folded Strips Through Hyperstatic Applied Rotations. Front. Mater. 6:41. doi: 10.3389/fmats.2019.00041
Received: 18 December 2018; Accepted: 22 February 2019;
Published: 19 March 2019.
Edited by:Nicola Maria Pugno, University of Trento, Italy
Reviewed by:Lorenzo Bardella, Universitá degli Studi di Brescia, Italy
Stefano Vidoli, Sapienza University of Rome, Italy
Copyright © 2019 Barbieri. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ettore Barbieri, email@example.com