Pillow Box Design

Abstract

This paper focuses on packaging design using origami techniques, specifically designs incorporating curves, known as pillow boxes. While conventional paper packaging boxes are typically cuboid, pillow box designs include curved surfaces, offering both aesthetic and practical advantages. This study analyzes the specific curved folds of pillow boxes, clarifying the fundamental geometric condition these curves must meet. Additionally, it proposes new design variations for pillow boxes based on the condition. The relationship between the shape of the folds and the volume of the final three-dimensional shape is also explored. This research extends the boundaries of functionality and aesthetics in origami design and explores new possibilities in packaging solutions.

1 Introduction

One of the notable applications of origami techniques in practical use can be observed in the design of packaging boxes. Boxes made from paper come in many variations, ranging from simplistic designs to those that are intricately crafted. Among these, some even have curved folds. The pillow box can be cited as one of the most simplistic yet popular packaging designs that incorporate curved folds.

Figure 1a shows an example of a pillow box, while Figure 1b presents its crease pattern, i.e. unfolded layout. Curved creases can also be observed in the French fry containers used at hamburger chain restaurants. The pillow box can be considered as a form composed of two of these containers joined together. Different three-dimensional shapes of boxes can be made depending on the shape of the curved folds of the pillow box shown in Figure 1b. However, there are specific conditions that these curves must fulfill for a design to function as a box. This paper aims to clarify the geometrical constraint of the curves used for pillow boxes, with the objective of expanding the design variations. Additionally, the latter part of the paper explores the relationship between the shape of the folds and the volume of the box.

The crease pattern of the pillow box targeted in this paper, as depicted in Figure 1b, exhibits symmetry both vertically and horizontally. After folding, the three-dimensional shape possesses symmetry along the z-axis and y-axis directions. Upon folding, the sides overlap, resulting in a two-layered structure.

As shown in Figure 1c (a quarter of the three-dimensional shape), the shape after folding consists of both the top and side surfaces being cylindrical. The rulings that constitute the top surface are horizontal, while those forming the side surface are vertical. By mirror-inverting a part of the top surface over a plane inclined at 45°, the side surface can be formed. This means that the three-dimensional folds of the pillow box are planar curves.

Refer to caption
Figure 1: Photo of a pillow box (a), its crease pattern (b) and one-fourth of the entire three-dimensional shape.

2 Related Work

2.1 Curved origami

The curved surfaces created by bending paper are limited to a specific type of ruled surface, where a straight line moves through space, characterized by a Gaussian curvature of zero; these are known as developable surfaces. Non-planar developable surfaces are confined to conical, cylindrical, and tangential surfaces. However, by incorporating folds into these surfaces, it is possible to create a wide variety of shapes.

Many studies have been conducted on the geometry of developable surfaces with curved folds. Fuchs and Tabachnikov [Fuchs and Tabachnikov 99] have clarified the relationship between three-dimensional spatial curves that become folds, their corresponding curves mapped onto a crease pattern, and the angles of the folds. Tachi streamlined this into clear equations [Tachi 11], paving the way for shape design with curved folds. Demaine et al. have explored the geometric properties of developable surfaces with curved folds, uncovering many facts [Demaine et al. 15]. When the fold is a planar curve, it is characterized by equal folding angles on both sides of the fold. This can be seen as dividing the developable surface into two parts with a plane and then mirror-inverting one side. This principle has been widely used in the creation of origami art by David Huffman [Demaine et al. 10]. Utilizing this principle of mirror inversion, Mitani and Igarashi developed a tool for interactively designing shapes that include planer curved folds [Mitani and Igarashi 11].

2.2 Volume of the pillow box

The volume of the pillow box is discussed in the latter part of this paper. Related to this is a discussion that took place in 1997-1998 on a web page named The Geometry Junkyard (which remains in the web archive [web 97]), in which various contributors discussed the geometrical problem of maximising the volume inside two square sheets of length one when sealed along their edges. At the end of the web page compiled by Andrew Kepert on these entitled ’Teabag problem’ (which also remains in the web archive [Kepert 98]), there is a description of the shape of the pillow box discussed in this paper, where the maximum volume is approximately 0.1703844172 when θ𝜃\thetaitalic_θ is approximately 1.047 radians under the condition that the surface is along an arc of angle 2θ2𝜃2\theta2 italic_θ (the value we have obtained with our optimisation method is larger than this.)

3 Variation of Shapes

3.1 Symmetric Shape with Respect to the Horizontal Plane

