**Debkumar Ghosh**^{1}
,
**Ibrahim A. Abbas**^{3}
,
**Abhijit Lahiri**^{2}

^{1, 2}Department of Mathematics, Jadavpur University, Kolkata, 700032, India

^{3}Department of Mathematics, Sohag University, Sohag, Egypt

^{1}Corresponding author

Mathematical Models in Engineering, Vol. 3, Issue 1, 2017, p. 27-40.
https://doi.org/10.21595/mme.2017.18236

Received 7 February 2017; accepted 11 March 2017; published 30 June 2017

Copyright © 2017 JVE International Ltd.

This paper concerns with the study of wave propagation in fibre reinforced anisotropic half space under the influence of temperature and hydrostatic initial stress. Lord-Shulman theory is applied to the heat conduction equation. The resulting equations are written in the form of vector matrix differential equation by using Normal Mode technique, finally which is solved by Eigen value approach.

**Keywords: ** eigenvalue, generalized thermoelasticity, normal mode, vector-matrix differential equation.

Fibre-reinforced composite(FRC) materials are usually low weight and high strength used in construction engineering. The physical property of FRC material is governed by the theory of elasticity for different materials with the direction along the direction of fibre. Green [1] studied wave propagation in anisotropic elastic plates. Abbas and Othman [2] discussed the distribution of wave propagation under hydrostatic initial stress of fibre-reinforced materials in anisotropic half-space. Baylies and Green [3] analyse the flexural waves in fibre-reinforced laminated plates. Rogerson [4] discussed effect of penetration in a six-ply composite laminates.

Most of the thermoelasticity and generalized thermoelasticity (coupled or uncoupled) problems have been solved by potential function approach. This method is not always suitable as discussed by Dhaliwal and Sherief [5] and Sherief and Anwar [6]. These may be summarized by the initial conditions and the boundary conditions for physical problems which are directly concern with the material quantities under consideration and not with the potential function. Also, the potential function representations are not convergent always while the physical problems in natural variables constitute convergent solution. So, the alternative method of potential function approach is eigenvalue approach. In this method, we obtain a vector-matrix differential equation from the basic equations which reduces finally to an algebraic eigenvalue problem and the solutions for the field variables are obtained by determining the eigenvalues and eigenvectors from the corresponding coefficient matrix. In this theory, body forces and/or heat sources are also accommodated as in Das and Lahiri [7], Bachher et al. [8]. Now, two different models of generalized thermoelastisities are extensively used. One is Lord and Shulman (L-S) [9] theory and the other is Green and Lindsay (G-L) [10] theory. Introducing one relaxtion time parameter in L-S theory the heat conduction equation becomes hyperbolic type without violating conventional Fourier’s law. Whereas the G-L theory modified the heat conduction equation as well as the equation of motion in coupled thermoelasticitywo relaxation time parameters. There are other three models (Model I, II and III by Green and Nagdhi [11-13]) for generalized thermoelasticity concerned to the theory of with or without energy dissipation.

The stress-strain relation and the governing equations of motion without body forces and heat sources are written as follow:

(1)

${\sigma}_{ij,j}-P{\omega}_{ij,j}=\rho {\ddot{u}}_{i},$
(2)

${\sigma}_{ij}=\lambda {e}_{kk}{\delta}_{ij}+2{\mu}_{T}{e}_{ij}+\alpha \left({a}_{k}{a}_{m}{e}_{km}{\delta}_{ij}+{a}_{i}{a}_{j}{e}_{kk}\right)$
$+2\left({\mu}_{L}-{\mu}_{T}\right)\left({a}_{i}{a}_{k}{e}_{kj}+{a}_{j}{a}_{k}{e}_{ki}\right)+\beta {a}_{k}{a}_{m}{e}_{km}{a}_{i}{a}_{j}-{\beta}_{ij}\left(T-{T}_{0}\right){\delta}_{ij},i,j,k,m=\mathrm{1,2},3,$

(3)

${e}_{ij}=\frac{1}{2}\left({u}_{i,j}+{u}_{j,i}\right),$
(4)

${\omega}_{ij}=\frac{1}{2}\left({u}_{j,i}-{u}_{i,j}\right),$
(5)

${K}_{ij}{T}_{ij}=\rho {c}_{e}\left(\dot{T}+{t}_{0}\ddot{T}\right)+{T}_{0}\left({\dot{u}}_{i,j}{\ddot{u}}_{i,j}\right),i,j=\mathrm{1,2},3.$
We consider the problem of a elastic half-space $(x\ge 0)$ in fibre-reinforced anisotropic material with $a\equiv ({a}_{1},{a}_{2},{a}_{3})$ where ${a}_{1}^{2}+{a}_{2}^{2}+{a}_{3}^{2}=1$ as in I. A. Abbas [14], where the displacements are given:

(6)

$u={u}_{x}=u\left(x,y,t\right),v={u}_{y}=v\left(x,y,t\right),w={u}_{z}=0.$
We consider the direction of fibre as $a\equiv (\mathrm{1,0},0)$ with $x$-axis as prefered direction, and Eqs. (1-5), reduces as given:

(7)

${\sigma}_{11}=\left(\lambda +2\alpha +4{\mu}_{L}-2{\mu}_{T}+\beta \right)\frac{\partial u}{\partial x}+\left(\lambda +\alpha \right)\frac{\partial v}{\partial y}-{\beta}_{11}\left(T-{T}_{0}\right),$
(8)

