Abstract
In this work we consider a problem in the context of the generalized theory of thermoelasticity for a half space. The material of the half space is functionally graded in which Lamé's modulii are functions of the vertical distance from the surface of the medium. The surface is traction free and subjected to a time dependent thermal shock. The problem was solved by using the Laplace transform method together with the perturbation technique. The obtained results are discussed and compared with the solution when Lamé's modulii are constants. Numerical results are computed and represented graphically for the temperature, displacement and stress distributions.
Keywords:
Half Space; Generalized Thermoelasticity; Perturbation method; FGM; Variable Lamé's modulii
1 INTRODUCTION
The theory of generalized thermoelasticity with one relaxation time was introduced by (Lord and Shulman, 1967Lord H., Shulman, Y., (1967). A Generalized Dynamical Theory Of Thermo- elasticity. J Mech Phys Solid 15: 299-309.). In this theory Cattaneo-Maxwell law of heat conduction replaces the conventional Fourier's law. The heat equation associated with this theory is a hyperbolic one and hence automatically eliminates the paradox of infinite speeds of propagation inherent in both the uncoupled and the coupled theories of thermoelasticity. For many problems involving steep heat gradients and when short time effects are sought this theory is indispensable. (Sherief and El-Maghraby, 2003Sherief, H., El-Maghraby, N., (2003). An internal pennyshaped crack in an infinite thermoelastic solid. J Thermal Stresses 26: 333-352, 2005Sherief, H., El-Maghraby, N., (2005). A mode-I crack problem for an infinite space in generalized thermoelasticity. J Thermal Stresses 28: 465-484) solved some crack problems for this theory. Sherief and Hamza has obtained the solution of axisymmetric problems in cylindrical regions in (1994Sherief, H., Hamza, F., (1994). Generalized thermoelastic problem of a thick plate under axisymmetric temperature distribution. J Thermal Stresses 17: 435-452) and in spherical regions in (1996Sherief, H., Hamza, F., (1996). Generalized two-dimensional thermoelastic problems in spherical regions under axisymmetric distributions. J Thermal Stresses 19: 55-76). (Sherief and Ezzat, 1994Sherief, H., Ezzat, M., (1994). Solution of the generalized problem of thermoelasticity in the form of series of functions. J Thermal Stresses 17: 75-95) have obtained the solution in the form of series. (Sherief and Dhaliwal, 1981Sherief, H., Dhaliwal, R., (1981). A generalized one-dimensional thermal shock problem for small times. J Thermal Stresses 4: 407-420.) used asymptotic expansions to obtain the solution of a 1D problem and to find the locations of the wave fronts and the speed of propagation of thermoelastic waves. This theory was extended to deal with thermoelastic diffusion in (2004Sherief, H., Hamza, F., Saleh, H., (2004). The theory of generalized thermoelastic diffusion. Int J Engng Sci. 42: 591-608.), micropolarity of the medium in (2005Sherief, H., Hamza, F., El-Sayed, A., (2005). Theory of generalized micropolar thermoelasticity and an axisymmetric half-space problem. J. Thermal Stresses 28: 409-437.), viscoelastic effects in (2011Sherief, H., Allam, M., El-Hagary, M., (2011). Generalized theory of thermoviscoelasticity and a half-space problem. Int J Thermophys 32:1271-1295). (Sherief et.al, 2010Sherief, H., El-sayed, A., Abd El-Latief, A. (2010) Fractional order theory of thermoelasticity. Int J Solids Structures 47: 269-275.) extended this theory using fractional derivatives. Under this theory, (Sherief and Abd-Ellatief, 2013Sherief, H., Abd El-Latief, A., (2013). Effect of variable thermal conductivity on a Half-Space under the Fractional order Theory of thermoelasticity. Int J Mech Sci 74: 185-189, 2014aSherief, H., Abd El-Latief, A., (2014a). Application of Fractional order Theory of thermoelasticity to a 1D Problem for a Half-Space. ZAMM 94: 509 - 515., 2014bSherief, H., Abd El-Latief, A., (2014b). Application of fractional order theory of thermoelasticity to a 2D problem for a half-space. Appl Math Computation 248: 584-592) have solved some problems.
Functionally graded materials (FGMs) are a new class of materials with the material properties varying continuously along specified directions. Due to their desirable properties, FGMs have been increasing used in modern engineering applications ((Wang, 2013Wang, H., (2013). An effective approach for transient thermal analysis in a functionally graded hollow cylinder. Int J Heat Mass Transfer 67: 499-505.)). For example, the functionally graded metal-ceramic composites have been used as thermal barriers or thermal shields in various applications ((Lee et.al, 1996Lee, W., Stinton, D., Berndt, C., Erdogan, F., Lee, Y., Mutasim, Z., (1996). Concept of functionally graded materials for advanced thermal barrier coating applications. J Am Ceram Soc. 79: 3003-3012.), (Tsukamoto, 2010Tsukamoto, H., (2010). Design of functionally graded thermal barrier coatings based on a nonlinear micromechanical approach. Comput Mater Sci. 50: 429-436.)). Especially, (Praveen et.al., 1999Praveen, G., Chin, C., Reddy, J., (1999). Thermoelastic analysis of functionally graded ceramic-metal cylinder. J Eng Mech. 125: 1259-1267.) in severe temperature environments, such as extremely high temperature and thermal shock, widely potential applications are opening for FGMs.
In recent years, (Asgari and Akhlaghi, 2009Asgari, M., Akhlaghi, M., (2009). Transient heat conduction in two-dimensional functionally graded hollow cylinder with finite length. Heat Mass Transfer 45: 1383 1392.) composition of several different materials is often used in structural components in order to optimize the thermal resistance and temperature distribution of structures subjected to thermal loading. Other works that consider FGM can be found in ((Nahvi et.al., 2008Nahvi, H., Salehipour, H., Shahidi, A., (2008). Generalized coupled thermoelasticity of functionally graded annular disk considering the Lord-Shulman theory. Composite Structures 83:168-179.), (Zafarmand et. al., 2015Zafarmand, H., Salehi, M., Asemi, K., (2015). Theree dimensional free vibtation transient analysis of two directional functionally graded thick cylinderical panels under impact loading. Lat Am J Solids Stru. 12(2):205-225), (Dey et.al., 2015Dey, S., Adhikari, S., Karmakar, A., (2015). Impact response of functionally graded conical shells. Lat Am J Solids Stru. 12(1):133-152.)).
2 FORMULATION OF THE PROBLEM
In this work, we consider a homogeneous isotropic thermoelastic solid occupying the region x ≥ 0 composed of a FGM material whose Lamé's parameters depend on the vertical distance '' x '' from the surface. The surface of the half-space is taken to be traction free and is subjected to a thermal shock that is a function of time. We assume also that there are no body forces or heat sources inside the medium. The initial conditions are assumed to be homogenous.
From the physics of the problem, it is clear that all the variables depend on x and t only. The displacement vector will thus have the form.
The equation of motion is given by
Where σ = σxx is the normal component of the stress tensor given by
where λ, μ are Lamé's modulii. αt is the coefficient of thermal expansion and θ = T - T0 where T is the absolute temperature, T 0 is the temperature of the medium in its normal state such that |θ| << 1.
Substituting from Equation (3) into Equation (2), we get
where ( is the density. The equation of the heat conduction has the form
where CE is the specific heat per unit mass in the absence of deformation, k is the thermal conductivity and (0 is the relaxation time. From now on, we shall take λ,μ in the form
where (0 , (0 and "a" are constants. Thus equations (3-5) take the form
We shall use the following non-dimensional variables
x' = cηx, u' = cηu, t' = c 2ηt, τ'0 = c 2ητ0, a' = a/cη
Using the above non-dimensional variables, equations (7)-(9) take the form
The remaining components of the stress tensor are given by:
where ε = T 0(3λ0 + 2μ0)20 + 2μ0)kη, β2 = /(λ
We assume that the boundary conditions have the form
where f(t) is a known function of t
3 SOLUTION IN THE LAPLACE TRANSFORM DOMAIN
Applying the Laplace transform with parameter s defined by the relation
to both sides of equations (10)-(13), we get the following equations
In order to solve the above equations, we shall use the perturbation method. We expand the temperature, displacement and stress functions as follows:
θ = θ(0) + aθ(1) + a 2θ(2) + ...
u = u(0) + au (1) + a 2 u (2) + ...
σ = σ(0) + aσ(1) + a 2σ(2) + ...
where θ( i ) and u ( i ) are functions to be determined, i = 0, 1, ...
Equations (18-19) gives, upon equating the coefficients of "a " in both sides up to order 1
Using expansion of equation (17), we obtain
Eliminating 21) and (23), we get between equations (
where D =
The general solution of equation (27) which is bounded at infinity can be written as
where A1 and A2 are parameters depending on s only and k 1 and k 2 are the roots with positive real parts of the characteristic equation
k4 - k 2[s 2 + (1 + ε)(1 + τ0 s)s] + s 3(1 + τ0 s) = 0
From equations (21) and (28), we get
The boundary conditions (15-16) give
In order to determine A 1, A 2 we shall use the boundary conditions (30a), (30c) to obtain
Eliminating 22) and Eq. (24), we get between Eq. (
Substituting from equations (31-32) into the right hand side of equation (33) and factorizing the left hand side as before, we obtain
where
The solution of the above non homogenous differential equation takes the form
where B1 and B2 are unknown parameters depending on s while the remaining parameters, due to the particular solution, are given by
By the same manner the displacement differential equation takes the form
where
The solution of Eq. (36) is given by
where L1 and L2 are unknown parameters depending on s while the remaining parameters, due to the particular solution, are given by
Substituting form Eqs. (28), (35) and (37) into Eq. (24), comparing the coefficient of the exponentials in the resulting equation, we get
Substituting from Eqs. (35), (37) into the boundary conditions (30b), (30d), we obtain
Solving the above system, we obtain
The other constants can be easily obtained from equation (38a, b,c).
Thus we have the solution in the Laplace transform domain for the temperature and displacement functions. The stress σ can be determined from Eqs. (25-26). The solution in the physical domain can be obtained by using numerical inversion method outlined in ((Honig and Hirdes, 1984Honig, G., Hirdes, U., (1984). A Method For The Numerical Inversion of the Laplace Transform. J Comp Appl Math 10:113-132.)). This method has been successfully used in solving many problems in the theory of thermoelasticity (see e.g. ((Sherief and Megahed, 1999Sherief, H., Megahed, F., (1999). A two-dimensional thermoelasticity problem for a half space subjected to heat sources, Int. j. solids struct. 36: 1369-1382)) , ((Sherief and Anwar, 1988Sherief, H., Anwar, M., (1988). A problem in generalized thermoelasticity for an infinitely long annular cylinder, Int. j. engng. sci 26: 301-306))
4 NUMERICAL RESULTS AND DISCUSSION
The copper material was chosen for purposes of numerical evaluations. The constants of the problem are shown in table 1
The computations were carried out for different functions f(t). We have chosen the following two cases:
Case 1
Case 2
The temperature, displacement and stress components were calculated by using the numerical method for the inversion of Laplace transform illustrated in (Honig and Hirdes, 1984Honig, G., Hirdes, U., (1984). A Method For The Numerical Inversion of the Laplace Transform. J Comp Appl Math 10:113-132.). The FORTRAN programming language was used on a personal computer. The accuracy maintained was 5 digits for the numerical program.
4.2 FIGURES and Discussion
Figures (1) and (2) show the temperature distribution for t = 0.1 for cases 1 and 2, respectively. In these figures and subsequent ones the dotted lines denote the case of constant Lamé's modulii (CL) while solid lines represent the case of variable Lamé's modulii (VL). We note that changing of Lamé's modulii has very small effect on the temperature profile.
Figures (3) and (4) show the displacement distribution for t = 0.1 for cases 1 and 2, respectively. We note that an increase in the values of Lamé's modulii results in an increase in the value of the displacement.
Figures (5) and (6) show the stress distribution for t = 0.1 for cases 1 and 2, respectively. We note that an increase in the values of Lamé's modulii results in a decrease in the value of the stress component σxx .
Figures (7-12) show the distributions of θ, u and σ for different times. In these figures the solid lines and the dotted lines represent t = 0.08 and t = 0.12, respectively. We note that the increase of the time results in an increase in the thermal wave penetration into the medium. The same effect is noticed for the displacement and the stress distributions.
Comparison between the values of all the functions are shown in tables 2-4.
As usual in dealing with problems of the theory of generalized thermoelasticity, the finite speed of wave propagation is apparent. For t = 0.1 the solution is identically zero for x > 0.6341, approximately, for both cases.
We note that for case 2, the temperature and stress distributions have two discontinuities. The first discontinuity is at x = 0.973 approximately while the second discontinuity is at x = 0.6341 approximately. The first discontinuity in the temperature is very small in value and does not show in the figure. The displacement distribution is continuous but has discontinuous first derivatives at the above locations. These locations are the locations of the wave fronts. These discontinuities are due to the fact that the input thermal shock f 2(t) is a discontinuous function.
For case 1, where the thermal shock is represented by a continuous function f 1(t), all the considered functions are continuous. The temperature and stress have discontinuous first derivatives at the wave fronts.
References
- Asgari, M., Akhlaghi, M., (2009). Transient heat conduction in two-dimensional functionally graded hollow cylinder with finite length. Heat Mass Transfer 45: 1383 1392.
- Dey, S., Adhikari, S., Karmakar, A., (2015). Impact response of functionally graded conical shells. Lat Am J Solids Stru. 12(1):133-152.
- Honig, G., Hirdes, U., (1984). A Method For The Numerical Inversion of the Laplace Transform. J Comp Appl Math 10:113-132.
- Lee, W., Stinton, D., Berndt, C., Erdogan, F., Lee, Y., Mutasim, Z., (1996). Concept of functionally graded materials for advanced thermal barrier coating applications. J Am Ceram Soc. 79: 3003-3012.
- Lord H., Shulman, Y., (1967). A Generalized Dynamical Theory Of Thermo- elasticity. J Mech Phys Solid 15: 299-309.
- Nahvi, H., Salehipour, H., Shahidi, A., (2008). Generalized coupled thermoelasticity of functionally graded annular disk considering the Lord-Shulman theory. Composite Structures 83:168-179.
- Praveen, G., Chin, C., Reddy, J., (1999). Thermoelastic analysis of functionally graded ceramic-metal cylinder. J Eng Mech. 125: 1259-1267.
- Sherief, H., Abd El-Latief, A., (2013). Effect of variable thermal conductivity on a Half-Space under the Fractional order Theory of thermoelasticity. Int J Mech Sci 74: 185-189
- Sherief, H., Abd El-Latief, A., (2014a). Application of Fractional order Theory of thermoelasticity to a 1D Problem for a Half-Space. ZAMM 94: 509 - 515.
- Sherief, H., Abd El-Latief, A., (2014b). Application of fractional order theory of thermoelasticity to a 2D problem for a half-space. Appl Math Computation 248: 584-592
- Sherief, H., Allam, M., El-Hagary, M., (2011). Generalized theory of thermoviscoelasticity and a half-space problem. Int J Thermophys 32:1271-1295
- Sherief, H., Anwar, M., (1988). A problem in generalized thermoelasticity for an infinitely long annular cylinder, Int. j. engng. sci 26: 301-306
- Sherief, H., Dhaliwal, R., (1981). A generalized one-dimensional thermal shock problem for small times. J Thermal Stresses 4: 407-420.
- Sherief, H., El-Maghraby, N., (2003). An internal pennyshaped crack in an infinite thermoelastic solid. J Thermal Stresses 26: 333-352
- Sherief, H., El-Maghraby, N., (2005). A mode-I crack problem for an infinite space in generalized thermoelasticity. J Thermal Stresses 28: 465-484
- Sherief, H., El-sayed, A., Abd El-Latief, A. (2010) Fractional order theory of thermoelasticity. Int J Solids Structures 47: 269-275.
- Sherief, H., Ezzat, M., (1994). Solution of the generalized problem of thermoelasticity in the form of series of functions. J Thermal Stresses 17: 75-95
- Sherief, H., Hamza, F., (1994). Generalized thermoelastic problem of a thick plate under axisymmetric temperature distribution. J Thermal Stresses 17: 435-452
- Sherief, H., Hamza, F., (1996). Generalized two-dimensional thermoelastic problems in spherical regions under axisymmetric distributions. J Thermal Stresses 19: 55-76
- Sherief, H., Hamza, F., El-Sayed, A., (2005). Theory of generalized micropolar thermoelasticity and an axisymmetric half-space problem. J. Thermal Stresses 28: 409-437.
- Sherief, H., Hamza, F., Saleh, H., (2004). The theory of generalized thermoelastic diffusion. Int J Engng Sci. 42: 591-608.
- Sherief, H., Megahed, F., (1999). A two-dimensional thermoelasticity problem for a half space subjected to heat sources, Int. j. solids struct. 36: 1369-1382
- Tsukamoto, H., (2010). Design of functionally graded thermal barrier coatings based on a nonlinear micromechanical approach. Comput Mater Sci. 50: 429-436.
- Wang, H., (2013). An effective approach for transient thermal analysis in a functionally graded hollow cylinder. Int J Heat Mass Transfer 67: 499-505.
- Zafarmand, H., Salehi, M., Asemi, K., (2015). Theree dimensional free vibtation transient analysis of two directional functionally graded thick cylinderical panels under impact loading. Lat Am J Solids Stru. 12(2):205-225
Publication Dates
-
Publication in this collection
Apr 2016
History
-
Received
26 Apr 2015 -
Reviewed
10 Dec 2015 -
Accepted
06 Jan 2016