Closed Unstretchable Knotless Ribbons and the Wunderlich Functional

In 1962, Wunderlich published the article “On a developable Möbius band,” in which he attempted to determine the equilibrium shape of a free standing Möbius band. In line with Sadowsky’s pioneering works on Möbius bands of infinitesimal width, Wunderlich used an energy minimization principle, which asserts that the equilibrium shape of the Möbius band has the lowest bending energy among all possible shapes of the band. By using the developability of the band, Wunderlich reduced the bending energy from a surface integral to a line integral without assuming that the width of the band is small. Although Wunderlich did not completely succeed in determining the equilibrium shape of the Möbius band, his dimensionally reduced energy integral is arguably one of the most important developments in the field. In this work, we provide a rigorous justification of the validity of the Wunderlich integral and fully formulate the energy minimization problem associated with finding the equilibrium shapes of closed bands, including both orientable and nonorientable bands with arbitrary number of twists. This includes characterizing the function space of the energy functional, dealing with the isometry and local injectivity constraints, and deriving the Euler–Lagrange equations. Special attention is given to connecting edge conditions, regularity properties of the deformed bands, determination of the parameter space needed to ensure that the deformation is surjective, reduction in isometry constraints, and deriving matching conditions and jump conditions associated with the Euler–Lagrange equations.


