Dynamic problem in 3 D thermoelastic half-space with rotation in context of GN type II and type III

In this paper, comparison between G-N model of type II (without energy dissipation) and G-N model of type III (with energy dissipation) has been shown in a three dimensional thermoelastic half space with rotation subjected to time dependent heat source on the traction free boundary. Eigenvalue methodology has been adopted to solve the equations resulting from the application of the Normal mode analysis to the non-dimensional coupled equations. Variation of the numerically computed values of thermal stresses and temperature with and without rotation has been illustrated graphically.


Introduction
The inconsistency of heat conduction equation of classical uncoupled theory of thermoelasticity with the experimental results is due to the fact that i. no elastic term is included to account for elastic changes producing heat effects; ii.parabolic nature of heat conduction equation indicating infinite speed of propagation of heat waves means that thermal disturbances (with infinite speed) and elastic disturbances (with finite speed) from the classical theory of thermoelasticity, are coupled together.This suggests that every solution of the equations has a part which extends to infinity.
Biot [1] developed a theory of irreversible thermodynamics and gave a satisfactory derivation of the linear theory of coupled thermoelasticity.In order to obtain a wave type heat conduction equation the concept of generalized thermoelasticity was introduced modifying CCTE and later extended by Dhaliwal and Sherief [3] for anisotropic body, and the uniqueness of the solutions was proved by Ignaczak [2,4].Green-Naghdi proposed a generalized thermoelasticity by modifying the energy equation.There are three types of constitutive relations in G-N model [6][7][8].Type-I leads to classical heat conduction equation.Type-II provides solutions for thermal waves propagating finite speed without energy dissipation (TEWOED) and type-III also confirms propagation of thermal waves of finite speed with energy dissipation (TEWED).Several investigations with these extensions have been studied by Abd-Alla and Abo-Dahab [11], Kar and Kanoria [9] and Yousef [10].Pal et al. [5,14] studied the effect of homogeneity of the surface waves in anisotropic media.
In this paper, the stress distributions and temperature variation has been depicted in an anisotropic triclinic half space for G-N model II and III both with rotation.

Basic equations
For a linear thermoelastic anisotropic body subjected to rotation the field equations are as follows: The equation of motion in the absence of inner heat source: Heat -conduction equation in absence of body force: The Duhamel-Neumann stress-strain relations are: where = ( , , , = 1, 2, 3).Strain -displacement relation: where the differentiation with respect to space variable is denoted by and that with respect to time is denoted by notation.

Formulation of the problem
Let us consider a linear anisotropic thermoelastic half space within ( , , ): 0 ≤ < ∞, 0 ≤ < ∞, 0 ≤ < ∞ for a time dependent heat source on the boundary plane to the surface = 0.The surface = 0 is assumed to be traction free and the body is assumed initially at rest.The components of displacement vectors of three dimensional plane waves in anisotropic elastic medium, are given as: where is the time variable and , ( = 1, 2, 3) denotes the respective orthogonal Cartesian co-ordinate axes.The elastic medium is now considered as rotating uniformly with an angular velocity Ω = Ω , where is the unit vector along the direction of rotation.The equation of motion of the rotating frame contains two additional terms: (Ω×(Ω× ) representing the centripetal acceleration due to time varying motion only and (2Ω× ) representing the coriolis acceleration and = (1, 0, 0).Using Hook's law, the stress-strain-temperature relations can be written as follows: In absence of inner heat source and body force the equation of motion are given as: , + , + , = − Ω + 2Ω .
With the help of Eqs. ( 5) and ( 6), the equations of motions Eq. ( 7) become: The generalized heat-conduction Eq. ( 2) is written as: * + , The following non-dimensional variables are introduced to transform the above equations in non-dimensional form: where is some standard length.Using Eq. ( 10), the non dimensional forms of the Eqs.(8a)-( 9) reduces to (omitting primes for convenience): where: Non -dimensional stress components can be calculated as:

Solution of the vector-matrix differential equation: eigenvalue approach
Applying eigenvalue approach method as in Santra et al. [15] to solve the vector-matrix differential equation (25), we get the characteristic equation of the matrix as: The roots of the Eq. ( 26) are of the form: = ± , ( = 1,2,3,4).
The eigenvector corresponding to the eigenvalue can be calculated as: where , ( , = 1, 2, 3) are given in the Appendix I.
The eigenvector [ = 1(1)8] corresponding to the eigenvalue = [ = 1(1)8] can be calculated from equation ( 27).For our further reference, we use the following notations: As in Lahiri et al. [12], the general solution of Eq. ( 25) which is regular as → +∞ can be written as: where the arbitrary constants are to be determined from the boundary conditions of the problem and because of regularity condition of the solution at +∞ the terms containing exponential of growing nature in the space variables have been neglected.Thus, the field variables can be written from the Eq. ( 29) for ≥ 0 as: (36)

Boundary conditions
To determine the arbitrary constants , boundary conditions are considered at the surface on the half space = 0 as in Sarkar and Lahiri [13].

Numerical analysis
For the purpose of illustrating the problem, we now consider a numerical example for which computational results are presented.Since is complex, we take = + , for studying the effect of wave propagation, we use the following physical parameters in SI units given in the following.
The numerical constants are given by: (1) The nature of the stress component is contraction in nature for G-N model of type II and III with rotation Ω = 25 whereas without rotation stress is parallel to the axis.Stresses for G-N model type II and III are same and parallel in 0 < < 3 and then gradually decreasing.
(2) Fig. 2 shows that the stress component is extensive in nature.The maximum value occurs for G-N model of type III with rotation.Stresses are gradually increasing in 0 < < 0.2 and then decreasing.
(3) Nature of the stress is same as .(4) In Fig. 4 nature of the stress component is same for G-N model II and III without ritation.Small variation can be seen for the case of with rotation in the range 0 < < 0.2.and then stresses are gradually decreasing for all the four cases.
(5) In Fig. 5 nature of is same as for the case of without rotation.And for the case of with rotation the stresses are parallel to the axis except for a small variation in the range 0.3 < < 0.5 in case of G-N model type III. 2 for G-N model II and III respectively.In the range 0 < < 1, the stress component monotonically decreases without rotation whereas it increases initially starting from 0, attains a maximum value and then gradually decreases > 0.3 with rotation.(7) From Fig. 7 we can see that the temperature is gradually increasing in 0 < < 0.2 and then maximum value occurs at = 0.2 for the case of G-N model type III with rotation and after that temperature for all the four cases are gradually decreasing.

Conclusions
Thermal stresses and temperature on a traction free boundary in a half space for G-N model type II as well as type III due to time dependent heat source shows significant dependence on rotation as evident from the above curves.
Various stress components deviate under rotation from their rotation less counter parts.Certain stress component possesses non zero initial value when subject to rotation whereas for temperature a converse effect has been found.

Fig. 7 .
Fig. 7. Variation of with respect to with rotation and without rotation