Optimal Versus Equal Dimensions of Round Bales of Agricultural Materials Wrapped with Plastic Film—Conﬂict or Compliance?

: For the assumed bale volume, its dimensions (diameter, height), minimizing the consumption of the plastic ﬁlm used for bale wrapping with the combined 3D method, depend on ﬁlm and wrapping parameters. Incorrect selection of these parameters may result in an optimal bale diameter, which differs signiﬁcantly from its height, while in agricultural practice bales with diameters equal or almost equal to the height dominate. The aim of the study is to formulate and solve the problem of selecting such dimensions of the bale with a given volume that the ﬁlm consumption is minimal and, simultaneously, the bale diameter is equal or almost equal to its height. Necessary and sufﬁcient conditions for such equilibria of the optimal bale dimensions are derived in the form of algebraic equations and inequalities. Four problems of the optimal bale dimension design guaranteeing assumed equilibrium of diameter and height are formulated and solved; both free and ﬁxed bale volume are considered. Solutions of these problems are reduced to solving the sets of simple algebraic equations and inequalities with respect to two variables: integer number of ﬁlm layers and continuous overlap ratio in bottom layers. Algorithms were formulated and examples regarding large bales demonstrate that they can handle the optimal dimensions’ equilibria problems.


Introduction
Currently, the demand to limit consumption of the film used to wrap bales of agriculture materials has been receiving increasing attention to reduce both costs and damage to the environment caused by plastic waste [1][2][3][4]. For this purpose, both experimental studies comparing film usage for different wrapping conditions and methods [5][6][7][8][9] and model-based analytical approaches useful for estimation of the film usage [10][11][12][13][14] and its minimization [15][16][17][18] were used. In the last few years it has been shown that by appropriate optimal selection of the film width [15,16], overlaps between adjacent film strips [17], and film width and overlaps, together, [18] it is possible to reduce film consumption by up to 20%. The model-based optimization research concerned mainly conventional [15,17,18] and IntelliWrap [16] wrapping methods. The monograph [19] comprehends the issues of modelling of the film consumption for wrapping round bales and related optimization problems solved in the last few years.
Also, an appropriate choice of bale size dimensions (diameter, height) may guarantee decreasing film consumption [11,19]. In [11], where the dependence of the film consumption on the bale diameter for the conventional wrapping technique was investigated, the analytical analysis showed that the larger the bale diameter is, the lower is the film consumption per unit of bale volume, which led to the conclusion that the use of the bale with the largest permissible diameter ensures the smallest film consumption. These studies were based on a rough model which describes film usage as a continuous function of bale dimensions, film width and a number of wrapped film strips; however, the model did not take into account mechanical properties of the stretch film and the direct relation between the number of wrapped film strips and bale and film parameters.
The concept of optimal (from the point of view of film consumption) selection of bale dimensions for the assumed bale volume was first introduced in [15] for the conventional wrapping method. A more accurate model was used that describes the consumption of the film as a function of the bale and film dimensions, mechanical parameters of the film (Poisson ratio, unit deformation), overlap ratio and the number of bale rotations [12]. The film consumption per unit of the bale volume is used as a measure of film usage.
Since it is very difficult to find the optimal bale dimensions minimizing the original exact film usage index due to the discontinuity of this index, near-optimal parameters, being as important as optimal parameters for engineering applications, are sought. It was shown in [15] that in the case of using the conventional wrapping method, the optimal bale height is twice its optimal diameter. This optimality rule, which holds also for the IntelliWrap method [19], has only a theoretical character and is a consequence, among others, of multiple overlapping segments on the bale cylinder top and bottom, where there are 2-16 times more film layers than on the bale's lateral surface [10,19]. However, in the case of the combined 3D wrapping method [6,14], the optimal bale dimensions turn out to be useful and applicable from the engineering practice point of view.
In the 3D wrapping method, which offers the potential to minimize film usage [6,14,20] as well as to enhance the quality of silage [6], biaxial rotation of the film applicators results in two types of film layers wrapped perpendicularly: the bottom layers are wrapped around the bale's lateral surface and the upper layers are wrapped along the bale's longitudinal axis. For a detailed description of the wrapping process, see [6] and the producer's documentations [21][22][23]. The problem of the choice of the best bale dimensions to guarantee the minimal film consumption was solved in [20], where the optimality conditions were established in the form of algebraic cubic equations, which can easily be solved using both analytical and numerical methods. Analytical and numerical studies [19,20] have shown that the relation between the optimal diameter and height is not evident. The optimal diameter can be larger than the bale height [19] (Figure 6.5b), [20] (Figure 11c) or can be smaller than the bale height [19] (Figure 6.5d), [20] (Figure 11a), depending on the bale volume as well as film and wrapping parameters. But the optimal bale diameter can also be equal to its height [19] (Figure 6.5a,c), [20] (Figure 11b). The last case corresponds to typical large bales of 1.2 m diameter and height [10,24,25]. However, in agricultural practice typical bales are 1.2 to 1.3 m diameter and height [6,10,19], where diameter is equal [26][27][28][29][30] or almost equal [9,31,32], but not much bigger, than bale height. Other dimensions of large round bales, for example 1.2 m × 1.6 m Ø [33], or 1.2 m × 1.5 m Ø [34] or 1.5 m Ø [35], are used less frequently. Although round bales of the height greater than the diameter are also being investigated, for example 1.2 m × 0.9 m Ø [36], they are much less common in agriculture practice [37]. Standard 3D bale wrappers are designed to wrap cylindrical bales up to 1.6 m in diameter and up to 1.2 m in height [22], or up to 1.5 m diameter and 1.2 m height [21,23]. Therefore, the question becomes how to choose such film and wrapping process parameters for a round bale with a given volume wrapped using the 3D method so that the bale diameter minimizing film consumption is equal to or almost equal to its height. The aim of the paper is to solve the problem of the selection of the film and wrapping parameters (overlaps, numbers of bottom and upper film layers, etc.) that for a given bale volume, i.e., a given bale weight, its diameter and height are equal or near-equal and, at the same time, they minimize the film consumption. It is, in essence, the inverse problem in which the film and wrapping parameters are determined that the equilibria of optimal dimensions are guaranteed, while in the direct problem considered in [20] the optimal dimensions for a given parameters are sought. The thesis is that for any given bale volume and number of global film layers there exist film and wrapping parameters for which the optimal bale dimensions are equal or near-equal, with the pre-assumed proportions diameter/height. To prove the above, the necessary and sufficient conditions of equilibrium and near-equilibrium of the optimal bale dimensions were derived in the form of algebraic equations. The equilibrium conditions are dichotomous, as some depend on the given bale volume and others apply to any bale volume. Then, the problems of selecting such bale wrapping parameters for which the optimal bale dimensions are equal or near-equal were formulated and solved, separately, for a given and an arbitrary bale volume. Suitable design algorithms were proposed and numerically verified for large bales.