Introduction
For the past half-century, efforts to find the equilibrium shapes of Möbius bands and various other ribbon-like objects made from materials like paper, which can be modeled as unstretchable two-dimensional sheets, have relied mainly on numerical strategies for minimizing the Wunderlich (1962) functional, 1 the derivation of which can be viewed as beginning with the bending energy 2μ S H 2 da (1) for a developable surface S with bending modulus μ > 0, mean curvature H , and element of area da. Granted that S is a half-twist Möbius band with midline M of length and uniform width 2b, Wunderlich (1962) noted that S must lie on the rectifying developable of M and on that basis concluded that the surface integral (1) reduces without approximation to a line integral whereκ 0 is the restriction to M of the nonvanishing principal curvature of S, V is the point of regression of the generatrix of S, and ds is the element of arclength along M.
In (2), it is tacit that the curvature of S is at least piecewise continuous. Wunderlich (1962) worked with dimensionless integrals that differ from (1) and (2) by a factor of 2/μ. We nevertheless refer to (2) as the Wunderlich functional. The purpose of Wunderlich's (1962) paper was to find the equilibrium shape of a developable Möbius band made from a rectangular strip of paper. Following Sadowsky (1930) he promoted the idea of finding the equilibrium shape by using an energy minimization principle 2 . His derivation of (2) provided an important step in that direction. However, as he concluded, finding the final equilibrium shape "would have required considerable effort, which is why it must be left undone for the time being." The following comments provide a comprehensive understanding of this task and the effort needed to accomplish it.
• For Wunderlich (1962), the equilibrium shape under discussion is a solution to a variational problem equivalent to 2μ S H 2 da = min.
However, Wunderluch did not mention that a proper formulation of the problem (3) should begin with first specifying a reference configuration D of the strip of paper 1 See Todres (2015) for a translation of Wunderlich's (1962) paper into English.
2 See Hinz and Fried (2015) for a translation of Sadowsky's (1930) paper into English. and then finding the entire family F of kinematically admissible deformations χ of D. Since the mean curvature H depends on χ implicitly through the current configuration S = χ (D), the bending energy E is a functional of the deformation χ : The variational problem (3) therefore amounts to finding a deformation χ e that complies with the condition Specifying the family F of kinematically admissible deformations χ in (4) is essential in formulating the variational problem. • Kinematic admissibility here means, among other things, that the deformed surface must be isometric to the flat strip, a requirement that stems from the material being unstretchable. This means that for a deformation χ to be kinematically admissible, it must be an isometry, as defined, for instance, by do Carmo (1976, §4-2). Moreover, to disallow self-penetration of the surface, all deformations must be injective. • To ensure that the strip deforms to a closed ribbon, certain matching edge conditions must be imposed on the class of admissible deformations. As Wunderlich (1962) put it, a Möbius band is formed "by connecting the narrow edges of a long, rectangular strip of paper so that all diagonally opposed vertices meet." However, Wunderlich (1962) did not explicitly invoke any such conditions when deriving (2), and moreover, it has been shown by Chen et al. (2018) that for nonclosed ribbons, (2) does not accurately represent the bending energy (1). • Wunderlich (1962) used a "partially numerical, partially graphical procedure" to find an "approximate portrayal of the equilibrium model" but did not derive the Euler-Lagrange equations of an appropriately constrained version of the functional (2) whose solution determines the exact equilibrium shape of the Möbius band.
The study of unstretchable material surfaces has been taken up by others in the aftermath of the pioneering work of Wunderlich (1962). Starostin andvan der Heijden (2007, 2015) minimized the bending energy (1) over the class of Möbius bands of a fixed length and width that lie on a rectifying developable surface and dimensionally reduced the bending energy to obtain the Wunderlich energy. They then used a higher-order variational approach introduced by Gay- Balmaz et al. (2012) to obtain the corresponding Euler-Lagrange equations and specialized those equations to Möbius bands with certain symmetries. To construct a rectifying developable surface from an energy minimizing midline, they implicitly ruled out situations where the curvature of the midline vanishes over any open interval. However, Sadowsky's (1930) construction of a Möbius band shows that there exist bands for which that assumption is violated. Hence, there exist Möbius bands made from unstretchable strips that do not lie on rectifying developable surfaces. The problem of minimizing the bending energy (4) over the collection of rectifying developable Möbius bands is therefore not the same as that of minimizing over all Möbius bands of a given length and width. While it might well be that a Möbius band with minimum bending energy lies on a rectifying developable, it is illegitimate to assume so a priori. The approach of Starostin and van der Heijden also yields only the shapes of Möbius bands with various symmetries without determining the underlying isometric deformations from a reference configuration. Moreover, the matter of injectivity is not considered.
There are also a series of papers where Hornung (2011a, b, c) studied the structure of isometric deformations that are locally injective and are minimizers of the bending energy (4). In particular, these papers show that locally the surface corresponding to the deformation with minimum energy has a parameterization satisfying certain conditions which can be interpreted as Euler-Lagrange equations for the bending energy. Although these equations provide a basis for establishing regularity results for extrema, they do not constitute a boundary-value problem for determining the deformation from a reference configuration. The relevance of Hornung's (2011a;2011b;2011c) analysis to the equilibrium shapes of Möbius strips was discussed by Bartels and Hornung (2015). More broadly, the study of isometric immersions has been a subject of recent interest in the mathematics community, as exemplified by the works of Pakzad (2004), Müller and Pakzad (2005), and Kirchheim (2003). Dias and Audoly (2015) studied the closely related problem of finding the equilibrium shape of a thin, unstretchable, elastic ribbon that is obtained by isometrically deforming a flat reference ribbon which need not be rectangular, but which must have a unique, well-defined midline. They utilized a double parameterization to describe the deformation: one parameterization for the reference configuration and another for the deformed configuration. The parameterization for the deformed configuration allows for a dimensional reduction in the bending energy, leading to an energy similar in spirit to Wunderlich's but for a broader class of ribbons. The parameterizations utilized in this work do not necessarily completely cover the reference and current configurations. This results in a dimensionally reduced energy that may underestimate or overestimate the true bending energy of the ribbon, as demonstrated by Chen et al. (2018).
The present work is aimed at filling the gaps left by Wunderlich (1962), gaps that, as itemized above, have not yet been treated in the ensuing literature. Our work contains many novel contributions, the most salient of which we now list: 1. We show that a kinematically admissible deformation χ can be represented by two vector-valued functions of arclength: a directrix d and a generatrix g, which satisfy certain constraints. By this, we mean that each kinematically admissible deformation induces a directrix and generatrix that satisfy the constraints, and conversely, that given a directrix and generatrix satisfying these constraints, it is possible to construct a surface S and a kinematically admissible deformation from D to S. 2. We formulate explicit matching conditions that ensure that the short edges of the strip are smoothly joined to form a closed ribbon and consistently apply those conditions when deriving the dimensionally reduced bending energy. 3. We show that the bending energy (1) can be dimensionally reduced to a functional of the directrix d and the generatrix g and we establish the equivalence of that functional and the Wunderlich functional (2). 4. We introduce Lagrange and Karush-Kuhn-Tucker multipliers to enforce the constraints on the the directrix d and generatrix g, and we formulate an associated variational problem involving an augmented energy functional. 5. We use elementary methods to derive the Euler-Lagrange equations associated with our augmented energy functional, resulting in a system of ordinary differential equations amenable to numerical analysis.
We begin by identifying the family F of kinematically admissible deformations χ, each of which is a mapping from a reference configuration D, which we take to be a flat rectangle of length and width 2b, to a spatial configuration that must satisfy certain matching edge conditions. Each deformation χ is required to be an isometry, meaning that it preserves the length of material curves. We refer to such deformations as isometric deformations. Each such χ can be represented by a pair of ruled parameterizations, one for the reference configuration and the other for the spatial configuration, namelŷ where α and β are the parameters for the referential and spatial configurations, c and d are, respectively, referential and spatial directrices, and f and g are, respectively, referential and spatial generatrices. The stipulation that χ be an isometric deformation amounts to a collection of conditions involving the referential and spatial directrices and generatrices, along with the derivatives of those quantities with respect to arclength. We find, however, that those conditions can be recast entirely in terms of an equivalent collection of purely spatial conditions. These conditions are where a prime denotes differentiation with respect to the parameter α. Whereas (7) 1 and (7) 2 simply require that the spatial directrix be parametrized by arclength and that the spatial generatrix be of unit magnitude, (7) 3 requires that the normal vector of the spatial directrix be orthogonal to the spatial generatrix and (7) 4 requires that S be developable. Hence, contrary to some statements in the literature, dictating that M be an inextensible space curve is not sufficient to ensure that a minimizer of (2) corresponds to an isometric deformation. Restricting attention to situations where the deformation χ is such that the rulings which form the ribbon in any particular configuration do not intersect, we show that this requirement implies that where θ is the angle between d and g. Moreover, we show that the condition (8) ensures that χ is locally injective on the interior of D. While there might be reasons to relinquish (8) in favor of the stronger requirement that χ be globally injective on the interior of D, this restriction is difficult to implement so we settle here with imposing local injectivity.
Although Wunderlich (1962) considered only half-twist Möbius bands, the steps leading from (1) to (2) are independent of the number of twists that S possesses. Granted that (2) can be used to calculate the bending energy (1) of a half-twist Möbius band S, it should apply just as well if S is instead a closed orientable ribbon, with or without twist, or a Möbius band with more than one half-twist. Our choice of reference configuration accommodates all such possibilities. We find, however, that the explicit form of the referential directrix depends on whether S is orientable or nonorientable. In either case, the referential directrix includes not only the midline of the rectangle D but also segments of its short edges. If S is orientable, then those segments must lie on opposite sides of M. If S is nonorientable, then those segments must lie on the same side of M. These additional segments are needed to ensure that D and S are completely covered by rulings, and thus, as is crucial, that the deformation χ preserves the identity of material points.
The distinction between orientable and nonorientable closed ribbons is imposed through two alternative conditions that the deformation χ must satisfy on the short edges of D. For either of these alternatives, χ maps the short edges of D to a single curve on S and thus can at best be injective on the interior of D. The curve in question may be a ruling of S, but this is not generally so. If that curve is a ruling of S, then the reduction from (1) to (2) is trivial. We therefore focus on situations in which the short edges of D are not associated with a ruling.
Our approach affords a unified derivation of a dimensionally reduced version of the bending energy (1) regardless of the orientability of the closed ribbon S. The resulting functional involves integration over the midline M of S and depends on both the parametrizations d and g of the spatial directrix and spatial generatrix. In contrast to being independent, the variations of d and g must be consistent with the spatial constraints (7) to ensure that χ is an isometric deformation. Using the relations d · g = cos θ and V = sin 2 θ θ , we prove with a change of variables that our dimensional reduction is equivalent to the Wunderlich functional (2). We therefore provide an independent basis for that functional and, in so doing, confirm that it can, as suggested above, be used to calculate the bending energy (1) of a closed orientable ribbon, with or without twist, or a Möbius band with any number of half-twists. In that regard, the number of twists that can be sustained without violating the requirement that the curvature of S be piecewise continuous or necessitating points of self-contact is not unlimited but depends on the half-width to length aspect ratio b/ . To deduce the dimensionally reduced bending energy hinges on appropriate extensions of the parametrization of S, the precise features of which depend on the orientability of S. In both cases, however, the extensions make use of the requirements involving rulings that terminate on the short ends of the reference region D.
Having shown that (1) can be expressed as a functional of the functions d and g that parametrize the spatial directrix and spatial generatrix, we prove that knowledge of those quantities is sufficient to construct both a C 1 closed ribbon S with piecewise continuous curvature and an underlying isometric deformation χ from the rectangular strip D to S. This construction relies on assuming that d and g satisfy the constraints (7) and the local injectivity condition (8).
To assure that constraints (7) and (8) are met, we introduce four Lagrange multipliers together with a Karush-Kuhn-Tucker multiplier and consider a correspondingly augmented version of our dimensionally reduced energy functional. We then derive the Euler-Lagrange equations arising from that functional. In addition to ordinary differential equations for d and g, the Euler-Lagrange equations include matching conditions that apply at the point of the midline where the ends of D are joined and jump conditions that apply at points on M across which the curvature of S is discontinuous. The Euler-Lagrange equations must be augmented by the isometry constraints (7) and complementary slackness condition associated with the Karush-Kuhn-Tucker multiplier in conjunction with an additional condition which encodes the number of twists in S through a constraint on the self-linking number of its midline M.
The remainder of the paper is organized as follows. In Sect. 2, we focus on kinematical issues. After introducing the notion of an isometric deformation χ from a flat rectangular strip D to a closed ribbon S, we present the connecting edge conditions that decide whether or not S is orientable, explain the basis for treating D and S as ruled surfaces, consider the degree of regularity that is allowed if two rulings meet at a point on the boundary of S, and define the explicit parametrizations of the referential directrix that we utilize in the orientable and nonorientable cases. Having presented ruled parametrizations for D and S, we conclude this section by describing the extensions of the parametrization of S that are needed subsequently to achieve our dimensional reduction. In Sect. 3, we define the bending energy, conduct the dimensional reduction, and exhibit the calculations needed to convert our dimensional reduction to the Wunderlich functional. In Sect. 4, we consider questions related to the regularity of S and χ, derive versions of the constraint of isometry in terms of the functions d and g that parametrize the spatial directrix and spatial generatrix, and show how S and χ can be constructed given knowledge of d and g. In Sect. 5, we derive the Euler-Lagrange equations for the energy functional arising from augmenting our dimensionally reduced bending energy to account for the constraints that d and g must satisfy to ensure that χ is locally injective and an isometry. We also provide interpretations of the matching conditions that apply at the point on the midline M of S where the end points of the midline of D are joined consistent with the connecting edge conditions. Finally, we summarize and discuss our most salient results in Sect. 6.

