Next Article in Journal
Tissue Conditioner Incorporating a Nano-Sized Surface Pre-Reacted Glass-Ionomer (S-PRG) Filler
Next Article in Special Issue
Mechanical Behaviour of Aluminium-Timber Composite Connections with Screws and Toothed Plates
Previous Article in Journal
Method for High-Yield Hydrothermal Growth of Silica Shells on Nanoparticles
Previous Article in Special Issue
Investigation of the Effect of Pallet Top-Deck Stiffness on Corrugated Box Compression Strength as a Function of Multiple Unit Load Design Variables
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Efficient Design of Thin Wall Seating Made of a Single Piece of Heavy-Duty Corrugated Cardboard

Mechanical Engineering Department, Escuela Técnica Superior de Ingenieros Industriales, Universidad Politécnica de Madrid, C/José Gutiérrez Abascal, 2, 28006 Madrid, Spain
*
Author to whom correspondence should be addressed.
Materials 2021, 14(21), 6645; https://doi.org/10.3390/ma14216645
Submission received: 18 September 2021 / Revised: 27 October 2021 / Accepted: 30 October 2021 / Published: 4 November 2021
(This article belongs to the Special Issue Mechanics of Corrugated and Composite Materials)

Abstract

:
Corrugated cardboard has waved cores with small flutes that prevent the use of detailed numerical models of whole structures. Many homogenization methods in the literature overcome this drawback by defining equivalent homogeneous plates with the same mechanical behaviour at a macro-mechanical scale. However, few homogenization works have considered complete structures, focusing mainly on beams or plates. For the first time, this study explores the application of homogenization approaches to larger structures as an aid in their design process. We also considered triple-wall boards rather than single- and double-wall configurations commonly addressed in the literature. To this end, we adapted the homogenization methods proposed by Talbi and Duong to analyze thin-walled stools made of triple-wall corrugated cardboard. Using a progressive design process, we performed an efficient stool design by removing material zones with lower stresses, with 35% less material, 35% lower vertical deflections, and 66% lower stresses than the initial design. Unlike other corrugated cardboard stools, this design comprises just one folded piece instead of three, thus saving storage space. These results demonstrate the utility of homogenization techniques as an aid in the design process of whole structures made of corrugated cardboard. Further research will consider buckling analysis.

1. Introduction

Finite element analysis (FEA) greatly facilitates the design process of many products, avoiding the construction of failed prototypes. Concerning products made of corrugated cardboard, this advantage is not so evident since it is inexpensive and easy to handle, so that prototypes have low economic and time costs. In this paper, the authors aim to show that FEA can also be very useful when designing products made with this material. The main advantage is not to avoid prototyping, but to guide the design stages towards more efficient solutions. Likewise, it could help to choose the most suitable type of cardboard for each product, avoiding the need to gather an extensive assortment of materials to test different prototypes.
In this work, we applied FEA to a piece of furniture made of corrugated cardboard to achieve a more efficient design. To define the material properties, we adapted the homogenization methods proposed by Talbi [1] and Duong [2], as described in Section 2.3.
Conventional furniture designs often rely on traditional knowledge in handicraft manufacturing. Moreover, their structural elements are often intentionally oversized. However, FEA becomes an essential tool when dealing with unconventional furniture made of thin-wall structural elements. In [3,4,5,6,7,8,9], we can find some studies on the FEA of wood furniture. Other previous research studies also considered other materials, such as laminated bamboo [10], honeycomb cardboard [11], corrugated cardboard [12] or fibre-reinforced concrete [13].

1.1. Corrugated Cardboard

Corrugated cardboard is a material for everyday use, light, economical and sustainable. In addition to packaging, it can have other uses, such as construction and indoor furniture [14,15,16]. Its high strength-to-weight ratio makes it ideal for furniture manufacturing, though a careful design is needed to ensure rigidity.
It presents a sandwich structure with small waves in the intermediate layers (Figure 1), called fluting. Flutings are glued to flat sheets of paper, called liners, with a water-resistant starch-based adhesive [17]. Liners support bending loads, and flutings support transverse shear, helping to stabilize the former by resisting out-of-plane deformations [18,19]. In this way, the mechanical properties of liners and flutings are efficiently combined [20], providing a higher stiffness-to-weight ratio than an equivalent solid panel made of any of the individual constituent materials [21]. Liners are usually made of softwood kraft pulp to provide strength, with grammages ranging from 125 to 440 g/m2, while flutings have lower grammages, from 80 to 180 g/m2 [15,17,20]. Boards can present various wall configurations: single-sided, with only one fluting and one liner, and single-, double- and triple-wall (Figure 1), with the strength increasing with the number of plies.
Flutings are classified by their height and the number of flutes per unit length. Table 1 shows the most common flute types, designated as A, B, C, E or F, the C flute being the most commonly used for boxes. There are other less common flute types, such as D, with a height of 2 mm; G, thinner than 0.55 mm; K, thicker than 5.0 mm; and even a thinner flute, called O [15,22]. These letters were assigned according to their introduction into the market, having no relation to their size [17]. Larger flutes provide greater vertical strength and cushioning, while smaller flutes enhance graphic capabilities and structural integrity.
In paperboard manufacture, cellulose fibres tend to align in the flow direction, called machine direction (MD) [15]. The perpendicular direction on the paperboard surface is called cross direction (CD). Corrugated cardboard has the same manufacture direction that paperboard [24], MD being perpendicular to the principal axes of the corrugations and CD parallel to them (Figure 2). Then, both paper and corrugated cardboard are orthotropic materials, with better mechanical properties in MD than in CD [15,25].
Being a low-cost, lightweight, and environmentally friendly material, the use of corrugated cardboard for packaging has steadily increased in the past decade [19,26]. The global production of packaging paper and board increased from 193 to 256 million tons between 2008 and 2018 [27]. This effect was also influenced by the growth of online commerce [28]. In 2020, the global demand for containerboard was 69 million tons, 40% of the global demand for paper [29]. In 2018, the recycling rate for paper and cardboard packaging in the EU was 83% [30] of waste material. Waste cardboard can be used in its original form, but it can also be used in new composite materials [28].
Due to its great strength-to-weight ratio, excellent burst strength and resistance to crushing, corrugated cardboard is also suitable for furniture manufacture. However, a careful design is needed to ensure rigidity [31]. Thus, a good understanding of its mechanical behaviour is required to use it in an optimum way. Many previous studies have focused on the properties of corrugated cardboard and how the external environment affects its performance [20,32,33,34,35,36]. The mechanical properties of various types of liners and flutings in MD and CD can be found in [1,17,24,25,33,37,38,39,40,41,42,43,44,45,46,47,48,49].

1.2. Thin-Wall Furniture

Based on its thickness, we can classify the structural elements of furniture as ultrathin, below 10 mm; thin, from 10 to 15 mm; standard, from 16 to 19 mm; thick, from 20 to 40 mm; and ultra-thick, above 40 mm [50]. Thin-wall furniture, made of thin or ultrathin structural materials, is a current trend in furniture design [51]. It is usually made of wood composite panels, such as plywood, particleboard, or medium-density fibreboard (MDF), which can be laminated with other materials [50]. Due to its light weight, it can be considered a good alternative for trade shows and conventions. It can even be a suitable option for students or professionals with upward mobility, who will probably move often.
When dealing with the design of thin-wall furniture, a structural calculation is of particular relevance [50]. In addition to the strength requirements imposed on the materials, a second challenge lies in the joints between different panels [28,51]. Thin-walled structures can also exhibit buckling and warping problems, extensively studied in the scientific literature. Some analytical, numerical, and experimental studies on the buckling analysis of thin-wall beams can be found in [52,53]. Other studies on the buckling of corrugated cardboard structures can be found in [37,44,46,54,55].
In this work, our objective was to design a thin-wall furniture piece made of a different material, such as heavy-duty corrugated cardboard, whose sandwich structure could provide the required strength. Compared to wood composites, it has the advantage of being foldable. Thus, it requires fewer joints. Being low-cost and easy to transport and mount, in addition to the applications mentioned above, it can also be considered to meet the needs for accommodation in improvised shelters for emergencies [56].

Corrugated Cardboard Furniture

Corrugated cardboard furniture is usually made of pieces that could be flat-packed and assembled at home, using folds, slots and tabs. In the early 1960s, Peter Murdoch designed the Spotty chair [57], a flat-pack disposable chair that could be assembled simply by folding it in shape. In the early 1970s, Craig Hodgetts, Robert Mangurian, and Keith Godard designed Punch-Out [31], a low-cost furniture line made of heavy-duty corrugated cardboard, with flat pieces that even children could assemble to form their own tables and chairs. Today, many specialized companies [58,59,60,61,62,63,64,65,66,67,68,69] offer a great variety of corrugated cardboard furniture (such as chairs, armchairs, tables, shelves, beds, standing desks or podiums) [70], to be used at home, the office or trade shows. Many freelance designers also present their designs of corrugated cardboard furniture in design and architectural social media platforms or blogs [71,72,73].
As evidence of the growing interest in this type of furniture, the Japanese bedding company Airweave [74] provided 18,000 and 8,000 high-resistance cardboard beds for Olympians and Paralympians at the 2020 Tokyo Olympics [75,76]. They were conceived as a recycling initiative and were intended to be converted into other paper products. They will be reused for COVID-19 patients in a temporary medical facility in Osaka [77].
Another use of waste corrugated cardboard, as part of lightweight multi-layered panels with alternating plies of corrugated cardboard and veneer, was examined in [28]. Their study, considering different types of end corner joints between rigid panels, confirmed the suitability of this material for furniture and interior applications.

1.3. Homogenization Techniques

Different approaches can be used to analyze the strength of corrugated cardboard products: experimental [78]; analytical [79,80]; analytical-numerical [81,82,83] or purely numerical [33,84,85,86]. Due to the small size of the fluting, numerical methods are inadequate to analyze any structure made with this material on a micromechanical scale. Instead, we may use homogenization approaches. They allow considering its sandwich structure as a homogeneous plate [87,88], providing almost as accurate responses for homogenized models as for real structures [89].
Some homogenization techniques use analytic methods to obtain the engineering constants of the equivalent material [48,90,91,92,93]. Others apply the classical laminate theory (CLT) or the first-order shear deformation theory (FSDT) [1,2,45,94,95] to obtain the stiffness matrix of an equivalent plate [1,2,19,45,96,97,98,99]. Others use FEA of a representative volume element (RVE) to find an equivalent homogeneous plate [43,100,101,102,103,104].
Most homogenization studies centre on isolated flutings or single-wall corrugated boards, though some of them also consider double-walled corrugated panels [41,94,98,105,106,107]. Moreover, most of the existing literature on corrugated cardboard models focuses on homogenization methods, with few practical applications in actual designs.

1.4. Scope of the Study

This work aims to apply FEA for the structural calculation of corrugated cardboard furniture as an aid in its design process. As an example, we chose a stool made of this material to show the effectiveness of this method. This paper shows the process we followed to design the stool, performing a structural calculation of each intermediate design to assess its validity. In a future study, we also intend to consider a buckling analysis of the different design stages. However, this is beyond the scope of this work.

2. Materials and Methods

2.1. Design Stages