Materials and Methods
In this section, the notions of equal and near-equal optimal dimensions of a round bale are introduced and the related necessary and sufficient conditions are derived. The research methodology is also described.

Equal and Near-Equal Optimal Bale Dimensions
A complete mathematical model describing stretch film consumption for wrapping a bale with p b bottom and p u upper film layers using the combined 3D technique was derived in [14]; the main formula describing film usage is recalled in Appendix A with the corresponding assumptions. Symmetry of the bale is assumed. Mechanical properties of the plastic film are described by the Poisson's ratio v f and unit deformation ε l f ; thickness of the film is ignored, e.g., polyethylene film is 25 µm thick [38]. The main symbols are summarized in Nomenclature, Appendix C.
The problem of the choice of bale dimensions (diameter, height) minimizing the consumption of the film used to wrap a cylindrical bale by 3D method has been stated and solved for the first time in [20], where the necessary and sufficient condition of the existence of the unique optimal bale diameter D * b of the bale of pre-assumed volume V b0 was derived in the form of cubic equation with zero linear term coefficient [20]: where b f r is the width of stretched film described by Equation (A1), δ denotes the overlap of the extreme film strips in bottom film layers at the bases of the bale, k f b and k f are the overlap ratios determining the width of the contact between adjacent film strips in bottom and upper film layers; function Ω k f is defined by Equation (A3). The coefficients in the first and second terms and the free term of Equation (1) depend on the pre-assumed bale volume as well as all the film and wrapping parameters. Based on the optimality condition, Equation (1), analytical and numerical analysis of the influence of film width, pre-assumed bale volume and numbers of bottom and upper film layers on optimal bale dimensions and optimal film consumption were carried out in [20], where many detailed conclusions regarding the impact of these parameters were formulated. Additionally, the influence of the bottom layers overlaps was studied in [19]. Equation (1) has one real positive root [20]. The corresponding optimal bale height For given V b0 the optimal D * b and H * b depend on film and wrapping parameters, they are linearly dependent and the difference between them is described by From a quick inspection of Equation (3) it can be seen that the relation between the optimal diameter D * b and height H * b is not evident. The conducted research [19,20] has shown that the relations D * b > H * b and D * b < H * b , as well as D * b = H * b are possible, depending on the values of the film and wrapping parameters and volume V b0 ; compare [20] ( Figure 11), where the proportions diameter/height are depicted for large bales. Only the relation D * b ≥ H * b , desirable from the point of view of baling systems engineering, will be considered. The proportion diameter/height do not exceed 1.25 [22] (RW 1819, 150/120), or 1.2 [21] (BW 1850, 150/125), or [23] (WM 1851, 150/125); however, predominantly, it is of 1.05 order [21][22][23].
The goal of this paper is to study the selection of equal and nearly-equal optimal bale dimensions, therefore, the respective precise definitions are given.
where a sufficiently small ε 0 > 0. Then D * b and H * b are called nearly-equal with order ε 0 .