Isometric Deformations
Let D be a rectangular strip, boundary included, of length and width 2b that is embedded in a three-dimensional Euclidean point space E. Suppose that D is made from an unstretchable material, meaning that it can sustain only isometric deformations. Let χ be such a deformation, assumed to be C 2 on the interior of D, and denote the surface in E corresponding to the spatial configuration of D by S = χ(D).

Connecting Edge Conditions
We consider only isometric deformations χ that are C 2 on the interior of D and satisfy certain connecting edge conditions. To describe those conditions, we choose an origin o, a positively oriented orthonormal basis {ı 1 , ı 2 , ı 3 }, and associated Cartesian coordinates We then assume that S is a C 1 surface and that one of the connecting edge conditions or holds. For either (10) or (11), the surface S is what we call a "closed ribbon." Whereas S is orientable for (10), it is nonorientable for (11). We allow for the possibility of twisting D any number of times in conjunction with the imposition of (10) or (11). If S is orientable, then its boundary ∂S has two disjoint edges. If, otherwise, S is nonorientable, then ∂S has only a single edge. We restrict attention to situations where ∂S is free of both traction and couple traction regardless of the orientability of S.

Closed Ribbons as Ruled Surfaces
Since the Gaussian curvature of S must vanish for any isometric deformation χ of D, any closed ribbon S can be partitioned depending on whether or not its mean curvature H vanishes. The work of Hartman and Nirenberg (1959) shows that the subset of S within which H = 0 can be covered by straight line segments, or rulings, which are asymptotic curves. Due to the rectangular shape of D, the connected components of S where H = 0 must be a trapezoid, a triangle, or an isolated line segment. With such a partition, S can be represented as a ruled surface. Specifically, any curved part of S is ruled by the corresponding asymptotic curves. Furthermore, any trapezoidal or triangular subset of S can be ruled in a way that ensures a continuous transition to the adjacent curved portions of S and an isolated line segment of S upon which H = 0 automatically provides a smooth transition between the adjacent curved portions of S. A complete ruling of S induces an associated ruling of D via the underlying isometric deformation χ. Thus, since the connecting edge conditions (10) and (11) make it impossible for any particular ruling of D to span its entire length, each ruling of S must intersect the boundary ∂S of S at exactly two points.
Let the left-and right-hand edges, , of D be denoted by ∂D L and ∂D R , respectively, and choose a fixed value of x 2 belonging to [0, 2b]. Then, if (10) holds, so that S is orientable, the slopes of the rulings that terminate at the points x 2 ı 2 and ı 1 + x 2 ı 2 on ∂D L and ∂D R must be identical, as depicted in Fig. 1a. Moreover, the deformed images of those rulings must meet on S to form a single ruling that contains the point χ (0, x 2 , 0) = χ ( , x 2 , 0). Alternatively,  (10) and (11). The red curves indicate the directrices. Whereas the slopes of the yellow rulings containing the points c(α) and c(α(α)) are equal for orientable bands, they are of equal magnitude but opposite sign for nonorientable bands. Hence, although the rulings passing through c(α) and c(α(α)) are disjoint on D, their images under the deformation χ join smoothly to form a single ruling on the band S = χ(D). Restricting attention to deformations for which the rulings intersecting the right-hand edge of D have positive slopes if (11) holds, so that S is nonorientable, the slopes of the ruling that terminate at the points x 2 ı 2 and ı 1 + (2b − x 2 )ı 2 on ∂D L and ∂D R must differ only by their signs, as depicted in Fig. 1b. Moreover, the deformed images of those rulings must meet to form a single ruling on S that contains the point Let J denote the line segment on S along which ∂D L and ∂D R are joined via either of the connecting edge conditions (10) or (11). Then, J might or might not be a ruling of S and these alternatives are mutually exclusive. Since the situation in which J is a ruling of S is both atypical and straightforward to treat, we consider only the more prevalent and substantially more challenging alternative in which J is not a ruling of S. If J is not a ruling of S, then there is necessarily a family of rulings that cross J . That being so, no generality is lost by considering deformations for which the rulings that intersect the right-hand edge of D have positive slopes as depicted in Fig. 1a, b.

