Skip to main content


Front. Mater., 19 March 2019
Sec. Mechanics of Materials
Volume 6 - 2019 |

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.

1. Introduction

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.

The equilibrium equations for the Elastica are (Landau and Lifshitz, 1959; Frisch-Fay, 1962)

μ+H0Bsinθ̄V0Bcosθ=0    (1a)
θ=μ    (1b)
x=cosθ    (1c)
y=sinθ    (1d)

with boundary conditions

μ(0)=M0B    (2a)
θ(0)=θ0    (2b)
x(0)=x0    (2c)
y(0)=y0    (2d)

where (·)′ = d/dS, with 0 ≤ SL 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:

S̄=SL  μ̄=μL  x̄=xL  ȳ=yL    (3)

resulting in

H0̄=H0L2B  V0̄=V0L2B  M0̄=M0LB      (4)

and the following system of non-linear ordinary differential equations [the (·-) is removed for ease of notation]

μ+H0sinθ̄V0cosθ=0    (5a)
θ=μ    (5b)
x=cos θ    (5c)
y=sinθ    (5d)

with boundary conditions

μ(0)=M0    (6a)
θ(0)=θ0    (6b)
x(0)=x0    (6c)
y(0)=y0    (6d)

where (·)=d/dS-, with 0S-1. Following standard techniques (Bigoni, 2012) of solution for Equation (5), we get

θ(S)=θF0+2ϕ(S)    (7)

with θF0 such that

H0=F0cosθF0    (8a)
V0=F0sinθF0    (8b)
F0=H02+V02>0    (8c)


ϕ(S)=sign(M0)am(α(S)|k2)    (9)

where am is the Jacobi amplitude function and

α(S)=sign(M0)McS+f0    (10)


f0=F(ϕ0|k2)    (11)

with F is the elliptic integral of the first kind and

ϕ0=θ0̄θF02    (12)
Mc2=M02+2F0(1̄cos2ϕ0)4    (13)
k2=F0Mc2    (14)

The Cartesian coordinates of the deformation are

x(S)=x0+cosθF0Cϕ(S)̄sinθF0Sϕ(S)    (15)
y(S)=y0+sinθF0Cϕ(S)+cosθF0Sϕ(S)    (16)


Sϕ(S)=0Ssin(2ϕ(S))dS=2Mck2(dn(α(S)|k2)̄dn(f0|k2))    (17)

where dn is the Jacobi dn function

Cϕ(S)=0Scos(2ϕ(S))dS=S̄2Mc(signϕI(ϕ)signϕ0I(ϕ0)+2Imax nϕ(S))    (18)


I(ϕ)=1k2(F(ϕ|k2)̄E(ϕ|k2))    (19)

where E is the elliptic integral of the second kind, Imax = I(ϕmax) and nϕ is the number of times from 0 to S when ϕ′ changes sign. Furthermore,

ϕmax=arcsin(1k)    (20)
ϕ(S)=sign(M0)Mcdn(α(S)|k2)    (21)


μ(S)=2ϕ(S)    (22)

The strain energy of the rod is

U=1201μ2dS=201ϕ(S)2dS=2Mc2̄2F001sin2ϕdS   =2Mc2̄F0(1̄Cϕ(1))    (23)

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)

θ=θ(Sf+)̄θ(Sf̄)=M(Sf)κ+φ0    (24)

with Sf+ the curvilinear coordinate immediately to the right of Sf and Sf- 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

β0=π̄(φ0+M(Sf)κ)    (25)

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)

φ0~ϵp~σYE    (26)

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:

θ(0)=θ0  θ(1)=θ1  x(1)=x1  y(1)=y1    (27)

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 ≤ SSf: 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:

θ(Sf)=θ̄(Sf)+(φ0+μ̄(Sf)κ)    (28a)
μ(Sf)=μ̄(Sf)    (28b)
x(Sf)=x̄(Sf)    (28c)
y(Sf)=ȳ(Sf)    (28d)

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.

3. Results

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

y(1)=y1    (29)


H0=0  M0=0  V0=0  θ0=unknown    (30)

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:

y(1)=y1=0  θ(0)=θ0  θ(1)=θ1    (31)

We also impose that the end lies in the same plane as the start of the strip, therefore y1 = 0. Under boundary conditions (31),

H0=0  M0=unknown  V0=unknown    (32)

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 φ0=-60° and applied rotations θ0=20°, θ1=-20°.

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

U~Mc2~M02+ F0    (33)

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 φ0=-60° 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 φ0=-45° and variable rotational stiffness κ. The boundary conditions are given by a fixed y-coordinate at the end y1 = 0 and symmetrical applied rotations θ0=-θ1=-20°. 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 θ0=-θ1=-20°, with rest angle φ0=-45°.

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 θ0=θ1=-20°, with rest angle φ0=-45°.

4. Discussion

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

μ(S)=φ0δ(S̄Sf)    (34a)
θ(S)=θ0+φ0H(S̄Sf)    (34b)
x(S)=cosθ0S̄(cosθ0̄cos(θ0+φ0))S̄Sf    (34c)
y(S)=sinθ0S̄(sinθ0̄sin(θ0+φ0))S̄Sf    (34d)

