Two-Dimensional Fractional Order Generalized Thermoelastic Porous Material

In the work, a two-dimensional problem of a porous material is considered within the context of the fractional order generalized thermoelasticity theory with one relaxation time. The medium is assumed initially quiescent for a thermoelastic half space whose surface is traction free and has a constant heat flux. The normal mode analysis and eigenvalue approach techniques are used to solve the resulting non-dimensional coupled equations. The effect of the fractional order of the temperature, displacement components, the stress components, changes in volume fraction field and temperature distribution have been depicted graphically.


Latin American Journal of Solids and
diction of behavior of sound-absorbing materials and in the area of exploration geophysics, the steadily growing literature bearing witness to the importance of the subject Pecker and Deresiewiez (1973).
The problem of a fluid-saturated porous material has been studied for many years.A short list of papers pertinent to the present study includes Biot(1941Biot( , 1956)), Gassmann (1951), Biot and Willis (1957), Biot (1962), Deresiewicz and Skalak (1963), Mandl (1964), Nur and Byerlee (1971), Brown and Korringa (1975), Rice and Cleary (1976), Burridge and Keller (1981), Zimmerman et al. (1986Zimmerman et al. ( ,1994)), Berryman and Milton (1991), Thompson and Willis (1991)], Pride et al. (1992), Berryman and Wang (1995), Tuncay and Corapcioglu (1995), Alexander and Cheng (1991), Charlez, P. A., and Heugas, O. (1992), Abousleiman et al. (1998), Ghassemi and Diek (2002), Tod (2003).Eringen (1970) and Nowacki (1966)developed the linear theory of micropolar thermoelasticity which are known as micropolar coupled thermoelasticity to include thermal effects.Goodman and Cowin (1972) established a continuum theory for granular materials, whose matrix material (or skeletal) is elastic and interstices are voids and they introduced the concept of distributed body, which represents a continuum model for granular materials (sand, grain, powder, etc) as well as porous materials (rock, soil, sponge, pressed powder, cork, etc.).Nunziato and Cowin (1979), developed the non-linear theory of elastic materials with void, underlying the basic concept that the bulk density of the material is written as the product of two fields, the density field of the matrix material and the volume fraction field (the ratio of volume occupied by grains to the bulk volume at a point of the material) Kumar and Gupta (2010)].Othman (2007) studied the effect of rotation and relaxation time on a thermal shock problem for a half-space in generalized thermo-viscoelasticity and Othman and Singh (2005) studied the effect of rotation on generalized micropolar thermoelasticity for a half-space under five theories.Youssef (2007)constructed theory of generalized porothermoelasticity which describe the behavior of thermoelastic porous medium in the context of the theory of generalized thermoelasticity with one relaxation time (Lord-Shulman).The energy and the entropy equations have been derived also in general co-ordinates.The uniqueness of the solution for the complete system of the equations of the theorem has been proved by Kumar et al. (2013) and he discussed the plane deformation due to thermal source in fractional order thermoelastic media, while Abbas and Kumar (2014) studied the deformation due to thermal source in micropolar generalized thermoelastic half-space by finite element method.
Recently, a new formula of heat conduction has been considered in the context of the fractional integral operator definition by Youssef (2010).This new consideration generated the fractional order generalized thermoelasticity which was cited by Youssef who approved the uniqueness of its solutions.
Youssef solved one dimensional problem in the context of the fractional order generalized thermoelasticity and discussed the effects of the fractional order parameter on all the studied fields and with Al-Leheabi i(2010).Youssef (2012) solved two-dimensional thermal shock problem of fractional order generalized thermoelasticity with thermal shock.Povstenko (2005) solved a problem of fractional heat conduction equation and associated thermal stress.The counterparts of our problem in the contexts of the thermoelasticity theories have been considered by using analytical and numerical methods Abbas et al. (2002Abbas et al. ( , 2008Abbas et al. ( , 2009Abbas et al. ( , 2011Abbas et al. ( , 2012)).Structures 12 (2015) 1415-1431 In this paper, a two-dimensional problem of a porous material will be considered within the context of the fractional order generalized thermoelasticity theory with one relaxation time.The medium will be assumed initially quiescent for a thermoelastic half space whose surface is traction free and has a constant heat flux.The normal mode analysis and eigenvalue approach techniques will be used to solve the resulting non-dimensional coupled equations.The effect of the fractional order of the temperature, displacement components, the stress and components, changes in volume fraction field distribution will be depicted graphically.

GOVERNING EQUATIONS
For homogeneous, linear and thermally elastic medium with voids and temperature dependent mechanical properties, the basic equations in the context of the Lord and Shulman (1997) model and Cowin and Nunziato (1983) in absence of body forces and heat source are given by Kumar and Devi (2011).
The equations of motion Kumar and Devi (2011): The generalized heat conduction equation Youssef (2010) and Kumar and Devi (2011): where the fractional integral operator defined as follows Youssef (2010): and ( ) Γ α is the Gamma function.
The constitutive equations The cubical dilatation where ρ is the mass density, T the temperature change of a material particle, o T the reference uniform temperature of the body, i u the displacement vector components, ij e the strain tensor; ij σ the Latin American Journal of Solids and Structures 12 (2015) 1415-1431 stress tensor, E c the specific heat at constant strain, γ the thermal elastic coupling tensor in which ( ) β ξ and ψ are the material constants due to presence of voids and φ is the change in volume fraction field of voids.

Formulation and solution of the problem
We consider an isotropic, homogenous and elastic body with voids in two-dimensional fills the region subjected to a time-dependent heat source and traction free on the surface x 0 = .The governing equations will be written in the context of Lord and Shulman model when the body has no heat sources or anybody forces, and we will use the Cartesian co-ordinates ( ) , , x y z and the components of the displacement ( ) ,0, = i u u w to write them as follows: The equations of motion are in the forms and The heat conduction equation ( ) The heat flux equation in x-direction The constitutive relations are ( ) and The cubical dilatation For our convenience, the following non-dimensional variables and notations are used: In terms of the non-dimensional quantities defined above, the governing equations will be reduce to (dropping the dashed for convenience) ( ) where

SOLUTION OF THE VECTOR-MATRIX DIFFERENTIAL EQUATION
Let us now proceed to solve equation ( 30) by the eigenvalue approach proposed by Das et al. (2009).The characteristic equation of the matrix W takes the form where

NQ AFHR ADNR A HS ADMS F
The roots of the characteristic equation (33) which are also the eigenvalues of matrix W in the form Latin American Journal of Solids and Structures 12 (2015) 1415-1431 The eigenvector [ ] which are corresponding to eigenvalue λ can be calculated as From equations ( 36)-( 40) we can easily calculate the eigenvector j X , corresponding to eigenvalue , 1, 2,3, 4,5, 6, 7,8.= j j λ For further reference, we shall use the following notations: The solution of equation ( 30) can be given by: 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 2 3 , , A A A and 4 A are constants to be determined from the boundary condition of the problem.Thus, the field variables can be written for 0, 0, , ≥ > − ∞ ≤ ≤ ∞ x t z as: ( ) 4 5 1 ( , , ) ,  ( , , ) , ) , ) , To complete the solution we have to know the constants 1 2 3 , , A A A and 4 A , so we will use the following boundary conditions.

APPLICATION
We will consider that the bounding plane of the medium 0 = x traction free and has a constant heat flux with constant strength.
Thus, the appropriate boundary conditions are and which gives where o q is the strength of the heat flux and it is constant From the boundary conditions ( 50), ( 51) and ( 53), we obtain where the element of matrix rs H are given by: ( )( ) ) ) )

NUMERICAL RESULTS AND DISCUSSIONS
Following Kumar and Devi (2011), magnesium material was chosen for purposes of numerical evaluations.The physical data are given as In figure 1, the fractional order parameterα has a significant effect on the temperature distribution, where increasing onα causes increasing on T and the rate of change of T with respect to x also increases when α increases which is compatible with the definition of the thermal conductivity.
In figures 2 and 3, the fractional order parameterα has a significant effect on the displacement u and w distributions, where increasing onα causes increasing on the absolute values of u and w, and the rate of change of them with respect to x also increase when α increases which is compatible with the definition of the thermal conductivity.
Figure 4 shows the variation of change in volume fraction field respect to x with different value of the fractional order parameterα .It is seen that, the volume fraction starts with its maximum value at the origin and decreases until attaining zero.The fractional order parameterα has a significant effect on the change in volume fraction field of voids distributionφ , where decreases with the decrease in the value of fractional parameterα .In figure 5, the fractional order parameterα has a significant effect on strain distribution e, where increasing onα causes increasing on e , and the rate of change of e with respect to x also increases when α increases which is compatible with the definition of the thermal conductivity.
In figures 6-8, the fractional order parameterα has significant effects on all components of the stress distribution, where increasing onα causes increasing the absolute values of the stresses, and the rate of change of them with respect to x also increase when α increases which is compatible with the definition of the thermal conductivity.

CONCLUSION
In this work, the effect of the fractional order of the temperature, displacement components, the stress components, changes in volume fraction field and temperature distribution have been studying for a two-dimensional problem of a porous material is considered within the context of the fractional order generalized thermoelasticity theory with one relaxation time.We found that, the fractional order parameter has significant effects on all the studied fields and the results supporting the definition of the classification of the thermal conductivity of the materials to three types; weak, normal and super conductivity.
14)Latin American Journal ofSolids and Structures 12 (2015) 1415-1431 Figures 1-8 represent the temperature distribution, displacement u distribution, displacement w distribution, the change in volume fraction field of voids distribution, the strain distribution, the stress xx σ distribution, the stress xz σ distribution and the stress zz σ distribution respectively at constant time 2.5 = t and constant 0.5 = z with different values of the fractional parameter 0.5, 1.0, 1.5 = α which express for the weak thermal conductivity, normal thermal conductivity and super thermal conductivity respectively.

Figure 1 :
Figure 1: The temperature distribution with different value of the fractional parameter.

Figure 2 :
Figure 2: The displacement u distribution with different value of the fractional parameter.

Figure 3 :
Figure 3: The displacement w distribution with different value of the fractional parameter.

Figure 4 :
Figure 4: The change in volume fraction field of voids distribution with different value of the fractional parameter.

Figure 5 :
Figure 5: The strain distribution with different value of the fractional parameter.

Figure 6 :
Figure 6: The stress xx σ distribution with different value of the fractional parameter.

Figure 7 :
Figure 7: The stress xz σ distribution with different value of the fractional parameter.

Figure 8 :
Figure 8: The stress zz σ distribution with different value of the fractional parameter.