Methodology
In this paper, a model-based analytical approach was applied, which addressed the goals (derivation of the conditions for equal and near-equal optimal bale dimensions and solving the problems of wrapping parameter design guaranteeing such equilibria) by using mathematical tools. The solutions of the optimality conditions for the bale diameter were obtained by applying numerical tools.
Firstly, the necessary and sufficient conditions under which for the given bale volume the optimal bale diameter is near-equal to its optimal height with order ε 0 were derived using the optimality conditions expressed by Equations (1) and (2); Proposition 1 abstracts these results (Section 2.3.1). The conditions for optimal bale dimensions equilibria are given by algebraic equations. One of them uniquely relates bale volume to order ε 0 and film and wrapping parameters, i.e., this condition is satisfied for a given volume. The complementary condition determined by two simple equations relating only to ε 0 and film and wrapping parameter applies regardless of the bale volume. Based on these conditions, using differential calculus, the relationships between the volume of the bale with near-equal optimal dimensions and the order ε 0 , width of the film and overlaps were analyzed (Section 2.3.2). Next, laying ε 0 = 0 the conditions for exact equilibrium of the optimal bale dimensions was obtained directly from the conditions of their nearequilibrium (Section 2.3.3).
Knowing the conditions of equilibrium and near-equilibrium of optimal bale dimensions, the tasks of designing the wrapping process in such a way that these conditions would be satisfied were formulated and solved. Due to the dichotomous nature of the equilibrium conditions, design tasks for a given bale volume and those in which the bale volume is arbitrary were considered separately; these are covered in Section 3.1 and Section 3.2, respectively. Consequently, four design problems were solved; however, the solutions for equal optimal bale dimensions resulted directly from those for near-equilibrium by substituting ε 0 = 0 (Sections 3.3 and 3.4). It was assumed that the film parameters (width, mechanical parameters v f , ε l f ) are known, i.e., a practically available plastic film can be used, e.g., a commercial PE film used traditionally due to its mechanical properties and low costs [26,27,39]. Also the global number of film layers p l wrapped on the bale's lateral surface was taken as a given. Many studies have investigated the number of desired film layers for baled silage preservation, for example [7,29,40,41]. Mostly, four, six or eight layers of film are applied [6,24,35]; however, ten, twelve, and even sixteen film layers in which the silages are wrapped are also considered [7,22,23,29]. Therefore no specific assumptions were taken concerning the number of global film layers; the numerical studies were conducted for p l from four to sixteen. The standard overlaps 50%, 67% or 75% resulting in uniform film coverage [17] were assumed for upper film layers. Thus, the bottom layers overlaps δ, k f b and film layers decomposition p l = p b + p u had to be selected. The upper and lower limits were taken for the overlaps, the selection of which should take into account the knowledge and experience of baling as well as the possibilities of the wrappers available. Consequently, the considered design problems consisted in solving a set of algebraic equations (resulting from the conditions of the equilibria of optimal bale dimensions) and inequalities (the constraints of variables) with respect to three decision variables. Two variables-the overlaps δ, k f b -were continuous, while the third-the number of bottom film layers p b -was integer. These mixed (hybrid) sets of equations and inequalities were solved analytically, separately for fixed and free bale volume. The solutions are abstracted by Propositions 4-7. For volume-free problems the solution is given by two inequalities directly related to the integer p b and simple rules for computing the overlaps δ, k f b . In the case of fixed bale volume the solution was derived in the form of two inequalities related to the overlap ratio k f b , which must be verified for every considered p b and unique algebraic rule for the determining of the overlap δ. If equal optimal bale dimensions are sought then these inequalities and formulas take particularly simple forms which are obtained by substituting ε 0 = 0. Computational algorithms were developed which enable designation of the sought wrapping parameters in a few steps. The examples illustrated how to use them (Excel is enough) as well as the effectiveness for optimal bale dimensions selection for a standard large bale.
In sum, in the multistage process several problems of optimal equal and near-equal bale dimensions design for cylindrical bales of a given volume were formulated and solved. The research framework is graphically shown in Figure 1, which also illustrates the relations between these problems.

Near-Equal Optimal Dimensions
By Definition 2, the optimal bale dimensions are such that Equation (4) holds and the optimality condition expressed by Equation (1) is satisfied. For a given bale volume V b0 Equations (2) and (4) yield and The necessary and sufficient conditions of near-equal optimal dimensions are given by the following proposition, which is proved in Appendix B. Proposition 1. For given: bale volume V b0 , width b f of the film and its mechanical parameters v f , ε l f , numbers of film layers p b , p u and the overlaps δ, k f b , k f , such that the applicability condition expressed by Equation (A4) holds, the optimal bale diameter D * b is near-equal to its optimal height H * b with order ε 0 in the sense of Definition 2 if and only if one (and only one) of the conditions is satisfied: (i) order ε 0 and the bale, film, and wrapping parameters are related by the equation (ii) the equations and are satisfied, simultaneously. The optimal bale diameter D * b and height H * b are described by Equations (5) and (6), respectively; in case (i) for a given bale volume, in case (ii) for an arbitrary volume V b0 .
The necessary and sufficient conditions for near-equal optimal bale dimensions, case (ii), are independent of V b0 ; they are volume-free. However, those from case (i) hold for given V b0 , i.e., they are volume-fixed.
Note, that the inequality must hold to satisfy Equation (7); however, the expression 2δ − b f r k f b may be positive or negative, depending on the sign of the denominator. The necessary condition expressed by Equation (10) is independent on V b0 , while Equation (7) depends on the particular value of bale volume V b0 . Both conditions depend on ε 0 .