with δ being the Dirac delta function, H the Heaviside function and

S=S+|S|2    (35)

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

θ0=atan2(̄(1̄Sf)sinφ0Sf+cosφ0(1̄Sf))    (36)

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:

μ=V0    (37)

and the solutions for the beam are

μ(S)=M0+V0S+φ0δ(S̄Sf)    (38a)
θ(S)=θ0+M0S+V0S22+φ0H(S̄Sf)    (38b)
x(S)=S    (38c)
y(S)=θ0S+M0S22+V0S33+φ0S̄Sf    (38d)

Assigning the following boundary conditions

θ(0)=θ0  θ(1)=θ1  y(1)=y1=0    (39)

the two unknowns V0 and M0 are given by

V0=6(θ0+θ1+φ0(1̄2Sf))    (40a)
M0=̄4(θ0+θ12+φ0(1̄32Sf))    (40b)

For Sf = 1/2, the curvature will be uniform (V0 = 0) if

θ0=̄θ1    (41)

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

θ0=φ0(1̄Sf)  θ1=̄φ0Sf    (42)

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:

μ(0)>0   μ(1)>0  same positive curvature    (43a)
μ(0)<0   μ(1)<0  same negative curvature    (43b)
μ(0)μ(1)<0  opposite curvatures    (43c)

The values of μ(0) and μ(1) are

μ(0)=M0    (44a)
μ(1)=M0+V0    (44b)

Figure 10 shows the presence of the three regions for the signs of the curvatures, as observed in Figure 6; also, there is agreement between the Euler-Bernoulli theory and the Elastica model.


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

μ(S)=M0κ+V0κS+(φ0+M0κ+V0κSfκ)δ(S̄Sf)    (45a)
θ(S)=θ0+M0κS+V0κS22+(φ0+M0κ+V0κSfκ)H(S̄Sf)    (45b)
x(S)=S    (45c)
y(S)=θ0S+M0κS22+V0κS33+(φ0+M0κ+V0κSfκ)S̄Sf    (45d)


V0κ=V01+4κ(1̄3Sf+3Sf2)+12θ1+Sf(θ0̄θ1)1+4κ(1̄3Sf+3Sf2)    (46a)
M0κ=M01+4κ(1̄3Sf+3Sf2)+12Sfθ1+Sf(θ0̄θ1)1+4κ(1̄3Sf+3Sf2)    (46b)

with V0 given by Equation (40a) and M0 by Equation (40b). By substituting Equations (42) into (46), the values of V0κ and M0κ also become zero.

In addition, for Sf = 1/2, θ0 = −θ1, also V0κ=0, 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

κ+12Sf2̄12Sf+4>0    0<Sf<1 κ>0    (47)

preventing the appearance of any singularities.

The discontinuity in rotation at the fold is given by

θκ=θ(Sf+)̄θ(Sf̄)=M(Sf)κ+φ0             =φ0κ+θ0(6Sf̄4)+θ1(6Sf̄2)κ+12Sf2̄12Sf+4    (48)

and is showed in Figure 12 for θ0=-θ1=20° and φ0=-45°.


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

κc=4θ0φ0(1̄32Sf)+2θ1φ0(1̄3Sf)    (49)

For Sf = 1/2:

κc=θ0φ0̄θ1φ0    (50)

therefore, when the applied rotations are equal, then the critical rotational stiffness is zero, as in Figures 9, 12, when the applied rotations are opposite, then

κc=2θ0φ0    (51)

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.

5. Conclusions

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.

Data Availability

All datasets generated for this study are included in the manuscript and/or the supplementary files.

Author Contributions

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.


Audoly, B., and Pomeau, Y. (2010). Elasticity and Geometry: From Hair Curls to the Non-linear Response of Shells. Oxford: Oxford University Press.

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

CrossRef Full Text | Google Scholar

Bayart, E., Boudaoud, A., and Adda-Bedia, M. (2014). Tuning the ordered states of folded rods by isotropic confinement. Phys. Rev. E 89:012407. doi: 10.1103/PhysRevE.89.012407

PubMed Abstract | CrossRef Full Text | Google Scholar

Beharic, J., Lucas, T. M., and Harnett, C. K. (2014). Analysis of a compressed bistable buckled beam on a flexible support. J. Appl. Mech. 81:081011. doi: 10.1115/1.4027463

CrossRef Full Text | Google Scholar

Bigoni, D. (2012). Nonlinear Solid Mechanics: Bifurcation Theory and Material Instability. Cambridge: Cambridge University Press.

Google Scholar

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

CrossRef Full Text | Google Scholar

Cazzolli, A., and Dal Corso, F. (2018). Snapping of elastic strips with controlled ends. Int. J. Solids Struct. 162, 285–303. doi: 10.1016/j.ijsolstr.2018.12.005

CrossRef Full Text | Google Scholar

