Mathematical Model of the Criterion of Optimization by Compensation for Designing Commercial Bottles with Lateral Surfaces of Revolution and a Straight Section along Its Silhouette

In this article, the mathematical foundations of the so called Criterion of Optimization by Compensation for designing commercial bottles with a straight section along its silhouette and with lateral surfaces of revolution is presented. Such mathematical model uses as main tools, Lagrange polynomial interpolation and Newton’s Method for Nonlinear Systems being first necessary to formulate and demonstrate a theorem. It was redesigned and manufactured a bottle of a half-liter of Fanta soda of the well-known Coca Cola Company, which uses 18.86% more material that such criterion establishes. It was expected that the redesigned bottle use 4.91% more of material with respect to what is established by the Criterion of Optimization by Compensation. However, it was reported a 13% of mistake due to important limitations that must be overcome.


Introduction
During the production stages of a bottle until being finally sent to the consumers, Open Journal of Optimization, 6, 115-147.https://doi.org/10.4236/ojop.2017.63009 the manufacturers and merchants must confront a more exigent market and society every day.The bottle has to satisfy not only the necessity of containing, protecting, preserving, commercializing and distributing merchandise, if not also, of retraining, decrease of the ecological impacts and minimization of costs.
Therefore, it is necessary to design appropriate bottles according to necessities above indicated, making this evident the necessity of generating and transmitting the knowledge of science and technology, in which the concept of optimization is of great importance.For instance, Fletcher [1] and Pierre [2] describe optimization methods that are currently most valuable in solving real-life problems.
Up to now, there has not been found studies about a clear mathematical model of an optimization method of a bottle with a lateral surface of revolution and a straight section along its silhouette.In the case of PET bottles, there are some studies about optimizing and redesigning the whole or part of the body of them by using different software programs.Masood and KeshavaMurthy [3] have reported the process, design and optimization of a bottle shape using Pro/Engineer Parametric Modelling Software and Pro/Mechanica Finite Element Software (FES); Demeril and Dave [1] have used numerical modelling with finite element analysis (FEA) techniques to redesign the petaloid base of bottles to improve stress-crack resistance.An experimental design based on an algorithmic partial cubic method was employed.Quinchung et al. [4] have optimized the structure of the PET bottle in order to increase the buckling load, based on the Abaqus/ Explicit computer program.Moreover, according to the stress contour of PET bottle obtained by Abaqus, plastic distribution of PET bottle was optimized in order to improve the efficiency of PET material and reduce the weight of the PET bottle.In addition, Mohammad [5] described the formulation of the mathematical model of the PET thermoplastic material for FEA drop-test simulation.Finally, Silva et al. [6] have designed and optimized a PET bottle through parametric computer aided design software (solidworks) and finite element method analysis, allowing the simulation of the blowing process from data input to the process variables listed in the available literature.
On the other hand, for aluminum bottles, Han et al. [7] have used numerical simulation and mathematical programming to optimized a cylindrical shell body of an aluminum can (volume: 500 ml), which was triangulated as one of the expectant choices of crushable cans for being folded down easily and safely for recycling.At the same way, Han et al. [8] have applied the structural optimization technique, to aluminum beverage bottle design, based on nonlinear finite element analyses to know the influence of the design parameters on the buckling strength and the stiffness of the bottom under an axial load and internal pressure, respectively.Similarly, Han et al. [9] have also performed multi-objective optimization of a two-piece aluminum beverage bottle considering tactile sensation of heat and embossing formability.Karl et al. [10] have performed simulation of the filling of PET bottles with a volumetric swirl chamber valve on the basis of calculations models and experiments.Hopmann et al. [11] have provide a sism tiao mulative approach to determine a well-adapted preform and bottle design as well as corresponding process parameters.For this, a three-dimensional simulation of the stretch-blow molding process was used within an iterative optimization routine.At the same way, numerical modelling and optimization of the production process of glass bottles have been the topic of several papers [12] [13] [14] [15] [16] [17].
Therefore, since most containers whose shapes are determined by surface of revolution and manufactured from plastic or metal material, the enterprises that decide to apply criterion of optimization on its manufacture, could try to reach the following achievements: reduction of rubbish, saving energy, decrease of the negative environmental impact and hence a friendly environmental image.
In this work, we report a mathematical model of the Criterion of Optimization by Compensation proposed by Reyna and Moore [18] in order to design commercial bottles with lateral surfaces of revolution and a straight section along its silhouette.With the aim of designing and manufacturing bottles using the less amount of material in their fabrication, it can avoid so environmental pollution by solid residues.Since the expenditure of material is proportional to the area of a bottle, to design a bottle by using the lesser as amount of material as possible, means that the superficial area of the bottle have to be minimized.For example, to design a cylinder of minimal area is very easy since it is possible to find its area as a function of its radius and then in using differential calculus to get the radius and height of the cylinder that makes the area of the cylinder a minimum.In the case of any container (bottle in particular) with lateral surface of revolution and a predetermined silhouette, to get a mathematical relation similar to the case of the cylinder, in general is not possible.In this case, that the Criterion of Optimization by Compensation becomes appropriate.
Such Criterion of Optimization by Compensation tells us that in order to design a bottle with a minimal total superficial area, it must first design a cylinder with a minimal total superficial area.From the cylinder, it must be get the shape of the bottle by removing certain solid parts of the cylinder.With the solid parts removed, it must be formed a cylinder whose volume is equal to the sum of the volumes of the solid parts removed from the initial cylinder.The lateral surface area of the formed cylinder must be equal to or less than the sum of the outer surface areas of the solid parts removed from the initial cylinder.This new formed cylinder has to fit the straight section of the bottle that is being designing.Since de cylinder has a minimal area, it is optimized, and at a beginning, the superficial area of the bottle is optimized too.Indeed, theorem 2.1 (see below) shows that the area of the bottle results to be less than that of the cylinder.However, this is only a descriptive fact and a mathematical modelling turns out to be necessary.Then, this work presents a mathematical modelling as the one required, according to the Criterion of Optimization by Compensation.