Kinematically Imposed Smoothness of the Deformation on the Boundary of the Strip
Due to the restrictive nature of the class of isometric deformations, it can be shown that, in addition to being C 2 on the interior of the rectangular strip D, the deformation χ from D to a band S is also C 2 at all points on the boundary ∂D of D at which rulings do not intersect. To verify this assertion, consider a point x • on ∂D at which rulings do not intersect. Since the direction of the rulings is a continuous function of position, rulings that are sufficiently close to x • can be extended slightly outside of D without intersecting. Thus, there is a small two-dimensional neighborhood N of x • , in the x 3 = 0 plane, that is covered by rulings that do not intersect. Since an isometric deformation χ transforms rulings rigidly, 3 χ has a natural isometric extension to N . Moreover, by the construction of that extension, χ has the same regularity on N as it does on the interior of D. Thus, the isometric extension of χ is C 2 on N , and, consequently, the original deformation χ is itself C 2 at x • . To illustrate that χ need not be smooth at any point of the boundary ∂D of the strip D at which rulings intersect, let T c denote a curved portion of S covered by rulings that intersect at a point on ∂S, and let T r be the part of the rectangular strip D associated with T c . Notice that T r is a triangular region, as depicted in Fig. 2. Thus, T r can be parameterized by polar coordinates (φ, ρ), with the point where the rulings intersect serving as the origin, with φ denoting the angle measured counterclockwise from the x 1 -axis, and ρ measuring the radius from that origin. If p denotes the position on S where the rulings intersect, then the curved portion T c of S under consideration can be parameterized byr where g(φ) is a unit vector in the direction of the ruling associated with the angle φ.
A unit normal n to T c is then given by As defined, n is independent of position along each ruling and, thus, is independent of β. Since T c is a curved portion of S, n cannot be uniform on T c . We thus see that T c has rulings along which n takes different values. Let n 1 and n 2 denote the normal vectors on two such rulings. If χ were C 1 , it would then follow that the unit normal to S would be continuous up to ∂S, implying that n 1 = n 2 , which is incompatible with the curved nature of T c . We thus conclude, as claimed, that χ need not be C 1 at a point on ∂D where the rulings intersect.

Parameterizations of the Strip and the Closed Ribbon
In describing the rectangular strip D and the ribbon S as parametrized surfaces, it is necessary to consider whether the closed ribbon S is orientable or nonorientable. In either case, the parametrizations rest on a convenient identification of a curve C r that Fig. 3 Examples of (a) orientable and (b) nonorientable closed ribbons S. The spatial directrix C s of each ribbon S consists not only of the midline M of S but also of a geodesic that is the image of the subset ]} of the referential directrix C r and connects two points on ∂S serves as a directrix for D. If S is orientable, we take C r to be parameterized by If, alternatively, S is nonorientable, we take C r to be parameterized by Notice that (14) and (15) are defined such that the referential directrix C r is parametrized by arclength. In Fig. 1, C r is depicted in red. Regardless of the orientability of S, each ruling of D intersects C r exactly once.
The spatial directrix C s is the image of C r under the deformation χ from D to S and has arclength parameterization where c is determined by (14) or (15) depending on whether S is orientable or nonorientable. Notice that, as Fig. 3 shows, the midline M of S is parametrized by and, thus, constitutes but a portion of C s . For α ∈ [0, + 2b], let f (α) denote the unit vector in the direction of the ruling that passes through c(α) chosen so that (c (α) × f (α)) · ı 3 > 0, and let θ(α) denote the angle between c (α) and f (α). Then, {c , ı 3 × c , ı 3 } is a Darboux frame for the referential directrix of D and it follows that f has the representation where since the referential directrix consists of straight line segments and ı 3 is a fixed unit vector, For each α ∈ [0, + 2b], it can be shown that it is possible to find are the end points of the ruling that passes through the point c(α) on the referential directrix C r . It follows that there is a parameter set such thatx : P → D defined bŷ covers D. In a similar way, there is a functionr : P → S of the form which covers S, where g(α) = (∇χ (x(α, 0))) f (α) is a unit vector parallel to the ruling that passes through the point d(α) on the spatial directrix C s . Notice, from (21) and (22), that the parameterizationsx andr determine the deformation χ through the requirement that if and only if x =x(α, β) and r =r(α, β).
Due to the assumed smoothness of χ , d is C 2 except at α = b and α = + b. Also, by a result due to Hartman and Wintner (1951), f and g are both C 1 except possibly at α = b and α = + b and at values of α at which flat and curved portions of S meet. Chen et al. (2018) showed that the parameterizations (21) and (22) are subject to the following constraints, which are necessary and sufficient to guarantee that χ is an isometric deformation from D to S: Notice from (19) 1 and (24) 4,5 that If we set t = d , we observe that is a unit vector orthogonal to S along its directrix C s , and that is a unit vector that is tangent to S and normal to C s . The triad {t, m, n} is a Darboux frame for C s , from which it follows that where κ g , κ n , and τ g , respectively, denote the geodesic curvature, normal curvature, and geodesic torsion of C s . Since the referential directrix C r consists of three straightline segments and χ is an isometric deformation, C s is the union of three geodesic curves. Thus, the geodesic curvature of C s obeys except at b and + b, where it is undefined. From (24) 4 , we see that the angle between t and g must be identical to the angle θ between c and f , and thus, that g admits the representation g = cos θ t + sin θ m.
Using (31) and (18) in (24) 3 , we arrive at the requirement κ n cos θ + τ g sin θ = 0 ( 3 2 ) from which we see that (31) simplifies to Moreover, we find from (28) 2 , (29), and (32) that and from (28) 3 and (32), that Since c(α)+β − (α) f (α) and c(α)+β + (α) f (α) are the end points of the ruling that passes through c(α), we find that it is possible to determine explicit representations for the functions β ± entering the definition (20) of P. To express those representations, we first let α L ∈ [b, + b] denote the largest value of α for which the ruling containing c(α L ) intersects the left-hand edge ∂D L of D and let α R ∈ [b, +b] denote the smallest value of α for which the ruling containing c(α R ) intersects the right-hand edge ∂D R of D, as depicted in Fig. 1. Notice that these values of α must satisfy α L < α R . If S is orientable, we find that and that If, alternatively, S is nonorientable, we find that and that