Bale Volume of Near-Equilibrium
From Equation (7), provided that inequality from Equation (10) is satisfied, we have whence Thus, for given film and wrapping parameters, bale volume resulting in near-equal optimal dimensions is monotonically increasing function of ε 0 whenever and decreases with ε 0 , if 2δ < b f r k f b . If inequality from Equation (12) holds, then the volume V b0 decreases if the overlap k f b grows, while in the opposite case the influence of k f b is not so evident and depends on the sign of the expression in the square brackets of the numerator in the last fraction of the right hand side of the following equation If condition from Equation (12) holds, then V b0 given by Equation (11) grows with the increase of the overlap δ. In the opposite case the larger δ is, the smaller is the volume resulting in near-equal optimal dimensions. From the analysis of Equation (11) it also follows that the greater film width is, the smaller is V b0 , provided that inequality from Equation (12) holds. The influence of the bottom layers overlaps k f and δ is illustrated in Figure 2a,b, where near-equilibrium volume V b0 is depicted as a function of ε 0 for plastic film (e.g., polyethylene, PE) characterized by Poisson's ratio v f = 0.34 [38] and unit deformation ε l f = 0.7(−) [7,12,14] of popular width b f = 0.75 m [21][22][23]29,42]. Four bottom and upper film layers p b = p u = 4 and overlap ratio k f = 0.5 are assumed. Figure 3a,b illustrate the effect of the film width b f for commonly used overlap δ = 0.2 m [21][22][23] and two overlap ratios k f b = 0.3, 0.35; other parameters remain unchanged. The volume V b0 = 1.357 m 3 corresponding to a standard large bale of diameter and height [24,25] is also marked on these pictures. Finally, note that for standard overlap δ = 0.2 m [6,14,[20][21][22][23], PE film and since k f b < 1 2 , the condition from Equation (12)    Since in practically used wrapping systems the condition expressed by Equation (12) is usually satisfied and the near-equilibrium bale volume increases with the order ε 0 , a question arises about the conditions to be met by the film and wrapping parameters to guarantee for any 0 ≤ ε ≤ ε 0 the existence of a bale volume such that the optimal dimensions D * b , H * b are near-equal with order ε. The following result, derived in Appendix B based on Proposition 1, answers this question. The notation V b0 (ε), D * b (ε), and H * b (ε) is introduced, locally, for near-equal bale parameters to emphasize the relationship between ε and V b0 , D * b , and H * b .
Proposition 2. For given film parameters b f v f , ε l f , numbers of film layers p b , p u and the overlaps δ, k f b , k f , such that the applicability condition expressed by Equation (A4) holds, for any 0 ≤ ε ≤ ε 0 , where a sufficiently small ε 0 > 0, there exists bale volume V b0 (ε) > 0 such that the optimal bale diameter D * b is near-equal to its optimal height H * b with order ε not greater than ε 0 , i.e., if and only if and one of the following conditions is satisfied: (a) inequality from Equation (12) holds, hold, simultaneously. Then, for any 0 ≤ ε ≤ ε 0 and volume V b0 (ε) Equation (7) holds for (13) and (14) hold, but inequality expressed by Equation (15) not, then there exists positive ε such that 0 < ε < ε 0 and For any ε < ε ≤ ε 0 bale volume V b0 (ε) satisfying Equation (7) and optimal bale dimensions Inequality expressed by Equation (12) holds for narrow films, while inequality form Equation (14) is satisfied for wider films.

Equal Optimal Dimensions
In this case the optimal bale diameter is equal to bale height, i.e., by Equation (2) and the optimality condition expressed by Equation (1) is satisfied, simultaneously. The necessary and sufficient conditions for the equilibrium are given by the following proposition, which results directly from Proposition 1 for ε 0 = 0.
Proposition 3. For given: bale volume V b0 , width of the film b f and its mechanical parameters v f , ε l f , numbers of film layers p b , p u and the overlaps δ, k f b , k f , such that the applicability condition expressed by Equation (A4) holds, the optimal bale diameter D * b is equal to its optimal height H * b if and only if one (and only one) of the conditions is satisfied: (i) the bale, film and wrapping parameters are related by the equation (ii) the equation and Equation (9) are satisfied, simultaneously. The optimal bale dimensions are given directly by Equation (19) for fixed volume V b0 in case (i) or for arbitrary volume in case (ii).
Note, that in case (i) the necessary condition for bale dimensions equilibrium expressed by Equation (10) takes the form

Results and Discussion
In this section four problems of the selection of optimal bale dimensions being equal and near-equal are formulated, solved and illustrated by related examples.
It is known [17] that any overlap ratio of the form of irreducible fraction in which dividend is the divisor minus one, i.e., where q ∈ N , N denotes the set of positive integer numbers, results in the uniform film distribution on the bale's whole lateral surface and guarantee the same minimal film usage. An important special case for even p u is q = 2, which means 50% overlap between the successive film strips [6,17,24,43]. For k f = k f ,u function Ω k f = 1 and since p u = p l − p b , the applicability condition expressed by Equation (A4) takes the form From the practical perspective, only four or five smallest such overlap ratios are useful [6,17,22,24,43,44], i.e., 2 ≤ q ≤ 5 are worth considering. Bearing in mind the dichotomous nature of the equilibrium conditions specified by Propositions 1 and 3, the problems of volume-free and volume-fixed near-equal and equal bale optimal dimensions design will be considered, separately.

Volume-Free Near-Equal Optimal Bale Dimensions Design
Let us consider the following problem of the optimal choice of bale dimensions.

Problem 1.
Given film parameters b f , v f , ε l f , upper layers overlap ratio k f , Equation (20), number of global film layers p l , and the order of near-equilibrium ε 0 . Find film layers decomposition (p b , p u ) and the bottom layers overlaps where δ min and δ max are the smallest and largest admissible overlap δ, and where k f b,min and k f b,max are the smallest and largest admissible k f b , such that the applicability condition, Equation (21), holds and for any bale volume V b0 the optimal bale dimensions D * b , H * b are near-equal with given order ε 0 .
In view of Proposition 1, case (ii), the solution to the above problem exists if and only if there exist integer p b , p u and continuous parameters k f b , δ satisfying Equations (8) and (9), inequalities expressed by Equations (22) and (23), and the applicability condition, Equation (21). This set of equations and inequalities is solved in the Appendix B, this solution is summarized in the next result.
holds, then the solution to the Problem 1 exists if and only if there exists an integer p b such that inequalities: are satisfied together with the applicability condition, Equation (21). For any p b solving Problem 1 the overlap radio k f b is given by equation while by Equation (9) the overlap δ = 1 2 b f r k f b . For the assumed bale volume optimal bale dimensions D * b and H * b are given by Equations (5) and (6), respectively.

1.
Determine the set P 1 of all integer p b defined by the inequalities max p l for which there exists 2 ≤ q ≤ 5 satisfying applicability condition, Equation (21). 2.
If the set P 1 is empty, then the solution to Problem 1 does not exist-go to step 3. Otherwise, go to step 4.

3.
Change the lower k f b,min , δ min or upper k f b,max , δ max bounds of wrapping parameters, or the order ε 0 , or the film width b f and repeat the computations starting from step 1.

4.
For any p b ∈ P 1 compute the overlap ratio k f b according to Equation (27) and, next, If a bale volume V b0 is assumed, then for any p b ∈ P 1 compute film usage FC obtained for optimal D * b , H * b , Equations (5) and (6), using Equation (A2) and choose that p b which yields the minimal film consumption.
Application of Proposition 4 and Algorithm 1 is illustrated by an example.

Example 1
Plastic film of the parameters v f = 0.34, ε l f = 0.7(−) and b f = 0.75 m is assumed. Even numbers of global film layers 4 ≤ p l ≤ 16 and "uniform" overlap ratios k f , Equation (20), are considered. Firstly, the following bounds for bottom layers overlaps are assumed: Three orders ε 0 = 0.05, 0.1, 0.15 of near-equilibrium are considered. Inequality from Equation (24) is satisfied. Unfortunately, the set P 1 of integer solutions p b of the inequalities expressed by Equation (28) is empty for all ε 0 considered. When the lower constraint of δ is changed into δ min = 0.1 m, then for each ε 0 the set P 1 is composed of a few p b . They are listed in Table 1 together with respective overlaps k f b and δ. Bale volume V b0 = 1.357 m 3 is assumed, film usage for the optimal bale dimensions is computed and given in the last column of Table 1; for all p l for which the set P 1 is nonempty, there is only one p b ∈ P 1 . Note, that for given p l film usage obtained for D * b , H * b and decomposition (p b , p u ) are the same for all ε 0 considered. Thus, the smallest ε 0 can be chosen. For example, for V b0 = 1.357 m 3 taking ε 0 = 0.05 we obtain near-equal optimal dimensions D * b = 1.22 m and H * b = 1.162 m. The same regularity holds for wider film b f = 0.9 m, for which the solutions to Problem 1 are given in Table 2. Table 1. The numbers of global p l and bottom p b film layers and bottom layers overlaps k f b , δ for which Problem 1 of volume-free near-equal optimal bale dimensions design has solution for the assumed order of near-equilibrium ε 0 and bale from Example 1; film width b f = 0.75 m, the bounds for constraints from Equations (22) Table 2. The numbers of global p l and bottom p b film layers and bottom layers overlaps k f b , δ for which Problem 1 of volume-free near-equal optimal bale dimensions design has solution for given order of near-equilibrium ε 0 and bale from Example 1; film width b f = 0.9 m, the bounds for constraints from Equations (22) and (23) Parameter ε 0 influences the constraints expressed by Equations (25) and (26). If ε 0 grows, then for given p l both the lower and upper bounds for p b expressed by Equations (25) and (26) increase provided that the inequality form Equation (24) holds; however, these changes are not significant (see Table 1).