Theory
Flow diagram of the Optimization method for designing an optimized bottle of volume V.The starting point is a cylinder of minimal area of volume V. Following the Criterion of Optimization by Compensation according to steps 1, 2,…, where k is a finite number, the final silhouette of the optimized bottle can be found.From this, the three dimensional bottle is obtained.

Criterion of Optimization of Containers with a Straight Cylindrical Shape
To design a container with a straight cylindrical shape and of a certain volumetrically capacity in using a smaller amount of material, lead us to determine a minimum of the following function ( ) where V, A and r are the volume, total area and radius of the cylinder respectively.

Deriving the function ( )
A r with respect to r and equaling to zero we have that Solving this equation with respect to r Which is a critical point for the function given in Equation ( 1) But ( ) By replacing (2) in (3) Which is major than zero.This means that the value of r given by Equation (2) corresponds to a minimum and which let us determine the height h of the cylinder in using As it can be seen in this case, it is not difficult to get the values of r given by Equation ( 2) and h given by Equation (4) which let us minimize the area of a cylinder for a given value of its volume V. Furthermore Equation (3) tell us that Equation (1) is always concave above because means that  given by Equation (2) corresponds to a global minimum.

Criterion of Optimization by Compensation
In the case of designing containers with a non-cylindrical shape as it can be the general case, of a bottle, its silhouette could result whimsical due to aesthetic considerations as many others, so it is impossible to write an Equation (function), as in the case of the cylinder discussed in Section 2.1.Therefore, in order to solve this apparent difficult, we will use the Criterion of Optimization by Compensation proposed by Reyna and Moore [18].This optimization criterion consists in designing first a cylinder with a minimal area for a given volume in using Equations ( 2) and ( 4) and then redistributes the total volume and area that will be removed from de cylinder in order to get the desired shape of the bottle that we wish to manufacture.At the end of the process, the volume of the resultant bottle must be equal to that of the original cylinder, where the area of the resultant bottle could be less than that of the cylinder, as it is shown in theorem 2.1 below.Now, the unique proposition given in [19] tell us that the area of a cylinder of minimal area is among the area of a sphere and the area of a cube having all of them the same volume.But, the area of the bottle turns out to be less than that of the cylinder of minimal area, so we have the area of a bottle closer to the area of a sphere than the area of a cylinder of minimal area.Of course, the volumes of the bottle, cylinder of minimal area and sphere are the same.Since the sphere is considered as a geometric object with maximum volume and lesser area, the Criterion of Optimization by Compensation turns out to be a good method for optimizing.
An outline of the Criterion of Optimization by Compensation, in order to minimize the area of a bottle, keeping a constant volume is as follow: 1) Know the volumetric capacity V of the bottle that will be optimized; 2) Optimize the total area of a cylinder whose volumetric capacity is V; 3) In each zone of the optimized cylinder where will be modified (Figure 1(a)) in order to get the desired shape of the bottle, it must be removed a solid which is enclosed by an external surface, the surface of the cylinder, and internal surface, the surface of revolution that will give the shape to the bottle.See in the upper part of Figure 1  4) Fit the cylinder without bases, obtained in the step 3, in the cylindrical zone of the optimized cylinder.In this way, it is possible to obtain an optimized bottle with lateral surface of revolution (see Figure 1(c)).
The theorem below is based on the Criterion of Optimization by Compensation.This theorem tell us that it is possible to deform a cylinder of volume c V and minimal area c A to a solid of revolution of volume s V and area s A such that c s V V = and s c A A < .This theorem is formulated in order to show in par- ticular that the area of a bottle could be less than the corresponding minimal area of a cylinder, both enclosing the same volume, so that the optimization of the bottle is better.
Theorem 2.1 There exists a solid of revolution of the non-convex type whose surface is less than the surface of a right circular cylinder of minimal area, keeping both of them the same volume.