The shape targeted in this paper is discretized into a collection of elongated planar rectangles, as illustrated in Figure 2a.

Refer to caption
Figure 2: Discrete model of the pillow box.

Due to the geometric constraint that the area of \triangleSPQ in Figure 2b must be non-negative, the angle θ𝜃\thetaitalic_θ needs to fall within the range specified by the following equation.

π/4θπ/4𝜋4𝜃𝜋4-\pi/4\leq\theta\leq\pi/4- italic_π / 4 ≤ italic_θ ≤ italic_π / 4 (1)

This fact is explained as follows:

Figure 2b is an enlargement of the planar rectangles on both sides of edge PR from Figure 2a. At points P and R, edges are connected that are parallel to the z𝑧zitalic_z-axis and y𝑦yitalic_y-axis. The distance w𝑤witalic_w is the distance between the vertical lines passing through P and R. Q is the foot of the perpendicular dropped from P to the vertical edge passing through R. S is the foot of the perpendicular dropped from R to the horizontal edge passing through P. If RPQ𝑅𝑃𝑄\angle RPQ∠ italic_R italic_P italic_Q is denoted as θ𝜃\thetaitalic_θ, then PS and RQ are equal to ΔwtanθΔ𝑤𝜃{\Delta}w\tan\thetaroman_Δ italic_w roman_tan italic_θ. Since \triangleSPQ is a right triangle lying horizontally, the following equation holds based on the Pythagorean theorem:

PS2+SQ2𝑃superscript𝑆2𝑆superscript𝑄2\displaystyle PS^{2}+SQ^{2}italic_P italic_S start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_S italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT =PQ2absent𝑃superscript𝑄2\displaystyle=PQ^{2}= italic_P italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT (2)
(Δwtanθ)2+SQ2superscriptΔ𝑤𝜃2𝑆superscript𝑄2\displaystyle({\Delta}w\tan{\theta})^{2}+SQ^{2}( roman_Δ italic_w roman_tan italic_θ ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_S italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT =Δw2absentΔsuperscript𝑤2\displaystyle={\Delta}w^{2}= roman_Δ italic_w start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT
(tanθ)2superscript𝜃2\displaystyle(\tan{\theta})^{2}( roman_tan italic_θ ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT =1SQ2Δw2absent1𝑆superscript𝑄2Δsuperscript𝑤2\displaystyle=1-\frac{SQ^{2}}{{\Delta}w^{2}}= 1 - divide start_ARG italic_S italic_Q start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG roman_Δ italic_w start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG

Therefore

|tanθ|1𝜃1|\tan{\theta}|\leq 1\\ | roman_tan italic_θ | ≤ 1 (3)

and the following equation holds

π/4θπ/4𝜋4𝜃𝜋4-\pi/4\leq\theta\leq\pi/4- italic_π / 4 ≤ italic_θ ≤ italic_π / 4 (4)
Refer to caption
Figure 3: Out-of-the-ordinary pillow box design. Polyhedral model and curved surface model.

The value of θ𝜃\thetaitalic_θ corresponds to the angle of the tangent to the fold curve in the crease pattern shown on the right side of Figure 2. Consequently, the constraint for the curve of the fold is that the absolute value of the angle of its tangent must be 45 degrees or less.

Under this condition, when fold curves are created on the crease pattern, a pillow box can be made by folding them. Figure 3 shows two examples of pillow boxes created with designed folds while satisfying the constraints. The left example features a polyline fold, while the right one is designed with smooth curved folds.

The constraint on the shape of the curve is simple, and under this constraint, it was possible to create various shapes of pillow boxes as shown in Figure 4.

Refer to caption
Figure 4: Various pillow box designs.

3.2 Asymmetric Shape with Respect to the Horizontal Plane

In previous designs of pillow boxes, a constraint was imposed such that their shape was symmetrical with respect to the xy𝑥𝑦xyitalic_x italic_y-plane, meaning that the rulings composing the side surface were parallel to the z𝑧zitalic_z-axis. This section demonstrates that by relaxing this constraint, it is possible to increase the variety of designs further.

Figure 5 presents a sectional view of a pillow box parallel to the yz𝑦𝑧yzitalic_y italic_z-plane. As previously mentioned, the rulings that constitute the top surface, side surface, and bottom surface are all parallel to the yz𝑦𝑧yzitalic_y italic_z-plane. Therefore, in this sectional view, the reflection planes R1subscript𝑅1R_{1}italic_R start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and R2subscript𝑅2R_{2}italic_R start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT appear to reflect the rulings.

Refer to caption
Figure 5: Cross-sectional view of a pillow vox parallel to the yz𝑦𝑧yzitalic_y italic_z-plane, where rulings appear reflected by reflection planes

Assuming both the top and bottom surfaces are horizontal, the following relationships hold for the angles α1subscript𝛼1\alpha_{1}italic_α start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT, α2subscript𝛼2\alpha_{2}italic_α start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT, θ1subscript𝜃1\theta_{1}italic_θ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT, θ2subscript𝜃2\theta_{2}italic_θ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT depicted in the figure:

θ1+θ2=180\tcdegree,α1=θ22,α2=θ12.formulae-sequencesubscript𝜃1subscript𝜃2180\tcdegreeformulae-sequencesubscript𝛼1subscript𝜃22subscript𝛼2subscript𝜃12\theta_{1}+\theta_{2}=180\tcdegree,\,\alpha_{1}=\frac{\theta_{2}}{2},\,\alpha_% {2}=\frac{\theta_{1}}{2}.\\ italic_θ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_θ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = 180 , italic_α start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = divide start_ARG italic_θ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG , italic_α start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = divide start_ARG italic_θ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG . (5)

Thus, determining one of these values automatically determines the remaining values.

In the designs discussed in the previous section, both α1subscript𝛼1\alpha_{1}italic_α start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and α2subscript𝛼2\alpha_{2}italic_α start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT were set at 45 degrees, and both θ1subscript𝜃1\theta_{1}italic_θ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and θ2subscript𝜃2\theta_{2}italic_θ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT were 90 degrees. However, setting different values allows for the creation of pillow box shapes where the side surface is not perpendicular to the xy𝑥𝑦xyitalic_x italic_y-plane.

In the examples shown in Figures 3 and 4, curved creases were drawn on the development and then physically folded. However, calculating the shape of curves on the development using the approach in this section is challenging. Therefore, by initially creating a cylindrical top surface and applying two mirror transformations corresponding to the reflections at R1subscript𝑅1R_{1}italic_R start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT and R2subscript𝑅2R_{2}italic_R start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT in Figure 5, it is possible to design pillow boxes that are not symmetrical with respect to the xy𝑥𝑦xyitalic_x italic_y-plane as shown in Figure 6.

Refer to caption
Figure 6: Designs of pillow boxes not symmetrical with respect to the xy𝑥𝑦xyitalic_x italic_y-plane.

4 Relationship between Curve Shape and Volume

As seen in the previous section, the shape of the folds in a pillow box affects its form, and consequently, its volume also changes. An intriguing question arises: what is the shape of the fold curve that maximizes this volume, and what is the volume in that case? In the followings, this question will be discussed in detail. This exploration delves into the relationship between the design of the crease pattern and the functional capacity of the pillow box.

4.1 Preliminary

Here we revisit and organize the conditions for the development of the pillow box we are focusing on.

Firstly, whereas previously it was assumed that the side surfaces overlapped in two layers, we will now consider the side surface to consist of a single layer. Consequently, the development becomes a rectangle, as shown in the left side of Figure 7 (which is positioned as if rotated 90 degrees from the development in Figure 1). This implies that a pillow box can also be created by adding curved creases to a generally rectangular envelope.

Refer to caption
Figure 7: Curves serving as creases in a pillow box.

In the discussions that follow, we will assume that the curves c1subscript𝑐1c_{1}italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT to c4subscript𝑐4c_{4}italic_c start_POSTSUBSCRIPT 4 end_POSTSUBSCRIPT defining the shape of the pillow box are all of the same shape. Additionally, each curve is symmetrical left to right (with curve c1subscript𝑐1c_{1}italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT being symmetrical about the centerline of the figure). It can be easily understood from the following discussion that to find a curve that maximizes volume, it suffices to examine curves that are symmetric left to right; If the shape of the curve is not symmetrical, by comparing the volumes of the left and right sides of the solid and selecting the side that does not have a smaller volume, then adjusting the shape of the curve on the opposite side to mirror the selected side, the volume obtained after this operation will be equal to or greater than the volume of the original shape. Therefore, it is impossible for an asymmetrical shape to yield a unique maximum volume.

Furthermore, let’s assume as a premise that the size of the original envelope is a rectangle with a width of 1 and a height of 22\sqrt{2}square-root start_ARG 2 end_ARG. It’s important to note that if this ratio is different, the shape of the curve that maximizes the volume will also change. For example, if the height is infinitely large, the volume is maximized when the cross-sectional shape is a circle.

The fold curve on the development will be represented by the function f(u)𝑓𝑢f(u)italic_f ( italic_u ) when plotted in an orthogonal uv𝑢𝑣uvitalic_u italic_v-coordinate system, as depicted in the right side of Figure 7. This function f(u)𝑓𝑢f(u)italic_f ( italic_u ) is symmetrical with respect to the line u=12𝑢12u=\frac{1}{2}italic_u = divide start_ARG 1 end_ARG start_ARG 2 end_ARG and satisfies f(0)=f(1)=0𝑓0𝑓10f(0)=f(1)=0italic_f ( 0 ) = italic_f ( 1 ) = 0. Additionally, based on the condition derived in the previous section, the function must satisfy the condition |f(u)|1superscript𝑓𝑢1|f^{\prime}(u)|\leq 1| italic_f start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ( italic_u ) | ≤ 1.

In the followings, we will consider examples where the cross-section in the xz𝑥𝑧xzitalic_x italic_z-plane is a circle, a rectangle, or a rhombus, as these shapes allow for a simpler calculation of the volume. After that, we will consider the shape with an arched cross-section commonly seen in pillboxes.

4.2 Circular cross section

When the cross-section is a circle, given that the circumference of the circle is 2, the radius r𝑟ritalic_r is 1π1𝜋\frac{1}{\pi}divide start_ARG 1 end_ARG start_ARG italic_π end_ARG. The function f(u)𝑓𝑢f(u)italic_f ( italic_u ) is expressed as follows:

f(u)=sin(πu)π𝑓𝑢𝜋𝑢𝜋f(u)=\frac{\sin(\pi u)}{\pi}italic_f ( italic_u ) = divide start_ARG roman_sin ( italic_π italic_u ) end_ARG start_ARG italic_π end_ARG (6)

The graph in the left side of Figure 8 represents the fold curve and the figure in the right side represents quarter of the resulting three-dimensional shape.

Refer to caption
Figure 8: The shape of one quarter of a pillow box with a circular cross-section.

The volume of this shape can be calculated as a quarter of a cylinder with a radius of 1π1𝜋\frac{1}{\pi}divide start_ARG 1 end_ARG start_ARG italic_π end_ARG and a length of 22\sqrt{2}square-root start_ARG 2 end_ARG, from which the volumes of a part

02r(g(x))2𝑑xsuperscriptsubscript02𝑟superscript𝑔𝑥2differential-d𝑥\int_{0}^{2r}\left(g(x)\right)^{2}\ dx∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 italic_r end_POSTSUPERSCRIPT ( italic_g ( italic_x ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d italic_x (7)

is subtracted. Since g(x)𝑔𝑥g(x)italic_g ( italic_x ) represents a semicircle centered at x=r𝑥𝑟x=ritalic_x = italic_r with radius r𝑟ritalic_r, it can be expressed as follows:

g(x)=r2(rx)2=2rxx2𝑔𝑥superscript𝑟2superscript𝑟𝑥22𝑟𝑥superscript𝑥2g(x)=\sqrt{r^{2}-(r-x)^{2}}=\sqrt{2rx-x^{2}}italic_g ( italic_x ) = square-root start_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - ( italic_r - italic_x ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG = square-root start_ARG 2 italic_r italic_x - italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG (8)

Therefore, the volume of a pillow box with a circular cross section can be calculated as follows.

2r2π402r(g(x))2𝑑x0.278152superscript𝑟2𝜋4superscriptsubscript02𝑟superscript𝑔𝑥2differential-d𝑥0.27815\sqrt{2}r^{2}\pi-4\int_{0}^{2r}\left(g(x)\right)^{2}\ dx\approx 0.27815square-root start_ARG 2 end_ARG italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_π - 4 ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 italic_r end_POSTSUPERSCRIPT ( italic_g ( italic_x ) ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d italic_x ≈ 0.27815 (9)

4.3 Rectangular cross section

Assuming the cross-section is rectangular, by adding fold lines at a 45-degree angle from both endpoints as shown in Figure 9, a rectangular box shape can be created in the manner of a caramel box. As depicted in the figure, if the distance from the endpoints to the rectangular prism is denoted as hhitalic_h, then the height, width, and depth of the rectangular prism are 2h22h2 italic_h, 12h121-2h1 - 2 italic_h, and 22h22\sqrt{2}-2hsquare-root start_ARG 2 end_ARG - 2 italic_h respectively, thus the volume can be expressed by the following equation.

Refer to caption
Figure 9: The shape of one quarter of a pillow box with a rectangular cross-section.
2h(12h)(22h)212222h(1-2h)(\sqrt{2}-2h)2 italic_h ( 1 - 2 italic_h ) ( square-root start_ARG 2 end_ARG - 2 italic_h ) (10)

This expression reaches its maximum value

1/27(422622+632)0.243691274226226320.243691/27(4-\sqrt{2}-2\sqrt{6-2\sqrt{2}}+6\sqrt{3-\sqrt{2}})\approx 0.243691 / 27 ( 4 - square-root start_ARG 2 end_ARG - 2 square-root start_ARG 6 - 2 square-root start_ARG 2 end_ARG end_ARG + 6 square-root start_ARG 3 - square-root start_ARG 2 end_ARG end_ARG ) ≈ 0.24369

when

h=16+132326.16132326h=\frac{1}{6}+\frac{1}{3\sqrt{2}}-\frac{\sqrt{3-\sqrt{2}}}{6}.italic_h = divide start_ARG 1 end_ARG start_ARG 6 end_ARG + divide start_ARG 1 end_ARG start_ARG 3 square-root start_ARG 2 end_ARG end_ARG - divide start_ARG square-root start_ARG 3 - square-root start_ARG 2 end_ARG end_ARG end_ARG start_ARG 6 end_ARG .

4.4 Rhombic cross section

If the cross section is a rhombus as shown in Figure 10, the volume can be expressed by the following equation, where hhitalic_h is the height of the triangle formed by the fold line and the u𝑢uitalic_u-axis,

Refer to caption
Figure 10: The shape of one quarter of a pillow box with a rhombic cross-section.
4(122wh23h2w)4122𝑤23superscript2𝑤4(\frac{1}{2}\sqrt{2}wh-\frac{2}{3}h^{2}w)4 ( divide start_ARG 1 end_ARG start_ARG 2 end_ARG square-root start_ARG 2 end_ARG italic_w italic_h - divide start_ARG 2 end_ARG start_ARG 3 end_ARG italic_h start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_w )

where

w=14h2.𝑤14superscript2w=\sqrt{\frac{1}{4}-h^{2}}.italic_w = square-root start_ARG divide start_ARG 1 end_ARG start_ARG 4 end_ARG - italic_h start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG .

Since the maximum volume could not be analytically determined even by using symbolic computation software, it was confirmed that the maximum volume, approximately 0.24351, occurs when the value of hhitalic_h is approximately 0.30547 by a numerical approach.

4.5 Arched cross section

To address a more general arch-shaped top surface, here the segment of the crease curve within the range u[0,1/2]𝑢012u\in[0,1/2]italic_u ∈ [ 0 , 1 / 2 ] is represented by a quadratic Bézier curve. Given that the crease curve is symmetric, the segment for u[1/2,1]𝑢121u\in[1/2,1]italic_u ∈ [ 1 / 2 , 1 ] is the reflection of this quadratic Bézier curve about u=1/2𝑢12u=1/2italic_u = 1 / 2. The shape of this Bézier curve is determined by the positions of control points P0subscript𝑃0P_{0}italic_P start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, P1subscript𝑃1P_{1}italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT, and P2subscript𝑃2P_{2}italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT, as depicted in Figure 11. Since the endpoints of the crease curve coincide with the origin of the uv𝑢𝑣uvitalic_u italic_v-coordinates, the position of P0subscript𝑃0P_{0}italic_P start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is fixed at the origin. Additionally, the u𝑢uitalic_u-coordinate of P2subscript𝑃2P_{2}italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT is 1/2121/21 / 2.

Refer to caption
Figure 11: Crease curve defined by Bézier curve and discrete model.

Under this setup, the coordinates of each control point can be represented as follows, with the v𝑣vitalic_v-coordinate value of P2subscript𝑃2P_{2}italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT being denoted as hhitalic_h:

P0(0,0),P1(a,b),P2(12,h)subscript𝑃000subscript𝑃1𝑎𝑏subscript𝑃212P_{0}(0,0),P_{1}(a,b),P_{2}(\frac{1}{2},h)italic_P start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( 0 , 0 ) , italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_a , italic_b ) , italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( divide start_ARG 1 end_ARG start_ARG 2 end_ARG , italic_h ) (11)

In essence, the shape of the curve is determined by three variable a𝑎aitalic_a, b𝑏bitalic_b, and hhitalic_h. In this case, the Bézier curve is expressed by the following equation using the parameter t[0,1]𝑡01t\in[0,1]italic_t ∈ [ 0 , 1 ]:

{u(t)=2a(1t)t+12t2v(t)=2b(1t)t+ht2cases𝑢𝑡2𝑎1𝑡𝑡12superscript𝑡2otherwise𝑣𝑡2𝑏1𝑡𝑡superscript𝑡2otherwise\begin{cases}u(t)=2a(1-t)t+\frac{1}{2}t^{2}\\ v(t)=2b(1-t)t+ht^{2}\end{cases}{ start_ROW start_CELL italic_u ( italic_t ) = 2 italic_a ( 1 - italic_t ) italic_t + divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_CELL start_CELL end_CELL end_ROW start_ROW start_CELL italic_v ( italic_t ) = 2 italic_b ( 1 - italic_t ) italic_t + italic_h italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_CELL start_CELL end_CELL end_ROW (12)

The three-dimensional shape obtained from folding along this curve is discretized, as shown in the right side of Figure 11. By calculating the sum of the volumes of each discrete element, the total volume of the shape can be determined. When the discretization interval ΔtΔ𝑡\Delta troman_Δ italic_t is set to 0.001, it was found that the maximum volume, approximately 0.294436, is achieved when a=b=0.2731𝑎𝑏0.2731a=b=0.2731italic_a = italic_b = 0.2731, and c=0.2544𝑐0.2544c=0.2544italic_c = 0.2544. To obtain these values, the SLSQP algorithm, which is widely used to solve constrained optimisation problems, was used. The optimise function optimize. minimize from the scipy library for Python was used for this purpose.

Furthermore, as an attempt to achieve larger values, experiments were conducted to optimize the crease curve using two representations: a cubic Bézier curve (as shown in Figure 12(a)) and a fine polyline (as shown in Figure 12(b)).

Refer to caption
Figure 12: Crease curve defined by a cubic-Bézier curve and a fine polyline.

For the cubic Bézier curve, the shape is defined by four control points from P0subscript𝑃0P_{0}italic_P start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT to P3subscript𝑃3P_{3}italic_P start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT using five variables: a𝑎aitalic_a, b𝑏bitalic_b, c𝑐citalic_c, d𝑑ditalic_d, and hhitalic_h, each represented by the following coordinate values:

P0(0,0),P1(a,b),P2(c,d),P3(12,h).subscript𝑃000subscript𝑃1𝑎𝑏subscript𝑃2𝑐𝑑subscript𝑃312P_{0}(0,0),P_{1}(a,b),P_{2}(c,d),P_{3}(\frac{1}{2},h).italic_P start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( 0 , 0 ) , italic_P start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_a , italic_b ) , italic_P start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_c , italic_d ) , italic_P start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT ( divide start_ARG 1 end_ARG start_ARG 2 end_ARG , italic_h ) . (13)

The curve is expressed by the following equation using the parameter t[0,1]𝑡01t\in[0,1]italic_t ∈ [ 0 , 1 ]:

{u(t)=3(1t)2ta+3(1t)t2c+12t3v(t)=3(1t)2tb+3(1t)t2d+t3hcases𝑢𝑡3superscript1𝑡2𝑡𝑎31𝑡superscript𝑡2𝑐12superscript𝑡3otherwise𝑣𝑡3superscript1𝑡2𝑡𝑏31𝑡superscript𝑡2𝑑superscript𝑡3otherwise\begin{cases}u(t)=3(1-t)^{2}ta+3(1-t)t^{2}c+\frac{1}{2}t^{3}\\ v(t)=3(1-t)^{2}tb+3(1-t)t^{2}d+t^{3}h\end{cases}{ start_ROW start_CELL italic_u ( italic_t ) = 3 ( 1 - italic_t ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_t italic_a + 3 ( 1 - italic_t ) italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_c + divide start_ARG 1 end_ARG start_ARG 2 end_ARG italic_t start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_CELL start_CELL end_CELL end_ROW start_ROW start_CELL italic_v ( italic_t ) = 3 ( 1 - italic_t ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_t italic_b + 3 ( 1 - italic_t ) italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d + italic_t start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_h end_CELL start_CELL end_CELL end_ROW (14)

By discretizing this curve in the same manner as before and calculating the volume, and solving it as an optimization problem with five variables, the maximum value of 0.295448 was obtained when a=b=0.1125𝑎𝑏0.1125a=b=0.1125italic_a = italic_b = 0.1125, c=d=0.2526𝑐𝑑0.2526c=d=0.2526italic_c = italic_d = 0.2526, and h=0.25430.2543h=0.2543italic_h = 0.2543.

For the fine polyline, a polyline with 1000 segments was used (setting the N𝑁Nitalic_N value in Figure 12(b) to 1000). Starting from the origin as the 0-th point, the u𝑢uitalic_u-coordinate of the i𝑖iitalic_i-th point pisubscript𝑝𝑖p_{i}italic_p start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT was set as i/2000𝑖2000i/2000italic_i / 2000, and the v𝑣vitalic_v-coordinate was determined through optimization. This approach treated the problem as an optimization with 1000 variables. As a result, a value of 0.295449 was obtained for the shape shown in Figure 13.

4.6 Summary of evaluation

Throughout this study, the maximum volumes for various cross-sectional shapes were determined. The results are summarized in Table 1. Among the easily creatable shapes—circle, rectangle, and diamond—the circular cross-section was found to have the largest volume.

For the general pillow box shape with an arch-shaped top surface, using a cubic Bézier curve yielded a maximum volume approximately 0.3% larger than when using a quadratic Bézier curve. Furthermore, employing a 1000-segment polyline allowed for a slight increase in volume, but the difference compared to using a cubic Bézier curve was extremely small, at only 0.0003%.

These findings suggest that the values obtained here are likely very close to the upper limit of the volume that can be realized with a pillow box made from an envelope of width 1 and length 22\sqrt{2}square-root start_ARG 2 end_ARG.

Using this method, the maximum volume value was also determined when the sheet was a square with a length of one side. This was 0.174628, as shown in the table. This value is 2.49% greater than the value of 0.1703844172 calculated by Andrew Kepert as presented in Section 2.2.

Figure 13 depicts the 1000-segment polyline that maximizes volume. The cubic Bézier curve that maximizes volume is almost identical in shape to this polyline, to the extent that if drawn on the same graph, the lines would overlap and appear as a single line, making it impossible to distinguish between them.

Table 1: Cross-sectional shapes and maximum volumes
Cross-sectional shape Maximum volume
rhombus 0.243507
rectangle 0.243692
circle 0.278150
arch (quad-Bézier) 0.294436
arch (cubic-Bézier) 0.295448
arch (polyline with 1000 segments) 0.295449 (0.174628 )

* The value when the vertical dimension is set to 1.

Refer to caption
Refer to caption
Figure 13: Volume-maximising curve (1000-segment polyline) and the 3d model (100-segment polyline).

5 Conclusions

In this paper, we focused on the fold curves, a critical element in determining the shape of a pillow box, and clarified their geometric constraints, demonstrating that a variety of shapes can be created within these constraints. We presented various design variations under the constraint of top-bottom symmetry, showcasing the diversity achievable in pillow box designs. Additionally, we demonstrated the feasibility of asymmetric top-bottom designs, expanding the scope of creativity in pillow box aesthetics. This approach enhances the ability to freely control the shape of the curves, thereby broadening the range of design possibilities. As a result, it becomes possible to create more complex and aesthetically pleasing pillow boxes.

Another significant aspect of this research is the focus on optimizing the volume of pillow boxes. We addressed this by examining cases with circular, rectangular, and rhombus cross-sectional shapes, providing insights into determining the optimal dimensions for maximizing volume within design constraints, using both analytical and numerical methods. This approach allowed for understanding of the relationship between design parameters and the functional capacity of the pillow boxes.

Furthermore, for more general pillow box designs with an arch-shaped top surface, we applied quadratic Bézier curve, cubic Bézier curve and fine-polyline to the creation of folding curves. Through numerical calculations, we determined the shape of the curve that maximizes the volume. The results of this study confirm that the common pillow box shape with an arch-shaped top surface can be larger in volume than those with a circular, rectangular or rhombic cross-section. The shape of the curve that maximises the volume was also obtained by optimisation calculations. For the value of the maximum volume when the sheet is square, it was confirmed that the proposed method can yield a shape with a value 2.49% larger than that calculated by Andrew Kepert [Kepert 98].

An approach to analytically determine the curve that maximises the volume is currently being worked on by Miyuki Koiso and the curve has been elucidated to be an elastic curve [Koiso 23]. The specific type of elastic curve has not yet been clarified, but it is not expected to depart significantly from the shape and maximum volume values obtained in this paper.

In conclusion, this study not only deepens the understanding of pillow box design but also expands the possibilities of origami-based packaging design by discussing both aesthetic appeal and practical utility.

Although we have focused on exploring the maximum volume of pillow boxes, if we consider the broader context of creating a three-dimensional object from two glued rectangular sheets of paper, the pillow box design might not always be the optimal choice. A cushion made of a square bag filled to the brim with cotton, or an inflated plastic balloon, will have a concave shape on all four sides while the centre is inflated. This shape cannot be easily made from paper and requires intricate folds. Regarding the maximum volume achievable without limiting the design to pillow boxes, information can be found on the English Wikipedia under the paper bag problem or teabag problem. This problem investigates the maximum volume that can be achieved with a given amount of material, a fundamental question in packaging design. The maximum volume is proposed to be represented by an approximate formula.

V=w3(h/(πw)0.142(110(h/w)))𝑉superscript𝑤3𝜋𝑤0.1421superscript10𝑤V=w^{3}(h/(\pi w)-0.142(1-10^{(-h/w)}))italic_V = italic_w start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT ( italic_h / ( italic_π italic_w ) - 0.142 ( 1 - 10 start_POSTSUPERSCRIPT ( - italic_h / italic_w ) end_POSTSUPERSCRIPT ) ) (15)

When the values of h=22h=\sqrt{2}italic_h = square-root start_ARG 2 end_ARG and w=1𝑤1w=1italic_w = 1 are substituted into the formula, the calculated volume is approximately 0.313629. The maximum volume of 0.295449 for the pillow box design demonstrated in this paper corresponds to about 94.2% of this theoretical upper limit. This shows that the pillow box design is quite efficient and close to the theoretical limit, even though it can be realised with simple folds.

From a mathematical perspective, future work could include analytically determining the curve that maximizes volume. However, from a practical standpoint, verifying the strength and ease of assembly of these designs is crucial. Additionally, from a design aspect, loosening shape constraints and increasing the number of folds could introduce more variations in designs.

Given the significance of packaging as a practical application of origami, research should not be limited to the pillow box form but should encompass a broader range of shapes. This expanded focus could lead to the development of more innovative and versatile packaging solutions.

Acknowledgement

This work was supported by JST CREST Grant Number JPMJCR1911, Japan.

References

  • [Demaine et al. 10] Erik D. Demaine, Martin L. Demaine, and Duks Koschitz. “Reconstructing David Huffman’s Legacy in Curved-Crease Folding.” In Origami 5: Fifth International Meeting of Origami Science, Mathematics, and Education, edited by Patsy Wang-Iverson, Robert J. Lang, and Mark YIM, p. 39–52. A K Peters/CRC Press, 2010.
  • [Demaine et al. 15] Erik D. Demaine, Martin L. Demaine, David A. Huffman, Duks Koschitz, and Tomohiro Tachi. “Characterization of Curved Creases and Rulings: Design and Analysis of Lens Tessellations.” In Origami 6, p. 209–230, 2015.
  • [Fuchs and Tabachnikov 99] Dmitry Fuchs and Serge Tabachnikov. “More on Paperfolding.” The American Mathematical Monthly 106:1 (1999), 27–35.
  • [Kepert 98] Andrew Kepert. “Teabag Problem.”, 1998. Accessed via the Wayback Machine. Available online (https://web.archive.org/web/20050404002125/https://frey.newcastle.edu.au/~andrew/teabag/#ONE).
  • [Koiso 23] Miyuki Koiso. “Problem solving using differential geometry.” Mathematical Sciences 60:9 (2023), 14–20. (in Japanese).
  • [Mitani and Igarashi 11] Jun Mitani and Takeo Igarashi. “Interactive Design of Planar Curved Folding by Reflection.” In Pacific Conference on Computer Graphics and Applications, 2011.
  • [Tachi 11] Tomohiro Tachi. “One-DOF Rigid Foldable Structures from Space Curves.” In Proceedings of the IABSE-IASS Symposium 2011, 2011.
  • [web 97] “The Geometry Junkyard.”, 1997. Accessed via the Wayback Machine. Available online (https://web.archive.org/web/20050407162522/https://www.ics.uci.edu/~eppstein/junkyard/teabag.html).

 Jun Mitani
University of Tsukuba, 1-1-1 Tennohdai Tsukuba Ibaraki, Japan, e-mail: [email protected]