${\sigma}_{22}=\left(\lambda +2{\mu}_{T}\right)\frac{\partial v}{\partial y}+\left(\alpha +\lambda \right)\frac{\partial u}{\partial x}-{\beta}_{22}\left(T-{T}_{0}\right),$
(9)

${\sigma}_{12}={\mu}_{L}\left(\frac{\partial v}{\partial x}+\frac{\partial u}{\partial y}\right),$
(10)

${A}_{11}\frac{{\partial}^{2}u}{\partial {x}^{2}}+\left({A}_{12}+{\mu}_{L}-\frac{P}{2}\right)\frac{{\partial}^{2}v}{\partial x\partial y}+\left({\mu}_{L}+\frac{P}{2}\right)\frac{{\partial}^{2}u}{\partial {y}^{2}}-{\beta}_{11}\frac{\partial T}{\partial x}=\rho \frac{{\partial}^{2}u}{\partial {t}^{2}},$
(11)

${A}_{22}\frac{{\partial}^{2}v}{\partial {y}^{2}}+\left({A}_{12}+{\mu}_{L}-\frac{P}{2}\right)\frac{{\partial}^{2}u}{\partial x\partial y}+\left({\mu}_{L}+\frac{P}{2}\right)\frac{{\partial}^{2}v}{\partial {x}^{2}}-{\beta}_{22}\frac{\partial T}{\partial y}=\rho \frac{{\partial}^{2}v}{\partial {t}^{2}},$
(12)

${K}_{11}\frac{{\partial}^{2}T}{\partial {x}^{2}}+{K}_{22}\frac{{\partial}^{2}T}{\partial {y}^{2}}=\left(\frac{\partial}{\partial t}+{t}_{0}\frac{{\partial}^{2}}{\partial {t}^{2}}\right)\left(\rho {c}_{e}T+{T}_{0}{\beta}_{11}\frac{\partial u}{\partial x}+{T}_{0}{\beta}_{22}\frac{\partial v}{\partial y}\right),$
with:

${A}_{11}=\lambda +2\left(\alpha +{\mu}_{L}\right)+4\left({\mu}_{L}-{\mu}_{T}\right)+\beta ,{A}_{12}=\alpha +\lambda ,{A}_{22}=\lambda +2{\mu}_{T},$

${\beta}_{11}=\left(2\lambda +3\alpha +4{\mu}_{L}-2{\mu}_{T}+\beta \right){\alpha}_{11}+\left(\lambda +\alpha \right){\alpha}_{22},{\beta}_{22}=\left(2\lambda +\alpha \right){\alpha}_{11}+\left(\lambda +2{\mu}_{T}\right){\alpha}_{22},$

where ${\alpha}_{11}$, ${\alpha}_{22}$ are linear thermal expansion coeeficients.

To transform the above governing equations in non-dimensional forms, we introduce the non-dimensional variables as follows:

(13)

