Effect of Initial Stress on a Thermoelastic Functionally Graded Material with Energy Dissipation ()
1. Introduction
Thermoelasticity has many applications and one of them is thermoelastic stress analysis (TSA) which consists of qualities appraisal of a model of radar antennae support column subjected to the bending; the magneto elasticity has practical applications in geophysics, optics, plasma physics etc. The linear theory of elasticity is important for engineering when the structure or material is under external or internal loadings. When the structure or the material’s temperature changes this implies a thermal loading. A thermal shock is a sudden variation in temperature from hot to cold or vice versa, when the interval of thermal shock is equal to the order of the lowest natural frequency of structure, this type of problems classified under thermoelasticity which assumes that temperature is expressed from the 1st law of thermodynamic. Biot [1] presented the uncoupled thermo-elasticity theory. There was a shortcoming in this theory and the temperature resulting from the equation of heat follows a parabolic behavior which means that there is an infinite speed of propagation for it, which has a paradox with the physical observations.
Two relaxations times were firstly introduced by Muller in generalized theory of thermoelasticity. With the aid of his entropy inequality he inserted a close approach by considering restrictions on a class off constitutive equations. Green and Naghdi [2] introduced three completely different models of thermoelasticity, which are called thermoelasticity of type I, thermoelasticity of type II and thermoelasticity of type III. Green and Naghdi created a concept in generalized thermoelasticity. He postulated that there is no wasted energy in thermo-elastic process i.e. thermo-elasticity without energy dissipation. This principal means that the heat flow does not consist of energy dissipation. When Fourier conductivity is dominant, then the temperature relation has undammed the solutions of thermal wave without wasted energy or energy dissipation. Many researchers devoted to carry out different articles in thermoelasticity, in the light of the GN theories of type II or/and of type III. Chandrasekharaiah [3] showed uniqueness theorems using energy method. Ezzat et al. [4] [5] studied the usage of GN theories a mathematical model of magneto-thermoelasticity. Kumar et al. [6] and Sharma et al. [7] studied several cases in GN theories.
Functionally graded materials was introduced first as a thermal barrier by a group of material scientists [8], other than homogenous materials (FGM) is different class of advances composite material which were developed in recent decades and used in many engineering applications. The properties of (FGM) is gradually change in material properties with respect to spatial coordinates and has the ability to reduce thermal stresses when a thermal shock occur. For these properties (FGM) has a widely applications in industry and engineering applications as we mentioned such as aerospace, nuclear reactors, pressure vessels and pipes. Kumar et al. discussed the influence of rotation, magnetic on (FGM) subjected to mechanical load [9], Mohamed H. Hendy discussed the application of fractional order theory to a functionally graded perfect conducting thermoelastic half space with variable Lame’s Modulii [10], Vel and Batra, and Qian and Batra studied the three dimensional steady or transient thermal stress problems of functionally graded rectangular plate [11] “sadek”, Javaheri and Eslami studied thermal bulking of functionally graded plates [12], Sherief and Abd El-Latief discussed modeling of variable Lamé’s modulii for a FGM generalized thermoelastic half Space [13]. Abd-Alla et al. [14] discussed Propagation of Rayleigh waves in magneto-thermoelastic half-space of a homogenous orthotropic material under the effect of rotation, initial stress and gravity field. Abd-Alla et al. [15] studied Rayleigh waves in magneto elastic half-space of orthotropic material under influence of initial stress and gravity field. Abd-Alla et al. [16] studied Influence of rotation, magnetic field, initial stress, and gravity on Rayleigh waves in a homogeneous orthotropic elastic half-space. Othman et al. [17] studied the effect of initial stress and Hall current on a magneto-thermoelastic porous medium with microtemperatures. Kumar et al. [18] studied the effect of initial stress on the propagation characteristics of waves in fiber-reinforced transversely isotropic thermoelastic material under an inviscid liquid layer.
In our paper, we introduce the effect initial stress on a functionally graded thermoelastic solid. Normal mode technique is used to find the expressions for the variables considered. The comparison of different theories of thermo-elasticity, i.e. Green Naghdi type II, Green Lindsay, Chandrasekharaiah and Tzou (DPL) model is carried out. The distributions of variables are displayeded graphically, and to the best of our knowledge, this case not discussed before.
2. Basic Equations
The governing equation of an isotropic, homogenous elastic medium without heat source:
2.1. Equation of Motion
(1)
2.2. Constitutive Relations
(2)
2.3. Equation of Heat
(3)
2.4. Strain-Displacement Relation
(4)
(5)
3. Formulation of the Problem
By using the Cartesian coordinates the thermoelastic isotropic medium is considered in half space (
) and in the plane of xy, therefore
.
In non-homogenous medium the parameters
are not constant and depend on the position so we can replace them by
where
are considered to be constants and
is a dimensionless function in space variable.
From Equations (2), (4) and (5) in (1) we get
(6)
(7)
(8)
By assuming that:
, where r is a dimensionless parameter. By using dimensionless variables
By using the dimensionless quantities the Equations (6) and (7) take the form
(9)
(10)
where
,
,
,
,
,
,
, (11)
where
,
,
,
Equation (3) after dimensionless become
(12)
where
,
,
,
We will use the normal mode analysis method to get the exact solution as
(13)
where
is the complex time constant (frequency), i is the imaginary unit, b is the wave number in the x-direction and
are the amplitudes of the functions.
Applying Equation (16) in Equations (13)-(15) we get
(14)
(15)
(16)
where
Eliminating
from Equations (17)-(19) we obtain the sixth order differential equation
(17)
where
,
,
Equation (17) can be factored into
(18)
, (19)
, (20)
, (21)
where
are some constants,
To get the solution of stresses, substitute from Equations (23) and (24) in the dimensionless of Equation (2) we get
(22)
(23)
(24)
where
,
,
.
4. Applications
To determine the coefficients
we consider the boundary conditions at
;
1) The normal stress condition:
(25)
2) The tangential stress condition:
(26)
2) A thermal shock is applied at the boundary:
(27)
Apply these conditions we get
(28)
By using the inverse matrix or Cramer rule we can get
and
.
5. Numerical Results and Discussion
For the purpose of numerical results we choose copper as a thermoelastic material for which we take the following values of the physical constants. Following Dhaliwal and Singh [19] magnesium material was chosen for this purpose.
The previous numerical technique was used for the distribution of the real part of the temperature T, the displacement, components
and the stresses components
, for the problem in the following figures.
Figures 1-4 compare the values in absence and existence of initial stress. The red line denote for the absence and the blue line denote for the existence of initial stress.
Figure 1 displays the horizontal displacement and y, the curves decrease first as
than increase in
than converges to zero after that.
Figure 2 shows the vertical displacement and y, the curves increase till reach the value 0.8 then converges to zero next.
Figure 3 describes
with y, the graph sharply increase intially till y is 0.3 followed by decreasing as
then converges to zero as
.
Figure 4 investigates
with y, the energy dissipation values of
rapidly increasing when
then decreasing when
and get closer to zero next, without energy dissipation the values of
have small increment as
and decrease slowly and converges to zero as
.
Figures 5-10 are carried out in comparing different values for the parameter r, the blue line represent the case
, while the red line when
Figure 5 shows the horizontal displacement u with y, the curve decrease in the range
then converges to zero.
Figure 6 depicts the vertical displacement v withy, the graph increase as
and then it get closer to zero as
.
Figure 7 depicts the temperature T with y, at
the graph increase very slowly in
then slowly decrease then approach after that to zero, at
the temperature rapidly increase when
then converges to zero
Figure 8 displays shear stress
with y. The values of shear stress decrease as
then get closer to zero as y become greater.
Figure 9 illustrates tangential stress
with y. In this figure; the values of the tangential stress decrease in the range
and increase when
then it tends to zero.
Figure 10 illustrates normal stress
versus y, the curve decrease in the range
and then converges to zero, and we notice the values of normal stress are almost equal in the two different values of r.
Figure 1. Horizontal displacement u with y.
Figure 2. Vertical displacement v with y.
Figure 5. Horizontal displacement u with y.
Figure 6. Vertical displacement v with y.
Figure 9. Tangential stress
with y.
6. Conclusion
In this paper we studied the effect of initial stress on functionally graded material under the context of Green-Naghdi theory. We used Normal Mode Method to get an analytical solution and expressions for the distributions of stresses, strains, displacements, current density and temperature. It is observed that these distributions are affected by current and its magnetic field is more than its absence in other words the effect of hall current cannot be neglected. The temperature distribution is always positive and the phenomenon of infinite speed does not exist, according to different values of r the distributions differ.
Nomenclature
Components of stresses
Magnetic permeability
Density
,
Lame’s constants
uDisplacement vector
K Thermal conductivity
T Absolute temperature
Specific heat at constant strain