ISRN Mechanical EngineeringVolume 2014 (2014), Article ID 694259, 6 pageshttp://dx.doi.org/10.1155/2014/694259

Research Article

Exact Solution of Thermoelastic Problem for a One-Dimensional Bar without Energy Dissipation

A. M. Abd El-Latief¹ and S. E. Khader²

¹Department of Mathematics, Faculty of Science, University of Alexandria, Egypt²Department of Mathematical and Theoretical Physics, Atomic Energy Authority, Egypt

Received 3 November 2013; Accepted 26 November 2013; Published 4 February 2014

Academic Editors: J. K. Chen and X. H. Liu

Copyright © 2014 A. M. Abd El-Latief and S. E. Khader. 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 a homogeneous isotropic thermoelastic half-space in the context of the theory of thermoelasticity without energy dissipation. There are no body forces or heat source acting on the half-space. The surface of the half-space is affected by a time dependent thermal shock and is traction free. The Laplace transform with respect to time is used. The inverse transforms are obtained in an exact manner for the temperature, thermal stress, and displacement distributions. These solutions are represented graphically and discussed for several cases of the applied heating. Comparison is made between the predictions here and those of the theory of thermoelasticity with one relaxation time.

1. Introduction

In 1967, Lord and Shulman [1] introduced the theory of generalized thermoelasticity with one relaxation time for an isotropic body. This theory corrects the unrealistic conclusions of the older theories (the uncoupled and the coupled theories of thermoelasticity) that heat waves travel with infinite speeds. The Heat conduction law of this theory is the Cattaneo law which is different from Fourier’s law utilized in both the coupled and the uncoupled theories. Among the contribution to the subject are the works in [2–6]. In 1972 Green and Lindsay [7] developed the theory of generalized thermoelasticity with two relaxation times, based on a generalized inequality of thermodynamics. In this theory both the equations of motion and of heat conduction are hyperbolic. The heat conduction law is the same as Fourier’s law when the system has a centre of symmetry. Among the contributions to this theory are the works in [8, 9].

Green and Naghdi [10–12] have formulated three new models of thermoelasticity. In one of these models Green and Naghdi [12] predict that the internal rate of production of entropy is identically zero; that is, there is no dissipation of thermal energy. This theory (GN theory) is known as thermoelasticity without energy dissipation theory. In the development of this theory the thermal displacement gradient is considered as a constitutive variable, whereas in the conventional development of a thermoelasticity theory, the temperature gradient is taken as a constitutive variable [12]. A couple of uniqueness theorems have been proved in [13, 14], and one-dimensional waves in a half-space and in an unbounded body have been studied in [15, 16]. In view of some experimental evidence available in favour of finiteness of heat propagation speed, generalized thermoelasticity theories are considered to be more realistic than the conventional theory in dealing with practical problems involving very large heat fluxes at short intervals, like those occurring in laser units and energy channels [17].

2. Formulation of the Problem

In this work, we consider a homogeneous isotropic and thermoelastic half-space in the context of the theory of thermoelasticity without energy dissipation. There are no body forces or heat source acting on the half-space taken as . The surface of the half-space is affected by a time dependent thermal shock and is traction free.

From the physics of the problem, all functions will depend on the distance and the time only. The problem is thus one-dimensional. The displacement components have the form

The basic equations due to Green and Nagdhi [12] in the absence of body forces and heat sources for isotropic elastic medium are given by the following:

(i) equation of motion:

(ii) the equation of heat conduction:

(iii) the constitutive equation: where and are Lamé’s modulii, is the absolute temperature of the medium, and is a material constant given by , where is the coefficient of linear thermal expansion and is a reference temperature assumed to be such that . is a material constant, characteristic of the theory and is the specific heat at constant strain.

Let us introduce the following nondimension variables:

The governing equations (2)–(4) in nondimensional form become (dropping the asterisks for convenience) where , and.

3. Solution of the Problem in the Laplace Transform Domain

