Normal Mode Analysis to a Poroelastic Half-Space Problem under Generalized Thermoelasticity

The thermo-hydro-mechanical problems associated with a poroelastic half-space soil medium with variable properties under generalized thermoelasticity theory were investigated in this study. By remaining faithful to Biot’s theory of dynamic poroelasticity, we idealized the foundation material as a uniform, fully saturated, poroelastic half-space medium. We first subjected this medium to time harmonic loads consisting of normal or thermal loads, then investigated the differences between the coupled thermohydromechanical dynamic models and the thermo-elastic dynamic models. We used normal mode analysis to solve the resulting nondimensional coupled equations, then investigated the effects that non-dimensional vertical displacement, excess pore water pressure, vertical stress, and temperature distribution exerted on the poroelastic half-space medium and represented them graphically.

The effects of temperature on the behavior of soil is an important engineering concern.Temperature can impact the disposal of high-level radioactive waste, the extraction of oil or geothermal energy, the storage of hot fluids, road subgrades, and furnace foundations.The theory of thermoelasticity was developed in order to analyze these problems associated with thermodynamics.
The classical coupled thermoelasticity proposed by Biot (1956) predicts an infinite speed for heat propagating in elastic media, which is physically impossible.To eliminate this inherent paradox, subsequent generalized thermoelastic theories have been developed by Lord and Shulman (L-S) (1967), Green and Lindsay (G-L) (1972), and Green and Naghdi (G-N) (1991,1992,1993).An extensive theory of thermoelasticity was formulated by Lord and Shulman that introduced a heat flow rate term together with a thermal relaxation time into Fourier's law of heat conduction.By modifying both the energy equation and Duhamel-Neumann's law, Green and Lindsay introduced the concept of allowing two relaxation times.A form of thermoelastic theory that allowed for no energy dissipation was proposed by Green and Naghdi.Other generalized thermoelastic theories include the two-temperature generalized thermoelasticity formulated by Youssef (2006) and the dual-phase-lag thermoelasticity introduced by Tzou (1995).Hetnarski and Ignaczak (1999) examined five concerning theories of the coupled theory of thermoelasticity.Structures 14 (2017) 930-949 A great deal of research attention has been focused on the wave propagation problem in context of generalized thermoelasticity.Singh (2007) solved a two-dimensional homogeneous, isotropic generalized thermoelastic half-space problem in the context of Lord and Shulman theory of generalized thermoelasticity.The propagation of waves in a transversely isotropic micropolar medium possessing thermoelastic properties under three different theories are discussed by Kumar and Gupta (2010).Bijarnia and Singh (2012) investigated the propagation of four different plane waves in a solid half-space with diffusion in the context of the Lord and Shulman theory of generalized thermoelasticity.Lotfy (2014) investigated a two-temperature problem for an isotropic homogeneous elastic half-space under three thermoelastic theories.Lotfy and Hassan (2014) solved a thermoelastic wave propagation problem for a half-space media subject to thermal shock based on the two-temperature theory of thermoelasticity.Abo-Dahab et al. (2015) studied the effects of magnetic field, relaxation times, and rotation on the propagation of surface waves in the context of Green-Lindsay theory.Singh (2016) explored plane wave propagation in a solid half-space which is rotating, transversely isotropic, two-temperature thermoelastic without energy dissipation.Lotfy (2016) introduced the dual-phase-lag (DPL) heat transfer model to solve the problem for a isotropic generalized thermoelastic medium with an internal heat source.