Effect of the Film Width
Film width b f , or equivalently b f r , influence the inequality from Equation (24); if b f r grows, then this inequality is still satisfied. Also, the lower and upper bounds from inequalities expressed by Equation (26) Table 2.

Volume-Fixed Near-Equal Optimal Bale Dimensions Design
Based on Proposition 1, case (i), the following problem of bale dimensions optimal design was formulated.   (30) are satisfied together with the applicability condition, Equation (21), where . (33) For any such p b and k f b the overlap δ is given by equation the optimal bale dimensions are described by Equations (5) and (6) and film usage where x denotes ceiling function [45]. Integer p * b solving Problem 2 is such that It may be easily shown that if the following inequality holds, then and condition from Equation (30) can be neglected. In all the tested examples the above inequality has been satisfied, compare Table 3 below. The inequality from Equation (37), similar to the condition expressed by Equation (12), holds especially if the number of bottom film layers is such that The inequalities specified in the above proposition must be solved for two variables: integer p b and continuous k f b . For any fixed p b the overlap δ given by Equation (34) increases linearly with growing k f b .

Algorithm 2
Assume film parameters b f , v f , ε l f , overlap ratio k f , Equation (20), number of global film layers p l , the order of near-equilibrium ε 0 , and bale volume V b0 are given. Take minimal δ min , k f b,min and maximal δ max , k f b,max values of the overlaps δ, k f b to guarantee the appropriate tightness of the wrappings.

1.
Determine the set P 2 of all integer p b for which there exists 2 ≤ q ≤ 5 satisfying applicability condition, Equation (21), and the set of overlap ratios K 2 (p b ) defined by the inequalities expressed by Equations (29) and (30) is nonempty.

2.
If the set P 2 is empty, then the solution to Problem 2 does not exist-go to step 3. Otherwise, go to step 4.

3.
Change the lower k f b,min , δ min or upper k f b,max , δ max bounds of wrapping parameters, or the order ε 0 , or the film width b f and go to step 1.

4.
Solve in p * b the integer programming task expressed by Equation (36) for p b ∈ P 2 . 5.
For the best p * b choose practically reasonable k f b ∈ K 2 p * b , compute the overlap δ = δ p * b according to Equation (34) and optimal D * b , H * b using Equations (5) and (6). The minimal film consumption is equal to FC p * b computed in step 4.

Example 2
Film parameters and overlaps constraints from Example 1 are taken, again. Bale volume V b0 = 1.357 m 3 is assumed together with orders of near-equilibrium ε 0 = 0.05 , parameters k f b,ne and closed intervals K 2 of k f b are given in Table 3. In all tested examples k f b,ne <ǩ f b (p b ). The upper bound of the closed interval of k f b defined by Equation (29) is determined bŷ For any k f b ∈ K 2 the formulas from Equation (34) for computing δ are presented in the penultimate column. In the last column film usage FC(p b ), Equation (35), is given. For most p b ∈ P 2 the closed intervals K 2 (p b ) are wide enough to select a practically convenient overlap ratio k f b . In some cases this interval is very narrow, e.g., for p l = 10, p b = 6, we have K 2 = [0.200, 0.204]; however, k f b = 0.2 is acceptable from the engineering point of view. For ε 0 = 0.05 and p l = 4, 6, 8 set P 2 is composed of only one number p b of bottom film layers and this p b solves Problem 2. For p l = 10, 12, 14, 16 there are at least two p b for which the set K 2 (p b ) is nonempty, the optimal film usage FC p * b defined in Equation (36) is marked by bold in Table 3. For ε 0 = 0.05, according to Equations (5) and (6) Table 3, the solutions to Problem 2 can be found for the greater order ε 0 = 0.1. Table 3. The numbers of global p l and bottom p b film layers, lowerǩ f b (p b ), Equation (32), and upperk f b , = k f b bounds defined in Equations (33) and (31), parameters k f b,ne , Equation (30), the non-empty sets of overlap ratios K 2 (p b ) defined by the inequalities expressed by Equations (29) and (30) and the linear functions, Equation (34), describing overlap δ for k f b ∈ K 2 (p b ) determined to find the solution to Problem 2 of near-equal optimal bale dimensions design for fixed volume V b0 = 1.357 m 3 , bale from Example 2; orders of near-equilibrium ε 0 = 0.05, 0.1, film width b f = 0.75 m, bounds for the constraints from Equations (22) and (23): (35), for p b ∈ P 2 and k f b ∈ K 2 (p b ); the optimal FC p * b defined in Equation (36) is marked by bold.