As a starting point, we based the first design on the geometry of a commercial stool, the so-called Kenno Stool [108,109] (Figure 3), designed by the Finnish designer Heikki Ruoho [110,111]. We chose this model for its simplicity. It comprises three pieces assembled perpendicularly, forming a closed structure that can be used as either a stool or a low table. It has a trapezoidal shape, resting on the ground, indistinctly, either on the wide or narrow part of the trapezoid. It has two vertical sidewalls with a vertical groove in the middle of their upper side. They are placed parallel to each other and covered by a third piece, whose ends fit into the groove of the former pieces.
We slightly increased its dimensions, since the original stool was conceived for children. We also replaced the original honeycomb cardboard with heavy-duty triple-wall corrugated cardboard [112,113,114,115], with which we obtained excellent results in a previous study of cardboard seating [12]. A 1970s child’s chair design from the hplusf design lab was made with this material. It was called Punch-out [116] and was temporarily exhibited at the MoMA [117]. Today, some contemporary furniture manufacturers, such as Chairigami (USA) [118] or Konno Konpou (Japan) [119], also use this material.
We applied FEA to this design, using a homogenization approach to characterize the mechanical properties of corrugated cardboard. In Section 2.3 and Appendix A, we present a thorough description of the homogenization technique used in this work.
From the numerical analysis performed, we obtained the deflections and stresses of this stool under some applied loads, according to the European Standards EN 1728 [120] and EN 12520 [121], both applicable to seating designs.
We then modified this design by removing both side panels. Therefore, the second design consisted of a single piece that the final user could fold for storage (Figure 4).
We lengthened the ends of the cardboard panel towards the opposite face and crossed them to ensure the structural strength of the stool. To maintain the total width of the top/bottom face, we placed the crossing point near it. We also reduced the width of one end to insert it into a slot made at the opposite end. The stool should also have two grooves on the top/bottom surface for inserting both ends, preventing them from moving. We also analyzed this design under the same load conditions.
To achieve more significant savings in material and storage space, we even opened the stool downward by removing the lower face (Figure 5). We now crossed both ends at an intermediate height inside the stool. However, it could rest only on the edges that limit the open surface, having a single possible position, unlike the previous designs.
Next, we modified the design by cutting both ends of the stool directly from the front and rear walls, opening a hole in those walls and folding the cut material inward (Figure 6a). This design saves even more material and storage space, since its ends could be placed inside the cut walls again. The angle formed between the front/rear wall and the seating surface should be the same as the angle between the ends and the seating surface, since both pieces should have the same length. During preliminary simulations, high longitudinal displacements were found at the bottom edges. Hence, we closed it on the bottom side by extending the front/rear walls to the bottom (Figure 6b) and connecting them.
Taking into account the orthotropic behaviour of corrugated cardboard, we analyzed each design for two material orientations: orientation I, with MD (x-axis) parallel to the folding lines, supposed to provide higher bending stiffness, and orientation II, with CD (y-axis) parallel to the folding lines, supposed to ease the folding process. We also considered two different body orientations: the wide part of the trapezium facing up and down.
The results thus obtained clearly show the utility of FEA, even for products made of an inexpensive and easy-to-handle material such as corrugated cardboard.

2.2. Finite Element Models

To develop the FE models of the stool designs, we used commercial software that includes a specific module for the structural analysis of composite materials.
We modelled the stool as a layered linear elastic shell. To do so, we combined the shell elements with a layered linear elastic material suitable for orthotropic laminates. In this way, the program applied the FSDT formulation internally. As input data, we introduced the stiffness matrices of the inner liners, the outer liners, and the fluting, together with the thickness and material model of each layer of the sandwich panel. We also used solid elements to model the loading pad used to apply loads on the seating surface.
We defined the contact conditions between intersecting panels using a mapped mesh defined so that two intersecting panels share the shell nodes lying on their intersection line. To define the contact between the solid elements of the loading pad and the shell elements of the stool panels, we used a multiphysics coupling provided by the commercial software; specifically, we used a solid-thin structure connection for this purpose.
Finally, we performed a static analysis with each model.
In the following sections, we define the FE model in more detail.

2.2.1. Geometry

All the designs considered had a seating surface 380 mm long and 400 mm wide and a height of 400 mm (Figure 7).
We built finite element (FE) models of all stool designs using homogenized shells with a mapped mesh (Figure 8) made up of square elements approximately 5 mm long. The number of boundary elements used in the models shown in Figure 8 ranges from 1704 for model (f) to 3810 for model (b).
The angle α between the top and front/rear panels was modified from 70° to 90°, with a 5° step, preserving the length of the seating surface. We considered the fourth design with α = 90° just for comparison, since it could rotate around the edges formed by the top and front/rear panels, thus being unstable. Figure 9 gathers the geometry variations for the fourth design to show where the board ends intersect the seating surface.

2.2.2. Material

The material considered for all designs was a heavy-duty triple-wall A-flute corrugated cardboard. Its homogenized properties were defined in the FE model using a layered material with seven layers: 1 and 7 are outer liners, 3 and 5 inner liners, and 2, 4, and 6 flutings (Figure 10). For each layer, we introduced either the liner thickness or the fluting height, together with its homogenized stiffness matrix, previously computed as described in Appendix A.
We used the engineering constants of the constituent materials reported in [45] to compute the stiffness matrices, since they have high elastic moduli and would provide high bending stiffness. Table 2 shows the engineering constants, Ei, Gij, and νij. They are given in the lamina reference frame, with the 2-axis parallel to the CD, and the 1-axis parallel to the MD.
To model the height and period of the fluting, we took the values indicated in the Tri-Wall Pak patent [114] for A flutes. We also took the thicknesses stated in [114] for the liners. For the fluting, we considered the grammage of 150 g/m2 specified in [113], corresponding to a thickness of 0.25 mm. Table 2 also shows the thickness of the liners and fluting, t, and the height, h, and period, P, of the fluting, all taken from the references mentioned above.
We also considered two orientations: I, with MD (red x-axis) parallel to the folding lines, and II, with CD (green y-axis) parallel to the folding lines (Figure 11).

2.2.3. Loads and Constraints

We applied the load distribution defined in the Eurocode EN 1728 [120], using a cylindrical loading pad, placed 175 mm from the front edge of the seat and centred on the width of the seating surface (Figure 12).
We modelled the pad as a solid steel cylinder with 180 mm diameter, covered with a 10 mm layer of polyurethane foam, using a free tetrahedral mesh. We applied a vertical force of 1300 N, according to Eurocode EN 12520 [121], for domestic seats. It was uniformly distributed on the upper surface of the cylinder and transmitted to the shell through a multiphysics coupling.
We applied simply supported boundary conditions at the lower edges of the folded panels (Figure 13). We restricted the three displacements of the lower front edge, but only the lateral, y, and vertical, z, displacements of the lower back edge (shown in blue).

2.3. Homogenization Approach

In this study, we applied a homogenization approach based on the first-order shear deformation theory (FSDT). It is an evolution of our previous work [12], which was in turn based on previous research by Talbi [1] and Duong [2].
The stiffness matrix of any lamina of a laminate can be easily formulated in the lamina reference frame, 123. However, to use a common reference system, we need to express the stiffness matrices of all laminas in the global laminate reference frame, xyz. This process is straightforward for liners, since they are flat, but not for flutings. Due to their waved shape, the material parameters for each section differ from the laminate reference frame, xyz, to the lamina reference frame, 123, in which they are known [87] (see Figure 14). Thus, we need to change the reference system of the stiffness matrix of the flutings.
Berthelot [55] applied a similar method to composite materials by rotating around the z-axis, normal to the laminate. For corrugated materials, however, the rotation has to be performed around the y-axis, or CD. Talbi [1] and Duong [2] performed this change of reference system to formulate their homogenization methods for single- and double-wall corrugated cardboard panels, respectively. Once the stiffness matrix of the fluting was transformed, they applied the FSDT to simplify the constitutive equations. Then, they integrated the stresses through the whole laminate thickness to get the internal forces, N and T, and the bending moments, M. After the integration, the z coordinate disappeared from the formulation, reducing the problem’s dimensionality from 3D to 2D. Then, they performed a second integration along the MD over a fluting period to obtain the average values. In this way, they expressed the generalized constitutive law as follows.
N M T = A B 0 B D 0 0 0 H · ε m κ γ s  
ε m is the membrane strain vector, κ the curvature vector, and γ s the transverse shear strain vector. A is the extensional stiffness matrix, D the bending stiffness matrix, B the bending-extension coupling stiffness matrix and H the transverse shear stiffness matrix. These matrices can be used to model a homogenized shell. For small structures, such as beams or plates, FE analysis can be performed analytically. However, when dealing with larger structures, an FE code is needed. Some FE packages include the FSDT formulation and directly work with the A, B, D, and H matrices. If it is not included, we can use the expressions found in the literature for the engineering constants of the homogenized shell as functions of these matrices [19,122].
In a previous work [12], we also applied this homogenization method. We computed the A, B, D, and H matrices outside the FE model and introduced them into the FE model. However, no additional information concerning the thickness and number of laminas was needed to perform the analysis. Since the FE model had no information to undo the homogenization after the simulation, the results of the analyses were averaged over the laminate thickness, and we needed to post-process them.
In this work, we used a different approach to avoid this post-processing, thus facilitating the graphical representation of the simulation results. As before, we changed the reference system to express the stiffness matrices of the corrugated layers in the laminate reference frame. Unlike before, this time, we directly introduced these matrices into the FE model. However, since they depend on the x-coordinate, they need to be processed before being introduced into the FE model. Thus, we performed a similar integration to that made by Talbi and Duong, but not on the A, B, D, and H matrices, but on the stiffness matrix of the corrugated layers. To do so, we first averaged each matrix through the z-coordinate and then over the x-direction, or MD (see Appendix A).
We then introduced the stiffness matrix of each layer into the FE model. We used a specialized module for composite materials that includes a layered linear elastic material model, which internally performs a second homogenization through the thickness of the whole laminate. It is based on the FSDT, like the methods of Talbi and Duong. This time, the total number of laminas and their respective thickness had to be introduced into the FE model. Then, it had the necessary information to undo the homogenization after the simulations. In this way, the results directly show different stress fields for each lamina, instead of just an average value, with no further post-processing.
The main drawback of this method is that it cannot be performed with basic FE packages but only with specific modules for composite materials. In return, we could simplify the calculation of the stiffness matrices while increasing the precision of the results of the FE analysis. Unlike before, any change in the number of sandwich layers or their thickness can be made directly inside the FE model, keeping the same stiffness matrices. Only when we want to change the geometry of the corrugated layers, we would need to recalculate their stiffness matrices outside the FE model. Using an FE module specialized in composite materials, this methodology also allows one to change the orientation of the corrugated panel and even to consider different orientations for individual layers inside the panel. If desired, it is also possible to perform delamination studies.

3. Results and Discussion

3.1. Homogenized Material Properties

Table 3 gathers the nonzero elements of the stiffness matrices computed for each layer of the corrugated board in the laminate reference frame, using Voigt notation. The fluting has lower values than the liners, since it is mainly void.

3.2. Parametric Study for α = 70° to 90°

