Mathematical Problems in Engineering
Volume 2015 (2015), Article ID 150493, 13 pages
http://dx.doi.org/10.1155/2015/150493

Research Article

A Theorem for Finding Maximum Temperature in Wet Grinding

J. L. González-Santander and G. Martín

Departamento de Ciencias Experimentales y Matemáticas, Universidad Católica de Valencia “San Vicente Mártir”, C/Guillem de Castro 96, 46001 Valencia, Spain

Received 10 February 2015; Revised 18 May 2015; Accepted 28 May 2015

Academic Editor: Bogdan Dumitrescu

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.

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

Figure 1: Setup in surface grinding.

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⁻¹) 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² s⁻¹). 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⁻² K⁻¹) 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⁻¹ K⁻¹) is the thermal conductivity, and (W m⁻²) 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⁻¹) 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 profile

Sauer 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 have

Since

we may rewrite (21) as

Moreover, the following inequality is also satisfied:

Taking and in (24), we have

so, 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 as

which is positive since , , and the integrand is positive within the integration interval. Similarly, let us rewrite , defined in (15), as

which 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 have

Taking into account now the asymptotic expansion [21, Eq. 7.12.1]

up to first order, we have

so, 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 have

and, 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 , where

Applying 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) becomes

According 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 by

Since (38) is satisfied for all , according to (44), we obtain

Thus, performing the limit in (43), taking into account (45), we have

so

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.

Figure 2: Relationship between the coordinates systems, one fixed to the wheel and the other one to the workpiece.

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 have

where, according to (12),

and, according to (15),

Let us set now the following variables:

thus

since the integration interval in (50) means that . Notice also that implies that ; so taking into account (54), we have

On the one hand, rewriting now (51) in terms of and , we have

due 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 as

Applying the mean value theorem for derivatives [24, Theorem 4.5], we have

Therefore, , so and (58) is written as

Taking into account (59) in (57), we conclude

Finally, collecting the results (56) and (60) and taking into account (20), we may write

so 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 have

so it directly follows (63), since .
In order to prove (64), rewrite as follows, taking into account (11):

Define now the function

so

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,

so

The equation with the “+” sign leads to but this contradicts (73). The equation with the “−” sign leads to

which 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 since

Indeed

so, if

then

but 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

Figure 3: The crossing point lies between the and .

Indeed, taking into account (79) and (85) in (86), we arrive at

which 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 have

and 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 at

so 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 have

where

Taking into account that , rewrite as

so (97) is satisfied and (96) reads out as

Substituting (99) in (52), we have

Taking 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 have

Notice 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 is

and for a linear heat flux profile is

Figure 4: Temperature field for a constant heat flux profile in the stationary regime.

Figure 5: Temperature field for a linear heat flux profile in the stationary regime.

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 have

If 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 have

Since is a differentiable function, according to (A.7), we have

so, according to (A.12),

Therefore, changes sign in the neighborhood of , that is to say,

On the one hand, assume now that

However, 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 have

Assume now that

Therefore, from (A.20) and (A.21), within the interval a sign change occurs; so applying Bolzano’s theorem, it follows that

but 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 have

that 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 have

thus

where we have taken into account (A.31) and (A.34). Rewrite now (A.36) as

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

S. Malkin and C. Guo, “Thermal analysis of grinding,” CIRP Annals—Manufacturing Technology, vol. 56, no. 2, pp. 760–782, 2007.

C. Guo and S. Malkin, “Analysis of energy partition in grinding,” Journal of Engineering for Industry, vol. 117, no. 1, pp. 55–61, 1995.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

N. N. Lebedev, Special Functions and their applications, Dover Publications, New York, NY, USA, 1972.

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.

D. S. Mitrinović, Analytic Inequalities, Springer, New York, NY, USA, 1970.

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.