Volume-Free Optimal Equal Bale Dimensions Design
In this section the problem of volume-free design of equal optimal bale dimensions is considered being, in fact, a special case of Problem 1 for the order of near-equilibrium ε 0 = 0. Problem 3. Given: film parameters b f , v f , ε l f , overlap ratio k f , Equation (20), and number of global film layers p l . Find film layers decomposition (p b , p u ) and overlaps δ, k f b satisfying constraints expressed by Equations (22) and (23) such that the applicability condition, Equation (21), hold and for any bale volume V b0 the optimal dimensions D * b , H * b are equal.
The solution to the above problem is identical with that of Problem 1 for ε 0 = 0. Thus, the next result follows directly from Proposition 4 by lying ε 0 = 0 in inequalities expressed by Equations (25) and (26) and in Equations (5), (6) and (27). Inequality from Equation (24), here b f r > δ max , is identity.

Proposition 6. The solution to Problem 3 exists if and only if there exists an integer p b such that inequalities
are satisfied together with the applicability condition, Equation (21). For any such p b the bottom layers overlaps are uniquely determined by For a given bale volume optimal bale dimensions D * b and H * b are given by Equation (19).
Thus, to solve Problem 3 an integer p b fulfilling inequalities from Equation (39) must be found. The fractions dependent on b f r in "max" and "min" functions, which define the lower and upper bounds in this inequalities, increase with the film width. However, the upper and lower bounds may be determined by other arguments of "max" or "min", thus a larger b f r does not necessarily imply a greater p b . Proposition 6 yielded the following algorithm.

Algorithm 3
Assume film parameters b f , v f , ε l f , overlap ratio k f , Equation (20), and the number of global film layers p l are given. Take minimal δ min , k f b,min and maximal δ max , k f b,max values of the overlaps δ, k f b .

1.
Determine the set P 3 of all integer p b defined by the inequalities from Equation (39) for which there exists 2 ≤ q ≤ 5 satisfying applicability condition, Equation (21).

2.
If the set P 3 is empty, then the solution to Problem 3 does not exist-go to step 3. Otherwise, go to step 4.

3.
Change the lower k f b,min , δ min or upper k f b,max , δ max bounds of wrapping parameters or the film width b f and go to step 1.

4.
For any p b ∈ P 3 compute the overlap ratio k f b and, next, overlap δ according to Equation (40).

5.
For any p b ∈ P 3 and assumed bale volume V b0 compute film usage obtained for the optimal D * b , H * b , Equation (19), using Equation (A2) and choose that p b which yields the minimal film consumption.

Example 3
The same mechanical parameters of the film and wrapping parameters constraints as in the previous examples are assumed. Film widths b f = 0.75, 0.9 m and 4 ≤ p l ≤ 16 are considered. For b f = 0.75 m, similarly as in Example 1, the set of integer solutions p b of inequalities expressed by Equation (39) is empty for all p l considered, while for the wider film for p l = 14 solution to Problem 3 exists. When the lower constraint expressed by Equation (22) is changed into δ min = 0.1 m, then for some p l Problem 3 has a solution. These solutions are listed in Table 4. For every p l , nonempty P 3 is a singleton (a unit set), thus step 5 of the selection of best film layers composition p l = p b + p u can be omitted here; however film usage FC is added in the last column. Table 4. The numbers of global p l and bottom p b film layers and bottom layers overlaps k f b , δ for which Problem 3 of volume-free equal optimal bale dimensions design has a solution for bale from Example 3; film width b f = 0.75, 0.9 m, the bounds from Equations (22) and (23) (20), number of global film layers p l and bale volume V b0 . Find film layers decomposition (p b , p u ) and overlaps δ, k f b satisfying constraints from Equations (22) and (23) such that the optimal bale dimensions D * b , H * b resulting in the volume V b0 are equal and film usage is minimal.
the optimal bale dimensions are given by Equation (19) and film usage Integer p * b solving Problem 4 is defined by optimization task expressed by Equation (36).
Thus, the inequalities specified in the proposition must be solved for two variables: integer p b and continuous k f b . Next optimization task, Equation (36), must be solved and the overlap δ can be computed.

Example 4
Film parameters and overlap ratio k f as in the previous examples are given. Film widths b f = 0.75, 0.9 m, even p such that 4 ≤ p l ≤ 16 and bale volume V b0 = 1.357 m 3 are assumed. For any p l the sets of overlap ratios k f b satisfying inequalities from Equations (41) and (42) are given in the fourth column of Table 5. As previously, for any k f b from these closed intervals the linear formulas from Equation (45) for computing δ are given in the next column, in the last column film consumptions FC(p b ), Equation (46), are enclosed. For these p l for which there are more than one p b such that inequalities from Equations (41) and (42) are satisfied the minimal film usage, solving optimization task from Equation (36), is marked by bold. Note that in most cases film usage for b f = 0.9 m is smaller from that for b f = 0.75 m; this confirms the analysis of the film width influence on the optimal film usage from [14,20]-the broader film width is applied, the smaller minimal film usage is achieved.