For the four stool designs, we performed a parametric variation of α, from 70° to 90°, with a step of 5°. The influence of α found in the vertical deflections for the first, second, and third designs is very low. Similarly, its influence on the longitudinal deflections is also low for the first design. Figure 15 shows the longitudinal deflections for the second and third designs. For the second design, they decrease with increasing α, while for the third design, they increase with increasing α.
Figure 16 shows the vertical and longitudinal deflections found for the fourth design. The former have a minimum for α = 80° and the latter for α = 75°. However, since the vertical deflections are better for α = 80°, we consider this the best angle.
Figure 17 shows the stress distributions σxx and σyy in the global reference system for the fourth design. For both orientations, σxx and σyy also present a minimum for α = 80°.

3.3. Analysis of Designs with α = 80°

3.3.1. First Design for α = 80°

This paragraph shows the results obtained for the first design with α = 80° (see Figure 7), considered the best angle from the parametric analysis. Deflections u, v and w, are respectively aligned with the global x-, y- and z-axes (see Figure 13).
Figure 18 shows the vertical deflections, w, for designs with bottom and top discontinuities and both ply orientations.
For the designs with bottom discontinuity, the vertical deflections show a revolution geometry about the vertical axis, with a flat bottom. Their maximum values for top discontinuity are located on the seating surface panel closest to the load application area.
For the designs with bottom discontinuity, they are 11% lower for orientation II. However, for the designs with top discontinuity, they are 23% lower for orientation I, which provides a higher bending stiffness. They are lower for the designs with top discontinuity. They show an 82% reduction for orientation I from bottom to top discontinuity and a 74% reduction for orientation II. We can explain this reduction by the span length of the seating surface, which has a single panel for bottom discontinuity, but is divided into two panels with half the span length for top discontinuity.
For both discontinuities, the longitudinal deflections are lower for orientation II. In any case, the four designs analyzed show small values, below 1 mm.
Since the material is orthotropic, we should not use von Mises stresses. Figure 19 and Figure 20 respectively show the components σxx and σyy of the stress tensor in the laminate reference frame.
For the designs with bottom discontinuity, σxx and σyy are distributed mainly on the seating surface. For orientation I, σyy is also transmitted to the front and rear panels. On the contrary, for orientation II, σxx and σyy are transmitted to the side panels.
For the designs with top discontinuity, the maximum stresses were found on the panels covering the sidewalls, specifically at the vertical ends inserted into the side panels’ slots. Figure 21 shows the stress distribution for σxx and σyy for the designs with top discontinuity again, but now removing the front panel of the seating surface, thus revealing the stress distribution in such central panels, with the maximum stresses shown in dark red and dark blue.
According to the maximum stress criterion [123], applicable to orthotropic materials, the maximum values of σxx and σyy should be lower than the tensile strength of the constituent materials in the MD, σt,MD, and in the CD, σt,CD, respectively (see Figure 11).
The tensile strength of structural paper can vary from 17 to more than 75 MPa in MD and from 9 to 35 MPa in CD [124]. In this study, we considered as reference values the tensile strengths found in [124] for a base paper with similar elastic moduli that the constituent materials of the analyzed stool: σt,MD = 75.4 MPa in MD and σt,CD = 22.7 MPa in CD. For these limit values, all the configurations analyzed meet the maximum stress criterion. Moreover, even for other materials with tensile strengths quite close to the lower limit of the stress ranges indicated above, the stresses obtained would be above the limit values.

3.3.2. Second Design for α = 80°

This paragraph presents the results found for the second design with α = 80°. Figure 22 shows the vertical deflections for designs with bottom and top discontinuities and both ply orientations.
For orientation I, the distributions of vertical deflections show a geometry of revolution about the vertical axis, but they have an almost cylindrical shape for orientation II.
They are lower for orientation I and top discontinuity, showing a 92% reduction (from 14.44 to 1.07 mm). The improvement due to orientation is substantially more significant than for the first design, with 57% and 65% reductions for configurations with bottom and top discontinuities, respectively. Regarding the discontinuity location, for the top position, we found 82% and 78% improvements for orientations I and II, respectively.
Longitudinal deflections range from 0.3 to 1.4 mm. They are also lower for orientation I and top discontinuity. For both orientations, the highest values are in the middle-upper part of the front inner panel.
Figure 23 and Figure 24, respectively, show the stress distributions σxx and σyy.
For the bottom discontinuity, the maximum stresses concentrate in the central area of the seating panel. For orientation I, they are transmitted to the front and rear panels. However, for orientation II they are transmitted to the lateral edges of the seating surface.
Both σxx and σyy, are below σt,MD (75.4 MPa) and σt,CD (22.7 MPa), thus complying with the maximum stress criterion. Moreover, they would also be valid for any other structural paper, whose tensile strengths in MD and CD are, respectively, higher than 17 and 9 MPa.

3.3.3. Third Design for α = 80°

The stresses and vertical deflections are similar for the third design and the second design with bottom discontinuity. However, the longitudinal displacements are somewhat higher due to the removal of the lower panel. This effect was also shown in preliminary studies for the fourth design, with high longitudinal displacements of the lower rear edge (see Figure 5). Thus, we reintroduced the lower panel in the fourth design, since it prevents relative sliding between the front and rear lower edges.

3.3.4. Fourth Design for α = 80°

This paragraph presents the results found for the fourth design with α = 80°.
Figure 25 shows both the vertical and longitudinal deflections for ply orientation I.
The trend showing lower vertical deflections for orientation I than for orientation II is maintained. For orientation I, they show an almost trapezoidal shape. For orientation II, they show an almost cylindrical shape placed on the front side of the stool, with its axis oriented from side to side. Maximum vertical and longitudinal deflections for orientation I are 0.6 and 1.2 mm, respectively, these being quite low.
Figure 26 shows the stress distributions σxx and σyy in the global reference system.
Besides the seating surface, there are other higher stresses at the intersection of the inner panels with the seating surface and on the folding lines at the lower edge of the front and rear panels, which appear to act as stress concentrators. This effect can be seen as a consequence of using less material. However, the stresses in these zones are quite below the tensile stresses. So, they do not pose any problem, at least from a static point of view.
Both σxx and σyy are much lower than σt,MD (75.4 MPa) and σt,CD (22.7 MPa), thus fulfilling the maximum stress criterion. They would also be valid for any other structural paper, whose tensile strengths in MD and CD are, respectively, higher than 17 and 9 MPa.

3.3.5. Comparative Results for α = 80°

Figure 27 and Figure 28 show the vertical and longitudinal deflections for all designs.
The vertical deflections for orientation II are higher in the second than in the first design due to the elimination of the side panels. In contrast, for orientation I, they are of the same order of magnitude. In the worst case, the vertical deflections increase 206% (from 4.7 to 14.4 mm), from the first to the second design, while in the most favourable case, they increase by 22% (from 0.9 to 1.1 mm). We can explain this behaviour as the combination of the antagonistic effect of two factors: on the one hand, the negative effect of eliminating the side panels, especially in designs with orientation II, in which the stresses of the seating surface were transmitted to that panel, and on the other hand, the positive effect of introducing a new supporting system, with inner triangular structures. Although this second effect is beneficial for both orientations, it cannot overcome the negative effect of the other factor, especially in designs with orientation II.
For the best configurations of all the designs analyzed, the vertical deflections are lower for the fourth design, being 0.6 and 1.1 mm, respectively, for orientations I and II. Compared to the first design with top discontinuity, with vertical deflections of 0.9 and 1.2 mm, they are reduced by 33% and 8% for orientations I and II, respectively. Compared to the second design with top discontinuity, with deflections of 1.1 and 3.1 mm, they are reduced by 45% and 64%. Compared to the third design, with deflections of 6.0 and 13.8 mm, they are reduced by 90% and 92%. These reductions are lower with respect to the first and second designs because it was not possible to consider any configuration with top discontinuity in the third design.
It is remarkable that after removing a significant amount of material, the vertical deflections for the fourth design are even lower than for the first design with top discontinuity. This effect is due to the high efficiency of the inner panels added when removing the side panels from the first design, since they have a triangular structure, whose effectiveness is well known. Additionally, the results obtained for the fourth design are even better than those found for the second and third designs, which already included the inner panels. This behaviour is due to a better distribution of the intersection lines of the inner panels with the seating surface in the fourth design. The inner panels divide the seating surface into two halves, acting as intermediate supports in the second and third designs. However, they divide it into three zones in the fourth design, acting as two intermediate supports, thus reducing the span and, consequently, the maximum vertical deflections.
The longitudinal deflections are higher for the fourth than for the first and second designs. However, they are kept within reasonable limits of 1.2 mm for orientation I.
Figure 29 and Figure 30 respectively show the stress distributions σxx and σyy.
For orientation II, σxx and σyy are lower for the fourth design than for any other design. We can see the same trend for orientation I, except for the second design with top discontinuity.

3.4. Summary Results

For the best configurations of each design for α = 80°, Table 4 gathers the area, A, the vertical, w, and longitudinal, u, deflections and the stresses along MD, σxx, and CD, σyy.
Figure 31 shows the variation of these parameters compared to the best results of the first design, expressed as a percentage of the corresponding value in the first design.
In terms of deflections, the best results correspond to the fourth design, with orientation I and α = 80°, with maximum vertical and longitudinal deflections of 0.6 and 1.2 mm, respectively. There is a noticeable improvement compared to the first design, from which it evolved, whose best results are maximum vertical and longitudinal deflections of 0.9 and 0.2 mm, respectively. These results lead to a 33% reduction for the vertical deflections, but a 500% increase for the longitudinal deflections. Despite the high increase for the longitudinal deflection, its maximum value is just slightly above 1 mm for the fourth design. Moreover, the fourth design has an area 34% lower than the first design.
We can extract the following conclusions regarding different aspects of the possible configurations of the stool design:
  • Ply orientation. The vertical and longitudinal deflections are lower for orientation I, except for the first design, with slightly lower values for orientation I with bottom discontinuity.
  • Discontinuity location. The vertical and longitudinal deflections are lower for designs with top discontinuity, because the seating surface is divided into two different panels with half the span of the whole seating surface.
  • Bottom panel. We should keep the bottom panel, because it prevents longitudinal sliding between the lower edges of the front and rear panels.
  • α angle: In the first three designs, α has little influence on the vertical deflections. However, in the fourth design, the lowest vertical deflections correspond to an intermediate angle of 80°. The best results correspond to those angles leading to a more uniform distribution of the seating surface. That is, for those designs with the inner panels dividing the seating surface into three zones of equal length, so that none of them tends to present more significant deflections than the others (see Figure 8).

4. Conclusions

4.1. Main Findings

It is known that corrugated cardboard has higher bending stiffness for orientation I. In this work, we quantified this improvement for real applications. For vertical deflections, it ranges from 23% to 65%, finding with the best results for orientation I, except for the first design with bottom discontinuity, with slightly better results for orientation II.
For orientation I, the first design sidewalls show low stresses and can be removed.
As expected, the triangular structures inside the stool improve its static behaviour.
In the first and second designs with top discontinuity, the seating surface is divided into two parts with half the span of the seating surface. This division leads to a more favourable configuration than the corresponding designs with bottom discontinuity.
We should keep the bottom panel, since it prevents any longitudinal sliding between the lower edges of the front and rear panels.
The edges where the inner triangular structure contacts the seating panel act as intermediate supports. The seating surface has two intermediate supports in the fourth design, but only one in the second. Since configurations with more supports are most favourable, the fourth design has better static behaviour than the second and third designs.