$\left({x}^{\text{'}},{y}^{\text{'}},{u}^{\text{'}},{v}^{\text{'}}\right)={c}_{1}\chi \left(x,y,u,v\right),{t}^{\text{'}}={c}_{1}^{2}\chi t,{T}^{\text{'}}=\frac{{\beta}_{11}\left(T-{T}_{0}\right)}{\rho {c}_{{1}^{2}}},\chi =\frac{\rho {c}_{e}}{{K}_{11}},$
$\left({{\sigma}^{\text{'}}}_{11},{{\sigma}^{\text{'}}}_{12},{{\sigma}^{\text{'}}}_{22}\right)=\frac{1}{\rho {c}_{1}^{2}}\left({\sigma}_{11},{\sigma}_{12},{\sigma}_{22}\right),{c}_{1}^{2}=\frac{{A}_{11}}{\rho}.$

Using non-dimensional Eq. (13), the governing equations reduces to (eleminating primes for convenience):

(14)

${\sigma}_{11}=\frac{\partial u}{\partial x}+{B}_{1}\frac{\partial v}{\partial y}-T,$
(15)

${\sigma}_{22}={B}_{1}\frac{\partial u}{\partial x}+{B}_{2}\frac{\partial v}{\partial y}-{B}_{3}T,$
(16)

${\sigma}_{12}={B}_{4}\left(\frac{\partial v}{\partial x}+\frac{\partial u}{\partial y}\right),$
(17)

$\frac{{\partial}^{2}u}{\partial {x}^{2}}+\left({B}_{1}+{B}_{4}-\frac{{R}_{p}}{2}\right)\frac{{\partial}^{2}v}{\partial x\partial y}+\left({B}_{4}+\frac{{R}_{p}}{2}\right)\frac{{\partial}^{2}u}{\partial {y}^{2}}-\frac{\partial T}{\partial x}=\frac{{\partial}^{2}u}{\partial {t}^{2}},$
(18)

${B}_{2}\frac{{\partial}^{2}v}{\partial {y}^{2}}+\left({B}_{1}+{B}_{4}-\frac{{R}_{p}}{2}\right)\frac{{\partial}^{2}u}{\partial x\partial y}+\left({B}_{4}+\frac{{R}_{p}}{2}\right)\frac{{\partial}^{2}v}{\partial {x}^{2}}-{B}_{3}\frac{\partial T}{\partial y}=\frac{{\partial}^{2}v}{\partial {t}^{2}},$
(19)

$\frac{{\partial}^{2}T}{\partial {x}^{2}}+{\epsilon}_{1}\frac{{\partial}^{2}T}{\partial {y}^{2}}=\left(\frac{\partial}{\partial t}+{t}_{0}\frac{{\partial}^{2}}{\partial {t}^{2}}\right)\left(T+{\epsilon}_{2}\frac{\partial u}{\partial x}+{\epsilon}_{3}\frac{\partial v}{\partial y}\right),$
where:

$\left({B}_{1},{B}_{2},{B}_{3}\right)=\frac{1}{{A}_{11}}\left({A}_{12},{A}_{22},{\mu}_{1}\right),{B}_{3}=\frac{{\beta}_{22}}{{\beta}_{11}},{R}_{p}=\frac{P}{{A}_{11}},$

$\left({\epsilon}_{2},{\epsilon}_{3}\right)=\frac{{T}_{0}{\beta}_{11}}{{A}_{11}\rho {c}_{e}}\left({\beta}_{11},{\beta}_{22}\right),{\epsilon}_{1}=\frac{{K}_{11}}{{K}_{22}}.$

$\left({\epsilon}_{2},{\epsilon}_{3}\right)=\frac{{T}_{0}{\beta}_{11}}{{A}_{11}\rho {c}_{e}}\left({\beta}_{11},{\beta}_{22}\right),{\epsilon}_{1}=\frac{{K}_{11}}{{K}_{22}}.$

For the solution of the Eqs. (14-19), physical variables can be decomposed using normal modes Eq. (20) in the following form:

(20)

$\left[u,v,T,{\sigma}_{11},{\sigma}_{12},{\sigma}_{22}\right]\left(x,y,t\right)=\left[{u}^{*},{v}^{*},{T}^{*},{\sigma}_{11}^{*},{\sigma}_{12}^{*},{\sigma}_{22}^{*}\right]\left(x\right){e}^{\omega t+iay},$
where $i=\sqrt{-1}$, $\omega $ is the angular frequency and a is the wave number along $x$-axis.

Using Eq. (20), Eqs. (14-19) reduces to omitting ‘*’ for convenience:

(21)

${\sigma}_{11}=\frac{du}{dx}+\left({B}_{1}ia\right)v-T,$
(22)

${\sigma}_{22}={B}_{1}\frac{du}{dx}+\left({B}_{2}ia\right)v-{B}_{3}T,$
(23)

${\sigma}_{12}={B}_{4}\frac{dv}{dx}+\left({B}_{4}ia\right)u,$
(24)

$\frac{{d}^{2}u}{d{x}^{2}}=M41.u+0.v+0.T+0.{u}^{\text{'}}+M45.{v}^{\text{'}}+{T}^{\text{'}},$
(25)

$\frac{{d}^{2}v}{d{x}^{2}}=0.u+M52.v+M53.T+M54.{u}^{\text{'}}+0.{v}^{\text{'}}+0.{T}^{\text{'}},$
(26)

$\frac{{d}^{2}T}{d{x}^{2}}=0.u+M62.v+M63.T+M64.{u}^{\text{'}}+0.{v}^{\text{'}}+0.{T}^{\text{'}},$
where:

$M41={a}^{2}\left[{B}_{4}+\frac{{R}_{p}}{2}\right]+{\omega}^{2},M45=-ia\left[{B}_{1}+{B}_{4}-\frac{{R}_{p}}{2}\right],M52=\frac{{a}^{2}{B}_{2}+{\omega}^{2}}{{B}_{4}+\frac{{R}_{p}}{2}},$

$M53=\frac{\frac{{R}_{p}}{2}-{B}_{1}-{B}_{4}}{{B}_{4}+\frac{{R}_{p}}{2}},M54=\frac{ia{B}_{3}}{{B}_{4}+\frac{{R}_{p}}{2}},M62=ia{\epsilon}_{3}(\omega +{t}_{0}{\omega}^{2}),$

$M63=\omega +{t}_{0}{\omega}^{2}+{\epsilon}_{2}{a}^{2},M64={\epsilon}_{2}^{2}\left(\omega +{t}_{0}{\omega}^{2}\right).$

$M53=\frac{\frac{{R}_{p}}{2}-{B}_{1}-{B}_{4}}{{B}_{4}+\frac{{R}_{p}}{2}},M54=\frac{ia{B}_{3}}{{B}_{4}+\frac{{R}_{p}}{2}},M62=ia{\epsilon}_{3}(\omega +{t}_{0}{\omega}^{2}),$

$M63=\omega +{t}_{0}{\omega}^{2}+{\epsilon}_{2}{a}^{2},M64={\epsilon}_{2}^{2}\left(\omega +{t}_{0}{\omega}^{2}\right).$

Eqs. (24-26) can be written in the form of vector-matrix differential equation as [2, 8]:

(27)

$\frac{d\overrightarrow{W}}{dx}=\overrightarrow{A}\overrightarrow{W},$
where $\overrightarrow{W}$ =${\left[uvT{u}^{\text{'}}{v}^{\text{'}}{T}^{\text{'}}\right]}^{T}$ and $A=\left[\begin{array}{ll}{L}_{11}& {L}_{12}\\ {L}_{21}& {L}_{22}\end{array}\right]$. Where ${L}_{11}$ is null matrix and ${L}_{12}$ identity matrix of order 3×3 respectively and ${L}_{21}$ and ${L}_{22}$ are given by:

${L}_{11}=\left(\begin{array}{lll}0& 0& 0\\ 0& 0& 0\\ 0& 0& 0\\ & & \end{array}\right),{L}_{12}=\left(\begin{array}{lll}1& 0& 0\\ 0& 1& 0\\ 0& 0& 1\\ & & \end{array}\right),$

${L}_{21}=\left(\begin{array}{ccc}M41& 0& 0\\ 0& M52& M53\\ 0& M62& M63\\ & & \end{array}\right),{L}_{22}=\left(\begin{array}{ccc}0& M45& 1\\ M54& 0& 0\\ M64& 0& 0\\ & & \end{array}\right).$

${L}_{21}=\left(\begin{array}{ccc}M41& 0& 0\\ 0& M52& M53\\ 0& M62& M63\\ & & \end{array}\right),{L}_{22}=\left(\begin{array}{ccc}0& M45& 1\\ M54& 0& 0\\ M64& 0& 0\\ & & \end{array}\right).$

To solve the vector-matrix differential Eq. (27), we apply the method of eigenvalue approach,

The characteristic equation of the matrix $\overrightarrow{A}$ is given by:

(28)

$\left|\left.A-\lambda I\right|\right.=0.$
The roots of the characteristic Eq. (28) are $\lambda ={\lambda}_{i}$, $i=$ 1, 2, 3 which are of the form $\lambda =\pm {\lambda}_{1}$, $\lambda =\pm {\lambda}_{2}$ and $\lambda =\pm {\lambda}_{3}$ and they are also eigenvalues of the matrix.

The eigenvector, $\overrightarrow{W}$ corresponding to the eigenvalue $\lambda $ can obtained as:

(29)

${\overrightarrow{X}}_{\lambda}={\left[{\delta}_{1}{\delta}_{2}{\delta}_{3}\lambda {\delta}_{1}\lambda {\delta}_{2}\lambda {\delta}_{3}\right]}^{T},$
where ${\delta}_{1}={a}_{2}{b}_{3}-{a}_{3}{b}_{2}$, ${\delta}_{2}={a}_{3}{b}_{1}-{a}_{1}{b}_{3}$, ${\delta}_{3}={a}_{1}{b}_{2}-{a}_{2}{b}_{1}$.

As in Lahiri et al. [7], the general solution of Eq. (27) which is regular as can be written as:

(30)

$\overrightarrow{W}={\sum}_{i=1}^{3}\u200d{A}_{i}{X}_{i}{e}^{-{\lambda}_{i}x},x\ge 0.$
Hence the field variables can be written as the following:

$u={A}_{1}{x}_{11}{e}^{-{\lambda}_{1}x}+{A}_{2}{x}_{21}{e}^{-{\lambda}_{2}x}+{A}_{3}{x}_{31}{e}^{-{\lambda}_{3}x},$

$v={A}_{1}{x}_{12}{e}^{-{\lambda}_{1}x}+{A}_{2}{x}_{22}{e}^{-{\lambda}_{2}x}+{A}_{3}{x}_{32}{e}^{-{\lambda}_{3}x},$

$T={A}_{1}{x}_{13}{e}^{-{\lambda}_{1}x}+{A}_{2}{x}_{23}{e}^{-{\lambda}_{2}x}+{A}_{3}{x}_{33}{e}^{-{\lambda}_{3}x}.$

$v={A}_{1}{x}_{12}{e}^{-{\lambda}_{1}x}+{A}_{2}{x}_{22}{e}^{-{\lambda}_{2}x}+{A}_{3}{x}_{32}{e}^{-{\lambda}_{3}x},$

$T={A}_{1}{x}_{13}{e}^{-{\lambda}_{1}x}+{A}_{2}{x}_{23}{e}^{-{\lambda}_{2}x}+{A}_{3}{x}_{33}{e}^{-{\lambda}_{3}x}.$

The simplified form of Eqs. (21-23) can be written as:

$\begin{array}{l}{\sigma}_{11}={A}_{1}{R}_{11}\left(x\right)+{A}_{2}{R}_{12}\left(x\right)+{A}_{3}{R}_{13}\left(x\right),\\ {\sigma}_{22}={A}_{1}{R}_{21}\left(x\right)+{A}_{2}{R}_{22}\left(x\right)+{A}_{3}{R}_{23}\left(x\right),\\ {\sigma}_{33}={A}_{1}{R}_{31}\left(x\right)+{A}_{2}{R}_{32}\left(x\right)+{A}_{3}{R}_{33}\left(x\right),\end{array}$

where:

${R}_{11}\left(x\right)=\left[-{\lambda}_{1}{x}_{11}+{B}_{1}ia{x}_{12}-{x}_{13}\right]{e}^{-{\lambda}_{1}x},$

${R}_{12}\left(x\right)=\left[-{\lambda}_{2}{x}_{21}+{B}_{1}ia{x}_{22}-{x}_{23}\right]{e}^{-{\lambda}_{2}x},$

${R}_{13}\left(x\right)=\left[-{\lambda}_{3}{x}_{31}+{B}_{1}ia{x}_{32}-{x}_{33}\right]{e}^{-{\lambda}_{3}x},$

${R}_{21}\left(x\right)=\left[-{\lambda}_{1}{B}_{1}{x}_{11}+{B}_{2}ia{x}_{12}-{B}_{3}{x}_{13}\right]{e}^{-{\lambda}_{1}x},$

${R}_{21}\left(x\right)=\left[-{\lambda}_{2}{B}_{1}{x}_{21}+{B}_{2}ia{x}_{22}-{B}_{3}{x}_{23}\right]{e}^{-{\lambda}_{2}x},$

${R}_{21}\left(x\right)=\left[-{\lambda}_{3}{B}_{1}{x}_{31}+{B}_{2}ia{x}_{32}-{B}_{3}{x}_{33}\right]{e}^{-{\lambda}_{3}x},$

${R}_{31}\left(x\right)=\left[{B}_{4}ia{x}_{11}-{\lambda}_{1}{x}_{12}\right]{e}^{-{\lambda}_{1}x},$

${R}_{32}\left(x\right)=\left[{B}_{4}ia{x}_{21}-{\lambda}_{2}{x}_{22}\right]{e}^{-{\lambda}_{2}x},$

${R}_{33}\left(x\right)=\left[{B}_{4}ia{x}_{31}-{\lambda}_{3}{x}_{32}\right]{e}^{-{\lambda}_{3}x}.$

${R}_{12}\left(x\right)=\left[-{\lambda}_{2}{x}_{21}+{B}_{1}ia{x}_{22}-{x}_{23}\right]{e}^{-{\lambda}_{2}x},$

${R}_{13}\left(x\right)=\left[-{\lambda}_{3}{x}_{31}+{B}_{1}ia{x}_{32}-{x}_{33}\right]{e}^{-{\lambda}_{3}x},$

${R}_{21}\left(x\right)=\left[-{\lambda}_{1}{B}_{1}{x}_{11}+{B}_{2}ia{x}_{12}-{B}_{3}{x}_{13}\right]{e}^{-{\lambda}_{1}x},$

${R}_{21}\left(x\right)=\left[-{\lambda}_{2}{B}_{1}{x}_{21}+{B}_{2}ia{x}_{22}-{B}_{3}{x}_{23}\right]{e}^{-{\lambda}_{2}x},$

${R}_{21}\left(x\right)=\left[-{\lambda}_{3}{B}_{1}{x}_{31}+{B}_{2}ia{x}_{32}-{B}_{3}{x}_{33}\right]{e}^{-{\lambda}_{3}x},$

${R}_{31}\left(x\right)=\left[{B}_{4}ia{x}_{11}-{\lambda}_{1}{x}_{12}\right]{e}^{-{\lambda}_{1}x},$

${R}_{32}\left(x\right)=\left[{B}_{4}ia{x}_{21}-{\lambda}_{2}{x}_{22}\right]{e}^{-{\lambda}_{2}x},$

${R}_{33}\left(x\right)=\left[{B}_{4}ia{x}_{31}-{\lambda}_{3}{x}_{32}\right]{e}^{-{\lambda}_{3}x}.$

Considering the problem of a half-space $\varphi $, defined as follows:

$\varphi =\left(x,y,z\right):0\le x\le \infty ,-\infty \le y\le \infty ,-\infty \le z\le \infty .$

In order to determine the arbitrary constants ${A}_{i}^{\text{'}}s$, $i=$ 1, 2, 3, we consider the boundary conditions as follows.

a) Mechanical Boundary condition:

For stress-free surface $x=$0, ${\sigma}_{11}=0$, ${\sigma}_{12}=0.$

b) Thermal Boundary condition:

(31)

$\nu T-\frac{dT}{dx}=r,$
where $\nu $ is Biot’s number.

a) Mechanical Boundary condition:

For stress-free surface $x=$0, ${\sigma}_{11}=-{P}_{1}+{P}_{2}{e}^{\omega t+iay}$, ${\sigma}_{12}=0$.

b) Thermal Boundary condition:

(32)

$T={P}_{3}{e}^{\omega t+iay}.$
Fig. 1 represents distribution of normal stress ${\sigma}_{11}$ for $y=$0.3.

For fixed time $t$, ${\sigma}_{11}$ gradually increases as $x$ increases. For fixed $x$ numerical values of ${\sigma}_{11}$ gradually decreases as t increases.

Fig. 2 represents distribution of normal stress ${\sigma}_{12}$ for $y=$0.2.

For fixed time $t$, ${\sigma}_{12}$ gradually decreases as $x$ increases. For fixed x numerical values of ${\sigma}_{12}$ gradually increases as $t$ increases.

Fig. 3 represents distribution of normal stress ${\sigma}_{22}$ for $y=$0.5.

For fixed time $t$, ${\sigma}_{22}$ gradually decreases as $x$ increases. For fixed $x$ numerical values of ${\sigma}_{12}$ gradually increases as $t$ increases.

**Fig. 1.**
Stress component ${\sigma}_{11}$ at $y=$0.3 for different values of $t$ verses $x$

**Fig. 2.**
Stress component ${\sigma}_{12}$ at $y=$0.2 for different values of $t$ verses $x$

Fig. 4 represents distribution of normal stress ${\sigma}_{11}$ for different values of $x$ and $y$ for fixed $t=$0.1 and $\omega =$0.5.

$x$ numerical value of ${\sigma}_{11}$ gradually decreases as $y$ increases. For fixed $y$ the numerical value of ${\sigma}_{11}$ gradually increases as $x$ increases. ${\sigma}_{11}$ is maximum when $x=$1 and $y=$0.

Fig. 5 represent distribution of normal stress ${\sigma}_{12}$ for different values of $x$ and $y$ for fixed $t=$0.4 and $\omega =$5.

For fixed $x$ numerical value of ${\sigma}_{12}$ gradually decreases as $y$ increases. For fixed $y$ the numerical value of ${\sigma}_{12}$ gradually decreases as $x$ increases. Significant changes occur in the region 0.2$\le x\le $0.6 and 0.6$\le y\le $1.0.

Fig. 6 represent distribution of stress component ${\sigma}_{22}$ at for different values of $x$ and $y$ for fixed $t=$0.1 and 0.1.

**Fig. 4. **
Stress component ${\sigma}_{11}$ at $t=$0.1 and $\omega =$0.5 verses $x$ and $y$

**Fig. 5.**
The variation of stress component ${\sigma}_{12}$ at $t=$0.4 and $\omega =$3 verses $x$ and $y$

**Fig. 6.**
Stress component ${\sigma}_{22}$ at $t=$0.1 and $\omega =$0.1 verses $x$ and $y$

For fixed $x$ numerical value of ${\sigma}_{22}$ gradually increases as $y$ increases. For fixed $y$ numerical value of ${\sigma}_{22}$ gradually increases as $x$ increases.

Fig. 7 represent distribution of normal stress ${\sigma}_{12}$ for different values of $x$ and $t$ for fixed $y=$0.2 and 1.

For fixed $x$, nominal decreasing of numerical values of ${\sigma}_{12}$ has been seen as $t$ increases, while For fixed $t$, numerical values of ${\sigma}_{12}$ decreases gradually as $x$ increases. numerical values of ${\sigma}_{12}$ minimum at $x=$1 and 0.02$\le t\le $0.1

Fig. 8 represent distribution of normal stress ${\sigma}_{22}$ for different values of $x$ and $t$ for fixed $y=$0.5 and 2.

For fixed $x$, nominal decreasing of ${\sigma}_{22}$ has been seen as t increases. For fixed $t$, numerical values of ${\sigma}_{12}$ decreases as $x$ increases. Also, significant changes occur in the region 0.6$\le x\le $1.0 and 0$\le t\le $1.0.

**Fig. 7.**
Stress component ${\sigma}_{12}$ at $y=$0.2 and $\omega =$1 verses $x$ and $t$

**Fig. 8.**
The distribution of stress component ${\sigma}_{22}$ at $y=$0.5 and $\omega =$2 verses $x$ and $t$

Fig. 9 represent distribution of temperature, $T$ for different values of $x$ and $y$ for fixed $t=$0.3 and 2.

For fixed $x$ numerical value of $T$ gradually decreases as $y$ increases. For fixed $y$ numerical value of $T$ gradually increases as $x$ increases. $T$ in minimum at $x=$0 and significant changes occurs in the region 0.6$\le x\le $ 1.0 and 0 $\le y\le $1.0

Fig. 10 represent distribution of temperature, $T$ for different values of $x$ and $t$ for fixed $y=$0.1 and 1.5.

**Fig. 9. **
The variation of $T$ at $t=$0.3 and $\omega =$2 verses $x$ and $y$

For fixed $t$ numerical value of $T$ nominally increases as $x$ increases. For fixed $x$ numerical value of $T$ gradually increases as $t$ increases.

Fig. 11 represent distribution of temperature, $T$ for different values of $y$ and $t$ for fixed $x=$0.5 and 3.

For fixed $t$ numerical value of $T$ nominally increases as $y$ increases. For fixed $x$ numerical values of $T$ decreases as $t$ increases.

**Fig. 10.**
Variation of $T$ at $y=$0.1 and $\omega =$1.5 verses $x$ and $t$

**Fig. 11. **
Variation of $T$ at $x=$0.5 and $\omega =$3 verses $y$ and $t$

Fig. 12 represents distribution of normal stress ${\sigma}_{11}$ for $y=$0.3, $t=$0.01 and 5 for different numerical values of ${R}_{p}$.

For fixed ${R}_{p}$, ${\sigma}_{11}$ gradually decreases as $x$ increases. For fixed x numerical values of ${\sigma}_{11}$ gradually decreases as ${R}_{p}$ increases.

Fig. 13 represents distribution of normal stress ${\sigma}_{11}$ for $y=$0.3, $t=$0.01 and 5 for different fractional values of ${R}_{p}$.

For fixed time ${R}_{p}$, ${\sigma}_{11}$ gradually increases as $x$ increases. For fixed $x$ numerical values of ${\sigma}_{11}$ gradually decreases as ${R}_{p}$ increases. Significant changes occurred for 0$\le x\le $ 0.4.

Fig. 14 represents the distribution of stress component ${\sigma}_{12}$ at $y=$0.2, $t=$0.03 and 0.3 for different fractional values ${R}_{p}$ of verses $x$ for $P1=$1.

For fixed time ${R}_{p}$, ${\sigma}_{11}$ gradually decreases as $x$ increases. For fixed x numerical values of ${\sigma}_{11}$ increases as ${R}_{p}$ increases. Significant changes occurred for 0$\le x\le $ 0.4.

Fig. 15 represents the distribution of stress component ${\sigma}_{12}$ at $y=$0.3, $t=$0.1 and 4 for different integral values ${R}_{p}$ for $P1=$0.

For fixed $Rp$, ${\sigma}_{12}$ gradually decreases as $x$ increases. For fixed $x$ numerical values of ${\sigma}_{12}$ increases as ${R}_{p}$ increases.

Fig. 16 represents the distribution of stress component ${\sigma}_{22}$ at $y=$0.2, $t=$0.6 and 0.4 for different values ${R}_{p}$ for $P1=$0.

For fixed ${R}_{p}=$0.9, ${\sigma}_{22}$ gradually increases as $x$ increases but for other fixed values of ${R}_{p}$, ${\sigma}_{22}$ gradually decreases as $x$ increases. For fixed $x$, ${\sigma}_{22}$ gradually increases as $x$ increases for ${R}_{p}=$0.9 but for other fixed values of ${R}_{p}$, ${\sigma}_{22}$ gradually decreases as $x$ increases.

Fig. 17 represents the distribution of stress component ${\sigma}_{11}$ for different values of $t$ for fixed $y=$0.2 and 1.5.

**Fig. 12.**
Stress component ${\sigma}_{11}$ at $y=$0.3, $t=$0.01 and $\omega =$ 5 for different values ${R}_{p}$ of verses $x$

**Fig. 13.**
Stress component ${\sigma}_{11}$ at $y=$0.2 and $\omega =$ 0.2 for different values ${R}_{p}$ of verses $x$

**Fig. 14.**
Stress component ${\sigma}_{12}$ at $y=$0.2, $t=$0.03 and $\omega =$ 0.3 for different values ${R}_{p}$ of verses $x$ for $P1=$1

**Fig. 15.**
Stress component ${\sigma}_{12}$ at $y=$0.3, $t=$0.1 and $\omega =$ 4 for different values ${R}_{p}$ of verses $x$ for $P1=$0

For fixed $t$, the numerical value of ${\sigma}_{11}$ gradually increases as $x$ increases. For fixed $x$, the numerical value of ${\sigma}_{11}$ gradually increases as $t$ increases.

Fig. 18 The distribution of stress component ${\sigma}_{12}$ for fixed $y=$0.3 and 0.5 for different values of $t$.

For fixed $t$, the numerical value of ${\sigma}_{12}$ gradually decreases as $x$ increases. For fixed $x$, the numerical value of ${\sigma}_{12}$ gradually increases as $t$ increases.

**Fig. 16.**
Stress component ${\sigma}_{22}$ at $y=$0.2, $t=$0.6 and $\omega =$ 0.4 for different values ${R}_{p}$ of verses $x$ for $P1=$0

**Fig. 17. **
Stress component ${\sigma}_{11}$ at $y=$0.2 and $\omega =$ 1.5 for different values of $t$ verses $x$

**Fig. 18.**
Stress component ${\sigma}_{12}$ at $y=$0.3 and $\omega =$ 0.5 for different values of $t$ verses $x$

**Fig. 19.**
Stress component ${\sigma}_{22}$ at $y=$0.5 and $\omega =$ 0.2 for different values of $t$ verses $x$

Fig. 19 represents the distribution of stress component ${\sigma}_{22}$ at fixed $y=$0.5 and 0.2 for different values of $t$.

For fixed $t$, the numerical value of ${\sigma}_{12}$ gradually increases as $x$ increases, but for fixed $x$, the numerical value of ${\sigma}_{12}$ gradually decreases as t increases.

Fig. 20 represents distribution of normal stress ${\sigma}_{11}$ for different values of $x$ and for fixed $t=$0.3 and $y=$0.3.

For fixed $x$ numerical value of ${\sigma}_{11}$ remain constant as increases, but for fixed $\omega $ the numerical value of ${\sigma}_{11}$ gradually decreases as $x$ increases.

Fig. 21 represents distribution of normal stress ${\sigma}_{12}$ for different values of $x$ and for fixed $y=$0.3 and $t=$0.5.

For fixed $x$, nominal increasing of numerical values of ${\sigma}_{12}$ has been seen as increases, while for fixed $\omega $, numerical values of ${\sigma}_{12}$ increases gradually as $x$ increases after $x=$0.6 (approx.). Numerical values of ${\sigma}_{12}$ minimum at $x=$0.5 (approx.) and 1$\le t\le $3.

Fig. 22 represents distribution of stress component ${\sigma}_{22}$ at for different values of $x$ and for fixed $t=$0.1 and $y=$0.4.

For fixed $\omega $ numerical value of ${\sigma}_{22}$ gradually decreases as $x$ increases. Numerical values of ${\sigma}_{22}$ minimum at $x=$1 and 0$\le \omega \le $1.

**Fig. 20.**
Stress component ${\sigma}_{11}$ at $y=$0.3 and $t=$0.3 verses $x$ and $\omega $

**Fig. 21.**
Stress component ${\sigma}_{12}$ at $y=$0.3 and $t=$0.5 verses $x$ and $\omega $

**Fig. 22.**
Stress component ${\sigma}_{22}$ at $y=$0.4 and $t=$0.1 verses $x$ and $\omega $

Fig. 23 represents distribution of temperature, $T$ for different values of $x$ and $\omega $ for fixed $t=$ 0.1 and $y=$ 0.4.

For fixed $x$ numerical value of $T$ gradually increases in the region 0$\le \omega \le $0.3 (approx.) For fixed $\omega $ numerical value of $T$ nominally increases as $x$ increases.

**Fig. 23.**
The variation of $T$ at $y=$0.4 and $t=$0.1 verses $x$ and $\omega $

We consider the physical parameters in SI units given in Dhaliwal and Singh following below to obtain the numerical result to observe the effect of wave propagation:

$\rho =2660kg/{m}^{3},\lambda =5.65\times 1{0}^{10}N/{m}^{2},$

${\mu}_{T}=2.46\times 1{0}^{10}N/{m}^{2},{\mu}_{L}=5.66\times 1{0}^{10}N/{m}^{2},$

$\alpha =-1.28\times 1{0}^{10}N/{m}^{2},\beta =220.90\times 1{0}^{10}N/{m}^{2},$

${\alpha}_{11}=0.017\times 1{0}^{-4}de{g}^{-1},{\alpha}_{22}=0.015\times 1{0}^{-4}de{g}^{-1},$

$l=0.5,{T}_{0}=293K,{c}_{e}=0.787\times 1{0}^{3}JK{g}^{-1}de{g}^{-1}$

${K}_{11}=0.0921\times 1{0}^{10}J{m}^{-1}{s}^{-1}de{g}^{-1}$

${K}_{22}=0.0963\times 1{0}^{10}J{m}^{-1}{s}^{-1}de{g}^{-1}$

${P}_{1}=0or1,{P}_{2}=0.1,{P}_{3}=0.2.$

${\mu}_{T}=2.46\times 1{0}^{10}N/{m}^{2},{\mu}_{L}=5.66\times 1{0}^{10}N/{m}^{2},$

$\alpha =-1.28\times 1{0}^{10}N/{m}^{2},\beta =220.90\times 1{0}^{10}N/{m}^{2},$

${\alpha}_{11}=0.017\times 1{0}^{-4}de{g}^{-1},{\alpha}_{22}=0.015\times 1{0}^{-4}de{g}^{-1},$

$l=0.5,{T}_{0}=293K,{c}_{e}=0.787\times 1{0}^{3}JK{g}^{-1}de{g}^{-1}$

${K}_{11}=0.0921\times 1{0}^{10}J{m}^{-1}{s}^{-1}de{g}^{-1}$

${K}_{22}=0.0963\times 1{0}^{10}J{m}^{-1}{s}^{-1}de{g}^{-1}$

${P}_{1}=0or1,{P}_{2}=0.1,{P}_{3}=0.2.$

**Green W. A.**Bending waves in strongly anisotropic plates. Quarterly Journal of Mechanics and Applied Mathematics, Vol. 35, 1982, p. 485-507.**Abbas I. A., Othman M. I. A.**Generalized thermoelastic interaction in a fiber-reinforced anisotropic half-space under hydrostatic initial stress. Journal of Vibration Control, Vol. 18, Issue 2, 2011, p. 175-182.**Baylis E. R., Green W. A.**Flexural waves in fiber reinforced laminated plates. Journal of Sound and Vibration, Vol. 110, 1986, p. 1-26.**Rogerson G. A.**Penetration of impact waves in a six-ply fiber composite laminate. Journal of Sound and Vibration, Vol. 158, 1992, p. 105-120.**Dhaliwal R. S., Sherief H. H.**Generalized thermoelasticity for anisotropic media. Quarterly of Applied Mathematics, Vol. 33, 1980, p. 1-8.**Sherief H. H., El-Sayed A., El-Latief A.**Fractional order theory of thermoelasticity. International Journal of Solids and Structures, Vol. 47, 2010, p. 269-275.**Santra S., Das N. C., Kumar R., Lahir A.**Three dimensional fractional order generalized thermoelastic problem under effect of rotation in a half space. Journal of Thermal Stresses, Vol. 38, 2015, p. 309-324.**Bachher M., Sarkar N., Lahiri A.**Generalized thermoelastic infinite medium with voids subjected to a instantaneous heat sources with fractional derivative heat transfer. International Journal of Mechanical Sciences, Vol. 89, 2014, p. 84-91.**Lord H. W., Shulman Y.**A generalized dynamical theory of thermoelasticity. Journal of the Mechanics and Physics of Solids, Vol. 15, 1967, p. 299-309.**Green A. E., Lindsay K. A.**Thermoelasticity. Journal of Elasticity, Vol. 2, 1972, p. 1-7.**Green A. E., Naghdi P. M.**A re-examination of the basic results of thermomechanics. Proceedings of the Royal Society of London A, Vol. 432, 1991, p. 171-194.**Green A. E., Naghdi P. M.**On undamped heat waves in an elastic solid. Thermal Stresses, Vol. 15, 1992, p. 252-264.**Green A. E., Naghdi P. M.**Thermoelasticity without energy dissipation. Elasticity, Vol. 31, 1993, p. 189-208.**Ibrahim A. Abbas**Generalized magnetothermoelasticity in a fiber-reinforced anisotropic half-space. International Journal of Thermophysics, 2011, https://doi.org/10.1007/s10765-011-0957-3.