Proof
Figure 2 shows the scheme to be used in order to prove the theorem given above.
From scheme of Figure 2 we have: Due to the Criterion of Optimization by Compensation, x are values that must be freely given in such a way as to allow us to determine an acceptable value for a .

x rx Hx x x H r r h r x H r rx Hr x x x x x x a h Ha Rh a h R H H h x R H Ha Rh a h a h a h
Equation ( 54) is an equation in the unknown parameter a.
Since we have a theorem of existence, we solve numerically Equation (54) in the unknown parameter a in using the testing data Here, from Equations ( 2) and ( 4) and the value 8.86

Mathematical Model of the Criterion of Optimization by Compensation for Designing Commercial Bottles with a Straight Section along Its Silhouette
In this section, the mathematical model of the criterion of optimization by compensation for designing commercial bottles with lateral surfaces of revolution and a straight section along its silhouette is presented.The start point is the schema shown in Figure 3, which is based in the criterion of optimization by compensation seen in the xy-plane.
The mathematical model taking into account the Criterion of Optimization by Compensation is as follows: consider the schema (shown in Figure 3) corresponding to a general bottle with a straight section along its silhouette, where , are volumes generated when the regions below the curves represented by the non-constant polynomials degree rotates around the x-axis.The value of k depends on the number of polynomials needed to complete the shape of the silhouette of the bottle out of the rectangle with vertices ( ) 0, 0 , ( ) volumes generated when the regions above of the polynomials tates around the x-axis.0 V and 0 V are volumes below and above the constant polynomial ( ) ( ) Then, from Figure 3 ( ) where c V is the volume of the cylinder of minimal area with radius c r and height c h which is generated when the rectangle with vertices ( ) ( ) ( ) 0, 0 , 0, , , Putting Equation ( 58) in (57) Putting Equations ( 60)-( 66) in (59), we have finally, Equation ( 67) is a fundamental equation derived from the Criterion of Optimization by Compensation.
N o w l e t ( ) be, the polynomials that are shown in Figure 3, where the coefficients of these polynomials are unknowns.Thus, the Optimization problem can be formulated as follows: given the equation Subject to the following conditions: ( ) ( ) ( ) in Equations ( 70)-( 78) are constant ordinates given by the designer.On the other hand, in order to give the shape of the silhouette of the bottle, ( ) are also given by the designer such that each of them determines a set of points through which the polynomials obtained by interpolation, passes respectively.Similarly, the ( ) ; are also given by the designer, such that each of them determines a set of points through which the polynomials obtained by interpolation, passes respectively.These polynomials are such that their corresponding polynomials by translation according to, = fit with the remaining part of the bottle in the interval x n are constants given by the designer, while c r and c h are the radius and height of a cylinder of minimal area.We ,0 the area of the bottle of volume c V is less or equal than the area of a cylinder of minimal area according to the theorem given above.