4.2. Concluding Remarks

Corrugated cardboard has a great strength-to-weight ratio, excellent burst strength and resistance to crushing, thus being an ideal material for furniture manufacture. However, a careful design is needed to ensure rigidity [31]. This work aimed to apply numerical methods for the structural analysis of corrugated cardboard furniture, as an aid in their design process, to obtain efficient designs with the best resistance-to-cost ratio.
As an example, we chose a stool made of heavy-duty triple-wall A-flute corrugated cardboard. We performed static analyses on various stool designs, with a geometric evolution guided by the stresses found in previous design stages. The selection of this specific type of furniture has no particular relevance, being just a way to show the feasibility and benefits of numerical analysis in the design practice of corrugated cardboard furniture.
To define the mechanical properties of corrugated cardboard, we used a homogenization approach based on the first-order shear deformation theory (FSDT). It is an evolution of our previous work [12], which was in turn based on prior research by Talbi [1] and Duong [2]. Together with [12], a novelty of this work, is applying a homogenization technique to the numerical analysis of whole structures made of corrugated board, thus extending the scope of previous studies, usually limited to beams and plates. Although the analysis of simple structures, such as beams or plates, can be performed analytically, more complex structures, such as those considered in this work, should be studied with numerical techniques, such as FEA. A second novelty is the possibility of analyzing multiple-wall panels of any number of layers, in addition to the single- and double-wall configurations commonly addressed in the literature, also broadening the scope of previous works.
We computed the stiffness matrix of an equivalent homogeneous plate for each fluting, first averaging over the laminate thickness and then along the MD. To model the whole board, we inserted these stiffness matrices of flutings and liners into a FE model, using a layered material model based on the FSDT. Unlike other previous works, this methodology provides a way to easily model multiple-wall boards, since the homogenized matrices are independent of the number of plies of the laminate. In this way, the number of plies and the thickness of the liners can easily be changed inside the FE model.
We then performed a static analysis. The starting design of the stool evolved to three other designs, taking into account the deflections and stresses found in the FEA. Together, we analyzed four different designs under the load conditions defined by the Eurocodes EN 1728 [120] and EN 12520 [121] for seating. The first design, based on the geometry of a commercial stool made of three panels assembled in perpendicular directions, forming a closed structure, was chosen because of its simplicity. We found zones with lower stresses and progressively removed some of them. We also included an inner triangular structure to compensate for removing the side panels from the initial design.
The fourth design has higher strength than the others, showing the lowest vertical deflections and stresses, with reductions of 44% for w, 90% for σxx, and 67% for σyy compared to the starting design. It also requires 44% less material, thus reducing material costs. It is also made from a single foldable piece, requiring less storage space and reducing the possibility of losing pieces when stored. Therefore, it is significantly more efficient than the first design, based on its static behaviour, the amount of material needed and the required storage space. However, we do not discuss aesthetic or ergonomic aspects.
As expected, the results of this study demonstrate the utility of homogenization techniques as an aid in the design process of whole structures made of corrugated cardboard. The proposed methodology can be applied to the design process of any other piece of furniture, such as a shelf, a bed, a desk, or any other structural element made of corrugated cardboard. It can help to optimize its design by choosing an optimal geometry for a given material. It can also help to choose the most suitable material for a predefined geometry, by comparing panels with different numbers of walls. In both cases, it would lead to material savings. FE models can also be used to analyze delamination or buckling situations and take corrective actions when needed. These potential situations should be considered in future research. Comparative fatigue analyses would also be interesting [3].

Author Contributions

B.S.: conceptualization, methodology, software, formal analysis, investigation, resources, data curation, writing, writing—review & editing, visualization and supervision. L.M.M.: conceptualization, data curation, writing, writing–review & editing, funding acquisition; J.D.S.-B.: software, formal analysis, investigation, resources, writing. G.R.: methodology, investigation, resources, visualization, writing. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the raw/processed data required to reproduce these findings were presented in this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Homogenized Stiffness Matrix

As shown in (Figure 14), in woven laminas, the 2-axis is parallel to the y-axis. For the 1- and 3-axes, they are tangent and normal to the shape of the flute, respectively.
Flutings have a sinusoidal shape with height h, period P and thickness t (see Figure 14). Thus, in the xyz frame, the geometry of the k-th corrugated layer can be described as follows:
z = h k t k 2 s i n 2 π P k x .  
For the k-th corrugated layer, the angle θk between the local 1-axis, tangent to the flute, and the global x-axis is given by:
𝜃 𝑘 x = a t a n h k t k π P k c o s 2 π P k x .  
We can use the angle θk to transform the material properties of the constitutive layers from the lamina reference frame to the laminate reference frame. Thus, knowing the material parameters in the reference frame 123, they can be transformed into the reference frame xyz. For a rotation angle θk around the y-axis, the inverse strain transformation, from the reference frame xyz to the reference frame 123, can be written as [1,87] (see Appendix B.2):
ε 11 ε 22 ε 33 γ 23 γ 13 γ 12 = c 2 0 s 2 0 c s 0 0 1 0 0 0 0 s 2 0 c 2 0 c s 0 0 0 0 c 0 s 2 c s 0 2 c s 0 c 2 s 2 0 0 0 0 s 0 c · ε x x ε y y ε z z γ y z γ x z γ x y .
where:
c = c o s θ k ; s = s i n θ k .  
On the other hand, we can write the stress transformation (see Appendix B.3) as follows:
σ 11 σ 22 σ 33 σ 23 σ 13 σ 12 = c 2 0 s 2 0 2 c s 0 0 1 0 0 0 0 s 2 0 c 2 0 2 c s 0 0 0 0 c 0 s c s 0 c s 0 c 2 s 2 0 0 0 0 s 0 c · σ x x σ y y σ z z σ y z σ x z σ x y .
We can write Equations (A3) and (A5) in compact form as:
ε = T ε · ε  
σ = T σ · σ ,
ε and ε′ being the strain vectors in the global and local frames; σ and σ′ the stress vectors in the global and local frames; and Tε and Tσ the strain and stress transformation matrix, respectively.
The constitutive equations in the lamina and laminate reference frames are the following:
ε = S σ
ε = S · σ ,
S and S′ being the compliance matrices in the global and local reference frames, respectively.
Substituting Equations (A6) and (A7) in Equation (A8) and then comparing with Equation (A9), we get:
S = T ε 1 · S · T σ .  
Taking into account that (see Appendix B.3):
T ε 1 = T σ T  
we finally get:
S = T σ T · S · T σ .  
Equation (A12) transforms S′ from the lamina reference frame, where material properties are known, to the laminate reference frame, as a function of θ. For orthotropic materials, S′ is given by [55]:
S = 1 E 1 ν 12 E 1 ν 13 E 1 0 0 0 ν 12 E 1 1 E 2 ν 23 E 2 0 0 0 ν 13 E 1 ν 23 E 2 1 E 3 0 0 0 0 0 0 1 G 23 0 0 0 0 0 0 1 G 13 0 0 0 0 0 0 1 G 12 .
Using Equations (A12) and (A13), we get the following nonzero elements for S:
S 11 = c 4 E 1 + s 4 E 3 2 c 2 s 2 ν 13 E 1 + c 2 s 2 G 13 ;   S 12 = c 2 ν 12 E 1 s 2 ν 23 E 2 ; S 13 = c 4 ν 13 E 1 + c 2 s 2 E 3 + c 2 s 2 E 1 s 4 ν 13 E 1 c 2 s 2 G 13 ;   S 15 = 2 c 3 s ν 13 E 1 + 2 c s 3 E 3 2 c 3 s E 1 + 2 c s 3 ν 13 E 1 + c 3 s G 13 c s 3 G 13 ; S 22 = 1 E 2 ; S 23 = c 2 v 23 E 2 s 2 v 12 E 1 ; S 25 = 2 c s v 12 E 1 2 c s v 23 E 2 S 33 = c 4 E 3 + s 4 E 1 2 c 2 s 2 ν 13 E 1 + c 2 s 2 G 13 ;   S 35 = 2 c 3 s E 3 2 c s 3 ν 13 E 1 + 2 c 3 s ν 13 E 1 2 c s 3 E 1 c 3 s G 13 + c s 3 G 13 ; S 44 = c 2 G 23 + s 2 G 12 ;   S 46 = c s G 23 c s G 12 ; S 55 = c 4 G 13 + s 4 G 13 2 c 2 s 2 G 13 + 4 c 2 s 2 E 1 + 8 c 2 s 2 ν 13 E 1 + 4 c 2 s 2 E 3 ; S 66 = c 2 G 12 + s 2 G 23 .
After the transformation, we lose the typical structure of the stiffness matrix for orthotropic materials, with the elements of the fourth, fifth, and sixth rows and columns that lie outside its main diagonal equal to zero. The new structure is typical for monoclinic materials, with symmetry around y = 0 [125].
In the global reference frame, we can easily obtain the stiffness matrix, Q, by matrix inversion. Thence:
Q = S 1 .
Sij being functions of θ (through c and s), Q is also a function of θ, which is also a function of x.
For the liners, we can directly compute Qk from Equations (A13) and (A15). However, for the flutings, we need to average through the z-coordinate, taking into account that the thickness of a corrugated layer, obtained by a vertical cutting, is also a function of θ:
t z k x = t k c o s θ k x
tk being the uniform thickness of the fluting layer (see Figure 14). Thus, we can express the averaged stiffness through the layer thickness as:
Q ¯ k , z x = 1 h k 0 h k Q k x · d z = 1 h k Q k x 0 h k d z = 1 h k Q k x · t z k x = Q k x · t k h k · c o s θ k x .
Furthermore, to eliminate the dependence of the stiffness matrix on the x-coordinate, a second homogenization was performed along the x-direction, computing the average values of Qij over a fluting period [1,2,19,126]. This way, we got a constant value for each matrix element:
Q ¯ k , z x = 1 P 0 P Q ¯ k z x · d x = 1 P 0 P Q k x · t k h k · c o s θ k x · d x = t k P · h k 0 P Q k x c o s θ k x · d x .  

Appendix B. Transformation Matrices

Appendix B.1. Coordinate Transformation of a Generic Vector

The components of a vector r , in the laminate and lamina reference frames, can be related by:
r = A · r ,
where A is the transformation matrix:
A = c 0 s 0 1 0 s 0 c ,
being c = cosθ and s = sinθ.
The inverse transformation is:
r = A 1 · r = A T · r .  
From these relationships, the derivatives of x′, y′ and z′ with respect to x, y, and z are:
δ x δ x = c ;   δ x δ z = s ;   δ y δ y = 1 ;   δ z δ x = s ;   δ z δ z = c ;   δ x δ y = δ y δ x = δ y δ z = δ z δ y = 0 .
If, instead of a generic vector, we consider the displacement vector with components u, v, w:
u = c u + s w ;   v = v ;   w = s u + c w .
From these relationships, the derivatives of u, v, w with respect to x′, y′, and z′ are:
δ u δ x = c δ u δ x + s δ w δ x ;   δ u δ y = c δ u δ y + s δ w δ y ;   δ u δ z = c δ u δ z + s δ w δ z ; δ v δ x = δ v δ x ;   δ v δ y = δ v δ y ;   δ v δ z = δ v δ z ; δ w δ x = s δ u δ x + c δ w δ x ;   δ w δ y = s δ u δ y + c δ w δ y ;   δ w δ z = s δ u δ z + c δ w δ z .

