Axisymmetric deformation in transversely isotropic magneto-thermoelastic solid with Green–Naghdi III due to inclined load
Iqbal Kaur ORCID: orcid.org/0000-0002-2210-77011 &
Parveen Lata1
The axisymmetric problem in two-dimensional transversely isotropic magneto-thermoelastic (TIMT) solid due to inclined load with Green–Naghdi (GN)-III theory and two temperature (2T) has been studied. The Laplace and Hankel transform has been used to get the expressions of temperature distribution, displacement, and stress components with the horizontal distance in the physical domain. The effect of Green–Naghdi theories of type I, II, and III theories of thermoelasticity has been studied graphically on the resulting quantities. A special case for the magneto-thermoelastic isotropic medium has also been studied.
The study of deformation in a thermoelastic medium is one of the wide and dynamic domains of continuum dynamics. It is well known that all the rotating large bodies have angular velocity, as well as magnetism; therefore, the thermoelastic interactions in a rotating medium under magnetic field are of importance. The study of thermoelasticity is beneficial to analyze the deformation field such as geothermal engineering, advanced aircraft structure design, thermal power plants, composite engineering, geology, high-energy particle accelerators, and many developing technologies.
Eubanks and Sternberg (1954) discussed the axisymmetric issue of elasticity concept for a transversely isotropy medium. Vendhan and Archer (1978) electrostatically analyzed transversely isotropic (TI) finite stress-free cylinders with lateral surfaces using displacement potential. Green and Naghdi (1992, 1993) dealt with the linear and the nonlinear theories of the thermoelastic body with and without energy dissipation. Three new thermoelastic theories were proposed by them, based on entropy equality. Their theories are known as GN-I, GN-II, and GN-III theories of thermoelasticity. On linearization, type I becomes the classical heat equation, whereas on linearization, type II as well as type III theories predicts the finite speed of thermal wave. Savruk (1994) discussed the axisymmetric deformation of a TI body containing cracks.
Tarn et al. (2009) analyzed the axisymmetric and stress dispersal in a TI roundabout barrel-shaped body utilizing Hamiltonian variational definition through Legendre's change. Liang and Wu (2012) discussed the axisymmetric deformation of one TI cylinder with the Lure method. Mahmoud (2012) considered the impact of relaxation times, the rotation, and the initial stress on Rayleigh waves. Shi et al. (2016) presented the thermomagnetoelectroelastic field in a heterogeneous annular multi-ferric composite plate with thermal loadings which is uniformly distributed on the boundaries. Kumar et al. (2016a, 2016b) studied the conflicts of thermomechanical sources in a TI homogeneous thermoelastic rotating medium with magnetic effect as well as two temperature applied to the thermoelasticity GN-III theories. Li et al. (2016) presented a set of axisymmetric solutions of the thermoelastic field in a heterogeneous circular plate which could either simply reinforced or fastened persuaded by the external thermal load.
Including these, various researchers dealt with various theory of thermoelasticity such as Marin (Marin 1997a, 1997b; Marin 1998; Marin 1999), Marin (Marin 2008; Marin 1997a, 1997b), Ezzat et al. (2012), Atwa (2014), Marin et al. (2013), Marin (2016), Marin and Baleanu (2016), Bijarnia and Singh (2016), Sharma et al. (Sharma et al. 2015a, 2015b; 2016; 2017), Ezzat et al. (2017), Lata (2018), Lata et al. (2016), Marin and Öchsner (2017), Othman and Marin (2017), Ezzat et al. (2017), Chauthale and Khobragade (2017), Kumar et al. (2017), Lata and Kaur (2018), Shahani and Torki (2018), Lata and Kaur (2019a, 2019b, 2019c, 2019d), Kaur and Lata (2019a, 2019b), Bhatti and Lu (2019a, 2019b), and Marin et al. (2019).
Despite these, very less work has been done in thermomechanical interactions in TIMT rotating solid with GN-III theory, with two temperature in the axisymmetric medium. Remembering these contemplations, analytic expressions for the displacement components, stress components, and temperature distribution in two-dimensional homogeneous, TIMT rotating solids with GN-III theories, with two temperature have been derived.
Basic equations
The field equations with and without energy dissipation, without body forces and heat sources for an anisotropic thermoelastic medium following Lata and Kaur (2019d), are:
$$ {t}_{ij}={C}_{ij kl}{e}_{kl}-{\beta}_{ij}T, $$
$$ {K}_{ij}{\varphi}_{, ij}+{K}_{ij}^{\ast }{\overset{.}{\varphi}}_{, ij}={\beta}_{ij}{T}_0{\ddot{\mathrm{e}}}_{ij}+\rho {C}_E\ddot{T}. $$
and the equation of motion for a medium rotating uniformly and Lorentz force is
$$ {t}_{ij,j}+{F}_i=\rho \left\{{\ddot{u}}_i+\right(\varOmega \times {\left(\varOmega \times u\right)}_i+{\left(2\varOmega \times \overset{.}{u}\right)}_i\Big\}, $$
Ω = Ωn, n is a unit vector signifying the direction of the rotating axis.
$$ {F}_i={\mu}_0\left(\overrightarrow{j}\times {\overrightarrow{H}}_0\right), $$
$$ T=\varphi -{a}_{ij}{\varphi}_{, ij}, $$
$$ {\beta}_{ij}={C}_{ij kl}{\alpha}_{ij}, $$
$$ {e}_{ij}=\frac{1}{2}\left({u}_{i,j}+{u}_{j,i}\right).i=1,2,3 $$
$$ {\beta}_{ij}={\beta}_i{\delta}_{ij},{K}_{ij}={K}_i{\delta}_{ij},{K}_{ij}^{\ast }={K}_i^{\ast }{\delta}_{ij},i\; is\kern0.17em not\kern0.17em summed $$
$$ Here\;{C}_{ijkl}\; having\kern0.17em symmetry\;\left({C}_{ijkl}={C}_{klij}={C}_{jikl}={C}_{ijlk}\right). $$
Method and formulation of the problem
Consider a TIMT homogeneous medium with an initial temperature T0. We consider a cylindrical polar coordinate system (r, θ, z) with symmetry about the z-axis. For a plane axisymmetric problem, v = 0, and u, w, and φ are independent of θ. Additionally, we take
$$ \varOmega =\left(0,\varOmega, 0\right). $$
$$ {J}_2=0. $$
The density components J1and J3 are given as
$$ {J}_1=-{\varepsilon}_0{\mu}_0{H}_0\frac{\partial^2w}{\partial {t}^2}, $$
$$ {J}_3={\varepsilon}_0{\mu}_0{H}_0\frac{\partial^2u}{\partial {t}^2}. $$
Using the appropriate transformation following Slaughter (2002) on Eqs. (1)–(3) to determine the conditions for TI thermoelastic solid with 2T and with and without energy dissipation, we get
$$ {C}_{11}\left(\frac{\partial^2u}{\partial {r}^2}+\frac{1}{r}\frac{\partial u}{\partial r}-\frac{1}{r^2}u\right)+{C}_{13}\left(\frac{\partial^2w}{\partial r\partial z}\right)+{C}_{44}\frac{\partial^2u}{\partial {z}^2}+{C}_{44}\left(\frac{\partial^2w}{\partial r\partial z}\right)-{\beta}_1\frac{\partial }{\partial r}\left\{\varphi -{a}_1\left(\frac{\partial^2\varphi }{\partial {r}^2}+\frac{1}{r}\frac{\partial \varphi }{\partial r}\right)-{a}_3\frac{\partial^2\varphi }{\partial {z}^2}\right\}-{\mu}_0{J}_3{H}_0=\rho \left(\frac{\partial^2u}{\partial {t}^2}-{\varOmega}^2u+2\varOmega \frac{\partial w}{\partial t}\right), $$
$$ \left({C}_{11}+{C}_{44}\right)\left(\frac{\partial^2u}{\partial r\partial z}+\frac{1}{r}\frac{\partial u}{\partial z}\right)+{C}_{44}\left(\frac{\partial^2w}{\partial {r}^2}+\frac{1}{r}\frac{\partial w}{\partial r}\right)+{C}_{33}\frac{\partial^2w}{\partial {z}^2}-{\beta}_3\frac{\partial }{\partial z}\left\{\varphi -{a}_1\left(\frac{\partial^2\varphi }{\partial {r}^2}+\frac{1}{r}\frac{\partial \varphi }{\partial r}\right)-{a}_3\frac{\partial^2\varphi }{\partial {z}^2}\right\}+{\mu}_0{J}_1{H}_0=\rho \left(\frac{\partial^2w}{\partial {t}^2}-{\varOmega}^2w-2\varOmega \frac{\partial u}{\partial t}\right), $$
$$ \left({K}_1+{K}_1^{\ast}\frac{\partial }{\partial t}\right)\left(\frac{\partial^2\varphi }{\partial {r}^2}+\frac{1}{r}\frac{\partial \varphi }{\partial r}\right)+\left({K}_3+{K}_3^{\ast}\frac{\partial }{\partial t}\right)\frac{\partial^2\varphi }{\partial {z}^2}={T}_0\frac{\partial^2}{\partial {t}^2}\left({\beta}_1\frac{\partial u}{\partial r}+{\beta}_3\frac{\partial w}{\partial z}\right)+\rho {C}_E\frac{\partial^2}{\partial {t}^2}\left\{\varphi -{a}_1\left(\frac{\partial^2\varphi }{\partial {r}^2}+\frac{1}{r}\frac{\partial \varphi }{\partial r}\right)-{a}_3\frac{\partial^2\varphi }{\partial {z}^2}\right\}. $$
Constitutive relations are
$$ {\displaystyle \begin{array}{c}{t}_{rr}={c}_{11}{e}_{rr}+{c}_{12}{e}_{\theta \theta}+{c}_{13}{e}_{zz}-{\beta}_1T,\\ {}{t}_{zr}=2{c}_{44}{e}_{rz},\\ {}{t}_{zz}={c}_{13}{e}_{rr}+{c}_{13}{e}_{\theta \theta}+{c}_{33}{e}_{zz}-{\beta}_3T,\\ {}{t}_{\theta \theta}={c}_{12}{e}_{rr}+{c}_{11}{e}_{\theta \theta}+{c}_{13}{e}_{zz}-{\beta}_3T,\end{array}} $$
$$ {\displaystyle \begin{array}{c}{e}_{rz}=\frac{1}{2}\left(\frac{\partial u}{\partial z}+\frac{\partial w}{\partial r}\right),\\ {}{e}_{rr}=\frac{\partial u}{\partial r},\\ {}{e}_{\theta \theta}=\frac{u}{r},\\ {}{e}_{zz}=\frac{\partial w}{\partial z},\\ {}T=\varphi -{a}_1\left(\frac{\partial^2\varphi }{\partial {r}^2}+\frac{1}{r}\frac{\partial \varphi }{\partial r}\right)-{a}_3\frac{\partial^2\varphi }{\partial {z}^2},\\ {}{\beta}_1=\left({c}_{11}+{c}_{12}\right){\alpha}_1+{c}_{13}{\alpha}_3,\\ {}{\beta}_3=2{c}_{13}{\alpha}_1+{c}_{33}{\alpha}_3.\end{array}} $$
We consider that a primary medium is at rest. Therefore, the preliminary and symmetry conditions are assumed as
$$ {\displaystyle \begin{array}{c}u\left(r,z,0\right)=0=\overset{.}{u}\left(r,z,0\right),\\ {}\ w\left(r,z,0\right)=0=\overset{.}{w}\left(r,z,0\right),\\ {}\varphi \left(r,z,0\right)=0=\overset{.}{\varphi}\left(r,z,0\right)\ \mathrm{for}\ z\ge 0,-\infty <r<\infty, \\ {}u\left(r,z,t\right)=w\left(r,z,t\right)=\varphi \left(r,z,t\right)=0\ \mathrm{for}\ t>0\ \mathrm{when}\ z\to \infty .\end{array}} $$
To simplify the solution, the following dimensionless quantities are introduced
$$ {r}^{\hbox{'}}=\frac{r}{L},{z}^{\hbox{'}}=\frac{z}{L},{t}^{\hbox{'}}=\frac{c_1}{L}t,{u}^{\hbox{'}}=\frac{\rho {c}_1^2}{L{\beta}_1{T}_0}u,{w}^{\hbox{'}}=\frac{\rho {c}_1^2}{L{\beta}_1{T}_0}w,{T}^{\hbox{'}}=\frac{T}{T_0},{t}_{zr}^{\hbox{'}}=\frac{t_{zr}}{\beta_1{T}_0},{t}_{zz}^{\hbox{'}}=\frac{t_{zz}}{\beta_1{T}_0},{\varphi}^{\hbox{'}}=\frac{\varphi }{T_0},{a}_1^{\hbox{'}}=\frac{a_1}{L^2},{a}_3^{\hbox{'}}=\frac{a_3}{L^2},{h}^{\hbox{'}}=\frac{h}{H_0},{\varOmega}^{\hbox{'}}=\frac{L}{C_1}\Omega . $$
Applying the dimensionless quantities introduced in (13) on Eqs. (9)–(11) and subsequently suppressing the primes and using the following Laplace and Hankel transforms
$$ \overline{f}\left(r,z,s\right)=\underset{0}{\overset{\infty }{\int }}f\left(r,z,t\right){e}^{- st} dt, $$
$$ \tilde{f}\left(\xi, z,s\right)={\int}_0^{\infty}\overline{f}\left(r,z,s\right)r{J}_n\left( r\xi \right) dr. $$
on the resulting quantities, we obtain
$$ \left(-{\xi}^2+{\delta}_2{D}^2-{s}^2{\delta}_7+{\varOmega}^2\right)\tilde{u}+\left({\delta}_1 D\xi -2\varOmega \mathrm{s}\right)\tilde{w}+\left(-\xi \left(1+{a}_1{\xi}^2\right)+{a}_3\xi {D}^2\right)\tilde{\varphi}=0, $$
$$ \left({\delta}_1 D\xi +2\varOmega s\right)\tilde{u}+\left({\delta}_3{D}^2-{\delta}_2{\xi}^2-{s}^2{\delta}_7+{\varOmega}^2\right)\tilde{w}-\left(\frac{\beta_3}{\beta_1}D\left[\left(1+{\xi}^2{a}_1\right)-{a}_3{D}^2\right]\right)\tilde{\varphi}=0, $$
$$ -{\delta}_6{s}^2\xi \tilde{u}-\frac{\beta_3}{\beta_1}{\delta}_6{s}^2D\tilde{w}+\left(-{\delta}_8{s}^2\left(1+{\xi}^2{a}_1-{a}_3{D}^2\right)-{\xi}^2\left({K}_1+{\delta}_4s\right)+{D}^2\left({K}_3+{\delta}_5s\right)\right)\tilde{\varphi}=0. $$
$$ {\delta}_1=\frac{c_{13}+{c}_{44}}{c_{11}},{\delta}_2=\frac{c_{44}}{c_{11}},{\delta}_3=\frac{c_{33}}{c_{11}},{\delta}_4=\frac{K_1^{\ast }{C}_1}{L},{\delta}_5=\frac{K_3^{\ast }{C}_1}{L},{\delta}_6=-\frac{T_0{\beta}_1^2}{\rho },{\delta}_7=\frac{\varepsilon_0{\mu}_0^2{H}_0^2}{\rho }+1,{\delta}_8=-\rho {C}_E{C}_1^2. $$
The non-trivial solution of (16)–(18) by eliminating\( \overset{\sim }{\ u} \), \( \overset{\sim }{w} \), and \( \overset{\sim }{\varphi } \) yields
$$ A{D}^6+B{D}^4+C{D}^2+E=0 $$
$$ {\displaystyle \begin{array}{c}A={\delta}_2{\delta}_3{\zeta}_{11}-{\zeta}_9{\updelta}_2{\zeta}_7,\\ {}B={\delta}_2{\zeta}_3{\zeta}_{11}+{\delta}_3{\zeta}_1{\zeta}_{11}+{\delta}_2{\delta}_3{\zeta}_{10}-{\delta}_2{\zeta}_7{\zeta}_{10}-{\zeta}_7{\zeta}_1{\zeta}_9-{\zeta}_2^2{\zeta}_{11}+{\zeta}_2{\zeta}_4{\zeta}_9+{\zeta}_8{\zeta}_7{\zeta}_2-{\delta}_3{\zeta}_4{\zeta}_8,\\ {}C={\zeta}_1{\zeta}_3{\zeta}_{11}+{\delta}_2{\zeta}_{10}{\zeta}_5+{\delta}_3{\zeta}_3{\zeta}_{10}-{\zeta}_9{\zeta}_1{\zeta}_6-{\zeta}_2^2{\zeta}_{10}+{\zeta}_2{\zeta}_6{\zeta}_8+{\zeta}_3{\zeta}_2{\zeta}_9-{\zeta}_3{\zeta}_8{\delta}_3-{\zeta}_4{\zeta}_8{\zeta}_5+4{\varOmega}^2{s}^2{\zeta}_{11},\\ {}E={\zeta}_5{\zeta}_1{\zeta}_{10}-{\zeta}_8{\zeta}_5{\zeta}_3+4{\varOmega}^2{s}^2{\zeta}_{10}.\\ {}{\zeta}_1=-{\xi}^2-{s}^2{\delta}_7+{\varOmega}^2,\\ {}{\zeta}_2={\delta}_1\xi, \\ {}{\zeta}_3=-\xi \left({a}_1{\xi}^2+1\right),\\ {}{\zeta}_4={a}_3\xi, \\ {}{\zeta}_5=-{\delta}_2{\xi}^2-{s}^2{\delta}_7+{\varOmega}^2,\\ {}{\zeta}_6=-\frac{\beta_3}{\beta_1}\left(1+{a}_1{\xi}^2\right),\\ {}{\zeta}_7={a}_3\frac{\beta_3}{\beta_1},\\ {}{\zeta}_8=-\xi {\delta}_6{s}^2,\\ {}{\zeta}_9=-\frac{\beta_3}{\beta_1}{\delta}_6{s}^2,\\ {}{\zeta}_{10}=-{\delta}_8{s}^2\left(1+{a}_1{\xi}^2\right)-{\xi}^2\left({K}_1+{\delta}_4s\right),\\ {}{\zeta}_{11}=\left({K}_3+{\delta}_5s\right)+{\delta}_8{s}^2{\mathrm{a}}_3.\end{array}} $$
The solutions of Eq. (19) can be written as
$$ \hat{u}\left(\xi, z,s\right)={\sum}_{j=1}^3{A}_j{e}^{-{\lambda}_jz}, $$
$$ \hat{w}\left(\xi, z,s\right)={\sum}_{j=1}^3{d}_j{A}_j{e}^{-{\lambda}_jz}, $$
$$ \hat{\varphi}\left(\xi, z,s\right)=\sum \limits_{j=1}^3{l}_j{A}_j{e}^{-{\lambda}_jz}, $$
where Aj being arbitrary constants, ±λj represents the roots of Eq. (19) and dj and lj are given by
$$ {d}_j=\frac{{\delta}_{2}{\zeta}_{11}{\lambda}^{4}_{j}+(\zeta_{11}\zeta_{1}-\zeta_{4}\zeta_{8}+\delta_{2}\zeta_{10})\lambda^{2}_{j}+{\zeta}_{1}{\zeta}_{10}-{\zeta}_{8}{\zeta}_{3}} {(\delta_{3}\zeta_{11}-\zeta_{7}\zeta_{9})\lambda^{4}_{j}+(\delta_{3}\zeta_{10}+\zeta_{5}\zeta_{11}-\zeta_{9}\zeta_{6})\lambda^{2}_{j}+ \zeta_{5}\zeta_{10}}, $$
$$ {l}_j=\frac{{\delta}_{2}{\delta}_{3}{\lambda}^{4}_{j}+(\delta_{2}\zeta_{5}+\zeta_{1}\delta_{3}-\zeta^{2}_{2})\lambda^{2}_{j}+{\zeta}_{1}{\zeta}_{5}+4\Omega^{2}s^{2}} {(\delta_{3}\zeta_{11}-\zeta_{7}\zeta_{9})\lambda^{4}_{j}+(\delta_{3}\zeta_{10}+\zeta_{5}\zeta_{11}-\zeta_{9}\zeta_{6})\lambda^{2}_{j}+ \zeta_{5}\zeta_{10}}, $$
$$ \tilde{t_{zz}}=\sum {A}_j\left(\xi, s\right){\eta}_j{e}^{-{\lambda}_jz}, $$
$$ \tilde{t_{rz}}=\sum {A}_j\left(\xi, s\right){\mu}_j{e}^{-{\lambda}_jz}, $$
$$ \tilde{t_{rr}}=\sum {A}_j\left(\xi, s\right){Q}_j{e}^{-{\lambda}_jz}, $$
$$ {\eta}_j={\delta}_{11}\xi -{\delta}_3{\lambda}_j{d}_j-\frac{\beta_3}{\beta_1}\left(1+{a}_1{\xi}^2\right){\mathrm{l}}_{\mathrm{j}}+\frac{\beta_3}{\beta_1}{\mathrm{a}}_3{l}_j{\lambda}_j^2, $$
$$ {\mu}_j={\delta}_2\left(-{\lambda}_j+\xi {d}_j\right), $$
$$ {Q}_j=\left({\delta}_{10}+1\right)\xi -{\delta}_{11}{\lambda}_j{d}_j-{l}_j\left(1+{a}_1{\xi}^2\right)+{\mathrm{a}}_3{l}_j{\lambda}_{\mathrm{j}}^2,i,j=1,2,3. $$
The boundary conditions when normal force and tangential load are applied to the half-space (z = 0) are
$$ {t}_{zz}\left(r,z,t\right)=-{F}_1{\psi}_1(r)H(t), $$
$$ {t}_{rz}\left(r,z,t\right)=-{F}_2{\psi}_2(r)H(t), $$
$$ \frac{\partial \varphi \left(r,z,t\right)}{\partial z}+ h\varphi \left(r,z,t\right)=0. $$
where ψ1(r) and ψ2(r) are the vertical and the tangential load applied on along the r-axis and
$$ H(t)=\Big\{{\displaystyle \begin{array}{c}1t>0\\ {}0t<0\end{array}} $$
$$ {\displaystyle \begin{array}{c}\sum {A}_j\left(\xi, s\right){\eta}_j=-{F}_1{\psi}_1\left(\xi \right),\\ {}\sum {A}_j\left(\xi, s\right){\mu}_j=-{F}_2{\psi}_2\left(\xi \right),\\ {}\sum {A}_j\left(\xi, s\right){P}_j=0.\ \mathrm{where},{P}_j={l}_j\left(-{\lambda}_j+h\right).\end{array}} $$
Solving Eqs. (27)–(29) with the aid of (20)–(25), we obtain
$$ \tilde{u}=\frac{F_1\tilde{\psi_1}\left(\xi \right)}{\varLambda}\left[\sum \limits_{j=1}^3{\varLambda}_{1j}{\theta}_j\right]+\frac{F_2\tilde{\psi_2}\left(\xi \right)}{\varLambda}\left[\sum \limits_{j=1}^3{\varLambda}_{2j}{\theta}_j\right], $$
$$ \tilde{w}=\frac{F_1\tilde{\psi_1}\left(\xi \right)}{\varLambda}\left[\sum \limits_{j=1}^3{d}_j{\varLambda}_{1j}{\theta}_j\right]+\frac{F_2\tilde{\psi_2}\left(\xi \right)}{\varLambda}\left[\sum \limits_{j=1}^3{d}_j{\varLambda}_{2j}{\theta}_j\right], $$
$$ \tilde{\varphi}=\frac{F_1\tilde{\psi_1}\left(\xi \right)}{\varLambda}\left[\sum \limits_{j=1}^3{l}_j{\varLambda}_{1j}{\theta}_j\right]+\frac{F_2\tilde{\psi_2}\left(\xi \right)}{\varLambda}\left[\sum \limits_{i=1}^3{l}_j{\varLambda}_{2j}{\theta}_j\right], $$
$$ \tilde{t_{rr}}=\frac{F_1\tilde{\psi_1}\left(\xi \right)}{\varLambda}\left[\sum \limits_{j=1}^3{Q}_j{\varLambda}_{1j}{\theta}_j\right]+\frac{F_2\tilde{\psi_2}\left(\xi \right)}{\varLambda}\left[\sum \limits_{j=1}^3{Q}_j{\varLambda}_{2j}{\theta}_j\right], $$
$$ \tilde{t_{zr}}=\frac{F_1\tilde{\psi_1}\left(\xi \right)}{\varLambda}\left[\sum \limits_{j=1}^3{\mu}_j{\varLambda}_{1j}{\theta}_j\right]+\frac{F_2\tilde{\psi_2}\left(\xi \right)}{\varLambda}\left[\sum \limits_{j=1}^3{\mu}_j{\varLambda}_{2j}{\theta}_j\right], $$
$$ \tilde{t_{zz}}=\frac{F_1\tilde{\psi_1}\left(\xi \right)}{\varLambda}\left[\sum \limits_{j=1}^3{\eta}_j{\varLambda}_{1j}{\theta}_j\right]+\frac{F_2\tilde{\psi_2}\left(\xi \right)}{\varLambda}\left[\sum \limits_{j=1}^3{\eta}_j{\varLambda}_{2j}{\theta}_j\right], $$
$$ {\displaystyle \begin{array}{l}{\varLambda}_{11}=-{\mu}_2{P}_3+{P}_2{\mu}_3,\\ {}{\varLambda}_{12}={\mu}_1{P}_3-{P}_1{\mu}_3,\\ {}{\varLambda}_{13}=-{\mu}_1{P}_2+{P}_1{\mu}_2,\\ {}{\varLambda}_{21}={\eta}_2{P}_3-{P}_2{\eta}_3\\ {}{\varLambda}_{22}=-{\eta}_1{P}_3+{P}_1{\eta}_3\\ {}{\varLambda}_{23}={\eta}_1{P}_2-{P}_1{\eta}_2\\ {}\varLambda =-{\eta}_1{\varLambda}_{11}-{\eta}_2{\varLambda}_{12}-{\eta}_3{\varLambda}_{13},j=1,2,3.\end{array}} $$
Concentrated normal force (CNF)
The CNF applied on the half-space is taken as
$$ {\psi}_1(r)=\frac{\delta (r)}{2\pi r},{\psi}_2(r)=\frac{\delta (r)}{2\pi r}. $$
Applying Hankel transform, we get
$$ {\hat{\psi}}_1\left(\xi \right)=\frac{1}{2\pi \xi},{\hat{\psi}}_2\left(\xi \right)=\frac{1}{2\pi \xi}. $$
The solution of Eqs. (30)–(35) with CNF is obtained using (37).
Uniformly distributed force (UDF)
Let a uniform force F1/constant temperature F2 be applied over a uniform circular region of radius a. We obtained the solution with UDF applied on the half-space by taking
$$ {\psi}_1(r)={\psi}_2(r)=\frac{H\left(a-r\right)}{\pi {a}^2}, $$
where H(a − r) is a Heaviside function. The Hankel transforms of ψ1(r) and ψ2(r)are given by
$$ {\hat{\psi}}_1\left(\xi \right)={\hat{\psi}}_2\left(\xi \right)=\left\{\frac{J_1\left(\xi a\right)}{2\pi a\xi}\right\},\xi \ne 0. $$
The solution of Eqs. (30)–(35) with UDF is obtained using (39).
We considered an inclined load (F0/unit length) applied on a uniform circular region and its inclination with the z-axis is θ (see Fig. 1), we have
$$ {F}_1={F}_0\cos \theta\ \mathrm{and}\ {F}_2={F}_0\sin \theta $$
Geometry of the problem
Using Eq. (40) in Eqs. (30)–(35) and with the aid of Eqs. (37) and (39), we obtain displacement components, stress components, and conductive temperature with uniformly distributed force and concentrated force on the surface of TIMT body with and without energy dissipation.
Particular cases
If we take\( {K}_{ij}^{\ast}\ne 0 \), Eq. (2) is GN-III theory, and thus we obtain the solution of (30)–(35) for TIMT solid with rotation and GN-III theory.
Equation (2) becomes GN-II theory if we take \( {K}_{ij}^{\ast }=0 \), and we obtain the solution of (30)–(35) for TIMT solid with rotation and GN-II theory.
If we take Kij = 0, the equation of GN-III theory reduces to the GN-I theory, which is identical with the classical theory of thermoelasticity, and thus we obtain the solution of (30)–(35) for TIMT solid with rotation and GN-I theory.
If \( {C}_{11}={C}_{33}=\lambda +2\mu, {C}_{12}={C}_{13}=\lambda, {C}_{44}=\mu, \kern0.5em {\alpha}_1={\alpha}_3={\alpha}^{\prime },\kern0.5em {a}_1={a}_3=a,\kern0.5em {b}_1={b}_3=b,{K}_1={K}_3=K,{K}_1^{\ast }={K}_3^{\ast }={K}^{\ast } \), we obtain the solution of (30)–(35) for TIMT materials with rotation and with GN-III theory.
Inversion of the transformation
To obtain the solution to the problem in the physical domain following Sharma et al. (2015a, 2015b), invert the transforms in Eqs. (30)–(35) by inverting the Hankel transform using
$$ {f}^{\ast}\left(r,z,s\right)=\underset{0}{\overset{\infty }{\int }}\xi \overset{\sim }{f}\left(\xi, z,s\right){J}_n\left(\xi r\right) d\xi . $$
The integral in Eq. (41) is calculated using Romberg's integration through adaptive step size as defined in Press et al. (1986).
Numerical results and discussion
For determining the theoretical results and influence of GN-I, GN-II, and GN-III theories of thermoelasticity, the physical data for cobalt material has been considered from Dhaliwal and Singh (1980) as
$$ {\displaystyle \begin{array}{l}{c}_{11}=3.07\times {10}^{11}\ \mathrm{N}\ {\mathrm{m}}^{-2},{c}_{33}=3.581\times {10}^{11}\mathrm{N}\ {\mathrm{m}}^{-2},{c}_{13}=1.027\times {10}^{10}\ \mathrm{N}\ {\mathrm{m}}^{-2},{c}_{44}=1.510\times \\ {}{10}^{11}\ \mathrm{N}\ {\mathrm{m}}^{-2},{\beta}_1=7.04\times {10}^6\mathrm{N}\ {\mathrm{m}}^{-2}\ {\deg}^{-1},{\beta}_3=6.90\times {10}^6\mathrm{N}\ {\mathrm{m}}^{-2}\ {\deg}^{-1},\rho =8.836\times {10}^3\mathrm{kg}\ {\mathrm{m}}^{-3},\\ {}\ {C}_E=4.27\times {10}^2{\mathrm{jkg}}^{-1}{\deg}^{-1},{K}_1=0.690\times {10}^2\mathrm{W}\ {\mathrm{m}}^{-1}\ {\deg}^{-1},{K}_3=0.690\times {10}^2\mathrm{W}\ {\mathrm{m}}^{-1}\ {\mathrm{K}}^{-1},{T}_0=298\ \mathrm{K},\\ {}\ {H}_0=1\ \mathrm{J}\ {\mathrm{m}}^{-1}{\mathrm{nb}}^{-1},\kern0.5em {\varepsilon}_0=8.838\times {10}^{-12}\mathrm{F}\ {\mathrm{m}}^{-1},\kern0.5em L=1,\kern0.5em \mathrm{and}\ \varOmega =0.5.\end{array}} $$
Using these values, the graphical illustrations of displacement components (u and w), conductive temperature φ, normal force stress tzz, tangential stress tzr, and radial stress trr,for a TIMT solid with GN-III theory and with 2T due to inclined load, have been illustrated. The numerical calculations have been obtained by developing a FORTRAN program using the above values for cobalt material.
The black line with a square symbol relates to \( {K}_{ij}^{\ast}\ne 0 \) for TIMT solid with rotation and GN-III theory.
The red line with a circle symbol relates to \( {K}_{ij}^{\ast }=0 \) for TIMT solid with rotation and GN-II theory.
The blue line with a triangle symbol relates to Kij = 0 the GN theory of type I.
Case 1: Concentrated normal force
Figures 2, 3, 4, 5, 6, and 7 illustrate the deviations of u and w, conductive temperature φ, and trr, trz, and tzz for a TIMT medium with concentrated normal force, with rotation, and due to inclined load. From the graph, we find that that the displacement component (u) and conductive temperature φ decreases while stress components ( trr, trz, and tzz) and displacement component (w) show a sharp increase and then an oscillatory pattern with an amplitude difference.
Deviations of u with x (with CNF)
Deviations of w with x (with CNF)
Deviations of φ with x (with CNF)
Deviations of trr with x (with CNF)
Deviations of trz with x (with CNF)
Deviations of tzz with x (with CNF)
Case II: Uniformly distributed force (UDF)
Figures 8, 9, 10, 11, 12, and 13 illustrate the deviations of u and w, conductive temperature φ, and trr, trz, and tzz for a TIMT medium with UDF and with rotation, and due to inclined load. The displacement components (u), stress components ( trr, trz, and tzz), and temperature φ decrease sharply and then show the small oscillatory pattern, while the displacement component (w) first increases during the initial range of distance near the loading surface and follows the small oscillatory pattern for the rest of the values of distance.
Deviations of u with x (with UDF)
Deviations of w with x (with UDF)
Deviations of φ with x (with UDF)
Deviations of trr with x (with UDF)
Deviations of the trzwith x (with UDF)
Deviations of tzzwith x (with UDF)
In the above research, we conclude:
The components of displacement, stress, and temperature distribution for TIMT solid with GN-III theory, with 2T with inclined load, are calculated numerically.
The study motivates to consider magneto-thermoelastic materials as an inventive field of thermoelastic solids. The shape of curves demonstrates the effect of various GN theories and rotation on the body and fulfills the purpose of the study.
The outcomes of this research are extremely helpful in the 2-D problem with dynamic response of inclined load in TIMT medium with rotation which is beneficial to detect the deformation field such as geothermal engineering, advanced aircraft structure design, thermal power plants, composite engineering, geology, high-energy particle accelerators, geophysics, auditory range, and geomagnetism. The proposed model is significant to different problems in thermoelasticity and thermodynamics.
Name of Symbol
SI Unit
δ ij Kronecker delta, ω Frequency Hz
C ijkl Elastic parameters, Nm −2 β ij Thermal elastic coupling tensor, Nm −2 K −1
τ 0 Relaxation Time s Ω Angular Velocity of the Solid s −1
Fi Components of Lorentz force N T Absolute temperature, K
\( {\overrightarrow{H}}_0 \) Magnetic field intensity vector Jm−1nb−1 e ij Strain tensors, Nm −2
φ conductive temperature, Wm −1 K −1 \( \overrightarrow{j} \) Current Density Vector Am −2
t ij Stress tensors, Nm −2 \( \overrightarrow{u} \) Displacement Vector m
μ 0 Magnetic permeability Hm−1 T 0 Reference temperature, K
u i Components of displacement, m ε 0 Electric permeability Fm−1
ρ Medium density, Kgm −3 δ(t) Dirac's delta function
C E Specific heat, JKg −1 K −1 K ij Materialistic constant, Wm −1 K −1
α ij Linear thermal expansion coefficient, K −1 \( {K}_{ij}^{\ast } \) Thermal conductivity, Ns −2 K −1
H(t) Heaviside unit step function F1, F2 magnitude of applied forces N
For the numerical outcomes, the physical data of cobalt material has been considered from Dhaliwal and Singh (1980).
Atwa, S. Y. (2014). Generalized magneto-thermoelasticity with two temperature and initial stress under Green-Naghdi theory. Applied Mathematical Modelling, 38, 21–22.
Bhatti, M. M., & Lu, D. Q. (2019a). An application of Nwogu's Boussinesq model to analyze the head-on collision process between hydroelastic solitary waves. Open Physics, 233(17), 6135.
Bhatti, M. M., & Lu, D. Q. (2019b). Analytical study of the head-on collision process between hydroelastic solitary waves in the presence of a uniform current. Symmetry, 11(3), 333.
Bijarnia, R., & Singh, B. (2016). Propagation of plane waves in a rotating transversely isotropic two temperature generalized thermoelastic solid half-space with voids. International Journal of Applied Mechanics and Engineering, 21(1), 285–301.
Chauthale, S., & Khobragade, N. W. (2017). Thermoelastic Response of a Thick Circular Plate due to Heat Generation and its Thermal Stresses. Global Journal of Pure and Applied Mathematics., 7505-7527.
Dhaliwal, R., & Singh, A. (1980). Dynamic coupled thermoelasticity. New Delhi,India: Hindustan Publication Corporation.
Eubanks, R. A., & Sternberg, E. (1954). On the Axisymmetric Problem of Elasticity Theory for a Medium with Transverse Isotropy. Journal of Rational Mechanics and Analysis, 3, 89–101.
MathSciNet MATH Google Scholar
Ezzat, M. A., El-Karamany, A. S., & El-Bary, A. A. (2017). Two-temperature theory in Green–Naghdi thermoelasticity with fractional phase-lag heat transfer. Microsystem Technologies- Springer Nature, 24(2), 951–961.
Ezzat, M. A., El-Karamany, A. S., & Ezzat, S. M. (2012). Two-temperature theory in magneto-thermoelasticity with fractional order dual-phase-lag heat transfer. Nuclear Engineering and Design (Elsevier), 252, 267–277.
Green, A., & Naghdi, A. P. (1992). On undamped heat waves in an elastic solid. Journal of Thermal Stresses, 15(2), 253–264.
Green, A., & Naghdi, P. (1993). Thermoelasticity without energy dissipation. Journal of Elasticity, 31(3), 189–208.
Kaur, I., & Lata, P. (2019a). Effect of hall current on propagation of plane wave in transversely isotropic thermoelastic medium with two temperature and fractional order heat transfer. SN Applied Sciences, 1, 900.
Kaur, I., & Lata, P. (2019b). Transversely isotropic thermoelastic thin circular plate with constant and periodically varying load and heat source. International Journal of Mechanical and Materials Engineering, 14(10), 1–13.
Kumar, R., Sharma, N., & Lata, P. (2016a). Thermomechanical Interactions Due to Hall Current in Transversely Isotropic Thermoelastic with and Without Energy Dissipation with Two Temperatures and Rotation. journal of solid mechanics, 840-858.
Kumar, R., Sharma, N., & Lata, P. (2016b). Thermomechanical interactions in transversely isotropic magnetothermoelastic medium with vacuum and with and without energy dissipation with combined effects of rotation, vacuum and two temperatures. Applied Mathematical Modelling, 40, 6560–6575.
Kumar, R., Sharma, N., & Lata, P. (2017). Effects of Hall current and two temperatures in transversely isotropic magnetothermoelastic with and without energy dissipation due to ramp-type heat. Mechanics of Advanced Materials and Structures, 625-635.
Lata, P. (2018). Effect of energy dissipation on plane waves in sandwiched layered thermoelastic medium. Steel and Composite Structures, An Int'l Journal , 27(4).
Lata, P., & Kaur, I. (2018). Effect of hall current in Transversely Isotropic magnetothermoelastic rotating medium with fractional order heat transfer due to normal force. Advances in Materials Research, 7(3), 203–220.
Lata, P., & Kaur, I. (2019a). Transversely isotropic thick plate with two temperature and GN type-III in frequency domain. Coupled Systems Mechanics-Techno Press, 8(1), 55–70.
Lata, P., & Kaur, I. (2019b). Thermomechanical Interactions in Transversely Isotropic Thick Circular Plate with Axisymmetric Heat Supply. Structural Engineering and Mechanics, 69(6), 607–614.
Lata, P., & Kaur, I. (2019c). Transversely isotropic magneto thermoelastic solid with two temperature and without energy dissipation in generalized thermoelasticity due to inclined load. SN Applied Sciences, 1, 426.
Lata, P., & Kaur, I. (2019d). Effect of rotation and inclined load on transversely isotropic magneto thermoelastic solid. Structural Engineering and Mechanics, 70(2), 245–255.
Lata, P., Kumar, R., & Sharma, N. (2016). Plane waves in an anisotropic thermoelastic. Steel & Composite Structures, 22(3), 567–587.
Li, X.-Y., Li, P.-D., Kang, G.-Z., & Pan, D.-Z. (2016). Axisymmetric thermo-elasticity field in a functionally graded circular plate of transversely isotropic material. Mathematics and Mechanics of Solids, 18(5), 464–475.
Liang, J., & Wu, P. (2012). The Refined Analysis of Axisymmetric Transversely Isotropic Cylinder under Radial Direction Surface Loading. Applied Mechanics and Materials, 198-199, 212–215.
Mahmoud, S. (2012). Influence of rotation and generalized magneto-thermoelastic on Rayleigh waves in a granular medium under effect of initial stress and gravity field. Meccanica, Springer, 47, 1561–1579.
Marin, M. (1997a). Cesaro means in thermoelasticity of dipolar bodies. Acta Mechanica, 122(1-4), 155–168.
Marin, M. (1997b). On weak solutions in elasticity of dipolar bodies with voids. Journal of Computational and Applied Mathematics, 82(1-2), 291–297.
Marin, M. (1998). Contributions on uniqueness in thermoelastodynamics on bodies with voids. Revista Ciencias Matemáticas, 16(2), 101–109.
Marin, M. (1999). An evolutionary equation in thermoelasticity of dipolar bodies. Journal of Mathematical Physics, 40(3), 1391–1399.
Marin, M. (2008). Weak solutions in elasticity of dipolar porous materials. Mathematical Problems in Engineering, 1–8.
Marin, M. (2016). An approach of a heat flux dependent theory for micropolar porous media. Meccan., 51(5), 1127–1133.
Marin, M., Agarwal, R. P., & Mahmoud, S. R. (2013). Modeling a Microstretch Thermoelastic Body with Two Temperatures. Abstract and Applied Analysis, 2013, 1–7.
Marin, M., & Baleanu, D. (2016). On vibrations in thermoelasticity without energy dissipation for micropolar bodies. Boundary Value Problems, Springer, 111.
Marin, M., & Öchsner, A. (2017). The effect of a dipolar structure on the Hölder stability in Green–Naghdi thermoelasticity. Continuum Mechanics and Thermodynamics, 29, 1365–1374.
Marin, M., Vlase, S., & Bhatti, M. (2019). On the Partition of Energies for the Backward in Time Problem of Thermoelastic Materials with a Dipolar Structure. Symmetry, 11(7), 863.
Othman, M., & Marin, M. (2017). Effect of thermal loading due to laser pulse on thermoelastic porous medium under G-N theory. Results in Physics, 7, 3863–3872.
Press, W., Teukolshy, S. A., Vellerling, W. T., & Flannery, B. (1986). Numerical recipes in Fortran. Cambridge: Cambridge University Press.
Savruk, M. (1994). Axisymmetric deformation of a transversely isotropic body containing cracks. Materials Science, 29(4), 420–430.
Shahani, A. R., & Torki, H. S. (2018). Determination of the thermal stress wave propagation in orthotropic hollow cylinder based on classical theory of thermoelasticity. Continuum Mechanics and Thermodynamics, Springer, 30(3), 509–527.
Sharma, N., Kumar, R., & Lata, P. (2015a). Effect of two temperature and anisotropy in an axisymmetric problem in transversely isotropic thermoelastic solid without energy dissipation and with two temperature. American Journal of Engineering Research, 4(7), 176–187.
Sharma, N., Kumar, R., & Lata, P. (2015b). Effect Of Two Temperature On The Time Harmonic Behaviour Of An Axisymmetric Problem In Transversely Isotropic Thermoelastic Solid With Green-Naghdi Theory Of Type-II. Afro Asian J SciTech, 199–214.
Shi, T. F., Wang, C. J., Liu, C., Liu, Y., Dong, Y. H., & Li, A. X. (2016). Axisymmetric thermo-elastic field in an annular plate of functionally graded multiferroic composites subjected to uniform thermal loadings. Smart Materials and Structures, 25(3), 1–19.
Slaughter, W. S. (2002). The Linearized Theory of Elasticity. Birkhäuser.
Tarn, J.-Q., Chang, H.-H., & Tseng, W.-D. (2009). Axisymmetric Deformation of a Transversely Isotropic Cylindrical Body: A Hamiltonian State-Space Approach. Journal of Elasticity, 97(2), 131–154.
Vendhan, C. P., & Archer, R. R. (1978). Axisymmetric stresses in transversely isotropic finite cylinders. International Journal of Solids and Structures, 14(4), 305–318.
No fund/scholarship/grant has been taken for this research work.
Department of Basic and Applied Sciences, Punjabi University, Patiala, Punjab, India
Iqbal Kaur & Parveen Lata
Iqbal Kaur
Parveen Lata
The research work is done by IK under the supervision of PL. Both authors read and approved the final manuscript.
Correspondence to Iqbal Kaur.
Kaur, I., Lata, P. Axisymmetric deformation in transversely isotropic magneto-thermoelastic solid with Green–Naghdi III due to inclined load. Int J Mech Mater Eng 15, 3 (2020). https://doi.org/10.1186/s40712-019-0111-8
Transversely isotropic
Magneto-thermoelastic
Mechanical and thermal stresses
Axisymmetric deformation