Applying the Laplace transform with parameter defined by the relation [14] to both sides of (6), we obtain

Eliminating from (12) and (14), we get where .

Equation (11) can be factorized as where and are the roots of the characteristic equation and are given by where and .

The solutions of (12), bounded as tends to have the form where and are some parameters depending on only.

To calculate using (9), we obtain

Integrating both sides of (16) with respect to , we get

Substituting (15) and (17) into (10), we obtain

The boundary conditions of the problem can be written as follows:

Applying the Laplace transform with parameter to both sides of (19), we get

Substituting (15) and (18) into (20), we obtain the following system of linear equations in the unknown parameters and :

The solution of the above equations is where .

Substituting (22) into (15), (17), and (18), we obtain The above equations can be written as where This completes the solution of the problem in the Laplace transform domain.

4. Exact Solution in the Physical Domain

We shall derive the temperature and stress distributions in an exact manner; by applying the inverse Laplace transform to both sides of (24) and (26), we get the following expressions: In obtaining the above equations, we have used the following property of the Laplace Transform: where is the Heaviside unit step function; we have According to the definition of Heaviside unit step function, we have two cases:(i)if , then the integral right hand side of (31) becomes ;(ii)while if in this case the integral vanishes; thus (31) becomes

Now, using (32) and (25), we obtain the exact form of the displacement distribution in the form This completes the solution of the problem in the physical domain.

5. Numerical Results and Discussion

We will consider the solution of the problem for different choices for the function as shown in Table 1, where and are constants.

Table 1

Case

1

2

3

Substituting into (29), (30), and (33), we obtain the following.

Case 1. Consider the following:

Case 2. Consider the following:

Case 3. Consider the following:

Obviously the above exact solution verifies the boundary conditions introduced in (19) for all different cases.

We shall now compare, graphically, our results, where the theory of thermoelasticity without energy dissipation (GN theory) was used with the results of the solution of a similar problem in [18], where the generalized theory of thermoelasticity due to Lord-Shulman (LS theory) was used. Due to the nature of the characteristic roots of the LS theory, the authors of [18] could not obtain a closed form exact solution. They used asymptotic expansions valid for small values of time.

For purposes of numerical evaluations, the calculations were performed for a copper-like material for which [19]

The temperature, displacement, and stress distributions are plotted for one value of time, namely, for . In all figures, LS theory is represented by dotted lines, while GN theory is represented by solid lines.

For the three cases, the temperature distributions are shown in Figures 1, 4, and 7, the displacement distributions are shown in Figures 2, 5, and 8, while the stress distributions are shown in Figures 3, 6, and 9.

Figure 1: Temperature distribution for Case 1.

Figure 2: Displacement distribution for Case 1.

Figure 3: Stress distribution for Case 1.

Figure 4: Temperature distribution for Case 2.

Figure 5: Displacement distribution for Case 2.

Figure 6: Stress distribution for Case 2.

Figure 7: Temperature distribution for Case 3.

Figure 8: Displacement distribution for Case 3.

Figure 9: Stress distribution for Case 3.

At the locations and the temperature and stress distributions have two finite jumps, while the displacement is continuous everywhere having a discontinuous first derivative at the above locations. The first discontinuity is not shown in Figures 1, 4, and 7 because it is very small .

and are the locations of the wave fronts. We thus have two waves propagating into the medium from the surface of the half space. The first is mainly mechanical and has velocity equal to . The second is mainly thermal and has velocity equal to .

For the solution for all functions is identically zero. Thus, the solution is confined to a finite region of space and does not extend to infinity. This is not the case for the coupled and the uncoupled theories of thermoelasticity, where the solution extends to infinity instantaneously signifying an infinite speed of propagation of waves. The values of the velocities and locations of the wave fronts are as follows:

The difference between the predictions of the theories of LS and GN is most apparent in the graphs of the temperature distribution. In the LS theory the temperature decreases monotonically signifying continuous dissipation of heat energy. This is not the case for GN theory.

