Introduction
Recent trends and increasing demand in corrugated board packaging market require from manufacturers and researchers, the development of innovative solutions to provide ease of shaping and attractive appearance of the packaging, robustly supported by sufficient box strength. Such needs highlight, in the engineering community, specific computational and experimental research challenging subjects, both from a material standpoint[
1] and from a packaging point of view[
2].
Proper understanding of packaging structural behaviour requires, as a first step, a detailed knowledge of mechanical properties of the employed material, namely of cardboard. Corrugated board is built as a sandwich composite material with individual layers, alternatively structured by flat and corrugated papers, usually ranging from two to seven layers. A practical classification of the corrugated layer, called fluting, based on wave height is usually adopted; it is denoted by capital letters, typical wave heights are from A to F. Due to constitutive material properties and internal composite structure, the layered corrugated board typically displays two characteristic in-plane directions of orthotropy, namely, the machine direction (MD), perpendicular to the main axis of the fluting, and cross direction (CD), parallel to the fluting, directly affecting the mechanical response of the paperboard, both in elasticity range and for strength thresholds. In recent literature, numerous works are available to accurately model the mechanical constitutive behaviour of corrugated board, upon computational developments, accounting for anisotropic behaviour,[
3,
4] plastic behaviour,[
5] creep response,[
6] forming process,[
7,
8] creasing and folding conditions,[
9,
10,
11,
12,
13] as also effectively supported by experimental testing (see, e.g.,[
9,
14,
15]). Consistently, a crucial role in the design of corrugated board packaging is played also by the evaluation of paperboard strength and buckling resistance, e.g., in [
16,
17,
18].
Toward the goals of structural modelling in packaging research and applications, the above mentioned refined understanding and modelling capability for corrugated board requires to be associated with a detailed knowledge of various packaging solutions, both in boxes and trays, particularly with reference to functionality,[
19,
20,
21] numerical modelling,[
22,
23,
24,
25] experimental testing[
26,
27] and, specifically, to box strength.[
28,
29,
30]
In view of reliable structural applications, over the years, several approaches have been developed, specifically conceived as aiming at estimation of box strength of corrugated board packaging. In particular, various approaches devised analytical formulae, although restricted to typical design, subsequently improved by the adoption of numerical tools, such as finite element modelling, toward broader applications or solution of demanding specific issues, such as consideration of layered structure, possibly tackled by homogenisation methodologies. In order to afford such computational and experimental challenging tasks in suitable unified and standardised procedures, two tests have been widespread, namely the Box Compression Test (BCT) and the Edge Crush Test (ECT), as nowadays widely employed in industry. However, despite their measuring effectiveness, these two testing methods are not sufficient, in current practice, in providing reliable data for strength estimation by computational approaches.
Within such context, several authors, in recent years proposed experimental and computational advancements in box strength estimation of corrugated board packaging. In particular, to the aims of the present paper, it is worth to mention the consideration of material moisture level,[
31] the analysis of buckling effects,[
32] the estimation of compressive strength accounting for openings,[
33] perforations,[
34] shifted flaps[
35] and open-top configurations,[
36] the modification of compression test suitable for trays[
37] and the extension to drop tests and gable-top shapes.[
38]
Despite research innovative contributions and proposed advanced computational methods, in the field of box strength estimation, the demanding task of parameter identification continuously highlights challenging situations both from problem complexity and from time-computational cost viewpoints. Such difficulties may be overcome, in and efficient and reliable way, by the adoption of Artificial Intelligence (AI) strategies, suitable to mimic neural schemes and to reproduce the behaviour of complex structured systems. In the last decades, AI and, specifically, Artificial Neural Network (ANN) methodology have been successfully applied, with growing interest in various engineering and multidisciplinary research fields, such as structural engineering (see, e.g., [
39,
40,
41,
42]), biomedical engineering (see, e.g., [
43,
44,
45]), agricultural engineering (see, e.g.,[
46,
47,
48,
49]). In particular, in corrugated paperboard research and related applications, ANNs have been employed limited to calibration of mechanical constitutive parameters,[
50] estimation of edge crush resistance,[
51] evaluation of effects by hand and ventilation holes on box compressive strength.[
52]
The main objectives of the present work are devoted to the estimation of corrugated board box strength, in various material and structural configurations, relying on ANN models. Therefore, tackling the demanding task of box strength evaluation by an AI approach aims at providing beneficial contributions regarding wide applicability of the devised method joined with reliability of the computational approach, toward an innovative methodology, also suitable for practical engineering applications.
Following the current Introduction, the paper is organised in two main sections. Section 2 presents the adopted methods and the selected materials for the numerical analyses, particularly with reference to the overall research investigation methodology (Section 2.1), to the collected data to be employed as input dataset in ANN processing (Section 2.2), to the assessed load capacity of the packaging to be analysed (Section 2.3) and to the used ANN structure (Section 2.4). Consequently, Section 3 gathers and discusses the obtained results toward effective reduction of the number of ANN calibration input parameters. A final section, as Conclusions, highlights the innovative contributions and summarises the main steps of the study.
Results and Discussion
Due to the relatively small set of input data, the selection of input parameters of the ANN plays an important role. Therefore, in the paper, the attention was paid to the analysis of which set of input parameters would ensure the smallest ANN estimation error, in order to not exaggerate the number of input parameters. Several cases have been selected according to the theory of mechanics and subject knowledge.[
32,
59] Finally, the following four ANN cases with different sets of input parameters were selected:
ANN1: Packaging dimensions (, and ), representative properties of cardboard as an orthotropic composite (, , , and ), , effective thickness and critical forces of load-bearing panels (, , and )
ANN2: Packaging dimensions (, and ), and critical forces of load-bearing panels (, , and )
ANN3: Ratios of packaging dimensions ( and ), , and critical forces of load-bearing panels (, , and )
ANN4: Packaging dimensions (, and ) and critical forces of load-bearing panels in relation to (, , and )
For clarity, the summary of the input parameters for all types of ANNs considered is presented in
Table 1.
For data explained in Section 2.2 and grouped above, the training of four types of ANNs was executed 1000 times due to different starting point in minimization algorithm, various initial values of weights and different subdivision of training/testing sets. Each time for randomly taken 45 packaging out of total 54 box samples collected from our previous research papers. The rest, i.e., nine packaging out of 54 box samples were left for testing ANNs accuracy. The accuracy was measured by the root mean square error:
in which
is the predicted value, while
is the reference (tested) value of BCT and
is the box sample number,
for training set and
for test set.
Due to the training based on stochastic selectin of 45 samples and completely random distribution of the initial values of the network parameters, the multiple ANNs were obtained for which the RMSEs were computed. For each of four types of ANN, the best one (with lowest RMSE) was selected, see pink bars in
Figure 5. Then, for those four best ANNs, the RMSEs for test set was computed, see grey bars in
Figure 5.
Detailed results for predicting BCT for test set (9 box samples) were presented as box compression strengths (so called BCT indices) in
Table 2. In second row, the real values from press machine were shown, while in later rows the ANN estimations were demonstrated. In
Figure 6, the same results of best ANNs were divided to the reference strengths values (BCT), therefore, the bars show how far from the real values the best ANN estimations are. Bar equal 0.0 would mean the perfect fit to the reference (BCT) value.
Comparing various possible input parameter sets studied in this research, one can notice that the minimal value of the RMSE was obtained for ANN2, in which eight input parameters were taken into account, see
Figure 5. This is also confirmed while comparing the results for single cases in the test set, see
Figure 6, it is clearly visible that the value of predicted BCT for the ANN2 was closest to the reference values in most cases. This suggest that the number of inputs of the ANN for predicting the BCT can be limited to eight parameters: packaging dimensions (
,
and
),
and critical forces of load-bearing panels (
,
,
and
). The RMSE obtained for this case was equal to 8.2% in the training set and similar value, i.e., 8.6%, in the test set. Moreover, the errors for single test cases from 1 to 9 were 7.9%, 9.6%, 7.7%, 9.3%, 3.0%, 5.1%, 7.2%, 15.5% and 6.0%, respectively, which is a good result when estimating the compressive strength for such different types of packaging by a single ANN model.
The worst results (in the training set) were obtained for the ANN3, in which ratios of the packaging dimensions ( and ) were used as the inputs instead of all packaging dimensions separately (, and ). Furthermore, for the ANN4 the RMSE obtained in the training set was lower than for ANN1 and ANN3. However, one can notice that in this case the worst results were obtained in the test set. The RMSE in the test set is about two times greater than in the training set. It shows that the ANN4 was overfitted to the training data, the results are satisfactory in the training set, but are not sufficiently general for the other data.
The results presented in this study can be easily compared with the other results presented in the previous papers on estimation of the BCT in various cases (boxes with openings, perforations or without any changes in their construction).[
32,
33,
34] In these previous studies, in which an analytical-numerical approach was proposed for estimation of the BCT, the authors modified also parameters in the well-known McKee formula in order to obtain the optimal values for the specific case (boxes with openings or perforations or without both). For the boxes with openings,[
33] the authors obtained the mean error of 6.5% for the optimal parameters. In the case of boxes with perforations,34 the mean error was equal to 3.5%. For the boxes without holes or perforations,[
32] the authors obtained the mean error in the range of 5% to 8%. One can notice that the values of optimal parameters obtained in both cases were totally different, e.g., for boxes with openings the optimal parameters achieved were
k = 0.755 and
r = 435, while for boxes with perforations the optimal parameters obtained were
k = 0.4 and
r = 0.75. In the current study, the error measure obtained was equal to 8.2% in the training set and 8.6% in the test set, but in this study the ANN approximation is more universal (i.e. is valid for all kind of boxes e.g. with perforations, opening and without them). The formulas with their optimal parameters proposed in the papers[
32,
33,
34] were dedicated to specific cases of box designs, while here the ANN model is much more universal and can be used for any box geometry with possible holes and perforations.
The methodology presented in this study with the application of ANN models for BCT prediction in various cases can be repeated and the results can be improved in the future for bigger number of data, which should greatly generalize the propose methodology.. If the average error of prediction in the new test set will become bigger than in the current test set, then the training process can be repeated with representation of new data both in training and test sets. This approach leads to an obvious asymptote, which is a kind of limit to the possibility of adapting the model for the selected neural network architecture to new data, which can also be improved by rebuilding the network architecture when the amount of training data allows it.