Latin American Journal of Solids and
Wave propagation in a thermoelastic porous media is useful in various engineering fields such as petroleum engineering, chemical engineering, pavement engineering, and nuclear waste management, making it a popular research subject as well.Kumar and Devi (2008) investigated a porous, generalized thermoelastic medium subjected to thermomechanical boundary conditions permeated with various heat sources.Sherief and Hussein (2012) developed a set of governing equations that effectually create a mathematical model of generalized thermoelasticity in poroelastic materials, then, they used this model to solve a thermal shock problem regarding the use of half-space.Abbas and Youssef (2015) solved a two-dimensional problem of a porous material in the context of the fractional order generalized thermoelasticity theory with one relaxation time.Wei et al. (2016) studied the reflection and refraction phenomenon that occurs in an oblique incidence longitudinal wave at a plane interface between an isotropic, homogeneous, thermoelastic medium and a porous thermoelastic medium.Schanz (2009) provided an overview of several poroelastodynamic models and analytical solutions, and also discussed two different numerical methodologies including the finite element method and the boundary element method.
The role that fluid flow plays in thermoelasticity was neglected in the above studies.With this in mind, we conducted the present study because we believe that a thorough knowledge of hydrothermal, hydromechanical, and thermal-mechanical processes, as well as knowledge of the interdependence of these processes, is necessary to accurately model the coupled behavior of fluid-saturated media.The thermo-hydro-mechanical (THM) coupled behavior of unsaturated soil is a very common physical phenomenon in nature and very worthy of further research.
Several previous studies of fluid-saturated media informed our current study.The theoretical developments have matured from a simple isotropic poroelastic theory proposed by Biot into a general approach to managing coupling and material anisotropy (Abousleiman and Ekbote, 2005).Booker and Savvidou (1984) developed an analytic solution for the consolidation of soil around an impermeable rigid, spherical, heat source and a point source.Biot (1977) also formulated a new thermodynamics theory for elastic saturated porous solids.Bai (2006a) solved a double-layered po-Latin American Journal of Solids and Structures 14 (2017) 930-949 rous space subjected to exponential decaying thermal loading, and derived the thermal consolidation of layered, saturated porous half-space to various thermal loads with time.Bai (2006b) established coupled governing equations for one-dimensional thermal consolidation problems regarding saturated porous media subjected to cyclical thermal loading.Based on the thermodynamics of irreversible processes; Bai and Li (2009) investigated a saturated medium with a long cylindrical cavity subjected to variable thermal loading and variable hydrostatic pressure.Lu et al. (2010) investigated a porous elastic medium subjected to a normal force and a thermal source in the context of generalized thermoelastic theory with one relaxation time.Liu et al. (2009;2010a) developed a method to overcome one-dimensional problems for an isotropic saturated poroelastic medium including a cylindrical cavity and spherical cavity subjected to a time-dependent thermal/mechanical shock in the context of thermodynamics theory.That same year, Liu et al. (2010b) solved a thermo-hydroelastodynamic problem in a two-phase porous thermoelastic medium.Bai and Li (2013) developed a spherical space with a spherical cavity based on the irreversible governing equations of saturated porous thermoelastic media and subjected it to variable mechanical and thermal loads.Kumar et al. (2014) investigated the propagation of Lamb waves in solids with layers or half-spaces of inviscid liquid subjected to stress-free boundary conditions in the context of the Green and Lindsay theory of generalized thermoelasticity.Ailawalia and Singla (2016) solved an infinite homogeneous isotropic generalized thermoelastic problem under the Lord and Shulman theory of generalized thermoelasticity.
To date, there have been very few works devoted to investigating thermo-hydro-mechanical problems involving normal or thermal loads via normal mode analysis.In this study, we investigated these thermo-hydro-mechanical problems and the characteristic variable properties in a poroelastic half-space soil medium in the context of generalized thermoelastic theory.We subjected this poroelastic half-space to time harmonic loads comprised of normal or thermal loads.The foundation material, under Biot's dynamic poroelastic theory, was idealized as a uniform, fully saturated poroelastic half-space.The distributions of non-dimensional vertical displacement, excess pore water pressure, vertical stress, and temperature distribution were obtained by means of normal mode analysis.The numerical results were obtained by performing numerical inversion of the transform integrals.These results were then used to analyze the differences between the two methods of coupled thermo-hydro-mechanical dynamic model (THMD) and thermoelastic dynamic model (TMD).

BASIC EQUATIONS
We consider the problem of a saturated, homogenous, isotropic porous elastic half-space.Based on theory of three-dimensional consolidation which proposed by Biot (1941), we assume that solid grains are incompressible because the Biot's effective stress coefficient 1   .The dynamic equation for the motion of thermo-hydromechanical coupling in the absence of body forces can be written by Smith and Booker, Bai and Abousleiman, and Bai as follows (1993, 1997, 2006c): where .
Latin American Journal of Solids and Structures 14 (2017) 930-949 The equation associated with heat conduction for poroelastic soil can be written by Sherief andSaleh, andRam et al. as follows (2005, 2008): where Based on Darcy's law, the equation of mass conservation of water is: where In the absence of a body force or inner heat source, the constitutive governing equation of thermo-hydromechanical can be written by Smith and Booker as follows (1993): The strain-displacement relation equation is:

FORMULATION OF THE PROBLEM
We consider the problem of a saturated, homogenous, isotropic porous elastic half-space soil medium subjected to both normal force and thermal loading.The governing equations are written in the context of the Lord and Shulman model, where the body is under no anybody force.The Cartesian These displacement components have the following forms: Based on Eqs. ( 5) and ( 6), the strain components have the following forms: where e is the cubical dilatation, which is expressed as follows: Based on Eqs. ( 1), (4), and (5): Latin American Journal of Solids and Structures 14 (2017) 930-949 where In order to facilitate the equation, the following non-dimensional quantities are also introduced: where In terms of the non-dimensional quantities in Eq. ( 14), we can transform Eqs. ( 2), (3), (9), and (10) into the following forms (dropping primes for convenience): Accordingly, the constitutive non-dimensional equations are: Latin American Journal of Solids and Structures 14 (2017) 930-949 where Differentiating Eq. ( 15) with respect to x , and Eq. ( 16) with respect to z , then adding them together, we arrive at follows:

NORMAL MODE ANALYSIS
The solutions for the variables that we considered can be decomposed in terms of normal modes in the following form: where  is frequency and a is the wave number in the x -direction.
The roots 1 k  and 2 k  in Eq. ( 31) satisfy the following relations: The solution to Eq. ( 27) can be expressed as follows: where is the solution to the following: The solution of Eq. ( 30), which is bounded as z   , is provided by: where ( , ) i R a  are parameters depending on a and  .Thus, ( ) i e z  takes the form: Similarly: Accordingly, Eqs. ( 38) and ( 39) can rewritten in the following forms: The following relations can be used to obtain the displacement w in terms of Eqs. ( 23) and ( 16): where The solution of Eq. ( 45), which is bounded as z   , is given by: where  is a parameter dependent on a and  .
The boundary conditions at 0 z  are necessary to determine the parameters ( 1,2) i R i  and F .
The boundary conditions are as follows: (1) Surface stress conditions: (2) Thermal boundary condition, where the half-space surface is subjected to a thermal load: (3) Excess pore water pressure boundary condition: By substituting the expressions of the considered variables into the above boundary conditions, we arrive at: Latin American Journal of Solids and Structures 14 (2017) 930-949 where   , x t  denotes the load and thermal source distribution function along the x axis  .By applying Eq. ( 23),   , x t  is expressed as follows:

NUMERICAL RESULTS AND DISCUSSION
It is necessary to know the material constants to secure accurate calculation results.Most of the mechanical, thermal, and hydraulic parameters used in this paper are the same as those used by Bai (2006b): The other constants of the problem are taken as To calculate and verify both the THMD and TMD cases, the density of pore water w  and porosity 0 n are assumed to be zero; thus, the solution of the THMD model reduces to the TMD situa- tion.The equivalent density of the model We divided all 12 figures into three cases.In Case 1, based on THMD and TMD methods, we investigated the non-dimensional vertical displacement, excess pore water pressure, vertical stress, and temperature distribution with two different frequencies 1.6   and 2.5

 
, respectively as shown in Figs.1-4.In Case 2, we considered four different frequencies non-dimensional vertical displacement w w q   , non-dimensional excess pore water pressure p p q   , non-dimensional vertical stress q     , and non-dimensional temperature distribution q     on the surface of the half-space with a distance of x under a uniformly distributed normal force.(For convenience, we dropped the punctuation marks in the figures.)In Case 3 we investigated how the considered variables varied according to the four different frequencies  , respectively.Figure 1 shows where the two methods yield identical results when frequency is 1.6   .When frequency is high, the difference between the two methods are significant.Because our research medium is saturated porous soil, the difference between the two methods is that the THMD method considers the effects of the fluid.When the normal load is constant, low frequency gives the water more time to drain from the saturated porous soil and the difference between the two methods are not significant.When the frequency is high, there is insufficient time for water to drain from the soil and the difference between the two methods are significant.It is important to take into account non-dimensional excess pore water pressure of the soil because it has a considerable effect on ground response.Accordingly, we investigated the variations in excess pore water pressure with x distance under uni-formly distributed normal force, as shown in Fig. 2. The TMD values remained at zero because we did not consider the effects of the fluid in this case.This means that by only using the TMD method, we were unable to assess excess pore water pressure.In actuality, many practical engineering applications are dependent on soil or saturated porous soil such as waste landfill engineering and nuclear waste management.Therefore, non-dimensional excess pore water pressure is an important physical variable in many real life applications.Figure 3 is similar to Fig. 1: As frequency increases, the difference between the two methods becomes more significant.As shown in Fig. 4, frequency  has an increasing effect under both methods., respectively.Our research medium is saturated porous soil, where normal load mainly includes non-dimensional excess pore water pressure and non-dimensional vertical stress.As frequency  increases, less water drains from the medium while the non-dimensional excess pore water pressure increases.When the normal load is constant, non-dimensional vertical stress decreases as the non-dimensional excess pore water pressure increases.Similarly, non-dimensional vertical displacement decreases as frequency  in- creases.The frequency  has an increasing effect upon all the magnitudes of all physical variables., respectively.In this case, we assumed that only a uniformly distributed thermal source was applied to the surface of the half-space.As discussed above, our research medium is saturated porous soil, where normal load mainly includes non-dimensional excess pore water pressure and non-dimensional vertical stress.In this case, only a uniformly distributed thermal source is applied and the normal load is zero.Thus, the absolute value of non-dimensional excess pore water pressure and the absolute value of nondimensional vertical stress are equal.A plus or minus sign only denotes that the medium is under tension or compression.The thermal expansion coefficient of pore water is greater than the thermal expansion coefficient of solid grains.The extent to which water drains from the medium decreases as frequency  increases, while non-dimensional vertical displacement, excess pore water pressure, vertical stress, and temperature distribution all increase as frequency continues to increase.The frequency  has an increasing effect to all the magnitudes of all physical variables.In addition, as observed from the figures, the effect of frequency has the same change trend.

CONCLUSIONS
In this study, we investigated the thermo-hydro-mechanical problems of a poroelastic half-space soil medium with variable properties under the theory of generalized thermoelasticity.We were able to solve these problems using normal mode analysis.Based on the numerical results, we arrived at the following conclusions: (1) Thermo-hydro-mechanical problems for a porous elastic half-space soil medium with variable properties were solved by normal mode analysis, which has been successfully applied in dealing with such complex multi-field coupling problems.The proposed method ultimately yields the corresponding fourth-order characteristic equation and exact solutions for the desired physical variables.
Latin American Journal of Solids and Structures 14 (2017) 930-949 (2) When subjected to a uniformly distributed normal force, the values of non-dimensional excess pore water pressure in the TMD method are consistently zero because the effect of fluid is not considered.This means that the TMD method is not appropriate for determining nondimensional excess pore water pressure.Non-dimensional excess pore water pressure is an important physical variable in real life engineering.However, fluid presence attenuates the response in terms of non-dimensional vertical displacement and vertical stress, which is especially pronounced at higher frequencies.
(3) The frequency of a uniformly distributed thermal load has a pronounced effect on nondimensional vertical displacement, excess pore water pressure, vertical stress, and temperature distribution.Furthermore, as observed from the figures of Case 3, the effect of frequency has the same change trend.(4) Frequency plays an important role in examining the deformation of the body.
(5) As shown in all figures, it can be clearly observed that the figures reveal is consistent with simple harmonic oscillation, which is mainly dominated by the intrinsic nature of the solution obtained when the normal mode analysis method is applied to decouple the equation.

.
Other material parameters were kept the same as those given above to obtain the corresponding TMD solution.
the THMD method (Figs.5-8).The results of both cases indicate variations of Latin American Journal of Solids and Structures 14 (2017) 930-949 Figure 1: Vertical displacement for 1.0 z  with 0 Q  under two methods.

Figure 2 :
Figure 2: Excess pore water pressure for 1.0 z  with 0 Q  under two methods.

Figure 6 :
Figure 5: Vertical displacement for 1.0 z  with 0 Q  under THMD method.
Figure 9: Vertical displacement for 1.0 z  with 0 q  under THMD method.
Figure 11: Vertical stress for 1.0 z  with 0 q  under THMD method. 2 and i is the imaginary unit,