Optimal Film Usage
For the combined 3D wrapping technique, the analysis of the optimal film consumption reached for D * b , H * b was carried out in [19,20], where it was shown, in particular, that for the bottom film layers p b p u the optimal selection of bale diameter and height may result in a 3.37% to even 23.13% reduction in film usage depending on the number of global film layers and wrapping parameters, which means up to 23% film cost savings. Many conclusions regarding the impact of film width, pre-assumed bale volume and numbers of bottom and upper film layers on near-optimal bale dimensions and near-optimal film consumption were formulated [19,20].

Conclusions
It has been shown that for bale wrapping by a 3D combined method for given film parameters and given a global number of film layers it is possible, both for assumed and for an arbitrary bale volume, to select such wrapping parameters that guarantee the preassumed equilibria of optimal bale dimensions. Thus, the "compliance" is the answer to the question from the title.
The derived necessary and sufficient conditions for the balance of optimal bale dimensions allowed for the formulating of four practical problems of selecting the wrapping parameters, in particular the overlaps in the bottom film layers and decomposition of the global number of layers into bottom and upper layers. The mathematical form of the conditions of optimal bale dimensions equilibrium indicated the legitimacy of considering two separate design strategies, one for fixed and one for free bale volume. Design algorithms were derived. These algorithms only require a solution with respect to the number of film layers and overlap ratio of the bottom film layers of two or four simple algebraic inequalities. Then, on the basis of very simple formulas, the remaining wrapping parameters and the optimal bale diameter and height can be determined. Examples regarding standard large bales demonstrate that the algorithms can handle the optimal dimensions' equilibria problems.
Optimal dimensions mean a reduction in film consumption by up to 23%. Simultaneously, equal and near-equal bale dimensions are essential from the agriculture engineering practice and bale wrapping technique point of view. The algorithms developed can be applied for design of wrapping processes using a 2D method for cylindrical bales of arbitrary agricultural materials, e.g., lignocellulosic agricultural residues [46,47], bearing in mind the minimisation of the film usage and the desired equilibrium between bale dimensions. Potential applications of the algorithms include the baling of waste, such as municipal solid waste [48,49], for example.

. Film Consumption
A mathematical model describing the consumption of the stretch film to wrap a cylindrical bale using combined 3D technique was derived in [14]. It was assumed that the subsequent film strips, which are wrapped in p b bottom layers on the bale's lateral surface, overlap one another, creating the overlap k f b b f r , where k f b is a dimensionless relative ratio determining the width of the contact between adjacent film strips and the film width after stretching b f r is described by [12]: where b f is the width of un-stretched film, v f and ε l f are the Poisson's ratio and unit deformation of the film. The extreme film strips are overlapped at the bases of the bale for δ, as shown in [14] ( Figure 3). It was assumed that 0 < k f b < 1 2 , for which only one film layer results for one wrapping cycle. For the wrappings of p u upper layers it was assumed that the subsequent film strips are wrapped along the bale's longitudinal axis with the overlap ratio k f [14] (Figure 2).
The global film consumption measured by the index FC defined as the ratio of the surface area of un-stretched film used to wrap the bale to bale volume, is described by the function [14]: where x denotes ceiling function [45], function [20] is introduced for brevity of the notation; x is the floor function [45]. The above formula indicates the dependence of the film usage on film parameters ε l f , v f , b f , bale diameter D b and height H b , the overlaps δ, k f b and k f of the bottom and upper film layers, and the numbers of film layers p u , p b . Second summand of FC describe film used for upper layers wrappings provided that the following applicability condition [14]: holds, where N denotes the set of positive integer numbers. Global number of basic film layers wrapped on the bale's lateral surface is p l = p b + p u .

Appendix A.2. Film Consumption Minimization
As the non-continuity and non-differentiability of the original film usage index FC undermines our ability to directly analytically solve the problem of film usage minimization, special attention has been given to the simpler case of near-optimal bale dimensions design, this being as important as the optimal parameters for engineering applications. The problem of the selection of near-optimal bale dimensions has been constructed by minimizing continuous lower bound of the original film usage index FC, where the goal function of a non-linear optimization problem is convex and differentiable. The necessary and sufficient optimality condition for near-optimal bale diameter, called "optimal" in [20] and here, was established in the form of a standard cubic equation, Equation (1), which can easily be solved using both analytical and numerical methods. The results of the numerical experiments [19,20] demonstrated that for any four to sixteen even layers of the film there are such compositions of bottom and upper film layers that the relative near-optimality errors do not exceed 0.01% whenever the optimal bale dimensions are used.