﻿ 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

Open Journal of Optimization
Vol.06 No.03(2017), Article ID:79317,33 pages
10.4236/ojop.2017.63009

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

Lizandro B. R. Zegarra1*, Luis E. G. Armas2**, Aida D. Reyna3, Julio A. L. Vergara1, Fidel A. V. Obeso1

1Departamento de Matemática de la Universidad Nacional del Santa, Chimbote, Perú

2Grupo de Óptica Micro e Nanofabricação de Dispositivos (GOMNDI), Universidade Federal do Pampa, Campus Alegrete, Alegrete, Brazil

3Departamento de Matemática de la Universidad Nacional de Trujillo, Trujillo, Perú

Received: August 5, 2017; Accepted: September 23, 2017; Published: September 26, 2017

ABSTRACT

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.

Keywords:

Optimization, Modelling, Bottle Design, Environmental Pollution, Solid Residues

1. Introduction

During the production stages of a bottle until being finally sent to the consumers, 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.

2. 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,…, $s=3,k+1$ , where k is a finite number, the final silhouette of the optimized bottle can be found. From this, the three dimensional bottle is obtained.

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

$A\left(r\right)=2\text{π}{r}^{2}+\frac{2V}{r}$ (1)

where V, A and r are the volume, total area and radius of the cylinder respectively.

Deriving the function $A\left(r\right)$ with respect to r and equaling to zero we have that

$4\text{π}r-\frac{2V}{{r}^{2}}=0$

Solving this equation with respect to r

$r=\sqrt[3]{\frac{V}{2\text{π}}}$ (2)

Which is a critical point for the function given in Equation (1)

But

$\frac{{\text{d}}^{2}A\left(r\right)}{\text{d}{r}^{2}}=4\text{π}+\frac{4V}{{r}^{3}}$ (3)

By replacing (2) in (3)

$\frac{{\text{d}}^{2}A\left(\sqrt[3]{\frac{V}{2\text{π}}}\right)}{\text{d}{r}^{2}}=12\text{π}$

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 $V=\text{π}{r}^{2}h$ , equal to

$h=\sqrt[3]{\frac{4V}{\text{π}}}$ (4)

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 $4\text{π}+\frac{4V}{{r}^{3}}>0$ , $\forall r>0$ . This

means that given by Equation (2) corresponds to a global minimum.

2.2. 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(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));

(a) (b) (c)

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.

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 ${V}_{c}$ and minimal area ${A}_{c}$ to a solid of revolution of volume ${V}_{s}$ and area ${A}_{s}$ such that ${V}_{c}={V}_{s}$ and ${A}_{s}<{A}_{c}$ . This theorem is formulated in order to show in particular 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:

${V}_{c}={\sum }_{i=1}^{6}{V}_{i}$ (5)

Due to the Criterion of Optimization by Compensation,

${V}_{7}={\sum }_{i=4}^{6}{V}_{i}$ (6)

Figure 2. Schema according to the criterion of optimization by compensation. ${h}_{c}$ and ${r}_{c}$ correspond to the height and radius of a cylinder of volume ${V}_{c}$ and minimal area ${A}_{c}$ which is generated when the rectangle ODIK rotates around the x-axis. From Equations (2) and (4), we have that ${r}_{c}=\frac{{h}_{c}}{2}$ . $P\left(x\right)$ and $Q\left(x\right)$ are polynomials of first degree whose coefficients must be found. ${V}_{i},i=\stackrel{^}{1,7}$ are the volumes generated when the regions given by the correspondent closed polygons OABC, CBGT, TGJK, GFIJ, BEFG, ADEB and KJMN, rotate around the x-axis. ${A}_{i}$ , $i=\stackrel{^}{0,5}$ are the areas generated when the correspondent segments OA, AB, BG, GJ, JM and MN, rotate around the x-axis. The value of $a$ must be found. $H,r,R,{x}_{0}$ and ${x}_{1}$ are values that must be freely given in such a way as to allow us to determine an acceptable value for $a$ .

replacing Equation (6) in (5)

${V}_{c}={\sum }_{i=1}^{3}{V}_{i}+{V}_{7}$ (7)

We have as well that,

${V}_{c}=\text{π}{r}_{c}^{2}{h}_{c}$ (8)

${V}_{1}=\text{π}{r}^{2}{x}_{0}$ (9)

${V}_{2}=\text{π}{\int }_{{x}_{0}}^{{x}_{1}}{\left[P\left(x\right)\right]}^{2}\text{d}x$ (10)

${V}_{3}=\text{π}{H}^{2}\left({h}_{c}-{x}_{1}\right)$ (11)

${V}_{7}=\text{π}{\int }_{{h}_{c}}^{a}{\left[Q\left(x\right)\right]}^{2}\text{d}x$ (12)

$P\left(x\right)={\alpha }_{1}x+{\alpha }_{2}$ (13)

$Q\left(x\right)={\beta }_{1}x+{\beta }_{2}$ (14)

By replacing Equation (13) in (10)

${V}_{2}=\text{π}{\int }_{{x}_{0}}^{{x}_{1}}{\left[{\alpha }_{1}x+{\alpha }_{2}\right]}^{2}\text{d}x$

${V}_{2}=\text{π}\left[{\alpha }_{1}^{2}\frac{\left({x}_{1}^{3}-{x}_{0}^{3}\right)}{3}+{\alpha }_{1}{\alpha }_{2}\left({x}_{1}^{2}-{x}_{0}^{2}\right)+{\alpha }_{2}^{2}\left({x}_{1}-{x}_{0}\right)\right]$ (15)

Then,

${V}_{2}={V}_{2}\left({\alpha }_{1},{\alpha }_{2}\right)$ (16)

By replacing Equation (14) in (12),

${V}_{7}=\text{π}{\int }_{{h}_{c}}^{a}{\left[{\beta }_{1}x+{\beta }_{2}\right]}^{2}\text{d}x$

${V}_{7}=\text{π}\left[{\beta }_{1}^{2}\frac{\left({a}^{3}-{h}_{c}^{3}\right)}{3}+{\beta }_{1}{\beta }_{2}\left({a}^{2}-{h}_{c}^{2}\right)+{\beta }_{2}^{2}\left(a-{h}_{c}\right)\right]$ (17)

Then,

${V}_{7}={V}_{7}\left({\beta }_{1},{\beta }_{2},a\right)$ (18)

By replacing Equations. (8), (9), (11), (16) and (18) in (7),

$\text{π}{r}_{c}^{2}{h}_{c}-\text{π}{r}^{2}{x}_{0}-{V}_{2}\left({\alpha }_{1},{\alpha }_{2}\right)-\text{π}{H}^{2}\left({h}_{c}-{x}_{1}\right)-{V}_{7}\left({\beta }_{1},{\beta }_{2},a\right)=0$ (19)

Let ${A}_{c}$ be such that,

${A}_{c}={\sum }_{i=0}^{5}{A}_{i}$ (20)

Where,

${A}_{0}=\text{π}{r}^{2}$ (21)

${A}_{1}=2\text{π}r{x}_{0}$ (22)

${A}_{2}=2\text{π}{\int }_{{x}_{0}}^{{x}_{1}}P\left(x\right)\sqrt{1+{\left[{P}^{\prime }\left(x\right)\right]}^{2}}\text{d}x$ (23)

${A}_{3}=2\text{π}H\left({h}_{c}-{x}_{1}\right)$ (24)

${A}_{4}=2\text{π}{\int }_{{h}_{c}}^{a}Q\left(x\right)\sqrt{1+{\left[{Q}^{\prime }\left(x\right)\right]}^{2}}\text{d}x$ (25)

${A}_{5}=\text{π}{R}^{2}$ (26)

${A}_{c}=2\text{π}{r}_{c}{h}_{c}+2\text{π}{r}_{c}^{2}$ (27)

Deriving Equation (13),

${P}^{\prime }\left(x\right)={\alpha }_{1}$ (28)

Deriving Equation (14),

${Q}^{\prime }\left(x\right)={\beta }_{1}$ (29)

Replacing Equations (13) and (28) in (23) and integrating,

${A}_{2}=2\text{π}\sqrt{1+{\alpha }_{1}^{2}}\left[{\alpha }_{1}\frac{{x}_{1}^{2}}{2}+{\alpha }_{2}{x}_{1}-{\alpha }_{1}\frac{{x}_{0}^{2}}{2}-{\alpha }_{2}{x}_{0}\right]$ (30)

Then,

${A}_{2}={A}_{2}\left({\alpha }_{1},{\alpha }_{2}\right)$ (31)

Replacing Equations (14) and (29) in (25) and integrating,

${A}_{4}=2\text{π}\sqrt{1+{\beta }_{1}^{2}}\left[{\beta }_{1}\frac{{a}^{2}}{2}+{\beta }_{2}a-{\beta }_{1}\frac{{h}_{c}^{2}}{2}-{\beta }_{2}{h}_{c}\right]$ (32)

Then,

${A}_{4}={A}_{4}\left({\beta }_{1},{\beta }_{2},a\right)$ (33)

Replacing Equations (21), (22), (24), (26), (31) and (33) in (20),

$2\text{π}{r}_{c}{h}_{c}+2\text{π}{r}_{c}^{2}-\text{π}{r}^{2}-2\text{π}r{x}_{0}-{A}_{2}\left({\alpha }_{1},{\alpha }_{2}\right)-2\text{π}H\left({h}_{c}-{x}_{1}\right)-{A}_{4}\left({\beta }_{1},{\beta }_{2},a\right)-\text{π}{R}^{2}=0$ (34)

Now, the conditions that the polynomials $P\left(x\right)$ and $Q\left(x\right)$ must satisfy are:

$P\left({x}_{0}\right)=r$ (35)

$P\left({x}_{1}\right)=H$ (36)

$Q\left({h}_{c}\right)=H$ (37)

$Q\left(a\right)=R$ (38)

From Equations (35) and (13),

${\alpha }_{1}{x}_{0}+{\alpha }_{2}-r=0$ (39)

From Equations (36) and (13),

${\alpha }_{1}{x}_{1}+{\alpha }_{2}-H=0$ (40)

From Equations (37) and (14),

${\beta }_{1}{h}_{c}+{\beta }_{2}-H=0$ (41)

From Equations (38) and (14),

${\beta }_{1}a+{\beta }_{2}-R=0$ (42)

Taking into account Equations (19), (34), (39), (40), (41) and (42) we get the following system of equations,

${\alpha }_{1}{x}_{0}+{\alpha }_{2}-r=0$ (43)

${\alpha }_{1}{x}_{1}+{\alpha }_{2}-H=0$ (44)

${\beta }_{1}{h}_{c}+{\beta }_{2}-H=0$ (45)

${\beta }_{1}a+{\beta }_{2}-R=0$ (46)

$\text{π}{r}_{c}^{2}{h}_{c}-\text{π}{r}^{2}{x}_{0}-{V}_{2}\left({\alpha }_{1},{\alpha }_{2}\right)-\text{π}{H}^{2}\left({h}_{c}-{x}_{1}\right)-{V}_{7}\left({\beta }_{1},{\beta }_{2},a\right)=0$ (47)

$\begin{array}{l}2\text{π}{r}_{c}{h}_{c}+2\text{π}{r}_{c}^{2}-\text{π}{r}^{2}-2\text{π}r{x}_{0}-{A}_{2}\left({\alpha }_{1},{\alpha }_{2}\right)\\ -2\text{π}H\left({h}_{c}-{x}_{1}\right)-{A}_{4}\left({\beta }_{1},{\beta }_{2},a\right)-\text{π}{R}^{2}=0\end{array}$ (48)

Hence, we have six equations in five unknowns.

From Equations (43) and (44),

${\alpha }_{1}=\frac{H-r}{{x}_{1}-{x}_{0}}$ (49)

${\alpha }_{2}=\frac{r{x}_{1}-H{x}_{0}}{{x}_{1}-{x}_{0}}$ (50)

From Equations. (45) and (46),

${\beta }_{1}=\frac{R-H}{a-{h}_{c}}$ (51)

${\beta }_{2}=\frac{Ha-R{h}_{c}}{a-{h}_{c}}$ (52)

By replacing Equations (15) and (17) in (47),

$\begin{array}{l}\text{π}{r}_{c}^{2}{h}_{c}-\text{π}{r}^{2}{x}_{0}-\text{π}\left[{\alpha }_{1}^{2}\frac{\left({x}_{1}^{3}-{x}_{0}^{3}\right)}{3}+{\alpha }_{1}{\alpha }_{2}\left({x}_{1}^{2}-{x}_{0}^{2}\right)+{\alpha }_{2}^{2}\left({x}_{1}-{x}_{0}\right)\right]\\ -\text{\hspace{0.17em}}\text{π}{H}^{2}\left({h}_{c}-{x}_{1}\right)-\text{π}\left[{\beta }_{1}^{2}\frac{\left({a}^{3}-{h}_{c}^{3}\right)}{3}+{\beta }_{1}{\beta }_{2}\left({a}^{2}-{h}_{c}^{2}\right)+{\beta }_{2}^{2}\left(a-{h}_{c}\right)\right]=0\end{array}$ (53)

By replacing Equations (49), (50), (51) and (52) in (53) and simplifying,

$\begin{array}{l}\text{π}{r}_{c}^{2}{h}_{c}-\text{π}{r}^{2}{x}_{0}-\text{π}\left[{\left(\frac{H-r}{{x}_{1}-{x}_{0}}\right)}^{2}\frac{\left({x}_{1}^{3}-{x}_{0}^{3}\right)}{3}+\left(H-r\right)\left(r{x}_{1}-H{r}_{0}\right)\left(\frac{{x}_{1}+{x}_{0}}{{x}_{1}-{x}_{0}}\right)+\frac{{\left(r{x}_{1}-H{x}_{0}\right)}^{2}}{{x}_{1}-{x}_{0}}\right]\\ -\text{\hspace{0.17em}}\text{π}{H}^{2}\left({h}_{c}-{x}_{1}\right)-\text{π}\left[{\left(\frac{R-H}{a-{h}_{c}}\right)}^{2}\frac{\left({a}^{3}-{h}_{c}^{3}\right)}{3}+\left(R-H\right)\left(Ha-R{h}_{c}\right)\left(\frac{a+{h}_{c}}{a-{h}_{c}}\right)+\frac{{\left(Ha-R{h}_{c}\right)}^{2}}{a-{h}_{c}}\right]=0\end{array}$ (54)

Equation (54) is an equation in the unknown parameter a.

By replacing Equations (30) and (32) in (48)

$\begin{array}{l}2\text{π}{r}_{c}{h}_{c}+2\text{π}{r}_{c}^{2}-\text{π}{r}^{2}-2\text{π}r{x}_{0}-2\text{π}\sqrt{1+{\alpha }_{1}^{2}}\left[{\alpha }_{1}\frac{{x}_{1}^{2}}{2}+{\alpha }_{2}{x}_{1}-{\alpha }_{1}\frac{{x}_{0}^{2}}{2}-{\alpha }_{2}{x}_{0}\right]\\ -\text{\hspace{0.17em}}2\text{π}H\left({h}_{c}-{x}_{1}\right)-2\text{π}\sqrt{1+{\beta }_{1}^{2}}\left[{\beta }_{1}\frac{{a}^{2}}{2}+{\beta }_{2}a-{\beta }_{1}\frac{{h}_{c}^{2}}{2}-{\beta }_{2}{h}_{c}\right]-\text{π}{R}^{2}=0\end{array}$ (55)

By replacing Equations (49), (50), (51) and (52) in (55),

$\begin{array}{l}2\text{π}{r}_{c}{h}_{c}+2\text{π}{r}_{c}^{2}-\text{π}{r}^{2}-2\text{π}r{x}_{0}-2\text{π}\sqrt{1+{\left(\frac{H-r}{{x}_{1}-{x}_{0}}\right)}^{2}}\left[\left(\frac{H-r}{{x}_{1}-{x}_{0}}\right)\frac{{x}_{1}^{2}}{2}+\left(\frac{r{x}_{1}-H{x}_{0}}{{x}_{1}-{x}_{0}}\right){x}_{1}\\ -\left(\frac{H-r}{{x}_{1}-{x}_{0}}\right)\frac{{x}_{0}^{2}}{2}-\left(\frac{r{x}_{1}-H{x}_{0}}{{x}_{1}-{x}_{0}}\right){x}_{0}\right]-2\text{π}H\left({h}_{c}-{x}_{1}\right)-2\text{π}\sqrt{1+{\left(\frac{R-H}{a-{h}_{c}}\right)}^{2}}\\ ×\left[\left(\frac{R-H}{a-{h}_{c}}\right)\frac{{a}^{2}}{2}+\left(\frac{Ha-R{h}_{c}}{a-{h}_{c}}\right)a-\left(\frac{R-H}{a-{h}_{c}}\right)\frac{{h}_{c}^{2}}{2}-\left(\frac{Ha-R{h}_{c}}{a-{h}_{c}}\right){h}_{c}\right]-\text{π}{R}^{2}=0\end{array}$ (56)

Equation (56) as well as 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 $r=1.2u$ , ${r}_{c}=4.43u$ , ${h}_{c}=8.86u$ , ${x}_{0}=2u$ , $H=4.2u$ , $R=2.7u$ and ${x}_{1}=4.7u$ . Here, from Equa-

tions (2) and (4) ${r}_{c}=\frac{{h}_{c}}{2}$ and the value ${h}_{c}=8.86u$ is arbitrary. The rest of

values for $r,{x}_{0},H$ and $R$ are given heuristically by looking Figure 2 and according to the Criterion of Optimization by Compensation. This let us find $a=15.1384u$ . 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 Figure 2, rotates around x-axis, are equal to ${A}_{s}=360.684{u}^{2}$ and ${V}_{s}=546.249{u}^{3}$ respectively.

The 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 ${A}_{c}=369.921{u}^{2}$ , while its volume is equal to ${V}_{c}=546.249{u}^{3}$ which is the same value for the solid of revolution. So, we have ${V}_{c}={V}_{s}$ and ${A}_{s}<{A}_{c}$ which prove the theorem.

Equation (56) led us to a no wished solution since $a=15.5758u$ increase the volume of the solid of revolution to ${V}_{s}=562.862{u}^{3}$ such that the surface of minimal area of the cylinder is ${A}_{c}={A}_{s}=369.921{u}^{2}$ . Say, we have ${A}_{c}={A}_{s}$ and ${V}_{c}<{V}_{s}$ .

2.3. 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 ${V}_{i},i=\stackrel{^}{1,n+k}$ , are volumes generated when the regions below the curves represented by the non-constant polynomials ${P}_{i}\left(x\right),i=\stackrel{^}{1,n+k}$ of certain

Figure 3. Schema to get a mathematical model in using the Criterion of Optimization by Compensation. When the rectangle with vertices $\left(0,0\right),\left(0,{r}_{c}\right),\left({h}_{c},{r}_{c}\right)$ and $\left({h}_{c},0\right)$ rotates around the x-axis, the cylinder of minimal area is to be obtained. When the curve described by the polynomials that must be determined in using interpolation when a set of points is given, rotates around the x-axis, the bottle must be obtained. Here, $r,{r}_{c},{h}_{c}$ and ${x}_{i},i=\stackrel{^}{1,n}$ are known constants. ${P}_{i}\left(x\right),i=\stackrel{^}{0,n+k}$ are the polynomials that describe the silhouette of the bottle. ${V}_{i}$ and ${V}^{i},\stackrel{^}{0,n}$ are described on the text.

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 $\left(0,0\right)$ , $\left(0,{r}_{c}\right)$ , $\left({h}_{c},{r}_{c}\right)$ and $\left({h}_{c},0\right)$ . ${V}^{i},\stackrel{^}{1,n}$ are volumes generated when the regions above of the polynomials ${P}_{i}\left(x\right),\stackrel{^}{1,n}$ , rotates around the x-axis. ${V}_{0}$ and ${V}^{0}$ are volumes below and above the constant polynomial $P\left(x\right)={P}_{0}\left({x}_{n}\right)$ . Then, from Figure 3

${V}_{c}={\sum }_{i=0}^{n}\left({V}_{i}+{V}^{i}\right)$ (57)

where ${V}_{c}$ is the volume of the cylinder of minimal area with radius ${r}_{c}$ and height ${h}_{c}$ which is generated when the rectangle with vertices $\left(0,0\right),\left(0,{r}_{c}\right),\left({h}_{c},{r}_{c}\right)$ and $\left({h}_{c},0\right)$ rotates around the x-axis (see Figure 3). By Criterion of Optimization by Compensation

${\sum }_{i=0}^{n}{V}^{i}=\text{π}{\left({P}_{0}\left({x}_{n}\right)\right)}^{2}\left({x}_{n+1}-{h}_{c}\right)+{\sum }_{i=1}^{k}{V}_{n+i}$ (58)

Putting Equation (58) in (57)

${V}_{c}={\sum }_{i=0}^{n}{V}_{i}+\text{π}{\left({P}_{0}\left({x}_{n}\right)\right)}^{2}\left({x}_{n+1}-{h}_{c}\right)+{\sum }_{i=1}^{k}{V}_{n+i}$ (59)

But,

${V}_{c}=\text{π}{r}_{c}^{2}{h}_{c}$ (60)

${V}_{1}=\text{π}{\int }_{0}^{{x}_{1}}{\left[{P}_{1}\left(x\right)\right]}^{2}\text{d}x$ (61)

${V}_{2}=\text{π}{\int }_{{x}_{1}}^{{x}_{2}}{\left[{P}_{2}\left(x\right)\right]}^{2}\text{d}x$ (62)

$⋮$

${V}_{n}=\text{π}{\int }_{{x}_{n-1}}^{{x}_{n}}{\left[{P}_{n}\left(x\right)\right]}^{2}\text{d}x$ (63)

${V}_{0}=\text{π}{\left[{P}_{0}\left({x}_{n}\right)\right]}^{2}\left({h}_{c}-{x}_{n}\right)$ (64)

$\begin{array}{l}{V}_{n+1}=\text{π}{\int }_{{x}_{n+1}}^{{x}_{n+2}}{\left[{P}_{n+1}\left(x\right)\right]}^{2}\text{d}x\\ {V}_{n+2}=\text{π}{\int }_{{x}_{n+2}}^{{x}_{n+3}}{\left[{P}_{n+2}\left(x\right)\right]}^{2}\text{d}x\end{array}$ (65)

$⋮$

${V}_{n+k}=\text{π}{\int }_{{x}_{n+k}}^{{x}_{n+k+1}}{\left[{P}_{n+k}\left(x\right)\right]}^{2}\text{d}x$ (66)

Putting Equations (60)?(66) in (59), we have finally,

$\begin{array}{c}\text{π}{r}_{c}^{2}{h}_{c}=\text{π}{\int }_{0}^{{x}_{1}}{\left[{P}_{1}\left(x\right)\right]}^{2}\text{d}x+\text{π}{\sum }_{i=1}^{n-1}{\int }_{{x}_{i}}^{{x}_{i+1}}{\left[{P}_{i+1}\left(x\right)\right]}^{2}\text{d}x+\text{π}{\left[{P}_{0}\left({x}_{n}\right)\right]}^{2}\left({h}_{c}-{x}_{n}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\text{π}{\left({P}_{0}\left({x}_{n}\right)\right)}^{2}\left({x}_{n+1}-{h}_{c}\right)+\text{π}{\sum }_{i=1}^{k}{\int }_{{x}_{n+i}}^{{x}_{n+i+1}}{\left[{P}_{n+i}\left(x\right)\right]}^{2}\text{d}x\end{array}$ (67)

Equation (67) is a fundamental equation derived from the Criterion of Optimization by Compensation.

Now let ${P}_{1}\left(x\right)={\sum }_{i=0}^{{m}_{1}}{a}_{1}^{i}{x}^{i}$ , ${P}_{2}\left(x\right)={\sum }_{i=0}^{{m}_{2}}{a}_{2}^{i}{x}^{i}$ , $\cdots$ , ${P}_{n}\left(x\right)={\sum }_{i=0}^{{m}_{n}}{a}_{n}^{i}{x}^{i}$ , , ${P}_{n+2}\left(x\right)={\sum }_{i=0}^{{m}_{n+2}}{a}_{n+2}^{i}{x}^{i}$ , $\cdots$ , ${P}_{n+k}\left(x\right)={\sum }_{i=0}^{{m}_{n+k}}{a}_{n+k}^{i}{x}^{i}$ 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

$\begin{array}{c}\text{π}{r}_{c}^{2}{h}_{c}=\text{π}{\int }_{0}^{{x}_{1}}{\left[{P}_{1}\left(x\right)\right]}^{2}\text{d}x+\text{π}{\sum }_{i=1}^{n-1}{\int }_{{x}_{i}}^{{x}_{i+1}}{\left[{P}_{i+1}\left(x\right)\right]}^{2}\text{d}x+\text{π}{\left[{P}_{0}\left({x}_{n}\right)\right]}^{2}\left({h}_{c}-{x}_{n}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\text{π}{\left({P}_{0}\left({x}_{n}\right)\right)}^{2}\left({x}_{n+1}-{h}_{c}\right)+\text{π}{\sum }_{i=1}^{k}{\int }_{{x}_{n+i}}^{{x}_{n+i+1}}{\left[{P}_{n+i}\left(x\right)\right]}^{2}\text{d}x\end{array}$ (68)

Subject to the following conditions:

${P}_{1}\left({x}_{0}\right)=r$ (69)

${P}_{1}\left({x}_{1}\right)={P}_{2}\left({x}_{1}\right)={R}_{1}$ (70)

${P}_{2}\left({x}_{2}\right)={P}_{3}\left({x}_{2}\right)={R}_{2}$ (71)

${P}_{n-1}\left({x}_{n-1}\right)={P}_{n}\left({x}_{n-1}\right)={R}_{n-1}$

${P}_{n}\left({x}_{n}\right)={P}_{0}\left({x}_{n}\right)={R}_{n}={R}_{0}$ (72)

${P}_{0}\left({x}_{n}\right)={P}_{0}\left({x}_{n+1}\right)={P}_{n+1}\left({x}_{n+1}-{h}_{1}\right)={R}_{n}={R}_{0}={R}_{n+1}$ (73)

${P}_{n+1}\left({x}_{n+2}-{h}_{1}\right)={P}_{n+2}\left({x}_{n+2}-{h}_{2}\right)={R}_{n+2}$ (74)

${P}_{n+2}\left({x}_{n+3}-{h}_{2}\right)={P}_{n+3}\left({x}_{n+3}-{h}_{3}\right)={R}_{n+3}$ (75)

${P}_{n+3}\left({x}_{n+4}-{h}_{3}\right)={P}_{n+4}\left({x}_{n+4}-{h}_{4}\right)={R}_{n+4}$ (76)

$⋮$

${P}_{n+k-1}\left({x}_{n+k}-{h}_{k-1}\right)={P}_{n+k}\left({x}_{n+k}-{h}_{k}\right)={R}_{n+k}$ (77)

${P}_{n+k}\left({x}_{n+k+1}-{h}_{k}\right)={R}_{n+k+1}$ (78)

${P}_{1}\left({t}_{i}\right)={P}_{1}^{i},\text{\hspace{0.17em}}0\le {t}_{i}\le {x}_{1},\text{\hspace{0.17em}}0\le i\le {m}_{1}$ (79)

${P}_{2}\left({t}_{i}\right)={P}_{2}^{i},\text{\hspace{0.17em}}{x}_{1}\le {t}_{i}\le {x}_{2},\text{\hspace{0.17em}}0\le i\le {m}_{2}$ (80)

$⋮$

${P}_{n}\left({t}_{i}\right)={P}_{n}^{i},\text{\hspace{0.17em}}{x}_{n-1}\le {t}_{i}\le {x}_{n},\text{\hspace{0.17em}}0\le i\le {m}_{n}$ (81)

${P}_{n+1}\left({t}_{i}\right)={P}_{n+1}^{i},\text{\hspace{0.17em}}{x}_{n+1}\le {t}_{i}\le {x}_{n+2},\text{\hspace{0.17em}}0\le i\le {m}_{n+1}$ (82)

${P}_{n+2}\left({t}_{i}\right)={P}_{n+2}^{i},\text{\hspace{0.17em}}{x}_{n+2}\le {t}_{i}\le {x}_{n+3},\text{\hspace{0.17em}}0\le i\le {m}_{n+2}$ (83)

$⋮$

${P}_{n+k}\left({t}_{i}\right)={P}_{n+k}^{i},\text{\hspace{0.17em}}{x}_{n+k}\le {t}_{i}\le {x}_{n+k+1},\text{\hspace{0.17em}}0\le i\le {m}_{n+k}$ (84)

where ${R}_{i},i=1,\cdots ,n+k$ 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, $\left({t}_{i},{P}_{1}^{i}\right)$ , $0\le {t}_{i}\le {x}_{1}$ , $0\le i\le {m}_{1}$ ; $\left({t}_{i},{P}_{2}^{i}\right)$ , ${x}_{1}\le {t}_{i}\le {x}_{2}$ , $0\le i\le {m}_{2}$ ; $\cdots$ , $\left({t}_{i},{P}_{n}^{i}\right)$ , ${x}_{n-1}\le {t}_{i}\le {x}_{n}$ , $0\le i\le {m}_{n}$ are also given by the designer such that each of them determines a set of points through which the polynomials obtained by interpolation, ${P}_{i}\left(t\right),i=\stackrel{^}{1,n},$ passes respectively. Similarly, the $\left({t}_{i},{P}_{n+1}^{i}\right)$ , ${x}_{n+1}\le {t}_{i}\le {x}_{n+2}$ , $0\le i\le {m}_{n+1}$ ; $\left({t}_{i},{P}_{n+2}^{i}\right)$ , ${x}_{n+2}\le {t}_{i}\le {x}_{n+3}$ , $0\le i\le {m}_{n+2}$ ; $\cdots$ , $\left({t}_{i},{P}_{n+k}^{i}\right)$ , ${x}_{n+k}\le {t}_{i}\le {x}_{n+k+1}$ , $0\le i\le {m}_{n+k}$ ; are also given by the designer, such that each of them determines a set of points through which the polynomials obtained by interpolation, ${P}_{n+i}\left(t\right),i=\stackrel{^}{1,k},$ passes respectively. These polynomials are such that their corresponding polynomials by translation according to, ${P}_{n+i}\left(t-{h}_{i}\right),i=\stackrel{^}{1,k}$ fit with the remaining part of the bottle in the interval ${x}_{n+k}\le {t}_{i}\le {x}_{n+k+1},k\ge 1$ . The $r$ and ${x}_{i},\stackrel{^}{1,n},$ are constants given by the designer, while ${r}_{c}$ and ${h}_{c}$ are the radius and height of a cylinder of minimal area. We must find ${a}_{1}^{i},0\le i\le {m}_{1}$ ; ${a}_{2}^{i},0\le i\le {m}_{2}$ ; $\cdots$ ; ${a}_{n}^{i},0\le i\le {m}_{n}$ ; ${a}_{n+1}^{i},0\le i\le {m}_{n+1}$ ; ${a}_{n+2}^{i},0\le i\le {m}_{n+2}$ ; $\cdots$ ; ${a}_{n+k}^{i},0\le i\le {m}_{n+k}$ ; ${h}_{i},1\le i\le k$ ; ${x}_{n+i},1\le i\le k+1$ such that the area of the bottle of volume ${V}_{c}$ 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 ${a}_{1}^{i},0\le i\le {m}_{1}$ ; ${a}_{2}^{i},0\le i\le {m}_{2}$ ; $\cdots$ ; ${a}_{n}^{i},0\le i\le {m}_{n}$ ; ${a}_{n+1}^{i},0\le i\le {m}_{n+1}$ ; ${a}_{n+2}^{i},0\le i\le {m}_{n+2}$ ; $\cdots$ ; ${a}_{n+k}^{i},0\le i\le {m}_{n+k}$ , 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 ${x}_{n+1},{x}_{n+2}$ and ${h}_{1}$ can be found by solving the system of Equations (85), (86) and (87) given below, in using the Newton’s Method.

$\begin{array}{l}{\sum }_{i=1}^{n}\text{ }{I}_{i}+\text{π}{\left[{P}_{0}\left({x}_{n}\right)\right]}^{2}\left({h}_{c}-{x}_{n}\right)+\text{π}{\left({P}_{0}\left({x}_{n}\right)\right)}^{2}\left({x}_{n+1}-{h}_{c}\right)\\ +\text{\hspace{0.17em}}\text{π}{\int }_{{x}_{n+1}}^{{x}_{n+2}}{\left[{P}_{n+1}\left(x-{h}_{1}\right)\right]}^{2}\text{d}x={p}_{1}\left(\text{π}{r}_{c}^{2}{h}_{c}\right),\text{\hspace{0.17em}}0<{p}_{1}<1\end{array}$ (85)

where,

${I}_{1},{I}_{2},\cdots ,{I}_{n}$ are integrals that correspond to the polynomials

${P}_{1}\left(x\right),{P}_{2}\left(x\right),\cdots ,{P}_{n}\left(x\right)$ , and can be calculated as well as $\text{π}{\left[{P}_{0}\left({x}_{n}\right)\right]}^{2}\left({h}_{c}-{x}_{n}\right)$ .

${P}_{n+1}\left({x}_{n+1}-{h}_{1}\right)={R}_{n}={R}_{0}={R}_{n+1}$ (86)

${P}_{n+1}\left({x}_{n+2}-{h}_{1}\right)={R}_{n+2}$ (87)

Step 3 ${x}_{n+3}$ and1 can be found by solving the system of equations:

${P}_{n+2}\left({x}_{n+2}-{h}_{2}\right)={R}_{n+2}$

$\begin{array}{l}{\sum }_{i=1}^{n}\text{ }\text{ }{I}_{i}+\text{π}{\left[{P}_{0}\left({x}_{n}\right)\right]}^{2}\left({h}_{c}-{x}_{n}\right)+\text{π}{\left({P}_{0}\left({x}_{n}\right)\right)}^{2}\left({x}_{n+1}-{h}_{c}\right)+{I}_{n+1}\\ +\text{\hspace{0.17em}}\text{π}{\int }_{{x}_{n+2}}^{{x}_{n+3}}{\left[{P}_{n+2}\left(x-{h}_{2}\right)\right]}^{2}\text{d}x={p}_{2}\left(\text{π}{r}_{c}^{2}{h}_{c}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}0<{p}_{1}<{p}_{2}<1\end{array}$

${I}_{n+1}$ can be calculated in using the data calculated in step 2.

Step 4 ${x}_{n+4}$ and ${h}_{3}$ can be found by solving the system of equations:

${P}_{n+3}\left({x}_{n+3}-{h}_{3}\right)={R}_{n+3}$

$\begin{array}{l}{\sum }_{i=1}^{n}\text{ }{I}_{i}+\text{π}{\left[{P}_{0}\left({x}_{n}\right)\right]}^{2}\left({h}_{c}-{x}_{n}\right)+\text{π}{\left({P}_{0}\left({x}_{n}\right)\right)}^{2}\left({x}_{n+1}-{h}_{c}\right)+{I}_{n+1}+{I}_{n+2}\\ +\text{\hspace{0.17em}}\text{π}{\int }_{{x}_{n+3}}^{{x}_{n+4}}{\left[{P}_{n+3}\left(x-{h}_{3}\right)\right]}^{2}\text{d}x={p}_{3}\left(\text{π}{r}_{c}^{2}{h}_{c}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}0<{p}_{1}<{p}_{2}<{p}_{3}<1\end{array}$

${I}_{n+2}$ can be calculated in using the data calculated in step 3.

Step k ${x}_{n+k}$ and ${h}_{k-1}$ can be found by solving the system of equations:

${P}_{n+k-1}\left({x}_{n+k-1}-{h}_{k-1}\right)={R}_{n+k-1}$

$\begin{array}{l}{\sum }_{i=1}^{n}\text{ }{I}_{i}+\text{π}{\left[{P}_{0}\left({x}_{n}\right)\right]}^{2}\left({h}_{c}-{x}_{n}\right)+\text{π}{\left({P}_{0}\left({x}_{n}\right)\right)}^{2}\left({x}_{n+1}-{h}_{c}\right)+{\sum }_{i=1}^{k-2}\text{ }{I}_{n+i}\\ +\text{\hspace{0.17em}}\text{π}{\int }_{{x}_{n+k-1}}^{{x}_{n+k}}{\left[{P}_{n+k-1}\left(x-{h}_{k-1}\right)\right]}^{2}\text{d}x={p}_{k-1}\left(\text{π}{r}_{c}^{2}{h}_{c}\right),\text{\hspace{0.17em}}\text{\hspace{0.17em}}0<{p}_{1}<{p}_{2}<{p}_{3}<\cdots <{p}_{k-1}<1\end{array}$

${I}_{n+k-2}$ can be calculated in using the data calculated in step $k-1$ .

Step k+1 ${x}_{n+k+1}$ and ${h}_{k}$ can be found by solving the system of equations:

${P}_{n+k}\left({x}_{n+k}-{h}_{k}\right)={R}_{n+k}$

$\begin{array}{l}{\sum }_{i=1}^{n}\text{ }{I}_{i}+\text{π}{\left[{P}_{0}\left({x}_{n}\right)\right]}^{2}\left({h}_{c}-{x}_{n}\right)+\text{π}{\left({P}_{0}\left({x}_{n}\right)\right)}^{2}\left({x}_{n+1}-{h}_{c}\right)+{\sum }_{i=1}^{k-1}{I}_{n+i}\\ +\text{\hspace{0.17em}}\text{π}{\int }_{{x}_{n+k}}^{{x}_{n+k+1}}{\left[{P}_{n+k}\left(x-{h}_{k}\right)\right]}^{2}\text{d}x\\ ={p}_{k}\left(\text{π}{r}_{c}^{2}{h}_{c}\right),\text{\hspace{0.17em}}0<{p}_{1}<{p}_{2}<{p}_{3}<\cdots <{p}_{k}=1\end{array}$

${I}_{n+k-1}$ can 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

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.

nomials should be free of oscillations, we must use a computer program (optimizer of five sections).

2.4. 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 $\left({S}_{i},i=\stackrel{^}{1,5}\right)$ 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

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. Real bottle, (a) Silhouette of the real bottle according to the data of the Table 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.

computer program, the optimum area of the bottle according to the Criterion of Optimization by Compensation must be 365.96 cm2 when the base of the bottle is flat.

It is worth to emphasize that, the volume of 551.63 cm3 in Output Data, of Figure 6(a), for the bottle exceeds 537.5 cm3 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 such cavities, it is found a total volume of 14 cm3 approximately that have to be subtracted from 551.63 cm3. So we have 537.63 cm3 which is a good approximation of the volume of the real bottle. Similarly, the area 435.972 cm2 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 cm2, 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 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 ${x}_{4}$ 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 ${P}_{4}\left(x\right)$ . This number must be determined by moving section S5 to the right or left according to the sign of ${h}_{1}$ , positive or negative in ${P}_{4}\left(x-{h}_{1}\right)$ as it is shown in Figure 8. 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. Approximation by triangles and circles in order to determine the area of the petaloid base of the bottle.

Figure 8. Silhouette of the redesigned bottlein using a set of pointsof Section 1 (S1), Section 2 (S2)… Section 5 (S5) shown in Table 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.

according to the development of the model in Section 2.3, in general the polynomials ${P}_{n+i}\left(x\right),i=\stackrel{^}{1,k}$ are susceptible of changing, depending if some of ${h}_{1},{h}_{2},\cdots ,{h}_{k}$ are zero or different of zero. In this particular case we have only ${P}_{4}\left(x-{h}_{1}\right)$ from ${P}_{4}\left(x\right)$ and we will find that ${h}_{1}\approx 0$ , so that the Lagrange polynomials 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 ${P}_{4}\left(x-{h}_{1}\right)$ from ${P}_{4}\left(x\right)$ that fit correctly with the design of the silhouette of the bottle. To do this we put in Table 2, ${x}_{4}$ instead of 9.71 and ${x}_{5}$ instead of 12.01 and a new Table with new values for ${x}_{4}$ and ${x}_{5}$ must be found.

Equations (73) and (74) with ${R}_{4}=4.4$ , ${R}_{5}=3.8$ where $n=3$ , and in using the polynomial ${P}_{4}\left(x\right)$ become,

${P}_{4}\left({x}_{5}-{h}_{1}\right)-3.8=0$ , from Equation (74),

$-0.0406912{\left({x}_{5}-{h}_{1}\right)}^{3}+0.99403{\left({x}_{5}-{h}_{1}\right)}^{2}-7.40008\left({x}_{5}-{h}_{1}\right)+15.9863=0$ (88)

${P}_{4}\left({x}_{4}-{h}_{1}\right)-4.4=0$ , from Equation (73),

$-0.0406912{\left({x}_{4}-{h}_{1}\right)}^{3}+0.99403{\left({x}_{4}-{h}_{1}\right)}^{2}-7.40008\left({x}_{4}-{h}_{1}\right)+15.3863=0$ (89)

By looking Figure 8 we see that ${P}_{1}\left(x\right)=1.1$ , ${P}_{4}\left(x\right)=4.4$ , ${x}_{1}=1.5$ and $r=1.1$ , ${h}_{c}=8.81$ , ${r}_{c}=4.41$ . In using the data of Table 3, Equation (67) can be written as,

$\begin{array}{c}531.798009333735=\text{π}{\int }_{1.5}^{5.5}{\left[{P}_{2}\left(x\right)\right]}^{2}\text{d}x+\text{π}{\int }_{5.5}^{8.3}{\left[{P}_{3}\left(x\right)\right]}^{2}\text{d}x\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\text{π}{\left[{P}_{3}\left(8.3\right)\right]}^{2}\left({x}_{4}-8.3\right)+\text{π}{\int }_{{x}_{4}}^{{x}_{5}}{\left[{P}_{4}\left(x-{h}_{1}\right)\right]}^{2}\text{d}x\end{array}$

where ${P}_{3}\left(8.3\right)={P}_{4}\left({x}_{4}\right)=4.4$ .

By numerical integration of the first two integrals and taking into account that ${P}_{3}\left(8.3\right)=4.4$ , we have,

Table 3. Polynomials according to the points given in Table 2 corresponding to Section 1 (S1), 2 (S2), 3 (S3), 4 (S4) and 5 (S5) of the redesigned bottle.

$225.332507809904=\text{π}\left({x}_{4}-8.3\right){\left(4.4\right)}^{2}+\text{π}{\int }_{{x}_{4}}^{{x}_{5}}{\left[{P}_{4}\left(x-{h}_{1}\right)\right]}^{2}\text{d}x$

From which

$730.148748129941=60.8212337734984{x}_{4}+\text{π}{\int }_{{x}_{4}}^{{x}_{5}}{\left[{P}_{4}\left(x-{h}_{1}\right)\right]}^{2}\text{d}x$ (90)

Since

$\begin{array}{l}\text{π}{\int }_{{x}_{4}}^{{x}_{5}}{\left[{P}_{4}\left(x-{h}_{1}\right)\right]}^{2}\text{d}x\\ =\text{π}\left(0.00165577375744\frac{{\left({x}_{5}-{h}_{1}\right)}^{7}}{7}-0.080896547072\frac{{\left({x}_{5}-{h}_{1}\right)}^{6}}{6}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}+1.590331911492\frac{{\left({x}_{5}-{h}_{1}\right)}^{5}}{5}-16.32205962592\frac{{\left({x}_{5}-{h}_{1}\right)}^{4}}{4}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}+94.0975355844\frac{{\left({x}_{5}-{h}_{1}\right)}^{3}}{3}-292.840405808\frac{{\left({x}_{5}-{h}_{1}\right)}^{2}}{2}+391.49766769{x}_{5}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}-\text{π}\left(0.00165577375744\frac{{\left({x}_{4}-{h}_{1}\right)}^{7}}{7}-0.080896547072\frac{{\left({x}_{4}-{h}_{1}\right)}^{6}}{6}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}+1.590331911492\frac{{\left({x}_{4}-{h}_{1}\right)}^{5}}{5}-16.32205962592\frac{{\left({x}_{4}-{h}_{1}\right)}^{4}}{4}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}+94.0975355844\frac{{\left({x}_{4}-{h}_{1}\right)}^{3}}{3}-292.840405808\frac{{\left({x}_{4}-{h}_{1}\right)}^{2}}{2}+391.49766769{x}_{4}\right)\end{array}$

where integration was performed analytically, Equation (90) becomes finally,

$\begin{array}{l}60.8212337734984{x}_{4}+\text{π}\left(0.00165577375744\frac{{\left({x}_{5}-{h}_{1}\right)}^{7}}{7}-0.080896547072\frac{{\left({x}_{5}-{h}_{1}\right)}^{6}}{6}\\ +1.590331911492\frac{{\left({x}_{5}-{h}_{1}\right)}^{5}}{5}-16.32205962592\frac{{\left({x}_{5}-{h}_{1}\right)}^{4}}{4}\\ +94.0975355844\frac{{\left({x}_{5}-{h}_{1}\right)}^{3}}{3}-292.840405808\frac{{\left({x}_{5}-{h}_{1}\right)}^{2}}{2}+391.49766769{x}_{5}\right)\\ -\text{π}\left(0.00165577375744\frac{{\left({x}_{4}-{h}_{1}\right)}^{7}}{7}-0.080896547072\frac{{\left({x}_{4}-{h}_{1}\right)}^{6}}{6}\\ +1.590331911492\frac{{\left({x}_{4}-{h}_{1}\right)}^{5}}{5}-16.32205962592\frac{{\left({x}_{4}-{h}_{1}\right)}^{4}}{4}\\ +94.0975355844\frac{{\left({x}_{4}-{h}_{1}\right)}^{3}}{3}-292.840405808\frac{{\left({x}_{4}-{h}_{1}\right)}^{2}}{2}+391.49766769{x}_{4}\right)\\ -730.148748129941=0\end{array}$ (91)

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 ${f}_{1},{f}_{2}$ and ${f}_{3}$ , where:

${f}_{1}=-0.0406912{\left({x}_{5}-{h}_{1}\right)}^{3}+0.99403{\left({x}_{5}-{h}_{1}\right)}^{2}-7.40008\left({x}_{5}-{h}_{1}\right)+15.9863$

${f}_{2}=-0.0406912{\left({x}_{4}-{h}_{1}\right)}^{3}+0.99403{\left({x}_{4}-{h}_{1}\right)}^{2}-7.40008\left({x}_{4}-{h}_{1}\right)+15.3863$

$\begin{array}{c}{f}_{3}=60.8212337734984{x}_{4}+\text{π}\left(0.00165577375744\frac{{\left({x}_{5}-{h}_{1}\right)}^{7}}{7}-0.080896547072\frac{{\left({x}_{5}-{h}_{1}\right)}^{6}}{6}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+1.590331911492\frac{{\left({x}_{5}-{h}_{1}\right)}^{5}}{5}-16.32205962592\frac{{\left({x}_{5}-{h}_{1}\right)}^{4}}{4}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+94.0975355844\frac{{\left({x}_{5}-{h}_{1}\right)}^{3}}{3}-292.840405808\frac{{\left({x}_{5}-{h}_{1}\right)}^{2}}{2}+391.49766769{x}_{5}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\text{π}\left(0.00165577375744\frac{{\left({x}_{4}-{h}_{1}\right)}^{7}}{7}-0.080896547072\frac{{\left({x}_{4}-{h}_{1}\right)}^{6}}{6}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+1.590331911492\frac{{\left({x}_{4}-{h}_{1}\right)}^{5}}{5}-16.32205962592\frac{{\left({x}_{4}-{h}_{1}\right)}^{4}}{4}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+94.0975355844\frac{{\left({x}_{4}-{h}_{1}\right)}^{3}}{3}-292.840405808\frac{{\left({x}_{4}-{h}_{1}\right)}^{2}}{2}+391.49766769{x}_{4}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-730.148748129941\end{array}$

Additionally, we need the elements of the Jacobian matrix, which are:

${f}_{1x}=-0.1220736{\left({x}_{5}-{h}_{1}\right)}^{2}+1.98806\left({x}_{5}-{h}_{1}\right)-7.40008$

${f}_{1y}=0$

${f}_{1z}=0.1220736{\left({x}_{5}-{h}_{1}\right)}^{2}-1.98806\left({x}_{5}-{h}_{1}\right)+7.40008$

${f}_{2x}=0$

${f}_{2y}=-0.1220736{\left({x}_{4}-{h}_{1}\right)}^{2}+1.98806\left({x}_{4}-{h}_{1}\right)-7.40008$

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 ${f}_{1}$ , ${f}_{2}$ and ${f}_{3}$ (see Enter F). This figure also show the Enter Jacobian Matrix ( ${f}_{1x}$ , ${f}_{2x}$ , ${f}_{3x}$ , ${f}_{1y}$ , ${f}_{2y}$ , ${f}_{3y}$ , ${f}_{1z}$ , ${f}_{2z}$ and ${f}_{3z}$ ); Enter initial points ( ${x}_{o}$ , ${y}_{o}$ and ${z}_{o}$ ) as well as the Approximations to roots.

${f}_{2z}=0.1220736{\left({x}_{4}-h\right)}^{2}-1.98806\left({x}_{4}-h\right)+7.40008$

$\begin{array}{c}{f}_{3x}=\text{π}\left(0.00165577375744{\left({x}_{5}-{h}_{1}\right)}^{6}-0.080896547072{\left({x}_{5}-{h}_{1}\right)}^{5}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+1.590331911492{\left({x}_{5}-{h}_{1}\right)}^{4}-16.32205962592{\left({x}_{5}-{h}_{1}\right)}^{3}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+94.0975355844{\left({x}_{5}-{h}_{1}\right)}^{2}-292.840405808\left({x}_{5}-{h}_{1}\right)+391.49766769\right)\end{array}$

$\begin{array}{c}{f}_{3y}=-\text{π}\left(0.00165577375744{\left({x}_{4}-{h}_{1}\right)}^{6}-0.080896547072{\left({x}_{4}-{h}_{1}\right)}^{5}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+1.590331911492{\left({x}_{4}-{h}_{1}\right)}^{4}-16.32205962592{\left({x}_{4}-{h}_{1}\right)}^{3}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-94.0975355844{\left({x}_{4}-{h}_{1}\right)}^{2}-292.840405808\left({x}_{4}-{h}_{1}\right)+391.49766769\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+60.8212337734984\end{array}$

$\begin{array}{c}{f}_{3z}=\text{π}\left(-0.00165577375744{\left({x}_{5}-{h}_{1}\right)}^{6}+0.080896547072{\left({x}_{5}-{h}_{1}\right)}^{5}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-1.590331911492{\left({x}_{5}-{h}_{1}\right)}^{4}+16.32205962592{\left({x}_{5}-{h}_{1}\right)}^{3}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-94.0975355844{\left({x}_{5}-{h}_{1}\right)}^{2}+292.840405808\left({x}_{5}-{h}_{1}\right)\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\text{π}\left(-0.00165577375744{\left({x}_{4}-{h}_{1}\right)}^{6}+0.080896547072{\left({x}_{4}-{h}_{1}\right)}^{5}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-1.590331911492{\left({x}_{4}-{h}_{1}\right)}^{4}+16.32205962592{\left({x}_{4}-{h}_{1}\right)}^{3}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-94.0975355844{\left({x}_{4}-{h}_{1}\right)}^{2}+292.840405808\left({x}_{4}-{h}_{1}\right)\right)\end{array}$

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 $x={x}_{5},y={x}_{4}$ and $z={h}_{1}$ , we have

${f}_{1}=\frac{92501\ast z}{12500}-\frac{92501\ast x}{12500}{+}_{99403 *}\frac{{\left(x-z\right)}^{2}}{100000}-\frac{3179\ast {\left(x-z\right)}^{3}}{78125}+\frac{22498717100947801}{1407374883553280}$

${f}_{2}=\frac{92501\ast z}{12500}-\frac{92501\ast y}{12500}{+}_{99403 *}\frac{{\left(y-z\right)}^{2}}{100000}-\frac{3179\ast {\left(y-z\right)}^{3}}{78125}+\frac{21654292170815833}{1407374883553280}$

$\begin{array}{c}{f}_{3}=\frac{484\ast \text{π}y}{25}-\frac{31017640839698776609179259142809\ast \text{π}\ast \left(x-z\right)}{79228162514264337593543950336}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\frac{\left(31017640839698776609179259142809\ast \text{π}\ast \left(y-z\right)\right)}{79228162514264337593543950336}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\frac{515170290029661592553\ast \text{π}\ast {\left(x-z\right)}^{2}}{3518437208883200000}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{1379484460268183051814177\ast \text{π}\ast {\left(x-z\right)}^{3}}{43980465111040000000000}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\frac{\left(22432867934940286126331\ast \text{π}\ast {\left(x-z\right)}^{4}\right)}{5497558138880000000000}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{\left(397582977873\ast \text{π}\ast {\left(x-z\right)}^{5}\right)}{1250000000000}-\frac{\left(316002137\ast \text{π}\ast {\left(x-z\right)}^{6}\right)}{23437500000}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{\left(10106041\ast \text{π}\ast {\left(x-z\right)}^{7}\right)}{42724609375}+\frac{\left(515170290029661592553\ast \text{π}\ast {\left(y-z\right)}^{2}\right)}{3518437208883200000}\end{array}$

$\begin{array}{c}\text{\hspace{0.17em}}\text{\hspace{0.17em}}-\frac{\left(1379484460268183051814177\ast \text{π}\ast {\left(y-z\right)}^{3}\right)}{43980465111040000000000}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{\left(22432867934940286126331\ast \text{π}\ast {\left(y-z\right)}^{4}\right)}{5497558138880000000000}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\frac{\left(397582977873\ast \text{π}\ast {\left(y-z\right)}^{5}\right)}{1250000000000}+\frac{\left(316002137\ast \text{π}\ast {\left(y-z\right)}^{6}\right)}{23437500000}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\frac{\left(10106041\ast \text{π}\ast {\left(y-z\right)}^{7}\right)}{42724609375}-\frac{6422456308599681}{8796093022208}\end{array}$

${f}_{1x}=\frac{99403\ast x}{50000}-\frac{99403\ast z}{50000}{-}_{9537 *}\frac{{\left(x-z\right)}^{2}}{78125}-\frac{92501}{12500}$

${f}_{1y}=0$

${f}_{1z}=\frac{99403\ast z}{50000}-\frac{99403\ast x}{50000}{-}_{9537 *}\frac{{\left(x-z\right)}^{2}}{78125}-\frac{92501}{12500}$

${f}_{2x}=0$

${f}_{2y}=\frac{99403\ast y}{50000}-\frac{99403\ast z}{50000}{-}_{9537 *}\frac{{\left(y-z\right)}^{2}}{78125}-\frac{92501}{12500}$

${f}_{2z}=\frac{99403\ast z}{50000}-\frac{99403\ast y}{50000}{-}_{9537 *}\frac{{\left(y-z\right)}^{2}}{78125}-\frac{92501}{12500}$

$\begin{array}{c}{f}_{3x}=\frac{31017640839698776609179259142809\ast \text{π}}{79228162514264337593543950336}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\frac{515170290029661592553\ast \text{π}\ast \left(2x-2z\right)}{3518437208883200000}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{4138453380804549155442531\ast \text{π}\ast {\left(x-z\right)}^{2}}{43980465111040000000000}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\frac{22432867934940286126331\ast \text{π}\ast {\left(x-z\right)}^{3}}{1374389534720000000000}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{397582977873\ast \text{π}\ast {\left(x-z\right)}^{4}}{250000000000}-\frac{316002137\ast \text{π}\ast {\left(x-z\right)}^{5}}{3906250000}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{10106041\ast \text{π}\ast {\left(x-z\right)}^{6}}{6103515625}\end{array}$

$\begin{array}{c}{f}_{3y}=\frac{515170290029661592553\ast \text{π}\ast \left(2y-2z\right)}{3518437208883200000}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\frac{737094590335565475834206206607601\ast \text{π}}{1980704062856608439838598758400}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\frac{4138453380804549155442531\ast \text{π}\ast {\left(y-z\right)}^{2}}{43980465111040000000000}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{22432867934940286126331\ast \text{π}\ast {\left(y-z\right)}^{3}}{1374389534720000000000}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\frac{397582977873\ast \text{π}\ast {\left(y-z\right)}^{4}}{250000000000}+\frac{316002137\ast \text{π}\ast {\left(y-z\right)}^{5}}{3906250000}\end{array}$

$\text{\hspace{0.17em}}\text{\hspace{0.17em}}-\frac{10106041\ast \text{π}\ast {\left(y-z\right)}^{6}}{6103515625}$

$\begin{array}{c}{f}_{3z}=\frac{515170290029661592553\ast \text{π}\ast \left(2x-2z\right)}{3518437208883200000}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\frac{4138453380804549155442531\ast \text{π}\ast {\left(x-z\right)}^{2}}{43980465111040000000000}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{22432867934940286126331\ast \text{π}\ast {\left(x-z\right)}^{3}}{1374389534720000000000}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\frac{397582977873\ast \text{π}\ast {\left(x-z\right)}^{4}}{250000000000}+\frac{316002137\ast \text{π}\ast {\left(x-z\right)}^{5}}{3906250000}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\frac{10106041\ast \text{π}\ast {\left(x-z\right)}^{6}}{6103515625}-\frac{515170290029661592553\ast \text{π}\ast \left(2y-2z\right)}{3518437208883200000}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{4138453380804549155442531\ast \text{π}\ast {\left(y-z\right)}^{2}}{43980465111040000000000}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}-\frac{22432867934940286126331\ast \text{π}\ast {\left(y-z\right)}^{3}}{1374389534720000000000}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{397582977873\ast \text{π}\ast {\left(y-z\right)}^{4}}{250000000000}-\frac{316002137\ast \text{π}\ast {\left(y-z\right)}^{5}}{3906250000}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\frac{10106041\ast \text{π}\ast {\left(y-z\right)}^{6}}{6103515625}\end{array}$

Figure 9 shows the SNLEs, where it is possible to see the entry functions ${f}_{1}$ , ${f}_{2}$ and ${f}_{3}$ ; Jacobian matrix ${f}_{1x}$ , ${f}_{2x}$ , ${f}_{3x}$ , ${f}_{1y}$ , ${f}_{2y}$ , ${f}_{3y}$ , ${f}_{1z}$ , ${f}_{2z}$ and ${f}_{3z}$ ; initial points ${x}_{o}$ , ${y}_{o}$ and ${z}_{o}$ as well as the approximations to roots.

In Approximations to roots of the Solver for System of Nonlinear Equations, we can see that ${x}_{5}=12.0077\approx 12.01$ , ${x}_{4}=9.70766\approx 9.71$ and ${h}_{1}=-0.00231716\approx -0.0023$ . The values of ${x}_{5}=12.0077\approx 12.01$ , ${x}_{4}=9.70766\approx 9.71$ , ${x}_{3}=8.3$ , ${x}_{2}=5.5$ and ${x}_{1}=1.5$ , in Figure 9, are similar to values of Table 2; and ${P}_{4}\left(x-{h}_{1}\right)={P}_{4}\left(x+0.0023\right)\approx {P}_{4}\left(x\right)$ , ${P}_{3}\left(x\right)$ , $\cdots$ , ${P}_{1}\left(x\right)$ similar to values of 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 Figure 10, and which are shown in Table 4.

Therefore, with the data of Table 2 where ${x}_{5}=12.007\approx 12.01$ and ${x}_{4}=9.70766\approx 9.71$ 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 cm2 less than 365.96 cm2 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

Figure 10. Optimized redesigned bottle with flat base, (a) Silhouette of the optimized redesigned bottle according to the data of Table 2. 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 Table 4). 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.

Table 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, where ${x}_{5}=12.0077\approx 12.01$ and ${x}_{4}=9.70766\approx 9.7$ .

for the real bottle (see Figure 6), the area of the bottle is modified from 365.96 cm2 to 366.17 m2.

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.

3. 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] that establishes the inequality ${A}_{\text{sphere}}<{A}_{\text{cylinder}}<{A}_{\text{cube}}$ enclosing these areas the same volume. If we build a bottle from a cylinder of minimal area according to the Criterion of Optimization by Compensation we obtain the inequality ${A}_{\text{sphere}}<{A}_{\text{bottle}}<{A}_{\text{cyliner}}$ . 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 cm3 with a superficial area of 435.23 cm2 approximately. However, in using the Criterion of Optimization by Compensation for almost the same volume (537.63 cm3), the area can be reduced to 366.17 cm2 (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-

lation $\frac{\text{areaofthebottle}-\text{areaoftheoptimizedbottle}}{\text{areaoftheoptimizedbottle}}×100$ ) respecting to what

is established by the Criterion of Optimization by Compensation and, it has 69.06 cm2 more superficial area or it uses 3.995 g more of PET plastic. In indications such as * and ** are specified that the results 366.17 cm2 and 537.63 cm3 are

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

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 cm2 (This value is obtained by using

the relation $\frac{\text{areaofthebottle}-\text{areaoftheoptimizedbottle}}{\text{areaoftheoptimizedbottle}}×100=4.91$ , where

the area of the optimized bottle is 365.96 cm2). 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.916 cm2, while the area of the bottle in Optimum results is

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

Figure 12. Manufactured bottle, (a) Silhouette of the manufactured bottle according to the data of Table 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 cm2 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

$\frac{\text{areaofthemoldofthebottle}-\text{areaoftheoptimizedbottle}}{\text{areaoftheoptimizedbottle}}×100$ , where the area

of the optimized bottle is 365.96 cm2 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.

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

thickness of the bottles was changed, and the weight of the bottle reduced in average~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.

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

Acknowledgements

The authors would like to thank Universidad Nacional del Santa of Peru for its financial support in the development of the present research and the Brazilian Agency CNPq for its partial financial support. Process No 460733/2014?1.

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. Open Journal of Optimization, 6, 115-147. https://doi.org/10.4236/ojop.2017.63009

References

1. 1. Demeril, B. and Dave, F. (2009) Optimization of Poly (Ethylene Terephthalate) Bottles via Numerical Modeling: A Statistical Design of Experiment Approach. Journal of Applied Polymer Science, 114, 1126-1132. https://doi.org/10.1002/app.30644

2. 2. Pierre, D.A. (1987) Optimization Theory with Applications. Dover Publications, New York.

3. 3. Masood, S.H. and Keshava Murthy, V. (2005) Development of Collapsible PET Water Fountain Bottles. Journal of Materials Processing Technology, 162-163, 83-89. https://doi.org/10.1016/j.jmatprotec.2005.02.176

4. 4. Qingchung, H., Wenjian, S., Yanhui, L. and Yongsheng, W. (2012) Structural Optimization and Light-weight of PET Bottle Based on Abaqus. Advanced Materials Research, 346, 558-563.

5. 5. Mohammad, K.H. and Abunawas, A. (2010) Suitable Mathematical Model of PET for FEA Drop-Test Analysis. International Journal of Computer Science and Network Security, 10, No. 5.

6. 6. Silva de Miranda, C.A., Drummond Camera, J.J., Monken, O.P. and Gouvea, S. (2011) Design Optimization and Weight Reduction of 500 mL CSD PET Bottle through FEM Simulations. Journal of Materials Science and Engineering B., 1, 947-959.

7. 7. Han, J., Yamazaki, K. and Nishiyama, S. (2004) Optimization of the Crushing Characteristicsof Triangulated Aluminum Beverage Cans. Structural and Multidisciplinary Optimization, 28, 47-54. https://doi.org/10.1007/s00158-004-0418-8

8. 8. Han, J., Itoh, R., Nishiyama, S. and Yamasaki, K. (2005) Application of Structure Optimization Technique to Aluminum Beverage Bottle Design. Structural and Multidisciplinary Optimization, 29, 304-311. https://doi.org/10.1007/s00158-004-0485-x

9. 9. Han, J., Yamasaki, K., Itoh, R. and Nishiyama, S. (2006) Multiobjective Optimization of a Two-Piece Aluminum Beverage Bottle Considering Tactile Sensation of Heat and Embossing Formability. Structural and Multidisciplinary Optimization, 32, 141-151. https://doi.org/10.1007/s00158-005-0563-8

10. 10. Karl, H., WELS, H. and Martin, M. (2009) Simulation of the Filling of Polyethylene Terephthalate Bottles (PET) with a Volumetric Swirl Chamber Valve (Vodm 40335) on the Basis of Calculations Models and Experiments. Seventh International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia. http://www.cfd.com.au/cfd_conf09/PDFs/067HAI.pdf

11. 11. Hopmann, C.H., Rasche, S. and Windeck, C. (2015) Simulative Design and Process Optimization of the Two-Stage Stretch-Blow Molding Process. AIP Conference Proceedings, 1664, Article ID: 050011. https://doi.org/10.1063/1.4918415

12. 12. Marechal, C., Moreau, P. and Lochegnies, D. (2004) Numerical Optimization of a New Robotized Glass Blowing Process. Engineering with Computers, 19, 233-240.https://doi.org/10.1007/s00366-003-0263-1

13. 13. Thibault, F., Malo, A., Lanctot, B. and Diraddo, R. (2007) Preform Shape and Operating Condition Optimization for the Stretch Blow Molding Process. Polymer Engineering & Science, 47, 289-301. https://doi.org/10.1002/pen.20707

14. 14. Nocedal, J. and Wright, S.J. (1999) Numerical Optimization. Springer, New York.https://doi.org/10.1007/b98874

15. 15. McEvoy, J.P., Armstrong, C.G. and Crawford, R.J. (1998) Simulation of the Stretch Blow Molding Process of Pet Bottles. Advances in Polymer Technology, 17, 339-352.https://doi.org/10.1002/(SICI)1098-2329(199824)17:4<339::AID-ADV5>3.0.CO;2-S

16. 16. Haessly, W.P. and Ryan, M.P. (1993) Experimental Study and Finite Element Analysis of the Injection Blow Molding Process. Polymer Engineering & Science, 33, 1279-1287. https://doi.org/10.1002/pen.760331908

17. 17. Ha, S.-H., Choi, K.K. and Cho, S. (2010) Numerical Method for Shape Optimization Using T-Spline Based Isogeometric Method. Structural and Multidisciplinary Optimization, 42, 417-428. https://doi.org/10.1007/s00158-010-0503-0

18. 18. Reyna, L. and Moore, T. (2009) Optimization Criterion by Compensation for Designing Bottles with Surface of Revolution Used by Factories for Commercializing Foods and Drinks. Revista ECIPERU, 6, No. 1.

19. 19. Reyna, L. and Morales, H. (2010) Isometric Deformation of the Surface of Minimal Area of a Cylinder Increasing Its Volume. Revista ECIPERU, 7, No. 1.

20. 20. Hung, A.L., Minh, T.D., Hoai, T.L. and Duy, B.P. (2014) Structural Optimization of PET Bottle by Numerical Simulation. The 5th TSME International Conference on Mechanical Engineering, Chiang Mai, 2014.