Solving the Problem of Optimization
We can solve the optimization problem in using Lagrange polynomial interpolation and Newton's Method for Nonlinear Systems according to the following steps: Step 1 Find the constants ,0 , in using Lagrange interpolation, say, by determining the coefficients of the Lagrange interpolating polynomial in each case respectively.These coefficients are the constants we are trying to find.
Step 2 The 1 2 , and 1 h can be found by solving the system of Equa- tions (85), (86) and (87) given below, in using the Newton's Method.

c x I P x h x P x x h P x h x p r h p
where, , , , n I I I  are integrals that correspond to the polynomials , , , n P x P x P x  , and can be calculated as well as ( ) ( ) ( ) Step 3 3 n x + and1 can be found by solving the system of equations: ( ) x p r h p p I + can be calculated in using the data calculated in step 2.
Step 4 4 n x + and 3 h can be found by solving the system of equations: ( )

P x h x p r h p p p
I + can be calculated in using the data calculated in step 3.
Step k n k x + and 1 k h − can be found by solving the system of equations: ( ) I + − can be calculated in using the data calculated in step Step k+1 x + + and k h can be found by solving the system of equations: ( ) be calculated in using the data calculated in step k.
The number of unknown coefficients of the polynomials depends on the number of points that were arbitrary chosen, in order to get a Lagrange polynomial that must describe the silhouette of the bottle.Hence, we must be careful when the number of points are chosen, since the polynomial in that region could present slight oscillations and not describing the shape of the silhouette of the bottle in that region, as shown in Figure 4.If this were the case, we must try to choose other points by varying slightly some points of the initial set of points until the oscillations have vanished.Polynomial oscillation could be present when a number of chosen points are not appropriated in order to get a Lagrange polynomial interpolation.Therefore, in order to get a desirable solution for the problem of optimization and choose a set of points in such a way that the poly nomials should be free of oscillations, we must use a computer program (optimizer of five sections).

Application of the Criterion of Optimization by Compensation
As an application of the criterion by compensation, a bottle of Fanta soda of the Coca Cola Company was considered, as shown in Figure 5.The characteristics of the bottle are: Type of the product that it contains: carbonated drink; Factory: Coca Cola company; volume of the bottle: 537.5 ml; weight of the empty bottle without its cap: 25.18 g; material from which the bottle is manufactured: Polyethylene terephthalate (PET); base of the bottle: petaloid base of five cavities.
We first analyse the real bottle, in order to determine its superficial area mainly.In fact, after cutting along the bottle through the middle and putting half of it in a coordinate system, previously drawn on a millimetered paper, the x and y coordinates of the silhouette of the real bottle are shown in Table 1: Figure 6(a) shows for the real bottle, the interface of the computer program elaborated in the high-level programming language MATLAB.This computer program has five options: interpolate, volume, area, surface of revolution and optimize in order to perform the tasks.The calculations are performed by considering the bottle with a flat base.Adjustments are made in Input Data of each section.In this case, it was considered five sections  ( ) , 1, 5 i S i = for the bottle.In Optimizer, the volume of the bottle to be optimized must be entered.As a result, in using Equations ( 2)-( 4), we get optimum results for the cylinder, such as: radio, height and area.In Output Data of Figure 6(a) are shown the results of the volume and area of the bottle, according to the adjustments that were made in Input data of each section.When the optimization is being performed the value of the area varies until get approximately the value of the area obtained in Optimum results for the cylinder, keeping approximately constant the value of the volume in Optimizer (see Figure 6(a)).In using the computer program and with the data of the Table 1, the real bottle in three dimensions is shown in Figure 6(b).As, it can be seen from the results of the calculations in using the   1.It is shown the sections S1, S2…S5, according to data in Input Data of Each Section in Optimizer of Five Sections.In Optimizer, the volume V of the bottle is entered in order to obtain the optimum radius, optimum height and optimum area of the cylinder which are visualized in Optimum Results for the Bottle.In Output Data, for the bottle, the volume and area of the bottle are visualized when the curves in corresponding sections S1, S2…S5 rotates around the x-axis and which generate the bottle; (b) view of the real bottle in three dimensions.Table 1.Coordinates of the silhouette of the real bottle seen in a Cartesian coordinate system of the xy plane.Section 1 (S1), corresponds to the lip of the bottle, Section 2 (S2), section 3 (S3), section 4 (S4) and section 5 (S5) correspond to the form of the bottle.such cavities, it is found a total volume of 14 cm 3 approximately that have to be subtracted from 551.63 cm 3 .So we have 537.63 cm 3 which is a good approximation of the volume of the real bottle.Similarly, the area 435.972 cm 2 in Output Data for the bottle is due to fact that the program calculates the area of the bottle by considering it with flat base.By calculations made apart in the zone of the petaloid base in order to determine its area by approximation by triangles and circles as shown in Figure 7, we find that the area of the bottle with petaloid base is 435.23 cm 2 , so that the difference in areas between the bottle with flat base and petaloid base is not significant.
In order to redesign the bottle, it is necessary to give some points through which the silhouette of the redesigned bottle passes.These points are given by the designer and are shown in Table 2, being shown a set of them in Figure 8 2.
Table 2. Coordinates of the silhouette of the redesigned bottle seen in a Cartesian coordinate system of the xy plane.Section 1 (S1) corresponds to the lip of the bottle; Section 2 (S2), Section 3 (S3), and Section 5 (S5) correspond to the form of the bottle.P x and we will find that 1 0 h ≈ , so that the Lagrange po- lynomials of the redesigned bottle will be equal to the Lagrange polynomials of the optimized redesigned bottle.
It is worth to emphasize that, the coefficients of the polynomials, in Table 3, are expressed in the form of decimals in order to decrease distortions in the design of the bottle, optimization and during passing data to the CNC lathe.Greater the number of considered decimal places better will be the design.These criteria are taken into account in all the calculations that are made in all the process of optimization and must be presented in this way.
In order to get the optimized redesigned bottle, rest to find the polynomial ( )  .By numerical integration of the first two integrals and taking into account that ( ) , we have,