Appendix B.2. Strain Transformations

By definition, the strains in the laminate reference frame, xyz, are given by:
ε x x = δ u δ x ;   ε y y = δ v δ y ;   ε z z = δ w δ z ;   γ y z = δ v δ z + δ w δ y ;   γ x z = δ u δ z + δ w δ x ;   γ x y = δ u δ y + δ v δ x .  
In the same way, strains in the lamina reference frame, 123, are given by:
ε 11 = δ u δ x ;   ε 22 = δ v δ y ;   ε 33 = δ w δ z ;   γ 23 = δ v δ z + δ w δ y ;   γ 13 = δ u δ z + δ w δ x ;   γ 12 = δ u δ y + δ v δ x .
Using the chain rule, together with Equation (A22), the strains in the laminate frame can be computed as follows:
ε x x = δ u δ x = δ u δ x δ x δ x + δ u δ y δ y δ x + δ u δ z δ z δ x = δ u δ x c + δ u δ y 0 + δ u δ z s = c δ u δ x + s δ u δ z .  
Replacing in these expressions the derivatives of Equation (A24), we get:
ε x x = c 2 ε 11 + c s γ 13 + s 2 ε 33 .
Analogously:
ε y y = ε 22 ;   ε z z = s 2 ε 11 c s γ 13 + c 2 ε 33 ; γ y z = c γ 23 s γ 12 ;   γ x z = c 2 s 2 γ 13 + 2 c s ε 33 ε 11 ;   γ x y = c γ 12 + s γ 23 .  
In matrix form, using the Voigt notation, which gathers the strains in a vector, we can write:
ε x x ε y y ε z z γ y z γ x z γ x y c 2 0 s 2 0 c s 0 0 1 0 0 0 0 s 2 0 c 2 0 c s 0 0 0 0 c 0 s 2 c s 0 2 c s 0 c 2 s 2 0 0 0 0 s 0 c · ε 11 ε 22 ε 33 γ 23 γ 13 γ 12 .
We obtained the inverse transformation by replacing θ by –θ, that is, sinθ by -sinθ:
ε 11 ε 22 ε 33 γ 23 γ 13 γ 12 = c 2 0 s 2 0 c s 0 0 1 0 0 0 0 s 2 0 c 2 0 c s 0 0 0 0 c 0 s 2 c s 0 2 c s 0 c 2 s 2 0 0 0 0 s 0 c · ε x x ε y y ε z z γ y z γ x z γ x y .
This transformation and its inverse transformation can be written in compact form as follows:
ε = T ε · ε ;   ε = T ε 1 · ε ,  
being Tε the strain transformation matrix [1,87]:
T ε = c 2 0 s 2 0 c s 0 0 1 0 0 0 0 s 2 0 c 2 0 c s 0 0 0 0 c 0 s 2 c s 0 2 c s 0 c 2 s 2 0 0 0 0 s 0 c ;   T ε 1 = c 2 0 s 2 0 c s 0 0 1 0 0 0 0 s 2 0 c 2 0 c s 0 0 0 0 c 0 s 2 c s 0 2 c s 0 c 2 s 2 0 0 0 0 s 0 c .

Appendix B.3. Stresses Transformations

Cauchy’s law [123] gives the components t1, t2, t3 of the stress vector as functions of the stress tensor and the components n1, n2, n3 of the normal vector to the surface. In the lamina reference frame, 123:
t 1 t 2 t 3 = σ 11 σ 12 σ 13 σ 12 σ 22 σ 23 σ 13 σ 23 σ 33 · n 1 n 2 n 3 .
In the laminate reference frame, xyz, it can also be written as:
t x t y t z = σ x x σ x y σ x z σ x y σ y y σ y z σ x z σ y z σ z z · n x n y n z .
Or, in compact form:
t = σ · n ;   t = σ · n .
Using Equation (A21) to transform from the local to the global reference frames:
A · t = σ · A · n   t = A 1 · σ · A · n = A T · σ · A · n .
Identifying with Equation (A36), we have:
σ = A T · σ · A .
Thence, from Equation (A20):
σ = σ x x σ x y σ x z σ x y σ y y σ y z σ x z σ y z σ z z = c 0 s 0 1 0 s 0 c σ 11 σ 12 σ 13 σ 12 σ 22 σ 23 σ 13 σ 23 σ 33 c 0 s 0 1 0 s 0 c = c 2 σ 11 + 2 c s σ 13 + s 2 σ 33 c σ 12 + s σ 23 c 2 s 2 σ 13 + c s σ 33 σ 11 c σ 12 + s σ 23 σ 22 c σ 23 s σ 12 c 2 s 2 σ 13 + c s σ 33 σ 11 c σ 23 s σ 12 s 2 σ 11 + c 2 σ 33 2 c s σ 13 .
Using the Voigt notation again, we can rewrite the above expressions as:
σ x x σ y y σ z z σ y z σ x z σ x y = c 2 0 s 2 0 2 c s 0 0 1 0 0 0 0 s 2 0 c 2 0 2 c s 0 0 0 0 c 0 s c s 0 c s 0 c 2 s 2 0 0 0 0 s 0 c · σ 11 σ 22 σ 33 σ 23 σ 13 σ 12 .
The inverse transformation is obtained by replacing θ by −θ, that is, substituting sinθ by −sinθ:
σ 11 σ 22 σ 33 σ 23 σ 13 σ 12 = c 2 0 s 2 0 2 c s 0 0 1 0 0 0 0 s 2 0 c 2 0 2 c s 0 0 0 0 c 0 s c s 0 c s 0 c 2 s 2 0 0 0 0 s 0 c · σ x x σ y y σ z z σ y z σ x z σ x y .
We can write in compact form this transformation and its inverse transformation as follows:
σ = T σ · σ ;   σ = T σ 1 · σ .  
being Tσ the stress transformation matrix:
T σ = c 2 0 s 2 0 2 c s 0 0 1 0 0 0 0 s 2 0 c 2 0 2 c s 0 0 0 0 c 0 s c s 0 c s 0 c 2 s 2 0 0 0 0 s 0 c ;   T σ 1 = c 2 0 s 2 0 2 c s 0 0 1 0 0 0 0 s 2 0 c 2 0 2 c s 0 0 0 0 c 0 s c s 0 c s 0 c 2 s 2 0 0 0 0 s 0 c .
As can be seen, the stiffness and strain transformation matrices satisfy the following relations:
T ε 1 = T σ T  
T σ 1 = T ε T .  