Natural Extension of the Parameterization of S
From (22), we notice that the functionr : P → S defined in (22) possesses a natural extension to the set When using that extension, it is important, however, to keep in mind thatr(α, β) need not be a point on S unless (α, β) is an element of P. Defining parameter sets we observe that the rulings that intersect the left-and right-hand edges ∂D L and ∂D R of D induce a bijectionα : I → K. Moreover, we see that the bijectionα : I → K can be used to express a constraint onr : P e → S. The specific forms of the bijection and the constraint on the extension ofr differ depending on the orientability of S.

Orientable Closed Ribbons
For an orientable ribbon, the bijectionα : I → K takes the form as exemplified in Fig. 1a.

Nonorientable Closed Ribbons
For a nonorientable ribbon, the bijectionα takes the form as depicted in Fig. 1b.

Consequences of Requiring that Rulings Do Not Overlap
It is not generally possible to ensure that the rulings do not intersect in the interior of S without placing a suitable restriction on θ . To determine such a restriction, choose α 1 and α 2 satisfying b ≤ α 1 ≤ α 2 ≤ + b. The points d(α 1 ) + β + (α 1 )g(α 1 ) and d(α 2 ) + β + (α 2 )g(α 2 ) are on the edge of the ribbon and are obtained by starting at the points d(α 1 ) and d(α 2 ) and proceeding along the associated rulings in the direction of g. Consider the geodesic that passes through each of these points and is orthogonal to the midline of the ribbon. Those geodesics intersect the midline at the points corresponding to the values α 1 + b cot θ(α 1 ) and α 2 + b cot θ(α 2 ) of α . Since the rulings do not intersect in the interior of the ribbon, we must have Moving all the terms to the right-hand side of the inequality (53), dividing by α 2 − α 1 , and taking the limit as α 2 goes to α 1 , we find that Similarly, considering the points d(α 1 ) + β − (α 1 )g(α 1 ) and d(α 2 ) + β − (α 2 )g(α 2 ) on the edge of the ribbon obtained by starting at the points d(α 1 ) and d(α 2 ), proceeding along the associated ruling in the direction of −g, and taking steps analogous to those leading to (54), we find that Since csc θ > 0, we may express the conditions (54) and (55), which must hold for all α 1 , in the form 0 ≤ sin 2 θ ± bθ .
Recognizing that sin 2 θ + bθ and sin 2 θ − bθ cannot be simultaneously negative, we deduce that the inequalities in (56) are equivalent to the single inequality For situations in which the directrix is a line of curvature, a condition equivalent to (57) was derived by Hornung (2011c) and was shown to imply that, in the interior of the ribbon, rulings cannot intersect locally. 4

Bending Energy
The energy stored during an isometric deformation of the unstretchable rectangular strip D to a closed ribbon S is due purely to bending, and thus, the curvature tensor of S serves as a measure of strain. For simplicity, we assume that the bending energy density depends quadratically on the curvature tensor and is isotropic. Since the Gaussian curvature vanishes identically for any isometric deformation χ : D → S, it follows that the total bending energy E of S must be of the simple form with μ > 0 being a material constant known as the bending modulus.

Consequences of Requiring that the Total Bending Energy Be Bounded
It is possible for an isometric deformation χ of D to result in a surface S with unbounded bending energy E. This can occur if, for example, a collection of rulings that meet at a point on D cover a curved portion of S. To verify this assertion, consider a portion of the surface T c covered by rulings that intersect on the boundary of S as discussed in Sect. 2.3. Let φ 1 and φ 2 be the minimum and maximum angles associated with the referential preimage T r of T c . With reference to the parameterization (12), we find that the mean curvature H of T c can be expressed as and thus, that the bending energy of T c is given by which is unbounded unless H = 0. From the foregoing argument, we infer that the bending energy E of S cannot be bounded unless each portion of S between a pair of rulings that intersect at a point on ∂S is flat. We next show that if χ is C 2 on the interior of D and the bending energy E of S = χ(D) is bounded, then χ must be C 1 on the entirety of ∂D. To establish this claim, we first collect some preliminary results. Recall that in Sect. 2.2 we argued that D can be completely covered by a continuous family of rulings. This means that there is a functionf defined on the interior of D with unit vector values and with the property thatf (x) is parallel to the ruling that passes through x. Letf ⊥ be the unit vector chosen such that it is orthogonal tof and such that (f ×f ⊥ ) · ı 3 > 0. Then, using Q = ∇χ to denote the gradient of χ , we see that the vectorŝ must be tangent to S, withĝ(x) being parallel to the ruling on S that passes through χ(x). Since the χ is an isometric deformation, | Qf | = |f | = 1 and | Qf ⊥ | = |f ⊥ | = 1 and we see from (61) thatĝ(x) andĝ ⊥ (x) are of unit magnitude. From (61), we see further that Q admits the representation Q =ĝ ⊗f +ĝ ⊥ ⊗f ⊥ .
Sincef ,f ⊥ ,ĝ, andĝ ⊥ are uniform along each ruling, we deduce from (62) that Q must also be uniform along each ruling. Consider a point x • on ∂D where rulings intersect, so that the corresponding portion of S is flat. If T r denotes the corresponding triangular region in D that is covered by the rulings that intersect at x • , then Q is uniform on T r since χ (T r ) is flat. Thus, writing Q • for the value of Q on T r and invoking the continuity of χ, we find that for any point x f ∈ T r . Moreover, for any point x ∈ D sufficiently close to x • , choose We also have the relation Notice that |x − x f | → 0 as x → x • . Then, for ε > 0 and sufficiently small, there is a compact subset D ε of the interior of D such that any ruling that intersects the set also intersects D ε . Since D ε is compact and Q is continuous, Q is uniformly continuous on D ε ; moreover, since Q is uniform along each ruling, it must be uniformly continuous on a set of the form for some η > 0. From this, (63), and (64), we thus see that from which we deduce that ∇χ (x • ) = Q • , and, thus, that, as claimed, χ is C 1 at x • .