P x h x x h x h x h x h x h x h x
x h where integration was performed analytically, Equation (90) becomes finally, Now, must be solved the system of Equations ( 88), (89) and (91).Since a system of three nonlinear equations in three unknowns has been obtained, must be solved numerically in using the Newton's Method for Nonlinear Systems, whose interface Solver For System of Nonlinear Equations (SNLEs) shows us the results of the values for the unknowns after perform 5 iterations (see Figure 9).To solve the system, Equations (88), ( 89) and (91) are renamed as Additionally, we need the elements of the Jacobian matrix, which are: ) .32205962592 94.0975355844 292.840405808 In order to enter Equations ( 88), ( 89), (91) and the elements of the Jacobian matrix to the SNLEs, in Matlab language, we change to an equivalent form.Then, taking into account     and following the same procedure that was made for the real bottle, the results for the redesigned bottle, after applying criterion of optimization by compensation, are shown in Figure 10.
In Output data of this figure, we can see that we can manufacture a bottle with an area of 363.588 cm 2 less than 365.96 cm 2 specified in optimum area of optimum results for the bottle, which is in accordance with theorem 2.1 and similar results reported by Reyna and Morales [19].
When considering a petaloid base for the bottle of Figure 10, and according to the design of a new petaloid base following similar calculations to that made  4. Section 1 (S1), 2 (S2), 3 (S3), 4 (S4), and 5 (S5) of the optimized redesigned bottle with its corresponding Lagrange Polynomials that draw the silhouette of the bottle according to Table 2 for the real bottle (see Figure 6), the area of the bottle is modified from 365.96 cm 2 to 366.17 m 2 .Table 4 shows the Lagrange polynomials, of each section, of the optimized redesigned bottle.We can see similar results to that given in Table 3, which means that the optimization criterion by compensation works very well.

