Next Article in Journal
Effect of Blend Composition on Barrier Properties of Insulating Mats Produced from Local Wool and Waste Bast Fibres
Previous Article in Journal
Influence of Thickness of Opaque Porcelain and Alloy Color on the Final Color of Porcelain-Fused-to-Metal Restorations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simplified Modelling of the Edge Crush Resistance of Multi-Layered Corrugated Board: Experimental and Computational Study

by
Tomasz Garbowski
1,*,
Anna Knitter-Piątkowska
2 and
Piotr Winiarski
3
1
Department of Biosystems Engineering, Poznan University of Life Sciences, Wojska Polskiego 50, 60-627 Poznan, Poland
2
Institute of Structural Analysis, Poznan University of Technology, Piotrowo 5, 60-965 Poznan, Poland
3
Schumacher Packaging Sp. z o.o., Wrocławska 66, 55-330 Krępice, Poland
*
Author to whom correspondence should be addressed.
Materials 2023, 16(1), 458; https://doi.org/10.3390/ma16010458
Submission received: 3 December 2022 / Revised: 28 December 2022 / Accepted: 30 December 2022 / Published: 3 January 2023

Abstract

:
The edge crush test is the most popular laboratory test in the corrugated packaging industry. It measures the edge crush resistance of a sample in the cross-fiber direction (CD), also known as the ECT index. This parameter is widely used for the specification of the board by its producers. It is also utilized in most analytical formulas describing the load capacity of the packaging. On the other hand, the ECT value can be estimated from both analytical and numerical models based on the basic parameters of each constituent paper. Knowing the compressive strength in CD (commonly known as SCT) and the elastic properties of the individual layers, the sample geometry (i.e., the period and height of the corrugated layer), as well as the boundary conditions, the ECT value can be calculated. This is very useful as new boards can be virtually analyzed before being manufactured. In this work, both detailed numerical models based on finite elements (FE) methods and very simple analytical (engineering) models were used for the ECT calculations. All presented models were validated with experimental data. The surprising consistency and high precision of the results obtained with the simplest approach was additionally analyzed in the study.

1. Introduction

Corrugated cardboard plays a very important role in the economies of fast-developing countries. In a world where the social awareness of consumers related to, for example, environmental protection, a sustainable economy, and the closed product cycle is gradually growing, there is no longer any place for products that threaten or may endanger the goals of sustainable development. Earth’s resources are not inexhaustible and one should be aware of the dangers and threats associated with their careless exploitation. Producers of corrugated board and paper—although their business is based on the exploitation of one of the most important raw materials—wood obtained from logging does not affect the environment to the same extent as producers of packaging made of other, less eco-friendly materials. Since corrugated cardboard is a biodegradable material, it does not leave an imprint on the natural environment in the form of piles of hard to disposed-of material or plastic islands floating on seas and oceans. Moreover, through the use of recycled fibers in the production of paper, the factor related to deforestation is also gradually minimized. By extending the life cycle of recycled fibers in the paper production process, the amount of cut trees is reduced and, concurrently, more time is gained to rebuild forested areas.
Corrugated cardboard boxes protect their contents throughout the distribution chain and are currently used for the packaging of the vast majority of goods sold and conveyed commercially. The strength of the transport box is influenced by the weather conditions in which it is transported, resulting in the box moisture content [1,2], the method of transport itself, as well as the conditions and manner of storage [3,4,5]. The compressive strength of the box, usually measured while performing the box compression test (BCT) [6,7,8,9,10], is a key measure of the packaging resistance to mechanical stress during transport and storage.
Since the invention of corrugated cardboard and its first use in the production of packaging, engineers and scientists around the world have been looking for the appropriate methods to reliably and easily estimate the load-bearing capacity of the boxes in order to optimize the consumption of the raw material and, at the same time, adequately protect the products transported inside them. The simplest formulas, however, do not allow for a credible estimation of each type of packaging; on the other hand, more advanced methods demand specified knowledge of the designer and often also require the use of advanced numerical tools. The simplest and most frequently used method includes the McKee formula [11], which is based on the utmost popular parameters of corrugated board, that is, column crush resistance, commonly known as ECT, from the name of the test—edge crush test, and the thickness of the cardboard or bending stiffness, as well as the dimensions of the box in the base. The ECT is standardized—four methods, depending on the samples’ shapes, have been established, that is, the edge clamping method [12], the neck-down method [13], the rectangular test specimen method [13,14,15], and the edge-reinforced method [16,17]. In the literature, one can also find many modifications/adaptations aimed at increasing the effectiveness of McKee formulas provided by, for example, Frank [7], Garbowski et al. [18], Maltenfort [19], Allerby et al. [20], Schrampfer et al. [21], Kawanishi [22], and Batelka [23] and their adaptation to other, more demanding, packaging cases [8,9,10].
Nowadays, a huge emphasis is placed on ecology and optimization of the resources use (the so-called zero waste, fit-to-product, box-on-demand, etc.), therefore, the correct modeling of packaging is taking on increasing importance. The optimal design of the corrugated cardboard packages was examined by Mrówczyński et al. [24,25]. The use of recycled fibers in the production of paper and cardboard causes the fibers to become shorter and more brittle after many recovery cycles, and their ability to create a stable network of connections with other fibers keeps getting smaller over time. The utilization of starch and other agents, to halt the decline in the ability of fibers to bond to one another, causes the significant change in mechanical properties of the paper. The produced paper, although more ecological, becomes a much more heterogeneous material with a much larger number and size of initial imperfections. All geometrical and material imperfections affect the strength of paper and corrugated board; however, this is not yet a well-recognized process by engineers and scientists. Only a few works dealing with this problem can be found in the literature. Garbowski and Knitter-Piątkowska [26] took into account geometric imperfections in the analytical calculations of bending stiffness in machine direction (MD). This effect was previously investigated experimentally by Czechowski et al. [27].
Numerical techniques, based on the finite element method (FEM) [28,29] allow for accurate estimation of the load-bearing capacity of complex packaging structures [30], the determining of mechanical properties of cardboard with creases [31,32,33,34,35,36], as well as strength estimations of corrugated board packages [37,38,39,40]. FEM has also been applied to examine the torsional and transversal stiffness of orthotropic paper materials [41,42], and the bending stiffness [27,43] and buckling or post-buckling phenomena [44] of cardboard. Since the paper materials are anisotropic and the cardboard’s structure is layered, numerical simulations are demanding because the material parameters of each layer need to be known. By utilizing the method called homogenization, one can preserve the precision of the results with significant savings in computation time. In the analytical homogenization, the equations of the classical theory of strength of materials or the classical theory of laminates are used [45]. Numerical homogenization makes use of FEM, where first a numerical model of a representative volume element (RVE) is created [46]. This approach fits perfectly with cardboard issues where homogenization can be accomplished in two ways, namely homogenization to one layer or homogenization of fluting to the inner layer of the laminate. Due to its practical importance and aforementioned advantages, this procedure has been intensively developed and utilized in recent years, as reported in the literature [47,48,49,50,51,52,53,54,55,56,57,58].
FEM models can be created based on laboratory tests from which the material parameters of corrugated cardboard are obtained. Therefore, their correct definition becomes a priority. In order to comprehensively obtain all the necessary material parameters, laboratory devices can be used [59]. Through this innovative system for predicting the compressive strength of corrugated cardboard packaging, a material database is created on the basis of the bending, shear, twisting and ECT of the corrugated board samples. The BCT index is calculated for the selected project and the previously calculated material parameters of the corrugated board. However, this approach is based on research on corrugated board, which is limited only to the range of board already produced. The need to improve the corrugated board itself, driven mainly by the dynamically developing e-commerce industry, makes manufacturers outdo each other in optimizing their products by changing the composition of the produced cardboard.
Modeling from paper-to-packaging requires from the designer the knowledge of the material parameters of the constituent papers in the composition of the corrugated board. These parameters are commonly obtained by means of short-span compression tests (SCT), and by tensile tests of cardboard samples. These tests are used to determine the compressive strength and tensile strength, as well as the stiffness modulus in different directions with respect to the direction of the fibers, that is, machine direction (MD) or cross direction (CD), together with any chosen direction rotated relative to the MD by any angle. To obtain the data from the exterior surface of the specimen during the experiment, a video extensometry can be applied. This technique, comparable to digital image correlation (DIC), has been utilized by Garbowski et al. [60,61]. Such a method is distinguished by high accuracy of data capture, and is very productive in the field of experimental mechanics [42,62,63,64,65]. Correct determination of the corrugated board mechanical parameters and the knowledge of the corrugated layers geometry allows to build a model of the corrugated board. The simplest models for ECT estimation are based on empirical dependencies in which the load capacities of all papers are summed, taking into account the development factor of the fluting, namely, the so-called take-up factor. The sum of these load capacities is then scaled with the fit factor, the value of which is usually set at 0.7–0.8. In the literature, one can also find ECT models based on analytical [65,66], analytical-numerical [59,60] or experimental methods [67].
The review of several aspects of ECT testing methods, namely specimen height, test duration, and fixture-clamping effects, was drawn up by Popil [68]. It was demonstrated that the combination of the chosen testing procedures with the specific structural and strength characteristics of the cardboard being examined affects ECT values, and proved that the measurement of compression strength is sensitive to: the type of method, sample preparation, the effect of crushing, as well as the effect of test duration. In the test T 839 [12], a specimen with the dimensions 50 mm × 50 mm, irrespective of the board type, was held at its ends in the clamping fixture, which exerts pressure through its spring-loaded jaws that cover approximately two-thirds of the height of the test specimen preventing any bending of the board. According to Frank [69], T 839 may be inappropriate for some lightweight boards with relatively big calipers. Unsupported samples with waxed edges were examined in the test T 811 [16]. Their height was specifically defined for the common flute types A, B, and C. During this test, on both sides of the specimen, supporting guide blocks were positioned to ensure perpendicularity or vertical adjustment, and they were removed once the force reached the value of 22 N. Because the samples are less restrained than in the T 839 method, T 811 is expected to give lower ECT values in many cases of examining the corrugated board. Another common compression test is T 838 [13], in which the specimens with narrowed width are in the middle, thus concentrating the stress there. This technique eliminates the drawbacks of using wax to strengthen the edges in method T 811 and the consequences of exerting pressure with jaws in method T 839. In comparison to other methods, it might give the highest ECT value for lightweight boards. The influence of fluting grammage on the mechano-sorption creep of cardboard was investigated by Popil and Hojjatie [70] while using various ECT methods. Regarding the effect of sample height on ECT, the studies [68,71] have shown that ECT values decrease as specimen height increases due to the sample buckling rather than sample compression. The paper is a viscoelastic material, thus the rate of compression has an influence on ECT. The test methods indicate a compression platen running speed of 12.5 + 0.25 mm/min, however data for ECT can be received at different strain rates, and for this reason, to compare data sets, the quantification of the time relation of ECT is necessary. In general, strength properties decrease by about 7.5% per decade of strain rate change [68]. Edge crush test (ECT) measurements using the aforementioned tests on boxes crushed in varying amounts during the manufacturing process was discussed by Frank and Cash [72].
In this paper, attention is focused on the empirical ECT model based on the SCT parameters of the constituent papers of various three- and five-layer corrugated board compositions. For this purpose, SCT tests of the component papers were performed in the main directions of the material orthotropy, that is, in the MD, CD directions and in the direction rotated by 45 degrees in relation to the MD. Then, a series of ECT tests were performed on the cardboard built on these papers, loaded at different angles to the wave direction. The results of laboratory tests were used to construct an easy-to-calibrate and reliable ECT model of corrugated board. The empirical-analytical model proposed here is based on easy-to-obtain paper parameters and is much more practical than detailed models that require specific modeling knowledge, presented in our previous papers. The model is based on the theory of bearing capacity of a post-buckled plate, in which, due to the buckling of the central part, only a portion of the section is bearing the load, the so-called effective cross-section, similar to the theory used in thin-walled steel cross-sections. Due to the need to avoid laborious calculations of the critical force, the model has replaced the critical force with an empirical model that uses only the grammage of the paper (instead of its thickness) and the width of the analyzed section. The model is calibrated with only one parameter, which has a constant value for all single-walled boards and another constant value for all double-walled boards. The adopted simplification allows for an easy and very practical estimation of the ECT value based on the SCT and the grammage of the constituent papers. In addition, for comparative purposes, a numerical model of corrugated board was built on the basis of data obtained from mechanical laboratory tests of individual papers. The model was verified by laboratory tests of corrugated cardboard, and then the influence of various types and sizes of imperfections on the results of numerical analyses was checked. The observations clearly show that both material and geometric imperfections have a significant impact on the edge crush resistance of the corrugated board.

2. Materials and Methods

2.1. Laboratory Tests

In order to build the correct analytical or numerical model of the corrugated board, it was necessary to perform a series of laboratory tests of all component papers and obtain the necessary material parameters. Aiming to correctly define the computational models, the papers were first subjected to tensile tests (see Figure 1a) in three main directions, namely MD (machine direction), CD (cross direction), and at 45 degrees with respect to MD (see Figure 1b), and short-span compression tests (SCT) in the CD (i.e., direction transverse to the direction of the fibers). To find the tensile parameters of paper samples, a Testometric laboratory device was used in accordance with the ISO 9073-3 standard [73].
The parameters of the paper in compression were identified only in one direction, that is, in the CD, in which the corrugated board usually works. For this purpose, a Short Span Compression Tester from TMI, model 17–36 (Messmer Büchel-Industrial Physics, LLC, Veenendaal, The Netherlands), was used in accordance with the ISO 9895 standard [74]. The SCT device is shown in Figure 2a.
The reference values of the edge crush test (ECT) for various corrugated board samples were measured according to the PN-EN ISO 3037:2013-12 [15] standard. In such a test, the sample is 100 mm long and 25 mm high. The ECT tests were performed on the device from FEMAT (Poznań, Poland), model ECT-10-21 (see Figure 2b). As can be seen in Figure 2b, the samples were each time supported by metal blocks that prevented the samples from buckling or tilting during the test.
It is worth mentioning that in corrugated cardboard, all papers are arranged in such a way that the direction of the fibers coincides with the direction of the wave. This means that in the ECT test, all layers are loaded in the CD. Therefore, in the presented procedure, the compressive strength of papers was tested only in the CD direction. The tensile tests were only necessary to determine the orthogonal stiffnesses in the linear-elastic orthotropic model.
All tests, both paper and board tests, were carried out under standard laboratory conditions, that is, 22 °C and 50% humidity. In all laboratory tests, eight material samples were prepared (samples of paper cut in selected directions and samples of corrugated board). As the tests were performed in accredited laboratories, the statistical analysis of the obtained results was omitted in this paper and, therefore, all presented results are the mean values.

2.2. Corrugated Board

In order to develop a correct computational model of corrugated board in the ECT test, a representative group of different cardboards was selected. In total, six popular compositions of three-layer and five-layer corrugated board were chosen for the tests, the main parameters of which are listed in Table 1.
Table 2 compiles the grammage values of the tested cardboards and the measured cross-sectional heights. In order to systematize the presented results, the cardboard ID was adopted, consisting of a symbol describing the type of wave and the basis weight of individual corrugated cardboard samples.
The corrugated boards presented in Table 2 are widely used for the production of packaging and, at the same time, constitute a wide spectrum of its different types, from single-walled with a height of almost 1.6 mm to double-walled with a height of over 6.7 mm.

2.3. Component Papers

All paper samples were examined in two main orthotropic directions, that is, in the MD and CD, and at an angle of 45 degrees to the MD. Markings and grammages of individual component papers for all cardboard samples are listed in Table 3.

2.4. Numerical Model

In order to model the ECT test correctly and comprehensively, the finite element method was used. On the basis of the measured stiffness and compressive and tensile strength in the two main orthotropic directions, as well as in the direction rotated by 45 degrees, the remaining material parameters were determined. The modulus of shear stiffness in the plane was determined from the following formula:
1 E x α = cos 2 α E x + 1 G x y 2 ν x y E x sin 2 α cos 2 α + sin 4 α E y .
Because in the tests carried out, the examinations were performed in the direction of 45 degrees relative to the x direction, and the direction x = MD and y = CD, Formula (1) therefore takes the form:
G x y = 4 E 45 1 E x 1 E y + 2 ν x y E x   1 .
In the above formulas, E x is the stiffness modulus in the MD direction, E y is the stiffness modulus in the CD direction, E 45 is the stiffness in the direction rotated 45 degrees from the MD direction, and ν x y is the Poisson’s ratio in the xy plane, which can be determined by the empirical formula proposed by Baum [21]:
ν x y = 0.293 E y E x .
Due to the special case of loading the corrugated board sample at the edges along the CD direction, which is observed in the ECT test, only the SCT parameter in the CD direction was used to model the inelastic behavior. As a simplification, it is assumed that in the ECT test, the compressive strength of the paper in other directions is not activated, nor is the tensile strength. Therefore, the constitutive model consists of the orthotropic linear elastic part and the nonlinear isotropic inelastic part. As was already indicated, the numerical model was built while applying the finite element (FE) method. A shell element with four nodes and four Gaussian points with an approximate dimension of one side of the element of about 0.35 mm was used for the analysis. Calculations were made in the ABAQUS FEA 2021 commercial software.
In the numerical analysis of thin shell elements that can undergo large displacements or rotations, apart from material non-linearity, geometric non-linearity should also be taken into account. This allows for correct modelling when thin layers of paper start to buckle. Unfortunately, numerical models without initial imperfections are often characterized by a lack of convergence, thus geometric or material imperfections should be taken into account in the modeling. In order to account for imperfections in the FE model of the ECT sample, the geometrical imperfections were introduced by changing the position of nodes in the model, so that the shape of the geometric imperfection corresponded to the lowest eigenvectors buckling model (see Figure 3).
The numerical model shown in Figure 3 is only a small part of a whole model and is intended to show the shape of buckling mode applied to account for imperfections. The full model covers the entire sample, that is, the model with dimensions of 100 mm × 25 mm. All geometric features of the numerical models are presented in Table 1. The numerical model itself, the description of the boundary conditions, and the loading condition and convergence studies were carried out and shown in our previous work [60,61].

2.5. Empirical Model

In the definition of the simplest empirical ECT model of corrugated board, the SCT (in CD) of the component papers and the take-up factor of corrugated layers are used. The model can be defined as follows:
E C T = A i = 1 n S C T i α i .
In this model, the SCT values of individual layers are summed up and multiplied by the take-up factor, which for flat layers equals 1. The fitting factor A is the only parameter that is used to correlate the result obtained with the model and the actual result from the laboratory test. It is very impractical, because in this way we adjust the test results for both single and double-wall corrugated board with one parameter only, which appears to be different for almost any cardboard configuration.

2.6. Analytical-Empirical Model

In order to better represent the actual behavior of a cardboard sample, where imperfections often occur, an analytical–empirical model can be used. It is assumed that the model takes into account both the SCT parameter, which describes the compressive strength of individual paper segments (both liners between the tops of the waves and curvilinear segments of fluting) and the critical force, which describes the initiation of buckling. It does not require the correlation of the empirical parameter A that appears in Equation (4). Such a model can be written, for example, in the following form:
E C T = i = 1 n S C T i α i γ i ,
where the γ i parameter is a factor reducing the static load capacity on compression of individual layers. This coefficient takes into account the buckling of the layers and depends on the geometry of the flute and the grammage of the given component paper. It would be more correct to base this parameter on paper thickness, but this parameter is rarely determined in the laboratory and does not appear in paper specifications supplied by paper mills. The formula describing the parameter γ i can be written as:
γ i = S C T i 1 / 2 H a g i b i 1 / 2 1 ,
where a is a parameter that should be selected once and for all on the basis of experimental observations (empirical adjustment), H is the height of the cross-section, g i is the basis weight of a given layer, b i is the width of individual segments in the cross-section of corrugated board (details in Figure 4).
In the case of corrugated board with a double wall, a better solution is to omit the parameter H from the Equation (6), which then takes the form:
γ i = S C T i 1 / 2 1 a * g i b i 1 / 2 1 ,
where a * is again a parameter to be determined on the basis of experimental observations, g i is the basis weight of a given paper, while b i is the width of individual segments in the cross-section of corrugated board (details in Figure 5). The value of the width b 3 (inner flat layer) was assumed a priori as 0.7 of the smaller of the two wavelengths, that is:
b 3 = min P 1 , P 2 × 0.7 ,
which, in this case, leads to: b 3 = 0.7 P 1 .
It is worth noting that the parameter γ i is the square root of the ratio between the critical force and the compressive strength of the selected corrugated board segment. Instead of tedious analytical [75,76] or numerical [77,78] calculations of the critical load of the analyzed segment, its empirical form is assumed here. This factor somewhat reduces the width of each segment that actually carries the load during the edge crush test.

3. Results

3.1. Laboratory Tests Results

Table 4 summarizes the results of the SCT in CD and tensile stiffness tests in three directions (MD, CD and 45 degrees) obtained for all constituent papers of the six selected corrugated boards: (a) B-410; (b) C-590; (c) E-480; (d) BC-790; (e) BE-600; (f) BE-590. Each paper was tested 3–4 times, and when the results were widely dispersed, additional tests were performed. Average values are presented in the table.

3.2. The ECT Estimations by Numerical and Analytical-Empirical and Empirical Models

Table 5 presents the empirically determined coefficients presented in Equations (4), (6) and (7). These coefficients were determined using least squares minimization. It is worth noting that the coefficient A was selected individually for each quality of corrugated board and then averaged (columns 3 and 4 in Table 5). The coefficients a and a * were computed individually for single-walled and double-walled cardboards.
The simplest empirical model described by Equation (4) uses the coefficient A , which can be selected for each cardboard quality separately (column 2 of Table 5), and it may also be selected separately for each type of corrugated board (column 3 of Table 5) or commonly for all cardboards (column 4 of Table 5).
Finally, five models listed below were used for ECT estimation:
  • Model 1: proposed analytical–empirical model, Equation (5);
  • Model 2: numerical model FEM1 with imperfection 2%;
  • Model 3: numerical model FEM2 with imperfection 9%;
  • Model 4: empirical model described by Equation (4) with coefficient from Table 5 (column 3);
  • Model 5: empirical model described by Equation (4) with coefficient from Table 5 (column 4).
Table 6 presents the estimation results using the proposed Model 1 (see Equation (5)) and two numerical models—Model 2 and Model 3—which differ only in the assumed scale of imperfections.
Table 7 summarizes the estimation error obtained while using different models. The abbreviations M1 to M5 denote the model numbers that have already been described.
Figure 6 shows the absolute estimation errors obtained while utilizing different models. It is worth noting that the estimations obtained using the models based on Equation (4) clearly deviate from the estimations which make use of the other models.

4. Discussion

The test results presented in Table 4 made it possible to build a material model, which was then used in a numerical model based on the finite element method for ECT estimation. Both numerical models used for ECT calculations differ only in the size of imperfections that must be included in the model to obtain correct results. This is due to the fact that in the real sample, imperfections are built in a priori, therefore numerical models should also have initial geometrical imperfections. The problem is that it is not known a priori where exactly and how big these imperfections are. One way to incorporate the initial imperfections is to determine the first buckling mode and scale the obtained displacements to obtain a new geometry that takes these displacements into account. Here, two different magnitudes of these displacements have been chosen: 2% and 9%.
As can be seen in Table 6 (columns 4–5) and Table 7 (columns 3–4), the estimation results obtained while using numerical model M2 are more correct than those obtained using the model M3. In the M2 model, imperfections of 2% of the displacement values obtained from the buckling analysis for the first lowest mode were used. In the M3 model, imperfections at the level of 9% of the value of displacements on each wall of the corrugated board cross-section were applied. Unfortunately, the actual imperfections are not known, however, the obtained estimations show that they amount to about 3–5% of the buckling displacements. This is due to the fact that the values obtained using the M2 model were higher than the reference values in each analyzed case, while the M3 model underestimated the ECT value in each case.
Models M4 and M5 based on Equation (4), although quite effective in the case of single-walled cardboard, were completely ineffective in the case of double-walled corrugated board. This is shown very clearly by the results presented in Table 7 and Figure 6. Even the determination of the correlation coefficient A separately for two types of board did not help—the stimulation error for double-walled board still oscillated between +/−13%. This is, of course, related to the very simple form of the model (see Equation (4)).
By contrast, model M1 proposed here behaved very correctly in each case of the analyzed corrugated board. The results obtained using the M1 model are burdened with an estimation error not exceeding +/−5%. The numerical models M2 and M3 are characterized by similar accuracy. For the correct construction of the M1 model, the following ingredients are necessary: (a) the compressive strength of the constituent papers (SCT), (b) their grammage and (c) the geometrical characteristics of the corrugated layers. Correct numerical modeling requires additionally: (d) paper stiffness in each orthotropic direction (i.e., in MD, CD and 45 deg) and, of course, (e) specialized software, as well as (f) necessary expert knowledge in computer modeling. Therefore, the proposed analytical–empirical model allows to significantly speed up the estimation of ECT without losing the accuracy of the obtained results compared to numerical models. In addition, this model is much more precise than the currently used simplified models.

5. Conclusions

The paper presented an analytical–empirical model for the estimation of corrugated board resistance to edge crushing. On the basis of the obtained results, also comparing the accuracy of the estimation with currently applicable laboratory models and with advanced numerical models based on the finite element method, it is concluded that the model is: (1) as precise as numerical models; (2) much simpler to use than numerical models; (3) much more precise than simplified models; and (4) it requires only a few readily available parameters. As shown in the work on six different grades of corrugated board, simplified modeling provides a very practical tool that, in addition to being accurate, is also very easy to calibrate because the model is based on commonly known and easily obtained paper and corrugated board parameters. The common use of the proposed here model in the laboratory practice of companies producing corrugated board will not only significantly improve the process of selecting the appropriate inputs, but also through high precision, achieve significant savings of the raw material.

Author Contributions

Conceptualization, T.G.; methodology, T.G.; software, T.G.; validation, T.G., A.K.-P. and P.W.; formal analysis, T.G. and P.W.; investigation, T.G. and A.K.-P.; resources, P.W.; data curation, P.W.; writing—original draft preparation, T.G. and A.K.-P.; writing—review and editing, T.G. and A.K.-P.; visualization, T.G. and P.W.; supervision, T.G.; project administration, T.G. and P.W.; funding acquisition, T.G., A.K.-P. and P.W. All authors have read and agreed to the published version of the manuscript.

Funding

The APC was partially funded by the National Center for Research and Development, Poland, grant number POIR.01.01.01-00-1006/19.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This research was funded by the National Center for Research and Development, Poland, grant at Schumacher Packaging Sp. z o. o., grant number POIR.01.01.01-00-1006/19.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Stott, R.A. Compression and stacking strength of corrugated fibreboard containers. Appita Technol. Innov. Manuf. Environ. 2017, 70, 76–82. [Google Scholar]
  2. Junli, W.; Quancheng, Z. Effect of moisture content of corrugated box on mechanical properties. J. Lanzhou Jiaotong Univ. 2006, 25, 134–136. [Google Scholar]
  3. Zhang, Y.-L.; Chen, J.; Wu, Y.; Sun, J. Analysis of hazard factors of the use of corrugated carton in packaging low-temperature yogurt during logistics. Procedia Environ. Sci. 2011, 10, 968–973. [Google Scholar] [CrossRef] [Green Version]
  4. Gallo, J.; Cortés, F.; Alberdi, E.; Goti, A. Mechanical behavior modeling of containers and octabins made of corrugated cardboard subjected to vertical stacking loads. Materials 2021, 14, 2392. [Google Scholar] [CrossRef]
  5. Whitsitt, W.J.; Gander, J.W.; McKee, R.C. Effect of Box Dimensions and Combined Board Creep Life on Box Creep Life; Institute of Paper Chemistry: Appleton, WI, USA, 1968. [Google Scholar]
  6. Urbanik, T.J.; Frank, B. Box compression analysis of world-wide data spanning 46 years. Wood Fiber Sci. 2006, 38, 399–416. [Google Scholar]
  7. Frank, B. Corrugated Box Compression—A Literature Survey. Packag. Technol. Sci. 2013, 27, 105–128. [Google Scholar] [CrossRef]
  8. Garbowski, T.; Gajewski, T.; Grabski, J.K. Estimation of the compressive strength of corrugated cardboard boxes with various perforations. Energies 2021, 14, 1095. [Google Scholar] [CrossRef]
  9. Garbowski, T.; Gajewski, T.; Grabski, J.K. Estimation of the compressive strength of corrugated cardboard boxes with various openings. Energies 2021, 14, 155. [Google Scholar] [CrossRef]
  10. Mrówczyński, D.; Garbowski, T.; Knitter-Piątkowska, A. Estimation of the compressive strength of corrugated cardboard boxes with shifted creases on the flaps. Materials 2021, 14, 5181. [Google Scholar] [CrossRef]
  11. McKee, R.C.; Gander, J.W.; Wachuta, J.R. Compression strength formula for corrugated boxes. Paperboard Packag. 1963, 48, 149–159. [Google Scholar]
  12. TAPPI T 839 om-12; Edge Compression Test for Strength of Corrugated Fiberboard Using the Clamp Method (Short Column Test). TAPPI: Peachtree Corners, GA, USA, 2009.
  13. TAPPI T 838 cm-12; Edge Crush Test Using Neckdown. TAPPI: Peachtree Corners, GA, USA, 2009.
  14. FEFCO NO.8; Edgewise Crush Resistance of Corrugated Fiberboard. FEFCO: Brussels, Belgium, 1997.
  15. ISO 3037:2013; Corrugated Fibreboard—Determination of Edgewise Crush Resistance (Unwaxed Edge Method). ISO: Geneva, Switzerland, 2013.
  16. TAPPI T 811 om-11; Edgewise Compressive Strength of Corrugated Fibreboard (Short Column Test). TAPPI: Peachtree Corners, GA, USA, 2009.
  17. ISO 13821:2002; Corrugated Fibreboard—Determination of Edgewise Crush Resistance—Waxed Edge Method. ISO: Geneva, Switzerland, 2002.
  18. Garbowski, T.; Niziałek-Łukawska, M.; Kuca, M. Analytical verification of popular McKee’s formula. Pol. Pap. Rev. 2019, 75, 767–773. (In Polish) [Google Scholar]
  19. Maltenfort, G. Compression strength of corrugated containers. Fibre Contain. 1956, 41, 106–121. [Google Scholar]
  20. Allerby, I.M.; Laing, G.N.; Cardwell, R.D. Compressive strength—From components to corrugated containers. Appita Conf. Notes 1985, 1–11. [Google Scholar]
  21. Schrampfer, K.E.; Whitsitt, W.J.; Baum, G.A. Combined Board Edge Crush (ECT) Technology; Institute of Paper Chemistry: Appleton, WI, USA, 1987. [Google Scholar]
  22. Kawanishi, K. Estimation of the compression strength of corrugated fibreboard boxes and its application to box design using a personal computer. Packag. Technol. Sci. 1989, 2, 29–39. [Google Scholar] [CrossRef]
  23. Batelka, J.J.; Smith, C.N. Package Compression Model; Institute of Paper Science and Technology, Georgia Institute of Technology: Atlanta, GA, USA, 1993. [Google Scholar]
  24. Mrówczyński, D.; Knitter-Piątkowska, A.; Garbowski, T. Non-Local Sensitivity Analysis and Numerical Homogenization in Optimal Design of Single-Wall Corrugated Board Packaging. Materials 2022, 15, 720. [Google Scholar] [CrossRef] [PubMed]
  25. Mrówczyński, D.; Knitter-Piątkowska, A.; Garbowski, T. Optimal Design of Double-Walled Corrugated Board Packaging. Materials 2022, 15, 2149. [Google Scholar] [CrossRef] [PubMed]
  26. Garbowski, T.; Knitter-Piątkowska, A. Analytical determination of the bending stiffness of a five-layer corrugated cardboard with imperfections. Materials 2022, 15, 663. [Google Scholar] [CrossRef]
  27. Czechowski, L.; Kmita-Fudalej, G.; Szewczyk, W.; Gralewski, J.; Bieńkowska, M. Numerical and experimental study of five-layer non-symmetrical paperboard panel stiffness. Materials 2021, 14, 7453. [Google Scholar] [CrossRef]
  28. Zienkiewicz, O.C.; Taylor, R.L.; Zhu, J.Z. The Finite Element Method: Its Basis and Fundamentals, 1st ed.; Butterworth-Heinemann: Oxford, UK, 2013. [Google Scholar]
  29. De Borst, R.; Crisfield, M.A.; Remmers, J.J.C.; Verhoosel, C.V. Nonlinear Finite Element Analysis of Solids and Structures, 2nd ed.; Wiley: Chichester, UK, 2012. [Google Scholar]
  30. Fadiji, T.; Coetzee, C.J.; 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]
  31. Domaneschi, M.; Perego, U.; Borgqvist, E.; Borsari, R. An industry-oriented strategy for the finite element simulation of paperboard creasing and folding. Packag. Technol. Sci. 2017, 30, 269–294. [Google Scholar] [CrossRef]
  32. Awais, M.; Tanninen, P.; Leppänen, T.; Matthews, S.; Sorvari, J.; Varis, J.; Backfol, K. A computational and experimental analysis of crease behavior in press forming process. Procedia Manuf. 2018, 17, 835–842. [Google Scholar] [CrossRef]
  33. Thakkar, B.K.; Gooren, L.G.J.; Peerlings, R.H.J.; Geers, M.G.D. Experimental and numerical investigation of creasing in corrugated paperboard. Philos. Mag. 2008, 88, 3299–3310. [Google Scholar] [CrossRef] [Green Version]
  34. Beex, L.A.A.; Peerlings, R.H.J. An experimental and computational study of laminated paperboard creasing and folding. Int. J. Solids Struct. 2009, 46, 4192–4207. [Google Scholar] [CrossRef] [Green Version]
  35. Giampieri, A.; Perego, U.; Borsari, R. A constitutive model for the mechanical response of the folding of creased paperboard. Int. J. Solids Struct. 2011, 48, 2275–2287. [Google Scholar] [CrossRef] [Green Version]
  36. Leminen, V.; Tanninen, P.; Pesonen, A.; Varis, J. Effect of mechanical perforation on the press-forming process of paperboard. Procedia Manuf. 2019, 38, 1402–1408. [Google Scholar] [CrossRef]
  37. Garbowski, T.; Jarmuszczak, M. Numerical strength estimate of corrugated board packages. Part 1. Theoretical assumptions in numerical modeling of paperboard packages. Pol. Pap. Rev. 2014, 70, 219–222. (In Polish) [Google Scholar]
  38. Garbowski, T.; Jarmuszczak, M. Numerical strength estimate of corrugated board packages. Part 2. Experimental tests and numerical analysis of paperboard packages. Pol. Pap. Rev. 2014, 70, 277–281. (In Polish) [Google Scholar]
  39. Park, J.; Chang, S.; Jung, H.M. Numerical prediction of equivalent mechanical properties of corrugated paperboard by 3D finite element analysis. Appl. Sci. 2020, 10, 7973. [Google Scholar] [CrossRef]
  40. Park, J.; Park, M.; Choi, D.S.; Jung, H.M.; Hwang, S.W. Finite element-based simulation for edgewise compression behavior of corrugated paperboard for packing of agricultural products. Appl. Sci. 2020, 10, 6716. [Google Scholar] [CrossRef]
  41. Nordstrand, T.; Carlsson, L. Evaluation of transverse shear stiffness of structural core sandwich plates. Compos. Struct. 1997, 37, 145–153. [Google Scholar] [CrossRef]
  42. Słonina, M.; Dziurka, D.; Smardzewski, J. Experimental research and numerical analysis of the elastic properties of paper cell cores before and after impregnation. Materials 2020, 13, 2058. [Google Scholar] [CrossRef] [PubMed]
  43. Jamsari, M.A.; Kueh, C.; Gray-Stuart, E.M.; Dahm, K.; Bronlund, J.E. Experimental and numerical performance of corrugated fibreboard at different orientations under four-point bending test. Packag. Technol. Sci. 2019, 32, 555–565. [Google Scholar] [CrossRef]
  44. Urbanik, T.J.; Saliklis, E.P. Finite element corroboration of buckling phenomena observed in corrugated boxes. Wood Fiber Sci. 2003, 35, 322–333. [Google Scholar]
  45. Allaoui, S.; Benzeggagh, M.L.; Aboura, Z.; Talbi, N. Elastic behaviour of corrugated cardboard: Experiments and modeling. Compos. Struct. 2004, 63, 53–62. [Google Scholar]
  46. Biancolini, M.E. Evaluation of equivalent stiffness properties of corrugated board. Compos. Struct. 2005, 69, 322–328. [Google Scholar] [CrossRef]
  47. Garbowski, T.; Jarmuszczak, M. Homogenization of corrugated paperboard. Part 1. Analytical homogenization. Pol. Pap. Rev. 2014, 70, 345–349. (In Polish) [Google Scholar]
  48. Garbowski, T.; Jarmuszczak, M. Homogenization of corrugated paperboard. Part 2. Numerical homogenization. Pol. Pap. Rev. 2014, 70, 390–394. (In Polish) [Google Scholar]
  49. Garbowski, T.; Gajewski, T. Determination of transverse shear stiffness of sandwich panels with a corrugated core by numerical homogenization. Materials 2021, 14, 1976. [Google Scholar] [CrossRef]
  50. 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]
  51. Mrówczyński, D.; Knitter-Piątkowska, A.; Garbowski, T. Numerical Homogenization of Single-Walled Corrugated Board with Imperfections. Appl. Sci. 2022, 12, 9632. [Google Scholar] [CrossRef]
  52. Ramírez-Torres, A.; Penta, R.; Rodríguez-Ramos, R.; Merodio, J.; Sabina, F.J.; Bravo-Castillero, J.; Guinovart-Díaz, R.; Preziosi, L.; Grillo, A. Three scales asymptotic homogenization and its application to layered hierarchical hard tissues. Int. J. Solids Struct. 2018, 130–131, 190–198. [Google Scholar] [CrossRef]
  53. Ramírez-Torres, A.; Di Stefano, S.; Grillo, A.; Rodríguez-Ramos, R.; Merodio, J.; Penta, R. An asymptotic homogenization approach to the microstructural evolution of heterogeneous media. Int. J. Non-Linear Mech. 2018, 106, 245–257. [Google Scholar] [CrossRef] [Green Version]
  54. Hohe, J. A direct homogenization approach for determination of the stiffness matrix for microheterogeneous plates with application to sandwich panels. Compos. Part B Eng. 2003, 34, 615–626. [Google Scholar] [CrossRef]
  55. Buannic, N.; Cartraud, P.; Quesnel, T. Homogenization of corrugated core sandwich panels. Compos. Struct. 2003, 59, 299–312. [Google Scholar] [CrossRef] [Green Version]
  56. Abbès, B.; Guo, Y.Q. Analytic homogenization for torsion of orthotropic sandwich plates. Compos. Struct. 2010, 92, 699–706. [Google Scholar] [CrossRef]
  57. Marin, G.; Srinivasa, P.; Nygårds, M.; Östlund, S. Experimental and finite element simulated box compression tests on paperboard packages at different moisture levels. Packag. Technol. Sci. 2021, 34, 229–243. [Google Scholar] [CrossRef]
  58. 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]
  59. FEMat BSE Systems. Available online: http://fematsystems.pl/bse-system_en (accessed on 17 October 2022).
  60. Garbowski, T.; Grabski, J.K.; Marek, A. Full-field measurements in the edge crush test of a corrugated board—Analytical and numerical predictive models. Materials 2021, 14, 2840. [Google Scholar] [CrossRef]
  61. Garbowski, T.; Knitter-Piątkowska, A.; Marek, A. New edge crush test configuration enhanced with full-field strain measurements. Materials 2021, 14, 5768. [Google Scholar] [CrossRef]
  62. Hägglund, R.; Åslund, P.E.; Carlsson, L.A.; Isaksson, P. Measuring thickness changes of edgewise compression loaded corrugated board panels using digital image correlation. J. Sandw. Struct. Mater. 2010, 14, 75–94. [Google Scholar] [CrossRef]
  63. Viguié, J.; Dumont, P.J.J. Analytical post-buckling model of corrugated board panels using digital image correlation measurements. Comp. Struct. 2013, 101, 243–254. [Google Scholar] [CrossRef]
  64. Fadiji, T.; Coetzee, C.J.; Opara, U.L. Evaluating the displacement field of paperboard packages subjected to compression loading using digital image correlation (DIC). Food Bioprod. Process. 2020, 123, 60–71. [Google Scholar] [CrossRef]
  65. Maier, G.; Bolzon, G.; Buljak, V.; Garbowski, T.; Miller, B. Synergic combinations of computational methods and experiments for structural diagnoses. In Computer Methods in Mechanics; Kuczma, M., Wilmanski, K., Eds.; Advanced Structured Materials; Springer: Berlin/Heidelberg, Germany, 2010; Volume 1, pp. 453–476. [Google Scholar]
  66. Garbowski, T.; Imbierowicz, R. Sensitivity analysis of edge crush test. Przegląd Pap. 2014, 70, 559–564. (In Polish) [Google Scholar]
  67. 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. 2009, 71, 227–233. [Google Scholar] [CrossRef]
  68. Popil, R.E. Overview of recent studies at IPST on corrugated board edge compression strength: Testing methods and effects of interflute buckling. BioResources 2012, 7, 2553–2581. [Google Scholar] [CrossRef] [Green Version]
  69. Frank, B. Revisiting clamped ECT. Corrugat. Int. 2007, Summer, 7–11. [Google Scholar]
  70. Popil, R.E.; Hojjatie, B. Effects of component properties and orientation on corrugated container endurance. Packag. Technol. Sci. 2010, 23, 189–202. [Google Scholar] [CrossRef]
  71. Wilson, C.J.; Frank, B. An evaluation of ECT sample height for small flute board grades and Box Manufacturer’s Certification compliance. Tappi J. 2009, 8, 24–28. [Google Scholar] [CrossRef]
  72. Frank, B.; Cash, D. Edge crush testing methods and box compression modeling. Tappi J. 2022, 21, 418. [Google Scholar] [CrossRef]
  73. ISO 9073-3:1989; Textiles—Test methods for nonwovens—Part 3: Determination of tensile strength and elongation. ISO: Geneva, Switzerland, 1989.
  74. ISO 9895:2008; Paper and board—Compressive strength—Short-span test. ISO: Geneva, Switzerland, 2008.
  75. 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]
  76. Nyman, U.; Gustafsson, P.J. Material and structural failure criterion of corrugated board facings. Compos. Struct. 2000, 50, 79–83. [Google Scholar] [CrossRef]
  77. Daxner, T.; Flatscher, T.; Rammerstorfer, F.G. 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; pp. 349–358. [Google Scholar]
  78. Ma, Y.; Gong, Z.; Liang Zhao, L.; Han, Y.; Wang, Y. Finite Element Analysis of Buckling of Corrugated Fiberboard. Open Mech. Eng. J. 2014, 8, 257–263. [Google Scholar] [CrossRef]