Dimensionally Reduced Bending Energy
Using the parameterizationr of S from (22), we may express the bending energy E defined in (58) in the form To reduce E to an integral over the spatial directrix C s , we first use the decomposition (68) of [0, + 2b] to rewrite (67) as Next, utilizing (48), (49), and (37)-(38), we find that (69) can be rewritten as from which we see that Finally, noting from (22), (33), and (35) that where sin θ − βθ must be positive to ensure that rulings do not intersect in the interior of S. We infer that the innermost integral of (71) can be evaluated in closed form, and thus, that the bending energy E of S can be expressed as a single integral from α = b to α = + b: The term sin 2 θ + bθ must also be positive to ensure that rulings do not intersect in the interior of S, and thus, that the argument of the logarithm in (73) is always positive. The dimensional reduction (73) is valid for orientable and nonorientable closed ribbons since the identities (48) and (49) needed to pass from (69) to (71) hold regardless of the orientability of S, as noted at the end of Sect. 2. When using E in the form (73), it is essential to keep in mind that θ is not defined if d vanishes or does not exist. In addition, the bending energy density appearing in (73) is not defined if θ vanishes. Using the asymptotic properties of the logarithm, we find, however, that The expression on the right-hand side of (74) is the bending energy density that arises from using the dimensional reduction if the rulings are locally parallel, in which case θ = 0. Hence, the bending energy density appearing in (73) can be used even if θ vanishes with the understanding that it is continuous at any point of the directrix where that occurs.

Wunderlich Formula for the Bending Energy
Wunderlich (1962) was the first to carry out the dimensional reduction in the bending energy (58) for the problem of isometrically deforming an unstretchable rectangular strip D of length and finite width 2b to a Möbius band S with a single half-twist, leading to the expression 5 involving integration over the midline M of S, whereκ 0 is the principal normal curvature of M, as parameterized by (17), and V is defined by where θ denotes the angle between the tangent and generators of M. Invoking Euler's formula 6 for the curvature κ = |d |, we see that the principal normal curvature of M can be expressed asκ Using (76) and (77) in (75), we find that (75) is equivalent to (73).

Ribbons with Piecewise Continuous Curvature
Although we derived the dimensionally reduced bending energy (73) under the assumption that the deformation χ is C 2 , and hence, that the curvature of the closed ribbon S is continuous, (73) is well-defined and finite under weaker regularity conditions. To specify such conditions, it is first essential to select the independent variables over which (73) is to be minimized. From the structure of (73), it might seem natural to view (73) as a functional of d and θ . Recalling from (30) that cos θ = d · g and, thus, noting that θ and θ can be expressed in terms of d , d , g, and g , we instead choose d and g as the primary independent variables. On that basis, we next determine the conditions on d and g consistent with allowing the curvature of S to be piecewise continuous and recast the constraints (24) in terms of an equivalent set of constraints involving only d and g. Moreover, we demonstrate that, given d and g satisfying our relaxed regularity assumptions and our alternative constraints, it is possible to construct a closed ribbon S and an isometric deformation χ from D to a ribbon S with piecewise continuous curvature. Due to that construction, we infer that the problem of minimizing (73) over all d and g is equivalent to the problem of minimizing (58) over all isometric deformations.

Relaxed Regularity of the Spatial Directrix and Generatrix
The regularity properties of d and g, which stem from the assumed regularity of the deformation χ, are too strong to capture all physically relevant isometric deformations of a rectangular planar strip. The need for relaxed regularity is demonstrated by Sadowsky's (1930) construction of a Möbius band with a single half-twist, which yields a surface that is C 1 but only piecewise C 2 . For this reason, we hereafter relax our smoothness assumptions by assuming that d is C 1 and d is piecewise C 1 and g is continuous and piecewise C 1 . (78) We say that a function u defined on exist for 1 ≤ k ≤ n.

Reformulation of the Constraints Ensuring that the Deformation Is an Isometry
Our ultimate goal is to provide a framework for constructing isometric deformations that result in closed ribbons S that are C 1 and minimize the bending energy (58), which we have shown simplifies to (73). This dimensionally reduced expression involves the vector-valued functions d and g. When working directly with (73), it is essential to ensure that those quantities generate an isometric deformation of D, meaning that they must satisfy the constraints (24). Recalling that c is given a priori by (14) or (15), θ defined by (18) is the angle between c and f , and the parameterization (21) provides a complete covering of D, then it transpires that the constraints (24) hold if and only if d and g satisfy the alternative collection of conditions When using the conditions in (80) to establish the constraints (24) the function f is defined by (18) with θ being the angle between d and g. To verify the foregoing assertion, we begin by showing that (24) implies (80). It is clear that (80) 1,2 follow directly from (24) 1,2 . Additionally, we notice that (80) 3 follows from (25) and that (80) 4 follows from (27) and (33). To establish the reverse implication, we begin by noticing that (80) 1,2 together with the definitions of c and f imply (24) 1,2 . Moreover, we see from (18), (80) 1 , and the foregoing identification of θ with the angle between d and g that (24) 4 holds. Differentiating (24) 4 and utilizing (80) 3 , we arrive at (24) 5 . It remains only to establish (24) 3 . Notice that, to ensure the reference domain D is completely covered, β ± must be given by (36)-(37) if S is orientable or by (38)- (39) if S is nonorientable. For these expressions to be well defined, we must have Constructing the Darboux frame {t, m, n} for d as in Sect. 2.4, so that, in particular, (28) holds, we find from the definition of θ and (80) 2 that g has the form g = cos θ t + sin θ m.

Construction of a Closed Ribbon with Piecewise Continuous Curvature
We next show how to construct a closed ribbon S using vector-valued functions d and g under the assumption that they satisfy the relaxed regularity conditions (78), the constraints (80), and the inequality (57). To describe the construction, we first fix d and g satisfying (78), (80), and (57) and define S by To ensure that this defines a ribbon-like surface, it suffices to show that the parameterization defining S is locally injective on the interior of its domain. By the inverse function theorem, this is accomplished by showing that Using (57), in the form of the two equivalent inequalities (54) and (55), we can then show that (86) holds. The unit normal to S defined according to is continuous by (78). Although (78) guarantees only that S as defined by (85) is piecewise C 1 , a consequence of the underlying structure is that S is everywhere C 1 . To establish this claim, consider a point p • = d(α • ) + β • g(α • ) on S. Notice that if g is C 1 at α • , then S is C 1 at p • . Thus, we need only consider situations in which g is discontinuous at α • . Let S ± be sets defined such that Notice that S ± is, at p • , locally diffeomorphic to a half-disc. To specify the underlying diffeomorphisms, consider the plane Since S + and S − share the same normal at p • , A is tangent to both of these surfaces at p • . For any ε > 0, let B ε denote the closed ball of radius ε in A centered at p • . The ruling L • that passes through p • lies on A and, hence, divides B ε into two halves, say B + ε and B − ε , that correspond to whichever of S + or S − is closest. There is an ε > 0 and C 1 functions ϕ ± : Here, ϕ ± represents S ± near p • using the tangent space A of S ± at p • . An analytical consequence of the geometric requirement that A be tangent to S ± along L • is that ϕ ± and the natural injection of A into E have the same gradient along L • . Thus, Next, defining ϕ : B ε → S by we see that, due to the regularity of ϕ ± and (90), ϕ is C 1 . Thus, we conclude that S is C 1 at p • as claimed.

Construction of an Isometric Deformation from a Rectangular Strip to a Closed Ribbon with Piecewise Continuous Curvature
Having shown how a closed ribbon S with piecewise continuous curvature can be constructed, we next show how the underlying isometric deformation from the rectangular strip D to S can be constructed. In so doing, we must distinguish whether S is orientable or nonorientable, as determined by which alternate form, (14) or (15), of the referential directrix c is used. Granted that d and g satisfy the relaxed regularity conditions (78) and the constraints (57) and (80) (14) or (15) and that f is defined by (18), with θ being the angle between d and g, we can define the parameterization (21). The constraints (24) are satisfied since (80) hold. Thus, they induce an isometric deformation χ from D to S. We can use the inequality (57) to show that the parameterizationx of D is locally injective, just as we previously did for the parametrizationr of S and, hence, we may conclude that the deformation χ is also locally injective. Since D and S are both C 1 , a result of Myers and Steenrod (1939) ensures that χ is C 1 . It is important to note here that the deformation χ constructed as described need not be injective, but only locally injective. In the mathematics literature, χ would be called an isometric immersion. The properties of such immersions have been studied by Pakzad (2004), Müller and Pakzad (2005), and Hornung (2011a, c), among others. To ensure that the constructed χ is injective, it would be necessary to require that d and g satisfy an appropriate global constraint, the formulation and analysis of which we leave to a future work. In a related work, Halpern and Weaver (1977) discussed the restrictions on b/ that guarantee the existence of an isometric immersion from a rectangular strip of length and width 2b to a half-twist Möbius band along with the additional requirements on b/ that guarantee the existence of an isometric embedding-that is, an injective isometric emersion.

Necessary Conditions for Equilibrium
We next specify the variational problem associated with finding an equilibrium shape S of the rectangular strip D that can undergo only isometric deformations induced by one of the connecting edge conditions (10) or (11) while leaving ∂S free of traction and couple traction. To begin, we recall from the previous section that, to construct S, it suffices to find vector-valued functions d and g defined on [b, + b] and consistent with the constraints (57) and (80) that minimize the dimensionally reduced bending energy (73).
To incorporate the constraints (57) and (80), we consider the augmented bending energy functional (92) where F and G are given by 2θ sin 2 θ log sin 2 θ + bθ sin 2 θ − bθ , and Here λ, χ , γ , and ν are arclength dependent Lagrange multipliers and ω is an arclength dependent Karush-Kuhn-Tucker multiplier. The connecting edge conditions (10) and (11) must be formulated in terms of d and g, and we do this in the next two subsections, considering the orientable and nonorientable cases separately. In addition, the extent to which the ribbon S is twisted must also be specified. The edge condition (10) can be met in the absence of twist, in which case the equilibrium shape of S would be cylindrical and, by the regularity assumptions (78), smooth. It is also possible, however, to satisfy (10) in the presence of any number of full twists, assuming, as stated previously, that the half-width-to-length aspect ratio b/ is sufficiently small. An analogous statement involving half-twists applies to (11). To quantify the number of twists present in either case, we will use the linking number between the midline M of S and a perturbation of M along S. Although knotted strips are also compatible with the connecting edge conditions (10) and (11), we disregard that possibility in the present work.

Orientable Closed Ribbons
If D is isometrically deformed by holding one of its short ends fixed and subjecting the other short end to any number k full twists before satisfying (10), then the resulting orientable surface S is not a closed ribbon unless d and g satisfy Consider a perturbation M ε of the midline M of S, in the direction of g, with parameterization where, to ensure that d ε is on the surface generated by d and g, ε > 0 satisfies Then, the linking number Lk(M, M ε ) between the midline M parameterized by d and its perturbation M ε parameterized by d ε must obey The value of Lk(M, M ε ) determined by (98) does not depend on the choice of ε.
Since there are two directions in which the twist can be performed, k may be positive or negative. In the absence of twist, k = 0.
where on defining the various derivatives of F and G are given by and Setting (102) equal to zero gives a weak statement of the conditions that must hold in equilibrium.
Granted that d is piecewise C 3 and g is piecewise C 2 , we find upon integrating by parts several times that (102) can be expressed as where ϕ, ψ, ξ , η, and ζ are given by From the requirement that the first variation δ E * [d, g](δd, δ g) of E * [d, g] must vanish for all admissible variations δd and δ g, we obtain Euler-Lagrange equations in the form of ordinary differential equations which apply on the interval (b, l + b), matching conditions Notice the sign on the right-hand side of the matching condition (109) 3 is dictated by the orientability of S-the plus sign applies if S is orientable and the minus sign applies otherwise.
In addition to the equilibrium conditions (108)-(110) and the constraints (80), it is also necessary to satisfy the conditions and ω ≥ 0 (111) associated with the Karush-Kuhn-Tucker multiplier, together with either the consequence (95) of the edge condition (10) and the specification (98) of the linking number k ∈ Z, which apply if S is orientable, or the consequence (99) of the edge condition (11) and the specification (101) of the linking number k satisfying k ∈ Z\{0}, which apply if S is nonorientable. It is worth emphasizing that the equilibrium equations we obtain are written directly in terms of two vector-valued functions: the directrix d and the generatrix g, the two quantities that are necessary and sufficient to construct a surface S and a locally injective isometric deformation from D to S. This is in contrast to other equilibrium equations found in the literature. Whereas Starostin and van der Heijden (2015) write their equilibrium equations in terms of the Frenet-Serret frame of the midline and its derivatives, Dias and Audoly (2015) write their equilibrium equations using the curvature and torsion of the midline. Hornung's (2011b) equations are very different in structure since they are not traditional boundary value problems, but rather are in a form amenable to establishing regularity results.

Comments on the Matching Conditions
We next interpret the matching conditions (109).
• (109) 1 : This condition expresses the requirement that the forces at the points d(b) and d( + b) of S be equal. Moreover, these forces act to ensure that d(b) = d( + b), and hence, that the midline M is closed. • (109) 2 : Bearing in mind that d and g are orthogonal as a consequence of (80) 3 , we first observe that (109) 2 is equivalent to If the condition (112) 1 is applied in conjunction with either (95) 3 or (99) 3 , it follows that the Lagrange multiplier γ satisfies either (95) 3 or (99) 3 , respectively, with g replaced by γ . Condition (112) 2 expresses the requirement that the bending moments at the points d(b) and d( + b) are equal. Moreover, these bending moments act to ensure that d (b) = d ( + b), and hence, that the tangent to M is smooth at the point of closure. • (109) 3 : Upon noticing that F g (d , g, g ) and G g (d , g, g ) are tangent to M, we see that (109) 3 expresses the requirement that the twisting moments and at the points d(b) and d( + b) of S are of equal magnitude and orientation consistent with the orientability of S. Moreover, these twisting moments act to ensure that g(b) = g( + b), and hence, that the short edges of D are continuously joined.

Nonuniqueness of the Rulings on Flat Portions of the Ribbon
We next show that on intervals where d = 0, namely the intervals associated with the flat subregions of S, no solution of (108) and (110) and −2λt + (γ g) + (ν g × g ) + ωG d = 0, Resolving the components of (116) 2 relative to the basis {t, g, n}, we obtain three equations which can be viewed as determining ν, χ , and ω. In particular, we find from the component along g that χ = 0. By (116) 1 , the condition must hold for all α ∈ (α 1 , α 2 ) and any choice of the multiplier λ compatible with the jump conditions (117) 1 and (117) 3 . Since λ is arbitrary on the open interval (α 1 , α 2 ), g cannot be determined uniquely from (118). This, however, should not be surprising since the bending energy density in (73) vanishes when d = 0, and thus, the distribution of g on a flat portion of M does not influence the dimensionally reduced bending energy (73).

Summary and Discussion
Our analysis is the first study to rigorously demonstrate that the Wunderlich functional (2) is a valid dimensional reduction in the bending energy energy (1) for any closed ribbon S obtained by an isometric deformation of a flat rectangular strip S of length and width 2b, regardless of the orientability of S and the number of twists that S possesses. That finding is, however, predicated on the imposition of the connecting edge conditions (10) and (11), which determine whether or not S is orientable, and by the corresponding conditions (98) and (101) on the self-linking number of the midline M of S.
An essential feature of our approach to establishing the validity of the Wunderlich functional (2) involves the referential and spatial directrices. The referential directrix involves not only the midline of the rectangle D but also segments of its short edges. Hence, the spatial directrix, which is parametrized by d, involves not only the midline M of S but also the curve on S that is created by joining the short edges of D consistent with whichever of the connecting edge conditions (98) and (101) applies. This ensures that both D and S are completely covered, and thus, that χ is surjective. Choosing M alone as the spatial directrix, as done by Dias and Audoly (2015), does not generally suffice to ensure the satisfaction of that essential requirement.
We find that a closed ribbon S obtained by minimizing the Wunderlich functional (2) represents the shape determined by a locally injective isometric deformation χ of a flat rectangular strip D only if certain additional constraints are met. Those constraints can be expressed in terms of four scalar conditions and one inequality involving the directrix d and generatrix g of S and their derivatives with respect to arclength along M. This leads us to a constrained variational problem, incorporating four scalar Lagrange multipliers and one Karush-Kuhn-Tucker multiplier, each of which is generally dependent on arclength along M, for the augmented energy functional E * defined in (92). The local injectivity condition is an important ingredient in the literature on isometric immersions with finite bending energy. However, the works in that literature do not include derivations of dimensionally reduced bending energy or the associated boundary value problems for determining the deformation. On the other hand, the local injectivity condition is mostly overlooked in the remainder of the pertinent literature.
Our approach is contingent upon the requirement that the bending energy be finite. Although we allow rulings to intersect on the boundary ∂S of the ribbon, we do not consider situations where rulings intersect on the interior of S. Granted these provisions, we demonstrate that knowledge of the spatial directrix d and spatial generatrix g suffices to ensure that we can construct a C 1 surface S and corresponding deformation χ that is sufficiently regular to guarantee that the curvature of S is piecewise continuous.
The Euler-Lagrange equations corresponding to our constrained variational problem consist not only of ordinary differential equations for the directrix d and generatrix g of S but also of matching conditions that apply at the point on the midline M of S where the midline of D is joined consistent with whichever of the connecting edge conditions (98) and (101) applies and of jump conditions which apply at points on M where the curvature of S is discontinuous. Sadowsky's (1930) construction of a Möbius band with a single half twist, which yields a surface that is C 1 but only piecewise C 2 , demonstrates that there do exist configurations of half-twist ribbons that are C 1 but only piecewise C 2 . This justifies our regularity assumptions on d and g. We make no claims about the regularity of the equilibrium configurations, which may even be C 2 . The derivation of our Euler-Lagrange equations uses standard techniques and, in particular, does not require recourse to more sophisticated methods such as the higher-order variational approach of Gay-Balmaz et al. (2012).
To obtain the final form of their Euler-Lagrange equations for the equilibrium shapes of half-twist Möbius bands, Starostin andvan der Heijden (2007, 2015) restricted attention to shapes that possess 180-degree rotational symmetry about a particular axis. 7 They also considered the equilibrium shapes of Möbius bands with multiple half-twists and, to write down the final equations, assumed that the associated shapes posses certain specific symmetries. Although some physical constructions of half-twist Möbius bands appear to possess such symmetry, we are not aware of a proof demonstrating that such symmetry is always achieved in equilibrium. Since our framework does not involve symmetry assumptions, it provides, in particular, a basis for developing numerical methods for investigating the existence and stability of shapes with and without symmetry.
The problem considered here is a special version of a more general one in which displacements and dead loads, constituted by tractions and couple tractions, are applied on the boundary of the ribbon. In such problems, the underlying variational principle must include the potential energy of the dead loads. The task of generalizing the dimensional reduction argument leading to (73) when dead loads are present is the subject of a forthcoming work. Another important direction, which is also in progress, involves extending our results to deal with reference shapes that are not rectangular but which can be smoothly ruled. The most challenging extension of our work arises when a smooth ruling is not possible. This can occur even if the reference shape is rectangular, but its four corners are subjected to mild displacements perpendicular to the plane in which it resides.