Results and Discussions
On the basis of the Criterion of Optimization by Compensation it was obtained a mathematical model that let us optimize (minimize) the superficial area of a bottle with a straight section along its silhouette.The minimal area that this mathematical model let us obtain is in the sense that it tends to the area of the sphere which is considered as geometric object with maximal volume and less area.Say, in general the area of the bottle is between the area of the sphere and the area of the cylinder.This is not difficult to show in using proposition presented in [11]  .In this sense, we are optimizing.We really don't know the absolute minimum for the area of the bottle being this an open problem.
Results of the application of the Criterion of Optimization by Compensation for the case of the half-liter bottle of the Fanta soda, are summarized in Table 5.It can be seen that the real bottle is made in using 25.18 g of PET plastic to enclose a volume of 537.5 cm 3 with a superficial area of 435.23 cm 2 approximately.However, in using the Criterion of Optimization by Compensation for almost the same volume (537.63 cm 3 ), the area can be reduced to 366.17 cm 2 (with petaloid base) so that it is necessary only 21.185 g of PET plastic keeping the original thickness of the wall of the bottle.It is worth to emphasize that 21.185 g was obtained using a rule of three with the data of Table 5 keeping constant the thickness of the bottle, which is not considered as a parameter of design.It was optimized the given geometric shape of the bottle.Consequently, the real bottle of Fanta soda has a mistake of 18.86 % (This value was obtained by using the re- due to the petaloid base of the bottle.The difference in areas between the bottlewith flat base and petaloid base is not significant in the present work.The cost that it is paid when the optimization is achieved in using the Criterion of Optimization by Compensation is that the height of the bottle diminishes and a widening happens.The view at scale of the real bottle and the redesigned bottle optimized in extreme is shown in Figure 11.It can be seen that the height of the redesigned bottle is almost half of the real bottle, but fatter in order to compensate the loss of volume, keeping so the same volume of the real bottle. Since no necessarily, an optimization in extreme must be achieved, it was attempted to manufacture a bottle with a mistake of 4.91% considered arbitrarily by the designer, say, with an area of 383.916 cm 2 (This value is obtained by using the relation area of the bottle area of the optimized bottle 100 4.91 area of the optimized bottle − × = , where the area of the optimized bottle is 365.96cm 2 ).Such attempt of manufacturing the bottle was achieved with certain limitations such as: 1) It was impossible to find preforms of PET plastic of 23.25 g approximately, as it was required by calculations.
2) It was impossible to find preforms of PET plastic whose design must be according to the length of the lip of the bottle as it is required in the design.
3) The mechanism of exportation of data, to the numerical control lathe distorted the data.So, volume and area of the bottle were slightly altered.
Table 6 shows coordinates of the silhouette of the manufactured bottle in a Cartesian coordinate system of the  plane.The process to obtain this Table is the same fallowed to obtain Table 2 and Table 4.
In Figure 12 we can see in Output data that the area of the bottle to be manufactured is 383.916cm 2 , while the area of the bottle in Optimum results is   6.It is shown the sections S1, S2…S7, according to data in Input Data of Each Section in Optimizer of Seven Sections.In Optimizer, the volume V of the bottle is entered in order to obtain the optimum radius, optimum height and optimum area of the cylinder which are visualized in Optimum Results for the Bottle.In Output Data, for the bottle, the volume and area of the bottle are visualized when the curves in corresponding sections S1, S2…S7 rotates around the x-axis and which generate the bottle.The coefficients of the polynomials that represent each curve of each section S1, S2…S5 are placed on the Lagrange Polynomial Coefficients (shown in Table 7).These results were obtained following the Flow diagram of the Optimization method for designing an optimized bottle of volume V; (b) view of the manufactured bottle in three dimensions.Table 6.Coordinates of the silhouette of the manufactured bottle seen in a Cartesian coordinate system of the xy plane.Section 1 (S1) corresponds to the lip of the bottle, Setion 2 (S2), Section 3 (S3)…Section 7 (S7) correspond to the form of the bottle.365.960 cm 2 being in both cases when the base of the bottle is flat.This means that there is a mistake of 4.91% with respect to optimum results.As a consequence of inserting to the bottle a petaloid base plus the limitations stated above the mistake increase as it is indicated below.Table 7 shows the Lagrange polynomials, of each section, of the manufactured bottle which were obtained following the same procedure to obtain Table 4.
The above last limitation (number 3) allowed us to obtain a bottle whose geometric design had a mistake of 13% (This value was obtained by direct measurement of the area of the mold of the bottle and by using the formula area of the mold of the bottle area of the optimized bottle 100 area of the optimized bottle − × , where the area of the optimized bottle is 365.96cm 2 according to Optimum results for the bottle of Figure 12).The limitation, numbered as 1, caused that the thickness of the wall of the bottle resulted be 75% of that of the real bottle because a preform of 22 g of PET plastic was used.In spite of the limitations, a good bottle was obtained as shown in Figure 13.In this Figure we can see: a) the corresponding mold that was fabricated to manufacture the redesigned Fanta soda bottle of the Coca Cola company, b) comparison between the real (higher) and redesigned (smaller) bottle.Our results are in accordance with the literature, for instance, Silva et al. [6] have reported a reduction of 21% of PET material (weight of 4.6 g) of a bottle of 22 g, using simulations based on finite element method (FEM).On the other hand, Hung et al. [20] have reported a reduction of the weight of a bottle, with a volume of 500 ml, from 27 g to 22.3 g using numerical simulations, where thickness and pattern of the bottle were changed in order to reduce the weight of the bottle.At the same way, Hopmann et al. [11] have also reported a reduction of the weight of a 0.5 liter PET bottle of Krones AG, Neutraubling, Germany, from 18.5 g to 15.5 g.For this aim, they have developed a simulative approach to determine a well-adapted preform and bottle design with its corresponding process parameters.Therefore, a three-dimensional simulation of the stretch-blow molding process was used within an iterative optimization routine.Their bottles were manufactured on a stretch-blow molding machine LB1 (Krones AG, Neutraubling, Germany), with a changed wall thickness.Therefore, according to this reports, we can note that, in all cases the wall Table 7. Lagrange polynomials of the manufactured bottle obtained via interpolation according to data of Table 6 and in using Optimizer of seven sections.thickness of the bottles was changed, and the weight of the bottle reduced in av-erage~4.0 g, similarly to the found results in this work.However, in this work the wall thickness of the bottle have to remain unchanged.