References

  1. Talbi, N.; Batti, A.; Ayad, R.; Guo, Y. An analytical homogenization model for finite element modelling of corrugated cardboard. Compos. Struct. 2009, 88, 280–289. [Google Scholar] [CrossRef]
  2. Duong, P.T.M. Modeling and Numerical Simulation for the Double Corrugated Cardboard under Transverse Loading by Homogenization Method. Int. J. Eng. Sci. 2017, 6, 16–25. [Google Scholar]
  3. Ceylan, E.; Güray, E.; Kasal, A. Structural analyses of wooden chairs by finite element method (FEM) and assessment of the cyclic loading performance in comparison with allowable design loads. Maderas-Cienc. Tecnol. 2021, 23. [Google Scholar] [CrossRef]
  4. Smardzewski, J.; Gawronski, T. FEM Algorithm for chair optimisation. Electron. J. Pol. Agric. Univ. Ser. Wood Technol. 2001, 4, 2. [Google Scholar]
  5. Kasal, A. Determination of the strength of various sofa frames with finite element analysis. Gazi Univ. J. Sci. 2006, 19, 191–203. [Google Scholar]
  6. Hu, W.; Liu, N.; Guan, H. Optimal design of a furniture frame by reducing the volume of wood. Drewno 2019, 62, 85–97. [Google Scholar]
  7. Gustafsson, S.-I. Optimising ash wood chairs. Wood Sci. Technol. 1997, 31, 291–301. [Google Scholar] [CrossRef]
  8. Smardzewski, J. Numerical analysis of furniture constructions. Wood Sci. Technol. 1998, 32, 273–286. [Google Scholar] [CrossRef]
  9. Tankut, N.; Tankut, A.N.; Zor, M. Finite Element Analysis of Wood Materials. Drv. Ind. 2014, 65, 159–171. [Google Scholar] [CrossRef]
  10. Laemlaksakul, V. Innovative design of laminated bamboo furniture using finite element method. Int. J. Math. Comput. Simul. 2008, 2, 274–284. [Google Scholar]
  11. Smardzewski, J.; Kramski, D. Modelling stiffness of furniture manufactured from honeycomb panels depending on changing climate conditions. Thin-Walled Struct. 2019, 137, 295–302. [Google Scholar] [CrossRef]
  12. Suarez, B.; Muneta, M.L.M.; Sanz-Bobi, J.D.; Romero, G. Application of homogenization approaches to the numerical analysis of seating made of multi-wall corrugated cardboard. Compos. Struct. 2021, 262, 113642. [Google Scholar] [CrossRef]
  13. Baioni, E.; Alessi, R.; Corinaldesi, V.; Lancioni, G.; Rizzini, R. Feasibility Study of a Table Prototype Made of High-Performance Fiber-Reinforced Concrete. Technologies 2017, 5, 41. [Google Scholar] [CrossRef] [Green Version]
  14. Abdel-Mohsen, M.; Faggal, A.; El-Metwally, Y. Efficiency of Corrugated Cardboard as a Building Material. 2012. Available online: https://www.researchgate.net/publication/351323568 (accessed on 15 September 2021). [CrossRef]
  15. Latka, J.F. Paper in architecture: Research by design, engineering and prototyping. A+ BE Archit. Built Environ. 2017, 19, 1–532. [Google Scholar]
  16. El Damatty, A.; Mikhail, A.; Awad, A. Finite element modeling and analysis of a cardboard shelter. Thin-Walled Struct. 2000, 38, 145–165. [Google Scholar] [CrossRef]
  17. Zhang, Z.; Qiu, T.; Song, R.; Sun, Y. Nonlinear Finite Element Analysis of the Fluted Corrugated Sheet in the Corrugated Cardboard. Adv. Mater. Sci. Eng. 2014, 2014, 1–8. [Google Scholar] [CrossRef] [Green Version]
  18. Barbero, E.J. Finite Element Analysis of Composite Materials Using Abaqus TM; CRC-Press: Boca Raton, FL, USA, 2013. [Google Scholar]
  19. Aboura, Z.; Talbi, N.; Allaoui, S.; Benzeggagh, M. Elastic behavior of corrugated cardboard: Experiments and modeling. Compos. Struct. 2004, 63, 53–62. [Google Scholar] [CrossRef] [Green Version]
  20. Twede, D.; Selke, S.; Kamdem, D.; Shires, D. Cartons, Crates and Corrugated Board: Handbook of Paper and Wood Packaging Technology; DEStech Publications: Lancaster, PA, USA, 2014. [Google Scholar]
  21. Dayyani, I.; Shaw, A.D.; Flores, E.I.S.; Friswell, M.I. The Mechanics of Composite Corrugated Structures: A Review with Applications in Morphing Aircraft; Elsevier: Amsterdam, The Netherlands, 2015; pp. 358–380. [Google Scholar]
  22. Kirwan, M.J. Handbook of Paper and Paperboard Packaging Technology, 2nd ed.; John Wiley and Sons: Hoboken, NJ, USA, 2012. [Google Scholar]
  23. Jay Enterprises Limited. Cardboard Thickness and Material. Available online: https://www.jay.co.jp/page/Cardboard_thickness_and_material (accessed on 21 October 2021).
  24. Allaoui, S.; Aboura, Z.; Benzeggagh, M.L. Phenomena governing uniaxial tensile behaviour of paper-board and corrugated cardboard. Compos. Struct. 2009, 87, 80–92. [Google Scholar] [CrossRef] [Green Version]
  25. Luong, V.D.; Abbès, F.; Abbès, B.; Duong, P.T.M.; Nolot, J.-B.; Erre, D.; Guo, Y.-Q. Finite Element Simulation of the Strength of Corrugated Board Boxes Under Impact Dynamics. In Proceedings of the International Conference on Advances in Computational Mechanics, Phu Quoc, Vietnam, 2–4 August 2017; pp. 369–380. [Google Scholar] [CrossRef]
  26. Statista. Production Volume of Paper and Cardboard Worldwide from 2008 to 2018. 2020. Available online: https://www.statista.com/statistics/270314/global-paper-and-cardboard-production/ (accessed on 21 October 2021).
  27. Statista. Global Paper Production Volume from 2008 to 2018 by Type. 2020. Available online: https://www.statista.com/statistics/270317/production-volume-of-paper-by-type/ (accessed on 21 October 2021).
  28. Jivkov, V.; Simeonova, R.; Antov, P.; Marinova, A.; Petrova, B.; Kristak, L. Structural Application of Light-Weight Panels Made of Waste Cardboard and Beech Veneer. Materials 2021, 14, 5064. [Google Scholar] [CrossRef] [PubMed]
  29. Statista. Paper Demand Worldwide in 2020, by Type. 2021. Available online: https://www.statista.com/statistics/1089092/global-paper-consumption-by-type/ (accessed on 21 October 2021).
  30. Eurostat. Recycling Rate of Packaging Waste by Type of Packaging. 2021. Available online: https://ec.europa.eu/eurostat/databrowser/view/cei_wm020/default/bar?lang=en%2F (accessed on 21 October 2021).
  31. Vigna, F. Materiality. Cardboard. Kinder. A Journal Dedicated to Child Design, Past, Present and Future. Available online: http://thekinderjournal.com/003-materiality (accessed on 15 September 2021).
  32. Lu, T.; Chen, C.; Zhu, G. Compressive behaviour of corrugated board panels. J. Compos. Mater. 2001, 35, 2098–2126. [Google Scholar] [CrossRef]
  33. Nordstrand, T. Analysis and testing of corrugated board panels into the post-buckling regime. Compos. Struct. 2004, 63, 189–199. [Google Scholar] [CrossRef]
  34. Nyman, U.; Gustafsson, P.J. Material and structural failure criterion of corrugated board facings. Compos. Struct. 2000, 50, 79–83. [Google Scholar] [CrossRef]
  35. Daum, M.; Darby, D.; Batt, G.; Campbell, L. Application of the stress-energy method for generating corrugated board cushion curves. J. Test. Eval. 2013, 41, 590–601. [Google Scholar] [CrossRef]
  36. Viguié, J.; Dumont, P.; Orgéas, L.; Vacher, P.; Desloges, I.; Mauret, E. Surface stress and strain fields on compressed panels of corrugated board boxes: An experimental analysis by using digital image stereocorrelation. Compos. Struct. 2011, 93, 2861–2873. [Google Scholar] [CrossRef]
  37. Fadiji, T.; Ambaw, A.; Coetzee, C.J.; Berry, T.M.; Opara, U.L. Application of finite element analysis to predict the mechanical strength of ventilated corrugated paperboard packaging for handling fresh produce. Biosyst. Eng. 2018, 174, 260–281. [Google Scholar] [CrossRef]
  38. Allaoui, S.; Aboura, Z.; Benzeggagh, M.L. Effects of the environmental conditions on the mechanical behaviour of the corrugated cardboard. Compos. Sci. Technol. 2009, 69, 104–110. [Google Scholar] [CrossRef]
  39. Patel, P.; Nordstrand, T.; Carlsson, L.A. Local buckling and collapse of corrugated board under biaxial stress. Compos. Struct. 1997, 39, 93–110. [Google Scholar] [CrossRef]
  40. Allansson, A.; Svärd, B. Stability and Collapse of Corrugated Board: Numerical and Experimental Analysis. Master’s Thesis, Lund University, Lund, Sweden, 2001. [Google Scholar]
  41. Park, J.; Kim, G.; Kwon, S.; Chung, S.; Kwon, S.; Choi, W.; Mitsuoka, M.; Inoue, E.; Okayasu, T.; Choe, J. Finite Element Analysis of Corrugated Board under Bending Stress. J. Fac. Agric. Kyushu Univ. 2012, 57, 181–188. [Google Scholar] [CrossRef]
  42. Fadiji, T.; Berry, T.; Coetzee, C.; Opara, L. Investigating the Mechanical Properties of Paperboard Packaging Material for Handling Fresh Produce Under Different Environmental Conditions: Experimental Analysis and Finite Element Modelling. 2017. Available online: https://scholarworks.rit.edu/japr/vol9/iss2/3 (accessed on 21 October 2021).
  43. Daxner, T. Optimum Design of Corrugated Board under Buckling Constraints. In Proceedings of the 7th World Congress on Structural and Multidisciplinary Optimization, Seoul, Korea, 21–25 May 2007. [Google Scholar]
  44. Biancolini, M.; Brutti, C.; Porziani, S. Corrugated board containers design methods. Int. J. Comput. Mater. Sci. Surf. Eng. 2010, 3, 143–163. [Google Scholar] [CrossRef]
  45. Minh, D.P.T. Analysis and simulation for the double corrugated cardboard plates under bending and in-plane shear force by homogenization method. Int. J. Mech. 2017, 11, 176–181. [Google Scholar]
  46. Fadiji, T.; Coetzee, C.; Opara, U.L. Compression strength of ventilated corrugated paperboard packages: Numerical modelling, experimental validation and effects of vent geometric design. Biosyst. Eng. 2016, 151, 231–247. [Google Scholar] [CrossRef]
  47. Isaksson, P.; Krusper, A.; Gradin, P. Shear correction factors for corrugated core structures. Compos. Struct. 2007, 80, 123–130. [Google Scholar] [CrossRef]
  48. Nordstrand, T.; Carlsson, L.A.; Allen, H.G. Transverse shear stiffness of structural core sandwich. Compos. Struct. 1994, 27, 317–329. [Google Scholar] [CrossRef]
  49. Kalyankar, N.; Mahakalkar, S.; Giri, J. A Review on Optimization of Corrugated Sheet Box Size for an Industrial Part. Int. J. Res. Advent Technol. ICATEST 2015, 16–19. [Google Scholar]
  50. Jivkov, V.; Elenska-Valchanova, D. Mechanical Properties of Some Thin Furniture Structural Composite Materials. In Proceedings of the 30th International Conference on Wood Science and Technology, Zagreb, Croatia, 12–13 December 2019. [Google Scholar]
  51. Jivkov, V.; Petrova, B. Challenges for furniture design with thin structural materials. In Proceedings of the VI International Furniture Congress, Trabzon, Turkey, 2–4 November 2020. [Google Scholar]
  52. Ruta, G.; Pignataro, M.; Rizzi, N. A direct one-dimensional beam model for the flexural-torsional buckling of thin-walled beams. J. Mech. Mater. Struct. 2006, 1, 1479–1496. [Google Scholar] [CrossRef] [Green Version]
  53. Piana, G.; Lofrano, E.; Carpinteri, A.; Ruta, G. Effect of local stiffeners and warping constraints on the buckling of symmetric open thin-walled beams with high warping stiffness. Meccanica 2021, 56, 2083–2102. [Google Scholar] [CrossRef]
  54. Jiménez-Caballero, M.A.; Conde, I.; García, B.; Liarte, E. Design of different types of corrugated board packages using finite element tools. In Proceedings of the SIMULIA Customer Conference, London, UK, 19–21 May 2009. [Google Scholar]
  55. Berthelot, J. Composite Materials-Mechanical Behavior and Structural Analysis; Springer: Berlin/Heidelberg, Germany, 1999. [Google Scholar]
  56. Spinillo, C.G.; Fujita, P.T. Do-it-yourself (DIY) furniture for emergency situations: A study on assembling a cardboard bench in Brazil. Theor. Issues Ergon. Sci. 2012, 13, 121–134. [Google Scholar] [CrossRef]
  57. Cartonlab, La Silla de Cartón. Una Historia de Diseño Sostenible. Available online: https://cartonlab.com/blog/silla-carton/ (accessed on 21 October 2021).
  58. Chairigami. Available online: https://www.chairigami.com/ (accessed on 21 October 2021).
  59. Quart De Poil. Cardboard. Available online: https://quartdepoil.fr/en/16,cardboard.html (accessed on 21 October 2021).
  60. Cartonlab. Available online: https://cartonlab.com/en/ (accessed on 21 October 2021).
  61. Kartelier. Available online: https://www.facebook.com/kartelier/about/ (accessed on 21 October 2021).
  62. Kartent. Available online: https://shop.kartent.com/en/ (accessed on 21 October 2021).
  63. Stange. Available online: https://www.pappmoebelshop.de/home.html (accessed on 21 October 2021).
  64. Origami Furniture. Collections. Available online: http://www.origamifurniture.com/collections/?lang=en (accessed on 21 October 2021).
  65. Litencarton. Available online: https://www.litencarton.ch/ (accessed on 21 October 2021).
  66. Cardboard. Cardboard Furniture and Projects. Available online: http://cardboard.es/portfolio/#productos-de-carton (accessed on 21 October 2021).
  67. Cartone Design. Available online: https://cartonedesign.com.br/produtos/ (accessed on 21 October 2021).
  68. Igreen. Cool Green Gadgets. Available online: https://www.igreengadgets.com/ (accessed on 21 October 2021).
  69. Danbaul x Style. Available online: https://www.danbaul.com/ (accessed on 21 October 2021).
  70. Architecture Art Designs. 30 Amazing Cardboard DIY Furniture Ideas. Available online: https://www.architectureartdesigns.com/30-amazing-cardboard-diy-furniture-ideas/ (accessed on 21 October 2021).
  71. Behance. Available online: https://www.behance.net/ (accessed on 21 October 2021).
  72. Designboom. Available online: https://www.designboom.com/ (accessed on 21 October 2021).
  73. Inhabitat. Available online: https://inhabitat.com/ (accessed on 21 October 2021).
  74. Airweave. Available online: https://www.airweave.com/ (accessed on 21 October 2021).
  75. Gleeson, S. Athletes to Sleep on Recyclable ‘Cardboard’ Beds during 2020 Olympics in Tokyo. USA Today Sports. 2019. Available online: https://eu.usatoday.com/story/sports/olympics/2019/09/24/olympics-2020-athletes-sleep-recyclable-cardboard-beds-tokyo/2429695001/ (accessed on 21 October 2021).
  76. Finney, A. Airweave Creates Cardboard Beds and Modular Mattresses for Tokyo 2020 Olympics. 2021. Available online: https://www.dezeen.com/2021/07/11/cardboard-beds-modular-mattresses-airweave-tokyo-2020-olympics/ (accessed on 21 October 2021).
  77. Lloyd, O. Tokyo 2020 Cardboard Beds to Be Made Available for COVID-19 Patients in Osaka Inside the Games. 2021. Available online: https://www.insidethegames.biz/articles/1112915/tokyo-2020-cardboard-beds-covid19-osaka (accessed on 21 October 2021).
  78. Dimitrov, K.; Heydenrych, M. Relationship between the edgewise compression strength of corrugated board and the compression strength of liner and fluting medium papers. South. For. J. For. Sci. 2009, 71, 227–233. [Google Scholar] [CrossRef]
  79. Magnucka-Blandzi, E.; Magnucki, K.; Wittenbeck, L. Mathematical modeling of shearing effect for sandwich beams with sinusoidal corrugated cores. Appl. Math. Model. 2015, 39, 2796–2808. [Google Scholar] [CrossRef]
  80. Magnucka-Blandzi, E.; Magnucki, K. Transverse shear modulus of elasticity for thin-walled corrugated cores of sandwich beams. Theoretical study. J. Theor. Appl. Mech. 2014, 52, 971–980. [Google Scholar] [CrossRef] [Green Version]
  81. Garbowski, T.; Gajewski, T.; Grabski, J.K. The Role of Buckling in the Estimation of Compressive Strength of Corrugated Cardboard Boxes. Materials 2020, 13, 4578. [Google Scholar] [CrossRef]
  82. Garbowski, T.; Gajewski, T.; Grabski, J.K. Estimation of the Compressive Strength of Corrugated Card-board Boxes with Various Openings. Energies 2020, 14, 155. [Google Scholar] [CrossRef]
  83. Garbowski, T.; Gajewski, T.; Grabski, J.K. Estimation of the Compressive Strength of Corrugated Card-board Boxes with Various Perforations. Energies 2021, 14, 1095. [Google Scholar] [CrossRef]
  84. Garbowski, T.; Gajewski, T.; Grabski, J.K. Role of Transverse Shear Modulus in the Performance of Corrugated Materials. Materials 2020, 13, 3791. [Google Scholar] [CrossRef] [PubMed]
  85. Nordstrand, T.; Carlsson, L. Evaluation of transverse shear stiffness of structural core sandwich plates. Compos. Struct. 1997, 37, 145–153. [Google Scholar] [CrossRef]
  86. Luong, V.D.; Bonnin, A.-S.; Abbès, F.; Nolot, J.-B.; Erre, D.; Abbès, B. Finite Element and Experimental Investigation on the Effect of Repetitive Shock in Corrugated Cardboard Packaging. J. Appl. Comput. Mech. 2021, 7, 820–830. [Google Scholar]
  87. Marek, A.; Garbowski, T. Homogenization of sandwich panels. Comput. Assist. Methods Eng. Sci. 2015, 22, 39–50. [Google Scholar]
  88. Fadiji, T.; Coetzee, C.; Berry, T.M.; Ambaw, A.; Opara, U.L. The efficacy of finite element analysis (FEA) as a design tool for food packaging: A review. Biosyst. Eng. 2018, 174, 20–40. [Google Scholar] [CrossRef]
  89. Buannic, N.; Cartraud, P.; Quesnel, T. Homogenization of corrugated core sandwich panels. Compos. Struct. 2003, 59, 299–312. [Google Scholar] [CrossRef] [Green Version]
  90. Bartolozzi, G.; Orrenius, U.; Pratellesi, A.; Pierini, M. An Equivalent Orthotropic Plate Model for Sinusoidal Core Sandwich Panels in Optimization Processes. In Proceedings of the INTER-NOISE and NOISE-CON Congress and Conference Proceedings, Sorrento, Italy, 16 November 2012; pp. 169–180. [Google Scholar]
  91. Bartolozzi, G.; Pierini, M.; Orrenius, U.; Baldanzini, N. An equivalent material formulation for sinusoidal corrugated cores of structural sandwich panels. Compos. Struct. 2013, 100, 173–185. [Google Scholar] [CrossRef]
  92. Bartolozzi, G.; Baldanzini, N.; Pierini, M. Equivalent properties for corrugated cores of sandwich structures: A general analytical method. Compos. Struct. 2014, 108, 736–746. [Google Scholar] [CrossRef]
  93. Shaban, M.; Alibeigloo, A. Three-dimensional elasticity solution for sandwich panels with corrugated cores by using energy method. Thin-Walled Struct. 2017, 119, 404–411. [Google Scholar] [CrossRef]
  94. Minh, D.P.T.; Khoa, N.N. An analytic homogenization model in traction and bending for orthotropic composite plates with the type of double corrugated cardboard. Vietnam. J. Mech. 2016, 38, 205–213. [Google Scholar] [CrossRef] [Green Version]
  95. Nguyen-Minh, N.; Tran-Van, N.; Bui-Xuan, T.; Nguyen-Thoi, T. Free vibration analysis of corrugated panels using homogenization methods and a cell-based smoothed Mindlin plate element (CS-MIN3). Thin-Walled Struct. 2018, 124, 184–201. [Google Scholar] [CrossRef]
  96. Ishikawa, T.; Chou, T.-W. Stiffness and strength behaviour of woven fabric composites. J. Mater. Sci. 1982, 17, 3211–3220. [Google Scholar] [CrossRef]
  97. Biancolini, M. Evaluation of equivalent stiffness properties of corrugated board. Compos. Struct. 2005, 69, 322–328. [Google Scholar] [CrossRef]
  98. Duong, P.; Abbès, B.; Djilali, H.A.; Hammou, A.; Makhlouf, M.; Guo, Y. An analytic homogenisation model for shear–torsion coupling problems of double corrugated core sandwich plates. J. Compos. Mater. 2012, 47, 1327–1341. [Google Scholar] [CrossRef]
  99. Garbowski, T.; Gajewski, T.; Grabski, J.K. Torsional and transversal stiffness of orthotropic sandwich panels. Materials 2020, 13, 5016. [Google Scholar] [CrossRef]
  100. Sharma, A.; Sankar, B.V.; Haftka, R.T. Homogenization of Plates with Microstructure and Application to Corrugated Core Sandwich Panels. In Proceedings of the 51st AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics and Materials Conference, Orlando, FL, USA, 12–15 April 2010. [Google Scholar]
  101. Atashipour, S.R.; Al-Emrani, M. A realistic model for transverse shear stiffness prediction of composite corrugated-core sandwich elements. Int. J. Solids Struct. 2017, 129, 1–17. [Google Scholar] [CrossRef]
  102. Heller, D. A Nonlinear Multiscale Finite Element Model for Comb-Like Sandwich Panels. Ph.D. Thesis, Technische Universität Darmstadt, Darmstadt, Germany, 2015. [Google Scholar]
  103. Montemurro, M.; Catapano, A.; Doroszewski, D. A multi-scale approach for the simultaneous shape and material optimisation of sandwich panels with cellular core. Compos. Part B Eng. 2016, 91, 458–472. [Google Scholar] [CrossRef] [Green Version]
  104. Smardzewski, J.; Jasińska, D. Mathematical models and experimental data for HDF based sandwich panels with dual corrugated lightweight core. Holzforschung 2017, 71, 265–273. [Google Scholar] [CrossRef]
  105. Hernández-Pérez, A.; Hägglund, R.; Carlsson, L.A.; Avilés, F. Analysis of twist stiffness of single and double-wall corrugated boards. Compos. Struct. 2014, 110, 7–15. [Google Scholar] [CrossRef]
  106. Garbowski, T.; Knitter-Piątkowska, A.; Mrówczyński, D. Numerical Homogenization of Multi-Layered Corrugated Cardboard with Creasing or Perforation. Materials 2021, 14, 3786. [Google Scholar] [CrossRef] [PubMed]
  107. Gajewski, T.; Garbowski, T.; Staszak, N.; Kuca, M. Crushing of Double-Walled Corrugated Board and Its Influence on the Load Capacity of Various Boxes. Energies 2021, 14, 4321. [Google Scholar] [CrossRef]
  108. Philadelphia Museum of Art. “Kenno” Stool. Available online: https://www.philamuseum.org/collections/permanent/325729.html (accessed on 15 September 2021).
  109. Archiproducts. Showroom Finland Presents the Novelties for 2013. Available online: https://www.archiproducts.com/es/noticias/showroom-finland-presents-the-novelties-for-2013_31596 (accessed on 15 September 2021).
  110. Archilovers. Heikki Ruoho. Designer. Available online: https://www.archilovers.com/heikki-ruoho/ (accessed on 15 September 2021).
  111. Archilovers. Järvi & Ruoho. Design Firm. Available online: https://www.archilovers.com/teams/119359/jarvi-ruoho.html#info (accessed on 15 September 2021).
  112. Harta Packaging Industries Sdn Bhd Heavy-Duty Corrugated Board. Available online: https://hartapack.com/home/wp-content/uploads/2019/03/HEAVY-DUTY-CORRUGATED-BOARD-BROCHURE.pdf (accessed on 15 September 2021).
  113. Tri-Wall Pack. Tri-Wall Pack. Available online: http://www.tri-wall.co.th/images/column_1272621251/brochuse3.pdf (accessed on 15 September 2021).
  114. Anderson, H.R. Triple-Wall Corrugated Board. U.S. Patent 2,985,553, 23 May 1961. [Google Scholar]
  115. ADF Packaging Committee. Part 15: Packaging Specifications and Classification Systems; ADF Packaging Committee: Sydney, Australia, 2000. [Google Scholar]
  116. Chairigami. Cardboard Stool. Available online: https://www.chairigami.com/product-page/cardboard-stool (accessed on 15 September 2021).
  117. Good Design Award. Cardboard (Tri-Wall) Tables. Available online: https://www.g-mark.org/award/describe/42743?locale=en (accessed on 15 September 2021).
  118. Hplusf. PunchOut Furninture. Available online: https://hplusf.com/projects/punch-out-furniture/ (accessed on 15 September 2021).
  119. MoMA. A Child’s Chair: Delight in Ownership. 2012. Available online: https://www.moma.org/explore/inside_out/2012/07/05/a-childs-chair-delight-in-ownership/ (accessed on 15 September 2021).
  120. European Committee for Standardization (CEN). EN 1728:2013 Furniture—Seating—Test Methods for the Determination of Strength and Durability; European Committee for Standardization: Brussels, Belgium, 2013. [Google Scholar]
  121. European Committee for Standardization (CEN). EN 12520:2015—Furniture—Strength, Durability and Safety—Requirements for Domestic Seating; European Committee for Standardization: Brussels, Belgium, 2015. [Google Scholar]
  122. Garbowski, T.; Marek, A. Homogenization of corrugated boards through inverse analysis. In Proceedings of the 1st International Conference on Engineering and Applied Sciences Optimization, Kos Island, Greece, 4–6 June 2014. [Google Scholar]
  123. Barbero, E.J. Introduction to Composite Materials Design; CRC Press: Boca Raton, FL, USA, 2010. [Google Scholar]
  124. Pohl, A. Strengthened Corrugated Paper Honeycomb for Application in Structural Elements. Ph.D. Thesis, ETH Zurich, Zurich, Germany, 2009. [Google Scholar]
  125. Jones, R. Mechanics of Composite Materials; Taylor & Francis: Abingdon, UK, 1998. [Google Scholar]
  126. Marek, A. Homogenization Techniques and Constitutive Modeling of Sandwich Panels. Master’s Thesis, Poznań University of Technology, Poznań, Poland, 2014. [Google Scholar]
