Thermoelastic interactions in an isotropic unbounded medium due to moving heat source using GNIII model 1

In this work, we study a problem of thermoelastic interaction due to moving heat source in an isotropic infinite medium under Green and Naghdi model of type III (GNIII). The form of vector-matrix differential equation in the Laplace transform domain, the basic equations have been written, which is then solved by an eigenvalue technique. The analytical solution in the Laplace transforms domain with eigenvalue approach has been obtained. Numerical results for the displacement, temperature and the stress distributions are represented graphically. Some comparisons have been shown in figures to estimate the effect of heat source velocity and time in the physical quintets.

Many investigators have treated the non-isothermal problems of the theory of elasticity, and so it become important.This is due to their many applications in widely diverse fields.In the extremely high temperature, the nuclear field and temperature gradients originating inside nuclear reactors influence their design and operations.In the high velocities of modern aircraft give rise to aerodynamic heating, which produces intense thermal stresses that reduce the strength of the aircraft structure.The counterparts of our problem in the contexts of the thermoelasticity theories have been considered by using numerical and analytical methods (Mukhopadhyay, 2006;Mukhopadhyay and Kumar, 2009;Jiangong and Tonglong, 2010;Kumar and Chawla, 2010;Abbas et al., 2011;Abo-Dahab and Abbas, 2011;Abbas, 2011;2012;2013;2014a;2014b;2014c;Abbas and Abo-Dahab, 2014;Abbas and Zenkour, 2014).In which, Abbas solved different problems by eigenvalue approach in the Laplace transformation domain.Abbas and his collogues solved one and twodimension problems by finite element method.Kumar and Chawla studied the wave propagation at the boundary surface of elastic layer overlaying a thermoelastic without energy dissipation halfspace.Mukhopadhyay and his collogues used the state space approach for several problems.Chandrasekharaiah andSrinath (1998a, 1998b) studied the Thermoelastic interactions without energy dissipation due to a point and line heat source.He and Cao (2009) considered generalized magneto-thermoelastic problem in thin slim strip subjected to a moving heat source.Youssef (2009Youssef ( , 2010) ) established the thermoelastic interactions in a unbounded medium with cylindrical and spherical cavity subjected to moving heat source.
In the present paper we have applied the technique of eigenvalue approach developed in Das et al. (1997) to solve generalized thermoelastic interaction problem subjected to a moving heat source using GNIII model.The eigenvalue approach gives exact solution in the Laplace domain without any assumed restrictions on the actual physical quantities.The governing equations of the mathematical model is presented when the beam is quiescent first.Laplace transforms techniques with eigenvalue approach are used to get the general solution for any set of boundary conditions.Numerical results are represented graphically.The moving heat source velocity have a significant effect on all distributions.

BASIC EQUATION AND FORMULATION OF THE PROBLEM
Following Othman and Abbas (2012), the system of equations that include the displacement, the stress, the strain and the temperature for a linear, homogenous and isotropic thermoelastic continuum without body forces take the following form: The equations of motion The equation of heat conduction   , , 0 .
The constitutive equations are given by Latin American Journal of Solids and Structures 12 (2015) 1061-1073 where ,   are the Lame's constants; K is the thermal conductivity;  is the density of the medium; and t  is the coefficient of linear thermal expansion; e c is the specific heat at constant strain; t is the time; 0 T is the reference temperature; K  is the material constant T is the temperature; characteristic of the theory; ij  is the Kronecker symbol; i u are the components of displacement vector and Q is the moving heat source; ij  are the components of stress tensor.Let us consider a homogeneous isotropic thermoelastic solid at a uniform reference temperature 0 T occupying the region 0 x  where the x -axis is taken perpendicular to the bounding plane of the half-space pointing inwards.For one-dimensional problem the displacement vector u and temperatures field T can be expressed in the following form: Then the equations ( 1) to (3) take the following form For convenience, we introduce the following non-dimensional variables where 2 1 2 and .
We consider that the medium is subjected to a moving heat source in the following nondimensional form where o Q is constant and  is the delta function.

APPLICATION
We assume that the medium is initially at rest.The undisturbed state is maintained at reference temperature.Then we have We consider boundary conditions of two types: where ( ) H t denotes the Heaviside unit step function and 1 T is a constant.Applying the Laplace transform define by the formula Hence, we obtain the following system of differential equations Latin American Journal of Solids and Structures 12 (2015) 1061-1073 Equations ( 16) and ( 17) can be written in a vector-matrix differential equation as follows Das et al. (1997) , where Now, the eigenvalue approach are used to solve the equation ( 21) as Das et al. (1997) The characteristic equation of the matrix A takes the form The roots of the characteristic equation ( 22) which are also the eigenvalues of matrix A are of the form . The eigenvector , corresponding to eigenvalue R can be calculated as: From equations ( 23), we can easily calculate the eigenvector j X  , corresponding to eigenvalue For further reference, we can write Thus, the complementary solution of equation ( 21) take the form where the terms containing exponentials of growing nature in the space variable x have been discarded due to the regularity condition of the solution at infinity, 1 B and 2 B are constants to be determined from the boundary condition of the problem.
The general solutions V  of the nonhomogeneous system (21) are the sum of the complementary solution c V  of the associated homogeneous system and a particular solution p V  of the nonhomogeneous system.The inhomogeneous terms in (21) contain the exponential function , therefore, the particular solution p V  should be sought in the form of a vector quasi-polynomial where A  is a constant vector.From ( 25), ( 26) and ( 21), the general solutions of the field variables can be written for x and s as: To complete the solution we have to know the constants 1 B and 2 B , by using the boundary conditions (19) for case (I) while the boundary conditions (20) for case (II).

NUMERICAL INVERSION OF THE LAPLACE TRANSFORMS
For the final solution of the temperatures, the displacement, the concentration, the stress and chemical potential distributions in the time domain, we adopt a numerical inversion method based on the Riemann-sum approximation method is used to obtain the numerical results.In this method, any function in the Laplace domain can be inverted to the time domain as Re[ ( , )] Re ( 1) , , where Re is the real part and i is the imaginary number unit.For faster convergence, numerical experiments have shown that the value that satisfies the above relation is Tzou, 1996).

NUMERICAL RESULTS AND DISCUSSION
In the present work, the thermoelastic interactions due to moving heat source under Green and Naghdi of type III model is analyzed by considering an isotropic unbounded medium.The material parameters are given as following Abbas ( 2009) 8.954 10 (kg)(m) , 17.8 10 (K) , 0.5, 1.
Using this data set, the temperature T , displacement u and stress xx  are numerically computed for different values of the distance x and their graphical representation is presented in figures 1-12.As expected, in both the cases, the velocity of moving heat source has a great effect on the distribution of field quantities.The time have a great effect on all distributions.
Case (I): The figures 1-6 is investigating the variation of the non-dimensional temperature, displacement and stress when the traction free and subjected to a thermal shock on the surface .It can be found that the temperature, magnitude of displacement and the magnitude of the stress decreases as the velocity increases before the intersection of the three curves.However, after the intersection, its increases as the velocity increases.Figure 4-6 show the variations of nondimensional temperature, displacement and stress with distance for different value of time   0.2, 0.6, 1.0 t  when the moving heat source velocity   0.2 v  remains constant.It can be found that the temperature, magnitude of displacement and the absolute of the stress increases as the time increases.From figures 1-6, the temperature starts with   1 1 T T   at the origin and increases due to the moving heat source then decreases until attaining zero beyond a wave front for the generalized theory, which agree with the boundary conditions.The displacement component attains maximum negative values and gradually increases until it attains a peak value at a particular location in close proximity to the surface and then continuously decreases to zero.The stress, always starts from the zero value and terminates at the zero value to obey the boundary conditions.
Case (II): The figures 7-12 is investigating the variation of the non-dimensional temperature, displacement and stress when the surface thermally insulation and fixed.Figures 7-9 show the effects of velocity of moving heat source when time remains constant while, figures 10-12 show the effects of time when the velocity of moving heat source remains constant.

CONCLUSIONS
Two cases have been considered in our application.The first one for the traction free and subjected to a thermal shock on the surface while the second case for the surface thermally insulation and fixed.The eigenvalue approach gives exact solution in the Laplace domain without any assumed restrictions on the actual physical quantities.The velocity of moving heat source have a significant effect in on all distributions.

Figure 1 -
Figure 1-3 display the effects of velocity of moving heat source   0.2, 0.4, 0.6 v  when 0.5 t 

Figure 1 :
Figure 1: The variation of temperature with distance for different values of v (Case I).

Figure 2 :
Figure 2: The variation of displacement with distance for different values of v (Case I).

Figure 3 :Figure 4 :
Figure 3: The variation of stress with distance for different values of v (Case I).

Figure 5 :
Figure 5: The variation of displacement with distance for different values of t (Case I).

Figure 6 :Figure 7 :
Figure 6: The variation of stress with distance for different values of t (Case I).Latin American Journal of Solids and Structures 12 (2015) 1061-1073

Figure 8 :
Figure 8: The variation of displacement with distance for different values of v (Case II).

Figure 9 :Figure 10 :
Figure 9: The variation of stress with distance for different values of v (Case II).

Figure 11 :
Figure 11: The variation of displacement with distance for different values of t (Case II).

Figure 12 :
Figure 12: The variation of stress with distance for different values of t (Case II).Latin American Journal of Solids and Structures 12 (2015) 1061-1073