Conti, S., and Maggi, F. (2007). Confining thin elastic sheets and folding paper. Arch. Ration. Mech. Anal. 187, 1–48. doi: 10.1007/s00205-007-0076-2

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Deboeuf, S., Katzav, E., Boudaoud, A., Bonn, D., and Adda-Bedia, M. (2013). Comparative study of crumpling and folding of thin sheets. Phys. Rev. Lett. 110:104301. doi: 10.1103/PhysRevLett.110.104301

PubMed Abstract | CrossRef Full Text | Google Scholar

Frisch-Fay, R. (1962). Flexible Bars. London: Butterworths.

Gomez, M., Moulton, D. E., and Vella, D. (2017). Critical slowing down in purely elastic snap-through instabilities. Nat. Phys. 13, 142–145. doi: 10.1038/nphys3915

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Holmes, D., and Crosby, A. J. (2010). Draping films: a wrinkle to fold transition. Phys. Rev. Lett. 105:038303. doi: 10.1103/PhysRevLett.105.038303

PubMed Abstract | CrossRef Full Text | Google Scholar

Hong, W., Zhao, X., and Suo, Z. (2009). Formation of creases on the surfaces of elastomers and gels. Appl. Phys. Lett. 95:111901. doi: 10.1063/1.3211917

CrossRef Full Text | Google Scholar

Huang, Z., Hong, W., and Suo, Z. (2004). Evolution of wrinkles in hard films on soft substrates. Phys. Rev. E 70:030601. doi: 10.1103/PhysRevE.70.030601

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, Z., Hong, W., and Suo, Z. (2005). Nonlinear analyses of wrinkles in a film bonded to a compliant substrate. J. Mech. Phys. Solids 53, 2101–2118. doi: 10.1016/j.jmps.2005.03.007

CrossRef Full Text | Google Scholar

Jin, L., Takei, A., and Hutchinson, J. W. (2015). Mechanics of wrinkle/ridge transitions in thin film/substrate systems. J. Mech. Phys. Solids 81, 22–40. doi: 10.1016/j.jmps.2015.04.016

CrossRef Full Text | Google Scholar

Jules, T., Lechenault, F., and Adda-Bedia, M. (2019). Local mechanical description of an elastic fold. Soft Matter. 15, 1619–1626. doi: 10.1039/C8SM01791C

PubMed Abstract | CrossRef Full Text | Google Scholar

Landau, L. D., and Lifshitz, E. M. (1959). Course of Theoretical Physics Vol 7: Theory and Elasticity. London: Pergamon Press.

Lechenault, F., Thiria, B., and Adda-Bedia, M. (2014). Mechanical response of a creased sheet. Phys. Rev. Lett. 112:244301. doi: 10.1103/PhysRevLett.112.244301

PubMed Abstract | CrossRef Full Text | Google Scholar

Lobkovsky, A., Gentges, S., Li, H., Morse, D., and Witten, T. A. (1995). Scaling properties of stretching ridges in a crumpled elastic sheet. Science 270, 1482–1485.

Google Scholar

Lobkovsky, A. E. (1996). Boundary layer analysis of the ridge singularity in a thin plate. Phys. Rev. E 53:3750. doi: 10.1103/PhysRevE.53.3750

PubMed Abstract | CrossRef Full Text | Google Scholar

Lobkovsky, A. E., and Witten, T. (1997). Properties of ridges in elastic membranes. Phys. Rev. E 55:1577. doi: 10.1103/PhysRevE.55.1577

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Pocivavsek, L., Dellsy, R., Kern, A., Johnson, S., Lin, B., Lee, K. Y., et al. (2008). Stress and fold localization in thin elastic membranes. Science 320, 912–916. doi: 10.1126/science.1154069

PubMed Abstract | CrossRef Full Text | Google Scholar

Press, W. H., Teukolsky, S. A., Vetterling, W. T., and Flannery, B. P. (2007). Numerical Recipes 3rd Edition: The Art of Scientific Computing. Cambridge: Cambridge university press.

Google Scholar

Sano, T. G., and Wada, H. (2018). Snap-buckling in asymmetrically constrained elastic strips. Phys. Rev. E 97:013002. doi: 10.1103/PhysRevE.97.013002

PubMed Abstract | CrossRef Full Text | Google Scholar

Tanaka, T., Sun, S.-T., Hirokawa, Y., Katayama, S., Kucera, J., Hirose, Y., et al. (1987). Mechanical instability of gels at the phase transition. Nature 325:796. doi: 10.1038/325796a0

CrossRef Full Text | Google Scholar

Vliegenthart, G. A., and Gompper, G. (2006). Forced crumpling of self-avoiding elastic sheets. Nat. Mater. 5:216. doi: 10.1038/nmat1581

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Q., and Zhao, X. (2015). A three-dimensional phase diagram of growth-induced surface instabilities. Sci. Rep. 5:8887. doi: 10.1038/srep08887

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, T., and Hanna, J. (2019). Bifurcations of buckled, clamped anisotropic rods and thin bands under lateral end translations. J. Mech. Phys. Solids 122, 657–685. doi: 10.1016/j.jmps.2018.01.015

CrossRef Full Text | Google Scholar

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,