Conflict of Interests

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

References

H. Lord and Y. Shulman, “A generalized dynamical theory of thermo-elasticity,” Journal of the Mechanics and Physics of Solids, vol. 15, no. 5, pp. 299–309, 1967.

H. Sherief and N. M. El-Maghraby, “An internal penny-shaped crack in an infinite thermoelastic solid,” Journal of Thermal Stresses, vol. 26, no. 4, pp. 333–352, 2003.

H. Sherief and F. A. Hamza, “Generalized two-dimensional thermoelastic problems in spherical regions under axisymmetric distributions,” Journal of Thermal Stresses, vol. 19, no. 1, pp. 55–76, 1996.

H. Sherief and M. Anwar, “A problem in generalized thermoelasticity for an infinitely long annular cylinder composed of two different materials,” Acta Mechanica, vol. 80, no. 1-2, pp. 137–149, 1989.

M. Anwar and H. Sherief, “Boundary integral equation formulation of generalized thermoelasticity in a Laplace-transform domain,” Applied Mathematical Modelling, vol. 12, no. 2, pp. 161–166, 1988.

H. Sherief and M. A. Ezzat, “Solution of the generalized problem of thermoelasticity in the form of series of functions,” Journal of Thermal Stresses, vol. 17, no. 1, pp. 75–95, 1994.

A. E. Green and K. A. Lindsay, “Thermoelasticity,” Journal of Elasticity, vol. 2, no. 1, pp. 1–7, 1972.

H. Sherief, “State space approach to thermoelasticity with two relaxation times,” International Journal of Engineering Science, vol. 31, no. 8, pp. 1177–1189, 1993.

H. Sherief, “Problem in electromagneto thermoelasticity for an infinitely long solid conducting circular cylinder with thermal relaxation,” International Journal of Engineering Science, vol. 32, no. 7, pp. 1137–1149, 1994.

A. E. Green and P. M. Naghdi, “A re-examination of the basic postulate of thermomechanics,” Proceedings of the Royal Society A, vol. 432, no. 1885, pp. 171–194, 1991.

A. E. Green and P. M. Naghdi, “On undamped heat waves in an elastic solid,” Journal of Thermal Stresses, vol. 15, no. 2, pp. 253–264, 1992.

A. E. Green and P. M. Naghdi, “Thermoelasticity without energy dissipation,” Journal of Elasticity, vol. 31, no. 3, pp. 189–208, 1993.

D. S. Chandrasekharaiah, “A uniqueness theorem in the theory of thermoelasticity without energy dissipation,” Journal of Thermal Stresses, vol. 19, no. 3, pp. 267–272, 1996.

D. S. Chandrasekharaiah, “A note on the uniqueness of solution in the linear theory of thermoelasticity without energy dissipation,” Journal of Elasticity, vol. 43, no. 3, pp. 279–283, 1996.

D. S. Chandrasekharaiah, “One-dimensional wave propagation in the linear theory of thermoelasticity without energy dissipation,” Journal of Thermal Stresses, vol. 19, no. 8, pp. 695–710, 1996.

D. S. Chandrasekharaiah, “Thermoelastic plane waves without energy dissipation,” Mechanics Research Communications, vol. 23, no. 5, pp. 549–555, 1996.

D. S. Chandrasekharaiah, “Thermoelasticity with second sound: a review,” Applied Mechanics Reviews, vol. 39, no. 3, pp. 355–376, 1986.

R. Dhaliwal and H. Sherief, “A reciprocity theorem and integral representation for generalized thermoelasticity,” Journal of Mathematical and Physical Sciences, vol. 15, pp. 537–552, 1981.

H. Sherief and H. Saleh, “A problem for an infinite thermoelastic body with a spherical cavity,” International Journal of Engineering Science, vol. 36, no. 4, pp. 473–487, 1998.