energies
Article Estimation of the Compressive Strength of Corrugated Cardboard Boxes with Various Openings
Tomasz Garbowski 1 , Tomasz Gajewski 1 and Jakub Krzysztof Grabski 2,*
1 Institute of Structural Analysis, Poznan University of Technology, Piotrowo Street 5, 60-965 Pozna´n,Poland; [email protected] (T.G.); [email protected] (T.G.) 2 Institute of Applied Mechanics, Poznan University of Technology, Jana Pawła II Street 24, 60-965 Pozna´n,Poland * Correspondence: [email protected]
Abstract: This paper presents mixed analytical/numerical method for estimating the static top-to- bottom compressive strength of corrugated packaging with different ventilation openings and holes, in which the torsional and shear stiffness of corrugated cardboard as well as the panel depth-to-width ratio are included. Analytical framework bases on Heimerls assumption with a modification to a critical force, which is here computed by a numerical algorithm. The proposed method is compared herein with the successful McKee formula and is verified with the large number of experiment results of various packaging designs made of different qualities of corrugated cardboard. The results show that, for various hole dimensions or location of openings in no-flap and flap boxes, the estimation error may be reduced up to three times than in the simple analytical approach.
Keywords: corrugated board buckling; orthotropic panels; laboratory tests; box strength; openings; ventilation holes; McKee formula
1. Introduction Citation: Garbowski, T.; Gajewski, T.; The packaging industry plays an important role in the modern world. The leading Grabski, J.K. Estimation of the material used in this field is a corrugated cardboard. Corrugated boxes are used in many Compressive Strength of Corrugated branches of various industries, mainly to transport goods or products from one place Cardboard Boxes with Various to another. During the transportation and storage, the boxes are subject to considerable Openings. Energies 2021, 14, 155. loading. In order to ensure appropriate strength of the boxes, it is necessary to use ap- https://doi.org/10.3390/en14010155 propriate estimation method and/or to conduct a series of experiments to determine a
Received: 4 December 2020 behavior of the corrugated material and the whole construction of the box. The strength Accepted: 26 December 2020 estimation becomes challenging when box contains, e.g., hand holes and/or openings. Published: 30 December 2020 Correct characterization of material properties of the corrugated board and precise esti- mation technique of the box strength allows to design the carton boxes in line with the Publisher’s Note: MDPI stays neu- sustainable development strategy. In the material such as corrugated board, there are two tral with regard to jurisdictional clai- characteristic in-plane directions of orthotropy related to fiber orientation. The first one ms in published maps and institutio- is perpendicular to the main axis of the corrugation and it is called the machine direction nal affiliations. (MD). It is parallel to the paperboard fiber alignment. The second one is parallel to it and called cross direction (CD). One possibility for analysis of strength of the corrugated boxes is the execution of several iterative experiments, in the loop in which testing, validation and optimization Copyright: © 2020 by the authors. Li- of the box design take place. In the packaging industry the most commonly used are censee MDPI, Basel, Switzerland. compressive, bending or bursting strength tests of corrugated board. Another important This article is an open access article distributed under the terms and con- material property to examine is also a humidity, which has a big impact on the box strength. ditions of the Creative Commons At- However, the most important and practical are the box compression test (BCT) and the tribution (CC BY) license (https:// edge crush test (ECT). creativecommons.org/licenses/by/ Another possibility is to use analytical formulas for prediction compressive strength 4.0/). of boxes made from corrugated boards available in the literature. The formulas should
Energies 2021, 14, 155. https://doi.org/10.3390/en14010155 https://www.mdpi.com/journal/energies Energies 2021, 14, 155 2 of 20
be as simple as possible to make it useful for practical applications. Over the years, many different approaches has been proposed. The parameters used in these formulas can be categorized intro three groups: paper, board and box parameters [1]. In the first group, one can distinguish two tests—ring crush test (RCT), Concora liner test (CLT), liner type, weights of liner and fluting, corrugation ratio and a constant related to fluting. The parameters of board include thickness, flexural stiffnesses in MD and CD, ECT and moisture content. The parameters of boxes are dimensions, perimeter, applied load ratio, stacking time, buckling ratio and printed ratio. In 1952, Kellicutt and Landt used the parameters of paper (RCT, flute constant) and the box (perimeter, box constant) to propose one of the first formula for prediction of compressive strength of boxes [2]. Maltenfort related the critical force in the BCT to the parameters of the paper (CLT, liner type) and box dimensions [3]. The most popular and widely used formula has been proposed by McKee, Gander and Wachuta in 1963 [4]. They took into account the parameters of the paperboard (ECT, flexural stiffnesses) and the box perimeter. The formula is applicable for some relatively simple containers. Over the years, many researchers tried to extend the applicability of the McKee’s approach. Probably one of the first is the work of Allerby et al. [5] from 1985, in which the authors modified constant and exponents of the McKee’s formula. Two years later, Schrampfer et al. extended applicability of the McKee relationship for a wider range of cutting methods and equipment [6]. Batelka et al. included in this relationship also the dimensions of the box [7], while Urbanik et al. took into account also the Poisson’s ratio [8]. In the original paper of McKee et al. [4], an arbitrary chosen constant value was used, which makes their approach limited only for some typical boxes. This constant was later analyzed in more complicated cases by Garbowski et al. [9]. If one would like to include also the transverse shear, the additional tests are required and therefore updated formulation should be used, which was recently modified by Aviles et al. [10] and later by Garbowski et al. [11,12]. On the contrary, in the analysis of corrugated boxes numerical methods are very often used to compute the strength of the boxes. The most common method is the finite element method (FEM), which was successfully applied in many research on numerical analysis of boxes made from corrugated cardboards, including investigations of buckling phenomena [13] or transverse shear stiffness of corrugated cardboards [14,15]. In numerical simulations, the corrugated paperboards are commonly modelled as sandwich panels [16]. The main goal in the numerical analysis is to obtain results as accurate as possible with computational costs as small as possible. Thus, there is a need of using models, which are simplification of the real objects without losing adequacy to reality. In the numerical analysis of layered materials [17–19] and constructions, made from such materials, it is useful to perform a process called homogenization [20,21]. As a consequence, instead of structure of layers made from different materials one can obtain one layer characterized by effective properties of the composite. Application of the external loads to the homogenized plate should give similar effect like for the composite material with heterogeneous core. One of the homogenization methods, which was based on strain energy and can be applied for sandwich panels, was proposed by Hohe [22]. In this approach, an equivalence of a representative element of the heterogeneous and homogenized elements was assumed. A period homogenization method, in which an equivalent membrane and pure bending characteristic of period plates is obtained, has been also proposed by Buannic et al. [23]. They applied a modified version of this approach to take into account the transfer shear effect in the analysis of sandwich panels. Biancolini used the FEM for analysis of a micromechanical part of the considered plate [24]. Appling the energy equivalency between the model and the equivalent plate, he obtained the stiffness properties of the considered plate. Abbès and Guo decomposed the plates into two beams in directions of the plate to obtain the torsion rigidity of the orthotropic sandwich plates [25]. Two approaches to homogenization have been compared by Marek and Garbowski. These considerations were based on the classical laminated plate theory and the deformation energy equivalence method [26] or based on inverse analysis [27]. Energies 2021, 14, x FOR PEER REVIEW 3 of 20
These considerations were based on the classical laminated plate theory and the defor- Energies 2021, 14, 155 mation energy equivalence method [26] or based on inverse analysis [27]. 3 of 20 The strength of the boxes made from corrugated cardboards can be reduced by many factors, including moisture, prints, perforations, etc. However, one of the most common factor for loss of the strength of boxes is the presence of openings on the walls of the pack- The strength of the boxes made from corrugated cardboards can be reduced by many age, as most of the boxes contain ventilation or at least hand holes [28]. This issue was a factors, including moisture, prints, perforations, etc. However, one of the most common subject of some recently published papers. Fadji et al. [29] investigated an effect of the factor for loss of the strength of boxes is the presence of openings on the walls of the geometrical features of the holes, including height, area, orientation, shape and number package, as most of the boxes contain ventilation or at least hand holes [28]. This issue was of holes on the strength of boxes. The influence of ventilation holes on the stability and a subject of some recently published papers. Fadji et al. [29] investigated an effect of the strength of the boxes using FEM was examined by Fadiji et al. [30]. They considered cer- geometrical features of the holes, including height, area, orientation, shape and number tain type of boxes, which are commonly used in the fresh fruit industry, in which the of holes on the strength of boxes. The influence of ventilation holes on the stability and ventilation holes play an important role. strength of the boxes using FEM was examined by Fadiji et al. [30]. They considered certain typeTo of boxes,the best which knowledge are commonly of the authors, used in there the fresh is no fruit research industry, in which in which the effect the ventilation of losing strengthholes play caused an important by holes role. is taken into account in the analytical prediction of the box strength.To the In this best paper, knowledge we extend of the the authors, applicability there isof nothe research approach in presented which the in effect [9] for of corrugatedlosing strength boxes caused with holes. by holes is taken into account in the analytical prediction of the box strength. In this paper, we extend the applicability of the approach presented in [9] for 2.corrugated Materials boxesand Methods with holes. 2.1. Ultimate Compressive Strength of a Single Panel 2. Materials and Methods The most popular and the simplest approach used these days for computation of ul- timate2.1. Ultimate compressive Compressive load of Strength corrugated of a Single cardboar Paneld boxes was proposed by McKee et al. [4]. The startingThe most point popular in their and study the was simplest the determination approach used of thesethe ultimate days for plate computation load of ofa rectangularultimate compressive panel with load dimensions of corrugated × cardboard, loaded vertically boxes was on proposed the upper by McKeeedge, in et which al. [4]. theThe plate starting under point consideration in their study is a wasseparated the determination panel of a box of (Figure the ultimate 1a), despite plate the load factPf thatof a therectangular box is loaded panel as with a three-dimensional dimensions a × b structure, loaded vertically during the on compression the upper edge, test in(Figure which 1b). the Allplate four under edges consideration of the each separated is a separated plate panel are pinned of a box in (Figure the out-of-plane1a), despite direction, the fact thatand thethe bottom box is loaded edge is as additionally a three-dimensional constrained structure along during in the the vertical compression direction. test The (Figure basis1 forb). theAll derivation four edges given of the by each McKee separated et al. is plate the empirical are pinned formula in the out-of-planepresented in direction,Heimerl [31], and the bottom edge is additionally constrained along b in the vertical direction. The basis for where the load at failure [N/mm] is the result of a combination of buckling force and compressivethe derivation strength, given by namely: McKee et al. is the empirical formula presented in Heimerl [31], where the load P at failure [N/mm] is the result of a combination of buckling force and f compressive strength, namely: = r , (1) Pf ECT = k , (1) Pcr Pcr where is an exponent, ∈(0,1), is a certain constant, is a critical load of a corru- gatedwhere panelr is an [N/mm] exponent, andr ∈ (0, 1is), ktheis aedge-loaded certain constant, compressivePcr is a critical strength load of of a a corrugated corrugated boardpanel [N/mm]in cross direction and ECT [N/mm].is the edge-loaded compressive strength of a corrugated board in cross direction [N/mm].
(a) (b) (c)
Figure 1. (a) A single panelFigure of width 1. (a)b Aand single height panela separatedof width from and cardboardheight separated packaging from as acardboard supported packaging plate under as a compression; (b) cardboardsupported packaging inplate press under during compression; box compression (b) cardboard testing; (packagingc) box compression in press pressduring used box in compression laboratory tests [32]. testing; (c) box compression press used in laboratory tests [32].
Energies 2021, 14, 155 4 of 20
To include the nonlinear material effect when no elastic buckling occurs the above formula was later modified by Urbanik et al. [8,13]:
un Pf P = k cr , (2) ECT ECT
where n = 1 − r. Integer u plays an activating role and if it is equal 1, it corresponds to elastic buckling (ECT > Pcr), while if it is equal 0, it corresponds to inelastic buckling, so the ultimate load Pf is proportional to ECT, namely:
Pf = k ECT. (3)
In both cases, if the selected panel of a corrugated box does not have any openings, see Figure1a, the analytical formula describing the critical load for a rectangular orthotropic plate [9,33–35] can be utilized: √ π2 D D Pb = k 11 22 , (4) cr cr b2
b where Pcr is the critical load of the panel with dimensions a × b (a—panel height, b—panel width [mm]), kcr is a dimensionless buckling coefficient that depends, among others, on the ratio a/b, boundary conditions applied to the plate edges, the buckling mode and a material characteristics, D11 is the effective bending stiffness in the MD [N/mm] and D22 is the effective bending stiffness in the CD [N/mm]. In this study, the first buckling mode was used for computations in all cases. In this buckling mode, the displacements of panel with and without hole usually represent a shape of the half-wave in both directions. All fibrous materials, including corrugated cardboard, are characterized by orthotropy, which means that the mechanical properties in such materials change depending on the direction, here CD/MD. The buckling coefficient for such panels takes the following form:
s 2 s D11 mb (D12 + 2D33) D22 a 2 kcr = + 2 √ + , (5) D22 a D11D22 D11 mb
where m defines a buckling mode and counts the half-waves along the load direction, D33 is an effective torsion stiffness [N/mm] and D12 = D11ν21 = D22ν12, ν12 and ν21 are in-plane Poisson’s ratios.
2.2. Buckling—McKee’s Formula If the analyzed panel does not contain any holes, a simplified formula for critical load computation can be utilized. Such a formula was suggested by McKee et al. [4], in order to maximally simplify the complicated equations describing the critical load of compressed orthotropic plates so that the use of these equations could be common and practical. Therefore, the authors replaced Equation (5) with a simpler formula:
mb 2 θa 2 k = + + 2K, (6) cr θa mb q where θ = D−1D and K is a plate parameter assumed to be equal to 0.5, which leads 11 22√ to: D12 + 2D33 = 0.5 D11D22. They also assumed that the plate has equal sides, i.e., b = c and therefore b = Z/4 (where Z is the box perimeter), which led to further simplification of Equation (5) to the form: √ (4π)2 D D Pb = 11 22 k . (7) cr Z2 cr Energies 2021, 14, 155 5 of 20
Energies 2021, 14, x FOR PEER REVIEW 5 of 20
2.3. Buckling—Finite Element Method 2.3. Buckling—FiniteHowever, if the Element analyzed Method panel contains any openings, the analytical approach be- comesHowever, impractical if the due analyzed to the lack panel of an contains universal any solution openings, for athe panel analytical with different approach holes. be- Incomes such impractical a case, the numerical due to the method lack of worksan universal best, although solution it for is morea panel difficult with different to implement, holes. itIn is such very a general case, the and numerical allows to method determine works the best, critical although load of it panels is more with difficult any shape to imple- and boundaryment, it is conditions, very general including and allows holes to with determine complex the shapes. critical In load this of study, panels the with buckling any shape was computedand boundary for each conditions, separated including panel in theholes 3D with space complex having allshapes. external In this edges study, fixed the in out buck- of planeling was direction, computed additionally for each the separated bottom edge panel was in fixedthe 3D in thespace vertical having direction, all external other edges BC’s werefixed inactive. in out of plane direction, additionally the bottom edge was fixed in the vertical di- The critical load can be calculated from the formula: rection, other BC’s were inactive. The critical load can be calculated fromb the formula: Pcr = λq, (8) ̅ = , (8) where λ is a critical loading multiplier, i.e., the smallest eigenvalue and q is a distributed loadwhere on a̅ panel is a critical upper loading edge [N/mm]. multiplier, i.e., the smallest eigenvalue and is a distributed loadTo on calculate a panel upper the critical edge force[N/mm]. multiplier, the following two-step algorithm should be used.To In thecalculate first step, the critical we need force to: multiplier, the following two-step algorithm should be •used.discretize In the first the step, panel we with need finite to: elements (FE), see Figure2a, •• computediscretize the the stiffness panel with matrix finiteKe elementsof each element, (FE), see Figure 2a, •• assemblecompute thethe globalstiffness stiffness matrix matrix of Keach0 of element, a whole panel, •• ∗ computeassemble nodal the global forces stiffness representing matrix initial of loading a whole configuration panel, P , i.e., for loading ∗ • multipliercompute nodalλ = 1 forces (one-parameter representing loading initial assumed loadingP configuration= λP ), ∗, i.e., for loading • takemultiplier boundary =1 conditions (one-parameter into account, loading assumed = ∗), ∗ ∗ •• solvetake boundary equation 0conditionsd = P to into obtain account, nodal displacements in pre-buckling state: • ∗ ∗ solve equation = to obtain nodal displacements in pre-buckling state: d∗ = K−1·P∗, (9) ∗ 0 ∗ = · , (9) • extract element displacements d(e)∗ (from displacements of the system d∗) and com- • extract element displacements ( )∗ (from displacements of the system ∗) and com- pute in each element the generalized stresses s(e)∗. pute in each element the generalized stresses ( )∗. In the second step, we need to: In the second step, we need to: • generate initial geometrical stress matrices for each element Ke (s∗e) and the whole • generate initial geometrical stress matrices for each element σ ( ∗ ) and the whole panel K (s∗), panel σ ( ∗), • formulate initial buckling problem: • formulate initial buckling problem: [K + λK ]v = 0, (10) + 0 = ,σ (10) • • solvesolve thethe eigenproblemeigenproblem to to determine determine the the pairs pairs(λ (1 , v , 1) ,)...,…,, ((λ N ,, v N)),, where N isis aa = numbernumber ofof degrees degrees of of freedom, freedom,λi is i this eigenvalue, th eigenvalue,vi ∆ d=∆ i is ith eigenvectoris th eigenvector (post- buckling(post-buckling deformation deformation mode). mode).
(a) (b) (c)
FigureFigure 2. 2.( a(a)) triangular triangular mesh mesh on on a a singlesingle panelpanel withwith circularcircular opening;opening; ((bb)) eighteeneighteen degreesdegrees ofof freedomfreedom triangulartriangular Reissner-Reissner- MindlinMindlin plate plate FE; FE; ( c()c) isoperimetric isoperimetric formulation formulation of of 6-node 6-node triangle triangle FE. FE.
Energies 2021, 14, 155 6 of 20
In order to compute stiffness matrix of the plate, the proper formulation of the FE must be used. In the classical plate FE, the field variables are approximated through shape functions according to the associated node values as follows:
Nn d = ∑ Nj(x, y)dj, (11) j=1
where Nj(x, y) is a shape function associated with node j and Nn is a number of nodes in the element: T dj = uj vj wj θxj θyj , (12) is the displacement vector in the j-th node. The membrane, bending and shear strains associated to the displacement in Equation (11) can be therefore obtained as follows:
m ε = ∑ Bj dj, (13) j
b κ = ∑ Bj dj, (14) j s γ = ∑ Bj dj, (15) j where: Nj,x 0 0 0 0 m Bj = 0 Nj,y 0 0 0 , (16) Nj,y Nj,x 0 0 0 0 0 0 Nj,x 0 b Bj = 0 0 0 0 Nj,y , (17) 0 0 0 Nj,y Nj,x s 0 0 Nj,x Nj 0 Bj = . (18) 0 0 Nj,y 0 Nj The energy functional of the plate can be finally obtained as follows:
Z Z Z 1 T T Π = d B DBdA + S Ds SdA d − q w dA d, (19) 2 A A A
T T m T b T s T where B = (B ) , B , S = (B ) ; D and Ds are stiffnesses of the plates [9]. By using Lagrange’s equations for the Equation (19), the characteristic equation of the system is obtained as follows: K d = P, (20)
where: Z T T K = B DB + S Ds S dA, (21) A
and: Z P = q N dA. (22) A As the geometry of each panel does not have to be regular, especially in the case of panels with openings, it was assumed that the FE mesh will be made of shell triangular elements. Most Reissner-Mindlin triangular shell FE described in the literature were developed using the assumed transverse shear strain techniques. Here, the 6-noded quadratic triangle FE with assumed linear shear strains was utilized. Zienkiewicz et al. [36] developed such a 6-noded plate triangle with a standard quadratic interpolation for the Energies 2021, 14, 155 7 of 20
rotations and the deflection, and the assumed linear transverse shear strain field in the following form: γ = α + α ξ + α η, ξ 1 2 3 (23) γη = α4 + α5ξ + α6η, i.e.,: 1 ξ η 0 0 0 A = . (24) 0 0 0 1 ξ η
Figure2b,c show the position of the six shear sampling points and the ξi directions. Following the procedure presented in AppendixA, we find: 1 ξ1 η1 0 0 0 1 ξ_2 η2 0 0 0 −a −aξ −aη a aξ aη P = 3 3 3 3 , (25) −a −aξ −aη a aξ aη 4 4 4 4 0 0 0 1 ξ5 η5 0 0 0 1 ξ6 η6
100000000000 001000000000 0 0 0 0 −a a 0 0 0 0 0 0 T = , (26) 0 0 0 0 0 0 −a a 0 0 0 0 000000000100 000000000001
J1 0 J2 C = . , (27) .. 0 J6 B1 2 ^ B Bs = . , (28) . B6 where: √ 2 a = . (29) 2
Matrix Bs is computed by Equation (A24) in AppendixA. A full 3-point Gauss quadra- ture is used here for the numerical integration of element transverse shear stiffness matrix: Z (e) T (e) Ks = S Ds S dA . (30) A(e) This element performs well for most of applications, although sometimes is too flexible [36,37]. Improvements can be introduced by imposing a linear variation of the normal rotation θn along all three triangle sides as:
1 θ − θ + θ , (31) ni 2 ni−2 ni+2 where i = 4, 5, 6 are the mid-side nodes (Figure2a). Equation (31) can be later used for eliminating the normal rotation at the midsize nodes [36,38]. Energies 2021, 14, 155 8 of 20
2.4. Box Compression Strength—McKee’s Formula To compute failure load of a single panel we substitute the buckling force derived from Equation (7) into Equation (1), after [4] we have:
1−r rp −2(1−r) 2(1−r) 1−r Pf = k ECT D11D22 Z (4π) kcr , (32)
which can be simplified to:
1−r rp −2(1−r) Pf = k1ECT D11D22 Z , (33)
1−r where k1, under the assumption [4] that kcr is constant if a/Z > 1/7 and equals approxi- mately 1.33, while r = 0.746 and k = 0.4215, is equal to:
2(1−r) ∼ k1 = 1.33 k (4π) = 2.028. (34)
To calculate the failure load of the entire box, it is enough to multiply Equation (33) by the Z, which is the perimeter of the box. The original assumption [4] was that the width and length of a box are equal (i.e., b = c) thus the compressive strength of a box, known as the long McKee formula, reads:
1−r rp 2r−1 BCTMK1 = k1ECT D11D22 Z . (35) √ ∼ 2 The empirical observations [4] that D11D22 = 66.1 ECT h , led to a further simplifi- cation and finally to the most known short form of McKee’s equation:
2(1−r) 2r−1 BCTMK2 = k2ECT h Z , (36)
where k2, with previously assumed value of r = 0.746, equals:
1−r ∼ k2 = k1(66.1) = 5.874, (37)
So, Equation (36) takes the final form:
0.508 0.492 BCTMK2 = 5.874 ECT h Z . (38)
This is the most commonly used formula. However, its accuracy is acceptable only for the simplest boxes. The McKee formula should not be used if the length to width ratio or the height to length ratio of the box is too large [9]. The same problem appear when box contains any openings or ventilation holes. To overcome this limitations we have derived an analytical formulation for box strength from Equation (1), in which the critical force is computed using numerical methods instead of analytical formulas. Both will be shown in the following sections.
2.5. Box Compression Strength—General Case In order to calculate the compressive strength of a single panel of a box, see Figure1a, the force Pf described by the formula (1) should be multiplied by the i-th panel width b or c. To get the failure load of an entire box (BCT), it is required to sum up the compressive strength of all panels. Thus, in the general case, we get the following: ! 2 1−r 2 r b b h c 1−r b i BCT = k ECT Pcr γp γbb + (Pcr) γp γcc , (39) ∑ i i ∑ i i i=1 i=1
b c where (·)i indicates the i-th panel of width b or c (i = 1, ... , 2). Pcr and Pcr are the critical b c forces of panels of width b and c, respectively (see Figure1a). γp and γp are the reduction Energies 2021, 14, 155 9 of 20
factors taking into account the ratio of the compressive strength of the plate with and without a hole for plate width b and c, respectively. γb and γc are the reduction coefficients due to in-plate aspect ratio of the box dimensions, defined as: √ γb = b/c, γc = 1 i f b ≤ c, (40) √ γc = c/b, γb = 1 i f b > c. (41)
2.6. Practical Considerations Equation (39) uses different definitions of the critical force, i.e., Equations (4) and (8), depending on whether the i-th panel contains any hole or not. In any case, the values of constants k and r in Equation (39) are required. As suggested by other authors [4,39–41], k = 0.5 and r = 0.75 were assumed. The compressive strength of rectangular boxes with various openings (with base dimensions b × c, as shown in Figure1a) after taking the typical values of k and r in Equation (39) takes the form: ! 2 0.25 2 0.75 b b h c 0.25 b i BCT = 0.5 ECT Pcr γp γbb + (Pcr) γp γcc , (42) ∑ i i ∑ i i i=1 i=1
b c where Pcr and Pcr are defined in: (a) Equation (4)—for orthotropic rectangular plates without holes (analytical), (b) Equation (8)—for orthotropic rectangular plates with various openings (numerical – FEM). c In all of the above equations describing the critical force Pcr, all definitions of width b should be replaced by c, e.g., Equation (4) should read: √ π2 D D Pc = k 11 22 . (43) cr cr c2
3. Results 3.1. Box Compression Strength—Experiment vs. Estimation In the following section, experimental and computational examples regarding boxes are presented and the applicability of the analytical-numerical approach is validated. In the paper, the analyzed cases represent the simplest slotted-type boxes with holes varied in size, position at the panel side and location on the longer or shorter wall. Two typical designs of FEFCO boxes with an addition of different holes were considered, i.e., FEFCO F200—packaging without upper flaps, see Figure3, and FEFCO F201—packaging with flaps, see Figure4. Two box dimensions were considered, namely, 300 × 200 × 200 mm and 300 × 200 × 300 mm (length, width and height, respectively). Moreover, three positions of the opening at the panel side were assumed, i.e., in the upper corner (O1), in the middle of the package height at the edge of the panel (O2) and in the center of the panel (O3). Also, two cases were considered, in which the opening was at the shorter wall (B), and at the longer wall (L) of the box. In most cases, the openings were square, but rectangular openings were also considered. The hole dimensions for cases O1B1S, O2B1S, O3B1S, O1L1S, O2L1S and O3L1S were the same, namely, 100 × 100 mm. For O3B1S1a, O3B1S1b, O3B1S1c and O3B1S1d they varied with dimensions of 85 × 85, 120 × 120, 150 × 150 and 200 × 200 mm, respectively. For O3L1S1a, O3L1S1b, O3L1S1c, O3L1S1d and O3L1S1e they varied with dimensions of 85 × 85, 120 × 120, 150 × 150, 200 × 200 and 250 × 250 mm, respectively. For O3L1Ra and O3L1Rb the dimensions were 300 × 200 and 300 × 250 mm. EnergiesEnergies2021 ,202114, 155, 14, x FOR PEER REVIEW 10 of10 20 of 20 Energies 2021, 14, x FOR PEER REVIEW 10 of 20
(a) (b) (c) (a) (b) (c)
(d) (e) (f) (d) (e) (f) FigureFigure 3. Designs 3. Designs of boxes of boxes without without upper upper flaps flaps with openingswith openings considered considered in the study—modificationsin the study – modifica- Figure 3. Designs of boxes without upper flaps with openings considered in the study – modifica- tions of FEFCO F200 with box height of 300 mm: (a) O1L1S, (b) O2L1S, (c) O3L1S and box height oftions FEFCO of FEFCO F200 with F200 box with height box height of 300 mm: of 300 (a) O1L1S,mm: (a) ( bO1L1S,) O2L1S, (b ()c O2L1S,) O3L1S ( andc) O3L1S box height and box of 200 heightmm of 200 mm (d) O1B1S, (e) O2B1S, and (f) O3B1S. (ofd) 200 O1B1S, mm (e ()d O2B1S,) O1B1S, and (e) ( fO2B1S,) O3B1S. and (f) O3B1S.
(a) (b) (c) (a) (b) (c)
(d) (e) (f) (g) (d) (e) (f) (g)
(h) (i) (j) (k) (l) (h) (i) (j) (k) (l)
(m) (n) (o) (m) (n) (o) Figure 4. Designs of flap boxes with openings considered in the study—modifications of FEFCO F201 with box height of FigureFigure 4.4. DesignsDesigns ofof flapflap boxesboxes withwith openingsopenings consideredconsidered inin thethe study—modificationsstudy—modifications ofof FEFCOFEFCO F201F201 withwith boxbox heightheight ofof 300 mm: (a) O1L1S, (b) O2L1S, (c) O3L1S, (d), O3B1S1a (e) O3B1S1b, (f) O3B1S1c, (g) O3B1S1d, (h) O3L1S1a, (i) O3L1S1b, 300 mm: (a) O1L1S, (b) O2L1S, (c) O3L1S, (d), O3B1S1a (e) O3B1S1b, (f) O3B1S1c, (g) O3B1S1d, (h) O3L1S1a, (i) O3L1S1b, 300 (mm:j) O3L1S1c, (a) O1L1S, (k) O3L1S1d, (b) O2L1S, (l ()c )O3L1S1e O3L1S, (andd), O3B1S1a with box ( heighte) O3B1S1b, of 200 (f )mm: O3B1S1c, (m) O1B1S, (g) O3B1S1d, (n) O2B1S, (h) O3L1S1a,and (o) O3B1S. (i) O3L1S1b, ((jj)) O3L1S1c,O3L1S1c, ((kk)) O3L1S1d,O3L1S1d, ((ll)) O3L1S1eO3L1S1e andand withwith boxbox heightheight ofof 200200 mm: mm: (m ()m O1B1S,) O1B1S, (n )(n O2B1S,) O2B1S, and and (o )(o O3B1S.) O3B1S. For the purposes of this study, the selected boxes were tested in a box compression For the purposes of this study, thethe selectedselected boxesboxes werewere testedtested inin aa boxbox compressioncompression testing machine, namely BCT-19T10 from FEMat [32]; see Figure 1c. According to the spec- testingtesting machine, machine, namely namely BCT-19T10 BCT-19T10 from from FEMat FEMat [32]; [32 see]; see Figure Figure 1c. 1Accordingc. According to the to spec- the ifications, the machine measures the force up to 10 kN with a resolution of 0.1 N. The specifications,ifications, the machine the machine measures measures the force the force up to up 10 to kN 10 with kN witha resolution a resolution of 0.1 of N. 0.1 The N. displacement accuracy is 0.001 mm. In order to get reliable results, 4–6 samples of each Thedisplacement displacement accuracy accuracy is 0.001 is 0.001mm. In mm. order In to order get reliable to get reliable results,results, 4–6 samples 4–6 samples of each packaging design were tested. In the study, for all cases the number of tested samples packaging design were tested. In the study, for all cases the number of tested samples
Energies 2021, 14, 155 11 of 20
of each packaging design were tested. In the study, for all cases the number of tested samples exceeded 180. The boxes were manufactured from different three-layer corrugated cardboard, namely, 350, 380, 400, and 480 g/m2. All those cardboards had B or E flute, their thickness was 1.49, 1.52, 1.59, 2.80 and 2.82 mm. All cases tested are presented in Tables1 and2 for FEFCO type boxes of F200 and F201, respectively. In Tables1 and2 “opening type” columns contain the information about the opening, for instance O1B1S, means that the opening was located in the corner (O1) of the shorter wall (B). “Cardboard quality” columns contain the code of cardboard, for instance 3B400-1, means that three layer cardboard and B flute was used with 400 g/m2, 1st material production.
Table 1. Experimental data for FEFCO F200 boxes for different: opening type, box height, cardboard quality and box compression strength.
No. Opening Type Box Height [mm] Cardboard Quality Box Strength [N] 1 O3B1S 200 3E350-3 901.4 2 O2B1S 200 3E350-3 742.3 3 O2B1S 200 3B400A1-1 1777.9 4 O1B1S 200 3B400A1-1 1842.7 5 O3B1S 200 3B400A1-1 2101.6 6 O1L1S 300 3B400A1-3 1819.4 7 O2L1S 300 3B400A1-3 1903.0 8 O3L1S 300 3B400A1-3 2184.3
Table 2. Experimental data for FEFCO F201 boxes for different: opening type, box height, cardboard quality and box compression strength.
No. Opening Type Box Height [mm] Cardboard Quality Box Strength [N] 9 O1B1S 200 3E350-3 694.3 10 O3B1S 200 3E350-5 783.8 11 O1L1S 300 3E350-5 699.3 12 O3L1S 300 3E350-5 842.9 13 O2B1S 200 3E350-5 697.5 14 O3B1S 200 3E380A2-1 944.6 15 O2L1S 300 3E380A2-1 935.7 16 O3L1S 300 3E380A2-1 1085.3 17 O1L1S 300 3E380A2-2 854.8 18 O1B1S 200 3B400A1-1 1630.0 19 O3B1S 200 3B400A1-1 1902.9 20 O2B1S 200 3B400A1-1 1629.0 21 O1L1S 300 3B400A1-3 1647.5 22 O2L1S 300 3B400A1-3 1701.9 23 O3L1S 300 3B400A1-3 2133.6 24 O3B1S1a 300 3B480-1 1841.1 25 O3B1S1b 300 3B480-1 1717.9 26 O3B1S1c 300 3B480-1 1606.6 27 O3L1S1e 300 3B480-1 1362.9 28 O3L1S1b 300 3B480-1 1743.5 29 O3L1S1c 300 3B480-1 1772.4 30 O3L1S1d 300 3B480-1 1591.3 31 O3L1S1a 300 3B480-1 1782.3
During the tests, the boxes were folded by one operator manually, and the flaps of the packaging were taped at the bottom and top, if applicable. All packaging structures were cut on a plotter, so the material was not converted. Also, each box was inspected for visual damage of the cardboard before it was submitted to displacement control pro- tocol tests. TAPPI standard T402 was applied, the samples were laboratory conditioned, i.e., temperature 23 ◦C ± 1 ◦C and relative humidity of 50% ± 2% was hold. Figure5 demonstrates examples of laboratory measurements for selected designs of packaging with Energies 2021, 14, x FOR PEER REVIEW 12 of 20
visual damage of the cardboard before it was submitted to displacement control protocol tests. TAPPI standard T402 was applied, the samples were laboratory conditioned, i.e., Energies 2021, 14, 155 12 of 20 temperature 23°C ± 1°C and relative humidity of 50% ± 2% was hold. Figure 5 demon- strates examples of laboratory measurements for selected designs of packaging with holes, namely F200–O1B1S, F201–O3B1S1c and F201–O3L1S1d. The curves of force versus displacementholes, namely are F200–O1B1S, shown; a good F201–O3B1S1c repeatability and of F201–O3L1S1d. the measurement The results curves is ofobserved. force versus The samplesdisplacement were arenot shown;inspected a good after repeatabilitythe measurem ofents the in measurement order to determine results isthe observed. deformation The orsamples damage were mechanism. not inspected after the measurements in order to determine the deformation or damage mechanism.
(a) (b) (c)
Figure 5.5. ExampleExample measurementsmeasurements from from a boxa box compression compression test test machine machine [32] [32] for different for different boxes boxes and holes: and holes: (a) F0200—O1B1S, (a) F0200— (O1B1S,b) F201—O3B1S1c (b) F201—O3B1S1c and (c) F201O3L1S1d.and (c) F201O3L1S1d.
Moreover,Moreover, the measured strength of the boxes with openings selected in this study serveserve as as a a validation validation data data for for those those obtained obtained with with analytical-numerical analytical-numerical approach approach proposed proposed here,here, see Equation (39). In In order order to to use use Equation Equation (39) (39) the critical buckling loads of box panels (with and without openin openings)gs) were required. The material data used for computing criticalcritical loads byby FEMFEM are are presented presented in in Table Table3. The 3. The cardboard cardboard qualities qualities were were also presented also pre- sentedin Tables in1 Table and2 .1 All and input Table data 2. All for input Equation data (39) for usedEquation for the (39) computations used for the were computations attached werein Supplementary attached in Supplementary Materials, separately Materials, for F200separately and F201. for F200 and F201.
TableTable 3. Material data for the cardboard qualities used in FEM computations of critical loadsloads [[42,43].42,43].
No. Cardboard Quality Thickness [mm] ECT D [N/mm] D [N/mm] D [N/mm] R [N/mm] R [N/mm] No. Cardboard Quality Thickness [mm] ECT 11 [N/mm] 22 33 44 [N/mm] 55 [N/mm] [N/mm] [N/mm] 1 3B400A1-1 2.82 5.79 3484.8 1789.8 2262.1 7.73 12.78 12 3B400A1-1 3B400A1-3 2.80 2.82 5.50 5.79 3443.53484.8 1789.8 1565.5 2262.1 2115.5 7.73 6.09 12.78 11.30 23 3B400A1-3 3B480-1 2.82 2.80 5.29 5.50 3491.23443.5 1565.5 1820.1 2115.5 2359.9 6.09 5.66 11.30 12.40 34 3B480-1 3E350-3 1.49 2.82 3.96 5.29 3491.2 958.9 1820.1 431.5 2359.9 870.2 5.66 2.49 12.40 2.79 45 3E350-3 3E350-5 1.49 1.49 4.68 3.96 878.6958.9 431.5 376.9 870.2 904.0 2.49 2.86 2.79 2.97 56 3E380A2-13E350-5 1.59 1.49 5.41 4.68 1272.1878.6 376.9 505.4 1084.5904.0 2.86 3.83 2.97 3.86 68 3E380A2-1 3E380A2-2 1.52 1.59 5.31 5.41 1042.91272.1 505.4 445.1 1084.5 963.2 3.83 3.39 3.86 3.70 8 3E380A2-2 1.52 5.31 1042.9 445.1 963.2 3.39 3.70 In Figure6, the data from the tests were presented graphically (black plot), also separatelyIn Figure for 6, F200 the (Figuredata from6a) the and tests F201 were (Figure presented6b) boxes. graphically Table1, Table (black2 andplot), Figure also sep-6 aratelyshow the for mean F200 (Figure values of6a) the and ultimate F201 (Figure load measured 6b) boxes. at Table the box 1, Table failure 2 and (BCT). Figure In Figure 6 show6, thealso mean the values values obtained of the ultimate from the load analytical-numerical measured at the approachbox failure proposed (BCT). In were Figure presented 6, also the(blue values plot), obtained see Equation from (39).the analytical-num In the analytical-numericalerical approach approach proposed presented, were presentedk and r (blueused areplot), the see typical Equation values (39). taken In the from analytical-numerical the literature, namely approachk = 0.5 presented,and r = 0.75 and, see usedEquation are the (42). typical values taken from the literature, namely =0.5 and =0.75, see Equation (42).
Energies 2021, 14, 155 13 of 20 Energies 2021, 14, x FOR PEER REVIEW 13 of 20
(a)
(b)
Figure 6. MeasuredFigure box compression 6. Measured boxstrengths compression (black line), strengths computed (black by line), a method computed proposed by a method with typical proposed (blue with line) and optimal (red line)typical parameters (blue of line) and and optimal, and McKee (red line) formula parameters (dashed of magentak and r, andline) McKee for FEFCO formula type (dashed boxes of magenta (a) F200 and (b) F201. line) for FEFCO type boxes of (a) F200 and (b) F201.
The errorThe of the error method of the proposed method proposed may be calculated may be calculated by the comparing by the comparing the computed the computed values withvalues the ultimatewith the loadsultimate measured, loads measured, namely: namely: