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

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.


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 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.

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. 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 [74].
(a) (b) 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 [75]. 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.

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.  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.

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.

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.

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: 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: 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 ν xy is the Poisson's ratio in the xy plane, which can be determined by the empirical formula proposed by Baum [21]: 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 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).
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 [61,62].

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: 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].

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: 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.

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: 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: 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).

=1
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 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.

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 that appears in Equation (4). Such a model can be written, for example, in the following form: where the 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 can be written as: where is a parameter that should be selected once and for all on the basis of experimental observations (empirical adjustment), is the height of the cross-section, is the basis weight of a given layer, 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 from the Equation (6), which then takes the form: 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: 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: which, in this case, leads to: b 3 = 0.7P 1 .
which, in this case, leads to: 3 = 0.7 1 . It is worth noting that the parameter 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 [76,77] or numerical [78,79] 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.  Table 4     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. Table 4 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).

The ECT Estimations by Numerical and Analytical-Empirical and Empirical Models
Finally, five models listed below were used for ECT estimation: • Model 1: proposed analytical-empirical model, Equation (5) 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.  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.

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 singlewalled 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.

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.