Figure 1. Board styles: (a) single-wall; (b) double-wall; (c) triple-wall. Reprinted with permission from ref. [12]. Copyright 2021 Elsevier.
Figure 1. Board styles: (a) single-wall; (b) double-wall; (c) triple-wall. Reprinted with permission from ref. [12]. Copyright 2021 Elsevier.
Materials 14 06645 g001
Figure 2. Machine direction (MD), cross direction (CD) and through-thickness direction (ZD). Reprinted with permission from ref. [12]. Copyright 2021 Elsevier.
Figure 2. Machine direction (MD), cross direction (CD) and through-thickness direction (ZD). Reprinted with permission from ref. [12]. Copyright 2021 Elsevier.
Materials 14 06645 g002
Figure 3. 1st design with: (a) bottom discontinuity; (b) top discontinuity.
Figure 3. 1st design with: (a) bottom discontinuity; (b) top discontinuity.
Materials 14 06645 g003
Figure 4. 2nd design with: (a,b) bottom discontinuity; (c,d) top discontinuity.
Figure 4. 2nd design with: (a,b) bottom discontinuity; (c,d) top discontinuity.
Materials 14 06645 g004
Figure 5. 3rd design: (a) perspective view; (b) front view.
Figure 5. 3rd design: (a) perspective view; (b) front view.
Materials 14 06645 g005
Figure 6. 4th design: (a) preliminary open design; (b) final closed design; (c) side view.
Figure 6. 4th design: (a) preliminary open design; (b) final closed design; (c) side view.
Materials 14 06645 g006
Figure 7. Model main dimensions.
Figure 7. Model main dimensions.
Materials 14 06645 g007
Figure 8. Mapped meshes: (a,b) 1st design with bottom/top discontinuity; (c,d) 2nd design with bottom/top discontinuity; (e) 3rd design; (f) 4th design.
Figure 8. Mapped meshes: (a,b) 1st design with bottom/top discontinuity; (c,d) 2nd design with bottom/top discontinuity; (e) 3rd design; (f) 4th design.
Materials 14 06645 g008
Figure 9. Geometry variations for the 4th design.
Figure 9. Geometry variations for the 4th design.
Materials 14 06645 g009
Figure 10. Layered material.
Figure 10. Layered material.
Materials 14 06645 g010
Figure 11. 1st design with top discontinuity (red: MD; green: CD). Ply orientation: (a) I; (b) II.
Figure 11. 1st design with top discontinuity (red: MD; green: CD). Ply orientation: (a) I; (b) II.
Materials 14 06645 g011
Figure 12. Loads applied to the 1st design with bottom (a) and top (b) discontinuity.
Figure 12. Loads applied to the 1st design with bottom (a) and top (b) discontinuity.
Materials 14 06645 g012
Figure 13. Boundary conditions applied to the 1st design with bottom (a) and top (b) discontinuity.
Figure 13. Boundary conditions applied to the 1st design with bottom (a) and top (b) discontinuity.
Materials 14 06645 g013
Figure 14. Corrugated lamina showing the laminate, xyz, and lamina, 123, reference frames. Reprinted with permission from ref. [12]. Copyright 2021 Elsevier.
Figure 14. Corrugated lamina showing the laminate, xyz, and lamina, 123, reference frames. Reprinted with permission from ref. [12]. Copyright 2021 Elsevier.
Materials 14 06645 g014
Figure 15. Evolution of maximum longitudinal deflections, u, with α, for designs: (a) 2nd; (b) 3rd.
Figure 15. Evolution of maximum longitudinal deflections, u, with α, for designs: (a) 2nd; (b) 3rd.
Materials 14 06645 g015
Figure 16. Evolution of maximum deflections with α, for the 4th design: (a) w; (b) u.
Figure 16. Evolution of maximum deflections with α, for the 4th design: (a) w; (b) u.
Materials 14 06645 g016
Figure 17. Evolution of maximum stresses [MPa] with α, for the 4th design: (a) σxx; (b) σyy.
Figure 17. Evolution of maximum stresses [MPa] with α, for the 4th design: (a) σxx; (b) σyy.
Materials 14 06645 g017
Figure 18. 1st design. Vertical deflection [mm] for bottom (a,b) and top (c,d) discontinuity.
Figure 18. 1st design. Vertical deflection [mm] for bottom (a,b) and top (c,d) discontinuity.
Materials 14 06645 g018
Figure 19. 1st design. Stress σxx [MPa] for bottom (a,b) and top (c,d) discontinuity.
Figure 19. 1st design. Stress σxx [MPa] for bottom (a,b) and top (c,d) discontinuity.
Materials 14 06645 g019
Figure 20. 1st design. Stress σyy [MPa] for bottom (a,b) and top (c,d) discontinuity.
Figure 20. 1st design. Stress σyy [MPa] for bottom (a,b) and top (c,d) discontinuity.
Materials 14 06645 g020
Figure 21. 1st design. Stresses under the seating surface for top discontinuity: (a,b) σxx [MPa], (c,d) σyy [MPa].
Figure 21. 1st design. Stresses under the seating surface for top discontinuity: (a,b) σxx [MPa], (c,d) σyy [MPa].
Materials 14 06645 g021
Figure 22. 2nd design. Vertical deflection [mm] for bottom (a,b) and top (c,d) discontinuity.
Figure 22. 2nd design. Vertical deflection [mm] for bottom (a,b) and top (c,d) discontinuity.
Materials 14 06645 g022
Figure 23. 2nd design. Stress σxx [MPa] for bottom (a,b) and top (c,d) discontinuity.
Figure 23. 2nd design. Stress σxx [MPa] for bottom (a,b) and top (c,d) discontinuity.
Materials 14 06645 g023
Figure 24. 2nd design. Stress σyy [MPa] for bottom (a,b) and top (c,d) discontinuity.
Figure 24. 2nd design. Stress σyy [MPa] for bottom (a,b) and top (c,d) discontinuity.
Materials 14 06645 g024
Figure 25. 4th design. Deflections for ply orientation I [mm]: (a) vertical, w; (b) longitudinal, u.
Figure 25. 4th design. Deflections for ply orientation I [mm]: (a) vertical, w; (b) longitudinal, u.
Materials 14 06645 g025
Figure 26. 4th design. Stress for ply orientation I [MPa] (a,c) σxx; (b,d) σyy.
Figure 26. 4th design. Stress for ply orientation I [MPa] (a,c) σxx; (b,d) σyy.
Materials 14 06645 g026
Figure 27. Vertical deflections, w.
Figure 27. Vertical deflections, w.
Materials 14 06645 g027
Figure 28. Longitudinal deflections, u.
Figure 28. Longitudinal deflections, u.
Materials 14 06645 g028
Figure 29. Stresses σxx.
Figure 29. Stresses σxx.
Materials 14 06645 g029
Figure 30. Stresses σyy.
Figure 30. Stresses σyy.
Materials 14 06645 g030
Figure 31. Comparative results for the best configurations found for each design.
Figure 31. Comparative results for the best configurations found for each design.
Materials 14 06645 g031
Table 1. Common flute types [15,23]. Reprinted with permission from ref. [12]. Copyright 2021 Elsevier.
Table 1. Common flute types [15,23]. Reprinted with permission from ref. [12]. Copyright 2021 Elsevier.
DesignationPictureHeight (in)Height (mm)Flutes/mPitch (mm)Take-Up Factor
A flute Materials 14 06645 i0011/4″4.8108 ± 108.0–9.5 ≈1.50
B flute Materials 14 06645 i0021/8″3.2154 ± 105.5–6.5≈1.40
C flute Materials 14 06645 i00311/64″4.0128 ± 106.8–7.9 ≈1.45
E flute Materials 14 06645 i0041/16″1.6295 ± 133.0–3.5≈1.25
F flute Materials 14 06645 i0051/32″0.8420 ± 131.9–2.6≈1.25
Table 2. Material properties: elastic moduli, Ei, shear moduli, Gij, and Poisson ratios, νij. Reprinted with permission from ref. [12]. Copyright 2021 Elsevier.
Table 2. Material properties: elastic moduli, Ei, shear moduli, Gij, and Poisson ratios, νij. Reprinted with permission from ref. [12]. Copyright 2021 Elsevier.
ParameterUnitHeavy Duty
Outer LinerInner LinerFluting
E1MPa825081804500
E2MPa290031204500
E3MPa290031203000
G23MPa707035
G13MPa773.5
G12MPa189019501500
ν12-0.430.430.40
ν13-0.010.010.01
ν23-0.010.010.01
tmm0.750.400.25
hmm--4.8
Pmm--8.5
Table 3. Elements of the stiffness matrix for each layer, in Voigt notation.
Table 3. Elements of the stiffness matrix for each layer, in Voigt notation.
QijUnitOuter LinerInner LinerFluting
Q11[MPa]8824.28801.4146.2
Q12[MPa]1334.2144459.807
Q13[MPa]44.36148.01145.44
Q22[MPa]31023357.2361.6
Q23[MPa]35.7139.0859.755
Q33[MPa]2900.53120.6146.14
Q44[MPa]70704.5198
Q55[MPa]770.90365
Q66[MPa]189019505.9147
Table 4. Comparative results for the best configurations found for each design, for α = 80°.
Table 4. Comparative results for the best configurations found for each design, for α = 80°.
DesignArea (m2)w (mm)u (mm)σxx (MPa)σyy (MPa)
1st0.870.90.215.38.9
2nd0.831.10.31.42.6
3rd0.796.01.93.28.4
4th0.570.61.21.52.9
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Suarez, B.; Muneta, L.M.; Romero, G.; Sanz-Bobi, J.D. Efficient Design of Thin Wall Seating Made of a Single Piece of Heavy-Duty Corrugated Cardboard. Materials 2021, 14, 6645. https://doi.org/10.3390/ma14216645

AMA Style

Suarez B, Muneta LM, Romero G, Sanz-Bobi JD. Efficient Design of Thin Wall Seating Made of a Single Piece of Heavy-Duty Corrugated Cardboard. Materials. 2021; 14(21):6645. https://doi.org/10.3390/ma14216645

Chicago/Turabian Style

Suarez, Berta, Luisa M. Muneta, Gregorio Romero, and Juan D. Sanz-Bobi. 2021. "Efficient Design of Thin Wall Seating Made of a Single Piece of Heavy-Duty Corrugated Cardboard" Materials 14, no. 21: 6645. https://doi.org/10.3390/ma14216645

APA Style

Suarez, B., Muneta, L. M., Romero, G., & Sanz-Bobi, J. D. (2021). Efficient Design of Thin Wall Seating Made of a Single Piece of Heavy-Duty Corrugated Cardboard. Materials, 14(21), 6645. https://doi.org/10.3390/ma14216645

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop