## Recent Advances in Solution Methods for Nonlinear Evolution Equations, Fluid Flow, and Heat and Mass Transfer

View this Special IssueResearch Article | Open Access

# A Theorem for Finding Maximum Temperature in Wet Grinding

**Academic Editor:**Bogdan Dumitrescu

#### Abstract

We consider the solutions found in the literature for heat transfer in surface grinding, assuming a constant heat transfer coefficient for the coolant acting on the workpiece surface and a constant or linear heat flux profiles entering into the workpiece. From the integral form of the time-dependent temperature field reached in the workpiece, assuming the previous conditions, we prove that the maximum temperature always occurs in the stationary regime on the workpiece surface within the contact zone between the wheel and the workpiece. This result assures a very rapid method for the theoretical computation of the maximum temperature.

#### 1. Introduction

Surface grinding is an industrial machining process aiming to remove excess material of a workpiece by means of an abrasive wheel which rotates at high speed over its surface (see Figure 1). Most of the energy used in the grinding process is converted into heat, and it is accumulated within the contact zone between the wheel and the workpiece [1]. The high temperatures reached can thermally damage the quality of the workpiece, burning the workpiece, causing metallurgical phase transformations, softening (tempering) the surface layer with possible rehardening, onsetting residual tensile stresses, and causing cracks [1, 2]. Therefore, the determination of the temperature field evolution inside the workpiece is of great industrial importance [3–7]. In order to avoid thermal damage, coolant is usually delivered to the porous grinding zone at a high velocity, so friction is reduced and cooling by convection occurs. Some recommendations for the optimization of the coolant usage are found in [8].

##### 1.1. Heat Transfer Model in Surface Grinding

In surface grinding, the heat transfer inside the workpiece is usually modeled [9–11] by a strip heat source infinitely long and of width (m in SI units), which moves at a speed (m s^{−1}) over a semi-infinite solid surface (see Figure 1). Notice that, without losing generality, we can consider because reversing the motion of the workpiece relative to the grinding wheel, , is equivalent to change the direction of the -axis, .

Setting the Cartesian coordinate system fixed to wheel, as shown in Figure 1, the temperature field of the workpiece must satisfy the convective heat equation (notice that if represents the velocity of the motion of the workpiece, then ) [12, ]:where is the thermal diffusivity (m^{2} s^{−1}). Since initially the workpiece is at room temperature (K), (2) is subjected to the following initial condition:

From a theoretical point of view, heat transfer in wet surface grinding has been modelled assuming a constant heat transfer coefficient (W m^{−2} K^{−1}) on the workpiece surface [11, 13]. Despite the fact that this is not a very realistic assumption, because outside and inside the grinding zone are expected to have different values of [3], it is known that most of the energy removed by convection due to coolant occurs within the grinding zone [14]. Therefore, considering that the actual inside the grinding zone is equal to the one outside of it (constant assumption), it is expected that the temperature field does not differ very much from a nonconstant assumption. Moreover, the assumption of constant makes the problem to be analytically tractable. Assuming as well a dimensionless heat flux profile within the contact area between wheel and workpiece, we have the following boundary condition:where denotes the Heaviside function, (W m^{−1} K^{−1}) is the thermal conductivity, and (W m^{−2}) is the average heat flux entering into the workpiece along the contact width ; thus

##### 1.2. The Solutions

This paper will focus on the issue from a theoretical point of view by using the solution of (2)–(4) given by DesRuisseaux [10] for a constant heat flux profile and the solution given by Sauer [11] for a linear heat flux profile. In dry grinding, experiments suggest that the heat flux profile is linear (triangular with its apex at the inlet of the wheel-workpiece contact arc) in the case of upcutting [15]. Moreover, in wet grinding, according to [14], a constant heat flux profile is expected to occur in downcutting, while a linear heat flux profile is presumed to occur in upcutting. However, a linear heat flux profile may prevail in the literature, although a constant profile is used as well [3].

For both heat flux profiles, Green’s method is used to arrive to the desired solution. Assuming a Cartesian coordinate system fixed to the wheel (see Figure 1), the Green function is interpreted in this case as the temperature field rise with respect to the room temperature, evaluated at in a semi-infinite body () in which, at the position over the surface , an instantaneous point heat source of energy (J) appears at the instant , where the surface () is subject to a constant heat transfer coefficient : where the* radiation coefficient* (m^{−1}) is defined as

A formal derivation of (6) is given in [16].

The derivation of the solution of (2)–(4) by using (6) follows three steps:(i)Superposition in space to give temperature due to instantaneous line source, which acts on the surface parallel to the -axis.(ii)Superposition of line sources to give temperature due to a instantaneously acting strip source on the surface.(iii)Superposition in time to give temperature due to continuously acting strip source.

Let us consider now on the following dimensionless variables, which, taking into account (1), satisfy

In the literature, is termed as the Peclet number and as the Biot number. For a constant heat flux profile satisfying (5),DesRuisseaux in [10] provides the following solution: where

Similarly, for a linear heat flux profileSauer in [11] provides the following solution:where

It is worth noting that an equivalent solution to (14) can be found in [16] within the framework of the Samara-Valencia model [17].

*Notation.* We will drop the superscripts or when we refer to the temperature field regardless of the heat flux profile entering into the workpiece, so that we will write simultaneously (10) and (14) as

##### 1.3. Maximum Temperature Search Approach

In order to avoid thermal damage, we need to control the maximum temperature in grinding, calculated from the temperature fields expressions given in (10) and (14). Therefore, let us state the following theorem, which will be very useful for the numerical computation of .

Theorem 1. *The maximum temperature of the time-dependent temperature field is**where*

The proof of Theorem 1 follows three steps:(1)The maximum temperature is reached in the stationary regime, that is to say, when .(2)The maximum temperature is reached on the workpiece surface, that is to say, at .(3)The maximum temperature is located within the grinding zone, that is to say, .

As far as author’s knowledge, the occurrence of the maximum temperature in the stationary regime is taken for granted in the literature and we did not find any formal proof of it. Moreover, we have found the following considerations about the maximum temperature:(i)Jaeger calculates in [9] the temperature field in the case of dry grinding, considering a constant heat flux profile. He takes for granted that the maximum temperature occurs in the stationary regime on the workpiece surface. Moreover, he provides some approximations for large in order to calculate the surface temperature in the stationary regime, from which he estimates the maximum temperature.(ii)Malkin and Guo in [3] follow Jaeger approximated formula for maximum temperature.(iii)DesRuisseaux in [13] presents the temperature field in wet grinding for a constant heat flux profile. He does not consider his solution in the transient regime. Also, he plots the temperature for and for different depths , in order to illustrate that the maximum temperature occurs on the surface and within the grinding zone.(iv)Sauer in [11] calculates the temperature field in wet grinding for a linear heat flux profile, but again he considers his solution only in the stationary regime. He says that the region in which the maximum temperature is expected lies on Therefore, it might be expected that the maximum temperature occurs outside the grinding zone and at certain depth from the surface.

However, for the case of a constant heat flux profile, we find in the literature the following partial results regarding to Theorem 1:(i)The maximum temperature occurs when [18, Sect. 3.1].(ii)When , the temperature field reaches the stationary regime; that is, [16, Sect. 5.2].(iii)When the temperature field reaches the stationary regime, the maximum temperature occurs at the boundary, that is, on the workpiece surface [18, Sect. 3.2].(iv)In order to locate the maximum temperature in the stationary regime on the workpiece surface, it can be proved for small Biot numbers that always this maximum occurs within the grinding zone [18, Sect. 3.3].

It is worth noting also that the numerical simulation of the time-dependent field temperature by using FEM analysis [19] agrees with Theorem 1.

The principal scope of this paper is just to generalize the above results for any Biot number, for both a constant and a linear heat flux profiles. We will see later on that this is a challenging mathematical problem.

This paper is organized as follows. Section 2 is devoted to prove that the maximum temperature is found on the surface in the stationary regime, for both a linear and a constant heat flux profiles. Section 3 proves that the maximum temperature occurs within the grinding zone, regardless as well of the heat flux profile considered: constant or linear. Section 4 provides a numerical example of the temperature field in the stationary regime, considering a constant and a linear heat flux profile. Our conclusions are summarized in Section 5. In the Appendix we collect some auxiliary lemmas in order to support the results given in body of the paper.

#### 2. Maximum Temperature on the Surface in the Stationary Regime

In this section we will prove first that is a monotonically increasing function of (Theorem 4), so that must be asymptotically reached when . This result prevents the occurrence of a temperature peak during the transient regime. Once this is proved, we will prove that when , the temperature field does not evolve over time; that is, the stationary regime is asymptotically reached when (Proposition 8). This fact is the key result to prove that the temperature field in the stationary regime can be interpreted as a harmonic function, so that, according to the maximum principle of harmonics functions, the maximum temperature is found to be on the boundary, that is, on the workpiece surface (Theorem 9).

Proposition 2. *The function satisfies *

*Proof. *According to [20, Sect. 2.2], we haveSincewe may rewrite (21) asMoreover, the following inequality is also satisfied: Taking and in (24), we haveso, according to (11), it follows (20), as we wanted to prove.

Proposition 3. *The functions satisfies*

*Proof. *According to the definition of given in (12), we may rewrite this function aswhich is positive since , , and the integrand is positive within the integration interval. Similarly, let us rewrite , defined in (15), aswhich is also positive because , , and the integrand is positive within the integration interval.

Theorem 4. *The time-dependent temperature field is a monotonically increasing function of :*

*Proof. *Performing the derivative with respect to in (16) and taking into account (20) and (26), we have

Corollary 5. *The maximum temperature is asymptotically reached when .*

Let us prove now that the stationary regime is asymptotically reached at . For this purpose, let us start proving the following propositions.

Proposition 6. *The function satisfies*

*Proof. *According to the definition of the given in (11), we haveTaking into account now the asymptotic expansion [21, Eq. 7.12.1] up to first order, we haveso, substituting (34) in (32), we conclude (31), as we wanted to prove.

Proposition 7. *The functions satisfy*

*Proof. *According to the integral representation (27), we have since the integral given in (36) is finite for all . Similarly, from (28), we have

Proposition 8. *The stationary regime is asymptotically reached, when :*

*Proof. *See that, from (16), we haveand, according to (31) and (35), we conclude (38).

The above proposition allows us to prove now the following result.

Theorem 9. *The maximum temperature is reached at the surface of the workpiece, .*

*Proof. *Let us denote the Cartesian coordinates fixed to the workpiece as (see Figure 2), where, at , both coordinates systems are overlapped. Consider as well the temperature field referred to as , whereApplying chain’s rule, it is easy to see that the convective term in (2) disappears:Taking now similar dimensionless variables as in (8), that is to say,(2) becomesAccording to (40) and (42), it is clear that the relationship between the temperature field referred to the workpiece and the one referred to the wheel is given bySince (38) is satisfied for all , according to (44), we obtainThus, performing the limit in (43), taking into account (45), we haveso is a harmonic function. Therefore, according to the maximum principle of harmonics functions [22, Chap. VI], the maximum temperature of (47) occurs on the boundary, that is, on the workpiece surface, , as we wanted to prove.

#### 3. Maximum Temperature within the Grinding Zone

According to Corollary 5 and Theorem 9, the maximum temperature must be found on the surface of the workpiece in the stationary regime, . Therefore, let us define the function as the surface temperature in the stationary regime:

The idea is to prove that is a monotonically increasing (Theorem 14) or decreasing (Theorem 10) function outside the grinding zone. Therefore, must be found within the grinding zone, , as it is stated in Theorem 1.

Theorem 10. *The surface temperature in the stationary regime is a monotonically decreasing function on the right-hand side of the leading edge:*

*Proof. *From (48), we havewhere, according to (12), and, according to (15), Let us set now the following variables:thussince the integration interval in (50) means that . Notice also that implies that ; so taking into account (54), we haveOn the one hand, rewriting now (51) in terms of and , we havedue to (55) and the fact that the function is a monotonically decreasing function for positive arguments [23, Sect. 1.3. Def. 1].

On the other hand, (52) in terms of and reads out asApplying the mean value theorem for derivatives [24, Theorem 4.5], we haveTherefore, , so and (58) is written asTaking into account (59) in (57), we concludeFinally, collecting the results (56) and (60) and taking into account (20), we may writeso integrating (61) over the interval , according to [25, Theorem 3.3.4], then (49) follows directly, as we wanted to prove.

In the above theorem we have just seen that the sign of for is determined by the sign of . This is not the case for ; so before arriving to Theorem 14, we have to prove the following three propositions.

Proposition 11. *Let **Then is a positive and monotonically decreasing function for positive arguments:*

*Proof. *Particularizing (20) to , we haveso it directly follows (63), since .

In order to prove (64), rewrite as follows, taking into account (11): Define now the functionso Considering that [20, Sect. 2.2]we can say Since according to (70), it follows (64), as we wanted to prove.

Proposition 12. *Let **If **then consider the following.*(1)*There is an unique crossing point such that*(2)*The following inequalities are satisfied:*

*Proof. *Notice that so the values or satisfy (74). Nevertheless, these values do not belong to the dominion of . Therefore, let us substitute (72) in (74), that is to say,soThe equation with the “+” sign leads to but this contradicts (73). The equation with the “−” sign leads towhich is positive, according to (73), so it falls within the dominion of .

Notice that (79) is an actual crossing point between the graphs of and as functions of sinceIndeedso, if thenbut this contradicts (73). Therefore, there is a unique crossing point satisfying (74), as we wanted to prove.

In order to prove (75), see that the function for given satisfies the assumptions of Lemma A.1 of the appendix:According to (72), , for all .According to (76), .Also, such that . Indeed, it is easy to prove that the only value for satisfying is given by

Then, according to (A.4), is the location of the absolute maximum of for given . Moreover, as it is shown in Figure 3, the crossing point satisfies

Indeed, taking into account (79) and (85) in (86), we arrive atwhich is trivially satisfied and is coherent with (73).

Notice as well that substituting (85) in (72), we obtain

Now, since , by definition of absolute maximum, we haveand taking into account (88), we obtain

See now that the functions and , where , satisfy the assumptions of Lemma A.2 of the Appendix, that is to say,(1) such that , according to (74);(2), according to (80);(3)taking , according to (86) , and (90) reads as .

Therefore, the conclusion of Lemma A.2 given in (A.9)-(A.10) is equivalent to (75), as we wanted to prove.

Proposition 13. *The following improper integral is independent of , if :*

*Proof. *Using the definition (72) of the function, we have Applying the following integral [26, Eq. 3.325]we arrive atso it follows directly (91), as we wanted to prove.

Theorem 14. *The surface temperature in the stationary regime is a monotonically increasing function on the left-hand side of the trailing edge:*

*Proof. *Substituting in (58) the values of and given in (53) and taking into account (54), we havewhereTaking into account that , rewrite asso (97) is satisfied and (96) reads out asSubstituting (99) in (52), we haveTaking into account (51) and (100) and recalling the definitions of the and functions, (72) and (62), respectively, let us express (50) as follows:where or and depending on whether we are in the constant or linear case, respectively. Therefore, according to (98) and (102), we haveNotice that the integral given in (101) satisfies the assumptions of Lemma A.3 given in the Appendix:(1)According to Proposition 13 and taking into account (103), we have(2)According to Proposition 11, we have (3)According to Proposition 12 and taking into account (103), we have

Therefore, since , we conclude (95), as we wanted to prove.

#### 4. Numerical Example

For the numerical evaluation of the temperature field, we will use the thermal properties of a VT20 titanium alloy workpiece [27]:

Also, we will take the following values for the grinding regime:

Figures 4 and 5 show the temperature field in the stationary regime for a constant and a linear heat flux profile, respectively. Since the temperature gradient nearby the grinding zone is very high, note that the -axis is zoomed with respect to the -axis in order to obtain a more detailed view of the contour plots. Notice as well that, although the temperature fields are quite different in both graphs, the maximum temperature is found to be on the surface within the grinding zone. For a constant heat flux profile, the value and location on the surface of the maximum temperature isand for a linear heat flux profile is

Note that the maximum temperature is higher in the linear case because the heat flux is skewed to the leading edge, but in the constant case the heat flux is uniform along the grinding zone. Nonetheless, the temperature field is more skewed to the trailing edge in the constant case because of the motion of the workpiece, being compensated the effect of the motion of the workpiece with the skewness of the heat flux profile in the linear case.

It is worth noting that the principal result of this paper (Theorem 1) assures a very rapid method for searching the maximum temperature. Indeed, the numerical search of the maximum in a three-dimensional field within a unbounded region is reduced to the numerical computation of the maximum in a one-dimensional function within a particular interval; that is, .

#### 5. Conclusions

We have considered the solutions found in the literature to model the heat transfer in surface grinding, assuming a constant and a linear heat flux profile entering into the workpiece and a constant heat transfer coefficient for the coolant acting on the workpiece surface. From these time-dependent temperature fields in dimensionless variables , we have analyzed the location of the maximum temperature in the workpiece, proving that the latter must be located on the workpiece surface, in the stationary regime and within the grinding zone, that is, Theorem 1. As long as the heat transfer model stated in Section 1.1 can be applied to the real grinding process, Theorem 1 has the following consequences:(i)The proof that the maximum temperature is found in the stationary regime prevents to encounter a temperature peak during the transient regime.(ii)Also, the proof about the location of the maximum temperature on the workpiece surface prevents to find an internal thermal damage in the workpiece when its surface seems not to be damaged.(iii)Moreover, the proof that the maximum temperature is located in the grinding zone on the surface and in the stationary regime converts the search of the maximum temperature in a three-dimensional temperature field within an unbounded region into the search of the maximum for a one-dimensional function within a particular interval, which is a much more simple problem from a computational point of view.

These theoretical results are intended to be useful in order to avoid thermal damage in surface wet grinding.

#### Appendix

#### Auxiliary Lemmas

Lemma A.1. *Let be differentiable in and continuous in that satisfies **If **then is an absolute maximum: *

*Proof. *Since the function is continuous in , according to the Extreme-Value Theorem of Continuous Functions [24, Theorem 3.12], the absolute maximum belongs to the interval . Suppose that this absolute maximum is not located in , say in . Thus, by definition of absolute maximum and according to (A.2), we haveIf or , then, according to (A.1), , but this contradicts (A.5). If , then, since is differentiable in , is also a local maximum, so , but this contradicts (A.3). Therefore, the absolute maximum must be located in , as we wanted to prove in (A.4).

Lemma A.2. *Let be two continuous and differentiable functions which have a single crossing point, that is to say, **If**then*

*Proof. *Let us consider the following differentiable function:According to (A.6), we haveSince is a differentiable function, according to (A.7), we haveso, according to (A.12),Therefore, changes sign in the neighborhood of , that is to say, On the one hand, assume now thatHowever, by hypothesis (A.8),Therefore, from (A.16) and (A.17), within the interval a sign change occurs; so applying Bolzano’s theorem, it follows that but this contradicts the hypothesis that the crossing point is unique (A.6). Thus such that , that is, as we wanted to prove in (A.9).

On the other hand, notice that from (A.15) and (A.19) we haveAssume now thatTherefore, from (A.20) and (A.21), within the interval a sign change occurs; so applying Bolzano’s theorem, it follows thatbut this contradicts that the crossing point is unique (A.6). Thus such that , that is, as we wanted to prove in (A.10).

Lemma A.3. *Let be three real functions that satisfy the following.*(1)*The following improper integrals of and are convergent and satisfy *(2)*The function is positive and monotonically decreasing for positive arguments *(3)*There is an unique crossing point for and such that *

Proposition A.4. *Then*

*Proof. *From (A.24), we havethat is,Since is a positive (A.25) and decreasing (A.26) function, then, due to (A.27), we have According to [25, Theorem 3.3.4] we may integrate in (A.33) preserving the sense of the inequality:Similarly, we havethuswhere we have taken into account (A.31) and (A.34). Rewrite now (A.36) aswhich is just (A.29), as we wanted to prove.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### References

- S. Malkin,
*Grinding Technology: Theory and Application of Machining with Abrasives*, Ellis Horwood, John Wiley & Sons, 1989. - J. B. J. W. Hegeman, J. T. M. de Hosson, and G. de With, “Grinding of WC-Co hardmetals,”
*Wear*, vol. 248, no. 1-2, pp. 187–196, 2001. View at: Publisher Site | Google Scholar - S. Malkin and C. Guo, “Thermal analysis of grinding,”
*CIRP Annals—Manufacturing Technology*, vol. 56, no. 2, pp. 760–782, 2007. View at: Publisher Site | Google Scholar - C. Guo and S. Malkin, “Analysis of energy partition in grinding,”
*Journal of Engineering for Industry*, vol. 117, no. 1, pp. 55–61, 1995. View at: Publisher Site | Google Scholar - S. Malkin and R. B. Anderson, “Thermal aspects of grinding: part 1—energy partition,”
*Journal of Engineering for Industry*, vol. 96, no. 4, pp. 1177–1183, 1974. View at: Publisher Site | Google Scholar - A. S. Lavine and B. F. von Turkovich, “Thermal aspects of grinding: the effects of heat generation at the shear planes,”
*CIRP Annals—Manufacturing Technology*, vol. 40, no. 1, pp. 343–345, 1991. View at: Publisher Site | Google Scholar - A. S. Lavine, “An exact solution for surface temperature in down grinding,”
*International Journal of Heat and Mass Transfer*, vol. 43, no. 24, pp. 4447–4456, 2000. View at: Publisher Site | Google Scholar | Zentralblatt MATH - M. N. Morgan, A. R. Jackson, H. Wu, V. Baines-Jones, A. Batako, and W. B. Rowe, “Optimisation of fluid application in grinding,”
*CIRP Annals—Manufacturing Technology*, vol. 57, no. 1, pp. 363–366, 2008. View at: Publisher Site | Google Scholar - J. C. Jaeger, “Moving sources of heat and the temperature at sliding contracts,”
*Proceedings of the Royal Society of New South Wales*, vol. 76, pp. 204–224, 1942. View at: Google Scholar - N. R. DesRuisseaux,
*Thermal aspects of grinding processes [Ph.D. dissertation]*, University of Cincinnati, 1968. - W. J. Sauer,
*Thermal aspects of grinding [Ph.D. thesis]*, Carnegie-Mellon University, 1971. - H. S. Carslaw and J. C. Jaeger,
*Conduction of Heat in Solids*, Oxford Science Publications, Oxford, UK, 1947. View at: MathSciNet - N. R. DesRuisseaux and R. D. Zerkle, “Temperature in semi-infinite and cylindrical bodies subjected to moving heat surfaces and surface Cooling,”
*Journal of Heat Transfer*, vol. 92, pp. 456–464, 1970. View at: Google Scholar - S. Ohishi and Y. Furukawa, “Analysis of workpiece temperature and grinding burn in creep feed grinding,”
*Bulletin of the JSME*, vol. 28, no. 242, pp. 1775–1781, 1985. View at: Publisher Site | Google Scholar - I. Zarudi and L. C. Zhang, “A revisit to some wheel—workpiece interaction problems in surface grinding,”
*International Journal of Machine Tools and Manufacture*, vol. 42, no. 8, pp. 905–913, 2002. View at: Publisher Site | Google Scholar - J. L. González-Santander, J. M. Isidro, and G. Martín, “An analysis of the transient regime temperature field in wet grinding,”
*Journal of Engineering Mathematics*, vol. 90, pp. 141–171, 2015. View at: Publisher Site | Google Scholar | MathSciNet - D. L. Skuratov, Y. L. Ratis, I. A. Selezneva, J. Pérez, P. Fernández de Córdoba, and J. F. Urchueguía, “Mathematical modelling and analytical solution for workpiece temperature in grinding,”
*Applied Mathematical Modelling*, vol. 31, no. 6, pp. 1039–1047, 2007. View at: Publisher Site | Google Scholar | Zentralblatt MATH - J. L. González-Santander and G. Martín, “Closed form expression for the surface temperature in wet grinding: application to maximum temperature evaluation,”
*Journal of Engineering Mathematics*, vol. 90, pp. 173–193, 2015. View at: Publisher Site | Google Scholar | MathSciNet - M. Mahdi and L. Zhang, “The finite element thermal analysis of grinding processes by ADINA,”
*Computers and Structures*, vol. 56, no. 2-3, pp. 313–320, 1995. View at: Publisher Site | Google Scholar - N. N. Lebedev,
*Special Functions and their applications*, Dover Publications, New York, NY, USA, 1972. View at: MathSciNet - F. W. J. Olver, D. W. Lozier, R. F. Boisvert, and C. W. Clark,
*NIST Handbook of Mathematical Functions*, Cambridge University Press, New York, NY, USA, 2010. - J. B. Conway,
*Functions of One Complex Variable*, vol. 11 of*Graduate Texts in Mathematics*, Springer, Berlin, Germany, 2nd edition, 1978. View at: MathSciNet - D. S. Mitrinović,
*Analytic Inequalities*, Springer, New York, NY, USA, 1970. View at: MathSciNet - T. M. Apostol,
*Calculus*, vol. 1, John Wiley & Sons, New York, NY, USA, 1967. - W. F. Trench,
*Introduction to Real Analysis*, Books and Monographs, Book 7, 2013, http://digitalcommons.trinity.edu/mono/7. - I. S. Gradsthteyn and M. Ryzhik,
*Table of Integrals, Series and Products*, Academic Press, New York, NY, USA, 7th edition, 2007. - V. I. Murav'ev, A. V. Yakimov, and A. V. Chernyshev, “Effect of deformation, welding, and electrocontact heating on the properties of titanium alloy VT20 in pressed and welded structures,”
*Metal Science and Heat Treatment*, vol. 45, no. 11-12, pp. 419–422, 2003. View at: Publisher Site | Google Scholar

#### Copyright

Copyright © 2015 J. L. González-Santander and G. Martín. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.