Conclusion
In summary, we propose the mathematical foundations of the so-called Criterion of Optimization by Compensation and in using this, an optimized and redesigned half-liter bottle of Fanta soda of the well-known Coca Cola Company, was manufactured.This manufactured bottle was designed with a mistake of 4.91% with respect to what such criterion of optimization stablishes.However, it was reported a mistake of 13% in each manufactured bottle due to important technical limitations listed above that must be overcome.In general, in spite of such limitations a good bottle was obtained with such 13% of mistake, resulting the thickness of the wall of the bottle 75% of that of the real bottle because a preform of 22 g of PET plastic was used, instead of a preform of 23.25 g as it was required.
How to cite this paper: Zegarra, L.B.R., Armas, L.E.G.Reyna, A.D., Vergara, J.A.L. and Obeso, F.A.V. (2017) Mathematical Model of the Criterion of Optimization by Compensation for Designing Commercial Bottles with Lateral Surfaces of Revolution and a Straight Section along Its Silhouette.
(b) signed with the arrow.Calculate the difference of areas between the cylindrical surface and the surface of revolution that both enclose the solid that was removed.With this difference must be formed a cylinder without bases whose volume must be equal approximately to the volume that was removed (see Figure 1(b));

Figure 1 .
Figure 1.(a) Optimized cylinder with a volume V from which an optimized bottle will be obtained; (b) Cylinder with an upper region that will be lost and recovered in the lower region; (c) Bottle after the application of the optimization criterion.Adapted with permission from Reyna and Moore [18].Copyright 2016 ECI.

Figure 2 .
Figure 2. Schema according to the criterion of optimization by compensation.c h and and R are given heuristically by looking Figure 2 and ac- cording to the Criterion of Optimization by Compensation.This let us find 15.1384 a u =.With this data, the area and volume of the solid of revolution that is generated when the region enclosed by the polygonal OABGJMN and x-axis of Figure2, rotates around x-axis, are equal to minimal area of the cylinder of revolution that is generated when the region enclosed by the rectangle ODIK rotates around the x-axis, is equal to

Figure 3 .
Figure 3. Schema to get a mathematical model in using the Criterion of Optimization by Compensation.When the rectangle with vertices ( ) ( ) ( ) 0, 0 , 0, , , c c c r x-axis (see Figure3).By Criterion of Optimization by Compensation ( )

Figure 4 .
Figure 4. Schema showing the polynomial oscillations when the chosen points are not appropriated to get the Lagrange interpolation.This happens when the number of points chosen is great.

Figure 5 .
Figure 5. Bottle of Fanta soda, from the Coca-Cola Company, used to apply the application of the Criterion of Optimization by Compensation.

Figure 6 .
Figure 6.Real bottle, (a) Silhouette of the real bottle according to the data of theTable1.It is shown the sections S1, S2…S5, according to data in Input Data of Each Section in Optimizer of Five Sections.In Optimizer, the volume V of the bottle is entered in order to obtain the optimum radius, optimum height and optimum area of the cylinder which are visualized in Optimum Results for the Bottle.In Output Data, for the bottle, the volume and area of the bottle are visualized when the curves in corresponding sections S1, S2…S5 rotates around the x-axis and which generate the bottle; (b) view of the real bottle in three dimensions.

Figure 6 (
Figure 6(a), for the bottle exceeds 537.5 cm 3 because the program calculates the volume of the bottle with flat base.By measuring the volume of the cavities of the petaloid base (base of the real bottle) in laboratory by using water for filling

4 P 1 h
by circles in sections S1, S2, S3, S4 and S5.The task is to determine by Lagrange interpolation the polynomials that pass through these points.The number 4x in the x-axis, is the right extreme of the section S4 which is a straight line, and left extreme of section S5 represented by the polynomial ( ) x .This number must be determined by moving section S5 to the right or left according to the sign of , positive or negative in ( )4 1 P x h −as it is shown in Figure8.These polynomials are the starting points to get the final polynomials that determine the silhouette of the redesigned bottle and hence optimized.As it can be seen

Figure 7 .
Figure 7. Approximation by triangles and circles in order to determine the area of the petaloid base of the bottle.

Figure 9 .
Figure 9. Solver for System of Nonlinear Equations based on Newton-Raphson numerical method to find the unknowns parameters of Equations.(88), (89) and (91).These Equations were renamed as 1 f , 2 f and 3 f (see Enter F).This figure also show the

Figure 9
Figure9shows the SNLEs, where it is possible to see the entry functions 1 f ,

Figure 10 .
Figure 10.Optimized redesigned bottle with flat base, (a) Silhouette of the optimized redesigned bottle according to the data of Table2.It is shown the sections S1, S2… S5, according to data in Input Data of Each Section in Optimizer of five Sections.In Optimizer, the volume V of the bottle is entered in order to obtain the optimum radius, optimum height and optimum area of the cylinder which are visualized in Optimum Results for the Bottle.In Output Data for the bottle the volume and area of the bottle are visualized when the curves in corresponding sections S1, S2…S5 rotates around the x-axis and which generate the bottle.The coefficients of the polynomials that represent each curve of each section S1, S2…S5 are placed on the Lagrange Polynomial Coefficients (shown in Table4).Theseresults were obtained following the Flow diagram of the Optimization method for designing an optimized bottle of volume V; (b) view of the optimized redesigned bottle in three dimensions.

Figure 11 .
Figure 11.View at scale of the real and redesigned bottle.Adapted with permission from Reyna and Moore[18].Copyright 2016 ECI.

Figure 12 .
Figure 12.Manufactured bottle, (a) Silhouette of the manufactured bottle according to the data of Table6.It is shown the sections S1, S2…S7, according to data in Input Data of Each Section in Optimizer of Seven Sections.In Optimizer, the volume V of the bottle is entered in order to obtain the optimum radius, optimum height and optimum area of the cylinder which are visualized in Optimum Results for the Bottle.In Output Data, for the bottle, the volume and area of the bottle are visualized when the curves in corresponding sections S1, S2…S7 rotates around the x-axis and which generate the bottle.The coefficients of the polynomials that represent each curve of each section S1, S2…S5 are placed on the Lagrange Polynomial Coefficients (shown in Table7).These results were obtained following the Flow diagram of the Optimization method for designing an optimized bottle of volume V; (b) view of the manufactured bottle in three dimensions.

Figure 13 .
Figure 13.(a) Mold used to manufacture the redesigned bottle of Fanta soda, and (b) comparison between the real (higher) and redesigned (smaller) bottle.

Table 3 .
Polynomials according to the points given inTable 2 corresponding to Section 1

Table 3 .
With these results, and in using the Optimizer of Five Sections, we get the polynomials for the optimized redesigned bottle shown in Lagrange polynomial coefficients of Figure10, and which are shown in Table4.Therefore, with the data of Table2where 5 12.007 12.01 that establishes the inequality sphere

Table 5 .
Comparison between the results of the real and the optimized bottle.