Figure 1. Paper tensile stiffness test: (a) laboratory device; (b) samples cut in the directions: MD, CD and 45 degrees.
Figure 1. Paper tensile stiffness test: (a) laboratory device; (b) samples cut in the directions: MD, CD and 45 degrees.
Materials 16 00458 g001
Figure 2. Laboratory testing machines: (a) short-span compression tester; (b) edge crush tester.
Figure 2. Laboratory testing machines: (a) short-span compression tester; (b) edge crush tester.
Materials 16 00458 g002
Figure 3. Model of BC-790 five-layered corrugated board: (a) without imperfections; (b) with imperfections (scale ×10).
Figure 3. Model of BC-790 five-layered corrugated board: (a) without imperfections; (b) with imperfections (scale ×10).
Materials 16 00458 g003
Figure 4. Cross-section of three-layer corrugated cardboard.
Figure 4. Cross-section of three-layer corrugated cardboard.
Materials 16 00458 g004
Figure 5. Cross-section of five-layer corrugated cardboard.
Figure 5. Cross-section of five-layer corrugated cardboard.
Materials 16 00458 g005
Figure 6. Absolute estimation error obtained using five different models.
Figure 6. Absolute estimation error obtained using five different models.
Materials 16 00458 g006
Table 1. Component papers and geometry of the corrugated boards.
Table 1. Component papers and geometry of the corrugated boards.
WaveGrammageComponentHeightPeriodTake-Up
Type(g/m2)Papers(mm)(mm)Factor
B410TL3125/WS120/TL31252.556.341.337
C590KLB170/S.C.175/KLB1703.637.951.427
E480TLWC160/WS 135/TLW1601.163.501.236
BC790KLB170/WS135/WS80/WS135/KLB170---
BE600TLW140/WS95/WS80/WS95/TL3125---
BE590TL3125/WS95/WS80/WS95/TL3125---
Table 2. Basic parameters of cardboards.
Table 2. Basic parameters of cardboards.
Wave TypeGrammage (g/m2)Cardboard IDHeight (mm)
B410B-4102.912
C590C-5904.110
E480E-4801.586
BC790BC-7906.740
BE600BE-6004.150
BE590BE-5904.120
Table 3. Designation of the constituent papers of the tested corrugated board.
Table 3. Designation of the constituent papers of the tested corrugated board.
Cardboard IDType of PaperGrammage (g/m2)Layer
B-410TL3125Liner ext.
WS120Fluting
TL3125Liner intern.
C-590KLB170Liner ext.
S.C.175Fluting
KLB170Liner intern.
E-480TLWC160Liner ext.
WS135Fluting
TLW160Liner intern.
BC-790KLB170Liner ext.
WS135Fluting (B)
WS80Liner
WS135Fluting (C)
KLB170Liner intern.
BE-600TLW140Liner ext.
WB95Fluting (E)
WB80Liner
WB95Fluting (B)
TL3125Liner intern.
BE-590TL3125Liner ext.
WS95Fluting (E)
WS80Liner
WS95Fluting (B)
TL3125Liner intern.
Table 4. Average SCT values in CD and average tensile stiffness values of all constituent papers in the MD, CD and in 45 degree directions.
Table 4. Average SCT values in CD and average tensile stiffness values of all constituent papers in the MD, CD and in 45 degree directions.
Cardboard
ID
Paper
ID
SCT-CD
(N/mm)
TS-CD
(N/mm)
TS-MD
(N/mm)
TS-45
(N/mm)
B-410TL3-1252.14373.331012.7572.68
WS-1202.09365.061024.6516.50
TL3-1252.09381.151058.3595.08
C-590KLB-1703.28527.771472.1929.12
SC-1754.18686.081476.1924.73
KLB-1703.19568.051445.1956.22
E-480TLWC-1602.75412.141043.6635.01
WS-1352.13365.021067.5533.52
TLW-1602.43443.861102.1667.03
BC-790KLB-1703.39618.881534.2990.02
WS-1352.19369.081113.5572.41
WS-801.50317.09699.14445.02
WS-1352.23385.911147.4623.54
KLB-1703.30592.741418.9838.15
BE-600TLW-1402.61505.95999.95622.60
WS-951.69331.94872.70498.63
WS-801.42273.16812.83424.82
WS-951.52290.66885.54508.44
TL3-1252.13440.631082.2623.26
BE-590TL3-1202.26412.94961.28586.98
WS-951.50294.34756.40427.40
WS-801.47343.54696.51459.84
WS-951.75332.01854.73474.40
TL3-1252.32413.77883.13588.29
Table 5. Empirically determined coefficients A (column 2–4), a (column 5) and a * (column 6).
Table 5. Empirically determined coefficients A (column 2–4), a (column 5) and a * (column 6).
Cardboard
ID
A
Equation (4)
a
Equation (6)
a *
Equation (7)
(1)(2)(3)(4)(5)(6)
B-4100.7810.7910.81952-
C-5900.778
E-4800.815
BC-7900.7280.8480.819-18
BE-6000.869
BE-5900.947
Table 6. Reference (measured) ECT values and estimation values with empirical–analytical and numerical models.
Table 6. Reference (measured) ECT values and estimation values with empirical–analytical and numerical models.
Cardboard
ID
ECT (ref)
(N/mm)
ECT (M1)
(N/mm)
ECT (M2)
(N/mm)
ECT (M3)
(N/mm)
B-4105.485.545.595.27
C-5909.689.359.819.19
E-4806.376.326.476.17
BC-79010.4110.6510.759.89
BE-6008.959.129.078.67
BE-5909.689.229.899.24
Table 7. Estimation error obtained using five different models.
Table 7. Estimation error obtained using five different models.
Cardboard
ID
Error (M1)
(%)
Error (M2)
(%)
Error (M3)
(%)
Error (M4)
(%)
Error (M5)
(%)
B-4101.092.01−3.834.981.39
C-590−3.411.34−5.065.211.61
E-480−0.781.57−3.140.45−2.99
BC-7902.313.27−5.0012.5116,49
BE-6001.901.34−3.13−5.76−2.43
BE-590−4.752.17−4.55−13.54−10.48
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Garbowski, T.; Knitter-Piątkowska, A.; Winiarski, P. Simplified Modelling of the Edge Crush Resistance of Multi-Layered Corrugated Board: Experimental and Computational Study. Materials 2023, 16, 458. https://doi.org/10.3390/ma16010458

AMA Style

Garbowski T, Knitter-Piątkowska A, Winiarski P. Simplified Modelling of the Edge Crush Resistance of Multi-Layered Corrugated Board: Experimental and Computational Study. Materials. 2023; 16(1):458. https://doi.org/10.3390/ma16010458

Chicago/Turabian Style

Garbowski, Tomasz, Anna Knitter-Piątkowska, and Piotr Winiarski. 2023. "Simplified Modelling of the Edge Crush Resistance of Multi-Layered Corrugated Board: Experimental and Computational Study" Materials 16, no. 1: 458. https://doi.org/10.3390/ma16010458

APA Style

Garbowski, T., Knitter-Piątkowska, A., & Winiarski, P. (2023). Simplified Modelling of the Edge Crush Resistance of Multi-Layered Corrugated Board: Experimental and Computational Study. Materials, 16(1), 458. https://doi.org/10.3390/ma16010458

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